特許 6362744 機能性光源を用いた検査装置と検査方法及び機能性光源とその設計方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6362744

(24)【登録日】2018年7月6日

(45)【発行日】2018年7月25日

(54)【発明の名称】機能性光源を用いた検査装置と検査方法及び機能性光源とその設計方法

(51)【国際特許分類】

   G01N 21/84 20060101AFI20180712BHJP

   G01N 21/27 20060101ALI20180712BHJP

   G01J 3/36 20060101ALI20180712BHJP

   G01N 21/01 20060101ALI20180712BHJP

【ＦＩ】

   G01N21/84 E

   G01N21/27 A

   G01J3/36

   G01N21/01 D

【請求項の数】5

【全頁数】32

(21)【出願番号】特願2017-148493(P2017-148493)

(22)【出願日】2017年7月31日

(62)【分割の表示】特願2013-514001(P2013-514001)の分割

【原出願日】2012年5月2日

(65)【公開番号】特開2017-191116(P2017-191116A)

(43)【公開日】2017年10月19日

【審査請求日】2017年8月28日

(31)【優先権主張番号】特願2011-105474(P2011-105474)

(32)【優先日】2011年5月10日

(33)【優先権主張国】JP

(73)【特許権者】

【識別番号】304027349

【氏名又は名称】国立大学法人豊橋技術科学大学

(72)【発明者】

【氏名】中内 茂樹

(72)【発明者】

【氏名】西野 顕

(72)【発明者】

【氏名】姫野 徹

【審査官】 小野寺 麻美子

(56)【参考文献】

【文献】 特開２００２−４８７３３（ＪＰ，Ａ）

【文献】 特開２００１−１８３３０４（ＪＰ，Ａ）

【文献】 特開２０００−３５６５５３（ＪＰ，Ａ）

【文献】 特表２０１０−５２３９８４（ＪＰ，Ａ）

【文献】 国際公開第２００８／１０５４６０（ＷＯ，Ａ１）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ０１Ｎ ２１／００ − Ｇ０１Ｎ ２１／０１

Ｇ０１Ｎ ２１／１７ − Ｇ０１Ｎ ２１／６１

Ｇ０１Ｎ ２１／８４ − Ｇ０１Ｎ ２１／９５８

Ｇ０１Ｊ ３／００ − Ｇ０１Ｊ ４／０４

Ｇ０１Ｊ ７／００ − Ｇ０１Ｊ ９／０４

ＷＰＩ

ＪＳＴＰｌｕｓ／ＪＭＥＤＰｌｕｓ／ＪＳＴ７５８０（ＪＤｒｅａｍＩＩＩ）

(57)【特許請求の範囲】

【請求項1】

ピーク波長が互いに異なる３つ以上の狭帯域光源を有して前記狭帯域光源が発する光が混合されることにより光源色を白色とする機能性光源の設計方法であって、
２つの弁別対象の複数のサンプルの反射率スペクトルを取得する反射率スペクトル取得ステップと、
ピーク波長が互いに異なる狭帯域光源のスペクトル（以下、「狭帯域スペクトル」という。）を４つ以上用意して、前記狭帯域スペクトルのうちの３つ以上を放射輝度に対する第１の重み付けを行って組み合わせることにより白色となるスペクトル（以下、「白色スペクトル」という。）を、複数設定する光源色スペクトル設定ステップと、
複数の前記白色スペクトルを放射輝度に対する第２の重み付けを行って組み合わせることにより構成される設計光源スペクトルと、前記各サンプルの反射率スペクトルとを用いて、前記設計光源スペクトルを有する設計光源で照明しつつ所定の計測装置により前記各サンプルからの反射光を計測したときの前記各サンプルの計測値を算出し、当該計測値を用いて前記２つの弁別対象のサンプル間の計測値の差を表す評価量を算出する評価量算出ステップと、
を有し、
前記評価量が表す前記計測値の差が略最大となるように前記第２の重み付けを最適化し、
前記第１の重み付けと最適化した前記第２の重み付けとに基づいて、前記機能性光源を設計することを特徴とする設計方法。

【請求項2】

ピーク波長が互いに異なる３つ以上の狭帯域光源を有して前記狭帯域光源が発する光が混合されることにより光源色を白色とする機能性光源の設計方法であって、
２つの弁別対象の複数のサンプルの反射率スペクトルを取得する反射率スペクトル取得ステップと、
ピーク波長が互いに異なる狭帯域光源のスペクトル（以下、「狭帯域スペクトル」という。）を４つ以上用意して、前記狭帯域スペクトルのうちの３つ以上を放射輝度に対する第１の重み付けを行って組み合わせることにより白色となるスペクトル（以下、「白色スペクトル」という。）を、複数設定する光源色スペクトル設定ステップと、
複数の前記白色スペクトルを放射輝度に対する第２の重み付けを行って組み合わせることにより構成される設計光源スペクトルと、前記各サンプルの反射率スペクトルとを用いて、前記設計光源スペクトルを有する設計光源で照明しつつ前記各サンプルを観察したと仮定したときの前記各サンプルの色の見えを表すパラメータを算出し、当該パラメータを用いて前記２つの弁別対象間の色差を表す評価量を算出する評価量算出ステップと、
を有し、
前記評価量が表す色差が略最大となるように前記第２の重み付けを最適化し、
前記第１の重み付けと最適化した前記第２の重み付けとに基づいて、前記機能性光源を設計することを特徴とする設計方法。

【請求項3】

請求項２記載の設計方法であって、
前記評価量算出ステップにおいて、カラーアピアランスモデルに基づいて前記パラメータを算出することを特徴とする設計方法。

【請求項4】

所望の光源色を有する機能性光源の設計方法であって、
弁別の対象となる２つの観察対象についてそれぞれ複数のサンプルの反射率スペクトルを取得する反射率スペクトル取得ステップと、
ピーク波長が互いに異なる狭帯域光源のスペクトル（以下、「狭帯域スペクトル」という。）を４つ以上用意して、前記狭帯域スペクトルのうちの３つ以上を放射輝度に対する第１の重み付けを行って組み合わせることにより前記所望の光源色となるスペクトル（以下、「光源色スペクトル」という。）を、複数設定する光源色スペクトル設定ステップと、
複数の前記光源色スペクトルを放射輝度に対する第２の重み付けを行って組み合わせることにより構成される設計光源スペクトルと、前記各サンプルの反射率スペクトルとを用いて、前記設計光源スペクトルを有する設計光源の下で前記各サンプルを観察したと仮定したときの前記各サンプルの色の見えを表すパラメータを算出し、当該パラメータを用いて前記設計光源下における前記２つの観察対象間の色差を表す評価量を算出する評価量算出ステップと、
を有し、
前記評価量が表す前記色差が略最大となるように前記第２の重み付けを最適化し、
前記第１の重み付けと最適化した前記第２の重み付けとに基づいて、前記機能性光源を設計することを特徴とする設計方法。

【請求項5】

所望の光源色を有する機能性光源の設計方法であって、
２つの弁別対象についてそれぞれ複数のサンプルの反射率スペクトルを取得する反射率スペクトル取得ステップと、
ピーク波長が互いに異なる狭帯域光源のスペクトル（以下、「狭帯域スペクトル」という。）を４つ以上用意して、前記狭帯域スペクトルのうちの３つ以上を放射輝度に対する第１の重み付けを行って組み合わせることにより前記所望の光源色となるスペクトル（以下、「光源色スペクトル」という。）を、複数設定する光源色スペクトル設定ステップと、
複数の前記光源色スペクトルを放射輝度に対する第２の重み付けを行って組み合わせることにより構成される設計光源スペクトルと、前記各サンプルの反射率スペクトルとを用いて、前記設計光源スペクトルを有する設計光源で照明しつつ所定の計測装置により前記各サンプルからの反射光を計測したと仮定したときの前記各サンプルの計測値を算出し、
当該計測値を用いて前記２つの弁別対象のサンプル間の計測値の差を表す評価量を算出する評価量算出ステップと、
を有し、
前記評価量が表す前記計測値の差が略最大となるように前記第２の重み付けを最適化し、
前記第１の重み付けと最適化した前記第２の重み付けとに基づいて、前記機能性光源を設計することを特徴とする設計方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、異物検出、不良品検査、医療診断等における色弁別を補助するための機能性光源を用いた検査装置と検査方法及び機能性光源とその設計方法に関する。

【背景技術】

【0002】

従来の異物検出装置又は異物検出方法として、下記特許文献１〜３に記載されているように、近赤外線や赤外線を用いたものがある。また、色弁別を補助するための多層膜光学フィルタの製造方法として、出願人は下記特許文献４記載の方法を出願している。さらに、下記特許文献５記載の異物検出装置では、撮像カメラを用いて異物検出を行っている。

【先行技術文献】

【特許文献】

【0003】

【特許文献1】特開２００１−９９８７３号公報

【特許文献2】特開２００４−３０１６９０号公報

【特許文献3】特開２００５−８３８８９号公報

【特許文献4】国際公開２００７／０９４３３８号パンフレット

【特許文献5】特開２００８−１５１６１３号公報

【発明の概要】

【発明が解決しようとする課題】

【0004】

