JP2012204240A

JP2012204240A - 発光領域推定システム、発光領域推定装置、発光領域推定方法および発光領域推定プログラム

Info

Publication number: JP2012204240A
Application number: JP2011069378A
Authority: JP
Inventors: Kazuhiro Wakao; 一広若生; Tadashi Kishimoto; 匡史岸本; Kuniomi Kaga; 国臣加賀
Original assignee: Aomori Support Center for Industrial Promotion
Current assignee: Aomori Support Center for Industrial Promotion
Priority date: 2011-03-28
Filing date: 2011-03-28
Publication date: 2012-10-22

Abstract

【課題】有機ＥＬデバイスにおける発光層内の発光領域を推定することができるという発光位置領域システムを得ること。
【解決手段】有機ＥＬデバイス１の出射角ごとの分光放射強度である角度スペクトルを測定することにより測定角度スペクトルを取得する回転ステージ２および分光放射輝度計４と、有機ＥＬデバイス１の積層方向の複数の発光領域ごとに、有機ＥＬデバイス１の光学特性と仮定した発光領域モデル（発光分布）とに基づいて角度スペクトルを計算した計算値データベースと測定角度スペクトルとに基づいて推定発光領域を求める機器制御装置５と、を備える。
【選択図】図１

Description

本発明は、例えば有機ＥＬ（Electro−Luminescence）デバイス等の自発光デバイスの発光層内の発光領域を推定する発光領域推定システム、発光領域推定装置、発光領域推定方法および発光領域推定プログラムに関する。

有機ＥＬデバイスでは、発光層を含む様々な機能を有する薄膜層が積層されている（例えば、下記特許文献１参照）。有機ＥＬデバイスでは、発光層において電子と正孔が再結合することによって発光が促される。ここで、発光層内において発光する領域は電子と正孔とのバランスから決まり、この領域を発光領域（平面方向ではなく、積層方向の発光分布）と呼ぶ。有機ＥＬデバイスの発光領域は、材料やデバイス構造、製造プロセスの状態等に依存して変化することが知られている。この発光領域の変化は輝度、色度、発光効率や寿命等のデバイス性能に大きく影響する。

有機ＥＬデバイスの輝度・色度調整および製造プロセス最適化のための解析結果のフィードバック、さらには高効率化、長寿命化を図るためには発光層における発光領域を推定することが重要である。

特開２０１１−３５２８６号公報

しかしながら、有機ＥＬデバイスの発光層は数１０ｎｍのオーダーと非常に薄く、有機ＥＬデバイスにおける発光層内の発光領域を直接的に測定するのは困難である。このため、発光層内の発光領域を把握することができない、という問題があった。また能動的に発光領域を決定する手法も確立されておらず、所望の発光位置となるよう材料やデバイス構造、製造プロセス手順等を調整することができない、という問題があった。

本発明は、上記に鑑みてなされたものであって、有機ＥＬデバイスにおける発光層内の発光領域を推定することができる発光領域推定システム、発光領域推定装置、発光領域推定方法および発光領域推定プログラムを得ることを目的とする。

上述した課題を解決し、目的を達成するために、本発明は、発光層を有する積層構造の発光デバイスの出射角ごとの分光放射強度である角度スペクトルを取得する角度スペクトル測定装置と、前記発光デバイスの前記発光層内の積層方向の複数の発光領域ごとに、前記発光デバイスの光学特性と仮定した発光領域モデル（発光分布）とに基づいて角度スペクトルを計算した計算結果と前記測定角度スペクトルとに基づいて推定発光領域を求める発光領域推定装置と、を備えることを特徴とする。

本発明によれば、有機ＥＬデバイスにおける発光層内の発光領域を推定することができるという効果を奏する。また、発光領域を推定することにより、材料や製造プロセス手順等と発光領域の関係を把握することができ、所望の発光領域となるよう材料や製造プロセス手順等を調整することができる。

図１は、発光領域推定システムの構成例を示す図である。図２は、一般的な有機ＥＬデバイスの構成例を示す図である。図３は、発光領域推定方法の概念を示す図である。図４は、発光領域推定手順の一例を示すフローチャートである。図５は、発光領域推定装置（機器制御装置）の機能構成例を示す図である。図６は、評価関数フィルタの一例を説明するための図である。図７は、発光領域をデルタ関数と仮定した場合とガウス分布とした場合との計算結果の一例を示す図である。

