JP6443089B2 - 屋外グレア評価装置、屋外グレア評価方法及び屋外グレア評価プログラム - Google Patents

Description

本発明は、屋外の不快グレアを評価する屋外グレア評価装置、屋外グレア評価方法及び屋外グレア評価プログラムに関する。

従来、屋内の不快グレアを評価するグレア評価システムが知られている（例えば、特許文献１及び２参照）。このグレア評価システムでは、不快グレアの値を、光源輝度、背景輝度、立体角、位置指数をパラメータとしたＵＧＲ(Unified Glare Rating)として求めている。

特開平５−７２０３２号公報 特開２０００−１５５０４９号公報

屋外のグレアを評価する手法は確立されておらず、歩行者空間のグレア評価に、投光照明された屋外スポーツ施設の不快グレアを数値化するＧＲ(Glare Rating)が流用されている。
しかしながら、発明者は、ＧＲから算出されるＧＲ値と歩行者が主観評価した結果とに相関が取れないという知見を得た。
本発明は、上述した事情に鑑みてなされたものであり、屋外の不快グレアを評価可能な屋外グレア評価装置、屋外グレア評価方法及び屋外グレア評価プログラムを提供することを目的とする。

上述した目的を達成するために、本発明の屋外グレア評価システムは、少なくとも１つの屋外用照明器具を含めて撮影した輝度分布データから、屋外用照明器具の輝度分布を順応輝度に応じて重み付け加算した有効グレア輝度Ｌｅｆｆ、立体角ω及び位置指数Ｐを求め、前記有効グレア輝度Ｌｅｆｆ、前記立体角ω、前記位置指数Ｐおよび定数α、β、γに基づいて、屋外グレア評価値ｍ＿ＤＧＩを次式によって算出するグレア算出部を備え、前記有効グレア輝度Ｌｅｆｆを、輝度分布データに対して輝度の範囲を指定して求めたことを特徴とする。
ｍ＿ＤＧＩ＝ｌｏｇ（Ｌｅｆｆα×ωβ／Ｐγ）

上述の構成において、前記範囲は、輝度分布データにおける最大輝度Ｌｍａｘの１／１０以上としてもよい。

また、本発明の屋外グレア評価方法は、少なくとも１つの屋外用照明器具を含めて撮影した輝度分布データから、屋外用照明器具の輝度分布を順応輝度に応じて重み付け加算した有効グレア輝度Ｌｅｆｆ、立体角ω及び位置指数Ｐを求め、前記有効グレア輝度Ｌｅｆｆ、前記立体角ω、前記位置指数Ｐおよび定数α、β、γに基づいて、屋外グレア評価値ｍ＿ＤＧＩを次式によって算出し、前記有効グレア輝度Ｌｅｆｆを、輝度分布データに対して輝度の範囲を指定して求めることを特徴とする。
ｍ＿ＤＧＩ＝ｌｏｇ（Ｌｅｆｆα×ωβ／Ｐγ）

また、本発明の屋外グレア評価プログラムは、少なくとも１つの屋外用照明器具を含めて撮影した輝度分布データから、屋外用照明器具の輝度分布を順応輝度に応じて重み付け加算した有効グレア輝度Ｌｅｆｆ、立体角ω及び位置指数Ｐを求め、前記有効グレア輝度Ｌｅｆｆ、前記立体角ω、前記位置指数Ｐおよび定数α、β、γに基づいて、屋外グレア評価値ｍ＿ＤＧＩを次式によって算出し、前記有効グレア輝度Ｌｅｆｆを、輝度分布データに対して輝度の範囲を指定して求めることを特徴とする。
ｍ＿ＤＧＩ＝ｌｏｇ（Ｌｅｆｆα×ωβ／Ｐγ）

本発明によれば、屋外用照明器具の輝度分布を考慮した有効グレア輝度に基づいて屋外グレア評価値を算出したため、屋外の不快グレアを評価できる。