従来の近赤外線や赤外線を用いる異物検出装置又は異物検出方法は、コストが高く、導入が容易でないのみならず、目視検査には適用できないという問題があった。また、上記多層膜光学フィルタは、目視検査に応用可能であるが、安価ではなく、やはり導入が容易でないという問題があった。また、目視検査の場合、検査光源の光源色によっては、観察者が光源色に順応することが困難であったり疲労し易くなったりして、目視検査が容易でないという問題があった。また、撮像カメラ等の計測装置を用いて異物検出や不良品検査を行う場合においても、導入が容易な装置を用いるとともに検査精度を向上させたいという要望があった。

【0005】

本発明は、上述した問題を解決するものであり、導入が容易で検査精度を向上させ得る機能性光源を用いた検査装置と検査方法及び機能性光源とその設計方法を提供することを目的とする。

【課題を解決するための手段】

【0006】

本発明の機能性光源を用いた検査装置は、反射光を計測してその計測値を出力可能な計測装置と、ピーク波長が互いに異なる３つ以上の狭帯域光源を選択し、前記狭帯域光源が発する光の発光条件を制御して、それらの光が混合されることにより光源色を白色とする機能性光源であって、前記機能性光源で照明しつつ前記計測装置により２つの弁別対象からの反射光を計測したときの弁別対象間の計測値の差が、略最大となるように構成された前記機能性光源と、前記機能性光源で照明しつつ前記計測装置により検査対象からの反射光を計測したときの前記検査対象の複数の部分の計測値を取得して、当該計測値に基づいて検査結果を判定する判定部と、を備えることを特徴とする。

【0007】

これによれば、白色を中心として色の変化を計測することとなるので、赤、緑及び青の３つの色の計測値を出力するＲＧＢカメラ等の計測装置を用いたときに、それら３つの色の計測値の全てを有効に検査結果の判定に用いることができるとともに、検査対象に含まれる弁別対象間の計測値の差が略最大となるので、検査精度を向上可能である。また、近赤外線、赤外線や多層膜光学フィルタを用いる必要がないので、導入容易である。

【0008】

他の本発明の機能性光源を用いた検査装置は、反射光を計測してその計測値を出力可能な計測装置と、ピーク波長が互いに異なる３つ以上の狭帯域光源を有して前記狭帯域光源が発する光が混合されることにより光源色を白色とする機能性光源であって、前記機能性光源で照明しつつ前記計測装置により２つの弁別対象からの反射光を計測したときの弁別対象間の計測値の差が、前記機能性光源と同等の輝度を有する等エネルギー白色光源で照明しつつ前記計測装置により計測したときよりも拡大するように構成された前記機能性光源と、前記機能性光源で照明しつつ前記計測装置により検査対象からの反射光を計測したときの前記検査対象の複数の部分の計測値を取得して、当該計測値に基づいて検査結果を判定する判定部と、を備えることを特徴とする。

【0009】

これによれば、白色を中心として色の変化を計測することとなるので、赤、緑及び青の３つの色の計測値を出力するＲＧＢカメラ等の計測装置を用いたときに、それら３つの色の計測値の全てを有効に検査結果の判定に用いることができるとともに、同等の輝度を有する等エネルギー白色光源で照明したときよりも弁別対象間の計測値の差が拡大するので、検査精度を向上可能である。また、近赤外線、赤外線や多層膜光学フィルタを用いる必要がないので、導入容易である。

【0010】

ここで、前記計測装置が前記計測値としてＲＧＢ値を出力可能な撮像装置であることが好ましい。白色光源で照明しつつＲＧＢ値を出力可能な撮像装置で計測すれば、計測値の全てを有効に検査結果の判定に用いることができるからである。

【0011】

本発明の機能性光源を用いた検査装置は、前記２つの弁別対象が前記検査対象における異物混入部分と異物無し部分であり、前記判定部が検査結果として異物混入の有無を判定することとすれば、前記検査対象に混入した異物を検出するための検査装置とすることができる。

【0012】

また、本発明の機能性光源を用いた検査装置は、前記検査対象が皮膚であって、前記２つの弁別対象が皮膚における疾患部分と疾患無し部分であり、前記判定部が検査結果として皮膚疾患の診断結果を判定することとすれば、皮膚疾患を診断するための検査装置とすることができる。

【0013】

また、本発明の機能性光源を用いた検査装置は、前記検査対象が塗装面であって、前記２つの弁別対象が塗装むら有り部分と塗装むら無し部分であり、前記判定部が検査結果として塗装むらの有無を判定することとすれば、塗装むらを評価するための検査装置とすることができる。

【0014】

本発明の機能性光源は、ピーク波長が互いに異なる３つ以上の狭帯域光源を有して前記狭帯域光源が発する光が混合されることにより光源色を白色とする機能性光源であって、前記機能性光源で照明しつつ所定の計測装置により２つの弁別対象からの反射光を計測したときの弁別対象間の計測値の差が、前記機能性光源と同等の輝度を有する等エネルギー白色光源で照明しつつ前記計測装置により計測したときよりも拡大するように構成されたことを特徴とする。

【0015】

これによれば、同等の輝度を有する等エネルギー白色光源で照明したときよりも弁別対象間の計測値の差が拡大するので、計測装置による弁別精度を向上可能である。

【0016】

本発明の機能性光源の設計方法は、前記各弁別対象の複数のサンプルの反射率スペクトルを取得する反射率スペクトル取得ステップと、ピーク波長が互いに異なる狭帯域光源のスペクトル（以下、「狭帯域スペクトル」という。）を４つ以上用意して、前記狭帯域スペクトルのうちの３つ以上を放射輝度に対する第１の重み付けを行って組み合わせることにより白色となるスペクトル（以下、「白色スペクトル」という。）を、複数設定する光源色スペクトル設定ステップと、複数の前記白色スペクトルを放射輝度に対する第２の重み付けを行って組み合わせることにより構成される設計光源スペクトルと、前記各サンプルの反射率スペクトルとを用いて、前記設計光源スペクトルを有する設計光源で照明しつつ前記計測装置により前記各サンプルからの反射光を計測したと仮定したときの前記各サンプルの計測値を算出し、当該計測値を用いて前記２つの弁別対象のサンプル間の計測値の差を表す評価量を算出する評価量算出ステップと、を有し、前記評価量が表す前記計測値の差が略最大となるように前記第２の重み付けを最適化し、前記第１の重み付けと最適化した前記第２の重み付けとに基づいて、前記機能性光源を設計することを特徴とする。

【0017】

本発明の検査方法は、ピーク波長が互いに異なる３つ以上の狭帯域光源を有して前記狭帯域光源が発する光が混合されることにより光源色を白色とする機能性光源であって、前記機能性光源で照明しつつ２つの弁別対象を観察したときの弁別対象間の色差が、前記機能性光源と同等の輝度を有する等エネルギー白色光源で照明しつつ観察したときよりも拡大するように構成された前記機能性光源を用い、前記機能性光源で照明しつつ検査対象の複数の部分を観察して前記部分同士の色差に基づいて検査結果を判定することを特徴とする。

【0018】

これによれば、機能性光源により色差が拡大されるとともに、機能性光源の光源色が、通常の環境における光源（例えば、蛍光灯）と同様の白色であるので、観察者が光源色に順応する必要がなくなり、観察者に与える違和感や疲労を軽減できるので、目視検査が容易となって、検査精度が向上し、また、近赤外線、赤外線や多層膜光学フィルタを用いる必要がないので、導入容易である。

【0019】

他の本発明の機能性光源は、ピーク波長が互いに異なる３つ以上の狭帯域光源を有して前記狭帯域光源が発する光が混合されることにより光源色を白色とする機能性光源であって、前記機能性光源で照明しつつ２つの弁別対象を観察したときの弁別対象間の色差が、前記機能性光源と同等の輝度を有する等エネルギー白色光源で照明しつつ観察したときよりも拡大するように構成されたことを特徴とする。

【0020】

これによれば、色差が拡大されるとともに、観察者が光源色に順応する必要がなくなり、観察者に与える違和感や疲労を軽減できるので、目視検査が容易となって、検査精度が向上し、また、近赤外線、赤外線や多層膜光学フィルタを用いる必要がないので、検査への導入が容易である。

【0021】

他の本発明の機能性光源の設計方法は、前記各弁別対象の複数のサンプルの反射率スペクトルを取得する反射率スペクトル取得ステップと、ピーク波長が互いに異なる狭帯域光源のスペクトル（以下、「狭帯域スペクトル」という。）を４つ以上用意して、前記狭帯域スペクトルのうちの３つ以上を放射輝度に対する第１の重み付けを行って組み合わせることにより白色となるスペクトル（以下、「白色スペクトル」という。）を、複数設定する光源色スペクトル設定ステップと、複数の前記白色スペクトルを放射輝度に対する第２の重み付けを行って組み合わせることにより構成される設計光源スペクトルと、前記各サンプルの反射率スペクトルとを用いて、前記設計光源スペクトルを有する設計光源で照明しつつ前記各サンプルを観察したと仮定したときの前記各サンプルの色の見えを表すパラメータを算出し、当該パラメータを用いて前記２つの弁別対象間の色差を表す評価量を算出する評価量算出ステップと、を有し、前記評価量が表す色差が略最大となるように前記第２の重み付けを最適化し、前記第１の重み付けと最適化した前記第２の重み付けとに基づいて、前記機能性光源を設計することを特徴とする。

【0022】

なお、前記評価量算出ステップにおいて、カラーアピアランスモデルに基づいて前記パラメータを算出することができる。

【0023】

また、本発明の機能性光源では、前記狭帯域光源としてＬＥＤ（発光ダイオード）を用いることが好ましい。製造コストを低減でき、導入容易となるからである。

【0024】

