JP2020139963A - 分光装置を校正する方法および校正された分光装置を生産する方法 - Google Patents

Abstract

【課題】分光装置に備えられるセンサーの分光感度を精度よく適切に求める。【解決手段】分光装置は、複数の波長成分の各々のエネルギー量を示す信号を出力するセンサーを備える。複数のセンサーの各々について、センサーの基準分光感度が求められ、センサーの基準分光感度を示す指標を求められ、センサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれが分光装置の機械的誤差を示す指標の１次関数により表される。センサーの分光感度の集合がセンサーにより出力される信号の集合に適合するように分光装置の機械的誤差を示す指標が求められる。複数のセンサーの各々について、センサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれが求められ、センサーの分光感度が求められる。【選択図】図１

Description

本発明は、分光装置の校正に関する。

ポリクロメーターを備える分光装置は、複数の波長成分の各々のエネルギー量を示す信号を出力するセンサーを備える。複数のセンサーにより出力される複数の信号からは、分光反射率、分光放射輝度、分光透過率等の分光特性が求められる。

センサーの分光感度は、受光センサーの配置、形状、大きさ等により変化し、ポリクロメーターに備えられるスリット板、回折格子等の配置、形状、大きさ等により変化する。このため、分光特性の基礎となるセンサーの分光感度は、分光装置ごとに校正されなければならない。

特許文献１に記載された分光装置の波長校正方法においては、分光感度の中心波長の基準値からの差が画素番号の１次関数で与えられ（段落００２６）、分光感度の半値幅の基準値に対する比が画素番号の１次関数で与えられ（段落００２８）、分光感度の中心波長および分光感度の半値幅の基準値に対する比を用いて修正分光感度テーブルが作成され（段落００２９）、修正分光感度テーブルおよび輝線波長から算出される算出相対出力と実測相対出力とが互いに最も近い値となるように画素番号の１次関数の各々に含まれる係数が決められる（段落００３３）。

特許第４６６０６９４号明細書

特許文献１に記載された分光装置の波長校正方法においては、分光感度の中心波長の基準値からの差が画素番号の１次関数で与えられ、分光感度の半値幅の基準値に対する比が画素番号の１次関数で与えられる。しかし、実際の分光装置においては、当該差が画素番号の１次関数とは著しく異なる場合または当該比が画素番号の１次関数とは著しく異なる場合がある。これらの場合は、当該差を画素番号の１次関数で与え当該比を画素番号の１次関数で与えることは良好な近似ではなく、分光感度が精度よく求められない。当該差を画素番号の２次関数以上の高次関数で与え当該比を画素番号の２次関数以上の高次関数で与えることも考えられるが、高次関数が用いられた場合は、係数が修正された場合に算出相対出力が大きく変化し算出相対出力と実測相対出力とが互いに最も近い値となる係数が適切に求められず、分光感度が適切に求められない。

下記の発明は、この問題を解決するためになされる。下記の発明が解決しようとする課題は、分光感度を精度よく適切に求めることである。

校正の対象にされる分光装置は、光学系および受光センサーを備える。光学系は、被測定光をスペクトルに変換する。受光センサーは、スペクトルに含まれる複数の波長成分の各々のエネルギー量を示す信号を出力するセンサーを備えることにより複数の信号を出力する複数のセンサーを備える。

複数のセンサーの各々について、センサーの基準分光感度が求められ、求められたセンサーの基準分光感度を示す指標が求められ、求められたセンサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれが分光装置の機械的誤差を示す指標の１次関数により表されるモデルが作成される。センサーの分光感度を示す指標により示されるセンサーの分光感度がセンサーにより出力される信号に適合するように分光装置の機械的誤差を示す指標が求められる。複数のセンサーの各々について、作成されたモデルおよび求められた分光装置の機械的誤差を示す指標を用いて、センサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれが求められ、求められたセンサーの基準分光感度およびセンサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれを用いて、センサーの分光感度が求められる。

センサーの分光感度を示す指標が分光装置の機械的誤差を示す指標の１次関数により表されるため、分光装置の機械的誤差を示す指標を変化させた場合に分光感度の集合の信号の集合への適合性が大きく変化しない。このため、分光装置の機械的誤差を示す指標が適切に求められ、センサーの分光感度が精度よく適切に求められる。

分光装置および校正装置を示す模式図である。 分光ユニットを示す模式図である。 リニアアレイセンサーを示す模式図である。 信号処理機構を示すブロック図である。 センサーの分光感度および輝線スペクトルの例を示すグラフである。 基準分光感度と校正後の分光感度との関係を示す模式図である。 基準分光感度の中心波長と中心波長のずれ量との関係を示すグラフである。 基準分光感度の中心波長と半値幅比との関係を示すグラフである。 第３のモデルを作成する手順を示すフローチャートである。 機械誤差a1に起因するセンサーの各チャンネルの中心波長のずれ量を示すグラフである。 機械誤差a1に起因するセンサーの各チャンネルの半値幅のずれ量を示すグラフである。 機械誤差a2に起因するセンサーの各チャンネルの中心波長のずれ量を示すグラフである。 機械誤差a2に起因するセンサーの各チャンネルの半値幅のずれ量を示すグラフである。 機械誤差a3に起因するセンサーの各チャンネルの中心波長のずれ量を示すグラフである。 機械誤差a3に起因するセンサーの各チャンネルの半値幅のずれ量を示すグラフである。 機械誤差a4に起因するセンサーの各チャンネルの中心波長のずれ量を示すグラフである。 機械誤差a4に起因するセンサーの各チャンネルの半値幅のずれ量を示すグラフである。 機械誤差a5に起因するセンサーの各チャンネルの中心波長のずれ量を示すグラフである。 機械誤差a5に起因するセンサーの各チャンネルの半値幅のずれ量を示すグラフである。 分光装置の校正の手順を示すフローチャートである。 センサーの分光感度および輝線スペクトルの例を示すグラフである。 輝線光源測定時のセンサーのチャンネルと測定信号との関係を示すグラフである。 校正された分光装置を生産する手順を示すフローチャートである。

１ 分光装置
図１の模式図は、分光装置および校正装置を示す。図１には、分光装置の断面が示される。

分光装置１００は、図１に示されるように、分光ユニット１０４および信号処理機構１０５を備える。

分光装置１００が測定を行う場合は、分光ユニット１０４が、被測定光１０８を受光し、受光した被測定光１０８に含まれる第１の波長成分、第２の波長成分、・・・、第４０の波長成分のエネルギー量をそれぞれ示す第１の信号、第２の信号、・・・、第４０の信号を出力する。信号処理機構１０５は、第１の信号、第２の信号、・・・、第４０の信号の信号を用いて分光特性を求める。

分光装置１００が分光測色計である場合は、測定される物体に光が照射された場合に測定される物体により反射される光が被測定光１０８になり、分光反射率などが分光特性として求められる。分光装置１００が分光輝度計である場合は、測定される光源により放射される光が被測定光１０８になり、分光放射輝度などが分光特性として求められる。測定された物体に光が照射される場合に測定される物体を透過する光が被測定光１０８になってもよく、分光透過率などが分光特性として求められてもよい。測色値が分光特性として求められてもよい。測色値は、マンセル表色系、Ｌ^＊ａ^＊ｂ^＊表色系、Ｌ^＊Ｃ^＊ｈ表色系、ハンターＬａｂ表色系、ＸＹＺ表色系等で表現される。