本発明の実施形態に係る屋外グレア評価システムの構成を概略的に示す図である。 実験の条件を示す説明図である。 ＧＲ値と主観評価結果との関係を示すグラフである。 屋外グレア評価値と主観評価結果との関係を示すグラフである。 屋外グレア評価装置の機能的構成を示すブロック図である。 輝度分布画像を示す図である。 屋外グレア評価処理を示すフローチャートである。

以下、図面を参照して本発明の実施形態について説明する。
図１は、本実施形態に係る屋外グレア評価システム１００の構成を概略的に示す図である。
屋外グレア評価システム１００は、カメラ７と、屋外グレア評価装置１１と、表示装置１３とを備えている。
カメラ７は、少なくとも１つの評価対象を含む光環境３を撮像し、撮像画像を生成する撮像手段である。撮像画像は、歩行者の視線で、すなわち、所定の歩行者の目の高さ（約１．５ｍ）且つ歩行方向の角度（水平又は斜め下方視）で、評価対象を含めて撮影される。撮像画像は、カメラ７によりデジタルデータ化され、或いは、別途のスキャナ装置などを用いてデジタルデータ化されて撮像画像データ５として屋外グレア評価装置１１に入力される。

屋外グレア評価装置１１は、撮像画像データ５に基づいて屋外グレア評価値ｍ＿ＤＧＩを算出し、算出した屋外グレア評価値ｍ＿ＤＧＩを表示装置１３に出力する。
表示装置１３は、屋外グレア評価値ｍ＿ＤＧＩの出力手段として用いられ、屋外グレア評価値ｍ＿ＤＧＩを表示する表示手段である。なお、出力手段としては、表示装置１３の代わりに、プリンタ装置などを用いてもよい。

本実施形態の評価対象である屋外用照明器具１は、図示しない固定手段により、路面Ｓから所定の高さ（例えば、４．５ｍ）に支持され、路面Ｓを照明する照明器具であり、屋外用照明器具１の光源はＬＥＤ等の発光素子やランプ等、種々の光源が用いられる。本実施形態の屋外用照明器具１は、歩行者の目の高さＴ２（約１．５ｍ）よりも高い高さＴ１に配置されて所定照度で路面Ｓを照明する街路灯として構成されている。所定照度は、街路灯に必要な照度に設定される。ここで、高さＴ１は、路面Ｓから屋外用照明器具１の発光部の中心までの高さとする。また、屋外用照明器具１の発光部は、光源を含む発光部分とし、反射鏡やグローブ等の光学素子を有する場合には当該光学素子も含むものとする。

屋外グレア評価システム１００の詳細説明に先立って、屋外のグレア（眩しさ）について説明する。
屋外用照明器具に対しては、眩しさを軽減することが望まれているが、屋外のグレアを評価する手法は確立されておらず、歩行者空間のグレア評価に、投光照明された屋外スポーツ施設の不快グレアを評価する評価式であるＧＲ(Glare Rating)が流用されている。このＧＲは、投光照明された屋外スポーツ施設で行われた主観評価実験から得られた実験式である。
しかしながら、発明者は、実験により、ＧＲから得られたＧＲ値と歩行者が主観評価した結果（主観評価結果）とに相関が取れないという知見を得た。具体的には、発明者は、どのような場合に人が眩しさを感じるかを調べるため、被験者に、複数の評価対象の眩しさについて、複数の観測位置において、所定の歩行者の目の高さ（約１．５ｍ）且つ歩行方向の角度（水平視）で主観評価させる実験を行った。この実験では、評価対象をランプ又はＬＥＤを含む７つの屋外用照明器具（Ｌｉｇｈｔ−１〜７）として、これら屋外用照明器具の眩しさを被験者１６人で評価した。実験の条件を図２に示す。なお、図２（Ａ）は屋外用照明器具の配置位置及び観測位置を示す図であり、図２（Ｂ）は観測位置の観測距離及び鉛直角を示す図であり、図２（Ｃ）は主観評価スケールを示す図である。なお、図２（Ａ）中、符号２は、屋外用照明器具を支持する支柱である。