更に他の本発明の機能性光源の設計方法は、所望の光源色を有する機能性光源の設計方法であって、弁別の対象となる２つの観察対象についてそれぞれ複数のサンプルの反射率スペクトルを取得する反射率スペクトル取得ステップと、ピーク波長が互いに異なる狭帯域光源のスペクトル（以下、「狭帯域スペクトル」という。）を４つ以上用意して、前記狭帯域スペクトルのうちの３つ以上を放射輝度に対する第１の重み付けを行って組み合わせることにより前記所望の光源色となるスペクトル（以下、「光源色スペクトル」という。）を、複数設定する光源色スペクトル設定ステップと、複数の前記光源色スペクトルを放射輝度に対する第２の重み付けを行って組み合わせることにより構成される設計光源スペクトルと、前記各サンプルの反射率スペクトルとを用いて、前記設計光源スペクトルを有する設計光源の下で前記各サンプルを観察したと仮定したときの前記各サンプルの色の見えを表すパラメータを算出し、当該パラメータを用いて前記設計光源下における前記２つの観察対象間の色差を表す評価量を算出する評価量算出ステップと、を有し、前記評価量が表す前記色差が略最大となるように前記第２の重み付けを最適化し、前記第１の重み付けと最適化した前記第２の重み付けとに基づいて、前記機能性光源を設計することを特徴とする。

【0025】

これによれば、任意の光源色を有するとともに、その光源下で観察したときの観察対象間の色差が拡大される機能性光源を設計できる。

【0026】

また更に他の本発明の機能性光源の設計方法は、所望の光源色を有する機能性光源の設計方法であって、２つの弁別対象についてそれぞれ複数のサンプルの反射率スペクトルを取得する反射率スペクトル取得ステップと、ピーク波長が互いに異なる狭帯域光源のスペクトル（以下、「狭帯域スペクトル」という。）を４つ以上用意して、前記狭帯域スペクトルのうちの３つ以上を放射輝度に対する第１の重み付けを行って組み合わせることにより前記所望の光源色となるスペクトル（以下、「光源色スペクトル」という。）を、複数設定する光源色スペクトル設定ステップと、複数の前記光源色スペクトルを放射輝度に対する第２の重み付けを行って組み合わせることにより構成される設計光源スペクトルと、前記各サンプルの反射率スペクトルとを用いて、前記設計光源スペクトルを有する設計光源で照明しつつ所定の計測装置により前記各サンプルからの反射光を計測したと仮定したときの前記各サンプルの計測値を算出し、当該計測値を用いて前記２つの弁別対象のサンプル間の計測値の差を表す評価量を算出する評価量算出ステップと、を有し、前記評価量が表す前記計測値の差が略最大となるように前記第２の重み付けを最適化し、前記第１の重み付けと最適化した前記第２の重み付けとに基づいて、前記機能性光源を設計することを特徴とする。

【0027】

これによれば、任意の光源色を有するとともに、その光源下で計測装置により計測したときの弁別対象間の色差が拡大される機能性光源を設計できる。

【発明の効果】

【0028】

本発明は、検査への導入が容易で検査精度を向上させ得る。

【図面の簡単な説明】

【0029】

【図1】（ａ）は第１実施形態に係る照明装置の概略斜視図、（ｂ）は同照明装置の主要部の一部破断斜視図である。

【図2】CIECAM02における観察者の視野の定義を示す図である。

【図3】第１実施形態に係る設計方法のフローチャートである。

【図4】（ａ）は実施例で用いたサンプルを入れたシャーレの斜視図、（ｂ）は同シャーレの平面図である。

【図5】各サンプルの反射率スペクトルを示すグラフである。

【図6】モデル化されたＬＥＤの放射輝度特性を示すグラフである。

【図7】実施例の設計光源の放射輝度スペクトルを示すグラフである。

【図8】実施例の設計光源及び等エネルギー白色光による平均色差を示すグラフである。

【図9】（ａ）は等エネルギー白色光源を用いた場合、（ｂ）は設計光源を用いた場合のマルチスペクトルイメージによるシミュレーションの結果を示す画像である。

【図10】等エネルギー白色光を照射した場合の画像である。

【図11】多波長可変光源による設計光を照射した場合の画像である。

【図12】実施例の一対比較実験の結果を示すグラフである。

【図13】第２実施形態の検査装置の概略構成図である。

【図14】第２実施形態の画像処理の流れを示す図である。

【図15】第２実施形態のデータ処理の流れを示すフローチャートである。

【図16】第３実施形態の検査装置の概略構成図である。

【図17】第４実施形態の検査装置の概略構成図である。

【発明を実施するための形態】

【0030】

〈第１実施形態〉
第１実施形態に係る照明装置１は、図１に示すように、機能性光源４と図示しない電源部とを備え、機能性光源４は、狭帯域光源２であるＬＥＤが片面に多数配置された基板３と、拡散板５とから構成されている。基板３上に配置されるＬＥＤの種類、配置個数、及び、各ＬＥＤに対する印加電流の決定方法については、後述する。拡散板５は、基板３の前方に配置されている。なお、基板３では、ＬＥＤが配置されている側を前とする。電源部は、基板３及び外部電源に接続されている。基板３上のＬＥＤは、スイッチ６を入れることにより、電源部を介して外部電源から電気の供給を受けて発光し、全ＬＥＤから放射された光は、拡散板５で完全に混合されて、拡散板５の前方に放射される。

【0031】

実施形態に係る機能性光源４の設計方法（すなわち、基板３に配置するＬＥＤの種類、個数、及び、各ＬＥＤに対する印加電流の決定方法）について、次に説明する。

【0032】

機能性光源４の設計方法は、図２に示すように、弁別の対象となる２つの観察対象（弁別対象）についてそれぞれ複数のサンプルの反射率スペクトルを取得する反射率スペクトル取得ステップ（ステップＳ０１）と、背景の反射率スペクトル及び順応領域の輝度を設定しておくステップ（ステップＳ０２）と、ピーク波長の異なるＬＥＤのスペクトルを４つ以上用意して、それらのスペクトルのうちの３つを放射輝度に対する重み付けを行って組み合わせることにより所望の光源色となる光源色スペクトルを、複数設定する光源色スペクトル設定ステップ（ステップＳ０３）と、複数の光源色スペクトルを放射輝度に対する重み付けを行って組み合わせることにより構成される設計光源スペクトルと、各サンプルの反射率スペクトルとを用いて、設計光源スペクトルを有する設計光源の下で各サンプルを観察したと仮定したときの各サンプルの色の見えを表すパラメータを算出し、それらのパラメータを用いて、設計光源下における２つの観察対象間の色差を表す評価量を算出する評価量算出ステップ（ステップＳ０４〜０６）と、を有し、評価量が略最大となるように光源色スペクトルを組み合わせるときの重み付けを最適化する。すなわち、図２のステップＳ０７において評価量が略最大となる重み付けが得られたと判断するまで、光源色スペクトルを組み合わせるときの重み付けを変えながら評価量算出ステップを繰り返す。なお、ステップＳ０１〜０３の順序は問わない。そして、最適化した重み付けで光源色スペクトルを組み合わせることにより構成される設計光源スペクトルに基づいて、機能性光源を設計する。具体的には、設計光源スペクトルと略同等のスペクトルを、機能性光源のスペクトルとする。なお、「略最大」とは最大である場合と最大に近い場合を含み、「略同等」とは同等である場合と近似している場合を含む。

【0033】

また、所望の光源色は、任意の光源色を選択可能であるが、実施形態では白色とする。したがって、以下、「光源色スペクトル」という代わりに「白色スペクトル」という場合がある。なお、「白色」には、蛍光灯の色の分類に用いられる「電球色」「温白色」「白色」「昼白色」「昼光色」のいずれの色も含まれるものとする。

【0034】

以下、実施形態に係る設計方法について詳細に説明する。なお、弁別の対象とする２つの観察対象を、第１観察対象及び第２観察対象と表記する。

【0035】

〈反射率スペクトル取得ステップ〉
第１観察対象と第２観察対象について、それぞれN個（Nは１以上の整数）のサンプルを用意して、分光反射率を計測可能な計測装置で各サンプルを計測することにより、第１観察対象の各サンプルの反射率スペクトルR_1j(λ)と、第２観察対象の各サンプルの反射率スペクトルR_2j(λ)を取得する（但し、j=1,2,…,N）。λは波長であり、可視光域である420〜700nmの範囲とする（図２のステップＳ０１）。

【0036】

また、観察対象を観察するときの背景（Background）の反射率スペクトルR_b(λ)、及び、観察対象を観察するときの順応領域（Adapting field）の輝度L_Aを、観察環境に応じて設定する（図２のステップＳ０２）。

【0037】

〈光源色スペクトル設定ステップ〉
ピーク波長が互いに異なるＬＥＤをM個（Mは４以上の整数）用意する。ピーク波長は、380〜780nmの範囲内にあるものとする。用意したＬＥＤの分光放射輝度スペクトルを、I₁(λ),I₂(λ),I₃(λ),…,I_M(λ)とする。そして、これらのスペクトルのうちの３つを放射輝度に対する重み付けを行って組み合わせることにより、所望の光源色である白色となるスペクトルを、すべて求める（図２のステップＳ０３）。

【0038】

