JPH06201673A - Method for detecting peak of chromatogram - Google Patents

Method for detecting peak of chromatogram

Info

Publication number
JPH06201673A
JPH06201673A JP36001792A JP36001792A JPH06201673A JP H06201673 A JPH06201673 A JP H06201673A JP 36001792 A JP36001792 A JP 36001792A JP 36001792 A JP36001792 A JP 36001792A JP H06201673 A JPH06201673 A JP H06201673A
Authority
JP
Japan
Prior art keywords
section
point
peak
value
data
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.)
Pending
Application number
JP36001792A
Other languages
Japanese (ja)
Inventor
Tetsushi Okuyama
哲史 奥山
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shimadzu Corp
Original Assignee
Shimadzu Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shimadzu Corp filed Critical Shimadzu Corp
Priority to JP36001792A priority Critical patent/JPH06201673A/en
Publication of JPH06201673A publication Critical patent/JPH06201673A/en
Pending legal-status Critical Current

Links

Abstract

PURPOSE:To detect the starting point, vertex, and ending point of a peak by one time of scanning by fixing a noticeable point on chromatogram data and calculating whether the measured values in a fixed section after the noticeable point are larger or smaller than the data at the noticeable point, and then, recognizing a point of inflection by moving the noticeable point to the maximum value in the section when the measured values are larger or to the minimum value in the section when the measured values are smaller. CONSTITUTION:A section width W and continuous section value (n) are decided as peak detecting references. The section width W is set at 1/4 of the width at the half maximum to be detected. After resetting the minimum section value, maximum section value, and mean section value of a memory area, the measured values in the section corresponding to the section width W from a noticeable point are found and, when the mean value of the section is larger than the data at the noticeable point, the section is discriminated as a rising section and, when smaller, as a falling section. When, for example, the noticeable point moves from P0 to P2 through P1 and the means value of the section W starting from P2 is smaller than the data at P2, and then, the position of the minimum value in the section is PS, the noticeable point is moved to PS and, when the peak again starts to rise, the point of inflection PS is discriminated as the starting point of the peak and recorded.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】本発明は液体クロマトグラフやガ
スクロマトグラフにおけるデータ処理装置で、クロマト
グラムのピークの開始点、頂点及び終了点をみつけるピ
ーク検出方法に関するものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a peak detection method for finding a starting point, an apex and an ending point of a peak of a chromatogram in a data processing device in a liquid chromatograph or a gas chromatograph.

【0002】[0002]

【従来の技術】クロマトグラムのデータからピークの位
置を検出する方法としては、サンプリングデータを一定
区間ごとにサンプリングし、順次比較をしていってデー
タが増加から減少に転じた点でピークの頂点を検出す
る。ピーク頂点がみつかると、ピーク頂点を中心にもう
一度データを走査することによって、ピーク頂点からサ
ンプリングデータの変化量を順次計算していって所定の
条件が満たされた点をピークの開始点又は終了点として
検出している。
2. Description of the Related Art As a method of detecting the position of a peak from chromatogram data, the sampling data is sampled at regular intervals and compared sequentially, and the peak of the peak changes at the point where the data changes from increasing to decreasing. To detect. When the peak apex is found, the data is scanned again around the peak apex to calculate the amount of change in the sampling data from the peak apex one by one, and the point where the specified condition is satisfied is the start or end point of the peak. Is detected as.

【0003】[0003]

【発明が解決しようとする課題】従来の方法によりピー
クを検出しようとすると、クロマトグラムデータを2回
走査する必要があり、処理時間が長くなる問題がある。
また、ピークの開始点と終了点はピーク検出基準により
求めるが、ピークの鋭さによってその検出基準を変化さ
せる必要がある。鋭いピークと鈍いピークが混ざり合っ
ている場合には、ピーク検出基準を正確に定めなければ
いずれかのピークを検出できないこともある。
When a peak is to be detected by the conventional method, it is necessary to scan the chromatogram data twice, which causes a problem that the processing time becomes long.
Further, the start point and the end point of the peak are obtained by the peak detection reference, but the detection reference needs to be changed depending on the sharpness of the peak. When the sharp peaks and the blunt peaks are mixed, it may not be possible to detect any of the peaks unless the peak detection standard is accurately defined.