以下に、本発明にかかる発光領域推定システム、発光領域推定装置、発光領域推定方法および発光領域推定プログラムの実施の形態を図面に基づいて詳細に説明する。なお、この実施の形態によりこの発明が限定されるものではない。

図１は、本発明にかかる発光領域推定システムの構成例を示す図である。図１に示すように、本実施の形態の発光領域推定システムは、有機ＥＬデバイス１を保持する回転ステージ２と、偏光素子３と、分光放射輝度計４と、機器制御装置（発光領域推定装置）５と、コントローラー６と、電源７と、を備える。

回転ステージ２は、有機ＥＬデバイス１の積層方向１０が回転面に垂直な方向となるよう有機ＥＬデバイス１を固定する機構を有する。有機ＥＬデバイス１には、直流安定化電源等である電源７から電力が供給される。

機器制御装置（発光領域推定装置）５は、例えばパーソナルコンピュータ（ＰＣ）等の制御装置であり、コントローラー６経由で回転ステージ２の回転を制御し、分光放射輝度計４に対して測定の開始や測定条件を設定し、分光放射輝度計４から測定値を取得する等の動作を行うことにより、発光領域推定システム内の各機器を制御する。また、発光領域推定装置５は、自身が計算した、または他の計算機等により計算された計算値データベースを有し、計算値データベースには、Ｓ／Ｐ偏光、発光領域の分布モデルの種類等を条件として計算した分光放射強度の角度依存（角度スペクトル）計算結果が格納されている。

コントローラー６は、機器制御装置５からの指示に基づいて回転ステージ２を回転されることにより、回転ステージ２に固定されている有機ＥＬデバイス１を回転させる。

なお、本実施の形態では、発光領域の推定対象のサンプルとして有機ＥＬデバイス１を例に説明するが、発光領域の推定対象のサンプルは、例えば無機発光デバイス等でもよく、発光層を有する積層された発光デバイスであればよい。

図２は、一般的な有機ＥＬデバイスの構成例を示す図である。図２では、図１の積層方向１０に相当する方向を紙面の上下方向としている。図２は、一般的な例であり、本実施の形態の発光領域の測定対象とする有機ＥＬデバイス２の構成は図２の構成例に限定されない。ここでは、一例として測定対象の有機ＥＬデバイス１が図２の構成例であるとして説明する。

図２の例では、有機ＥＬデバイスは、ガラス基板１１、ＩＴＯ（インジウム・スズ酸化物）陽極１２、ホール注入層１３、ホール輸送層１４、発光層１５、電子輸送層１６、電子注入層１７および金属陰極１８で構成される。

本実施の形態では、発光層１５内の発光領域を推定する。以下では、発光領域を、発光層１５の電子注入側（紙面の下側）を０とし、発光層１５のホール注入側（紙面の上側）を１として正規化して説明する。

図３は、本実施の形態の発光領域推定方法の概念を示す図である。図３の左に示すように、機器制御装置５は、計算値データベースを保持する。そして、図３の右に示すように、機器制御装置５からの制御に基づいて回転ステージ２および分光放射輝度計４を制御することにより、有機ＥＬサンプル１について出射角θごとの分光放射強度を測定する。なお、偏光素子３のＳ偏光とＰ偏光の切替えは、手動で行ってもよいし、機器制御装置５からの指示により自動で切替えられるような機構を備えるようにしてもよい。

ここで、計算値データベースに格納されている分光放射強度の角度依存計算結果の算出方法について説明する。有機ＥＬデバイス２からの出射光の出射角に応じた分光放射強度は、発光層１５内の発光領域に依存する。有機ＥＬデバイス２からの出射光の出射角に応じた分光放射強度の算出は、一般的な多層膜に関する光学シミュレーションソフトウェア等を用いて算出することができる。

光学シミュレーションソフトウェアの入力項目は、例えば、各層の光学定数（屈折率等）、有機ＥＬデバイス１の構造、各層の膜厚、発光層の発光の分光特性（一例として、ＰＬ（Photoluminescence）スペクトル）、発光領域の分布モデル（変化させるパラメータ）等である。発光領域の分布モデルとしては、例えば、デルタ関数、ガウス分布、指数関数分布等、あるいはそれらの組み合わせのうち適切と思われるものを１つ、または複数設定して、それぞれ計算を行っておく。