この求め方について、以下、詳説する。スペクトルI₁(λ),I₂(λ),I₃(λ),…,I_M(λ)から任意に３つ選択したスペクトルを、ピーク波長が長い順にI_l(λ)、I_m(λ)、I_s(λ)と置くと、これらのスペクトルを放射輝度に対する重み付けを行って組み合わせることによって得られる、所望の光源色（Ｘ_ｗ，Ｙ_ｗ，Ｚ_ｗ）と輝度L_Wを持つ白色光スペクトルＷ（λ）は、次式[数１]で表せる。なお、a_l,a_m,a_sは、重みである。

【0039】

【数1】

【0040】

ここで、所望の白色光の光源色（Ｘ_ｗ，Ｙ_ｗ，Ｚ_ｗ）ならびに組み合わせる各ＬＥＤの光源色（Ｘ_ｌ，Ｙ_ｌ，Ｚ_ｌ），（Ｘ_ｍ，Ｙ_ｍ，Ｚ_ｍ），（Ｘ_ｓ，Ｙ_ｓ，Ｚ_ｓ）を、次式[数２]、[数３]で定義する。W(λ)は、等エネルギー白色光のスペクトルを用いる。関数f_cは、光源のスペクトルI(λ)を入力として等色関数から光源色の三刺激値ＸＹＺを求める関数である。また、[数３]におけるL，L_Wの単位はcd/m²であり、V(λ)は既知の分光視感度特性を表す。

【0041】

【数2】

【0042】

【数3】

【0043】

上記式[数２]の各光源色より、次式[数４]及び[数５]を用いて、所望の白色光を得るための各ＬＥＤのスペクトルの配合強度、すなわち、重みa_l,a_m,a_sを決定する。

【0044】

【数4】

【0045】

【数5】

【0046】

上記配合強度^T[a_la_ma_s]が一意に定まり（すなわち、行列Ａのランクが３であり）、更にa_l,a_m,a_sがすべて正の値をとるときのみ、選択した３つのＬＥＤによる所望の白色光出力が可能である。ゆえに、ＬＥＤのスペクトルの総数Mに対して、３つを取り出して作る組み合わせ_MC₃から、行列Ａのランクが３であり、かつ、重みa_l,a_m,a_sがすべて正の値をとるものを、すべて抽出することにより、総数Lの白色スペクトルW₁(λ),W₂(λ),W₃(λ),…,W_L(λ)が求められる。

【0047】

〈評価量算出ステップ〉
上記白色スペクトルW₁(λ),W₂(λ),W₃(λ),…,W_L(λ)を、放射輝度に対する重み付けを行って組み合わせることにより、設計光源スペクトルを構成する。W₁(λ),W₂(λ),W₃(λ),…,W_L(λ)の放射輝度に対する重みを、それぞれ、b₁,b₂,b₃,…,b_Lとすると、設計光源スペクトルW_F(λ)は、次式[数６]で表される。なお、b₁,b₂,b₃,…,b_Lはそれぞれ０以上であり、b₁＋b₂＋b₃＋…＋b_L＝１とする。

【0048】

【数6】

【0049】

設計光源スペクトルW_F(λ)、第１観察対象及び第２観察対象の各サンプルの反射率スペクトルR_1j(λ)，R_2j(λ)、背景の反射率スペクトルR_b(λ)、及び、順応領域の輝度L_Aを用いて、カラーアピアランスモデル(Color Appearance Model)CIECAM02に基づいて、スペクトルW_F(λ)を有する設計光源の下で第１観察対象及び第２観察対象の各サンプルを観察したと仮定したときの各サンプルの色の見えを表すパラメータを算出する（図２のステップＳ０５）。

【0050】

ここで、CIECAM02による各サンプルの色の見えを表すパラメータの算出方法について説明するが、この算出方法は下記参考文献２に詳しく記載されているので、詳説はしない。CIECAM02では、図３及び表１に示すように、注視点から２度視野内を刺激（Stimulus）２０、刺激２０の外側の１０度視野内を背景（Background）２１、背景２１の外側の視野を周辺領域（Surround field）２２、背景２１と周辺領域２２を合わせた領域を順応領域（Adapting field）と定義する。

【0051】

設計光源下での各サンプルの色の見えをCIECAM02で計算する場合、各サンプルの三刺激値（Ｘ_ｉｊ，Ｙ_ｉｊ，Ｚ_ｉｊ）、設計光源の三刺激値（Ｘ_ｗ，Ｙ_ｗ，Ｚ_ｗ）、順応領域の輝度L_A[cd/m²]、設計光源下で白色を観測したときの相対輝度Ｙ_ｗ、及び、背景の相対輝度Ｙ_ｂが必要である。Ｙ_ｗは上記式[数２]で求めたものを用いる。また、b₁＋b₂＋b₃＋…＋b_L＝１であることから、上記式［数２］で求めたＸ_ｗ，Ｙ_ｗ，Ｚ_ｗを設計光源の三刺激値とすることができる。L_Aは、上記ステップＳ０２で設定したものを用いる。Ｙ_ｂは、設計光源スペクトルW_F(λ)、及び、上記ステップＳ０２で設定した背景の反射率スペクトルR_b(λ)から、次式[数７]により算出する。

【0052】

【数7】

【0053】

また、各サンプルの三刺激値Ｘ_ｉｊ，Ｙ_ｉｊ，Ｚ_ｉｊを、次式［数８］により算出する。

【0054】

【数8】

【0055】

観察する周囲環境条件によるパラメータF，c，Ncを、表２から選択する。観察環境に応じて周辺視野を３段階(Average, Dim, Dark)から選択することで、より正確に見えを評価可能である。実施形態では、平均（Average）とする。また、次式［数９］に従って、パラメータk，F_L，n，N_bb，zを計算する。

【0056】

【数9】

【0057】

そして、次式［数１０］に従って、各サンプルの三刺激値Ｘ_ｉｊ，Ｙ_ｉｊ，Ｚ_ｉｊをＲ，Ｇ，Ｂに変換する。設計光源の三刺激値Ｘ_ｗ，Ｙ_ｗ，Ｚ_ｗについても同様にＲ_ｗ，Ｇ_ｗ，Ｂ_ｗに変換する。

【0058】

【数10】

【0059】

次式［数１１］に従って順応の度合いＤを計算する。Ｄは、０から１の範囲値で示され、選択された白色点（Ｘ_ｗ，Ｙ_ｗ，Ｚ_ｗ）に対し、無順応の場合には０、完全順応の場合には１となる。

【0060】

【数11】

【0061】

計算したＤから、各サンプルのＲ，Ｇ，Ｂを、次式［数１２］に従って、Ｒ_ｃ，Ｇ_ｃ，Ｂ_ｃに変換する。なお、［数１２］はＲについての式であるが、Ｇ，Ｂについても同様である。Ｒ_ｗ，Ｇ_ｗ，Ｂ_ｗも同様にＲ_ｃｗ，Ｇ_ｃｗ，Ｂ_ｃｗに変換する。

【0062】

【数12】

【0063】

次式［数１３］に従って、Ｒ_ｃ，Ｇ_ｃ，Ｂ_ｃをＲ'，Ｇ'，Ｂ'に変換する。

【0064】

【数13】

【0065】

次式［数１４］に従って、Ｒ'，Ｇ'，Ｂ'をＲ_ａ'，Ｇ_ａ'，Ｂ_ａ'に変換する。

【0066】

【数14】

【0067】

そして、これらから、予備の座標ａ，ｂ及びｔを、次式［数１５］に従って計算する。

【0068】

【数15】

【0069】

色相ｈを、次式［数１６］により計算する。

【0070】

【数16】

【0071】

計算したｈを、表３の各色（Red、Yellow、Green、Blue）のｈと比較し、最も近い色の表３におけるｅ，ｈ，Ｈを、ｅ_１，ｈ_１，Ｈ_１とし、最も離れている色の表３におけるｅ，ｈを、ｅ_２，ｈ_２として、次式［数１７］に従って、ｅ，Ｈを計算する。

【0072】

【数17】

【0073】

そして、次式［数１８］に従って、相対的な明るさ（Lightness）Ｊ、及び、絶対的な明るさ（Brightness）Ｑを計算する。

【0074】

【数18】

【0075】

また、次式［数１９］に従って、相対的な彩度（Chroma）Ｃ、及び、絶対的な彩度(Colorfulness)Ｍを計算する。

【0076】

【数19】

【0077】

このようにして、サンプルの三刺激値Ｘ_ｉｊ，Ｙ_ｉｊ，Ｚ_ｉｊ、背景の三刺激値Ｘ_ｂ，Ｙ_ｂ，Ｚ_ｂ、設計光源の三刺激値Ｘ_ｗ，Ｙ_ｗ，Ｚ_ｗ、及び、順応領域の輝度L_Aから、そのサンプルの色の見えを表すパラメータＪ_ｉｊ，Ｃ_ｉｊ，ｈ_ｉｊ，Ｑ_ｉｊ，Ｍ_ｉｊが得られる（［数２０］参照）。

【0078】

【数20】

【0079】

次に、各サンプルの色の見えを表すパラメータから、第１観察対象及び第２観察対象の対応するサンプル間の色差ΔＥ'を、CIECAM02に基づいて計算する。対応するサンプルとは、同時に計測したサンプルであり、番号ｊが同じである。以下、CIECAM02に基づく色差ΔＥ'の算出方法について説明するが、この算出方法については参考文献３に詳しく記載されているので、ここでは詳説しない。

【0080】