【0004】本発明はクロマトグラムデータを1回走査
するだけでピークの開始点、頂点及び終了点を検出する
ことができるようにして処理時間を短かくするととも
に、鋭いピークから鈍いピークの広い範囲にわたってピ
ーク検出を可能にする方法を提供することを目的とする
ものである。
The present invention makes it possible to detect the start point, apex and end point of a peak by scanning the chromatogram data only once, thereby shortening the processing time and in a wide range from a sharp peak to a dull peak. It is an object of the present invention to provide a method that enables peak detection over a period of time.

【0005】[0005]

【課題を解決するための手段】本発明ではクロマトグラ
ムデータ上に注目点を定め、その点から後の一定区間内
のデータが注目点に対し上昇しているか下降しているか
を計算し、上昇であればその区間内の最大値に注目点を
移動させ、下降であれば区間内の最小値に移動させるこ
とによりピークの変曲点を捉えてピークを抽出する。そ
のため、本発明ではクロマトグラムデータについてサン
プリングデータを複数個含む一定区間幅Wと、データの
上昇又は下降を判定するための連続区間値nを定め、ク
ロマトグラム上に注目点を定め、その注目点から一定区
間幅内のデータから区間最小値、区間最大値及び区間平
均値を求め、区間平均値が注目点の値より大きい場合は
その区間を上昇、小さい場合は下降と判定し、上昇区間
では注目点をその区間の最大値の点へ移動させ、下降区
間では注目点をその区間の最小値の点へ移動させること
によりピークの変曲点を捉え、n区間以上連続して下降
した後の変曲点をピーク開始点もしくはピーク終了点と
し、n区間以上連続して上昇した後の変曲点をピーク頂
点とする。
In the present invention, a point of interest is set on the chromatogram data, and whether the data within a certain section after that point is rising or falling with respect to the point of interest is calculated, If it is, the point of interest is moved to the maximum value in the section, and if it is descending, it is moved to the minimum value in the section to catch the inflection point of the peak and extract the peak. Therefore, in the present invention, a fixed section width W including a plurality of sampling data for chromatogram data and a continuous section value n for judging rise or fall of the data are set, and a point of interest is set on the chromatogram. From the data within a certain section width, the section minimum value, the section maximum value, and the section average value are obtained.If the section average value is larger than the value of the point of interest, the section is raised, and if it is smaller, it is determined to be descending. The point of inflection is detected by moving the point of interest to the point of maximum value in the section and moving the point of interest to the point of minimum value in the section in the descending section, and after descending continuously for n sections or more. The inflection point is defined as the peak start point or the peak end point, and the inflection point after continuously rising for n sections or more is defined as the peak apex.

【0006】[0006]

【実施例】図1から図3により実施例を説明する。図1
と図2は一実施例の手順をPAD(Problem Analysis D
iagram)により記述したものである。図1に示されるよ
うに、ピーク検出基準として区間幅Wと、その区間が上
昇区間であるのか、又は下降区間であるのかを判定する
ための連続区間値nを定める。これらはノイズの周期と
検出したい最小ピークを基にして、例えば次のように求
めることができる。 W=(検出したいピークの半値幅の1/4以下の値) n=(ノイズの周期)/(サンプリングレート)程度の
値 Wとnをともに小さめに設定することにより、ノイズも
微小ピークとして検出し、その後、ピーク面積の小さい
ピークを左右の連続したピークに含めてしまう方法を採
ることにより、精度よくピーク検出することができる。
一例としてWをデータ数20に該当する区間、nを4と
設定する。
EXAMPLES Examples will be described with reference to FIGS. 1 to 3. Figure 1
2 and FIG. 2 show the procedure of one embodiment in PAD
iagram). As shown in FIG. 1, a section width W and a continuous section value n for determining whether the section is an ascending section or a descending section are determined as peak detection criteria. These can be obtained, for example, as follows based on the noise cycle and the minimum peak to be detected. W = (value less than or equal to ¼ of half-width of peak to be detected) n = (noise period) / (sampling rate) value By setting both W and n small, noise is also detected as a minute peak Then, after that, by adopting a method of including peaks having a small peak area in continuous peaks on the left and right, it is possible to detect peaks with high accuracy.
As an example, W is set to a section corresponding to the number of data 20 and n is set to 4.

【0007】(現在位置)、(ピーク最小位置)、(ピ
ーク最大位置)、(区間最小値)、(区間最大値)及び
(区間平均値)のためのメモリ領域を用意する。初期状
態としてはデータはn回連続して下降状態にあるものと
し、最初のデータを(現在位置)に記録し、メモリ領域
の(ピーク最小位置)と(ピーク最大位置)の記録を消
す。
Memory areas for (current position), (peak minimum position), (peak maximum position), (section minimum value), (section maximum value) and (section average value) are prepared. As an initial state, the data is assumed to be in a descending state for n consecutive times, the first data is recorded at (current position), and the recording of (peak minimum position) and (peak maximum position) in the memory area is erased.