この光学シミュレーションソフトウェア等としては、どのようなアルゴリズムを用いたものを使用してもよいが、例えば以下のような計算方法を用いることができる。

Ｎ（Ｎは１以上の整数）層からなる多層薄膜のｊ（ｊ＝１，…，Ｎ）層におけるマトリクス成分をＭ_jとすると全層のマトリクス成分Ｍは、Ｍ_jを以下の式（１）で表す場合に、以下の式（２）で表すことができる。ここで、λ_jはｊ層の波長、ｎ_jはｊ層の屈折率、ｄ_jはｊ層の膜厚である。

基板（有機ＥＬデバイス１）のマトリクス成分をη_subとすると、多層膜のフレネル反射係数ρは、Ｂ，Ｃを以下の式（３）のように仮定して、式（４）で表すことができる。

そして、位相変化Φは以下の式（５）で表すことができる。

以上から基板から空気中に放射される発光強度は以下の式（６）で表わすことができる。

ここで、|ρ₀|²、|ρ₁|²は発光層１５の両端の各層における界面のエネルギー反射率、Φ₀、Φ₁は両端の各層における反射時の位相シフト、δ₀、δ₁はそれぞれ発光層全体における位相シフトと発光点と発光層最下面における位相シフトである。

上記の式（６）中のδ₀とδ₁は有機ＥＬデバイスからの発光の出射角θを用いて式（７）で表わすことができる。なお、zは発光層における発光点の位置を示している。

なお、以上述べた出射角θに応じた分光放射強度の算出方法は一例であり、上記の算出方法に限らず薄膜発光デバイスの光学特性を考慮した算出方法であればどのような算出方法を用いてもよい。

図３の左側の計算値データベースの図、および右側の測定値の図では、横軸は出射角（Ａｎｇｌｅ）θを示し、縦軸は正規化した分光放射強度を示している。ここでは、出射角θが０［ｄｅｇ］のときの波長５９０ｎｍの分光放射強度を１として正規化している。なお、正規化の基準は、一例であり、他の波長等を基準として正規化してもよい。

また、図３の左側の計算値データベースの図では、上側がＰ偏光についての計算結果を示し下側がＳ偏光についての計算結果を示している。また、図３の左側の計算値データベースの図では、発光領域を分布の無い、点で示される位置とするデルタ関数で仮定し、発光層１５内の電子注入側を０としホール注入側を１とした場合の正規化した発光位置（発光領域）を０．１刻みで計算しておくとし、発光位置０．１，０．２，…，０．９をそれぞれＮＯ１，ＮＯ２，…，ＮＯ９と番号を付け、ＮＯ１，ＮＯ５，ＮＯ９についてグラフ化している。なお、図３では、図の簡略化のため３つの発光位置のみを示しているがデータベースとしてはこれら以外の番号の位置のデータも格納されている。また、波長については、５６０、５７０、５９０、６１０、６４０ｎｍについてそれぞれグラフ化している。計算値データベースとして計算しておく波長範囲は、測定対象の有機ＥＬデバイス１の発光する波長を含むように余裕をもって適切に設定しておく。なお、ここでは、有機ＥＬデバイス１の発光波長が可視光の例を示しているが、有機ＥＬデバイス１の発光波長は可視光に限定されない。また、計算値データベースとして計算しておく波長刻み、および発光位置の刻みは、データベースの容量や要求される発光位置の推定精度等に応じて設定すればよく、図３に例示した値に限定されない。

また、図３の右側の測定値の図は、上側がＳ偏光についての測定結果を示し下側がＰ偏光についての測定結果を示しており、５７０、５８０、６００、６１０、６２０、６３０、６４０ｎｍの分光放射強度を正規化した値を示している。なお、図３では、横軸を出射角とした図を示しているが、実際には、測定結果としては、まず、各出射角において、波長ごとの分光放射強度が得られる。そして、出射角を変えて、それぞれ分光放射強度を測定し、それらの測定結果から図３に示したように、出射角ごとの分光放射強度を波長ごとに求める。

そして、計算値データベースのうち、発光位置（領域）推定のために最適な波長／偏光成分に対応する値を抽出し、抽出した値と測定結果とを比較してフィッティングすることにより、発光位置（領域）を推定する。すなわち、抽出した値と測定結果とが最も近くなるような発光位置（領域）を求める。