屋外用照明器具は、図２（Ａ）に示すように、幅員５ｍの道路の路肩に、路面から４．５ｍの高さに設置し、観測位置は、車線軸（２．５ｍ）上のＡ〜Ｆの６カ所とする。観測位置Ａにおける観測距離ＬＡは３４．３ｍ、鉛直角は８５°であり、観測位置Ｂにおける観測距離ＬＢは２４．４ｍ、鉛直角は８３°であり、観測位置Ｃにおける観測距離ＬＣは１７．０ｍ、鉛直角は８０°であり、観測位置Ｄにおける観測距離ＬＤは１１．２ｍ、鉛直角は７５°であり、観測位置Ｅにおける観測距離ＬＥは６．４ｍ、鉛直角は６５°であり、観測位置Ｆにおける観測距離ＬＦは３．６ｍ、鉛直角は５０°である。ここで、観測距離Ｌ（ＬＡ〜ＬＦ）は、器具中心から歩行者の目までの距離を示す。また、この実験では、図２（Ｃ）に示すように、９段階の評価スケールを用いており、１はまぶしくない、３はややまぶしい、５はまぶしい、７は非常にまぶしい、９はまぶしさが耐えがたい、を表す。

実験によれば、観測位置Ａから観測位置Ｄでは、観測距離が短くなるほど被験者がグレアを強く感じる傾向にあった。また、屋外用照明器具に最も近い観測位置Ｆでは、被験者の視界からグレア源となる屋外用照明器具が外れて被験者がグレアを感じず、次に近い観測位置Ｅにおいても、同様の理由で一部の被験者がグレアを感じなかった。したがって、データから観測位置Ｅ，Ｆの結果を除外した。
また、同一の評価対象及び条件（観測位置はＡ〜Ｄ）において、ＧＲを用いてＧＲ値を算出した。ＧＲは、撮像画像から写真測光法により得た輝度分布図から算出したグレア光による等価光幕輝度Ｌｖｌと反射光による等価光幕輝度Ｌｖｅとに基づいて、次式（１）と表せる。
ＧＲ＝２７＋２４×ｌｏｇ（Ｌｖｌ／Ｌｖｅ^０．９）・・・（１）
実験及びＧＲ値の結果を図３に示す。

図３は、ＧＲ値と主観評価結果との関係を示すグラフである。図３では、横軸にＧＲ値、縦軸に主観評価結果を示す。ここで、ＧＲ値のスケールは１０〜９０の９段階、主観評価結果のスケールは１〜９の９段階であり、図３では、ＧＲ値は２０〜６０を、主観評価結果は１〜６を示している。なお、本実施形態では、説明の便宜上、評価式ＧＲで得られるＧＲ値にも符号ＧＲを付すものとする。
図３に示すように、ＧＲ値と主観評価結果とに相関が取れておらず、ｙ（主観評価結果）＝ｘ（ＧＲ値）の一次関数を示す直線Ｈ上に結果が現れていない。また、データは、発光素子を有する等して発光部の輝度分布が比較的不均一なグループＧ１と、ランプや拡散グローブを有する等して発光部の輝度分布が比較的均一なグループＧ２と、に分かれて分布している。また、発光部の輝度分布が比較的不均一な屋外用照明器具は眩しいと評価され、発光部の輝度分布が比較的均一な屋外用照明器具ほど眩しさ評価が低くなっている。この結果は、ＧＲ値が、発光部の輝度分布が比較的不均一な光源と、発光部の輝度分布が比較的均一な光源とに共通して使える指標ではないことを示している。