２ 分光ユニット
図２の模式図は、分光ユニット１０４を示す斜視図である。

分光ユニット１０４は、図１および図２の各々に示されるように、光学系１１１およびリニアアレイセンサー１１２を備える。光学系１１１は、ポリクロメーターであり、スリット板１１５および凹面回折格子１１６を備える。

図３の模式図は、リニアアレイセンサーを示す平面図である。

リニアアレイセンサー１１２は、図３に示されるように、センサー１１９−１，１１９−２，・・・，１１９−４０を備える。

分光装置１００が測定を行う場合は、被測定光１０８がスリット板１１５に形成された矩形のスリット１２２に導かれる。スリット１２２に導かれた被測定光１０８は、スリット１２２を通過する。スリット１２２を通過した被測定光１０８は、スリット１２２から凹面回折格子１１６の回折面１２４まで進み、回折面１２４に反射される。被測定光１０８は、回折面１２４に反射されることにより、スペクトルに変換される。スペクトルに変換された被測定光１０８は、凹面回折格子１１６の回折面１２４からリニアアレイセンサー１１２の受光面１２５まで進み、受光面１２５において結像し、リニアアレイセンサー１１２に受光される。センサー１１９−１，１１９−２，・・・，１１９−４０は受光面１２５において方向１２８に直線状に配列される。スリット板１１５、凹面回折格子１１６およびリニアアレイセンサー１１２は、受光面１２５に結像する光の波長が方向１２８の位置に応じて変化するように配置される。このため、リニアアレイセンサー１１２がスペクトルに変換された被測定光１０８を受光した場合は、センサー１１９−１，１１９−２，・・・，１１９−４０が、互いに異なる第１の波長成分、第２の波長成分、・・・、第４０の波長成分のエネルギー量をそれぞれ示す第１の信号、第２の信号、・・・、第４０の信号をそれぞれ出力する。出力された第１の信号、第２の信号、・・・、第４０の信号は、信号処理機構１０５に入力される。信号処理機構１０５は、入力された第１の信号、第２の信号、・・・、第４０の信号を用いて分光特性を求める。

光学系１１１は、スリット１２２から回折面１２４に至る光軸１３４および回折面１２４から受光面１２５に至る光軸１３５を有する。

光学系１１１が他の種類の光学系に置き換えられてもよい。例えば、凹面回折格子１１６が平面回折格子および凹面鏡に置き換えられてもよい。また、スリット板１１５および凹面回折格子１１６が、円形のスリットが形成されたスリット板、シリンドリカルレンズおよびリニアバリアブルフィルターに置き換えられてもよい。円形のスリットを通過した被測定光は、シリンドリカルレンズを通過する。被測定光がシリンドリカルレンズを通過することにより、被測定光の断面形状が円形状から線状に変換される。シリンドリカルレンズを通過した被測定光は、リニアバリアブルフィルターを通過する。被測定光は、リニアバリアブルフィルターを通過することにより、スペクトルに変換される。

リニアアレイセンサー１１２が他の種類の受光センサーに置き換えられてもよい。例えば、リニアアレイセンサー１１２が、３９個以下または４１個以上のセンサーを備えるリニアアレイセンサーに置き換えられてもよい。光学系によっては、リニアアレイセンサー１１２がエリアセンサーに置き換えられてもよい。

３ 信号処理機構
図４のブロック図は、信号処理機構を示す。

信号処理機構１０５は、図４に示されるように、Ａ／Ｄ変換機構１３８および演算機構１３９を備える。

第１の信号、第２の信号、・・・、第４０の信号が信号処理機構１０５に入力された場合は、第１の信号、第２の信号、・・・、第４０の信号がＡ／Ｄ変換機構１３８に入力される。Ａ／Ｄ変換機構１３８に入力された第１の信号、第２の信号、・・・、第４０の信号は、第１の信号値、第２の信号値、・・・、第４０の信号値にそれぞれアナログ／デジタル変換される。第１の信号値、第２の信号値、・・・、第４０の信号値は、演算機構１３９に入力される。演算機構１３９は、入力された第１の信号値、第２の信号値、・・・、第４０の信号値および演算機構１３９に格納されたセンサー１１９−１の分光感度１４０−１、センサー１１９−２の分光感度１４０−２、・・・、センサー１１９−４０の分光感度１４０−４０を用いて分光特性を求める。分光感度１４０−１，１４０−２，・・・，１４０−４０に代えて、分光感度１４０−１，１４０−２，・・・，１４０−４０から導かれる情報であって分光特性を求めるために必要な情報が演算機構１３９に格納されてもよい。

演算機構１３９は、組み込みコンピューターであり、インストールされたプログラムにしたがって動作する。演算機構１３９が行う処理の全部または一部がプログラムを伴わない電子回路により行われてもよい。演算機構１３９が行う処理の全部または一部が手作業で行われてもよい。

４ センサーの分光感度の例
図５のグラフは、センサーの分光感度および輝線校正用の被測定光の輝線スペクトルの例を示す。

センサー１１９−１の分光感度１４１−１の中心波長、センサー１１９−２の分光感度１４１−２の中心波長、・・・、センサー１１９−４０の分光感度１４１−４０の中心波長は、図５に示されるように、互いに異なり、それぞれ約３５２ｎｍ、約３６３ｎｍ、・・・、約７４０ｎｍである。これにより、センサー１１９−１，１１９−２，・・・，１１９−４０は、互いに異なる第１の波長成分、第２の波長成分、・・・、第４０の波長成分のエネルギー量をそれぞれ示す第１の信号、第２の信号、・・・、第４０の信号をそれぞれ出力する。

５ 分光装置の校正の必要性
分光感度１４１−１，１４１−２，・・・，１４１−４０の各々は、スリット板１１５、凹面回折格子１１６およびリニアアレイセンサー１１２の配置、形状、大きさ等により変化する。このため、分光特性が精度よく求められるためには、演算機構１３９に格納された分光感度１４０−１，１４０−２，・・・，１４０−４０を、スリット板１１５、凹面回折格子１１６およびリニアアレイセンサーの配置、形状、大きさ等に応じて変更し、真の分光感度１４１−１，１４１−２，・・・，１４１−４０にそれぞれ近づけなければならない。以下では、演算機構１３９に格納された分光感度１４０−１，１４０−２，・・・，１４０−４０を真の分光感度１４１−１，１４１−２，・・・，１４１−４０にそれぞれ近づけることを分光装置１００の校正という。

６ 校正装置
校正装置１４２は、図１に示されるように、ＨｇＣｄランプ１４３および制御演算機構１４４を備える。