図３の中央に示した図では、Ｐ偏光について計算値データベースの値として波長５９０ｎｍのＮＯ１，ＮＯ３，ＮＯ５，ＮＯ７，ＮＯ９の発光位置に対応する値をグラフ化しており、点線は５９０ｎｍの測定結果を示している。なお、図３では、一例としては発光分布をデルタ関数で示される位置として計算例を示している。図３の例では、測定結果は、上から３番目の計算値であるＮＯ５の発光位置とした計算値に最も近く、この例では、発光位置としてＮＯ５に対応する位置が発光位置（領域）の推定結果として求められる。

図４は、本実施の形態の発光領域推定手順の一例を示すフローチャートである。また、図５は、機器制御装置５の機能構成例を示す図である。図５に示すように、機器制御装置５は、記憶部５１と、ＣＰＵ（Central Processing Unit）等の制御部５２と、を備える。記憶部５１には、計算値データベース５３が格納される。制御部５２は、計算値データベースを生成するデータベース生成部５４と、回転ステージ２や分光放射輝度計４等を制御する測定制御部５５と、最適波長抽出部５６と、フィッティング部５７と、発光領域推定部５８と、を備える。データベース生成部５４、測定制御部５５、最適波長抽出部５６、フィッティング部５７および発光領域推定部５８は、外部記憶やネットワーク回線等によりプログラムがインストールされることにより、制御部５２上で機能する。

以下、図４および図５を用いて本実施の形態の発光領域推定手順を説明する。図４に示すように、機器制御装置５のデータベース生成部５４または他の計算機等は、有機ＥＬデバイス１の発光領域をパラメータとした出射角に応じた分光放射強度（角度スペクトル特性）を計算し、計算結果を計算値データベースに格納する（ステップＳ１）。他の計算機等で計算を行った場合は、データベース生成部５４は計算結果を取得して計算値データベースに格納する。

次に、機器制御装置５の測定制御部５５は、回転ステージ２および分光放射輝度計４等を制御することにより測定値（Ｓ／Ｐ偏光ごとの角度スペクトル特性）を取得する（ステップＳ２）。そして、最適波長抽出部５６は、計算値データベースに基づいて、フィッティング処理で用いる最適な波長および偏光成分を、評価関数フィルタを用いて抽出する（ステップＳ３）。この際、最適な波長および偏光成分は１つ抽出してもよいし複数抽出してもよい。なお、ステップＳ２とステップＳ３の順序は逆であってもよい。また、ステップＳ３をステップＳ１の前に実施してもよい。

ここで、ステップＳ３の抽出処理に用いる評価関数フィルタについて説明する。図６は、評価関数フィルタの一例を説明するための図である。ここでは、評価関数フィルタを、計算値データベースに格納されている角度スペクトル特性の分散値を求め、分散値の積分値が最も大きい波長を選択するようなフィルタとする。

図６では、波長４６０ｎｍ、発光領域をデルタ関数で仮定した場合の発光位置０．１，０．５，０．９（ＮＯ１，ＮＯ５，ＮＯ９）の例を図示している。図６に示すように、Ｓ偏光成分（左側の上図）について、角度（出射角）ごとに分散値を求める（左側の下図）。同様に、Ｐ偏光成分（右側の上図）について、角度（出射角）ごとに分散値を求める（右側の下図）。そして、分散値の角度に関する積分値（面積）を求める。偏光成分ごとに各波長について、この積分値を求め、積分値が大きくなる波長および偏光成分の組み合わせを最適な波長および偏光成分として抽出（選択）する。この際、この積分値が大きくなる１つの波長および偏光成分の組み合わせを最適な波長および偏光成分として選択してもよいが、例えば積分値の大きい順に複数の波長や偏光成分を最適な波長および偏光成分として選択してもよい。また、最適な波長および偏光成分を選択する際に、必ずしも積分値の大きい順に選択する必要はない。例えば、最適な波長および偏光成分として選択された波長で分光放射強度の値が小さい場合には測定値の精度が低くなるため、発光領域の推定精度も低下する可能性がある。従って例えば、分光放射強度における測定値の絶対値精度も考慮して最適な波長および偏光成分を選択してもよい。

なお、ここでは、分散値の積分値に基づいて最適波長を選択する例を説明したが、発光領域の差の影響を受けやすい波長を最適波長として選択する方法であれば、これに限らずどのような方法で選択してもよい。