この結果に対し、発明者は、有効グレア輝度Ｌｅｆｆ、立体角ω、位置指数Ｐをパラメータとした式によって屋外グレア評価値ｍ＿ＤＧＩを算出することで、屋外グレア評価値ｍ＿ＤＧＩと主観評価結果とがｙ（主観評価結果）＝ｘ（屋外グレア評価値ｍ＿ＤＧＩ）の一次関数で近似できるという知見を得た。この屋外グレア評価値ｍ＿ＤＧＩは、有効グレア輝度Ｌｅｆｆ、立体角ω、位置指数Ｐ、定数α、β、γに基づいて、次式（２）の評価式ｍ＿ＤＧＩと表せる。なお、本実施形態では、評価式ｍ＿ＤＧＩで得られる屋外グレア評価値にも符号ｍ＿ＤＧＩを付すものとする。
ｍ＿ＤＧＩ＝ｌｏｇ（Ｌｅｆｆ^α×ω^β／Ｐ^γ）・・・（２）

ここで、有効グレア輝度Ｌｅｆｆとは、被験者が照明器具を直視しないで不快グレアを評価した結果と近似する値であり、有効グレア輝度Ｌｅｆｆは、順応輝度に応じて、発光部の輝度分布データを重み付け加算することで求まる。
順応輝度０．１ｃｄ／ｍ^２のときの重み付け関数ＷＦ_０．１は、評価対象を撮像した撮像画像の各ピクセルの位置をｉ行目及びｊ列目で示したピクセル位置ijとした場合、ピクセル位置ijにおけるピクセル輝度Ｌ_ｉｊに基づいて次式（３）と表せる。
ＷＦ_０．１＝７．０×１０^−４×Ｌ_ｉｊ ^{０．６０５}・・・（３）
順応輝度１．０ｃｄ／ｍ^２のときの重み付け関数ＷＦ_１．０は、次式（４）と表せる。
ＷＦ_１．０＝６．０×１０^−４×Ｌ_ｉｊ ^{０．５８８}・・・（４）
順応輝度１０ｃｄ／ｍ^２のときの重み付け関数ＷＦ_１０は、次式（５）と表せる。
ＷＦ_１０＝４．０×１０^−４×Ｌ_ｉｊ ^{０．４７１}・・・（５）
そして、有効グレア輝度Ｌｅｆｆは、総ピクセル数Ａｅｆｆを用いて、次式（６）と表せる。
Ｌｅｆｆ＝Σ（Ｌ_ｉｊ×ＷＦ_Ｌｂ（Ｌ_ｉｊ））／Ａｅｆｆ・・・（６）
ここで、重み付け関数ＷＦ_Ｌｂ（Ｌ_ｉｊ）はそのときの順応輝度に応じて重み付け関数ＷＦ_０．１，ＷＦ_１．０，ＷＦ_１０のいずれかから選択される。

図４は、屋外グレア評価値ｍ＿ＤＧＩと主観評価結果との関係を示すグラフである。図４では、横軸に屋外グレア評価値ｍ＿ＤＧＩ、縦軸に主観評価結果を示す。また、図４の屋外グレア評価値ｍ＿ＤＧＩは、有効グレア輝度Ｌｅｆｆとして、輝度データのうち最大輝度Ｌｍａｘの１／１０以上の輝度値を対象に求めた有効グレア輝度Ｌｅｆｆ＿１０を用い、定数αを１．５、定数βを０．５、定数γを０．８とした次式（７）によって算出された値である。定数α，β，γは、上述の実験の結果から求まる値である。
ｍ＿ＤＧＩ＝ｌｏｇ（Ｌｅｆｆ＿１０^１．５×ω^０．５／Ｐ^０．８）・・・（７）
図４に示すように、屋外グレア評価値ｍ＿ＤＧＩと主観評価結果とは、相関係数Ｒが０．９３となり、高い正の相関関係を有する。
そこで、本実施形態の屋外グレア評価装置１１は、評価式ｍ＿ＤＧＩを用いて屋外グレア評価値ｍ＿ＤＧＩを算出するように構成されている。