【0008】メモリ領域の(区間最小値)、(区間最大
値)及び(区間平均値)がリセットされた後、(現在位
置)から(現在位置+区間幅W)までの区間における、
区間最小値、区間最大値及び区間平均値を求められ、そ
れぞれのメモリ領域に記録される。区間平均が現在位置
のデータより大きい場合はその区間は上昇区間と判断さ
れ、小さい場合はその区間は下降区間と判断される。上
昇の場合は図2(A)に従って処理され、下降の場合は
図2(B)に従って処理される。
In the section from (current position) to (current position + section width W) after (section minimum value), (section maximum value) and (section average value) of the memory area are reset,
A section minimum value, a section maximum value, and a section average value are obtained and recorded in each memory area. If the section average is larger than the data of the current position, the section is judged to be an ascending section, and if it is smaller, the section is judged to be a descending section. In case of ascending, it is processed according to FIG. 2A, and in case of descending, it is processed according to FIG. 2B.

【0009】図3のピークを検出する場合を例にして説
明する。注目点(現在位置)がP0,P1,P2と移動し
てきて、P2から始まる区間では区間平均値が現在位置
2のデータより小さく、その区間の最小値の位置がP
sであったとする。次の区間では注目点(現在位置)が
Psに移動する。Psから始まる区間では、区間平均値
がPsよりも大きくなって上昇に転じたとすると、図2
(A)の手順に従って処理される。この場合、まだ(ピ
ーク最小位置)も(ピーク最大位置)も記録されていな
いので、(現在位置)Psが(ピーク最小位置)に記録
された後、(現在位置)には(区間最大値)の位置P3
が記録され、現在位置がその区間最大値の位置P3へ移
動する。これにより、変曲点Psが捉えられたことにな
る。
The case of detecting the peak in FIG. 3 will be described as an example. Point of interest (the current position) been moved P 0, P 1, P 2, interval average value in a section starting from P 2 is smaller than the data of the current position P 2, the position of the minimum value of the interval P
s. In the next section, the point of interest (current position) moves to Ps. In the section starting from Ps, assuming that the section average value becomes larger than Ps and starts to rise, FIG.
It is processed according to the procedure of (A). In this case, since neither (peak minimum position) nor (peak maximum position) has been recorded yet, after (current position) Ps is recorded in (peak minimum position), (current position) shows (section maximum value). Position P 3
Is recorded, and the current position moves to the position P 3 having the maximum value in the section. As a result, the inflection point Ps is captured.

【0010】再びメモリ領域(区間最小値)、(区間最
大値)及び(区間平均値)がリセットされた後、(現在
位置)P3から(現在位置P3+区間幅W)までの区間に
おける、区間最小値、区間最大値及び区間平均値を求め
られ、それぞれのメモリ領域に記録される。上昇区間が
続くかぎり、(現在位置)がその区間の(区間最大値)
の位置へ移動しながら処理が繰り返されていく。
After the memory area (section minimum value), (section maximum value) and (section average value) are reset again, in the section from (current position) P 3 to (current position P 3 + section width W) , Section minimum value, section maximum value and section average value are obtained and recorded in respective memory areas. As long as the climb section continues, (current position) is the (section maximum value) of that section.
The process is repeated while moving to the position of.

【0011】ピークの頂点Ptを含む上昇区間になる
と、注目点(現在位置)がPtに移動し、次の区間は
(現在位置)Ptから始まる区間幅Wの区間になる。そ
の区間では下降に転じ、図2(B)に従って処理され、
(現在位置)Ptが(ピーク最大位置)に記録される。
これにより、変曲点Ptが捉えられたことになる。その
後、(現在位置)にはその区間の(区間最小値)の位置
が記録され、現在位置がその区間最小値の位置へ移動す
る。そして、下降区間が続くかぎり、(現在位置)がそ
の区間の(区間最小値)の位置へ移動しながら処理が繰
り返されていく。
At the rising section including the peak Pt, the point of interest (current position) moves to Pt, and the next section becomes a section having a section width W starting from (current position) Pt. In that section, it starts to descend and is processed according to FIG.
(Current position) Pt is recorded in (peak maximum position).
As a result, the inflection point Pt is captured. After that, the position of (section minimum value) of the section is recorded in (current position), and the current position moves to the position of the section minimum value. Then, as long as the descending section continues, the process is repeated while the (current position) moves to the (section minimum value) position of the section.