図４の説明に戻り、最適な波長および偏光成分（以下、偏光成分も含めて最適波長という）を抽出した後、フィッティング部５７は、測定値に対する計算値のフィッティングを行い（測定値に最も近い計算値を選び）、最適波長ごとの発光位置（領域）を推定する（ステップＳ４）。この際のフィッティングの具体的方法としては、例えば、最適波長ごとに、計算値データベースの各発光位置（計算で仮定した発光領域モデル）について測定値と計算値との相関値等の評価係数を求めて、例えば相関値の大きい（最も相関が強い）発光位置を当該最適波長の推定発光位置とする。なお、フィッティングの具体的方法はこれに限定されない。また、本実施の形態では、計算した発光位置の刻みを最小分解能として発光位置を求めているが、補間等の適切な方法により計算した刻みより小さい分解能で発光位置を求めてもよい。

次に、フィッティング部５７は、全ての最適波長についてフィッティング処理を実施したか否かを判断し（ステップＳ５）、フィッティングを実施していない最適波長がある場合（ステップＳ５Ｎｏ）、ステップＳ４へ戻る。

発光領域推定部５８は、全ての最適波長についてフィッティング処理を実施した場合（ステップＳ５Ｙｅｓ）、最適波長ごとの推定した発光位置（領域）に基づいて推定発光位置を求める（ステップＳ６）。具体的には、例えば、多数決により推定発光位置を求める。多数決により推定発光位置をする場合、具体的には例えば３つの最適波長を抽出するとしたとき、最適波長ごとに推定された３つの発光位置のうち２つ以上が一致していれば、その一致した発光位置を推定発光位置とする。なお、最適波長ごとの推定した発光位置に基づいて推定発光位置を求める方法はこれに限らず、例えば、平均値を用いる等どのような方法としてもよい。

以上では、発光領域をデルタ関数として分布の無い発光位置のみを推定する例を説明したが、発光位置だけでなく分布を持つモデルについても推定を行ってもよい。この場合、計算値データベースとして、分布を持つ発光領域モデルについても複数のモデルについて計算しておく。

図７は、発光領域を点とした（デルタ関数）場合とガウス分布とした場合との計算結果の一例を示す図である。図７では、上側に発光分布を幅５ｎｍのガウス分布（ｅｘｐ（−（（ｚ−ｚ₀）／Ｗ）²／２）：ｚ＝発光層における発光点の位置、ｚ₀＝発光分布中心、Ｗ＝５ｎｍ）とした場合の角度スペクトルの計算結果を示し、下側に発光分布を点とした（デルタ関数状の分布）場合の計算結果を示している。いずれの場合も、波長５９０ｎｍ、発光位置を０．１刻みで計算した例を示している。このように、発光位置は同じでも発光分布が異なると角度スペクトルも異なる。従って、例えば、発光分布の幅を固定して発光（中心）位置を推定した後に、分布モデルが異なる計算データベースの計算値に基づいて測定値に対してさらにフィッティングを行い、測定値に最も近い発光領域を求めるようにしてもよい。

また、さらに分布の幅を変更して、測定値に対してさらにフィッティングを行い、測定値に最も近い幅を求めるようにしてもよい。

なお、本実施の形態では、回転ステージ２、偏光素子３、分光放射輝度計４およびコントローラー６を用いて、有機ＥＬデバイス１の角度スペクトルを取得している。すなわち、有機ＥＬデバイスの角度スペクトル（出射角ごとの分光放射強度）を取得する角度スペクトル測定装置が、回転ステージ２、偏光素子３、分光放射輝度計４およびコントローラー６で構成されている。角度スペクトル測定装置は、この構成に限らず、角度スペクトルが測定できる構成であればどのような構成でもよい。

以上のように、本実施の形態では、発光領域をパラメータとして発光領域ごとに角度スペクトルを計算した結果を計算値データベースとして保持し、一方で、測定対象の有機ＥＬデバイス１について偏光成分ごとに角度スペクトルを測定して測定結果を取得する。そして、計算値データベースから評価関数フィルタを用いて最適波長を抽出して、抽出した波長に対応する計算値データベース上の計算値と測定値とをフィッティングすることにより、推定発光領域を求めるようにした。そのため、有機ＥＬデバイスにおける発光層内の発光領域を推定することができる。また、評価関数フィルタを用いて最適波長を抽出して、抽出した波長に対応する計算値を用いていることにより、計算値データベースの全てのケースの計算結果と測定値との比較を実施しなくてすみ、計算負荷を抑制することができる。