図５は、屋外グレア評価装置１１の機能的構成を示すブロック図である。図６は、輝度分布画像を示す図である。
屋外グレア評価装置１１は、例えば、ＣＰＵやＲＡＭ、ＲＯＭ等の記憶部を備えたコンピュータを用いて構成され、図５に示すように、画像入力部２１と、輝度分布画像生成部２３と、パラメータ算出部２５と、グレア算出部２７と、グレア評価部２８と、結果出力部２９とを備え、これらの各機能ブロックは、屋外グレア評価プログラム１５（図１）をコンピュータが実行することで実現される。
画像入力部２１は、カメラ７や別途のスキャナ装置等から撮像画像データ５が入力される入力手段である。

輝度分布画像生成部２３は、撮像画像データ５に基づいて、光環境３の輝度分布を示す輝度分布画像（輝度分布データ）４１を生成する。輝度分布画像生成部２３による輝度分布画像４１の生成としては、例えば写真測光法（photographic photometry）を用いることができる。すなわち、輝度分布画像生成部２３は、露出条件を異ならせて光環境３を撮像した複数の撮像画像データ５を、それぞれ画像の階調値を輝度値に変換して輝度データを生成し、各撮像画像データ５の輝度データを合成して輝度分布画像４１（図６）を生成する。

パラメータ算出部２５は、輝度分布画像４１に基づいて、有効グレア輝度Ｌｅｆｆ、立体角ω及び位置指数Ｐを含むパラメータを算出する算出部であり、有効グレア輝度Ｌｅｆｆを求める範囲を設定する範囲設定部２６を備えている。
範囲設定部２６は、輝度分布画像４１に対して、有効グレア輝度の計算領域として所定範囲を設定する設定手段であり、ユーザから入力部３１を介して所定範囲の指定を受け、当該所定範囲を記憶する。所定範囲は、最大輝度Ｌｍａｘを基準に設定される場合と、所定輝度を基準に設定される場合とがある。

例えば、最大輝度Ｌｍａｘの１／１０以上の輝度値を計算対象とした場合には、計算領域は図６（Ｂ）のようになる。
また、例えば、最大輝度Ｌｍａｘの１／１００以上の輝度値を計算対象とした場合には、計算領域は図６（Ｃ）のようになる。
さらに、例えば、１０００ｃｄ／ｍ^２以上の輝度値を算出対象とした場合、対象領域は図６（Ｄ）のようになる。
図６の例（Ｌｉｇｈｔ−６）は最大輝度Ｌｍａｘが３８万ｃｄ／ｍ^２であるので、閾値は、図６（Ｂ）では３８，０００ｃｄ／ｍ^２、図６（Ｃ）では３，８００ｃｄ／ｍ^２となる。
なお、図６（Ａ）〜図６（Ｃ）は、計算領域の説明のために、それぞれ最大輝度Ｌｍａｘの１／１０、最大輝度Ｌｍａｘの１／１００、１０００ｃｄ／ｍ^２を閾値として輝度分布画像４１を二値化した画像であり、有効グレア輝度を求める際には二値化は行わない。

計算対象の最小輝度値を高くするほど、屋外グレア評価値ｍ＿ＤＧＩと歩行者の主観評価結果との相関をより高めることができ、屋外の不快グレアをより正確に評価できる。本実施形態では、輝度分布画像４１の最大輝度Ｌｍａｘを基準とし、最大輝度Ｌｍａｘの１／１０以上の輝度値を計算対象として有効グレア輝度Ｌｅｆｆ＿１０を求めている。