校正装置１４２を用いて分光装置１００の校正が行われる場合は、制御演算機構１４４がＨｇＣｄランプ１４３に被測定光１０８として輝線校正用の被測定光を放射させる。放射された輝線校正用の被測定光は、分光装置１００により測定される。輝線校正用の被測定光が分光装置１００により測定された場合は、第１の信号値、第２の信号値、・・・、第４０の信号値が演算機構１３９に入力される。入力された第１の信号値、第２の信号値、・・・、第４０の信号値は、演算機構１３９から制御演算機構１４４へ転送される。制御演算機構１４４は、転送されてきた第１の信号値、第２の信号値、・・・、第４０の信号値を用いてセンサー１１９−１の分光感度、センサー１１９−２の分光感度、・・・、センサー１１９−４０の分光感度を求める。求められたセンサー１１９−１の分光感度、センサー１１９−２の分光感度、・・・、センサー１１９−４０の分光感度は、制御演算機構１４４から演算機構１３９へ転送され、それぞれ演算機構１３９に新たに格納される分光感度１４０−１，１４０−２，・・・，１４０−３になる。これにより、分光装置１００の校正が行われた後は、演算機構１３９が第１の信号、第２の信号、・・・、第４０の信号および新たに求められたセンサー１１９−１の分光感度、センサー１１９−２の分光感度、・・・、センサー１１９−４０の分光感度を用いて分光特性を求めることができる。

制御演算機構１４４は、コンピューターであり、インストールされたプログラムにしたがって動作する。制御演算機構１４４が行う処理の全部または一部がプログラムを伴わない電子回路により行われてもよい。制御演算機構１４４が行う処理の全部または一部が手作業で行われてもよい。制御演算機構１４４が分光装置１００に内蔵されてもよい。

７ 輝線成分の波長
輝線校正用の被測定光は、図５に示されるように、輝線成分１４５−１，１４５−２，１４５−３，１４５−４，１４５−５および１４５−６を含む。輝線成分１４５−１，１４５−２，１４５−３，１４５−４，１４５−５および１４５−６の波長は、それぞれ４０４．５５ｎｍ，４３５．８４ｎｍ，５０８．５８ｎｍ，５４６．０７ｎｍ，５７８ｎｍおよび６４７．８５ｎｍである。輝線成分１４５−１，１４５−２，１４５−３，１４５−４，１４５−５および１４５−６は、分光装置１００の校正に用いられる。

輝線成分１４５−１，１４５−２，１４５−３，１４５−４，１４５−５および１４５−６以外の輝線成分が分光装置１００の校正に用いられてもよい。５個以下または７個以上の輝線成分が分光装置１００の校正に用いられてもよい。ＨｇＣｄランプ１４３以外の輝線光源から放射される光が輝線校正用の被測定光として用いられてもよい。輝線校正でない校正が行われ、輝線光源でない光源から放射される光が校正用の被測定光として用いられる場合もある。

８ 分光装置の校正において用いられるモデル
分光装置１００の校正においては、第１のモデル、第２のモデルまたは第３のモデルが用いられる。第１のモデルおよび第２のモデルは、参考例である。

９ センサーの分光感度を示す指標
第１のモデル、第２のモデルおよび第３のモデルの各々においては、センサー１１９−１，１１９−２，・・・，１１９−４０が位置ｉにより識別される。位置ｉは、互いに異なる４０個の値ｉ_１，ｉ_２，・・・，ｉ_４０のいずれかをとる。センサー１１９−１，１１９−２，・・・，１１９−４０が位置以外の指標により識別されてもよい。例えば、センサー１１９−１，１１９−２，・・・，１１９−４０が基準分光感度の中心波長、画素番号等により識別されてもよい。

第１のモデル、第２のモデルおよび第３のモデルの各々においては、位置ｉにあるセンサーの分光感度が、位置ｉにあるセンサーの分光感度の中心波長λ_Ｇ（ｉ）および半値幅ＦＷＨＭ（ｉ）により示される。中心波長λ_Ｇ（ｉ）および半値幅ＦＷＨＭ（ｉ）の各々は、位置ｉの関数である。

位置ｉにあるセンサーの分光感度は、独立変数が波長であり従属変数が感度であるガウス関数により良好に近似される。当該ガウス関数の形状は、中心波長および半値幅により決定される。このため、中心波長λ_Ｇ（ｉ）および半値幅ＦＷＨＭ（ｉ）は、位置ｉにあるセンサーの分光感度を示す指標として好適である。ただし、位置ｉにあるセンサーの分光感度を示す指標が変更されてもよい。位置ｉにあるセンサーの分光感度を示す指標が１個の変数または３個以上の変数からなることも許される。

１０ 第１のモデル（参考例）
第１のモデルが用いられる場合は、中心波長λ_Ｇ（ｉ）が位置ｉのｎ次関数で表される式（１）が作成される。

また、半値幅ＦＷＨＭ（ｉ）が位置ｉのｍ次関数により表される式（２）が作成される。

第１のモデルは、式（１）および（２）からなる。係数ａ_ｎ，・・・，ａ_０は、説明変数になり、位置ｉにあるセンサーの分光感度を決める校正パラメーターになる。係数ｂ_ｎ，・・・，ｂ_０は、説明変数になり、位置ｉにあるセンサーの分光感度を決める校正パラメーターになる。

第１のモデルを用いて分光装置１００の校正が行われる場合は、中心波長λ_Ｇ（ｉ）および半値幅ＦＷＨＭ（ｉ）により示されるセンサーの分光感度の２個以上のセンサーについての集合である分光感度の集合が、センサーにより出力される信号の２個以上のセンサーについての集合である信号の集合に適合するように、係数ａ_ｎ，・・・，ａ_０が求められ、係数ｂ_ｎ，・・・，ｂ_０が求められる。

続いて、作成された式（１）および求められた係数ａ_ｎ，・・・，ａ_０を用いて中心波長λ_Ｇ（ｉ）が求められ、作成された式（２）および求められた係数ｂ_ｎ，・・・，ｂ_０を用いて半値幅ＦＷＨＭ（ｉ）が求められる。

続いて、求められた中心波長λ_Ｇ（ｉ）および半値幅ＦＷＨＭ（ｉ）を用いて位置ｉにあるセンサーの分光感度が求められる。位置ｉにあるセンサーの分光感度は、求められた中心波長λ_Ｇ（ｉ）および半値幅ＦＷＨＭ（ｉ）により示される分光感度である。

第１のモデルを用いて分光装置１００の校正が行われた場合、輝線成分の波長λ_ＨｇＣｄ（１），λ_ＨｇＣｄ（２）、・・・，λ_ＨｇＣｄ（Ｋ_０）のいずれかに近い波長領域においては分光感度が精度よく求められるが、それ以外の波長領域においては分光感度が精度よく求められない。特に、最も短波長側にある波長領域および最も長波長側にある波長領域においては、分光感度が精度よく求められない。

１１ 第２のモデル（参考例）
第２のモデルは、第１のモデルの問題を解決するために提供される。

第２のモデルが用いられる場合は、スリット板１１５、凹面回折格子１１６およびリニアアレイセンサー１１２の配置、形状、大きさ等が設計狙いに一致する理想的な分光装置１００が想定され、想定された分光装置１００が備える位置ｉにあるセンサーの基準分光感度が光学シミュレーションにより求められる。位置ｉにあるセンサーの基準分光感度は、位置ｉの関数である。

続いて、位置ｉにあるセンサーの基準分光感度の中心波長λ_Ｇ０（ｉ）および半値幅ＦＷＨＭ_０（ｉ）が求められる。