CIECAM02に基づく色差ΔＥ'の算出方法では、パラメータＪ_ｉｊ，ｈ_ｉｊ，Ｍ_ｉｊを用いる。また、求める色差の程度によって、表４からパラメータＫ_Ｌ，ｃ_１，ｃ_２を選択して用いる。CAM02-LCDは大きな色差を求める場合、CAM02-SCDは小さな色差を求める場合、CAM02-UCSは、それらの中間程度の色差を求める場合に用いるバージョンであり、実施形態では、バージョンCAM02-SCDを用いる。

【0081】

そして、次式［数２１］に従って、第１観察対象のサンプルのパラメータＪ_１ｊ，ｈ_１ｊ，Ｍ_１ｊから第１観察対象のサンプルについてＪ'，ａ'，ｂ'を計算し、第２観察対象のサンプルのパラメータＪ_２ｊ，ｈ_２ｊ，Ｍ_２ｊから第２観察対象のサンプルについてＪ'，ａ'，ｂ'を計算して、対応するサンプル間のＪ'，ａ'，ｂ'の差分ΔＪ’，Δａ'，Δｂ'に基づいて、色差ΔＥ'を、計算する。

【0082】

【数21】

【0083】

色差ΔＥ'を全てのサンプル対について求め、その平均値を評価量とする。すなわち、第１観察対象及び第２観察対象のｊ番目のサンプル間の色差をΔＥ_ｊ'とすると、色差ΔＥ_ｊ'は次式［数２２］で表せ、評価量は次式［数２３］で計算できる（図２のステップＳ０６）。なお、［数２２］のsrdは、CAM02-LCD、CAM02-UCS、又は、CAM02-SCDから選択したパラメータＫ_Ｌ，ｃ_１，ｃ_２を示す。［数２３］で計算される評価量は、第１観察対象と第２観察対象との間の色差を表している。

【0084】

【数22】

【0085】

【数23】

【0086】

そして、重み付けb₁,b₂,b₃,…,b_Lを変えながら上記評価量算出ステップを繰り返すことにより、上記評価量が最大となるように重みb₁,b₂,b₃,…,b_Lを最適化する。重み０が付された白色スペクトルは使用されないこととなるため、ＬＥＤの組み合わせと配合強度を同時に最適化することとなる。実時間内で優良解に到達する最適化法としては、山登り法（Hill Climbing法）、焼きなまし法（Simulated Annealing法）、遺伝的アルゴリズム、進化的アルゴリズム等がある。

【0087】

そして、最適化された重みb₁,b₂,b₃,…,b_Lで白色スペクトルを組み合わせることにより構成される設計光源スペクトルと同等のスペクトルを、機能性光源４のスペクトルとする。そのために、機能性光源４では、３つのＬＥＤを組み合わせて白色スペクトルを作るときの重みa_l,a_m,a_s、及び、白色スペクトルを組み合わせて設計光源スペクトルを作るときの最適化した重みb₁,b₂,b₃,…,b_Lに基づいて、基板３上に配置されるＬＥＤの種類（すなわち、どのピーク波長のＬＥＤを用いるか）、個数、及び、各ＬＥＤに対する印加電流を決定する。

【0088】

例えば、上述したように用意したM個（すなわち、M種類）のＬＥＤのうち、ピーク波長λｐのＬＥＤが白色スペクトルW₁,W₂,W₃を作るときにだけ採用されたとし、白色スペクトルW₁,W₂,W₃を作るときにピーク波長λｐのＬＥＤに付された重みをそれぞれa₁,a₂,a₃とする。そして、設計光源スペクトルを構成するときに白色スペクトルW₁,W₂,W₃に付された重みをb₁,b₂,b₃とすると、ピーク波長λ_ｐのＬＥＤに対して付された重みは、全体でa₁b₁＋a₂b₂＋a₃b₃となる。したがって、この重みになるように、ピーク波長λ_ｐのＬＥＤの基板３上への配置個数及び配置されたピーク波長λ_ｐのＬＥＤへの印加電流を決定する。なお、どの白色スペクトルを作るときにも採用されなかったＬＥＤ、及び、重み０が付された白色スペクトルにのみ使用されているＬＥＤは、基板３上に実装しない。

【0089】

この機能性光源４は、観察対象間の色差を、設定したＬＥＤを組み合わせて構成される同等の輝度を有する白色光源の中で、略最大とするものであり、観察対象間の色差を、同等の輝度を有する等エネルギー白色光源下よりも、拡大する。

【0090】

実施形態の設計方法によれば、光源色が白色であって、その光源下で観察したときの観察対象間の色差が、同等の輝度を有する等エネルギー白色光源下で観察したときよりも拡大される機能性光源４を設計できる。光源色が通常の環境における光源（例えば、蛍光灯）と同様の白色であるので、観察者が光源色に順応する必要がなくなり、観察者に与える違和感や疲労を軽減でき、観察対象間の色差が拡大されることと相まって、目視検査を容易とし、検査精度を向上させ得る光源を設計可能である。また、フィルタを通さないので、出力された光の強度が落ちることがなく、高効率の光源を設計可能である。

【0091】

機能性光源４は、目視検査に好適であり、機能性光源４で照明しつつ検査対象の複数の部分を目視で観察してその部分同士の色差に基づいて検査結果を判定する検査方法に用いることができる。この方法によれば、機能性光源４により色差が拡大されるとともに、機能性光源４１の光源色が白色であるので、目視検査が容易となって、検査精度が向上し、また、近赤外線、赤外線や多層膜光学フィルタを用いる必要がないので、検査への導入が容易である。

【0092】

また、狭帯域光源２としては、半値幅１０〜２０ｎｍ程度のものが機能性光源４の設計に適しており（このことは、後述する第２〜４実施形態でも同様である。）、また、必ずしもＬＥＤを用いる必要はなく、例えば、有機ＥＬ（Electronic Luminescence）を用いてもよいが、ＬＥＤは安価であるため、狭帯域光源２としてＬＥＤを用いることにより、設計光源スペクトルと同等のスペクトルを有する機能性光源４を、安価に作成できる。

【0093】

また、カラーアピアランスモデルに基づいて色の見えを表すパラメータを算出すれば、カラーアピアランスは観察環境から色の見えを予測できるので、実際の光源使用環境において有効な光源を設計することができる。

【0094】

なお、機能性光源４の光源色は白色としたが、上述した設計方法によれば、任意の色を光源色とする機能性光源を設計可能であり、光源色を周囲の環境に応じた色とすれば、観察対象間の色差が拡大されることと相まって、検査精度が向上させることができる。

【実施例】

【0095】

本設計方法の実施例として、ブルーベリージャムに混入した異物の視認性を向上させる白色光源の設計を行った。表５に実施例の条件設定を示す。異物が混入している異物混入部が第１観察対象、異物が混入していない異物無し部が第２観察対象に相当する。以下、詳細に説明する。

【0096】

１．設計方法
次の(1)〜(5)に示す手順で設計した。

【0097】

(1)観察環境（観察対象及び背景）のスペクトルの設定
弁別対象である異物混入部ならびに異物無し部の反射率スペクトルは、以下に示す環境の下で実測し、取得した。

【0098】

シャーレにブルーベリージャムを入れ、異物としてブルーベリーの枝１０ならびに葉１１、そして黒色の石１２を、図４（ａ）に示すようにジャム中に混入しサンプルを作成し、サンプルの上方から反射率スペクトルの計測を行った。図４（ｂ）のように、枝１０、葉１１、及び、石１２は、平面視においてシャーレを４分割したうちの３つの領域Ｐ１、Ｐ２、Ｐ３に、それぞれ混入されている。そして、異物はシャーレ底からの厚さ5.25mmのブルーベリージャム上に配置されており、さらに異物を被覆するジャムを1.75mm、3.5mm、5.25mm、7mmの４段階に変化させた。使用したジャムは一般的市販ブルーベリージャム三種である。計測装置は、ケンブリッジ・リサーチ・アンド・インストゥルメンテーション（Cambridge Research & Instrumentation）社製マルチスペクトルイメージング装置Nuance（VIS）、光源にはキノフロ(Kino Flo)社製大型蛍光照明Diva-Lite、白色校正にラブスフェア（Labsphere）社製Spectralon反射板を使用した。計測された分光画像から、異物が混入している部分の反射率スペクトルをそれぞれ抽出し、異物混入部のスペクトル群（第１群）R_1j(l)：j={1,2,…,N}とし、対応する異物無し部のスペクトル群（第２群）R_2j(l)：j={1,2,…,N}には同じシャーレ中の異物が混入されていない領域Ｑの反射率スペクトルを割り当てた。つまり、異物無し部のスペクトル群には、同じ反射率スペクトルを葉１１、枝１０、石１２に対応させて三度使用している。シャーレはジャムの各種類につき４個用意したので、N＝3（異物の種類）×3（ジャムの種類）×4（ジャムの深さ）×4（シャーレの個数）＝144である。計測された被覆ジャム層5.25mm時の反射率スペクトルを図５に示す。図５では、縦軸は反射率、横軸は波長であり、異物無し部をBlueberry Jam、葉１１、枝１０、石１２が混入された異物混入部をそれぞれLeaf,Twig,Stoneと表記している。

【0099】

図５に示される長波長帯域の急峻な反射特性は、レッドエッジと呼ばれる植物特有の反射特性である。そのため、葉１１や枝１０が混入した場合にはこの反射強度が高い。故に、この波長帯に強いエネルギーを持つ光源を用いれば異物の視認性が向上する可能性が高い。

【0100】

