(57)【要約】【課題】分光装置に備えられるセンサーの分光感度を精度よく適切に求める。【解決手段】分光装置は、複数の波長成分の各々のエネルギー量を示す信号を出力するセンサーを備える。複数のセンサーの各々について、センサーの基準分光感度が求められ、センサーの基準分光感度を示す指標を求められ、センサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれが分光装置の機械的誤差を示す指標の１次関数により表される。センサーの分光感度の集合がセンサーにより出力される信号の集合に適合するように分光装置の機械的誤差を示す指標が求められる。複数のセンサーの各々について、センサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれが求められ、センサーの分光感度が求められる。【選択図】図１
2020139963.tif 000002
【特許請求の範囲】
【請求項１】
a) 被測定光をスペクトルに変換する光学系と、スペクトルに含まれる複数の波長成分の各々のエネルギー量を示す信号を出力するセンサーを備えることにより複数の信号を出力する複数のセンサーを備える受光センサーと、を備える分光装置を校正の対象にする工程と、
b) 前記複数のセンサーの各々について、センサーの基準分光感度を求める工程と、
c) 前記複数のセンサーの各々について、工程b)において求められたセンサーの基準分光感度を示す指標を求める工程と、
d) 前記複数のセンサーの各々について、工程c)において求められたセンサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれが前記分光装置の機械的誤差を示す指標の１次関数により表されるモデルを作成する工程と、
e) 工程d)において作成されたモデルに含まれるセンサーの分光感度を示す指標により示されるセンサーの分光感度がセンサーにより出力される信号に適合するように前記分光装置の機械的誤差を示す指標を求める工程と、
f) 前記複数のセンサーの各々について、工程d)において作成されたモデルおよび工程e)において求められた前記分光装置の機械的誤差を示す指標を用いて、センサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれを求める工程と、
g) 前記複数のセンサーの各々について、工程b)において求められたセンサーの基準分光感度および工程f)において求められたセンサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれを用いて、センサーの分光感度を求める工程と、
を備える分光装置を校正する方法。
【請求項２】
工程c)において求められたセンサーの基準分光感度を示す指標は、センサーの基準分光感度の中心波長および半値幅であり、
工程d)において作成されたモデルに含まれるセンサーの分光感度を示す指標は、センサーの分光感度の中心波長および半値幅であり、
工程d)は、前記複数のセンサーの各々について、センサーの基準分光感度の中心波長からのセンサーの分光感度の中心波長のずれが前記分光装置の機械的誤差を示す指標の第１の１次関数により表される第１の式を作成し、センサーの基準分光感度の半値幅からのセンサーの分光感度の半値幅のずれが前記分光装置の機械的誤差を示す指標の第２の１次関数により表される第２の式を作成し、
工程d)において作成されるモデルは、第１の式および第２の式からなる
請求項１の分光装置を校正する方法。
【請求項３】
前記受光センサーは、受光面を有し、
前記複数のセンサーは、前記受光面において第１の方向に配列され、
前記光学系は、第２の方向に伸び前記受光面に至る光軸を有し、
工程d)において作成されたモデルに含まれる前記分光装置の機械的誤差を示す指標は、前記第１の方向についての前記受光センサーの配置誤差、前記第２の方向についての前記受光センサーの配置誤差および前記第１の方向と前記受光面に至る光軸とに直交する軸の周りを周る回転方向についての前記受光センサーの配置誤差の全部または一部を含む
請求項１または２の分光装置を校正する方法。
【請求項４】
前記光学系は、
スリットが形成されたスリット板と、
主断面および回折面を有する回折格子と、
を備え、
前記光学系は、前記スリットから前記回折面に至る光軸を有し、
工程d)において作成されたモデルに含まれる前記分光装置の機械的誤差を示す指標は、前記主断面と平行をなし前記スリットから前記回折面に至る光軸と垂直をなす方向についての前記スリットの幅の製作誤差を含む
請求項１から３までのいずれかの分光装置を校正する方法。
【請求項５】
前記光学系は、
主断面を有する回折格子
を備え、
工程d)において作成されたモデルに含まれる前記分光装置の機械的誤差を示す指標は、前記主断面が同一平面上に維持される回転方向についての前記回折格子の配置誤差を含む
請求項１から４までのいずれかの分光装置を校正する方法。
【請求項６】
工程e)は、
複数の輝線成分を含む輝線校正用の被測定光を前記分光装置に測定させ、
前記複数の輝線成分の各々について、輝線成分の波長において感度を有する第１のセンサーおよび第２のセンサーを選択し、第１のセンサーの分光感度の輝線成分の波長における感度である第１の感度、第２のセンサーの分光感度の輝線成分の波長における感度である第２の感度、第１のセンサーにより出力される信号である第１の信号、第２のセンサーにより出力される信号である第２の信号を用いて、第１の信号に対する第２の信号の比からの第１の感度に対する第２の感度の比のずれの前記複数の輝線成分についての合計が最小となるように前記分光装置の機械的誤差を示す指標を求めることにより、分光感度の集合を信号の集合に適合させる
請求項１から５までのいずれかの分光装置を校正する方法。
【請求項７】
工程e)は、
複数の輝線成分を含む輝線校正用の被測定光を前記分光装置に測定させ、
前記複数の輝線成分の各々について、輝線成分の波長において感度を有する第１のセンサーおよび第２のセンサーを選択し、第１のセンサーの分光感度の輝線成分の波長における感度である第１の感度、第２のセンサーの分光感度の輝線成分の波長における感度である第２の感度、第１のセンサーにより出力される信号である第１の信号、第２のセンサーにより出力される信号である第２の信号を用いて、第１の信号からの第１の感度のずれの絶対値が大きくなるほど絶対値が大きくなる因子と第２の信号からの第２の感度のずれの絶対値が大きくなるほど絶対値が大きくなる因子との和の前記複数の輝線成分についての合計が最小になるように前記分光装置の機械的誤差を示す指標を求めることにより、分光感度の集合を信号の集合に適合させる
請求項１から５までのいずれかの分光装置を校正する方法。
【請求項８】
前記輝線校正用の被測定光が、ＨｇＣｄランプにより放射され、波長が４０４．５５ｎｍである輝線成分、波長が４３５．８４ｎｍである輝線成分、波長が５０８．５８ｎｍである輝線成分、波長が５４６．０７ｎｍである輝線成分、波長が５７８ｎｍである輝線成分および波長が６４７．８５ｎｍである輝線成分を含む
請求項６または７の分光装置を校正する方法。
【請求項９】
工程b)は、センサーの基準分光感度を光学シミュレーションにより求め、
工程d)は、前記分光装置の機械的誤差を示す指標の単位量のずれに対するセンサーの分光感度を示す指標のずれ量を光学シミュレーションにより求め、工程d)において作成されたモデルに含まれる前記分光装置の機械的誤差を示す指標の１次関数において前記分光装置の機械的誤差を示す指標に求めた単位誤差量当たりのずれ量を加減算する
請求項１から８までのいずれかの分光装置を校正する方法。
【請求項１０】
h) 分光装置を準備する工程と、
i) 請求項１から９までのいずれかの分光装置を校正する方法により工程h)において準備された分光装置を校正する工程と、
を備える校正された分光装置を生産する方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、分光装置の校正に関する。
【背景技術】
【０００２】
ポリクロメーターを備える分光装置は、複数の波長成分の各々のエネルギー量を示す信号を出力するセンサーを備える。複数のセンサーにより出力される複数の信号からは、分光反射率、分光放射輝度、分光透過率等の分光特性が求められる。
【０００３】
センサーの分光感度は、受光センサーの配置、形状、大きさ等により変化し、ポリクロメーターに備えられるスリット板、回折格子等の配置、形状、大きさ等により変化する。このため、分光特性の基礎となるセンサーの分光感度は、分光装置ごとに校正されなければならない。
【０００４】
特許文献１に記載された分光装置の波長校正方法においては、分光感度の中心波長の基準値からの差が画素番号の１次関数で与えられ（段落００２６）、分光感度の半値幅の基準値に対する比が画素番号の１次関数で与えられ（段落００２８）、分光感度の中心波長および分光感度の半値幅の基準値に対する比を用いて修正分光感度テーブルが作成され（段落００２９）、修正分光感度テーブルおよび輝線波長から算出される算出相対出力と実測相対出力とが互いに最も近い値となるように画素番号の１次関数の各々に含まれる係数が決められる（段落００３３）。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特許第４６６０６９４号明細書
【発明の概要】
【発明が解決しようとする課題】
【０００６】
特許文献１に記載された分光装置の波長校正方法においては、分光感度の中心波長の基準値からの差が画素番号の１次関数で与えられ、分光感度の半値幅の基準値に対する比が画素番号の１次関数で与えられる。しかし、実際の分光装置においては、当該差が画素番号の１次関数とは著しく異なる場合または当該比が画素番号の１次関数とは著しく異なる場合がある。これらの場合は、当該差を画素番号の１次関数で与え当該比を画素番号の１次関数で与えることは良好な近似ではなく、分光感度が精度よく求められない。当該差を画素番号の２次関数以上の高次関数で与え当該比を画素番号の２次関数以上の高次関数で与えることも考えられるが、高次関数が用いられた場合は、係数が修正された場合に算出相対出力が大きく変化し算出相対出力と実測相対出力とが互いに最も近い値となる係数が適切に求められず、分光感度が適切に求められない。
【０００７】
下記の発明は、この問題を解決するためになされる。下記の発明が解決しようとする課題は、分光感度を精度よく適切に求めることである。
【課題を解決するための手段】
【０００８】
校正の対象にされる分光装置は、光学系および受光センサーを備える。光学系は、被測定光をスペクトルに変換する。受光センサーは、スペクトルに含まれる複数の波長成分の各々のエネルギー量を示す信号を出力するセンサーを備えることにより複数の信号を出力する複数のセンサーを備える。
【０００９】
複数のセンサーの各々について、センサーの基準分光感度が求められ、求められたセンサーの基準分光感度を示す指標が求められ、求められたセンサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれが分光装置の機械的誤差を示す指標の１次関数により表されるモデルが作成される。センサーの分光感度を示す指標により示されるセンサーの分光感度がセンサーにより出力される信号に適合するように分光装置の機械的誤差を示す指標が求められる。複数のセンサーの各々について、作成されたモデルおよび求められた分光装置の機械的誤差を示す指標を用いて、センサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれが求められ、求められたセンサーの基準分光感度およびセンサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれを用いて、センサーの分光感度が求められる。
【発明の効果】
【００１０】
センサーの分光感度を示す指標が分光装置の機械的誤差を示す指標の１次関数により表されるため、分光装置の機械的誤差を示す指標を変化させた場合に分光感度の集合の信号の集合への適合性が大きく変化しない。このため、分光装置の機械的誤差を示す指標が適切に求められ、センサーの分光感度が精度よく適切に求められる。
【図面の簡単な説明】
【００１１】
【図１】分光装置および校正装置を示す模式図である。
【図２】分光ユニットを示す模式図である。
【図３】リニアアレイセンサーを示す模式図である。
【図４】信号処理機構を示すブロック図である。
【図５】センサーの分光感度および輝線スペクトルの例を示すグラフである。
【図６】基準分光感度と校正後の分光感度との関係を示す模式図である。
【図７】基準分光感度の中心波長と中心波長のずれ量との関係を示すグラフである。
【図８】基準分光感度の中心波長と半値幅比との関係を示すグラフである。
【図９】第３のモデルを作成する手順を示すフローチャートである。
【図１０】機械誤差a1に起因するセンサーの各チャンネルの中心波長のずれ量を示すグラフである。
【図１１】機械誤差a1に起因するセンサーの各チャンネルの半値幅のずれ量を示すグラフである。
【図１２】機械誤差a2に起因するセンサーの各チャンネルの中心波長のずれ量を示すグラフである。
【図１３】機械誤差a2に起因するセンサーの各チャンネルの半値幅のずれ量を示すグラフである。
【図１４】機械誤差a3に起因するセンサーの各チャンネルの中心波長のずれ量を示すグラフである。
【図１５】機械誤差a3に起因するセンサーの各チャンネルの半値幅のずれ量を示すグラフである。
【図１６】機械誤差a4に起因するセンサーの各チャンネルの中心波長のずれ量を示すグラフである。
【図１７】機械誤差a4に起因するセンサーの各チャンネルの半値幅のずれ量を示すグラフである。
【図１８】機械誤差a5に起因するセンサーの各チャンネルの中心波長のずれ量を示すグラフである。
【図１９】機械誤差a5に起因するセンサーの各チャンネルの半値幅のずれ量を示すグラフである。
【図２０】分光装置の校正の手順を示すフローチャートである。
【図２１】センサーの分光感度および輝線スペクトルの例を示すグラフである。
【図２２】輝線光源測定時のセンサーのチャンネルと測定信号との関係を示すグラフである。
【図２３】校正された分光装置を生産する手順を示すフローチャートである。
【発明を実施するための形態】
【００１２】
１  分光装置
図１の模式図は、分光装置および校正装置を示す。図１には、分光装置の断面が示される。
【００１３】
分光装置１００は、図１に示されるように、分光ユニット１０４および信号処理機構１０５を備える。
【００１４】
分光装置１００が測定を行う場合は、分光ユニット１０４が、被測定光１０８を受光し、受光した被測定光１０８に含まれる第１の波長成分、第２の波長成分、・・・、第４０の波長成分のエネルギー量をそれぞれ示す第１の信号、第２の信号、・・・、第４０の信号を出力する。信号処理機構１０５は、第１の信号、第２の信号、・・・、第４０の信号の信号を用いて分光特性を求める。
【００１５】
分光装置１００が分光測色計である場合は、測定される物体に光が照射された場合に測定される物体により反射される光が被測定光１０８になり、分光反射率などが分光特性として求められる。分光装置１００が分光輝度計である場合は、測定される光源により放射される光が被測定光１０８になり、分光放射輝度などが分光特性として求められる。測定された物体に光が照射される場合に測定される物体を透過する光が被測定光１０８になってもよく、分光透過率などが分光特性として求められてもよい。測色値が分光特性として求められてもよい。測色値は、マンセル表色系、Ｌ＊ａ＊ｂ＊表色系、Ｌ＊Ｃ＊ｈ表色系、ハンターＬａｂ表色系、ＸＹＺ表色系等で表現される。
【００１６】
２  分光ユニット
図２の模式図は、分光ユニット１０４を示す斜視図である。
【００１７】
分光ユニット１０４は、図１および図２の各々に示されるように、光学系１１１およびリニアアレイセンサー１１２を備える。光学系１１１は、ポリクロメーターであり、スリット板１１５および凹面回折格子１１６を備える。
【００１８】
図３の模式図は、リニアアレイセンサーを示す平面図である。
【００１９】
リニアアレイセンサー１１２は、図３に示されるように、センサー１１９−１，１１９−２，・・・，１１９−４０を備える。
【００２０】
分光装置１００が測定を行う場合は、被測定光１０８がスリット板１１５に形成された矩形のスリット１２２に導かれる。スリット１２２に導かれた被測定光１０８は、スリット１２２を通過する。スリット１２２を通過した被測定光１０８は、スリット１２２から凹面回折格子１１６の回折面１２４まで進み、回折面１２４に反射される。被測定光１０８は、回折面１２４に反射されることにより、スペクトルに変換される。スペクトルに変換された被測定光１０８は、凹面回折格子１１６の回折面１２４からリニアアレイセンサー１１２の受光面１２５まで進み、受光面１２５において結像し、リニアアレイセンサー１１２に受光される。センサー１１９−１，１１９−２，・・・，１１９−４０は受光面１２５において方向１２８に直線状に配列される。スリット板１１５、凹面回折格子１１６およびリニアアレイセンサー１１２は、受光面１２５に結像する光の波長が方向１２８の位置に応じて変化するように配置される。このため、リニアアレイセンサー１１２がスペクトルに変換された被測定光１０８を受光した場合は、センサー１１９−１，１１９−２，・・・，１１９−４０が、互いに異なる第１の波長成分、第２の波長成分、・・・、第４０の波長成分のエネルギー量をそれぞれ示す第１の信号、第２の信号、・・・、第４０の信号をそれぞれ出力する。出力された第１の信号、第２の信号、・・・、第４０の信号は、信号処理機構１０５に入力される。信号処理機構１０５は、入力された第１の信号、第２の信号、・・・、第４０の信号を用いて分光特性を求める。
【００２１】
光学系１１１は、スリット１２２から回折面１２４に至る光軸１３４および回折面１２４から受光面１２５に至る光軸１３５を有する。
【００２２】
光学系１１１が他の種類の光学系に置き換えられてもよい。例えば、凹面回折格子１１６が平面回折格子および凹面鏡に置き換えられてもよい。また、スリット板１１５および凹面回折格子１１６が、円形のスリットが形成されたスリット板、シリンドリカルレンズおよびリニアバリアブルフィルターに置き換えられてもよい。円形のスリットを通過した被測定光は、シリンドリカルレンズを通過する。被測定光がシリンドリカルレンズを通過することにより、被測定光の断面形状が円形状から線状に変換される。シリンドリカルレンズを通過した被測定光は、リニアバリアブルフィルターを通過する。被測定光は、リニアバリアブルフィルターを通過することにより、スペクトルに変換される。
【００２３】
リニアアレイセンサー１１２が他の種類の受光センサーに置き換えられてもよい。例えば、リニアアレイセンサー１１２が、３９個以下または４１個以上のセンサーを備えるリニアアレイセンサーに置き換えられてもよい。光学系によっては、リニアアレイセンサー１１２がエリアセンサーに置き換えられてもよい。
【００２４】
３  信号処理機構
図４のブロック図は、信号処理機構を示す。
【００２５】
信号処理機構１０５は、図４に示されるように、Ａ／Ｄ変換機構１３８および演算機構１３９を備える。
【００２６】
第１の信号、第２の信号、・・・、第４０の信号が信号処理機構１０５に入力された場合は、第１の信号、第２の信号、・・・、第４０の信号がＡ／Ｄ変換機構１３８に入力される。Ａ／Ｄ変換機構１３８に入力された第１の信号、第２の信号、・・・、第４０の信号は、第１の信号値、第２の信号値、・・・、第４０の信号値にそれぞれアナログ／デジタル変換される。第１の信号値、第２の信号値、・・・、第４０の信号値は、演算機構１３９に入力される。演算機構１３９は、入力された第１の信号値、第２の信号値、・・・、第４０の信号値および演算機構１３９に格納されたセンサー１１９−１の分光感度１４０−１、センサー１１９−２の分光感度１４０−２、・・・、センサー１１９−４０の分光感度１４０−４０を用いて分光特性を求める。分光感度１４０−１，１４０−２，・・・，１４０−４０に代えて、分光感度１４０−１，１４０−２，・・・，１４０−４０から導かれる情報であって分光特性を求めるために必要な情報が演算機構１３９に格納されてもよい。
【００２７】
演算機構１３９は、組み込みコンピューターであり、インストールされたプログラムにしたがって動作する。演算機構１３９が行う処理の全部または一部がプログラムを伴わない電子回路により行われてもよい。演算機構１３９が行う処理の全部または一部が手作業で行われてもよい。
【００２８】
４  センサーの分光感度の例
図５のグラフは、センサーの分光感度および輝線校正用の被測定光の輝線スペクトルの例を示す。
【００２９】
センサー１１９−１の分光感度１４１−１の中心波長、センサー１１９−２の分光感度１４１−２の中心波長、・・・、センサー１１９−４０の分光感度１４１−４０の中心波長は、図５に示されるように、互いに異なり、それぞれ約３５２ｎｍ、約３６３ｎｍ、・・・、約７４０ｎｍである。これにより、センサー１１９−１，１１９−２，・・・，１１９−４０は、互いに異なる第１の波長成分、第２の波長成分、・・・、第４０の波長成分のエネルギー量をそれぞれ示す第１の信号、第２の信号、・・・、第４０の信号をそれぞれ出力する。
【００３０】
５  分光装置の校正の必要性
分光感度１４１−１，１４１−２，・・・，１４１−４０の各々は、スリット板１１５、凹面回折格子１１６およびリニアアレイセンサー１１２の配置、形状、大きさ等により変化する。このため、分光特性が精度よく求められるためには、演算機構１３９に格納された分光感度１４０−１，１４０−２，・・・，１４０−４０を、スリット板１１５、凹面回折格子１１６およびリニアアレイセンサーの配置、形状、大きさ等に応じて変更し、真の分光感度１４１−１，１４１−２，・・・，１４１−４０にそれぞれ近づけなければならない。以下では、演算機構１３９に格納された分光感度１４０−１，１４０−２，・・・，１４０−４０を真の分光感度１４１−１，１４１−２，・・・，１４１−４０にそれぞれ近づけることを分光装置１００の校正という。
【００３１】
６  校正装置
校正装置１４２は、図１に示されるように、ＨｇＣｄランプ１４３および制御演算機構１４４を備える。
【００３２】
校正装置１４２を用いて分光装置１００の校正が行われる場合は、制御演算機構１４４がＨｇＣｄランプ１４３に被測定光１０８として輝線校正用の被測定光を放射させる。放射された輝線校正用の被測定光は、分光装置１００により測定される。輝線校正用の被測定光が分光装置１００により測定された場合は、第１の信号値、第２の信号値、・・・、第４０の信号値が演算機構１３９に入力される。入力された第１の信号値、第２の信号値、・・・、第４０の信号値は、演算機構１３９から制御演算機構１４４へ転送される。制御演算機構１４４は、転送されてきた第１の信号値、第２の信号値、・・・、第４０の信号値を用いてセンサー１１９−１の分光感度、センサー１１９−２の分光感度、・・・、センサー１１９−４０の分光感度を求める。求められたセンサー１１９−１の分光感度、センサー１１９−２の分光感度、・・・、センサー１１９−４０の分光感度は、制御演算機構１４４から演算機構１３９へ転送され、それぞれ演算機構１３９に新たに格納される分光感度１４０−１，１４０−２，・・・，１４０−３になる。これにより、分光装置１００の校正が行われた後は、演算機構１３９が第１の信号、第２の信号、・・・、第４０の信号および新たに求められたセンサー１１９−１の分光感度、センサー１１９−２の分光感度、・・・、センサー１１９−４０の分光感度を用いて分光特性を求めることができる。
【００３３】
制御演算機構１４４は、コンピューターであり、インストールされたプログラムにしたがって動作する。制御演算機構１４４が行う処理の全部または一部がプログラムを伴わない電子回路により行われてもよい。制御演算機構１４４が行う処理の全部または一部が手作業で行われてもよい。制御演算機構１４４が分光装置１００に内蔵されてもよい。
【００３４】
７  輝線成分の波長
輝線校正用の被測定光は、図５に示されるように、輝線成分１４５−１，１４５−２，１４５−３，１４５−４，１４５−５および１４５−６を含む。輝線成分１４５−１，１４５−２，１４５−３，１４５−４，１４５−５および１４５−６の波長は、それぞれ４０４．５５ｎｍ，４３５．８４ｎｍ，５０８．５８ｎｍ，５４６．０７ｎｍ，５７８ｎｍおよび６４７．８５ｎｍである。輝線成分１４５−１，１４５−２，１４５−３，１４５−４，１４５−５および１４５−６は、分光装置１００の校正に用いられる。
【００３５】
輝線成分１４５−１，１４５−２，１４５−３，１４５−４，１４５−５および１４５−６以外の輝線成分が分光装置１００の校正に用いられてもよい。５個以下または７個以上の輝線成分が分光装置１００の校正に用いられてもよい。ＨｇＣｄランプ１４３以外の輝線光源から放射される光が輝線校正用の被測定光として用いられてもよい。輝線校正でない校正が行われ、輝線光源でない光源から放射される光が校正用の被測定光として用いられる場合もある。
【００３６】
８  分光装置の校正において用いられるモデル
分光装置１００の校正においては、第１のモデル、第２のモデルまたは第３のモデルが用いられる。第１のモデルおよび第２のモデルは、参考例である。
【００３７】
９  センサーの分光感度を示す指標
第１のモデル、第２のモデルおよび第３のモデルの各々においては、センサー１１９−１，１１９−２，・・・，１１９−４０が位置ｉにより識別される。位置ｉは、互いに異なる４０個の値ｉ１，ｉ２，・・・，ｉ４０のいずれかをとる。センサー１１９−１，１１９−２，・・・，１１９−４０が位置以外の指標により識別されてもよい。例えば、センサー１１９−１，１１９−２，・・・，１１９−４０が基準分光感度の中心波長、画素番号等により識別されてもよい。
【００３８】
第１のモデル、第２のモデルおよび第３のモデルの各々においては、位置ｉにあるセンサーの分光感度が、位置ｉにあるセンサーの分光感度の中心波長λＧ（ｉ）および半値幅ＦＷＨＭ（ｉ）により示される。中心波長λＧ（ｉ）および半値幅ＦＷＨＭ（ｉ）の各々は、位置ｉの関数である。
【００３９】
位置ｉにあるセンサーの分光感度は、独立変数が波長であり従属変数が感度であるガウス関数により良好に近似される。当該ガウス関数の形状は、中心波長および半値幅により決定される。このため、中心波長λＧ（ｉ）および半値幅ＦＷＨＭ（ｉ）は、位置ｉにあるセンサーの分光感度を示す指標として好適である。ただし、位置ｉにあるセンサーの分光感度を示す指標が変更されてもよい。位置ｉにあるセンサーの分光感度を示す指標が１個の変数または３個以上の変数からなることも許される。
【００４０】
１０  第１のモデル（参考例）
第１のモデルが用いられる場合は、中心波長λＧ（ｉ）が位置ｉのｎ次関数で表される式（１）が作成される。
【００４１】
【数１】
2020139963.tif 000003
【００４２】
また、半値幅ＦＷＨＭ（ｉ）が位置ｉのｍ次関数により表される式（２）が作成される。
【００４３】
【数２】
2020139963.tif 000004
【００４４】
第１のモデルは、式（１）および（２）からなる。係数ａｎ，・・・，ａ０は、説明変数になり、位置ｉにあるセンサーの分光感度を決める校正パラメーターになる。係数ｂｎ，・・・，ｂ０は、説明変数になり、位置ｉにあるセンサーの分光感度を決める校正パラメーターになる。
【００４５】
第１のモデルを用いて分光装置１００の校正が行われる場合は、中心波長λＧ（ｉ）および半値幅ＦＷＨＭ（ｉ）により示されるセンサーの分光感度の２個以上のセンサーについての集合である分光感度の集合が、センサーにより出力される信号の２個以上のセンサーについての集合である信号の集合に適合するように、係数ａｎ，・・・，ａ０が求められ、係数ｂｎ，・・・，ｂ０が求められる。
【００４６】
続いて、作成された式（１）および求められた係数ａｎ，・・・，ａ０を用いて中心波長λＧ（ｉ）が求められ、作成された式（２）および求められた係数ｂｎ，・・・，ｂ０を用いて半値幅ＦＷＨＭ（ｉ）が求められる。
【００４７】
続いて、求められた中心波長λＧ（ｉ）および半値幅ＦＷＨＭ（ｉ）を用いて位置ｉにあるセンサーの分光感度が求められる。位置ｉにあるセンサーの分光感度は、求められた中心波長λＧ（ｉ）および半値幅ＦＷＨＭ（ｉ）により示される分光感度である。
【００４８】
第１のモデルを用いて分光装置１００の校正が行われた場合、輝線成分の波長λＨｇＣｄ（１），λＨｇＣｄ（２）、・・・，λＨｇＣｄ（Ｋ０）のいずれかに近い波長領域においては分光感度が精度よく求められるが、それ以外の波長領域においては分光感度が精度よく求められない。特に、最も短波長側にある波長領域および最も長波長側にある波長領域においては、分光感度が精度よく求められない。
【００４９】
１１  第２のモデル（参考例）
第２のモデルは、第１のモデルの問題を解決するために提供される。
【００５０】
第２のモデルが用いられる場合は、スリット板１１５、凹面回折格子１１６およびリニアアレイセンサー１１２の配置、形状、大きさ等が設計狙いに一致する理想的な分光装置１００が想定され、想定された分光装置１００が備える位置ｉにあるセンサーの基準分光感度が光学シミュレーションにより求められる。位置ｉにあるセンサーの基準分光感度は、位置ｉの関数である。
【００５１】
続いて、位置ｉにあるセンサーの基準分光感度の中心波長λＧ０（ｉ）および半値幅ＦＷＨＭ０（ｉ）が求められる。
【００５２】
中心波長λＧ０（ｉ）および半値幅ＦＷＨＭ０（ｉ）は、位置ｉにあるセンサーの基準分光感度を示す指標として好適である。
【００５３】
続いて、中心波長λＧ０（ｉ）からの中心波長λＧ（ｉ）のずれΔλＧ（ｉ）が位置ｉの１次関数で表される式（３）が作成される。
【００５４】
【数３】
2020139963.tif 000005
【００５５】
また、半値幅ＦＷＨＭ０（ｉ）に対する半値幅ＦＷＨＭ（ｉ）の比ｒａｔｉｏ（ｉ）が位置ｉの１次関数で表される式（４）が作成される。
【００５６】
【数４】
2020139963.tif 000006
【００５７】
第２のモデルは、式（３）および（４）からなる。係数ａ１およびａ０は、説明変数になり、位置ｉにあるセンサーの分光感度を決める校正パラメーターになる。係数ｂ１およびｂ０は、説明変数になり、位置ｉにあるセンサーの分光感度を決める校正パラメーターになる。
【００５８】
中心波長λＧ（ｉ）は、式（３）から導かれる式（５）を用いて求められる。
【００５９】
【数５】
2020139963.tif 000007
【００６０】
第２のモデルを用いて分光装置１００の校正が行われる場合は、中心波長λＧ（ｉ）および半値幅ＦＷＨＭ（ｉ）により示されるセンサーの分光感度の２個以上のセンサーについての集合である分光感度の集合が、センサーにより出力される信号の２個以上のセンサーについての集合である信号の集合に適合するように、係数ａ１およびａ０が求められ、係数ｂ１およびｂ０が求められる。
【００６１】
続いて、作成された式（３）ならびに求められた係数ａ１およびａ０を用いてずれΔλＧ（ｉ）が求められ、作成された式（４）ならびに求められた係数ｂ１およびｂ０を用いて比ｒａｔｉｏ（ｉ）が求められる。
【００６２】
続いて、求められたずれΔλＧ（ｉ）および比ｒａｔｉｏ（ｉ）を用いて位置ｉにあるセンサーの分光感度が求められる。
【００６３】
図６の模式図は、位置ｉにあるセンサーの基準分光感度と位置ｉにあるセンサーの分光感度との関係を示す。
【００６４】
位置ｉにあるセンサーの分光感度は、図６に示されるように、位置ｉにあるセンサーの基準分光感度を中心波長λＧ０（ｉ）を中心として波長軸方向にｒａｔｉｏ（ｉ）倍に拡大し、拡大された基準分光感度を波長軸方向にΔλＧ（ｉ）だけ移動させたものである。
【００６５】
図７のグラフは、実際の分光装置１００における基準分光感度の中心波長λＧ０とずれΔλＧとの関係を示す実測データである。図８のグラフは、実際の分光装置１００における基準分光感度の中心波長λＧ０と比ｒａｔｉｏとの関係を示す実測データである。
【００６６】
実際の分光装置１００においては、図７に示されるように、ずれΔλＧが基準分光感度の中心波長λＧ０の１次関数といいがたい場合があり、図８に示されるように、比ｒａｔｉｏが中心波長λＧ０の１次関数とはいいがたい場合がある。このような場合は、式（３）および（４）の各々は良好な近似ではない。このため、第２のモデルを用いて分光装置１００の校正が行われた場合は、分光感度が精度よく求められない場合がある。
【００６７】
この問題を解決するためにずれΔλＧ（ｉ）が位置ｉの２次関数以上の高次関数で表され比ｒａｔｉｏ（ｉ）が位置ｉの２次関数以上の高次関数で表された場合は、係数ａ１，ａ０，ｂ１およびｂ０が変化した場合に分光感度の集合の信号の集合への適合性が大きく変化する。このため、係数ａ１，ａ０，ｂ１およびｂ０が適切に求められない場合がある。
【００６８】
１２  第３のモデル
１２．１  第３のモデルの作成
第３のモデルは、第２のモデルの問題を解決するために提供される。
【００６９】
図９のフローチャートは、第３のモデルを作成する手順を示す。
【００７０】
複数の分光装置の各々について分光装置１００の校正が行われる場合でも、複数の分光装置の各々について図９のステップS1およびS2に示される手順を実行する必要はなく、図９に示されるステップS1およびS2の手順をあらかじめ１回だけ実行すればよい。ただし、複数の分光装置の各々について分光装置１００の校正が行われる間に図９に示される手順が再実行され、第３のモデルが再作成されてもよい。
【００７１】
第３のモデルが作成される場合は、ステップＳ１において、スリット板１１５、凹面回折格子１１６およびリニアアレイセンサー１１２の配置、形状、大きさ等が設計狙いに一致する理想的な分光装置１００が想定され、想定された分光装置１００が備える位置ｉにあるセンサーの基準分光感度が光学シミュレーションにより求められる。位置ｉにあるセンサーの基準分光感度は、位置ｉの関数である。このため、ステップＳ１により、位置ｉ１にあるセンサー、位置ｉ２にあるセンサー、・・・、位置ｉ４０にあるセンサーの各々について、センサーの基準分光感度が求められる。
【００７２】
続いて、ステップＳ２において、位置ｉにあるセンサーの基準分光感度の中心波長λＧ０（ｉ）および半値幅ＦＷＨＭ０（ｉ）が求められる。中心波長λＧ０（ｉ）および半値幅ＦＷＨＭ０（ｉ）の各々は、位置ｉの関数である。このため、ステップＳ２により、位置ｉ１にあるセンサー、位置ｉ２にあるセンサー、・・・、位置ｉ４０にあるセンサーの各々について、センサーの基準分光感度の中心波長および半値幅が求められる。
【００７３】
中心波長λＧ０（ｉ）および半値幅ＦＷＨＭ０（ｉ）は、位置ｉにあるセンサーの基準分光感度を示す指標として好適である。
【００７４】
続いて、ステップＳ３において、中心波長λＧ０（ｉ）からの中心波長λＧ（ｉ）のずれがリニアアレイセンサー１１２の配置誤差ａ１，ａ２およびａ３、スリット１２２の幅の製作誤差ａ４ならびに凹面回折格子１１６の配置誤差ａ５の１次関数により表される式（６）が作成される。また、半値幅ＦＷＨＭ０（ｉ）からの半値幅ＦＷＨＭ（ｉ）のずれがリニアアレイセンサー１１２の配置誤差ａ１，ａ２およびａ３、スリット１２２の幅の製作誤差ａ４ならびに凹面回折格子１１６の配置誤差ａ５の１次関数により表される式（７）が作成される。中心波長λＧ０（ｉ）および半値幅ＦＷＨＭ０（ｉ）の各々は、位置ｉの関数である。誤差ａ１，ａ２，ａ３，ａ４およびａ５は、分光装置１００の機械的誤差を示す指標である。このため、ステップＳ３により、位置ｉ１にあるセンサー、位置ｉ２にあるセンサー、・・・、位置ｉ４０にあるセンサーの各々について、センサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれが分光装置１００の機械的誤差を示す指標の１次関数により表される第３のモデルが作成される。
【００７５】
【数６】
2020139963.tif 000008
【００７６】
【数７】
2020139963.tif 000009
【００７７】
第３のモデルは、式（６）および（７）からなる。誤差ａ１，ａ２，ａ３，ａ４およびａ５は、説明変数になり、位置ｉにあるセンサーの分光感度を決める校正パラメーターになる。
【００７８】
リニアアレイセンサー１１２の配置誤差ａ１は、センサー１１９−１，１１９−２，・・・，１１９−４０が配列される方向１２８についてのリニアアレイセンサー１１２の配置誤差である。リニアアレイセンサー１１２の配置誤差ａ２は、光軸１３５が伸びる方向１２９についてのリニアアレイセンサーの配置誤差である。リニアアレイセンサー１１２の配置誤差ａ３は、方向１２８と光軸１３５とに直交するセンサーの軸１３０ｃ周り１３０の回転誤差である。スリット１２２の幅の製作誤差ａ４は、凹面回折格子１１６の主断面１４６と平行をなし光軸１３４と垂直をなす方向１３１についてのスリット１２２の幅の製作誤差である。配置誤差ａ５は、凹面回折格子１１６の軸１３２ｃ周り１３２の回転誤差である。凹面回折格子１１６の主断面１４６は、回折面１２４に形成される刻線と垂直をなす面である。
【００７９】
誤差ａ１，ａ２，ａ３，ａ４およびａ５は小さいことが期待されるため、中心波長λＧ０（ｉ）からの中心波長λＧ（ｉ）のずれを誤差ａ１，ａ２，ａ３，ａ４およびａ５の１次関数により表すことは良好な近似であり、半値幅ＦＷＨＭ０（ｉ）からの半値幅ＦＷＨＭ（ｉ）のずれを誤差ａ１，ａ２，ａ３，ａ４およびａ５の１次関数により表すことも良好な近似である。このため、第３のモデルを用いて分光装置１００の校正が行われた場合は、分光感度が精度よく求められる。
【００８０】
誤差ａ１，ａ２，ａ３，ａ４およびａ５は、中心波長λＧ（ｉ）または半値幅ＦＷＨＭ（ｉ）への影響が大きく、位置ｉにあるセンサーの分光感度への影響が大きい。このため、誤差ａ１，ａ２，ａ３，ａ４およびａ５は、分光装置１００の機械的誤差を示す指標として好適である。ただし、分光装置１００の機械的誤差を示す指標が変更されてもよい。分光装置１００の機械的誤差を示す指標に含まれる説明変数または校正パラメーターの数は、輝線校正用の被測定光に含まれる輝線成分の数以下に制限されるが、４個以下または６個以上であってもよい。
【００８１】
式（６）に含まれる誤差ａ１，ａ２，ａ３，ａ４およびａ５の１次関数においては、誤差ａ１，ａ２，ａ３，ａ４およびａ５に、係数δλＧ１（ｉ），δλＧ２（ｉ），δλＧ３（ｉ），δλＧ４（ｉ）およびδλＧ５（ｉ）がそれぞれ乗じられる。係数δλＧ１（ｉ），δλＧ２（ｉ），δλＧ３（ｉ），δλＧ４（ｉ）およびδλＧ５（ｉ）は、それぞれ誤差ａ１，ａ２，ａ３，ａ４およびａ５が単位量発生した場合の中心波長λＧ（ｉ）のずれ量であり、光学シミュレーションにより求められる。係数δλＧ１（ｉ），δλＧ２（ｉ），δλＧ３（ｉ），δλＧ４（ｉ）およびδλＧ５（ｉ）の各々は、位置ｉの関数である。
【００８２】
式（７）に含まれる誤差ａ１，ａ２，ａ３，ａ４およびａ５の１次関数においては、誤差ａ１，ａ２，ａ３，ａ４およびａ５に、係数δＦＷＨＭ１（ｉ），δＦＷＨＭ２（ｉ），δＦＷＨＭ３（ｉ），δＦＷＨＭ４（ｉ）およびδＦＷＨＭ５（ｉ）がそれぞれ乗じられる（加減算される）。係数δＦＷＨＭ１（ｉ），δＦＷＨＭ２（ｉ），δＦＷＨＭ３（ｉ），δＦＷＨＭ４（ｉ）およびδＦＷＨＭ５（ｉ）は、それぞれ誤差ａ１，ａ２，ａ３，ａ４およびａ５が単位量発生した場合の半値幅ＦＷＨＭ（ｉ）のずれ量であり、光学シミュレーションにより求められる。係数δＦＷＨＭ１（ｉ），δＦＷＨＭ２（ｉ），δＦＷＨＭ３（ｉ），δＦＷＨＭ４（ｉ）およびδＦＷＨＭ５（ｉ）の各々は、位置ｉの関数である。
【００８３】
図１０から図１９までの各々におけるチャンネルＣｈは、センサー１１９−１，１１９−２，・・・，１１９−４０の識別番号である。
【００８４】
図１０の１４９は、配置誤差ａ１が０から＋１単位ずれた場合の各Ｃｈの分光感度の中心波長のずれ量δλＧを示し、１５０は配置誤差ａ１が０から−１単位ずれた場合の各Ｃｈの分光感度の中心波長のずれ量δλＧを示す。関係１４９および１５０は、光学シミュレーションにより求められる。係数δλＧ１（ｉ）は、例えば、関係１４９および１５０を用いて求められる。
【００８５】
図１１の１５３は、配置誤差ａ１が０から＋１単位ずれた場合のセンサーの各Ｃｈの分光感度の半値幅のずれ量δＦＷＨＭを示し、１５４は、配置誤差ａ１が０から−１単位ずれた場合のセンサーの各Ｃｈの分光感度の半値幅のずれ量δＦＷＨＭを示す。関係１５３および１５４は、光学シミュレーションにより求められる。係数δＦＷＨＭ１（ｉ）は、例えば、関係１５３および１５４を用いて求められる。
【００８６】
図１２の１５７は、配置誤差ａ２が０から＋１単位ずれた場合の各Ｃｈの分光感度の中心波長のずれ量δλＧを示し、１５８は配置誤差ａ２が０から−１単位ずれた場合の各Ｃｈの分光感度の中心波長のずれ量δλＧを示す。関係１５７および１５８は、光学シミュレーションにより求められる。係数δλＧ２（ｉ）は、例えば、関係１５７および１５８を用いて求められる。
【００８７】
図１３の１６１は、配置誤差ａ２が０から＋１単位ずれた場合のセンサーの各Ｃｈの分光感度の半値幅のずれ量δＦＷＨＭを示し、１６２は配置誤差ａ２が０から−１単位ずれた場合のセンサーの各Ｃｈの分光感度の半値幅のずれ量δＦＷＨＭを示す。関係１６１および１６２は、光学シミュレーションにより求められる。係数δＦＷＨＭ２（ｉ）は、例えば、関係１６１および１６２を用いて求められる。
【００８８】
図１４の１６５は、配置誤差ａ３が０から＋１単位ずれた場合の各Ｃｈの分光感度の中心波長のずれ量δλＧを示し、１６６は配置誤差ａ３が０から−１単位ずれた場合の各Ｃｈの分光感度の中心波長のずれ量δλＧを示す。関係１６５および１６６は、光学シミュレーションにより求められる。係数δλＧ３（ｉ）は、例えば、関係１６５および１６６を用いて求められる。
【００８９】
図１５の１６９は、配置誤差ａ３が０から＋１単位ずれた場合のセンサーの各Ｃｈの分光感度の半値幅のずれ量δＦＷＨＭを示し、１７０は配置誤差ａ３が０から−１単位ずれた場合のセンサーの各Ｃｈの分光感度の半値幅のずれ量δＦＷＨＭを示す。関係１６９および１７０は、光学シミュレーションにより求められる。係数δＦＷＨＭ３（ｉ）は、例えば、関係１６９および１７０を用いて求められる。
【００９０】
図１６の１７３は、製作誤差ａ４が０から＋１単位ずれた場合の各Ｃｈの分光感度の中心波長のずれ量δλＧを示し、１７４は製作誤差ａ４が０から−１単位ずれた場合の各Ｃｈの分光感度の中心波長のずれ量δλＧを示す。関係１７３および１７４は、光学シミュレーションにより求められる。係数δλＧ４（ｉ）は、例えば、関係１７３および１７４を用いて求められる。
【００９１】
図１７の１７７は、製作誤差ａ４が０から＋１単位ずれた場合のセンサーの各Ｃｈの分光感度の半値幅のずれ量δＦＷＨＭのチャンネルＣｈを示し、１７８は製作誤差ａ４が０から−１単位ずれた場合のセンサーの各Ｃｈの分光感度の半値幅のずれ量δＦＷＨＭを示す。関係１７７および１７８は、光学シミュレーションにより求められる。係数δＦＷＨＭ４（ｉ）は、例えば、関係１７７および１７８を用いて求められる。
【００９２】
図１８の１８１は、配置誤差ａ５が０から＋１単位ずれた場合の各Ｃｈの分光感度の中心波長のずれ量δλＧを示し、１８２は配置誤差ａ５が０から−１単位ずれた場合の各Ｃｈの分光感度の中心波長のずれ量δλＧを示す。関係１８１および１８２は、光学シミュレーションにより求められる。係数δλＧ５（ｉ）は、例えば、関係１８１および１８２を用いて求められる。
【００９３】
図１９の１８５は、配置誤差ａ５が０から＋１単位ずれた場合のセンサーの各Ｃｈの分光感度の半値幅のずれ量δＦＷＨＭを示し、１８６は配置誤差ａ５が０から−１単位ずれた場合のセンサーの各Ｃｈの分光感度の半値幅のずれ量δＦＷＨＭを示す。関係１８５および１８６は、光学シミュレーションにより求められる。係数δＦＷＨＭ５（ｉ）は、例えば、関係１８５および１８６を用いて求められる。
【００９４】
１２．２  第３のモデルを用いた分光装置の校正
図２０のフローチャートは、第３のモデルを用いた分光装置の校正の手順を示す。
【００９５】
複数の分光装置の各々について分光装置１００の校正が行われる場合は、複数の分光装置の各々について図２０に示される手順が実行される。
【００９６】
第３のモデルを用いて分光装置１００の校正が行われる場合は、ステップＳ１１において、中心波長λＧ（ｉ）および半値幅ＦＷＨＭ（ｉ）により示されるセンサーの分光感度が、センサーにより出力される信号に適合するように、誤差ａ１，ａ２，ａ３，ａ４およびａ５が求められる。当該分光感度を当該信号に適合させることは、位置ｉにあるセンサーの真の分光感度が中心波長λＧ（ｉ）および半値幅ＦＷＨＭ（ｉ）により示されるセンサーの分光感度である場合にセンサーにより出力されると想定される信号を、センサーにより実際に出力される信号に近づけることを意味する。適合性の評価には、目的変数が用いられる。ステップＳ１１により、センサーの分光感度を示す指標により示されるセンサーの分光感度が、センサーにより出力される信号に適合するように分光装置１００の機械的誤差を示す指標が求められる。
【００９７】
校正にはセンサー１１９−１，１１９−２，・・・，１１９−４０のうちの輝線成分の波長λＨｇＣｄ（１），λＨｇＣｄ（２），・・・，λＨｇＣｄ（Ｋ０）において感度を有するセンサーの出力を用いる。
【００９８】
続いて、ステップＳ１２において、作成された式（６）および求められた誤差ａ１，ａ２，ａ３，ａ４およびａ５を用いて、中心波長λＧ０（ｉ）からの中心波長λＧ（ｉ）のずれが求められる。また、作成された式（７）および求められた誤差ａ１，ａ２，ａ３，ａ４およびａ５を用いて、半値幅ＦＷＨＭ０（ｉ）からの半値幅ＦＷＨＭ（ｉ）のずれが求められる。ステップＳ１２により、位置ｉ１のセンサー、位置ｉ２のセンサー、・・・、位置ｉ４０のセンサーの各々について、作成された第３のモデルおよび求められた分光装置１００の機械的誤差を示す指標を用いて、センサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれが求められる。
【００９９】
続いて、ステップＳ１３において、求められた位置ｉにあるセンサーの基準分光感度、中心波長λＧ０（ｉ）からの中心波長λＧ（ｉ）のずれおよび半値幅ＦＷＨＭ０（ｉ）からの半値幅ＦＷＨＭ（ｉ）のずれを用いて、位置ｉにあるセンサーの分光感度が求められる。これにより、位置ｉ１にあるセンサー、位置ｉ２にあるセンサー、・・・、位置ｉ４０にあるセンサーの各々について、求められたセンサーの基準分光感度およびセンサーの基準分光感度を示す指標からのセンサーの分光感度を示す指標のずれを用いてセンサーの分光感度が求められる。
【０１００】
求められる位置ｉにあるセンサーの分光感度は、図６に示されるように、基準分光感度を中心波長λＧ０（ｉ）を中心として波長軸方向にｒａｔｉｏ（ｉ）倍に拡大し、拡大された基準分光感度を波長軸方向にΔλＧ（ｉ）だけ移動させたものである。
【０１０１】
比ｒａｔｉｏ（ｉ）は、式（８）により表される。
【０１０２】
【数８】
2020139963.tif 000010
【０１０３】
第３のモデルを用いて分光装置１００の校正が行われた場合は、中心波長λＧ０（ｉ）および半値幅ＦＷＨＭ０（ｉ）の各々が誤差ａ１，ａ２，ａ３，ａ４およびａ５の１次関数により表されるため、誤差ａ１，ａ２，ａ３，ａ４およびａ５を変化させた場合に分光感度の信号への適合性が大きく変化しない。このため、誤差ａ１，ａ２，ａ３，ａ４およびａ５が適切に求められ、センサーの校正後の分光感度が適切に求められる。
【０１０４】
また、第３のモデルを用いて分光装置１００の校正が行われた場合は、中心波長λＧ０（ｉ）および半値幅ＦＷＨＭ０（ｉ）の各々が共通の説明変数である誤差ａ１，ａ２，ａ３，ａ４およびａ５で表されるため、中心波長と半値幅を別の説明変数で表す場合に比べて説明変数の数が減少し、分光装置１００の校正に必要な輝線成分の数が減少する。
【０１０５】
さらに、第３のモデルを用いて分光装置１００の校正が行われた場合は、中心波長λＧ０（ｉ）および半値幅ＦＷＨＭ０（ｉ）が共通の説明変数である誤差ａ１，ａ２，ａ３，ａ４およびａ５で表されるため、中心波長λＧ０（ｉ）および半値幅ＦＷＨＭ０（ｉ）の関係が不適切なものにならず、中心波長λＧ０（ｉ）および半値幅ＦＷＨＭ０（ｉ）が適切に求められる。
【０１０６】
１３  分光感度を信号に適合させる手順
以下では、波長がλＨｇＣｄ（ｋ）である輝線成分が位置Ｉｋにあるセンサーおよび位置Ｉｋ＋１にあるセンサーにまたがって入射するとする。位置Ｉｋ＋１にあるセンサーの基準分光感度の中心波長は、位置Ｉｋにあるセンサーの基準分光感度の中心波長に隣接する。輝線成分の識別番号ｋは、１，２，・・・，Ｋ０のいずれかの値をとるとする。
【０１０７】
図２１のグラフは、センサーの分光感度および輝線校正用の被測定光の輝線スペクトルの例を示す。図２１のグラフは、図５のグラフの波長が３９０ｎｍから４２０ｎｍまでである範囲を拡大したものである。
【０１０８】
波長が４０４．５４ｎｍである輝線成分１４５−１は、中心波長が約３９６ｎｍである分光感度１４１−５を有するセンサーおよび中心波長が約４０７ｎｍである分光感度１４１−６を有するセンサーにまたがって入射する。このため、分光感度１４１−５を有するセンサーおよび分光感度１４１−６を有するセンサーの各々は、図２１に示されるように、輝線成分１４５−１に対して感度を有する。
【０１０９】
図２２のグラフは、センサーのチャンネルとセンサーにより出力される信号との関係を示す。
【０１１０】
チャンネル５のセンサーおよびチャンネル６のセンサーの各々は輝線成分１４５−１に対して感度を有し、チャンネル８のセンサーおよびチャンネル９のセンサーの各々は輝線成分１４５−２に対して感度を有し、チャンネル１５およびチャンネル１６のセンサーの各々は輝線成分１４５−３に対して感度を有し、チャンネル１９のセンサーおよびチャンネル２０のセンサーの各々は輝線成分１４５−４に対して感度を有し、チャンネル２２のセンサーおよびチャンネル２３のセンサーは輝線成分１４５−５に対して感度を有し、チャンネル２９のセンサーおよびチャンネル３０のセンサーの各々は輝線成分１４５−６に対して感度を有する。その結果、図２２に示されるような信号が得られる。
【０１１１】
波長がλＨｇＣｄ（１）である輝線成分、波長がλＨｇＣｄ（２）である輝線成分、・・・、波長がλＨｇＣｄ（Ｋ０）である輝線成分を用いて分光装置１００の校正が完全に行われた場合は、位置ｉにあるセンサーの分光感度Ｒｅｓｐｏｎｓｅ（ｉ，λ）および位置ｉにあるセンサーにより出力される信号をアナログ／デジタル変換することにより得られる信号値Ｃｏｕｎｔ（ｉ）は、ｋ＝１，２，・・・，Ｋ０の各々について、式（９）および（１０）に示される関係を満たす。
【０１１２】
【数９】
2020139963.tif 000011
【０１１３】
【数１０】
2020139963.tif 000012
【０１１４】
式（９）および（１０）は、ｋ＝１，２，・・・，Ｋ０の各々について、波長λＨｇＣｄ（ｋ）において感度を有する位置Ｉｋにあるセンサーおよび位置Ｉｋ＋１にあるセンサーを選択した場合に、位置Ｉｋにあるセンサーの分光感度の波長λＨｇＣｄ（ｋ）における感度Ｒｅｓｐｏｎｓｅ（Ｉｋ，λＨｇＣｄ（ｋ））が位置Ｉｋにあるセンサーにより出力される信号を示す信号値Ｃｏｕｎｔ（Ｉｋ）に一致し、位置Ｉｋ＋１にあるセンサーの分光感度の波長λＨｇＣｄ（ｋ）における感度Ｒｅｓｐｏｎｓｅ（Ｉｋ＋１，λＨｇＣｄ（ｋ））が位置Ｉｋ＋１にあるセンサーにより出力される信号を示す信号値Ｃｏｕｎｔ（Ｉｋ＋１）に一致することを示す。
【０１１５】
分光感度の集合を信号の集合に適合させることは、感度Ｒｅｓｐｏｎｓｅ（Ｉｋ，λＨｇＣｄ（ｋ））および信号値Ｃｏｕｎｔ（Ｉｋ）の関係を式（９）に示される関係に近づけ、感度Ｒｅｓｐｏｎｓｅ（Ｉｋ＋１，λＨｇＣｄ（ｋ））および信号値Ｃｏｕｎｔ（Ｉｋ＋１）の関係を式（１０）に示される関係に近づけることを意味する。
【０１１６】
このため、分光感度を信号に適合させる第１の方法は、式（１１）に示される目的変数Ｆを最小にする誤差ａ１，ａ２，ａ３，ａ４およびａ５を求めることである。
【０１１７】
【数１１】
2020139963.tif 000013
【０１１８】
式（１１）に示される目的変数Ｆは、信号値Ｃｏｕｎｔ（Ｉｋ）からの感度Ｒｅｓｐｏｎｓｅ（Ｉｋ，λＨｇＣｄ（ｋ））のずれの２乗と信号値Ｃｏｕｎｔ（Ｉｋ＋１）からの感度Ｒｅｓｐｏｎｓｅ（Ｉｋ＋１，λＨｇＣｄ（ｋ））のずれの２乗との和の、ｋ＝１，２，・・・，Ｋ０についての合計である。
【０１１９】
ずれの２乗が、ずれの絶対値が大きくなるほど絶対値が大きくなる他の因子に置き換えられてもよい。例えば、ずれの２乗がずれの絶対値に置き換えられてもよい。
【０１２０】
もしある輝線波長λHgCd(k)に対して感度を有するセンサーが３つ(Ik-1, Ik, Ik+1)ある場合には信号Count(Ik-1)と感度Response(Ik-1)とのずれ量を目的関数に加えてもよい。
【０１２１】
また、式（９）および（１０）からは式（１２）が導かれる。
【０１２２】
【数１２】
2020139963.tif 000014
【０１２３】
分光感度を信号に適合させることは、感度Ｒｅｓｐｏｎｓｅ（Ｉｋ，λＨｇＣｄ（ｋ））およびＲｅｓｐｏｎｓｅ（Ｉｋ＋１，λＨｇＣｄ（ｋ））ならびに信号値Ｃｏｕｎｔ（Ｉｋ）およびＣｏｕｎｔ（Ｉｋ＋１）の関係を式（１２）に示される関係に近づけることを意味する。
【０１２４】
このため、分光感度を信号に適合させる第２の方法は、式（１３）に示される目的変数Ｆを最小にする誤差ａ１，ａ２，ａ３，ａ４およびａ５を求めることである。
【０１２５】
【数１３】
2020139963.tif 000015
【０１２６】
式（１３）に示される目的関数Ｆは、信号値Ｃｏｕｎｔ（Ｉｋ）に対する信号値Ｃｏｕｎｔ（Ｉｋ＋１）の比からの感度Ｒｅｓｐｏｎｓｅ（Ｉｋ，λＨｇＣｄ（ｋ））に対する感度Ｒｅｓｐｏｎｓｅ（Ｉｋ＋１，λＨｇＣｄ（ｋ））の比のずれの、ｋ＝１，２，・・・，Ｋ０についての合計である。
【０１２７】
式（１３）に示される目的関数Ｆによれば、感度Ｒｅｓｐｏｎｓｅ（Ｉｋ，λＨｇＣｄ（ｋ））およびＲｅｓｐｏｎｓｅ（Ｉｋ＋１，λＨｇＣｄ（ｋ））を信号値Ｃｏｕｎｔ（Ｉｋ）およびＣｏｕｎｔ（Ｉｋ＋１）とそれぞれ対比できるようにするために分光感度Ｒｅｓｐｏｎｓｅ（ｉ，λ）を規格化する必要がない。
【０１２８】
もしある輝線波長λHgCd(k)に対して感度を有するセンサーが３つ(Ik-1, Ik, Ik+1)ある場合には、信号Count(Ik)とCount(Ik+1)との比と感度Response(Ik)とResponse(Ik+1)との比のずれ量に加えて、信号Count(Ik)とCount(Ik-1)との比と感度Response(Ik)とResponse(Ik-1)との比のずれ量とのずれ量を目的関数に加えてもよい。
【０１２９】
１４  校正された分光装置の生産
校正された分光装置の生産においては、図２３に示されるように、ステップＳ２１において分光装置１００が準備され、ステップＳ２２において準備された分光装置１００が校正される。ステップＳ２１においては、分光装置１００の校正を行う事業者が分光装置１００を製作することにより分光装置１００が準備されてもよいし、分光装置１００の校正を行う事業者が他の事業者から分光装置１００を仕入れることにより分光装置１００が準備されてもよい。
【符号の説明】
【０１３０】
１００  分光装置
１１１  光学系
１１２  リニアアレイセンサー
１１５  スリット板
１１６  凹面回折格子
１１９−１，１１９−２，・・・，１１９−４０  センサー
１２２  スリット
１４３  ＨｇＣｄランプ
【図１】
2020139963.tif 000016
【図２】
2020139963.tif 000017
【図３】
2020139963.tif 000018
【図４】
2020139963.tif 000019
【図５】
2020139963.tif 000020
【図６】
2020139963.tif 000021
【図７】
2020139963.tif 000022
【図８】
2020139963.tif 000023
【図９】
2020139963.tif 000024
【図１０】
2020139963.tif 000025
【図１１】
2020139963.tif 000026
【図１２】
2020139963.tif 000027
【図１３】
2020139963.tif 000028
【図１４】
2020139963.tif 000029
【図１５】
2020139963.tif 000030
【図１６】
2020139963.tif 000031
【図１７】
2020139963.tif 000032
【図１８】
2020139963.tif 000033
【図１９】
2020139963.tif 000034
【図２０】
2020139963.tif 000035
【図２１】
2020139963.tif 000036
【図２２】
2020139963.tif 000037
【図２３】
2020139963.tif 000038