中心波長λ_Ｇ０（ｉ）および半値幅ＦＷＨＭ_０（ｉ）は、位置ｉにあるセンサーの基準分光感度を示す指標として好適である。

続いて、中心波長λ_Ｇ０（ｉ）からの中心波長λ_Ｇ（ｉ）のずれΔλ_Ｇ（ｉ）が位置ｉの１次関数で表される式（３）が作成される。

また、半値幅ＦＷＨＭ_０（ｉ）に対する半値幅ＦＷＨＭ（ｉ）の比ｒａｔｉｏ（ｉ）が位置ｉの１次関数で表される式（４）が作成される。

第２のモデルは、式（３）および（４）からなる。係数ａ_１およびａ_０は、説明変数になり、位置ｉにあるセンサーの分光感度を決める校正パラメーターになる。係数ｂ_１およびｂ_０は、説明変数になり、位置ｉにあるセンサーの分光感度を決める校正パラメーターになる。

中心波長λ_Ｇ（ｉ）は、式（３）から導かれる式（５）を用いて求められる。

第２のモデルを用いて分光装置１００の校正が行われる場合は、中心波長λ_Ｇ（ｉ）および半値幅ＦＷＨＭ（ｉ）により示されるセンサーの分光感度の２個以上のセンサーについての集合である分光感度の集合が、センサーにより出力される信号の２個以上のセンサーについての集合である信号の集合に適合するように、係数ａ_１およびａ_０が求められ、係数ｂ_１およびｂ_０が求められる。

続いて、作成された式（３）ならびに求められた係数ａ_１およびａ_０を用いてずれΔλ_Ｇ（ｉ）が求められ、作成された式（４）ならびに求められた係数ｂ_１およびｂ_０を用いて比ｒａｔｉｏ（ｉ）が求められる。

続いて、求められたずれΔλ_Ｇ（ｉ）および比ｒａｔｉｏ（ｉ）を用いて位置ｉにあるセンサーの分光感度が求められる。

図６の模式図は、位置ｉにあるセンサーの基準分光感度と位置ｉにあるセンサーの分光感度との関係を示す。

位置ｉにあるセンサーの分光感度は、図６に示されるように、位置ｉにあるセンサーの基準分光感度を中心波長λ_Ｇ０（ｉ）を中心として波長軸方向にｒａｔｉｏ（ｉ）倍に拡大し、拡大された基準分光感度を波長軸方向にΔλ_Ｇ（ｉ）だけ移動させたものである。

図７のグラフは、実際の分光装置１００における基準分光感度の中心波長λ_Ｇ０とずれΔλ_Ｇとの関係を示す実測データである。図８のグラフは、実際の分光装置１００における基準分光感度の中心波長λ_Ｇ０と比ｒａｔｉｏとの関係を示す実測データである。

実際の分光装置１００においては、図７に示されるように、ずれΔλ_Ｇが基準分光感度の中心波長λ_Ｇ０の１次関数といいがたい場合があり、図８に示されるように、比ｒａｔｉｏが中心波長λ_Ｇ０の１次関数とはいいがたい場合がある。このような場合は、式（３）および（４）の各々は良好な近似ではない。このため、第２のモデルを用いて分光装置１００の校正が行われた場合は、分光感度が精度よく求められない場合がある。

この問題を解決するためにずれΔλ_Ｇ（ｉ）が位置ｉの２次関数以上の高次関数で表され比ｒａｔｉｏ（ｉ）が位置ｉの２次関数以上の高次関数で表された場合は、係数ａ_１，ａ_０，ｂ_１およびｂ_０が変化した場合に分光感度の集合の信号の集合への適合性が大きく変化する。このため、係数ａ_１，ａ_０，ｂ_１およびｂ_０が適切に求められない場合がある。

１２ 第３のモデル
１２．１ 第３のモデルの作成
第３のモデルは、第２のモデルの問題を解決するために提供される。

図９のフローチャートは、第３のモデルを作成する手順を示す。

複数の分光装置の各々について分光装置１００の校正が行われる場合でも、複数の分光装置の各々について図９のステップS1およびS2に示される手順を実行する必要はなく、図９に示されるステップS1およびS2の手順をあらかじめ１回だけ実行すればよい。ただし、複数の分光装置の各々について分光装置１００の校正が行われる間に図９に示される手順が再実行され、第３のモデルが再作成されてもよい。

第３のモデルが作成される場合は、ステップＳ１において、スリット板１１５、凹面回折格子１１６およびリニアアレイセンサー１１２の配置、形状、大きさ等が設計狙いに一致する理想的な分光装置１００が想定され、想定された分光装置１００が備える位置ｉにあるセンサーの基準分光感度が光学シミュレーションにより求められる。位置ｉにあるセンサーの基準分光感度は、位置ｉの関数である。このため、ステップＳ１により、位置ｉ_１にあるセンサー、位置ｉ_２にあるセンサー、・・・、位置ｉ_４０にあるセンサーの各々について、センサーの基準分光感度が求められる。

続いて、ステップＳ２において、位置ｉにあるセンサーの基準分光感度の中心波長λ_Ｇ０（ｉ）および半値幅ＦＷＨＭ_０（ｉ）が求められる。中心波長λ_Ｇ０（ｉ）および半値幅ＦＷＨＭ_０（ｉ）の各々は、位置ｉの関数である。このため、ステップＳ２により、位置ｉ_１にあるセンサー、位置ｉ_２にあるセンサー、・・・、位置ｉ_４０にあるセンサーの各々について、センサーの基準分光感度の中心波長および半値幅が求められる。

中心波長λ_Ｇ０（ｉ）および半値幅ＦＷＨＭ_０（ｉ）は、位置ｉにあるセンサーの基準分光感度を示す指標として好適である。

続いて、ステップＳ３において、中心波長λ_Ｇ０（ｉ）からの中心波長λ_Ｇ（ｉ）のずれがリニアアレイセンサー１１２の配置誤差ａ_１，ａ_２およびａ_３、スリット１２２の幅の製作誤差ａ_４ならびに凹面回折格子１１６の配置誤差ａ_５の１次関数により表される式（６）が作成される。また、半値幅ＦＷＨＭ_０（ｉ）からの半値幅ＦＷＨＭ（ｉ）のずれがリニアアレイセンサー１１２の配置誤差ａ_１，ａ_２およびａ_３、スリット１２２の幅の製作誤差ａ_４ならびに凹面回折格子１１６の配置誤差ａ_５の１次関数により表される式（７）が作成される。中心波長λ_Ｇ０（ｉ）および半値幅ＦＷＨＭ_０（ｉ）の各々は、位置ｉの関数である。誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５は、分光装置１００の機械的誤差を示す指標である。このため、ステップＳ３により、位置ｉ_１にあるセンサー、位置ｉ_２にあるセンサー、・・・、位置ｉ_４０にあるセンサーの各々について、センサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれが分光装置１００の機械的誤差を示す指標の１次関数により表される第３のモデルが作成される。

第３のモデルは、式（６）および（７）からなる。誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５は、説明変数になり、位置ｉにあるセンサーの分光感度を決める校正パラメーターになる。