グレア算出部２７は、パラメータに基づいて、式（７）によって屋外グレア評価値ｍ＿ＤＧＩを算出する算出部である。
グレア評価部２８は、屋外グレア評価値ｍ＿ＤＧＩから屋外グレアを評価するものである。グレア評価部２８は、例えば、図２（Ｃ）に示す９段階の評価スケール（０以上２未満はまぶしくない、２以上４未満はややまぶしい、４以上６未満はまぶしい、６以上８未満は非常にまぶしい、８以上はまぶしさが耐えがたい）を用いて評価する。
結果出力部２９は、屋外グレア評価値ｍ＿ＤＧＩ及び評価結果を上述の表示装置１３等の出力装置に出力する出力手段である。

次いで、屋外グレア評価処理を図７に示すフローチャートを参照して説明する。
この図に示すように、屋外グレア評価装置１１は、画像入力部２１を介して複数の撮像画像データ５が入力されると（ステップＳ１）、輝度分布画像生成部２３が写真測光法に基づいて輝度分布画像４１を生成する（ステップＳ２）。

次に、パラメータ算出部２５が輝度分布画像４１から有効グレア輝度Ｌｅｆｆ（本実施形態では、有効グレア輝度Ｌｅｆｆ＿１０）、立体角ω及び位置指数Ｐを含むパラメータを求める（ステップＳ３）。
次いで、グレア算出部２７が式（７）によって屋外グレア評価値ｍ＿ＤＧＩを算出し（ステップＳ４）、グレア評価部２８が屋外グレア評価値ｍ＿ＤＧＩに基づいて屋外グレアを評価し、評価結果を生成する（ステップＳ５）。
そして、結果出力部２９が屋外グレア評価値ｍ＿ＤＧＩ及び評価結果を表示装置１３に出力する（ステップＳ６）。この表示装置１３の表示により、光環境３の屋外グレアが定量的に評価可能になる。

以上説明したように、本実施形態によれば、少なくとも１つの屋外用照明器具１を含めて撮影した輝度分布画像４１から、屋外用照明器具１の輝度分布を順応輝度に応じて重み付け加算した有効グレア輝度Ｌｅｆｆ、立体角ω及び位置指数Ｐを含むパラメータを求め、当該パラメータから屋外グレア評価値ｍ＿ＤＧＩを算出するグレア算出部２７を備える構成とした。具体的には、有効グレア輝度Ｌｅｆｆ、立体角ω、位置指数Ｐ、実験によって求まる定数α、β、γに基づいて、屋外グレア評価値ｍ＿ＤＧＩを次式（２）によって求める構成とした。
ｍ＿ＤＧＩ＝ｌｏｇ（Ｌｅｆｆ^α×ω^β／Ｐ^γ）・・・（２）
この構成により、屋外用照明器具１の輝度分布を考慮した有効グレア輝度Ｌｅｆｆに基づいて屋外グレア評価値ｍ＿ＤＧＩを算出できるので、屋外の不快グレアを適正に評価できる。また、発光部の輝度分布が比較的不均一な屋外用照明器具１と均一な屋外用照明器具１とを区別することなく、不快グレアを数値化できる。

また、本実施形態によれば、有効グレア輝度Ｌｅｆｆを、輝度分布データに対して所定範囲を指定して有効グレア輝度Ｌｅｆｆ＿１０として求めたため、屋外グレア評価値ｍ＿ＤＧＩと歩行者の主観評価結果との相関をより高めることができるので、屋外の不快グレアをより正確に評価できる。

但し、上述の実施形態は本発明の一態様であり、本発明の趣旨を逸脱しない範囲において適宜変更可能であるのは勿論である。
例えば、上述の実施形態では、有効グレア輝度の計算領域を所定範囲に指定していたが、輝度分布画像４１の全体から有効グレア輝度を算出してもよい。
また、上述の実施形態では、最大輝度Ｌｍａｘの１／１０以上の輝度値を計算対象としていたが、計算対象はこれに限定されるものではない。例えば、計算対象を最大輝度Ｌｍａｘの１／１０００以上としてもよい。

