(57)【特許請求の範囲】
【請求項１】
周波数、又は波長、又は波数、又はエネルギーのいずれかのＸに対して、吸収強度Ｙがプロットされたスペクトル（Ｘ，Ｙ）において、Ｙの値が有効であるＸの有効範囲を決定し、前記スペクトル（Ｘ，Ｙ）のうち、前記有効範囲内のスペクトルを解析するスペクトル解析方法であって、
前記吸収強度Ｙ、前記吸収強度Ｙの一次微分値Ｙ'、及び前記吸収強度Ｙの二次微分値Ｙ"のすべてを用い、それぞれに対して決められた規定値を用いて、前記有効範囲を決定することを特徴とするスペクトル解析方法。
【請求項２】
前記吸収強度Ｙまたは前記吸収強度Ｙの絶対値｜Ｙ｜のいずれかが、前記Ｙの規定値を超えるデータポイント（Ｘ，Ｙ）におけるＸを、上限値Ｘmaxまたは下限値Ｘminとし、前記上限値Ｘmaxまたは下限値Ｘminのいずれかから規定の間隔だけ離れたＸに対する吸収強度Ｙまたは吸収強度Ｙの絶対値のいずれかが、前記Ｙの規定値を超えない側の範囲を第1の有効範囲と判定し、
前記一次微分Ｙ'または前記一次微分Ｙ'の絶対値のいずれかが、前記Ｙ'の規定値を超えるデータポイント（Ｘ，Ｙ'）におけるＸを、上限値Ｘmaxまたは下限値Ｘminとし、前記上限値Ｘmaxまたは下限値Ｘminのいずれかから規定の間隔だけ離れたＸに対する一次微分Ｙ'または一次微分Ｙ'の絶対値のいずれかが、前記一次微分Ｙ'の規定値を超えない側の範囲を第２の有効範囲と判定し、
前記二次微分Ｙ"または前記二次微分Ｙ"の絶対値のいずれかが、前記Ｙ"の規定値を超えるデータポイント（Ｘ，Ｙ"）におけるＸを、上限値Ｘmaxまたは下限値Ｘminとし、前記上限値Ｘmaxまたは下限値Ｘminのいずれかから規定の間隔だけ離れたＸに対する二次微分Ｙ"または二次微分Ｙ"の絶対値のいずれかが、前記二次微分Ｙ"の規定値を超えない側の範囲を第３の有効範囲と判定し、
前記第1の有効範囲と前記第２の有効範囲と前記第３の有効範囲の重なった範囲を、前記有効範囲と決定することを特徴とする請求項１記載のスペクトル解析方法。
【請求項３】
前記吸収強度Ｙまたは前記吸収強度Ｙの絶対値｜Ｙ｜のいずれかが、前記Ｙの規定値を超えるデータポイント（Ｘ１，Ｙ１）、（Ｘ２，Ｙ２）、（Ｘ３，Ｙ３）、及び（Ｘ４，Ｙ４）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、規定の間隔より狭い場合で、かつ前記Ｘ１より小さい側に前記規定の間隔だけ離れたＸにおいて、Ｙまたは｜Ｙ｜のいずれかがＹlimを超えない場合、
又は、
前記一次微分Ｙ'または前記一次微分Ｙ'の絶対値｜Ｙ'｜のいずれかが、前記Ｙ'の規定値を超えるデータポイント（Ｘ１，Ｙ１'）、（Ｘ２，Ｙ２'）、（Ｘ３，Ｙ３'）、及び（Ｘ４，Ｙ４'）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、規定の間隔より狭い場合で、かつ前記Ｘ１より小さい側に前記規定の間隔だけ離れたＸにおいて、Ｙ'または｜Ｙ'｜のいずれかがＹ'limを超えない場合、
又は、
前記二次微分Ｙ"または前記二次微分Ｙ"の絶対値｜Ｙ"｜のいずれかが、前記Ｙ"の規定値を超えるデータポイント（Ｘ１，Ｙ１"）、（Ｘ２，Ｙ２"）、（Ｘ３，Ｙ３"）、及び（Ｘ４，Ｙ４"）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、規定の間隔より狭い場合で、かつ前記Ｘ１より小さい側に前記規定の間隔だけ離れたＸにおいて、Ｙ"または｜Ｙ"｜のいずれかがＹ"limを超えない場合、
のいずれかのＸ１を前記有効範囲の上限値として決定することを特徴とする請求項１記載のスペクトル解析方法。
【請求項４】
前記吸収強度Ｙまたは前記吸収強度Ｙの絶対値｜Ｙ｜のいずれかが、前記Ｙの規定値を超えるデータポイント（Ｘ１，Ｙ１）、（Ｘ２，Ｙ２）、（Ｘ３，Ｙ３）、及び（Ｘ４，Ｙ４）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、規定の間隔より狭い場合で、かつ前記Ｘ４より大きい側に前記規定の間隔だけ離れたＸにおいて、Ｙまたは｜Ｙ｜のいずれかがＹlimを超えない場合、
又は、
前記一次微分Ｙ'または前記一次微分Ｙ'の絶対値｜Ｙ'｜のいずれかが、前記Ｙ'の規定値を超えるデータポイント（Ｘ１，Ｙ１'）、（Ｘ２，Ｙ２'）、（Ｘ３，Ｙ３'）、及び（Ｘ４，Ｙ４'）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、規定の間隔より狭い場合で、かつ前記Ｘ４より大きい側に前記規定の間隔だけ離れたＸにおいて、Ｙ'または｜Ｙ'｜のいずれかがＹ'limを超えない場合、
又は、
前記二次微分Ｙ"または前記二次微分Ｙ"の絶対値｜Ｙ"｜のいずれかが、前記Ｙ"の規定値を超えるデータポイント（Ｘ１，Ｙ１"）、（Ｘ２，Ｙ２"）、（Ｘ３，Ｙ３"）、及び（Ｘ４，Ｙ４"）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、規定の間隔より狭い場合で、かつ前記Ｘ４より大きい側に前記規定の間隔だけ離れたＸにおいて、Ｙ"または｜Ｙ"｜のいずれかがＹ"limを超えない場合、
のいずれかのＸ４を前記有効範囲の下限値として決定することを特徴とする請求項１記載のスペクトル解析方法。
【請求項５】
前記吸収強度Ｙまたは前記吸収強度Ｙの絶対値｜Ｙ｜のいずれかが、前記Ｙの規定値を超えるデータポイント（Ｘ１，Ｙ１）、（Ｘ２，Ｙ２）、（Ｘ３，Ｙ３）、及び（Ｘ４，Ｙ４）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、規定の間隔より狭い場合で、かつ前記Ｘ１より小さい側に前記規定の間隔だけ離れたＸにおいて、Ｙまたは｜Ｙ｜のいずれかがＹlimを超えない場合、
又は、
前記一次微分Ｙ'または前記一次微分Ｙ'の絶対値｜Ｙ'｜のいずれかが、前記Ｙ'の規定値を超えるデータポイント（Ｘ１，Ｙ１'）、（Ｘ２，Ｙ２'）、（Ｘ３，Ｙ３'）、及び（Ｘ４，Ｙ４'）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、規定の間隔より狭い場合で、かつ前記Ｘ１より小さい側に前記規定の間隔だけ離れたＸにおいて、Ｙ'または｜Ｙ'｜のいずれかがＹ'limを超えない場合、
又は、
前記二次微分Ｙ"または前記二次微分Ｙ"の絶対値｜Ｙ"｜のいずれかが、前記Ｙ"の規定値を超えるデータポイント（Ｘ１，Ｙ１"）、（Ｘ２，Ｙ２"）、（Ｘ３，Ｙ３"）、及び（Ｘ４，Ｙ４"）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、規定の間隔より狭い場合で、かつ前記Ｘ１より小さい側に前記規定の間隔だけ離れたＸにおいて、Ｙ"または｜Ｙ"｜のいずれかがＹ"limを超えない場合、
のいずれかのＸ１のうち、最小のＸ１を前記有効範囲の上限値とし、
前記吸収強度Ｙまたは前記吸収強度Ｙの絶対値｜Ｙ｜のいずれかが、前記Ｙの規定値を超えるデータポイント（Ｘ１，Ｙ１）、（Ｘ２，Ｙ２）、（Ｘ３，Ｙ３）、及び（Ｘ４，Ｙ４）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、規定の間隔より狭い場合で、かつ前記Ｘ４より大きい側に前記規定の間隔だけ離れたＸにおいて、Ｙまたは｜Ｙ｜のいずれかがＹlimを超えない場合、
又は、
前記一次微分Ｙ'または前記一次微分Ｙ'の絶対値｜Ｙ'｜のいずれかが、前記Ｙ'の規定値を超えるデータポイント（Ｘ１，Ｙ１'）、（Ｘ２，Ｙ２'）、（Ｘ３，Ｙ３'）、及び（Ｘ４，Ｙ４'）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、規定の間隔より狭い場合で、かつ前記Ｘ４より大きい側に前記規定の間隔だけ離れたＸにおいて、Ｙ'または｜Ｙ'｜のいずれかがＹ'limを超えない場合、
又は、
前記二次微分Ｙ"または前記二次微分Ｙ"の絶対値｜Ｙ"｜のいずれかが、前記Ｙ"の規定値を超えるデータポイント（Ｘ１，Ｙ１"）、（Ｘ２，Ｙ２"）、（Ｘ３，Ｙ３"）、及び（Ｘ４，Ｙ４"）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、規定の間隔より狭い場合で、かつ前記Ｘ４より大きい側に前記規定の間隔だけ離れたＸにおいて、Ｙ"または｜Ｙ"｜のいずれかがＹ"limを超えない場合、
のいずれかのＸ４のうち、最大のＸ４を前記有効範囲の下限値として決定することを特徴とする請求項１記載のスペクトル解析方法。
【請求項６】
前記周波数が０．１[THz]〜１０[THz]、前記波長が３[μｍ]〜３０[μｍ]、前記波数が３．３[cm−１]〜３３３．３[cm−１]、前記エネルギーが４．１３５７×１０−４［ｅＶ］〜４．１３５７×１０−２［ｅＶ］の範囲に適用されることを特徴とする請求項１から請求項５のうちいずれか１項に記載されたスペクトル解析方法。
【請求項７】
コンピュータを、
周波数、又は波長、又は波数、又はエネルギーのいずれかのＸに対して、吸収強度Ｙがプロットされたスペクトル（Ｘ，Ｙ）において、Ｙの値が有効であるＸの範囲を有効範囲と決定する手段として動作させるプログラムであって、
前記吸収強度Ｙまたは前記吸収強度Ｙの絶対値のいずれかが、Ｙの規定値を超えるデータポイント（Ｘ，Ｙ）におけるＸを、上限値Ｘmaxまたは下限値Ｘminとし、前記上限値Ｘmaxまたは下限値Ｘminのいずれかから規定の間隔だけ離れたＸに対する吸収強度Ｙまたは吸収強度Ｙの絶対値のいずれかが、前記Ｙの規定値を超えない側の範囲を第1の有効範囲と判定し、
前記一次微分Ｙ'または前記一次微分Ｙ'の絶対値のいずれかが、Ｙ'の規定値を超えるデータポイント（Ｘ，Ｙ'）におけるＸを、上限値Ｘmaxまたは下限値Ｘminとし、前記上限値Ｘmaxまたは下限値Ｘminのいずれかから規定の間隔だけ離れたＸに対する一次微分Ｙ'または一次微分Ｙ'の絶対値のいずれかが、前記一次微分Ｙ'の規定値を超えない側の範囲を第２の有効範囲と判定し、
前記二次微分Ｙ"または前記二次微分Ｙ"の絶対値のいずれかが、Ｙ"の規定値を超えるデータポイント（Ｘ，Ｙ"）におけるＸを、上限値Ｘmaxまたは下限値Ｘminとし、前記上限値Ｘmaxまたは下限値Ｘminのいずれかから規定の間隔だけ離れたＸに対する二次微分Ｙ"または二次微分Ｙ"の絶対値のいずれかが、前記二次微分Ｙ"の規定値を超えない側の範囲を第３の有効範囲と判定し、
前記第１の有効範囲と前記第２の有効範囲と前記第３の有効範囲の重なった範囲を、前記有効範囲と決定することを特徴とするプログラム。
【請求項８】
コンピュータを、
周波数、又は波長、又は波数、又はエネルギーのいずれかのＸに対して、吸収強度Ｙがプロットされたスペクトル（Ｘ，Ｙ）において、Ｙの値が有効であるＸの範囲を有効範囲と決定する手段として動作させるプログラムであって、
前記吸収強度Ｙまたは前記吸収強度Ｙの絶対値｜Ｙ｜のいずれかが、吸収強度Ｙの規定値を超えるデータポイント（Ｘ１，Ｙ１）、（Ｘ２，Ｙ２）、（Ｘ３，Ｙ３）、及び（Ｘ４，Ｙ４）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、規定の間隔より狭い場合で、かつ前記Ｘ１より小さい側に前記規定の間隔だけ離れたＸにおいて、Ｙまたは｜Ｙ｜のいずれかがＹlimを超えない場合、
又は、
前記一次微分Ｙ'または前記一次微分Ｙ'の絶対値｜Ｙ'｜のいずれかが、Ｙ'の規定値を超えるデータポイント（Ｘ１，Ｙ１'）、（Ｘ２，Ｙ２'）、（Ｘ３，Ｙ３'）、及び（Ｘ４，Ｙ４'）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、規定の間隔より狭い場合で、かつ前記Ｘ１より小さい側に前記規定の間隔だけ離れたＸにおいて、Ｙ'または｜Ｙ'｜のいずれかがＹ'limを超えない場合、
又は、
前記二次微分Ｙ"または前記二次微分Ｙ"の絶対値｜Ｙ"｜のいずれかが、Ｙ"の規定値を超えるデータポイント（Ｘ１，Ｙ１"）、（Ｘ２，Ｙ２"）、（Ｘ３，Ｙ３"）、及び（Ｘ４，Ｙ４"）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、規定の間隔より狭い場合で、かつ前記Ｘ１より小さい側に前記規定の間隔だけ離れたＸにおいて、Ｙ"または｜Ｙ"｜のいずれかがＹ"limを超えない場合、
のいずれかのＸ１を前記有効範囲の上限値として決定することを特徴とするプログラム。
【請求項９】
コンピュータを、
周波数、又は波長、又は波数、又はエネルギーのいずれかのＸに対して、吸収強度Ｙがプロットされたスペクトル（Ｘ，Ｙ）において、Ｙの値が有効であるＸの範囲を有効範囲と決定する手段として動作させるプログラムであって、
前記吸収強度Ｙまたは前記吸収強度Ｙの絶対値｜Ｙ｜のいずれかが、Ｙの規定値を超えるデータポイント（Ｘ１，Ｙ１）、（Ｘ２，Ｙ２）、（Ｘ３，Ｙ３）、及び（Ｘ４，Ｙ４）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、規定の間隔より狭い場合で、かつ前記Ｘ４より大きい側に前記規定の間隔だけ離れたＸにおいて、Ｙまたは｜Ｙ｜のいずれかがＹlimを超えない場合、
又は、
前記一次微分Ｙ'または前記一次微分Ｙ'の絶対値｜Ｙ'｜のいずれかが、Ｙ'の規定値を超えるデータポイント（Ｘ１，Ｙ１'）、（Ｘ２，Ｙ２'）、（Ｘ３，Ｙ３'）、及び（Ｘ４，Ｙ４'）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、規定の間隔より狭い場合で、かつ前記Ｘ４より大きい側に前記規定の間隔だけ離れたＸにおいて、Ｙ'または｜Ｙ'｜のいずれかがＹ'limを超えない場合、
又は、
前記二次微分Ｙ"または前記二次微分Ｙ"の絶対値｜Ｙ"｜のいずれかが、Ｙ"の規定値を超えるデータポイント（Ｘ１，Ｙ１"）、（Ｘ２，Ｙ２"）、（Ｘ３，Ｙ３"）、及び（Ｘ４，Ｙ４"）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、規定の間隔より狭い場合で、かつ前記Ｘ４より大きい側に前記規定の間隔だけ離れたＸにおいて、Ｙ"または｜Ｙ"｜のいずれかがＹ"limを超えない場合、
のいずれかのＸ４を前記有効範囲の下限値として決定することを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、スペクトル解析において、有効なスペクトルデータ範囲を特定してスペクトル演算処理を行うスペクトル解析方法、及び、プログラムに関する。
【背景技術】
【０００２】
従来、周波数、波長、波数、エネルギーのいずれかに対して吸収強度（光学密度ともいう）がプロットされているスペクトルには、周波数、波長、波数のいずれかの範囲を予め限定して測定されるデータもあるが、例えば、時間領域分光法のように、時間領域のデータをフーリエ変換し、周波数、波長、波数のいずれかに対する吸収強度スペクトルを得る手法もある。後者の場合、予め有効な周波数、波長、波数の範囲を限定してデータの測定を行うことは困難である。
【０００３】
スペクトル測定時には、ノイズや装置由来が原因のアーティファクトなど、無効なデータを含んだ結果が得られる場合があるため、そのスペクトルを用いた成分分析や定量分析、他の装置で得られたデータとの比較、混合スペクトルの予想などで必要となるスペクトルの足し算、差分、定数倍、定量解析に代表されるスペクトル演算処理を行う上で、信頼性の低いスペクトル範囲が含まれてしまうという問題があった。このため、有効なスペクトル範囲を限定するなどのスペクトルの加工処理が必要となるが、従来では、大雑把にスペクトルの有効範囲を決定する手法（非特許文献１）や、信号と雑音との強度比から有効範囲を決定する客観的な手法（非特許文献２）などが用いられた。
【非特許文献１】田隅三生著「ＦＴ−ＩＲの基礎と実際」、東京化学同人（１９９４）、Ｐ２１９．
【非特許文献２】P.Uhd, B.M.Fischer, Opt.Lett.30(2005)29.
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかしながら、従来のスペクトル範囲の限定処理のうち、大雑把にスペクトルの有効範囲を決定する場合は、客観性に乏しい処理が行われる場合があるという問題があった。また、信号と雑音との強度比から有効範囲を決定する場合、強い吸収ピークを持つスペクトルの場合は、吸収ピークをノイズとして誤判定してしまう場合があるという問題があった。
本発明は、以上の点に鑑みてなされたものであり、客観的な手法で、かつ、鋭いピークを有するスペクトルにおいても、吸収ピークをノイズと誤判定することなく、スペクトルの足し算、差分、定数倍、定量解析に代表されるスペクトル演算処理を行う上で必要な、スペクトルの有効範囲を取り出すスペクトル解析手法、及び、そのプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００５】
前述した目的を達成するために本発明は、周波数、波長、波数、エネルギー（Ｘ）のいずれかに対して、吸収強度（Ｙ）がプロットされているスペクトル（Ｘ，Ｙ）において、Ｙの値が有効であるとみなせるＸの範囲［Ｘmin〜Ｘmax］を決定することで有効なスペクトルの最大範囲を限定することを特徴とし、この有効な最大範囲を超えないＸの範囲に限定した、一つないし複数のスペクトルを用いてスペクトルの足し算、差分、定数倍、定量解析に代表されるスペクトルの演算処理を行うことを可能とする。
本発明では、吸収強度Ｙ、そのＹの一次微分値Ｙ'、Ｙの二次微分値Ｙ"のすべてを用い、それぞれに対して決められた規定値を用いて、スペクトルの有効範囲を決定することを特徴とする。
【０００６】
本発明では、吸収強度ＹまたはＹの絶対値（｜Ｙ｜）のいずれかが、ある規定値（Ｙlim）を越えるデータポイント（Ｘ，Ｙ）におけるＸを限界値ＸminまたはＸmaxとし、限界値ＸminまたはＸmaxに隣接するＸにおいて、Ｙまたは｜Ｙ｜が、前記Ｙlimを超えない側のＸの範囲を有効範囲であると判定することが好ましい。
【０００７】
また、本発明では、吸収強度Ｙの一次微分値Ｙ'またはＹ'の絶対値（｜Ｙ'｜）のいずれかが、ある規定値（Ｙ'lim）を越えるデータポイント（Ｘ，Ｙ）におけるＸを限界値ＸminまたはＸmaxとし、限界値ＸminまたはＸmaxに隣接するＸにおいて、Ｙ'または｜Ｙ'｜が、前記Ｙ'limを超えない側のＸの範囲を有効範囲であると判定することが好ましい。
【０００８】
また、本発明では、吸収強度Ｙの二次微分値Ｙ"またはＹ"の絶対値（｜Ｙ"｜）のいずれかが、ある規定値（Ｙ"lim）を越えるデータポイント（Ｘ，Ｙ）におけるＸを限界値ＸminまたはＸmaxとし、限界値ＸminまたはＸmaxに隣接するＸにおいて、Ｙ"または｜Ｙ"｜が、前記Ｙ"limを超えない側のＸの範囲を有効範囲であると判定することが好ましい。
こうして、Ｙ、Ｙ'、Ｙ"を用いて決定された限界値Ｘminのうち最大のものをＸの有効範囲の下限値、Ｙ、Ｙ'、Ｙ"を用いて決定された限界値Ｘmaxのうち最小のものをＸの有効範囲の上限値として決定する、即ち、上述のＹ、Ｙ'、Ｙ"を用いて決定されたそれぞれの有効範囲の重なった部分をスペクトルの有効範囲と決定することを特徴とする。
【０００９】
また、本発明では、吸収強度Ｙ或いはその絶対値｜Ｙ｜のいずれかがＹlimを越える４つのデータポイント（Ｘ１，Ｙ１）、（Ｘ２，Ｙ２）、（Ｘ３，Ｙ３）及び（Ｘ４，Ｙ４）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４とする）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、ある規定の間隔（Ｘgap）より狭い場合で、かつＸ１より小さい側に隣接するＸにおいて、Ｙまたは｜Ｙ｜のいずれかがＹlimを超えない場合、Ｘ１＝ＸmaxとしてＸの有効範囲を決定することが好ましい。
【００１０】
また、本発明では、吸収強度の一次微分値Ｙ'或いはその絶対値｜Ｙ'｜のいずれかがＹ'limを越える４つのデータポイント（Ｘ１，Ｙ１）、（Ｘ２，Ｙ２）、（Ｘ３，Ｙ３）及び（Ｘ４，Ｙ４）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４とする）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、ある規定の間隔（Ｘgap）より狭い場合で、かつＸ１より小さい側に隣接するＸにおいて、Ｙ'または｜Ｙ'｜のいずれかがＹ'limを超えない場合、Ｘ１＝ＸmaxとしてＸの有効範囲を決定することが好ましい。
【００１１】
また、本発明では、吸収強度の二次微分値Ｙ"或いはその絶対値｜Ｙ"｜のいずれかがＹ"limを越える４つのデータポイント（Ｘ１，Ｙ１）、（Ｘ２，Ｙ２）、（Ｘ３，Ｙ３）及び（Ｘ４，Ｙ４）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４とする）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、ある規定の間隔（Ｘgap）より狭い場合で、かつＸ１より小さい側に隣接するＸにおいて、Ｙ"または｜Ｙ"｜のいずれかがＹ"limを超えない場合、Ｘ１＝ＸmaxとしてＸの有効範囲を決定することが好ましい。
【００１２】
また、本発明では、吸収強度Ｙ或いはその絶対値｜Ｙ｜のいずれかがＹlimを越える４つのデータポイント（Ｘ１，Ｙ１）、（Ｘ２，Ｙ２）、（Ｘ３，Ｙ３）及び（Ｘ４，Ｙ４）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４とする）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、ある規定の間隔（Ｘgap）より狭い場合で、かつＸ４より大きい側に隣接するＸにおいて、Ｙまたは｜Ｙ｜のいずれかがＹlimを超えない場合、Ｘ４＝ＸminとしてＸの有効範囲を決定することが好ましい。
【００１３】
また、本発明では、吸収強度の一次微分値Ｙ'またはその絶対値｜Ｙ'｜のいずれかが規定値Ｙ'limを越える４つのデータポイント（Ｘ１，Ｙ１）、（Ｘ２，Ｙ２）、（Ｘ３，Ｙ３）及び（Ｘ４，Ｙ４）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４とする）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、ある規定の間隔（Ｘgap）より狭い場合で、かつＸ４より大きい側に隣接するＸにおいて、Ｙ'または｜Ｙ'｜のいずれかがＹ'limを超えない場合、Ｘ４＝ＸminとしてＸの有効範囲を決定することが好ましい。
【００１４】
また、本発明の一つの形態として、吸収強度の二次微分値Ｙ"またはその絶対値｜Ｙ"｜のいずれかが規定値Ｙ"limを越える４つのデータポイント（Ｘ１，Ｙ１）、（Ｘ２，Ｙ２）、（Ｘ３，Ｙ３）及び（Ｘ４，Ｙ４）（ただし、Ｘ１＜Ｘ２＜Ｘ３＜Ｘ４とする）において、Ｘ１とＸ２の間隔（Ｘ２−Ｘ１）、Ｘ２とＸ３の間隔（Ｘ３−Ｘ２）、Ｘ３とＸ４の間隔（Ｘ４−Ｘ３）の全てが、ある規定の間隔（Ｘgap）より狭い場合で、かつＸ４より大きい側に隣接するＸにおいて、Ｙ"または｜Ｙ"｜のいずれかがＹ"limを超えない場合、Ｘ４＝ＸminとしてＸの有効範囲を決定することが好ましい。
【００１５】
また、Ｙ、Ｙ'、Ｙ"を用いて決定されたＸmax及びＸminのうち、最小のＸmaxと最大のＸminに囲まれたＸの範囲を、Ｘの有効範囲であると判定することが好ましい。
また、Ｘの有効範囲は、例えば、Ｘが周波数である場合には０．１[THz]〜１０[THz]の範囲に、Ｘが波長である場合には３［μｍ］〜３０［μｍ］の範囲に、Ｘが波数の場合には３．３［ｃｍ−１］〜３３３．３［ｃｍ−１］の範囲に、Ｘがエネルギーの場合には、４．１３５７×１０−４［ｅＶ］〜４．１３５７×１０−２［ｅＶ］の範囲に適用されることが好ましい。
【発明の効果】
【００１６】
本発明によれば、スペクトルの鋭いピークも誤判定することなく、有効なデータ範囲を特定し、その有効範囲のスペクトルを取り出し、スペクトルを用いた成分分析や定量分析、他の装置で得られたデータとの比較、混合スペクトルの予想などにおいて、信頼性の高い解析が可能となる。
【発明を実施するための最良の形態】
【００１７】
以下に、添付図面を参照しながら、本発明に係るスペクトル解析方法の好適な実施形態について詳細に説明する。
市販のテラヘルツ時間領域分光測定装置（例えば、先端赤外社製ＴＨｚ−ＴＤＳ２００４と同程度の仕様）を用いて、糖の粉末結晶とポリエチレン粉末を混合し、ペレット状にした濃度既知の試料のテラヘルツ透過スペクトルを測定する場合を例に説明する。テラヘルツ波の発生には、例えば、市販のフェムト秒レーザ（例えば、Coherent社製、Vitesse、800[nm]、100[fs]、80[MHz]、或いは、Femtolaser社製、Integral Pro、800[nm]、10[fs]、80[MHz]）を光源に用いる。前述のペレット状の試料は、一般的な赤外吸収スペクトルの測定用のＫＢｒペレットと同様に油圧式プレス機を用いて作成することが可能であり、１[mg]〜５０[mg]程度のアミノ酸粉末を５０[mg]〜２００[mg]程度のポリエチレン粉末と混合した濃度範囲が好ましい。また、ペレットの大きさとしては、直径５[mm]〜１５[mm]程度、厚さは１[mm]〜２[mm]程度のものが好ましい。
【００１８】
ペレット状の試料にテラヘルツ波を照射し、テラヘルツ時間領域分光測定装置により測定された時間領域テラヘルツ透過スペクトルをフーリエ変換することにより、周波数領域のスペクトルを得た場合を考える。図１は、得られた周波数領域のスペクトルを示す図である。図１に示すように、約０．５[THz]〜３[THz]の低周波数領域では、いくつかのピークが観測される。多くの糖類やアミノ酸の透過スペクトルは、図１に示すように、０．５[THz]〜３[THz]の低周波数領域にシャープで強度の高いピークを有し、その強度は濃度に比例して強くなることが知られている。これらは、定性定量分析のマーカーバンドとなる。これらのマーカーバンドのほか、高周波数領域に向かって吸収強度が徐々に増加する、ベースライン状のブロードな吸収が見られるが、これは、単純なベースラインではなく、分子振動モードに由来する多くの特性吸収が重ね合わさって観測されたものと考えられ、一般的にほとんどの物質のスペクトルで観測される。
【００１９】
従って、図１に示すようなテラヘルツスペクトルを用いた定性定量分析には、マーカーバンドのピーク位置と吸収強度だけではなく、このブロードな吸収も含めたスペクトル形状全体を含む標準スペクトルを基準として用いる必要がある。しかし、測定領域の上限周波数は、光源の励気パルス幅のほか、試料の吸収強度にも関係することが知られており、例えば約３[THz]より高い周波数領域では、試料の強い吸収に加えて、テラヘルツ光源のパワーが弱いことから、測定装置で測定された光のパワーが弱くなり、ノイズが大きくなる。このようにノイズを含んだスペクトルを用いた場合、定性定量分析の信頼性が低くなる。このため、有効なスペクトル範囲だけを取り出すことが必要となる。
【００２０】
以下、有効スペクトル範囲（単に、有効範囲ともいう）の決定方法について説明する。
まず、図１に示すスペクトルの有効範囲の上限Ｘmaxを決定するため、吸収強度（光学密度）Ｙの絶対値｜Ｙ｜に対する規定値ＹlimをＹlim＝２．５、また、規定の間隔ＸgapをＸgap＝０．２[THz]として、Ｘmaxを求めたところ、Ｘmax＝３．７８[THz]を得た。以下、このＸmaxの決定方法の手順について説明する。
【００２１】
規定値Ｙlim＝２．５と絶対値｜Ｙ｜との交点Ｘｚを基準に、規定の間隔Ｘgap＝０．２[THz]だけ離れた両側の隣接する基準点Ｘａ、Ｘｂ（Ｘａ＜Ｘｂ）における絶対値｜Ｙ｜の値Ｙａ、Ｙｂを検出する。Ｙａ又はＹｂがともに規定値Ｙlim＝２．５を越えない場合は、ＸｚはＸminにもＸmaxにもならず、ＹａまたはＹｂのいずれかがＹlim＝２．５を超えた場合、Ｙlim＝２．５を超えない側のＸの範囲をＸの有効範囲と判断し、ＸｚをＸminまたはＸmaxとする。
【００２２】
図１に示すスペクトルには、２[THz]付近、３．３[THz]付近、３．８[THz]付近、及びそれ以上の高周波領域にＹlim＝２．５を超えるピークが観測される。図１に示すように、２[THz]付近、３．３[THz]付近のピークにおいて、規定値Ｙlim＝２．５と絶対値｜Ｙ｜との交点Ｘｚを基準に、規定の間隔Ｘgap＝０．２[THz]だけ離れた両側の隣接する基準点Ｘａ、Ｘｂ（Ｘａ＜Ｘｂ）における絶対値｜Ｙ｜の値Ｙａ及びＹｂは、ともにＹlim＝２．５より小さいため、ＸｚはＸminまたはＸmaxにならない。
【００２３】
図１に示す３．８[THz]付近のピークでは規定値Ｙlim＝２．５と絶対値｜Ｙ｜との交点ＸｚはＸｚ＝３．７８[THz]であり、交点Ｘｚを基準に、規定の間隔Ｘgap＝０．２[THz]だけ離れた両側の隣接する基準点Ｘａ、Ｘｂにおける絶対値｜Ｙ｜の値Ｙａ及びＹｂのうち、ＹａはＹlim＝２．５を超えないが、ＹｂはＹlim＝２．５を超えることが観測される。従って、Ｘｚ＝３．７８[THz]がＸmaxとなる。尚、規定値Ｙlimとピーク線との交点は、交点がピーク頂点である場合を除いて左右２点あるが、Ｘmaxは左右２点の交点のうち、Ｘの値が小さい方とする。
【００２４】
図１において、３．８[THz]付近以上の高周波領域では、ほぼ一定の狭い周波数間隔でＹの値が激しく上下しており、ノイズとみなせることからしても、この判定手順により決定されたＸmax＝３．７８[THz]という値はスペクトル形状に照らして妥当な値とみなされる。
【００２５】
次に、吸収強度Ｙの一次微分Ｙ'を計算する。図２は、図１に示す吸収強度Ｙの一次微分Ｙ'を示す図である。｜Ｙ'｜に対する規定値Ｙ'lim＝５０、規定の間隔Ｘgap＝０．２[THz]としてＸmaxを求めたところ、３．７８[THz]を得た。この決定方法は、前述のＹに対するＸmaxを求めた手法と同様であり、規定値Ｙ'lim＝５０と絶対値｜Ｙ'｜との交点Ｘｚ＝３．７８を基準に、規定の間隔Ｘgap＝０．２[THz]だけ離れた両側の隣接する基準点Ｘａ、Ｘｂ（Ｘａ＜Ｘｂ）における絶対値｜Ｙ'｜の値Ｙ'ａ、Ｙ'ｂを検出する。Ｙ'ａ又はＹ'ｂがともに規定値Ｙ'lim＝５０を越えない場合は、ＸｚはＸminにもＸmaxにもならず、Ｙ'ａまたはＹ'ｂのいずれかがＹ'lim＝５０を超えた場合、Ｙ'lim＝５０を超えない側のＸの範囲をＸの有効範囲と判断し、ＸｚをＸminまたはＸmaxとする。この場合、Ｙ'bがＹlim＝５０を超えるので、Ｘｚ＝３．７８がＸmaxとなる。
【００２６】
次に、吸収強度Ｙの二次微分Ｙ"を計算する。図３は、図１に示す吸収強度Ｙの二次微分Ｙ"を示す図である。｜Ｙ"｜に対する規定値Ｙ"lim＝４０００、規定の間隔Ｘgap＝０．２[THz]としてＸmaxを求めたところ、３．６０[THz]を得た。この決定方法は、前述のＹ及びＹ'に対するＸmaxを求めた手法と同様であり、規定値Ｙ"lim＝４０００と絶対値｜Ｙ"｜との交点Ｘｚ＝３．６を基準に、規定の間隔Ｘgap＝０．２[THz]だけ離れた両側の隣接する基準点Ｘａ、Ｘｂ（Ｘａ＜Ｘｂ）における絶対値｜Ｙ"｜の値Ｙ"ａ、Ｙ"ｂを検出する。Ｙ"ａ又はＹ"ｂがともに規定値Ｙ"lim＝４０００を越えない場合は、ＸｚはＸminにもＸmaxにもならず、Ｙ"ａまたはＹ"ｂのいずれかがＹ"lim＝４０００を超えた場合、Ｙ"lim＝４０００を超えない側のＸの範囲をＸの有効範囲と判断し、ＸｚをＸminまたはＸmaxとする。この場合、Ｙ"ｂがＹ"lim＝４０００を超えるので、Ｘｚ＝３．６をＸmaxと決定する。
【００２７】
以上のように、吸収強度Ｙ、その一次微分Ｙ'、二次微分Ｙ"のそれぞれを用いた手法により決定されたＸmaxが異なるのは、それぞれの物理的な意味の違いによるものであり、｜Ｙ｜を用いて限界値を決定する場合は、吸収の大きなところでノイズの影響が大きくなることに重みを置いており、|Ｙ'｜を用いた場合には、ノイズに由来するピークの傾斜が大きいことに重みを置いており、|Ｙ"｜を用いた場合には、ノイズに由来するピーク頂点での傾きの変化が大きいことに重みを置いている。
【００２８】
次に、Ｘの有効範囲の下限値Ｘminを求める。｜Ｙ｜に対する規定値Ｙlim＝０．２、規定の間隔Ｘgap＝０．２[THz]として、前述と同様の方法で、Ｘminを求めたところ、０．１４[THz]を得た。同様に、｜Ｙ'｜に対する規定値Ｙ'lim＝２５、規定の間隔Ｘgap＝０．２[THz]としてＸminを求めたところ、０．１１[THz]を得た。同様に、｜Ｙ"｜に対する規定値Ｙ"lim＝２０００、規定の間隔Ｘgap＝０．２[THz]としてＸminを求めたところ、０．１２[THz]を得た。
【００２９】
以上のように、吸収強度Ｙ、その一次微分Ｙ'、二次微分Ｙ"のそれぞれを用いた手法により決定されたＸmax、Ｘminのうち、最小のＸmaxと最大のＸminに囲まれたＸの範囲を、Ｘの有効範囲であると判定し、例えば、図１に示すスペクトルの有効範囲は、０．１４[THz]〜３．６０[THz]と決定される。
【００３０】
ここで、従来の方法、例えば、前述の非特許文献２の手法に基づいてＸmaxを求めた場合と比較する。図４は、図１に示すスペクトルに関して非特許文献２の手法に基づいて求めたＸmaxの位置を示す図である。非特許文献２に記載の方法では、Ｘ＝２．０[THz]付近の大きな吸収ピークをノイズと誤判定してしまうため、上限は１．９７[THz]と計算された。大きなピークは、その物質を特徴付けるためのスペクトルの最重要な情報の一つであるため、ここで有効限界と誤判定されることは、このスペクトルを用いた定性定量分析を行う上で非常に不利となる。また、Ｘの範囲が狭まるように誤判定されると、データポイントの数も格段に小さくなるため、情報量の損失が大きく、信頼性が損なわれることとなる。
【００３１】
このように、本実施の形態によれば、客観的な手法で、鋭いピークを有するスペクトルにおいても、吸収ピークをノイズと誤判定することなく、スペクトル演算処理を行う上で有効なスペクトル範囲を取り出すことが可能である。
なお、本実施の形態における規定値Ｙlim、Ｙ'lim、Ｙ"lim、及び、規定の間隔Ｘgapはここで述べた値に限らず、試料や実験条件などに基づいて任意に設定可能である。
【００３２】
また、上述したＸの有効範囲、即ち、スペクトルの有効範囲を決定する処理は、コンピュータに記憶されたアプリケーションプログラムによって実現される。このプログラムは、記憶媒体に記憶して流通することも可能であり、また、ネットワークを介して流通させることも可能である。
また、上述したＸの有効範囲は、エネルギーで規定することも可能である。すなわち、エネルギーが４．１３５７×１０−４［ｅＶ］〜４．１３５７×１０−２［ｅＶ］の範囲として適用することが可能である。この範囲は、上述の周波数、波長、波数等のエネルギー単位を換算することによって算出することが可能である。
【００３３】
以上、添付図面を参照しながら、本発明に係るスペクトル解析方法の好適な実施形態について説明したが、本発明はかかる例に限定されない。当業者であれば、本願で開示した技術的思想の範疇内において、各種の変更例又は修正例に想到し得ることは明らかであり、それらについても当然に本発明の技術的範囲に属するものと了解される。
【産業上の利用可能性】
【００３４】
このように、本発明によれば、客観的な手法で、鋭いピークを有するスペクトルにおいても、吸収ピークをノイズと誤判定することなく、スペクトル演算処理を行う上で有効なスペクトル範囲を取り出すことが可能であり、スペクトルに基づいた定性定量解析の精度を向上することが可能である。
【００３５】
また、上述のテラヘルツ時間領域分光測定装置の機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することによりスペクトル解析を行ってもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。
【００３６】
また、「コンピュータシステム」は、ＷＷＷシステムを利用している場合であれば、ホームページ提供環境（あるいは表示環境）も含むものとする。
また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含むものとする。また上記プログラムは、前述した機能の一部を実現するためのものであっても良く、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであってもよい。
【図面の簡単な説明】
【００３７】
【図１】測定された周波数領域のテラヘルツスペクトル（吸収強度Ｙ）を示す図
【図２】吸収強度Ｙの一次微分値を示す図
【図３】吸収強度Ｙの二次微分値を示す図
【図４】従来方法に基づいて求められた上限値Ｘmaxの位置を示す図
【符号の説明】
【００３８】
Ｘ………周波数
Ｙ………吸収強度
Ｙ'………吸収強度Ｙの一次微分
Ｙ"………吸収強度Ｙの二次微分
【図１】
0005121415.tif 000002
【図２】
0005121415.tif 000003
【図３】
0005121415.tif 000004
【図４】
0005121415.tif 000005