リニアアレイセンサー１１２の配置誤差ａ_１は、センサー１１９−１，１１９−２，・・・，１１９−４０が配列される方向１２８についてのリニアアレイセンサー１１２の配置誤差である。リニアアレイセンサー１１２の配置誤差ａ_２は、光軸１３５が伸びる方向１２９についてのリニアアレイセンサーの配置誤差である。リニアアレイセンサー１１２の配置誤差ａ_３は、方向１２８と光軸１３５とに直交するセンサーの軸１３０ｃ周り１３０の回転誤差である。スリット１２２の幅の製作誤差ａ_４は、凹面回折格子１１６の主断面１４６と平行をなし光軸１３４と垂直をなす方向１３１についてのスリット１２２の幅の製作誤差である。配置誤差ａ_５は、凹面回折格子１１６の軸１３２ｃ周り１３２の回転誤差である。凹面回折格子１１６の主断面１４６は、回折面１２４に形成される刻線と垂直をなす面である。

誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５は小さいことが期待されるため、中心波長λ_Ｇ０（ｉ）からの中心波長λ_Ｇ（ｉ）のずれを誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５の１次関数により表すことは良好な近似であり、半値幅ＦＷＨＭ_０（ｉ）からの半値幅ＦＷＨＭ（ｉ）のずれを誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５の１次関数により表すことも良好な近似である。このため、第３のモデルを用いて分光装置１００の校正が行われた場合は、分光感度が精度よく求められる。

誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５は、中心波長λ_Ｇ（ｉ）または半値幅ＦＷＨＭ（ｉ）への影響が大きく、位置ｉにあるセンサーの分光感度への影響が大きい。このため、誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５は、分光装置１００の機械的誤差を示す指標として好適である。ただし、分光装置１００の機械的誤差を示す指標が変更されてもよい。分光装置１００の機械的誤差を示す指標に含まれる説明変数または校正パラメーターの数は、輝線校正用の被測定光に含まれる輝線成分の数以下に制限されるが、４個以下または６個以上であってもよい。

式（６）に含まれる誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５の１次関数においては、誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５に、係数δλ_Ｇ１（ｉ），δλ_Ｇ２（ｉ），δλ_Ｇ３（ｉ），δλ_Ｇ４（ｉ）およびδλ_Ｇ５（ｉ）がそれぞれ乗じられる。係数δλ_Ｇ１（ｉ），δλ_Ｇ２（ｉ），δλ_Ｇ３（ｉ），δλ_Ｇ４（ｉ）およびδλ_Ｇ５（ｉ）は、それぞれ誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５が単位量発生した場合の中心波長λ_Ｇ（ｉ）のずれ量であり、光学シミュレーションにより求められる。係数δλ_Ｇ１（ｉ），δλ_Ｇ２（ｉ），δλ_Ｇ３（ｉ），δλ_Ｇ４（ｉ）およびδλ_Ｇ５（ｉ）の各々は、位置ｉの関数である。

式（７）に含まれる誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５の１次関数においては、誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５に、係数δＦＷＨＭ_１（ｉ），δＦＷＨＭ_２（ｉ），δＦＷＨＭ_３（ｉ），δＦＷＨＭ_４（ｉ）およびδＦＷＨＭ_５（ｉ）がそれぞれ乗じられる（加減算される）。係数δＦＷＨＭ_１（ｉ），δＦＷＨＭ_２（ｉ），δＦＷＨＭ_３（ｉ），δＦＷＨＭ_４（ｉ）およびδＦＷＨＭ_５（ｉ）は、それぞれ誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５が単位量発生した場合の半値幅ＦＷＨＭ（ｉ）のずれ量であり、光学シミュレーションにより求められる。係数δＦＷＨＭ_１（ｉ），δＦＷＨＭ_２（ｉ），δＦＷＨＭ_３（ｉ），δＦＷＨＭ_４（ｉ）およびδＦＷＨＭ_５（ｉ）の各々は、位置ｉの関数である。

図１０から図１９までの各々におけるチャンネルＣｈは、センサー１１９−１，１１９−２，・・・，１１９−４０の識別番号である。

図１０の１４９は、配置誤差ａ_１が０から＋１単位ずれた場合の各Ｃｈの分光感度の中心波長のずれ量δλ_Ｇを示し、１５０は配置誤差ａ１が０から−１単位ずれた場合の各Ｃｈの分光感度の中心波長のずれ量δλ_Ｇを示す。関係１４９および１５０は、光学シミュレーションにより求められる。係数δλ_Ｇ１（ｉ）は、例えば、関係１４９および１５０を用いて求められる。

図１１の１５３は、配置誤差ａ_１が０から＋１単位ずれた場合のセンサーの各Ｃｈの分光感度の半値幅のずれ量δＦＷＨＭを示し、１５４は、配置誤差ａ_１が０から−１単位ずれた場合のセンサーの各Ｃｈの分光感度の半値幅のずれ量δＦＷＨＭを示す。関係１５３および１５４は、光学シミュレーションにより求められる。係数δＦＷＨＭ_１（ｉ）は、例えば、関係１５３および１５４を用いて求められる。

図１２の１５７は、配置誤差ａ_２が０から＋１単位ずれた場合の各Ｃｈの分光感度の中心波長のずれ量δλ_Ｇを示し、１５８は配置誤差ａ_２が０から−１単位ずれた場合の各Ｃｈの分光感度の中心波長のずれ量δλ_Ｇを示す。関係１５７および１５８は、光学シミュレーションにより求められる。係数δλ_Ｇ２（ｉ）は、例えば、関係１５７および１５８を用いて求められる。

図１３の１６１は、配置誤差ａ_２が０から＋１単位ずれた場合のセンサーの各Ｃｈの分光感度の半値幅のずれ量δＦＷＨＭを示し、１６２は配置誤差ａ_２が０から−１単位ずれた場合のセンサーの各Ｃｈの分光感度の半値幅のずれ量δＦＷＨＭを示す。関係１６１および１６２は、光学シミュレーションにより求められる。係数δＦＷＨＭ_２（ｉ）は、例えば、関係１６１および１６２を用いて求められる。

図１４の１６５は、配置誤差ａ_３が０から＋１単位ずれた場合の各Ｃｈの分光感度の中心波長のずれ量δλ_Ｇを示し、１６６は配置誤差ａ_３が０から−１単位ずれた場合の各Ｃｈの分光感度の中心波長のずれ量δλ_Ｇを示す。関係１６５および１６６は、光学シミュレーションにより求められる。係数δλ_Ｇ３（ｉ）は、例えば、関係１６５および１６６を用いて求められる。

図１５の１６９は、配置誤差ａ_３が０から＋１単位ずれた場合のセンサーの各Ｃｈの分光感度の半値幅のずれ量δＦＷＨＭを示し、１７０は配置誤差ａ_３が０から−１単位ずれた場合のセンサーの各Ｃｈの分光感度の半値幅のずれ量δＦＷＨＭを示す。関係１６９および１７０は、光学シミュレーションにより求められる。係数δＦＷＨＭ_３（ｉ）は、例えば、関係１６９および１７０を用いて求められる。