また、上述の実施形態では、定数αを１．５、定数βを０．５、定数γを０．８としたが、定数α，β，γはこれに限定されるものではない。また、屋外グレア評価装置１１は、定数α，β，γの値を変更可能に構成されてもよい。

また、上述の実施形態では、屋外用照明器具１は、路面Ｓから４．５ｍの高さＴ１に支持されていたが、この高さに限定されるものではなく、例えば、３ｍ〜１２ｍ（８ｍ等）の高さに支持されてもよい。
また、本発明は、路面を照射する防犯灯や街路灯、アプローチライト等の道路灯だけではなく、種々の照明器具に適用可能である。

また、上述の実施形態では、屋外用照明器具を評価対象としていたが、屋外用照明器具の背景が比較的暗い場合には、屋内用照明器具を評価対象としてもよい。

１ 屋外用照明器具
１００ 屋外グレア評価システム
１１ 屋外グレア評価装置
２７ グレア算出部
４１ 輝度分布画像（輝度分布データ）
Ｌｅｆｆ，Ｌｅｆｆ＿１０ 有効グレア輝度
ｍ＿ＤＧＩ 屋外グレア評価値
Ｐ 位置指数
α，β，γ 定数
ω 立体角

Claims (4)

1. 少なくとも１つの屋外用照明器具を含めて撮影した輝度分布データから、屋外用照明器具の輝度分布を順応輝度に応じて重み付け加算した有効グレア輝度Ｌｅｆｆ、立体角ω及び位置指数Ｐを求め、前記有効グレア輝度Ｌｅｆｆ、前記立体角ω、前記位置指数Ｐおよび定数α、β、γに基づいて、屋外グレア評価値ｍ＿ＤＧＩを次式によって算出するグレア算出部を備え、
   前記有効グレア輝度Ｌｅｆｆを、輝度分布データに対して輝度の範囲を指定して求めたことを特徴とする屋外グレア評価装置。
   ｍ＿ＤＧＩ＝ｌｏｇ（Ｌｅｆｆ ^α ×ω ^β ／Ｐ ^γ ）
2. 前記範囲は、輝度分布データにおける最大輝度Ｌｍａｘの１／１０以上としたことを特徴とする請求項１に記載の屋外グレア評価装置。
3. 少なくとも１つの屋外用照明器具を含めて撮影した輝度分布データから、屋外用照明器具の輝度分布を順応輝度に応じて重み付け加算した有効グレア輝度Ｌｅｆｆ、立体角ω及び位置指数Ｐを求め、前記有効グレア輝度Ｌｅｆｆ、前記立体角ω、前記位置指数Ｐおよび定数α、β、γに基づいて、屋外グレア評価値ｍ＿ＤＧＩを次式によって算出し、
   前記有効グレア輝度Ｌｅｆｆを、輝度分布データに対して輝度の範囲を指定して求めることを特徴とする屋外グレア評価方法。
   ｍ＿ＤＧＩ＝ｌｏｇ（Ｌｅｆｆ ^α ×ω ^β ／Ｐ ^γ ）
4. 少なくとも１つの屋外用照明器具を含めて撮影した輝度分布データから、屋外用照明器具の輝度分布を順応輝度に応じて重み付け加算した有効グレア輝度Ｌｅｆｆ、立体角ω及び位置指数Ｐを求め、前記有効グレア輝度Ｌｅｆｆ、前記立体角ω、前記位置指数Ｐおよび定数α、β、γに基づいて、屋外グレア評価値ｍ＿ＤＧＩを次式によって算出し、
   前記有効グレア輝度Ｌｅｆｆを、輝度分布データに対して輝度の範囲を指定して求めることを特徴とする屋外グレア評価プログラム。
   ｍ＿ＤＧＩ＝ｌｏｇ（Ｌｅｆｆ ^α ×ω ^β ／Ｐ ^γ ）