この物体を観察する環境として、配置されるジャムの周囲（背景輝度）は波長特性が一様の２０％グレーを仮定した。但し、背景輝度は順応作用や明度対比効果に影響するため、異物検査環境に合わせて修正されることが望ましい。

【0101】

(2)ＬＥＤのスペクトル設定と組み合わせ構築
実施例では、狭帯域光源としてＬＥＤを仮定しているが、ここでは実測したＬＥＤの放射輝度ではなく、図６に示されるモデル化された分光特性を使用した。これは後に使用する多波長可変光源の特性を基にしている。図６の縦軸は放射輝度、横軸は波長である。

【0102】

ＬＥＤの分光特性のモデリングは、次式［数２４］に示すスプライン関数を用いた。詳しくは、中心波長λ_kを420〜700nmの範囲で5nm刻み、スプライン関数の幅を表すパラメータωを10nmとし、これを［数２５］に示すように両隣合わせて3つ重ねたものを使用した。これを言い換えると，ピーク波長が420〜700nmの波長範囲で15nmずつ異なり、各々およそ17.5nmの半値幅の分光放射輝度を持つＬＥＤが存在するという仮定である。なお、中心波長420nmと700nmをのぞくすべてのＬＥＤのピーク波長強度は最大0.37W/sr/m²/nmとし、両端の二つはそれぞれ0.279W/sr/m²/nm，0.152W/sr/m²/nmとした。これは後に出力光をシミュレートするために使用する多波長可変光源の分光放射輝度を元にした設定である。

【0103】

【数24】

【0104】

【数25】

【0105】

上記のように仮定したＬＥＤの分光放射輝度を用い、白色光を構成する３つのＬＥＤの組み合わせを、全て求めた。

【0106】

(3)観察対象ならびに背景の色の見えの計算
上記(2)で求めた３つのＬＥＤの組み合わせで構築される白色光の分光放射輝度、ならびに、観察対象の反射率スペクトルから、観察対象ならびに背景の色の見えを表すパラメータを求めた。算出方法は上記評価量算出ステップで示したとおりであり、設計光源の輝度が一定値となるように配合強度b_kの正規化を行っている。設計光源の輝度が一定であるため、背景R_b(l)の相対輝度Y_bもまた一定である。順応領域の輝度L_Aは全ＬＥＤフル発光時の輝度を基準にその20%に設定した。設計光源の輝度は、このフル発光した場合の輝度を基準としてその割合で設定した。ＬＥＤの最大出力に制限があるため、輝度が低い方が最適化の自由度は高く、また装置全体の消費エネルギーも抑えることができる。しかしながら、輝度が低い場合は順応が進まず色弁別能力が低下する。そこで本実施例では10%,20%,30%,50%,70%の5段階について設計を行い、それぞれの色差を評価した。

【0107】

(4)色弁別性の計算
上記(3)で得られた色の見えを表すパラメータから、CIECAM02を用いて色差を求めた。算出方法は上記評価量算出ステップで示したとおりである。

【0108】

(5)ＬＥＤの組み合わせ修正（最適化）
評価量すなわち異物混入部と異物無し部との平均色差が最大化されるように、配合強度b_kを最適化した。最適化の方法として山登り法を採用した。この手法の特徴は、ある程度短い時間で局所解に到達することができ、到達する解は初期値から一意に定まるということである。そのため、初期値を変化させながら何度か繰り返すことで良好な解を得ることができる。

【0109】

２．設計結果
以下に、実施例の結果を示す。

【0110】

(1)設計結果
輝度を10%,20%,30%,50%,70%とした場合のそれぞれについて、得られた設計光源スペクトルW_F(λ)を図７に示す。図７の縦軸は放射輝度、横軸は波長である。また、輝度を10%,20%,30%,50%,70%とした場合のそれぞれについて、等エネルギー白色光源下での平均色差と、設計光源下での平均色差とを、図８に示す。図８では、縦軸は平均色差、横軸は輝度であり、黒丸●は設計光源、黒四角■は等エネルギー白色光源を示し、いずれも左から10%,20%,30%,50%,70%の場合を示す。

【0111】

図８より、設計光源を使用する場合には、同等の輝度をもつ等エネルギー白色光源を使用する場合と比較して、色差が向上していることが分かる。さらに、等エネルギー白色光源を使用する場合は輝度が小さいほど弁別能力が低下するが、対して設計光源は輝度が低い条件において、色差が向上していることが分かる。さらに輝度を低く設定すれば色差が低下すると考えられるが、今回設定した輝度条件に関して言えば低輝度条件において良好な光源が設計可能であったことが分かる。

【0112】

(2)分光画像ならびに多波長可変光源によるシミュレーション
得られた設計光源スペクトルを用いて、光源照射時のブルーベリージャムの見えをシミュレートした。シミュレートの方法として、(i)マルチスペクトルイメージから計算的に求める方法と、(ii)多波長可変光源を用いて設計光源スペクトルと同等のスペクトルを有する光を生成し、実際のサンプルに照射する方法とを用いた。また、(i)では、マルチスペクトルイメージとして、１．(1)で計測したスペクトル収集用の画像を使用した。(ii)では、多波長可変光源として、株式会社ニコン製ELS-VISを使用し、この多波長可変光源で、設計光源スペクトルと同等のスペクトルを有する出力光を生成し、シャーレ中のジャムに石１２を混入した評価用のサンプルに照射して、RGBカメラで撮影した。

【0113】

図９は、(i)のマルチスペクトルイメージを用いてシミュレートした画像であり、(a)は等エネルギー白色光源下での見え、(b)は設計光源下での見えを示す。両光源の輝度は同じであるが、明らかに、設計光源下では等エネルギー白色光源下よりも、ジャム中の異物が視認しやすくなっていることが分かる。

【0114】

また、図１１に、(ii)の多波長可変光源を用いて生成した光を照射した結果を示し、図１０に、等エネルギー白色光を照射した結果を示す。両光源の輝度は同じであるが、図１１のほうが明るく映っている。これは、ジャムに吸収されずに反射してくる光が多いためである。結果として、ジャム中の石１２の視認性が大きく向上していることが分かる。なお、図９，図１０，図１１は、いずれも元の画像はカラーであるが、白黒に変換して示している。

【0115】

また、図１０，図１１においては、各サンプルの右横に、白色から濃い灰色までの５段階の無彩色を示す同じグレーチャートが配置され、さらに、グレーチャートの右横に、青、緑、赤、黄色、桃色の５色を示す同じカラーチャートが配置されている。グレーチャートの見えは図１０，図１１でほとんど差がないことから、図１１の多波長可変光源からの光は、図１０の等エネルギー白色光と略同じ白色であり、観察者に疲労感や違和感を与え難いことが分かる。

【0116】

(3)心理物理実験による光源評価
多波長可変光源により設計光源スペクトルと同等のスペクトルの光を上記５種の輝度条件のそれぞれについて生成し、すなわち、多波長可変光源により５種の機能性光源を構成し、これら５種の機能性光源と、同等の輝度を持つ等エネルギー白色光源５種を用いて、異物の検出しやすさを問う主観評価実験を行った。実験にはサーストンの一対比較実験を採用した。具体的には、以下のパラダイムを全組み合わせについて行った。

【0117】

(i)上記計１０種の光源から２種の光源を選択
(ii)選択された光源を5秒ずつ（間に消灯1.5sec）照射
(iii)どちらが見分けやすかったかを二肢強制選択
なお、実験には、枝１０、葉１１、及び石１２を同じシャーレのジャム中に混入した評価用のサンプルを使用した。結果は選択確率から求めたZスコアで評価するものとした。図１２は、本評価実験の結果を示すものであり、縦軸は尺度値（Zスコア）、横軸は輝度であり、黒丸●は機能性光源、黒四角■は等エネルギー白色光源を示す。

【0118】

図１２より、低輝度の機能性光源が最も高性能であり、さらにいずれの機能性光源も等エネルギー白色光源の尺度値を上回ることが分かる。これは図８に示した平均色差とも一致する。このことから、CIECAM02を元に求まる色差は、主観的な色弁別性を良好に数値化できていることが確認できる。以上の結果より、本設計方法による機能性光源の設計は、観察者の熟練を要する目視検査のサポートに非常に有効であるといえる。

【0119】

なお、実施例では、多波長可変光源により機能性光源を構成したが、例えば、株式会社東京インスツルメンツ製の波長プログラマブル光源も、任意のスペクトルの光を生成可能であり、機能性光源として使用可能である。このように、本発明の機能性光源は、ＬＥＤ等の狭帯域光源を用いなくても構成可能である。

【0120】

なお、光源色スペクトルを設定するときは、少なくとも３つの狭帯域光源のスペクトルを組み合わせることとする。２つ以下の狭帯域光源のスペクトルの組み合わせでは、光源色スペクトルの光源色を任意に設定できないからである。

【0121】

また、評価量を最適化するとき、評価量が最大となるように最適化することが望ましいが、最大に近くなるように最適化すればよい。すなわち、略最大とは、最大である場合と、最大に近い場合を含む意である。同様に、機能性光源のスペクトルは、設計上のスペクトルである設計光源スペクトルと同じとすることが望ましいが、少なくとも設計光源スペクトルに近似していればよい。すなわち、略同等とは、同じである場合と、近似している場合を含む意である。実用上問題にならない程度に色差が拡大すればよいからである。

【0122】

なお、評価量は、第１観察対象と第２観察対象との間の色差を示す値であればよく、［数２３］で計算される平均色差に限られない。

【0123】