図１６の１７３は、製作誤差ａ_４が０から＋１単位ずれた場合の各Ｃｈの分光感度の中心波長のずれ量δλ_Ｇを示し、１７４は製作誤差ａ_４が０から−１単位ずれた場合の各Ｃｈの分光感度の中心波長のずれ量δλ_Ｇを示す。関係１７３および１７４は、光学シミュレーションにより求められる。係数δλ_Ｇ４（ｉ）は、例えば、関係１７３および１７４を用いて求められる。

図１７の１７７は、製作誤差ａ_４が０から＋１単位ずれた場合のセンサーの各Ｃｈの分光感度の半値幅のずれ量δＦＷＨＭのチャンネルＣｈを示し、１７８は製作誤差ａ_４が０から−１単位ずれた場合のセンサーの各Ｃｈの分光感度の半値幅のずれ量δＦＷＨＭを示す。関係１７７および１７８は、光学シミュレーションにより求められる。係数δＦＷＨＭ_４（ｉ）は、例えば、関係１７７および１７８を用いて求められる。

図１８の１８１は、配置誤差ａ_５が０から＋１単位ずれた場合の各Ｃｈの分光感度の中心波長のずれ量δλ_Ｇを示し、１８２は配置誤差ａ_５が０から−１単位ずれた場合の各Ｃｈの分光感度の中心波長のずれ量δλ_Ｇを示す。関係１８１および１８２は、光学シミュレーションにより求められる。係数δλ_Ｇ５（ｉ）は、例えば、関係１８１および１８２を用いて求められる。

図１９の１８５は、配置誤差ａ_５が０から＋１単位ずれた場合のセンサーの各Ｃｈの分光感度の半値幅のずれ量δＦＷＨＭを示し、１８６は配置誤差ａ_５が０から−１単位ずれた場合のセンサーの各Ｃｈの分光感度の半値幅のずれ量δＦＷＨＭを示す。関係１８５および１８６は、光学シミュレーションにより求められる。係数δＦＷＨＭ_５（ｉ）は、例えば、関係１８５および１８６を用いて求められる。

１２．２ 第３のモデルを用いた分光装置の校正
図２０のフローチャートは、第３のモデルを用いた分光装置の校正の手順を示す。

複数の分光装置の各々について分光装置１００の校正が行われる場合は、複数の分光装置の各々について図２０に示される手順が実行される。

第３のモデルを用いて分光装置１００の校正が行われる場合は、ステップＳ１１において、中心波長λ_Ｇ（ｉ）および半値幅ＦＷＨＭ（ｉ）により示されるセンサーの分光感度が、センサーにより出力される信号に適合するように、誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５が求められる。当該分光感度を当該信号に適合させることは、位置ｉにあるセンサーの真の分光感度が中心波長λ_Ｇ（ｉ）および半値幅ＦＷＨＭ（ｉ）により示されるセンサーの分光感度である場合にセンサーにより出力されると想定される信号を、センサーにより実際に出力される信号に近づけることを意味する。適合性の評価には、目的変数が用いられる。ステップＳ１１により、センサーの分光感度を示す指標により示されるセンサーの分光感度が、センサーにより出力される信号に適合するように分光装置１００の機械的誤差を示す指標が求められる。

校正にはセンサー１１９−１，１１９−２，・・・，１１９−４０のうちの輝線成分の波長λ_ＨｇＣｄ（１），λ_ＨｇＣｄ（２），・・・，λ_ＨｇＣｄ（Ｋ_０）において感度を有するセンサーの出力を用いる。

続いて、ステップＳ１２において、作成された式（６）および求められた誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５を用いて、中心波長λ_Ｇ０（ｉ）からの中心波長λ_Ｇ（ｉ）のずれが求められる。また、作成された式（７）および求められた誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５を用いて、半値幅ＦＷＨＭ_０（ｉ）からの半値幅ＦＷＨＭ（ｉ）のずれが求められる。ステップＳ１２により、位置ｉ_１のセンサー、位置ｉ_２のセンサー、・・・、位置ｉ_４０のセンサーの各々について、作成された第３のモデルおよび求められた分光装置１００の機械的誤差を示す指標を用いて、センサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれが求められる。

続いて、ステップＳ１３において、求められた位置ｉにあるセンサーの基準分光感度、中心波長λ_Ｇ０（ｉ）からの中心波長λ_Ｇ（ｉ）のずれおよび半値幅ＦＷＨＭ_０（ｉ）からの半値幅ＦＷＨＭ（ｉ）のずれを用いて、位置ｉにあるセンサーの分光感度が求められる。これにより、位置ｉ_１にあるセンサー、位置ｉ_２にあるセンサー、・・・、位置ｉ_４０にあるセンサーの各々について、求められたセンサーの基準分光感度およびセンサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれを用いてセンサーの分光感度が求められる。

求められる位置ｉにあるセンサーの分光感度は、図６に示されるように、基準分光感度を中心波長λ_Ｇ０（ｉ）を中心として波長軸方向にｒａｔｉｏ（ｉ）倍に拡大し、拡大された基準分光感度を波長軸方向にΔλ_Ｇ（ｉ）だけ移動させたものである。

比ｒａｔｉｏ（ｉ）は、式（８）により表される。

第３のモデルを用いて分光装置１００の校正が行われた場合は、中心波長λ_Ｇ０（ｉ）および半値幅ＦＷＨＭ_０（ｉ）の各々が誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５の１次関数により表されるため、誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５を変化させた場合に分光感度の信号への適合性が大きく変化しない。このため、誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５が適切に求められ、センサーの校正後の分光感度が適切に求められる。

また、第３のモデルを用いて分光装置１００の校正が行われた場合は、中心波長λ_Ｇ０（ｉ）および半値幅ＦＷＨＭ_０（ｉ）の各々が共通の説明変数である誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５で表されるため、中心波長と半値幅を別の説明変数で表す場合に比べて説明変数の数が減少し、分光装置１００の校正に必要な輝線成分の数が減少する。

さらに、第３のモデルを用いて分光装置１００の校正が行われた場合は、中心波長λ_Ｇ０（ｉ）および半値幅ＦＷＨＭ_０（ｉ）が共通の説明変数である誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５で表されるため、中心波長λ_Ｇ０（ｉ）および半値幅ＦＷＨＭ_０（ｉ）の関係が不適切なものにならず、中心波長λ_Ｇ０（ｉ）および半値幅ＦＷＨＭ_０（ｉ）が適切に求められる。

１３ 分光感度を信号に適合させる手順
以下では、波長がλ_ＨｇＣｄ（ｋ）である輝線成分が位置Ｉ_ｋにあるセンサーおよび位置Ｉ_ｋ＋１にあるセンサーにまたがって入射するとする。位置Ｉ_ｋ＋１にあるセンサーの基準分光感度の中心波長は、位置Ｉ_ｋにあるセンサーの基準分光感度の中心波長に隣接する。輝線成分の識別番号ｋは、１，２，・・・，Ｋ_０のいずれかの値をとるとする。

図２１のグラフは、センサーの分光感度および輝線校正用の被測定光の輝線スペクトルの例を示す。図２１のグラフは、図５のグラフの波長が３９０ｎｍから４２０ｎｍまでである範囲を拡大したものである。

