JPS5895222A - ピーク波数位置検出方法 - Google Patents
ピーク波数位置検出方法Info
- Publication number
- JPS5895222A JPS5895222A JP19343381A JP19343381A JPS5895222A JP S5895222 A JPS5895222 A JP S5895222A JP 19343381 A JP19343381 A JP 19343381A JP 19343381 A JP19343381 A JP 19343381A JP S5895222 A JPS5895222 A JP S5895222A
- Authority
- JP
- Japan
- Prior art keywords
- wave number
- data
- number position
- peak
- peak wave
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000001514 detection method Methods 0.000 title claims description 7
- 238000000034 method Methods 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000004304 visual acuity Effects 0.000 abstract 3
- 238000003491 array Methods 0.000 abstract 1
- CLOMYZFHNHFSIQ-UHFFFAOYSA-N clonixin Chemical compound CC1=C(Cl)C=CC=C1NC1=NC=CC=C1C(O)=O CLOMYZFHNHFSIQ-UHFFFAOYSA-N 0.000 abstract 1
- 238000002834 transmittance Methods 0.000 description 9
- 238000005070 sampling Methods 0.000 description 5
- 238000010521 absorption reaction Methods 0.000 description 3
- 230000004069 differentiation Effects 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 230000001580 bacterial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/02—Details
- G01J3/06—Scanning arrangements arrangements for order-selection
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Measurement Of Radiation (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
本発明は、分光光度計等におけるピーク波数（波長）位
＠検出方法に関する。
＠検出方法に関する。
分光光度計、特に赤外分光光度計では、吸収強度のピー
ク波数位置を正龜に検出する。ことは非常に重要であシ
、このためのピーク波数位置検出方法の先行技術として
、各波数に対する透過度のテーク列を１次微分し、その
１次微分饋か負値から正値に変化したときの反数位置を
ピーク反数ｆｆｌはとするようにしたものが旋条されて
おシ、この提案においては、いわゆる５ａｖｉｔｚｋｙ
ａｎｆｆｌ Ｇｏｌａｙ ｃＡ−５ａｖｉｔｚｋｙ＋
Ｍ、Ｊ−Ｅ、 Ｇｏｌａｙ：Ａｎａｌ、 Ｃｈｅｍ、
、 Ｖｏｌ ８６゜１６２７ （１９６４）の文献参照
〕によって１次倣分１直を求める方法が最良とされてい
る。ところが、５ａｖｉｔｚｋｙ ａｎｄ Ｇｏｌａｙ
による１次微分値の方法にあっては、（１）その方
法の本質的なものであるが、１次微分の分解能に基因す
る波数位置の誤差を生じること、ならびに（２）波数−
透過度のグラフ１こよってあられされる曲線における比
較的平担な部分、例えばベースラインの近くや幅が広い
部分などではノイズの影響を受けやすく、したかって真
のピークでないものをピークであると判定してしまうこ
と、という問題点がある。ここで、前出の文献における
５ａｖｉｔｚｋｙ ａｎｄ Ｇｏｌａｙ による１次
微分法を砥略的に説明し、上記問題点の理解ｆこ供する
。
ク波数位置を正龜に検出する。ことは非常に重要であシ
、このためのピーク波数位置検出方法の先行技術として
、各波数に対する透過度のテーク列を１次微分し、その
１次微分饋か負値から正値に変化したときの反数位置を
ピーク反数ｆｆｌはとするようにしたものが旋条されて
おシ、この提案においては、いわゆる５ａｖｉｔｚｋｙ
ａｎｆｆｌ Ｇｏｌａｙ ｃＡ−５ａｖｉｔｚｋｙ＋
Ｍ、Ｊ−Ｅ、 Ｇｏｌａｙ：Ａｎａｌ、 Ｃｈｅｍ、
、 Ｖｏｌ ８６゜１６２７ （１９６４）の文献参照
〕によって１次倣分１直を求める方法が最良とされてい
る。ところが、５ａｖｉｔｚｋｙ ａｎｄ Ｇｏｌａｙ
による１次微分値の方法にあっては、（１）その方
法の本質的なものであるが、１次微分の分解能に基因す
る波数位置の誤差を生じること、ならびに（２）波数−
透過度のグラフ１こよってあられされる曲線における比
較的平担な部分、例えばベースラインの近くや幅が広い
部分などではノイズの影響を受けやすく、したかって真
のピークでないものをピークであると判定してしまうこ
と、という問題点がある。ここで、前出の文献における
５ａｖｉｔｚｋｙ ａｎｄ Ｇｏｌａｙ による１次
微分法を砥略的に説明し、上記問題点の理解ｆこ供する
。
第１図は、１次微分、直を求めることによって、透過光
の強度（透過度）がグラフ上、下向きのピークとしてあ
られされる技数（ピーク仮数）の装置を検出する場合の
基本的アルゴリズム（Ａｌｇｏｒｉｔｈｍ）を示すため
のものであって、第１１（１）は、横綱を波数（波長）
Ｖに、ｋ軸を透過度Ｔにとるグラフであシ、第１図（２
）は横軸を波数（波長）Ｉ／に、縦軸を１次微分値（ｄ
″ｒ／ｄｌ／）にとるグラフである。
の強度（透過度）がグラフ上、下向きのピークとしてあ
られされる技数（ピーク仮数）の装置を検出する場合の
基本的アルゴリズム（Ａｌｇｏｒｉｔｈｍ）を示すため
のものであって、第１１（１）は、横綱を波数（波長）
Ｖに、ｋ軸を透過度Ｔにとるグラフであシ、第１図（２
）は横軸を波数（波長）Ｉ／に、縦軸を１次微分値（ｄ
″ｒ／ｄｌ／）にとるグラフである。
分光光度計に入射された白色光のうち、反数／’。
における特定の光の透過度は、第１図（１）に示すよう
に波数ビΦの付近において他の光の透過度と比較して最
も小さくなシ、シたがって、波数（ｅ長））ｏは、ピー
ク波数位−をあられしている。第１図（２）は、第１図
（１）であられされるグラフに対応する１次微分値（ｄ
：１−／ｄｌ／ ）をあられすグラフであシ、ピーク波
数位置ヒ。のところにおいて１次微分値（ｄＴ／ｃ＋Ｈ
）が負値から正値に変化するところは、反数Ｖ。がピー
ク波数位置であることに対応している。
に波数ビΦの付近において他の光の透過度と比較して最
も小さくなシ、シたがって、波数（ｅ長））ｏは、ピー
ク波数位−をあられしている。第１図（２）は、第１図
（１）であられされるグラフに対応する１次微分値（ｄ
：１−／ｄｌ／ ）をあられすグラフであシ、ピーク波
数位置ヒ。のところにおいて１次微分値（ｄＴ／ｃ＋Ｈ
）が負値から正値に変化するところは、反数Ｖ。がピー
ク波数位置であることに対応している。
第２図は、第１図における１次微分値ｆこよる方法によ
ってピーク波数位置を検出する場合のフローチャートを
あられしておシ、ステップ５Ｔ１にお０て、（ｄＴ／ｄ
Ｊ／）がビービ、１１こおいて正であるか否かが判断さ
れ、正でない（ＮＯ）と判断されると、ステップＳＴ２
においてピーク波数位置でないと処理され、正である（
ＹＥＳ ）と判断されると、次のステップＳＴ３に進む
。ステップ５Ｔａ ｉこおいて、（ｄ Ｔ／ｄ Ｖ ）
がヒ＝ ｔｄ−ｉ ＃＋こおいて負でちるか否かか判断
され、負でない（ＮＯ）と判断されると、ステップＳＴ
２に進み、負である（ＹＥＳ）と判断されると、ステッ
プＳＴ４においてピーク波数位置であると処理される。
ってピーク波数位置を検出する場合のフローチャートを
あられしておシ、ステップ５Ｔ１にお０て、（ｄＴ／ｄ
Ｊ／）がビービ、１１こおいて正であるか否かが判断さ
れ、正でない（ＮＯ）と判断されると、ステップＳＴ２
においてピーク波数位置でないと処理され、正である（
ＹＥＳ ）と判断されると、次のステップＳＴ３に進む
。ステップ５Ｔａ ｉこおいて、（ｄ Ｔ／ｄ Ｖ ）
がヒ＝ ｔｄ−ｉ ＃＋こおいて負でちるか否かか判断
され、負でない（ＮＯ）と判断されると、ステップＳＴ
２に進み、負である（ＹＥＳ）と判断されると、ステッ
プＳＴ４においてピーク波数位置であると処理される。
次に、５ａｖｉｔｚｋｙ ａｎｄ Ｇｏｌａｙ ｌこ
おいて、１次微分値を求める場合につき下記の表１にし
たがって説明する。
おいて、１次微分値を求める場合につき下記の表１にし
たがって説明する。
表１の波数は、第１図のグラフの横軸上における波数に
対応しておシ、１次倣分１＠は表１の波数毎に対応して
求める。第１図のグラフから明らかなように、透過度Ｔ
が極小のピークとなる波数はど。であるが、５ａｖｉｔ
ｚｋｙ ａｎｄ Ｇｏｌａｙ においてこの波数Ｖ。
対応しておシ、１次倣分１＠は表１の波数毎に対応して
求める。第１図のグラフから明らかなように、透過度Ｔ
が極小のピークとなる波数はど。であるが、５ａｖｉｔ
ｚｋｙ ａｎｄ Ｇｏｌａｙ においてこの波数Ｖ。
での１次微分値を求めるためには、波数Ｐ。
を中心として、前後の波数ど−ｎ〜Ｐ ＋ｎまで合計（
２ｎ＋１）個（波数ビ０も含める）のサンプリ表 １ ンクポイントでの各透過度Ｔ −ｎ ＋ Ｔ−（ｎ−１
）ｌ・・・Ｈ］’（１，＋・・・、 Ｔ＋ｎを多項式で
ｆｉｔｔｉｎｇ Ｌ、それを１次微分した値と等価な変
数をａｊ（ｊ ＝−ｎ”−＋ｎ）と疋め、次式Ｔ＊＝（
１／Ｎ）カａｊ−Ｔｊをもって、波数２０におけＪゝｎ る１次倣分値（ｄＴ／ｄ ）’ ） とする。ここ
で、Ｎはｒ巳絢 規格化するための定数である。
２ｎ＋１）個（波数ビ０も含める）のサンプリ表 １ ンクポイントでの各透過度Ｔ −ｎ ＋ Ｔ−（ｎ−１
）ｌ・・・Ｈ］’（１，＋・・・、 Ｔ＋ｎを多項式で
ｆｉｔｔｉｎｇ Ｌ、それを１次微分した値と等価な変
数をａｊ（ｊ ＝−ｎ”−＋ｎ）と疋め、次式Ｔ＊＝（
１／Ｎ）カａｊ−Ｔｊをもって、波数２０におけＪゝｎ る１次倣分値（ｄＴ／ｄ ）’ ） とする。ここ
で、Ｎはｒ巳絢 規格化するための定数である。
このようにして、求めた係数ａｊを縦軸に、波数を横軸
にとるグラフを、第８図（１）に示す。横軸上は、波数
ビーｎ ＋ ）’Ｌ（ｎ ｌ）＋ ”’ ｒ Ｊ’ ２
＋と一１＋ｆ’ｏ＋）’＋１ ＋ ｌ’＋２ ｐ・・
・、と十□の合計（２ｎ＋１ ） イ固のサンプリン
グポイントの個々の位置をあられし、縦軸は、それらの
個々の位１における上記係数ａ−ｎ。
にとるグラフを、第８図（１）に示す。横軸上は、波数
ビーｎ ＋ ）’Ｌ（ｎ ｌ）＋ ”’ ｒ Ｊ’ ２
＋と一１＋ｆ’ｏ＋）’＋１ ＋ ｌ’＋２ ｐ・・
・、と十□の合計（２ｎ＋１ ） イ固のサンプリン
グポイントの個々の位置をあられし、縦軸は、それらの
個々の位１における上記係数ａ−ｎ。
ａ−（ｎｌ）＋”’＋ ａ２ ’＋ａ １ ＋
ａＱ、ａ＋ｌ ｌ ａ＋２＋”°＋ａ＋ｏの各
々の値をあられしている。ところで、第１図（１）のよ
うなグラフであられされる曲線でピーク波数位置を検出
する場合（こはこの方法は正しくその位置を求めること
ができるが、＠３図（２）のグラフに示すように、２個
のピークを有する曲線（実線部分）についてピーク波数
位置を検出する場合で、これらの２個のピークの間隔が
、１次微分関数の分解能Ｒ（＝−Ｍ・Δと１ここて、↑
、１ハｆンプリングポイント数（２ｎ ＋１ ）、Δは
淋サンプリングポイント間隔、即ち、波数間隔である〕
よシ小さい場合には、この５ａｖｉｔｚｋｙ ａｎｄ
Ｇｏｌａｙ による１次微分値を求める方法では、破
線に示すところに１個のピークが恰も存在するものとし
て、１次微分し、第３図（３）に示すような微分曲線に
対応した１次微分値を求めてしまい、正しくピーク波数
位置を検出することができない。
ａＱ、ａ＋ｌ ｌ ａ＋２＋”°＋ａ＋ｏの各
々の値をあられしている。ところで、第１図（１）のよ
うなグラフであられされる曲線でピーク波数位置を検出
する場合（こはこの方法は正しくその位置を求めること
ができるが、＠３図（２）のグラフに示すように、２個
のピークを有する曲線（実線部分）についてピーク波数
位置を検出する場合で、これらの２個のピークの間隔が
、１次微分関数の分解能Ｒ（＝−Ｍ・Δと１ここて、↑
、１ハｆンプリングポイント数（２ｎ ＋１ ）、Δは
淋サンプリングポイント間隔、即ち、波数間隔である〕
よシ小さい場合には、この５ａｖｉｔｚｋｙ ａｎｄ
Ｇｏｌａｙ による１次微分値を求める方法では、破
線に示すところに１個のピークが恰も存在するものとし
て、１次微分し、第３図（３）に示すような微分曲線に
対応した１次微分値を求めてしまい、正しくピーク波数
位置を検出することができない。
したがって、本発明の目的は、上述の問題点を解消する
ためになされたものであって、極めて効率的にかつ正確
にピーク波数位置を検出することのできる、分光光度計
等におけるピーク波数位置検出方法を提供することであ
る。
ためになされたものであって、極めて効率的にかつ正確
にピーク波数位置を検出することのできる、分光光度計
等におけるピーク波数位置検出方法を提供することであ
る。
本発明の構成を、その作用とともに記載すると、ピーク
波数位１ｄを検出するために横軸上の各波数位置に個別
的に対応する縦軸上の各データ列から、ある波数位置に
対するデータの１次微分１直を算出し、その１次微分値
が、負値から正値へ、または正値から負値へ変化したと
きにその１次微分埴に対応する波数位−を、データがピ
ークとなるピーク波数位置であると判定するピーク波数
位置検出方法において、大小関係を比較される２つのデ
ータは、前記１次倣分値の算出胸係の分解り目の１／２
の波数間隔だけ陥れた波数位置関係になるものを選定し
、ピーク波数位置からＮ ＝ （１） ＋ １ （ただ
し、Ｎはある波数位−での１次微分値を求めるときの前
記データの列の数、〔〕はＮの整数部である）だけ丙出
れた他の波数位置へ向って順法、前記波数位置関係にあ
る２つのデータの大小関係を比較して、ピーク波数位置
を検出する方法であって、この方法１こよって２つのデ
ータの大小−係を比較していって前記のピーク波数位置
が真のピーク波数位置であるか否かを正確に検出するこ
とができる。
波数位１ｄを検出するために横軸上の各波数位置に個別
的に対応する縦軸上の各データ列から、ある波数位置に
対するデータの１次微分１直を算出し、その１次微分値
が、負値から正値へ、または正値から負値へ変化したと
きにその１次微分埴に対応する波数位−を、データがピ
ークとなるピーク波数位置であると判定するピーク波数
位置検出方法において、大小関係を比較される２つのデ
ータは、前記１次倣分値の算出胸係の分解り目の１／２
の波数間隔だけ陥れた波数位置関係になるものを選定し
、ピーク波数位置からＮ ＝ （１） ＋ １ （ただ
し、Ｎはある波数位−での１次微分値を求めるときの前
記データの列の数、〔〕はＮの整数部である）だけ丙出
れた他の波数位置へ向って順法、前記波数位置関係にあ
る２つのデータの大小関係を比較して、ピーク波数位置
を検出する方法であって、この方法１こよって２つのデ
ータの大小−係を比較していって前記のピーク波数位置
が真のピーク波数位置であるか否かを正確に検出するこ
とができる。
次に、本発明の一実施例について詳細に説明する。下記
の−Ｈｘは、本件実施例の説明に供するものである。
の−Ｈｘは、本件実施例の説明に供するものである。
先ず、１次碩分１匣古求めるためには、一般に、ある仮
数ＶｏＯ位置をサンプリングポイントの中心蚤 として、その前後にｎ個、合計（２ｎ＋１） 個（以
下、Ｍという）のサンプリンクポイントを定め、このサ
ンプリングポイントの数Ｍに対応するデータが必要とさ
れる。この場合、１次微分値を求めるための微分関数の
分解能Ｒは−ｐＭ・Δｒ（ただし、Δとはデータ間の波
数間隔である）であられされるので、波数Ｊ’ｏの位置
の前後で１次微分値の符号が、例えば負から正になシ、
シかも波数りを中心として前後、１／２ Ｒ（前記分解
能の１／２ ） ＝１／４Ｍ・Δどの波数範囲内で波数
ｒθの位１に対応するデータよシも小さいデータがなけ
れば、波数７’ｏの位置に対応するデータをピークであ
ると判定してよい。このため、実際には、波数位置が−
に、−（ｋ−１）ｋ、に−１にそれぞれ対応するデータ
、例えば透過層Ｔ−に、Ｔ−（ｋ−リ＋ Ｔｋ＋ Ｔｋ
−１が、Ｔｋ≧Ｔ−＋Ｃ（ｋ−］）・Ｔｋ ≧Ｔ（ｋ
−１）、（ｋ−１〜ｎ）の菌体式を全て一足するという
条件を付加すればよいのであるが、この条件では、デー
タにノイズが含まれてぃなびものでなければ、はとんど
リジェクトされてしまうことになる工 本件実施例では、このようなリジェクトされるおそれの
ある上記条件を使用せず表■に示すような条件の判別用
フィルタを付加するものであシ、これによって後述する
ように、ノイズの影響が大きく緩和された、したがって
最も効率的にかつ正確（こピーク波数位置を検出するこ
とができるものである。本件実施例の判別用フィルタに
ついて説明する。本件実施例にあっては、ある波数位置
に対応するデータと、その隣シの波数位置に対応するデ
ータとの大小関係を比較するようにはしないで、１次微
分値の算出関数の分解能をＲとした場合に、その１／２
−Ｒだけ、即ち、−（２ｎ＋１）Δｒだけ離れた波数位
置間１．係にある２つのデータの大小関係を比較するよ
うにし、しかも、表■の比較式の全てが成立し、かつ、
波数と。の前後で１次微分値が負値から正値（その逆も
可）になったときに、その波数Ｋ。の位置に対応するデ
ータ（透過度）をピークであると判定して、ピーク波数
位置を検出するように、なしたものである。なお、本件
実施例にあっては、ピーク波数位置検出のアルゴリズム
による誤差は、分光器そのものによる波数誤差がないと
仮定した場合には、±Δととなる。また、検出能力はピ
ークとピークとの間隔が、Ｍ・Δ／’−（２ｎ＋１）・
Δど以上酷れていることが必要であシ、また、ピークの
先端における平坦部の幅か（（１）＋ＩＪ・２１以上に
なるものについてもピーク波数位置を検出することがで
きない。
数ＶｏＯ位置をサンプリングポイントの中心蚤 として、その前後にｎ個、合計（２ｎ＋１） 個（以
下、Ｍという）のサンプリンクポイントを定め、このサ
ンプリングポイントの数Ｍに対応するデータが必要とさ
れる。この場合、１次微分値を求めるための微分関数の
分解能Ｒは−ｐＭ・Δｒ（ただし、Δとはデータ間の波
数間隔である）であられされるので、波数Ｊ’ｏの位置
の前後で１次微分値の符号が、例えば負から正になシ、
シかも波数りを中心として前後、１／２ Ｒ（前記分解
能の１／２ ） ＝１／４Ｍ・Δどの波数範囲内で波数
ｒθの位１に対応するデータよシも小さいデータがなけ
れば、波数７’ｏの位置に対応するデータをピークであ
ると判定してよい。このため、実際には、波数位置が−
に、−（ｋ−１）ｋ、に−１にそれぞれ対応するデータ
、例えば透過層Ｔ−に、Ｔ−（ｋ−リ＋ Ｔｋ＋ Ｔｋ
−１が、Ｔｋ≧Ｔ−＋Ｃ（ｋ−］）・Ｔｋ ≧Ｔ（ｋ
−１）、（ｋ−１〜ｎ）の菌体式を全て一足するという
条件を付加すればよいのであるが、この条件では、デー
タにノイズが含まれてぃなびものでなければ、はとんど
リジェクトされてしまうことになる工 本件実施例では、このようなリジェクトされるおそれの
ある上記条件を使用せず表■に示すような条件の判別用
フィルタを付加するものであシ、これによって後述する
ように、ノイズの影響が大きく緩和された、したがって
最も効率的にかつ正確（こピーク波数位置を検出するこ
とができるものである。本件実施例の判別用フィルタに
ついて説明する。本件実施例にあっては、ある波数位置
に対応するデータと、その隣シの波数位置に対応するデ
ータとの大小関係を比較するようにはしないで、１次微
分値の算出関数の分解能をＲとした場合に、その１／２
−Ｒだけ、即ち、−（２ｎ＋１）Δｒだけ離れた波数位
置間１．係にある２つのデータの大小関係を比較するよ
うにし、しかも、表■の比較式の全てが成立し、かつ、
波数と。の前後で１次微分値が負値から正値（その逆も
可）になったときに、その波数Ｋ。の位置に対応するデ
ータ（透過度）をピークであると判定して、ピーク波数
位置を検出するように、なしたものである。なお、本件
実施例にあっては、ピーク波数位置検出のアルゴリズム
による誤差は、分光器そのものによる波数誤差がないと
仮定した場合には、±Δととなる。また、検出能力はピ
ークとピークとの間隔が、Ｍ・Δ／’−（２ｎ＋１）・
Δど以上酷れていることが必要であシ、また、ピークの
先端における平坦部の幅か（（１）＋ＩＪ・２１以上に
なるものについてもピーク波数位置を検出することがで
きない。
次に、本件実施例においては、あるピーク波数から他の
波数位置へ向って２つのデータの大小関係を比較してい
く場合、大小関係を比較されるべき最初の２つのデータ
は、先ずピーク波数Ｐ。の位置に対応するデータ１と、
その位置から＝（１）−１即ち、表■からいって−Ｌ−
１だけ離れた波数Ｐ−Ｌ−１の位置に対応するデータＴ
−Ｌ−１であシ、これらのデータｌとＴ−Ｌ−１との大
小関係を比較し、データＴ。の方が小さければ、データ
Ｔ。に対応する波数と。をピーク波数位置とすることは
このデータの比較に関しては正しいと判定される。次に
、ピーク波数へよシ１つ前の波数Ｆ−１の動１こ対応す
るデータＴ−１と、前記分解能Ｒの１／２だけ、このデ
ータＴ、よシ前にある波数と−Ｌ−２に対応ず−るデー
タＴ−Ｌ−２との大小関係を比較し、テークＴｌＯ方が
データＴＬ２よシも小さいときは、テ゛−タＴ。
波数位置へ向って２つのデータの大小関係を比較してい
く場合、大小関係を比較されるべき最初の２つのデータ
は、先ずピーク波数Ｐ。の位置に対応するデータ１と、
その位置から＝（１）−１即ち、表■からいって−Ｌ−
１だけ離れた波数Ｐ−Ｌ−１の位置に対応するデータＴ
−Ｌ−１であシ、これらのデータｌとＴ−Ｌ−１との大
小関係を比較し、データＴ。の方が小さければ、データ
Ｔ。に対応する波数と。をピーク波数位置とすることは
このデータの比較に関しては正しいと判定される。次に
、ピーク波数へよシ１つ前の波数Ｆ−１の動１こ対応す
るデータＴ−１と、前記分解能Ｒの１／２だけ、このデ
ータＴ、よシ前にある波数と−Ｌ−２に対応ず−るデー
タＴ−Ｌ−２との大小関係を比較し、テークＴｌＯ方が
データＴＬ２よシも小さいときは、テ゛−タＴ。
に対応する波数／１０をピーク波廠位置とするるとはこ
のデータの比較に関しては正しいと判定される。
のデータの比較に関しては正しいと判定される。
このようにして、順次、データＴ−２とデータＴ＋ｓと
の大小関係、・・・、デニタＴ−にとデータＴ−Ｌ−に
−１との大小関係が比較されてゆき、分解能の１／２の
範囲内の全てにおいて、表■の比較式が全て成立したと
き、データＴ。に対応する波数／１０は真のピーク波数
位置であると判定される。
の大小関係、・・・、デニタＴ−にとデータＴ−Ｌ−に
−１との大小関係が比較されてゆき、分解能の１／２の
範囲内の全てにおいて、表■の比較式が全て成立したと
き、データＴ。に対応する波数／１０は真のピーク波数
位置であると判定される。
なお、この実施例では、波数Ｐａよシ前のデータのみに
ついて行ったが、波数／’ｏよシ後のデータについて、
あるいはその両方のデータについ゛て大小関係の比較を
行うこともできる。また、上述の実施例において、Ｔ−
ｋ＜ Ｔ−Ｌ−に−１（またはＴｋ〉ＴＬ十に＋１
でも同じことである）であって、かつ、Ｔ＝ｋ＜ Ｔ−
Ｌ−ｋ（またはＴｋ＜ＴＬ十に＋１）であって−かつ、
ｋ＝０〜ｎ −Ｌとしてもよい１また、本件実施例では
ピークを下向きのものとしたが、土間きのものについて
も同様に芙施することができ、またデータを透一度とし
ているが、ピークがあられれるデータであれば同様に実
施することができる。
ついて行ったが、波数／’ｏよシ後のデータについて、
あるいはその両方のデータについ゛て大小関係の比較を
行うこともできる。また、上述の実施例において、Ｔ−
ｋ＜ Ｔ−Ｌ−に−１（またはＴｋ〉ＴＬ十に＋１
でも同じことである）であって、かつ、Ｔ＝ｋ＜ Ｔ−
Ｌ−ｋ（またはＴｋ＜ＴＬ十に＋１）であって−かつ、
ｋ＝０〜ｎ −Ｌとしてもよい１また、本件実施例では
ピークを下向きのものとしたが、土間きのものについて
も同様に芙施することができ、またデータを透一度とし
ているが、ピークがあられれるデータであれば同様に実
施することができる。
以上説明したように、本発明によれば、大小関係を比較
される２つのデータは、１次微分値の算出関数の分解能
の１／２に相当する波数間隔だけ離れた波数位置関係に
なるものを選定し、ピーク波数位置からＮ−（Ｔ、］
＋１だけ離れた他の波数位、 置へ向って順次、前記波
数位置関係を有する２つのデータの大小関係を比較゛し
てピーク波数位置を検出するように構成されであるので
、極めて効率的にかつ正確にピーク波数位置を検出する
ことができる。即ち、本発明によれば、上述の先行技術
の微分に基因する分解能の低下による誤った波数をピー
クとすることがなく、また幅広い吸収帯をピークとした
シ、ボトムアウト（ｂ、ｏｔｔｏｍ ｏｕｔ）シた吸収
帯のノイズによって、これを数個のピークと判定するこ
ともなく、丈にベースラインのような平坦な部分でノイ
ズによる多数の偽のピークを発生させることもない、等
の効果を発揮することができる。
される２つのデータは、１次微分値の算出関数の分解能
の１／２に相当する波数間隔だけ離れた波数位置関係に
なるものを選定し、ピーク波数位置からＮ−（Ｔ、］
＋１だけ離れた他の波数位、 置へ向って順次、前記波
数位置関係を有する２つのデータの大小関係を比較゛し
てピーク波数位置を検出するように構成されであるので
、極めて効率的にかつ正確にピーク波数位置を検出する
ことができる。即ち、本発明によれば、上述の先行技術
の微分に基因する分解能の低下による誤った波数をピー
クとすることがなく、また幅広い吸収帯をピークとした
シ、ボトムアウト（ｂ、ｏｔｔｏｍ ｏｕｔ）シた吸収
帯のノイズによって、これを数個のピークと判定するこ
ともなく、丈にベースラインのような平坦な部分でノイ
ズによる多数の偽のピークを発生させることもない、等
の効果を発揮することができる。
第１図は、１次微分値によるピーク検出の基本的アルゴ
リズムを示す図、第２図は、第１図によるアルゴリズム
に対応したフローチャート、第３図は先行技術の動作を
説明するためのグラフである。 代理人 弁理士 岡田和秀 第１図 第３図
リズムを示す図、第２図は、第１図によるアルゴリズム
に対応したフローチャート、第３図は先行技術の動作を
説明するためのグラフである。 代理人 弁理士 岡田和秀 第１図 第３図
Claims (3)
- （１）ピーク波数位置を検出するために横軸上の各波数
位置に個別的をこ対応する縦軸上の各データ列から、あ
る波数位置に対するデータの１次微分値を算出し、その
１次微分値が、負部から正値へ、または正−から負値へ
変化したときに、その１次微分値に対応する波数位置を
、データがピークとなるピーク波数位置であると判定す
るピーク波数位置検出方法において、大小関係を比較さ
れる２つのデータは、前記１次微分値の算出関数の分解
能の１／２の波数間隔だけ離れた波数位置関係になるも
のを選定し、ピーク波数位−からＮ −（、−）＋１（
ただし、Ｎはある波数位置での１次微分値を求めるとき
の前記データの列の数であシ、〔〕は小関係を比較して
データがピークとなる波数位置を便用するようにしたこ
とを特徴とする、分光光度計等におけるピーク波数位置
検出方法。 - （２）前記特許請求の範囲第１項に記載の分光光度計等
におけるピーク波数位置検出方法において、ピーク波数
位置から、一方向ｌこある他の波数位置へ向って、前記
波数位置関係を有する２つのデータの大小関係を比較し
ていくことを特徴とする、分光光度計等ｆこおけるビー
ク波数位＠検出方法。 - （３）前記特許請求の範囲第１項に記載の分光光度計等
におけるピーク波数位置検出方法において、ピーク波数
位置から、両方向にある他の波数位置へ向って、前記波
数位置関係を有する２つのデータの大小間係を比較して
いくことを特徴とする、分光光度計等におけるピーク波
数位置検出方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP19343381A JPS5895222A (ja) | 1981-11-30 | 1981-11-30 | ピーク波数位置検出方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP19343381A JPS5895222A (ja) | 1981-11-30 | 1981-11-30 | ピーク波数位置検出方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS5895222A true JPS5895222A (ja) | 1983-06-06 |
| JPH0332727B2 JPH0332727B2 (ja) | 1991-05-14 |
Family
ID=16307894
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP19343381A Granted JPS5895222A (ja) | 1981-11-30 | 1981-11-30 | ピーク波数位置検出方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS5895222A (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH09184761A (ja) * | 1995-12-29 | 1997-07-15 | Horiba Ltd | ２つのスペクトルの一致度算出方法 |
-
1981
- 1981-11-30 JP JP19343381A patent/JPS5895222A/ja active Granted
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH09184761A (ja) * | 1995-12-29 | 1997-07-15 | Horiba Ltd | ２つのスペクトルの一致度算出方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0332727B2 (ja) | 1991-05-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5552997A (en) | Method of calibrating an analytical instrument | |
| WO2017094170A1 (ja) | ピーク検出方法及びデータ処理装置 | |
| CN110108658B (zh) | 一种污染气体红外光谱识别方法及系统 | |
| CN115859080A (zh) | 基于非对称高斯模型的光纤布拉格光栅信号峰值检测算法 | |
| JPS5895222A (ja) | ピーク波数位置検出方法 | |
| JPS6285828A (ja) | 多重チヤンネルパイロメ−タ−を用いて非接触的に温度を測定する方法 | |
| JP3461284B2 (ja) | 赤外線ガス分析計における検量線の作成方法 | |
| JP3174186B2 (ja) | ガス識別装置 | |
| JPH06249810A (ja) | ガス識別装置 | |
| JPH0833345B2 (ja) | デ−タ処理方法 | |
| Zipser | Selectivity of sensor systems | |
| JPS62263404A (ja) | パタ−ン検査装置 | |
| US5945639A (en) | Cancellation of common-mode signals in digitizing tablet | |
| JPS6217716Y2 (ja) | ||
| JP3327625B2 (ja) | 濃度測定装置の校正方法 | |
| JPH02296162A (ja) | ケーブルの部分放電検出方法 | |
| JPS6011172A (ja) | 測定出力のピ−ク検出回路 | |
| JPS63262511A (ja) | パタ−ン検査装置 | |
| JPH06242038A (ja) | ガス識別装置 | |
| JP3235402B2 (ja) | ディジタル保護継電器 | |
| JP3295124B2 (ja) | デバイス評価試験装置及びその方法 | |
| JP2001311689A (ja) | 物体表面の光沢度判別方法 | |
| JPH07210681A (ja) | 回転位置検出装置 | |
| CN118570436A (zh) | 识别方法、装置、设备、存储介质及程序产品 | |
| JPH0339608A (ja) | プリント基板の欠陥検出方法 |