〈第２実施形態〉
次に、図１３、１４、１５に基づいて、第２実施形態に係る検査装置３０について説明する。検査装置３０は、ブルーベリージャムの異物検出装置として構成されており、２つの機能性光源３１と、計測装置に相当するＲＧＢカメラ３２と、ＲＧＢカメラ３２に接続されパーソナル・コンピュータからなる判定部３３と、判定部３３に接続された排除機構制御ユニット３４とを備えている。ＲＧＢカメラ３２は、ＲＧＢ値を出力可能な撮像装置である。判定部３３は、画像処理部３３１とデータ処理部３３２とを備えている。

【0124】

各機能性光源３１は、ピーク波長の異なる３つ以上のＬＥＤが配列された基板３１１と、ＬＥＤからの光を混合する拡散板（図示せず。）とを備え、外部電源から電気の供給を受けて白色光を発光する。

【0125】

各機能性光源３１は、弁別の対象をブルーベリージャムの異物混入部と異物無し部として、第１実施形態で説明した設計方法により設計されたものである。但し、第１実施形態では、評価量算出ステップにおいて、各サンプルの色の見えを表すパラメータをCIECAM02により算出したが、第２実施形態では、各サンプルに対しＲＧＢカメラ３２から出力されると予想されるＲＧＢ値を、ＲＧＢカメラ３２の分光感度特性に基づいて算出する。後述する第３実施形態及び第４実施形態でも同様である。なお、ＲＧＢカメラ３２の分光感度特性は、予め調べておくものとする。そして、異物混入部と異物無し部とのＲＧＢ値の差を表す評価量（例えば、両ＲＧＢ値間のユークリッド距離）を算出する。

【0126】

機能性光源３１、３１は、第１移動手段であるベルトコンベヤ３５上の所定の撮影領域を照明するように、互いに向き合う方向に少し傾斜されて配置されている。ＲＧＢカメラ３２は、真上から機能性光源３１、３１間を通して撮影領域を撮影するように配置されている。ベルトコンベヤ３５は、第２移動手段であるベルトコンベヤ３７とともに、検査対象であるブルーベリージャム３６の移動手段を構成し、ジャム３６は、ベルトコンベヤ３５上に乗せられて、撮影領域を通ってベルトコンベヤ３７上に落下する。撮影領域は、ジャム３６とその両側のベルトコンベヤ３５部分を含むように設定されている。ベルトコンベヤ３５とベルトコンベヤ３７との段差部分には、排除機構であるシュータ３８が配置されており、シュータ３８は排除機構制御ユニット３４と接続されている。シュータ３８は、上下方向に回動して開閉する扉部３８１を有し、扉部３８１が閉鎖されているときは、ジャム３６はベルトコンベヤ３５からベルトコンベヤ３７に落下するが、扉部３８１が開放されているとき、ジャム３６はベルトコンベヤ３５から扉部３８１を伝って排除物ボックス３９に排出され、ベルトコンベヤ３７に落下しないように構成されている。

【0127】

以上のように構成された検査装置３０の動作について説明する。ＲＧＢカメラ３２は、ベルトコンベヤ３５上のジャム３６の移動速度に合わせて、ジャム３６に未撮影部分が残らないように、所定の時間間隔で繰り返し撮影領域を撮影し、ＲＧＢ値を含む計測データである撮影画像を、画像処理部３３１に出力する。

【0128】

図１４（ａ）は、ＲＧＢカメラ３２の撮影画像の模式図であり、この撮影画像には、ジャム３６の画像部分３６’と、ベルトコンベヤ３５の画像部分３５’と、ジャム３６中の異物の画像部分１０’とが含まれている。画像処理部３３１は、まず、撮影画像の各部分の濃度値（強度）をＲＧＢ値から計算し、予め定められた閾値にしたがって、背景領域である画像部分３５’を、撮影画像から除外し、図１４（ｂ）に示すような画像とする。なお、図１４では、除外された部分を黒色で表している。閾値は、ジャム３６、ベルトコンベヤ３５及び異物の各画像部分が含まれ、それらの各々が画像中のどこに含まれているかが分かっている教師画像を用いて、事前に学習を行うことにより、決定しておく。学習方法として線形判別分析を用いるが、これに限らずサポートベクターマシンなど様々な公知の機械分類手法を使用可能である。

【0129】

次に、画像処理部３３１は、上記背景領域の除外方法と同様の手法を用いて、ジャム３６の画像部分３６’と異物の画像部分１０’とを閾値に従って判別する。この閾値も予め学習により決定しておく。画像処理部３３１は、この判別結果に基づいて、画像部分１０’に例えば白を、それ以外の部分である画像部分３６’、３７’に例えば黒をラベル付けする２値化処理を行い、図１４（ｃ）に示すような２値化処理画像を生成し、データ処理部３３２に出力する。

【0130】

データ処理部３３２は、図１５に示すように、２値化処理画像において異物の画像部分１０’が存在するか否かを判定し（ステップＳ１０１）、画像部分１０’が存在すれば、排除信号を排除機構制御ユニット３４に出力する（ステップＳ１０２）。排除信号を受け取った排除機構制御ユニット３４は、シュータ３８に扉部３８１を開放する制御信号を出力し、シュータ３８は扉部３８１を開放して、ジャム３６を排除物ボックス３９に排出する。なお、排除機構制御ユニット３４は、ジャム３６の異物混入と判定された部分が排出されるのに必要な所定の時間が経過すると、扉部３８１を閉鎖する制御信号を出力し、これにより扉部３８１は閉鎖する。

【0131】

例えば赤色の光源下で計測を行った場合、赤色方向への色の変化は検出し難いため、Ｒ値は有効に使えない等、白色以外の光源下で計測を行った場合には、ＲＧＢカメラ３２が出力するＲ値、Ｇ値及びＢ値を有効に使うことができないが、第２実施形態によれば、白色光源である機能性光源３１の下で計測を行うので、すなわち、白色を中心として色の変化を計測することとなるので、ＲＧＢカメラ３２が出力するＲ値、Ｇ値及びＢ値の全てを有効に検査結果の判定に用いることができ、検査精度を向上可能である。また、機能性光源３１は、同等の輝度を有する等エネルギー白色光源で照明したときよりも弁別対象間の計測値の差が拡大することからも、検査精度を向上可能である。また、検査装置３０は、近赤外線、赤外線や多層膜光学フィルタを用いる必要がないので、導入容易である。

【0132】

〈第３実施形態〉
次に、図１６に基づいて、第３実施形態に係る検査装置４０について説明する。検査装置４０は、検査対象が皮膚であって皮膚疾患を診断する皮膚疾患診断装置として構成されており、２つの機能性光源４１と、接写計測用レンズ４２１が取り付けられたＲＧＢカメラ４２と、ＲＧＢカメラ４２に接続されパーソナル・コンピュータからなる判定部４３と、判定部４３に接続された装置制御ユニット４４とを備えている。ＲＧＢカメラ４２は、ＲＧＢ値を出力可能な撮像装置であり、計測装置に相当する。判定部４３は、ノイズ・鏡面反射除去部４３１と、色素分布画像生成部４３２と、画像統計量計算部４３３と、症状推定部４３４と、表示部４３５とを備えている。

【0133】

各機能性光源４１は、ピーク波長の異なる３つ以上のＬＥＤが配列された基板４１１と、ＬＥＤからの光を混合する拡散板（図示せず。）とを備え、外部電源から電気の供給を受けて白色光を発光する。

【0134】

各機能性光源４１は、弁別の対象を人の皮膚における疾患部分と疾患無し部分として、第２実施形態で説明した設計方法により設計されたものである。機能性光源４１、４１は、所定の撮影領域を照明するように、互いに向き合う方向に少し傾斜されて配置されている。ＲＧＢカメラ４２は、機能性光源４１、４１間を通して撮影領域を撮影するように配置されている。なお、撮影領域には無色透明の透明板４５が配置され、ＲＧＢカメラ４２は、透明板４５を通して皮膚４６を撮影するように構成されている。

【0135】

以上のように構成された検査装置４０の動作について説明する。ＲＧＢカメラ４２は、撮影領域内の皮膚４６を撮影し、ＲＧＢ値を含むデータである撮影画像を、ノイズ・鏡面反射除去部４３１に出力する。

【0136】

ノイズ・鏡面反射除去部４３１は、例えばガウシアンフィルタや移動平均フィルタ，メディアンフィルタ等の一般的な画像処理用二次元フィルタを用いて、撮影画像からノイズを除去する。次に、ノイズ・鏡面反射除去部４３１は、撮影画像から、極端に肌色から逸脱し光源色を示す領域を、予め定めておいた閾値に従って除外する。判断には、撮影画像のＲＧＢ値を用いる。閾値は、鏡面反射有りの部分と鏡面反射無しの部分とが分かっている教師画像を用いて事前に学習を行うことにより、定めておく。学習手法として線形判別分析を用いるが、これに限らずサポートベクターマシンなど様々な公知の機械分類手法を使用可能である。ノイズ・鏡面反射除去部４３１は、ノイズ及び鏡面反射を除去した画像を、色素分布画像生成部４３２に出力する。

【0137】

色素分布画像生成部４３２は、画像のＲＧＢ値をＬ^＊ａ^＊ｂ^＊値に変換し、ａ^＊をヘモグロビン分布、Ｌ^＊をメラニン分布として、それぞれの色素分布画像を生成する（下記参考文献５参照）。なお、下記参考文献６に記載された色素分布算出手法を用いて、色素分布画像を生成することとしてもよい。この手法は、肌画像に対して独立成分分析を適用し、色素成分を取り出すというものである。色素分布画像生成部４３２は、生成した各色素分布画像を、画像統計量計算部４３３に出力する。