波長が４０４．５４ｎｍである輝線成分１４５−１は、中心波長が約３９６ｎｍである分光感度１４１−５を有するセンサーおよび中心波長が約４０７ｎｍである分光感度１４１−６を有するセンサーにまたがって入射する。このため、分光感度１４１−５を有するセンサーおよび分光感度１４１−６を有するセンサーの各々は、図２１に示されるように、輝線成分１４５−１に対して感度を有する。

図２２のグラフは、センサーのチャンネルとセンサーにより出力される信号との関係を示す。

チャンネル５のセンサーおよびチャンネル６のセンサーの各々は輝線成分１４５−１に対して感度を有し、チャンネル８のセンサーおよびチャンネル９のセンサーの各々は輝線成分１４５−２に対して感度を有し、チャンネル１５およびチャンネル１６のセンサーの各々は輝線成分１４５−３に対して感度を有し、チャンネル１９のセンサーおよびチャンネル２０のセンサーの各々は輝線成分１４５−４に対して感度を有し、チャンネル２２のセンサーおよびチャンネル２３のセンサーは輝線成分１４５−５に対して感度を有し、チャンネル２９のセンサーおよびチャンネル３０のセンサーの各々は輝線成分１４５−６に対して感度を有する。その結果、図２２に示されるような信号が得られる。

波長がλ_ＨｇＣｄ（１）である輝線成分、波長がλ_ＨｇＣｄ（２）である輝線成分、・・・、波長がλ_ＨｇＣｄ（Ｋ_０）である輝線成分を用いて分光装置１００の校正が完全に行われた場合は、位置ｉにあるセンサーの分光感度Ｒｅｓｐｏｎｓｅ（ｉ，λ）および位置ｉにあるセンサーにより出力される信号をアナログ／デジタル変換することにより得られる信号値Ｃｏｕｎｔ（ｉ）は、ｋ＝１，２，・・・，Ｋ_０の各々について、式（９）および（１０）に示される関係を満たす。

式（９）および（１０）は、ｋ＝１，２，・・・，Ｋ_０の各々について、波長λ_ＨｇＣｄ（ｋ）において感度を有する位置Ｉ_ｋにあるセンサーおよび位置Ｉ_ｋ＋１にあるセンサーを選択した場合に、位置Ｉ_ｋにあるセンサーの分光感度の波長λ_ＨｇＣｄ（ｋ）における感度Ｒｅｓｐｏｎｓｅ（Ｉ_ｋ，λ_ＨｇＣｄ（ｋ））が位置Ｉ_ｋにあるセンサーにより出力される信号を示す信号値Ｃｏｕｎｔ（Ｉ_ｋ）に一致し、位置Ｉ_ｋ＋１にあるセンサーの分光感度の波長λ_ＨｇＣｄ（ｋ）における感度Ｒｅｓｐｏｎｓｅ（Ｉ_ｋ＋１，λ_ＨｇＣｄ（ｋ））が位置Ｉ_ｋ＋１にあるセンサーにより出力される信号を示す信号値Ｃｏｕｎｔ（Ｉ_ｋ＋１）に一致することを示す。

分光感度の集合を信号の集合に適合させることは、感度Ｒｅｓｐｏｎｓｅ（Ｉ_ｋ，λ_ＨｇＣｄ（ｋ））および信号値Ｃｏｕｎｔ（Ｉ_ｋ）の関係を式（９）に示される関係に近づけ、感度Ｒｅｓｐｏｎｓｅ（Ｉ_ｋ＋１，λ_ＨｇＣｄ（ｋ））および信号値Ｃｏｕｎｔ（Ｉ_ｋ＋１）の関係を式（１０）に示される関係に近づけることを意味する。

このため、分光感度を信号に適合させる第１の方法は、式（１１）に示される目的変数Ｆを最小にする誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５を求めることである。

式（１１）に示される目的変数Ｆは、信号値Ｃｏｕｎｔ（Ｉ_ｋ）からの感度Ｒｅｓｐｏｎｓｅ（Ｉ_ｋ，λ_ＨｇＣｄ（ｋ））のずれの２乗と信号値Ｃｏｕｎｔ（Ｉ_ｋ＋１）からの感度Ｒｅｓｐｏｎｓｅ（Ｉ_ｋ＋１，λ_ＨｇＣｄ（ｋ））のずれの２乗との和の、ｋ＝１，２，・・・，Ｋ_０についての合計である。

ずれの２乗が、ずれの絶対値が大きくなるほど絶対値が大きくなる他の因子に置き換えられてもよい。例えば、ずれの２乗がずれの絶対値に置き換えられてもよい。

もしある輝線波長λ_HgCd(k)に対して感度を有するセンサーが３つ(I_k-1, I_k, I_k+1)ある場合には信号Count(I_k-1)と感度Response(I_k-1)とのずれ量を目的関数に加えてもよい。

また、式（９）および（１０）からは式（１２）が導かれる。

分光感度を信号に適合させることは、感度Ｒｅｓｐｏｎｓｅ（Ｉ_ｋ，λ_ＨｇＣｄ（ｋ））およびＲｅｓｐｏｎｓｅ（Ｉ_ｋ＋１，λ_ＨｇＣｄ（ｋ））ならびに信号値Ｃｏｕｎｔ（Ｉ_ｋ）およびＣｏｕｎｔ（Ｉ_ｋ＋１）の関係を式（１２）に示される関係に近づけることを意味する。

このため、分光感度を信号に適合させる第２の方法は、式（１３）に示される目的変数Ｆを最小にする誤差ａ_１，ａ_２，ａ_３，ａ_４およびａ_５を求めることである。

式（１３）に示される目的関数Ｆは、信号値Ｃｏｕｎｔ（Ｉ_ｋ）に対する信号値Ｃｏｕｎｔ（Ｉ_ｋ＋１）の比からの感度Ｒｅｓｐｏｎｓｅ（Ｉ_ｋ，λ_ＨｇＣｄ（ｋ））に対する感度Ｒｅｓｐｏｎｓｅ（Ｉ_ｋ＋１，λ_ＨｇＣｄ（ｋ））の比のずれの、ｋ＝１，２，・・・，Ｋ_０についての合計である。

式（１３）に示される目的関数Ｆによれば、感度Ｒｅｓｐｏｎｓｅ（Ｉ_ｋ，λ_ＨｇＣｄ（ｋ））およびＲｅｓｐｏｎｓｅ（Ｉ_ｋ＋１，λ_ＨｇＣｄ（ｋ））を信号値Ｃｏｕｎｔ（Ｉ_ｋ）およびＣｏｕｎｔ（Ｉ_ｋ＋１）とそれぞれ対比できるようにするために分光感度Ｒｅｓｐｏｎｓｅ（ｉ，λ）を規格化する必要がない。

もしある輝線波長λ_HgCd(k)に対して感度を有するセンサーが３つ(I_k-1, I_k, I_k+1)ある場合には、信号Count(I_k)とCount(I_k+1)との比と感度Response(I_k)とResponse(I_k+1)との比のずれ量に加えて、信号Count(I_k)とCount(I_k-1)との比と感度Response(I_k)とResponse(I_k-1)との比のずれ量とのずれ量を目的関数に加えてもよい。