【0012】やがて、Peを含む下降区間になると、注
目点(現在位置)がPeに移動し、次の区間は(現在位
置)Peから始まる区間幅Wの区間になる。その区間で
は上昇に転じ、図2(A)に従って処理される。この場
合、(ピーク最大位置)が記録されていて、かつ(ピー
ク最小位置)が(ピーク最大位置)より以前にあるとい
う条件が満たされるので、(ピーク最大位置)がピーク
頂点、(ピーク最小位置)がピーク開始点、(現在位
置)がピーク終了点であるピークとみなされてピークデ
ータが作成される。
When a descending section including Pe is finally reached, the point of interest (current position) moves to Pe, and the next section becomes a section having a section width W starting from (current position) Pe. In that section, it starts to rise and is processed according to FIG. In this case, since the condition that (peak maximum position) is recorded and (peak minimum position) is before (peak maximum position) is satisfied, (peak maximum position) is the peak apex and (peak minimum position). ) Is regarded as the peak starting point and (current position) is the peak ending point, and the peak data is created.

【0013】その後、(ピーク最大位置)の記録が消さ
れ、(ピーク最小位置)には(現在位置)Peが記録さ
れた後、(現在位置)にはその区間の(区間最大値)が
記録されて、現在位置が区間最大値の位置へ移動して次
の区間のデータ処理が開始される。前のピークの終了点
Peは次のピークの開始点Psとなり、同様にしてピー
ク検出が繰り返されていく。この条件で半値幅が約46
倍異なる2つのピークを同一区間幅で検出することがで
きた。
After that, the record of (peak maximum position) is erased, (current position) Pe is recorded in (peak minimum position), and (section maximum value) of that section is recorded in (current position). Then, the current position moves to the position having the maximum value in the section, and the data processing of the next section is started. The end point Pe of the previous peak becomes the start point Ps of the next peak, and the peak detection is repeated in the same manner. Under this condition, the half width is about 46
It was possible to detect two peaks having different times with the same interval width.

【0014】[0014]

【発明の効果】本発明ではクロマトグラムデータを1回
走査するだけでピークの開始点、頂点及び終了点を検出
することができるので、データ処理時間が短かくなる。
ピーク検出基準を正確に定めなくてもピークの検出が可
能であり、半値幅の異なる広い範囲のピークに対応する
ことができる。
According to the present invention, since the start point, the peak and the end point of the peak can be detected by scanning the chromatogram data once, the data processing time becomes short.
It is possible to detect peaks without accurately defining the peak detection reference, and it is possible to deal with a wide range of peaks having different half widths.

【図面の簡単な説明】[Brief description of drawings]

【図1】一実施例のピーク検出処理の全体を示すPAD
図である。
FIG. 1 is a PAD showing an overall peak detection process according to an embodiment.
It is a figure.

【図2】上昇区間(A)と下降区間(B)での処理を示
すPAD図である。
FIG. 2 is a PAD diagram showing processing in an ascending section (A) and a descending section (B).

【図3】一実施例の動作を示すクロマトグラムピークの
図である。
FIG. 3 is a chromatogram peak diagram showing the operation of one example.

Claims (1)

【特許請求の範囲】[Claims] 【請求項1】 クロマトグラムデータについてサンプリ
ングデータを複数個含む一定区間幅Wと、データの上昇
又は下降を判定するための連続区間値nを定め、クロマ
トグラム上に注目点を定め、その注目点から一定区間幅
W内のデータから区間最小値、区間最大値及び区間平均
値を求め、区間平均値が注目点の値より大きい場合はそ
の区間を上昇、小さい場合は下降と判定し、上昇区間で
は注目点をその区間の最大値の点へ移動させ、下降区間
では注目点をその区間の最小値の点へ移動させることに
よりピークの変曲点を捉え、n区間以上連続して下降し
た後の変曲点をピーク開始点もしくはピーク終了点と
し、n区間以上連続して上昇した後の変曲点をピーク頂
点とすることを特徴とするピーク検出方法。
1. A constant section width W containing a plurality of sampling data for chromatogram data, a continuous section value n for judging rise or fall of the data, and an attention point is set on the chromatogram. From the data within the fixed section width W, the section minimum value, the section maximum value and the section average value are obtained. If the section average value is larger than the value of the point of interest, it is determined that the section is up Then, move the point of interest to the point of maximum value in that section, and in the descending section, move the point of interest to the point of minimum value in that section to capture the inflection point of the peak, and after descending continuously for n sections or more The peak detection method is characterized in that the inflection point is defined as the peak start point or the peak end point, and the inflection point after continuously rising for n sections or more is defined as the peak apex.
JP36001792A 1992-12-30 1992-12-30 Method for detecting peak of chromatogram Pending JPH06201673A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP36001792A JPH06201673A (en) 1992-12-30 1992-12-30 Method for detecting peak of chromatogram

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP36001792A JPH06201673A (en) 1992-12-30 1992-12-30 Method for detecting peak of chromatogram