【0138】

画像統計量計算部４３３は、各色素分布画像から、平均、分散、歪度、尖度、空間周波数等の画像統計量を算出し、症状推定部４３４に出力する。空間周波数にはガボール・ウェーブレット変換を採用し、あらかじめ設定した振動方向についてそれぞれ周波数成分を計算して、症状推定の指標とする。

【0139】

症状推定部４３４は、予め、様々な症例について上記の画像統計量を計算し、機械学習を行っておくことで、症状の分類ができるように構成されている。学習手法には、公知の機械学習手法を用いればよく、例えば正準判別分析、サポートベクターマシン、ニューラルネットワークなどがある。症状推定部４３４は、画像統計量計算部４３３から受け取った画像統計量に基づいて、症状の分類を行い、その分類を診断結果（検査結果）として、表示部４３５に出力する。なお、症状推定部４３４は、画像統計量計算部４３３から色素分布画像を受け取って、表示部４３５に出力する。表示部４３５には、色素分布及び診断結果が表示される。

【0140】

なお、判定部４３は装置制御ユニット４４に計測開始信号を出力し、それに応じて装置制御ユニット４４はＲＧＢカメラ４２に画像取得命令信号を出力し、それに応じてＲＧＢカメラ４２は撮影を行い、画像取得完了信号を装置制御ユニット４４に出力するとともに、計測画像データである撮影画像を判定部４３に出力する。

【0141】

また、装置制御ユニット４４は、計測開始信号を受け取ったとき、ＬＥＤ制御信号を機能性光源４１、４１に出力し、これに応じて機能性光源４１、４１は各ＬＥＤを所定の輝度で発光する。なお、検出したい対象に応じて、ＬＥＤの選択及び輝度を切り替えられるように構成してもよく、その場合には、ＬＥＤの選択及び輝度を切り替えるためのＬＥＤ制御信号に応じて、ピーク波長が互いに異なるＬＥＤが３つ以上選択されるとともに、それらに対する印加電流が定められる。すなわち、ＬＥＤ制御信号を受信した機能性光源４１は、ピーク波長が互いに異なる３つ以上のＬＥＤを選択して、それらが発する光の発光条件を制御し、それらの光を混合することにより、白色光を発光する。

【0142】

第３実施形態によれば、白色光源である機能性光源４１の下で計測を行うので、すなわち、白色を中心として色の変化を計測することとなるので、ＲＧＢカメラ４２が出力するＲ値、Ｇ値及びＢ値の全てを有効に検査結果の判定に用いることができ、検査精度を向上可能である。また、機能性光源４１は、同等の輝度を有する等エネルギー白色光源で照明したときよりも弁別対象間の計測値の差が拡大することからも、検査精度を向上可能である。また、検査装置４０は、近赤外線、赤外線や多層膜光学フィルタを用いる必要がないので、導入容易である。

【0143】

なお、第３実施形態の参考文献として、先に挙げた参考文献５、６以外に、下記参考文献７、８、９がある。

【0144】

〈第４実施形態〉
次に、図１７に基づいて、第４実施形態に係る検査装置５０について説明する。検査装置５０は、検査対象が塗装面であって塗装むらを評価するための装置として構成されており、２つの機能性光源５１と、計測装置に相当するＲＧＢカメラ５２と、ＲＧＢカメラ５２に接続されパーソナル・コンピュータからなる判定部５３とを備えている。判定部５３は、ノイズ・鏡面反射除去部５３１と、不均一性指標計算部５３２と、均一性判定部５３３とを備えている。

【0145】

各機能性光源５１は、ピーク波長の異なる３つ以上のＬＥＤが配列された基板（図示せず。）と、ＬＥＤからの光を混合する拡散板５１２とを備え、外部電源から電気の供給を受けて白色光を発光する。

【0146】

各機能性光源５１は、弁別の対象を塗装面における塗装むら有り部分と塗装むら無し部分として、第２実施形態で説明した設計方法により設計されたものである。機能性光源５１、５１は、所定の撮影領域を照明するように、互いに向き合う方向に少し傾斜されて配置されている。ＲＧＢカメラ５２は、機能性光源５１、５１間を通して撮影領域を撮影するように配置されている。撮影領域には検査対象である塗装面５６が配置される。

【0147】

以上のように構成された検査装置５０の動作について説明する。ＲＧＢカメラ５２は、撮影領域に配置された塗装面５６を撮影し、ＲＧＢ値を含むデータである撮影画像を、ノイズ・鏡面反射除去部５３１に出力する。

【0148】

ノイズ・鏡面反射除去部５３１は、第３実施形態のノイズ・鏡面反射除去部４３１と同様にして、撮影画像からノイズ及び鏡面反射を除去し、その結果の画像を不均一性指標計算部５３２に出力する。

【0149】

不均一性指標計算部５３２は、ノイズ・鏡面反射除去部５３１から受け取った画像から、分散、歪度、尖度、空間周波数等の画像統計量を算出し、均一性判定部５３３に出力する。

【0150】

均一性判定部５３３は、均一性の有無を判定する。均一性判定のプロセスは、画像統計量から均一性指標を計算する手順と、その指標と閾値とから均一性の有無を判定する手順とからなる。なお、予め、複数のサンプルについて、均一性を主観的に数値評価した主観評価値とそのサンプルの画像統計量とを取得して、画像統計量から主観評価値を推定する推定式を構築しておく。推定には重回帰分析を使用すればよい。そして、均一性判定部５３３は、その推定式を用いて、不均一性指標計算部５３２から受け取った画像統計量に対し、均一性指標である主観評価値を推定する。若しくは、主観評価値を「均一／不均一」の二値とし、画像統計量に対して線形判別分析を適用し、それによって得られる判別スコアを均一性指標としてもよい。そして、均一性判定部５３３は、均一性指標に基づいて、事前に設定した閾値に従って、均一性の有無、すなわち、塗装むらの有無の判断を行い、その結果を指標とともに表示部（図示せず。）に表示する。この閾値は、教師データとして使用する主観評価値を得たときに、求めておいたものである。

【0151】

第４実施形態によれば、白色光源である機能性光源５１の下で計測を行うので、すなわち、白色を中心として色の変化を計測することとなるので、ＲＧＢカメラ５２が出力するＲ値、Ｇ値及びＢ値の全てを有効に検査結果の判定に用いることができ、検査精度を向上可能である。また、機能性光源５１は、同等の輝度を有する等エネルギー白色光源で照明したときよりも弁別対象間の計測値の差が拡大することからも、検査精度を向上可能である。また、検査装置５０は、近赤外線、赤外線や多層膜光学フィルタを用いる必要がないので、導入容易である。

【0152】

なお、機能性光源３１、４１、５１の光源色はいずれも白色としたが、上述した設計方法によれば、任意の色を有する機能性光源を設計可能であり、光源色を周囲の環境に応じた色とすれば、弁別対象間の計測値差が拡大されることと相まって、検査精度を向上させることができる。

【0153】

参考文献を、以下に挙げる。これらの参考文献の内容は、参照によってここに援用される。

【0154】

〈参考文献１〉F. J. M. Schmitt,“A method for the treatment of metamerism in colorimetry”，Journal of the Optical Society of America，1976年，Vol.66，p.601-608
〈参考文献２〉N. Moroney，M. D. Fairchild，R. W. G. Hunt，C. J. Li，M. R. Luo，T. Newman，“The CIECAM02 color appearance model”，IS&T 10th Color Imaging Conference，2002年，p.23‐27
〈参考文献３〉Changjun Li，M. Ronnier Luo，Guihua Cui，“Uniform colour spacesbased on CIECAM02 colour appearance model”，COLOR research and application，2006年，Vol.31，p.320-330
〈参考文献４〉CIE，“The CIE 1997 Interim Colour Appearance Model(Simple Version),CIECAM97s”，1998年4月
〈参考文献５〉Facial Skin Coloration Affects Perceived Health of Human Faces,Ian D. Stephen, et. al., Int J Primatol, 30:845-857, 2009.
〈参考文献６〉Tsumura N (2003) Image-based skin color and texture analysis/synthesis by extracting hemoglobin and melanin information in the skin. ACM SIGGRAPH 2003, 770-779. (computer graphics)
〈参考文献７〉Computer-Based Classification of Dermoscopy Images of Melanocytic Lesions on Acral Volar Skin. Iyatomi H, et al. Journal of Investigative Dermatology advance online publication, 2008.
〈参考文献８〉Moncrieff M., Cotton S., Claridge E. and Hall P. Spectrophotometric intracutaneous analysis: a new technique for imaging pigmented skin lesions.Brit. J. Derm. 146(3),448-457, 2002.
〈参考文献９〉Matts PJ, et al, Chromophore mapping: a new technique to characterize aging human skin, in vivo. J Am Acad Dermatol 53(3):86.

【符号の説明】

【0155】

１ 照明装置
２ 狭帯域光源
４，３１，４１，５１ 機能性光源
３０，４０，５０ 検査装置
３２，４２，５２ ＲＧＢカメラ（計測装置）
３３，４３，５３ 判定部

【0156】

【表1】

【0157】

【表2】

【0158】

【表3】

【0159】

【表4】

【0160】

【表5】

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】

【図15】

【図16】

【図17】