１４ 校正された分光装置の生産
校正された分光装置の生産においては、図２３に示されるように、ステップＳ２１において分光装置１００が準備され、ステップＳ２２において準備された分光装置１００が校正される。ステップＳ２１においては、分光装置１００の校正を行う事業者が分光装置１００を製作することにより分光装置１００が準備されてもよいし、分光装置１００の校正を行う事業者が他の事業者から分光装置１００を仕入れることにより分光装置１００が準備されてもよい。

１００ 分光装置
１１１ 光学系
１１２ リニアアレイセンサー
１１５ スリット板
１１６ 凹面回折格子
１１９−１，１１９−２，・・・，１１９−４０ センサー
１２２ スリット
１４３ ＨｇＣｄランプ

Claims (10)

1. a) 被測定光をスペクトルに変換する光学系と、スペクトルに含まれる複数の波長成分の各々のエネルギー量を示す信号を出力するセンサーを備えることにより複数の信号を出力する複数のセンサーを備える受光センサーと、を備える分光装置を校正の対象にする工程と、
   b) 前記複数のセンサーの各々について、センサーの基準分光感度を求める工程と、
   c) 前記複数のセンサーの各々について、工程b)において求められたセンサーの基準分光感度を示す指標を求める工程と、
   d) 前記複数のセンサーの各々について、工程c)において求められたセンサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれが前記分光装置の機械的誤差を示す指標の１次関数により表されるモデルを作成する工程と、
   e) 工程d)において作成されたモデルに含まれるセンサーの分光感度を示す指標により示されるセンサーの分光感度がセンサーにより出力される信号に適合するように前記分光装置の機械的誤差を示す指標を求める工程と、
   f) 前記複数のセンサーの各々について、工程d)において作成されたモデルおよび工程e)において求められた前記分光装置の機械的誤差を示す指標を用いて、センサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれを求める工程と、
   g) 前記複数のセンサーの各々について、工程b)において求められたセンサーの基準分光感度および工程f)において求められたセンサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれを用いて、センサーの分光感度を求める工程と、
   を備える分光装置を校正する方法。
2. 工程c)において求められたセンサーの基準分光感度を示す指標は、センサーの基準分光感度の中心波長および半値幅であり、
   工程d)において作成されたモデルに含まれるセンサーの分光感度を示す指標は、センサーの分光感度の中心波長および半値幅であり、
   工程d)は、前記複数のセンサーの各々について、センサーの基準分光感度の中心波長からのセンサーの分光感度の中心波長のずれが前記分光装置の機械的誤差を示す指標の第１の１次関数により表される第１の式を作成し、センサーの基準分光感度の半値幅からのセンサーの分光感度の半値幅のずれが前記分光装置の機械的誤差を示す指標の第２の１次関数により表される第２の式を作成し、
   工程d)において作成されるモデルは、第１の式および第２の式からなる
   請求項１の分光装置を校正する方法。
3. 前記受光センサーは、受光面を有し、
   前記複数のセンサーは、前記受光面において第１の方向に配列され、
   前記光学系は、第２の方向に伸び前記受光面に至る光軸を有し、
   工程d)において作成されたモデルに含まれる前記分光装置の機械的誤差を示す指標は、前記第１の方向についての前記受光センサーの配置誤差、前記第２の方向についての前記受光センサーの配置誤差および前記第１の方向と前記受光面に至る光軸とに直交する軸の周りを周る回転方向についての前記受光センサーの配置誤差の全部または一部を含む
   請求項１または２の分光装置を校正する方法。
4. 前記光学系は、
   スリットが形成されたスリット板と、
   主断面および回折面を有する回折格子と、
   を備え、
   前記光学系は、前記スリットから前記回折面に至る光軸を有し、
   工程d)において作成されたモデルに含まれる前記分光装置の機械的誤差を示す指標は、前記主断面と平行をなし前記スリットから前記回折面に至る光軸と垂直をなす方向についての前記スリットの幅の製作誤差を含む
   請求項１から３までのいずれかの分光装置を校正する方法。
5. 前記光学系は、
   主断面を有する回折格子
   を備え、
   工程d)において作成されたモデルに含まれる前記分光装置の機械的誤差を示す指標は、前記主断面が同一平面上に維持される回転方向についての前記回折格子の配置誤差を含む
   請求項１から４までのいずれかの分光装置を校正する方法。
6. 工程e)は、
   複数の輝線成分を含む輝線校正用の被測定光を前記分光装置に測定させ、
   前記複数の輝線成分の各々について、輝線成分の波長において感度を有する第１のセンサーおよび第２のセンサーを選択し、第１のセンサーの分光感度の輝線成分の波長における感度である第１の感度、第２のセンサーの分光感度の輝線成分の波長における感度である第２の感度、第１のセンサーにより出力される信号である第１の信号、第２のセンサーにより出力される信号である第２の信号を用いて、第１の信号に対する第２の信号の比からの第１の感度に対する第２の感度の比のずれの前記複数の輝線成分についての合計が最小となるように前記分光装置の機械的誤差を示す指標を求めることにより、分光感度の集合を信号の集合に適合させる
   請求項１から５までのいずれかの分光装置を校正する方法。
7. 工程e)は、
   複数の輝線成分を含む輝線校正用の被測定光を前記分光装置に測定させ、
   前記複数の輝線成分の各々について、輝線成分の波長において感度を有する第１のセンサーおよび第２のセンサーを選択し、第１のセンサーの分光感度の輝線成分の波長における感度である第１の感度、第２のセンサーの分光感度の輝線成分の波長における感度である第２の感度、第１のセンサーにより出力される信号である第１の信号、第２のセンサーにより出力される信号である第２の信号を用いて、第１の信号からの第１の感度のずれの絶対値が大きくなるほど絶対値が大きくなる因子と第２の信号からの第２の感度のずれの絶対値が大きくなるほど絶対値が大きくなる因子との和の前記複数の輝線成分についての合計が最小になるように前記分光装置の機械的誤差を示す指標を求めることにより、分光感度の集合を信号の集合に適合させる
   請求項１から５までのいずれかの分光装置を校正する方法。
8. 前記輝線校正用の被測定光が、ＨｇＣｄランプにより放射され、波長が４０４．５５ｎｍである輝線成分、波長が４３５．８４ｎｍである輝線成分、波長が５０８．５８ｎｍである輝線成分、波長が５４６．０７ｎｍである輝線成分、波長が５７８ｎｍである輝線成分および波長が６４７．８５ｎｍである輝線成分を含む
   請求項６または７の分光装置を校正する方法。
9. 工程b)は、センサーの基準分光感度を光学シミュレーションにより求め、
   工程d)は、前記分光装置の機械的誤差を示す指標の単位量のずれに対するセンサーの分光感度を示す指標のずれ量を光学シミュレーションにより求め、工程d)において作成されたモデルに含まれる前記分光装置の機械的誤差を示す指標の１次関数において前記分光装置の機械的誤差を示す指標に求めた単位誤差量当たりのずれ量を加減算する
   請求項１から８までのいずれかの分光装置を校正する方法。
10. h) 分光装置を準備する工程と、
    i) 請求項１から９までのいずれかの分光装置を校正する方法により工程h)において準備された分光装置を校正する工程と、
    を備える校正された分光装置を生産する方法。
    

Images