１有機ＥＬデバイス
２回転ステージ
３偏光素子
４分光放射輝度計
５機器制御装置
６コントローラー
７電源
１０積層方向
５１記憶部
５２制御部
５３計算値データベース
５４データベース生成部
５５測定制御部
５６最適波長抽出部
５７フィッティング部
５８発光領域推定部

Claims

発光層を有する積層構造の発光デバイスの出射角ごとの分光放射強度である角度スペクトルを測定することにより測定角度スペクトルを取得する角度スペクトル測定装置と、
前記発光デバイスの前記発光層内の積層方向の複数の発光領域ごとに、前記発光デバイスの光学特性と仮定した発光領域に基づいて角度スペクトルを計算した計算結果と前記測定角度スペクトルとに基づいて推定発光領域を求める発光領域推定装置と、
を備えることを特徴とする発光領域推定システム。
前記発光領域推定装置は、前記計算結果として計算した波長のうちから最適波長を選択し、前記最適波長に対応する前記計算結果に基づいて前記推定発光領域を求める、ことを特徴とする請求項１に記載の発光領域推定システム。
前記角度スペクトルを偏光成分ごとに測定し、前記計算結果を偏光成分ごとに計算する、ことを特徴とする請求項１に記載の発光領域推定システム。
前記発光領域推定装置は、前記計算結果として計算した波長および偏光成分のうちから最適な波長および偏光成分の組み合わせを最適波長として選択し、前記最適波長に対応する前記計算結果に基づいて前記推定発光領域を求める、ことを特徴とする請求項３に記載の発光領域推定システム。
前記最適波長を１つとする、ことを特徴とする請求項２または４に記載の発光領域推定システム。
前記最適波長を複数とする、ことを特徴とする請求項２または４に記載の発光領域推定システム。
前記発光領域推定装置は、波長ごとに前記計算結果の分散値を求め、前記分散値に基づいて前記最適波長を選択する、ことを特徴とする請求項２、４〜６のいずれか１つに記載の発光領域推定システム。
前記発光領域推定装置は、前記分散値の前記出射角についての積分値に基づいて前記最適波長を選択する、ことを特徴とする請求項７に記載の発光領域推定システム。
前記発光領域推定装置は、前記最適波長ごとに前記計算結果と前記測定角度スペクトルとに基づいて発光領域を求め、前記最適波長ごとに求めた発光領域のうち発光領域が一致した前記最適波長が多い発光領域を前記推定発光領域とする、ことを特徴とする請求項１〜８のいずれか１つに記載の発光領域推定システム。
前記発光領域推定装置は、前記計算結果と前記測定角度スペクトルとの差に基づいて前記推定発光領域を求める、ことを特徴とする請求項１〜９のいずれか１つに記載の発光領域推定システム。
前記発光領域推定装置は、前記計算結果と前記測定角度スペクトルとの相関値を求め、前記相関値に基づいて前記推定発光領域を求める、ことを特徴とする請求項１０に記載の発光領域推定システム。
発光層を有する積層構造の発光デバイスの出射角ごとの分光放射強度である角度スペクトルを測定することにより測定角度スペクトルを取得する角度スペクトル測定装置を制御し、前記角度スペクトルの測定値である測定角度スペクトルを前記角度スペクトル測定装置から取得する測定制御部と、
前記発光デバイスの前記発光層内の積層方向の複数の発光領域ごとに、前記発光デバイスの光学特性と仮定した発光領域モデル（発光分布）とに基づいて角度スペクトルを計算した計算結果を保持する記憶部と、
前記計算結果と前記測定角度スペクトルとに基づいて推定発光領域を求めるフィッティング部と、
を備えることを特徴とする発光領域推定装置。
発光層を有する積層構造の発光デバイスの出射角ごとの分光放射強度である角度スペクトルを測定することにより測定角度スペクトルを取得する角度スペクトル測定ステップと、
前記発光デバイスの前記発光層内の積層方向の複数の発光領域ごとに、前記発光デバイスの光学特性と仮定した発光領域モデル（発光分布）とに基づいて角度スペクトルを計算した計算結果と前記測定角度スペクトルとに基づいて推定発光領域を求める発光領域推定ステップと、
を含むことを特徴とする発光領域推定方法。
請求項１３に記載の発光領域推定方法をコンピュータに実行させることを特徴とする発光領域推定プログラム。