Publications (1)

Publication Number Publication Date
JPH06201673A true JPH06201673A (en) 1994-07-22

Family

ID=18467468

Family Applications (1)

Application Number Title Priority Date Filing Date
JP36001792A Pending JPH06201673A (en) 1992-12-30 1992-12-30 Method for detecting peak of chromatogram

Country Status (1)

Country Link
JP (1) JPH06201673A (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015056329A1 (en) * 2013-10-17 2015-04-23 株式会社島津製作所 Method and device for detecting peak endpoints in waveform
JPWO2017094170A1 (en) * 2015-12-03 2018-09-20 株式会社島津製作所 Peak detection method and data processing apparatus
WO2019092837A1 (en) * 2017-11-09 2019-05-16 富士通株式会社 Waveform analysis device
CN111595992A (en) * 2020-06-30 2020-08-28 浙江三青环保科技有限公司 Rapid peak searching method for online gas chromatographic peak
CN114047160A (en) * 2021-10-28 2022-02-15 中南大学 Second harmonic threading peak-searching method

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015056329A1 (en) * 2013-10-17 2015-04-23 株式会社島津製作所 Method and device for detecting peak endpoints in waveform
CN105637360A (en) * 2013-10-17 2016-06-01 株式会社岛津制作所 Method and device for detecting peak endpoints in waveform
JPWO2015056329A1 (en) * 2013-10-17 2017-03-09 株式会社島津製作所 Method and apparatus for detecting peak end point in waveform
JPWO2017094170A1 (en) * 2015-12-03 2018-09-20 株式会社島津製作所 Peak detection method and data processing apparatus
WO2019092837A1 (en) * 2017-11-09 2019-05-16 富士通株式会社 Waveform analysis device
CN111316095A (en) * 2017-11-09 2020-06-19 株式会社岛津制作所 Waveform analysis device
JPWO2019092837A1 (en) * 2017-11-09 2020-11-12 富士通株式会社 Waveform analyzer
CN111595992A (en) * 2020-06-30 2020-08-28 浙江三青环保科技有限公司 Rapid peak searching method for online gas chromatographic peak
CN114047160A (en) * 2021-10-28 2022-02-15 中南大学 Second harmonic threading peak-searching method

Similar Documents

Publication Publication Date Title
JP3378166B2 (en) Pulse signal classification device
JP5705395B2 (en) Signal analyzer
US20100299078A1 (en) Chromatography data processing method and system
US20050129292A1 (en) Method and arrangement for electronic recording of rolled fingerprints
JP5211753B2 (en) Chromatographic data processor
JPH06201673A (en) Method for detecting peak of chromatogram
JP3706031B2 (en) Target identification device
WO2005002209A1 (en) A method and a circuit for deriving a synchronisation signal from a video signal
EP0698860A2 (en) Edge line measuring method
JPH1172489A (en) Liquid chromatography
JP2005221276A (en) Data processor for chromatograph mass spectrometry
JP4826579B2 (en) Chromatographic data processor
JPH0835960A (en) Data processor for chromatograph mass analyzer
JP3498683B2 (en) Analysis equipment
JP2910290B2 (en) Chromatographic data processor
US7702158B2 (en) Image processing apparatus, image processing method and record medium for the image processing apparatus
JP2773296B2 (en) Chromatographic data processor
EP0348888A2 (en) Overflow speech detecting apparatus
CN112750458A (en) Touch screen sound detection method and device
EP2528019A1 (en) Apparatus and method for detecting objects in moving images
US7409075B2 (en) Image processing apparatus and method, program, and recording medium
CN111854603B (en) Laser detection method, device and equipment based on CCD sensor and storage medium
JPH06230001A (en) Data processing device for chromatograph
JPH04138327A (en) Processing method for measured spectrum data
JPH06148163A (en) Chromatographic data processing device