JP5737403B2 - Chromatographic data processor - Google Patents

Chromatographic data processor Download PDF

Info

Publication number
JP5737403B2
JP5737403B2 JP2013523875A JP2013523875A JP5737403B2 JP 5737403 B2 JP5737403 B2 JP 5737403B2 JP 2013523875 A JP2013523875 A JP 2013523875A JP 2013523875 A JP2013523875 A JP 2013523875A JP 5737403 B2 JP5737403 B2 JP 5737403B2
Authority
JP
Japan
Prior art keywords
peak
baseline
intersection
detected
width
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.)
Active
Application number
JP2013523875A
Other languages
Japanese (ja)
Other versions
JPWO2013008611A1 (en
Inventor
年伸 柳沢
年伸 柳沢
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 JP2013523875A priority Critical patent/JP5737403B2/en
Publication of JPWO2013008611A1 publication Critical patent/JPWO2013008611A1/en
Application granted granted Critical
Publication of JP5737403B2 publication Critical patent/JP5737403B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/86—Signal analysis
    • G01N30/8624—Detection of slopes or peaks; baseline correction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/86—Signal analysis
    • G01N30/8624—Detection of slopes or peaks; baseline correction
    • G01N30/8631—Peaks
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/86—Signal analysis
    • G01N30/8624—Detection of slopes or peaks; baseline correction
    • G01N30/8641—Baseline

Landscapes

  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Treatment Of Liquids With Adsorbents In General (AREA)

Description

本発明は、ガスクロマトグラフ、液体クロマトグラフ等のクロマトグラフ用のデータ処理装置に関する。   The present invention relates to a data processing apparatus for a chromatograph such as a gas chromatograph or a liquid chromatograph.

クロマトグラフ装置の性能や効率を判断する指標として理論段数や分離度などがある。このうち理論段数はカラムの分離性能を表す指標であり、クロマトグラム上の成分の保持時間とそのピーク幅とから求められる。   There are the number of theoretical plates and the degree of separation as indices for judging the performance and efficiency of the chromatographic apparatus. Of these, the theoretical plate number is an index representing the separation performance of the column, and is determined from the retention time of the components on the chromatogram and the peak width.

米国食品医薬品局（Food and Drug Administration: FDA）が管轄する米国薬局方（United States Pharmacopeia: USP）では、ピーク幅を次のように規定している（非特許文献１）。すなわち、図１に示すように、ピークトップＰを挟んでピークの前半部及び後半部でそれぞれの変曲点Ｃ１、Ｃ２を求め、この変曲点Ｃ１、Ｃ２においてピーク波形に接線を引く。この両接線とベースラインＢＬとの交点Ｂ１、Ｂ２を求め、両交点Ｂ１、Ｂ２間の距離をピーク幅Ｗとする。理論段数Ｎは、このピーク幅Ｗと保持時間Ｔｒとから、次の式(1)により求められる。
Ｎ＝16×（Ｔｒ／Ｗ）2 …(1)In the United States Pharmacopeia (USP), which is under the jurisdiction of the US Food and Drug Administration (FDA), the peak width is defined as follows (Non-patent Document 1). That is, as shown in FIG. 1, the inflection points C1 and C2 are obtained at the first half and the second half of the peak across the peak top P, and the peak waveform is drawn at the inflection points C1 and C2. The intersections B1 and B2 between the both tangents and the base line BL are obtained, and the distance between the intersections B1 and B2 is defined as the peak width W. The theoretical plate number N is obtained from the peak width W and the holding time Tr by the following equation (1).
N = 16 × (Tr / W) 2 (1)

また、厚生労働省が管轄する日本薬局方(Japanese Pharmacopeia: JP)では、図２に示すように、ピークトップＰから垂線を下ろし、ベースラインＢＬとの交点Ｑを求め、その高さＰＱの1/2となる点を通り、ベースラインＢＬに平行な直線とピークの前半部と後半部のそれぞれの交点Ｄ１、Ｄ２間の距離をピーク幅Ｗ0.5と規定している。この場合、理論段数Ｎは、次の式(2)により求められる（非特許文献２）。
Ｎ＝5.54×（Ｔｒ／Ｗ0.5）2 …(2)In addition, in the Japanese Pharmacopeia (JP) under the jurisdiction of the Ministry of Health, Labor and Welfare, as shown in FIG. The distance between the straight line parallel to the base line BL and the intersections D1 and D2 of the first half and the second half of the peak is defined as a peak width W 0.5 . In this case, the theoretical plate number N is obtained by the following equation (2) (Non-Patent Document 2).
N = 5.54 × (Tr / W 0.5 ) 2 (2)

式(1)及び(2)により求められる理論段数Ｎは、クロマトグラムのピーク形状が理想的なガウス分布（正規分布）である場合には一致する。ピーク幅Ｗ、Ｗ0.5は、理論段数以外にも例えば分離度を求める際にも用いられる。また、ベースラインＢＬから5%、10%の高さのピーク幅Ｗ0.05、Ｗ0.1を算出し、シンメトリー係数等の他の指標を求めることもある。
なお、日本薬局方によるピーク幅の算出方法は、上記のようにベースラインＢＬから50%の高さにおいて行うものであるが、5%や10%など異なる高さのピーク幅を算出する場合も同様の手順で行うことができる。そのため、以下では、これらをまとめて「日本薬局方のピーク幅算出方法」と呼ぶことにする。The theoretical plate number N obtained from the equations (1) and (2) is the same when the peak shape of the chromatogram is an ideal Gaussian distribution (normal distribution). The peak widths W and W 0.5 are used not only for the number of theoretical plates but also for obtaining the degree of separation, for example. In addition, peak widths W 0.05 and W 0.1 having a height of 5% and 10% from the baseline BL may be calculated to obtain other indices such as symmetry coefficients.
The method for calculating the peak width by the Japanese Pharmacopoeia is performed at a height of 50% from the baseline BL as described above. However, there are cases where peak widths of different heights such as 5% and 10% are calculated. The same procedure can be used. Therefore, in the following, these will be collectively referred to as the “Japanese Pharmacopoeia peak width calculation method”.

“レビュワー・ガイダンス（Reviewer Guidance） バリデイション・オブ・クロマトグラフィック・メソッズ（Validation of Chromatographic Methods）”、［Online］、センター・フォー・ドラッグ・エバリュエーション・アンド・リサーチ（Center for Drug Evaluation and Research (CDER)）、［平成２４年６月１１日検索］、インターネット〈URL: http://www.fda.gov/downloads/Drugs/GuidanceComplianceRegulatoryInformation/Guidances/ucm072974.pdf〉“Reviewer Guidance Validation of Chromatographic Methods” [Online], Center for Drug Evaluation and Research ( CDER)), [Searched on June 11, 2012], Internet <URL: http://www.fda.gov/downloads/Drugs/GuidanceComplianceRegulatoryInformation/Guidances/ucm072974.pdf> “第十五改正日本薬局方”、［Online］、厚生労働省、［平成２４年６月１１日検索］、インターネット〈URL: http://jpdb.nihs.go.jp/jp15/YAKKYOKUHOU15.pdf〉“15th Japanese Pharmacopoeia”, [Online], Ministry of Health, Labor and Welfare, [Searched on June 11, 2012], Internet <URL: http://jpdb.nihs.go.jp/jp15/YAKKYOKUHOU15.pdf>

クロマトグラフ用のデータ処理装置では、上記のように、得られたクロマトグラムの理論段数や分離度等の指標を計算するために、ピーク幅の算出が必要となることが多い。しかしながら、図３及び図４に示すように隣り合うピークが重なる場合、実際と異なるピーク幅を算出してしまったり、ピーク幅を算出できなくなったりすることがある。   In a chromatograph data processing apparatus, as described above, it is often necessary to calculate a peak width in order to calculate an index such as the number of theoretical plates and the resolution of the obtained chromatogram. However, when adjacent peaks overlap as shown in FIGS. 3 and 4, a peak width different from the actual peak width may be calculated, or the peak width may not be calculated.

図３は、図１に示した米国薬局方のピーク幅算出方法によりピーク幅Ｗを算出する例を示したものである。この図の隣り合う２つのピークのうち、左側のピーク１については前半部と後半部の両方で適切に変曲点が得られるため、ピーク幅Ｗを正しく算出することができる。一方、右側のピーク２については、その前半部がピーク１と重なっていることにより、前半部における本来の変曲点の位置とは異なる位置に変曲点が検出されている。このような場合、ピーク２のピーク幅Ｗを正しく算出することができない。   FIG. 3 shows an example in which the peak width W is calculated by the US Pharmacopoeia peak width calculation method shown in FIG. Of the two adjacent peaks in this figure, for the left peak 1, inflection points are appropriately obtained in both the first half and the second half, so that the peak width W can be calculated correctly. On the other hand, for the peak 2 on the right side, the inflection point is detected at a position different from the position of the original inflection point in the first half because the first half overlaps the peak 1. In such a case, the peak width W of peak 2 cannot be calculated correctly.

図４は、日本薬局方のピーク幅算出方法によりピーク幅Ｗ0.5を算出する例を示したものである。この図の２つのピークのうち、ピーク１については、ベースラインに平行な直線とピーク１との交点が前半部と後半部の両方で得られるため、ピーク幅Ｗ0.5を算出することができる。一方、ピーク２については、その前半部がピーク１と重なっていることにより、前半部では交点を得ることができず、ピーク幅Ｗ0.5を算出することができない。FIG. 4 shows an example in which the peak width W 0.5 is calculated by the Japanese Pharmacopoeia peak width calculation method. Of the two peaks in this figure, for peak 1, the intersection point between the straight line parallel to the base line and peak 1 is obtained in both the first half and the second half, so the peak width W 0.5 can be calculated. On the other hand, since the first half of peak 2 overlaps with peak 1, no intersection can be obtained in the first half and the peak width W 0.5 cannot be calculated.

以上の場合、デコンボリューション演算によりピークを分離してからピーク幅を算出することも可能ではあるが、分離のための処理に時間を要すると共に、求める値は推定値の域を越えない。   In the above case, although it is possible to calculate the peak width after separating the peaks by the deconvolution calculation, it takes time for the separation processing, and the obtained value does not exceed the range of the estimated value.

本発明が解決しようとする課題は、ピークが重なることにより、ピークの前半部と後半部のいずれか一方の変曲点又は交点が適切に得られない場合でも、複雑な処理を行うことなく、より確かなピーク幅を算出することができるクロマトグラフ用データ処理装置を提供することである。   The problem to be solved by the present invention is that even if inflection points or intersections of any one of the first half and the second half of the peak cannot be obtained appropriately due to overlapping peaks, complicated processing is not performed. It is an object of the present invention to provide a chromatographic data processing apparatus capable of calculating a more reliable peak width.

上記課題を解決するために成された本発明に係るクロマトグラフ用データ処理装置の第一の態様は、
a) クロマトグラムのベースラインを決定するベースライン決定手段と、
b) 前記クロマトグラムのピークトップを検出するピークトップ検出手段と、
c) 前記ピークトップを挟んでピークの前半部又は後半部において変曲点を検出し、該変曲点における接線と前記ベースラインとの交点と、該ピークトップから該ベースラインに下ろした垂線と該ベースラインとの交点と、をそれぞれ検出し、これら両交点間の距離を求め、その距離の２倍の値を該ピークのピーク幅として算出するピーク幅算出手段と、
を備えることを特徴とする。The first aspect of the chromatograph data processing apparatus according to the present invention made to solve the above problems,
a) a baseline determination means for determining a baseline of the chromatogram;
b) Peak top detecting means for detecting the peak top of the chromatogram;
c) detecting an inflection point in the first half or the latter half of the peak across the peak top, an intersection of the tangent line at the inflection point and the base line, and a perpendicular line extending from the peak top to the baseline. A peak width calculating means for detecting an intersection with the baseline, obtaining a distance between the two intersections, and calculating a value twice as large as the peak width of the peak;
It is characterized by providing.

また、上記課題を解決するために成された本発明に係るクロマトグラフ用データ処理装置の第二の態様は、
a) クロマトグラムのベースラインを決定するベースライン決定手段と、
b) 前記クロマトグラムのピークトップを検出するピークトップ検出手段と、
c) 前記ベースラインからピークトップまでの高さのＭ%の高さに該ベースラインに平行な直線を引き、該ピークトップを挟んでピークの前半部又は後半部における該直線との交点と、該ピークトップから該直線に下ろした垂線と該直線との交点と、をそれぞれ検出し、これら両交点間の距離を求め、その距離の２倍の値を該ピークのピーク幅として算出するピーク幅算出手段と、
を備えることを特徴とする。In addition, the second aspect of the chromatographic data processing apparatus according to the present invention made to solve the above problems,
a) a baseline determination means for determining a baseline of the chromatogram;
b) Peak top detecting means for detecting the peak top of the chromatogram;
c) A straight line parallel to the baseline is drawn at a height of M% of the height from the baseline to the peak top, and the intersection with the straight line in the first half or the second half of the peak across the peak top; A peak width that detects a perpendicular line drawn from the peak top to the straight line and an intersection of the straight line, obtains a distance between the two intersection points, and calculates a value twice the distance as the peak width of the peak A calculation means;
It is characterized by providing.

クロマトグラムにおける理論段数等の指標は、元々クロマトグラムのピーク形状がガウス分布に従うものとして、計算式が作られている。クロマトグラムのピーク形状が理想的なガウス分布である場合、このピークはピークトップを挟んで前半部と後半部で対称になる。本発明は、この対称性を仮定することにより、デコンボリューション演算などによるピークの分離を行わずに、ピークの前半部と後半部のいずれか一方の算出可能な幅の２倍で以てピーク幅を算出するものである。   Indices such as the number of theoretical plates in the chromatogram are calculated by assuming that the peak shape of the chromatogram originally follows a Gaussian distribution. When the peak shape of the chromatogram is an ideal Gaussian distribution, this peak is symmetric between the first half and the second half with the peak top in between. By assuming this symmetry, the present invention does not perform peak separation by a deconvolution operation or the like, and the peak width is doubled by a width that can be calculated from either the first half or the second half of the peak. Is calculated.

なお、ピーク同士が重なっていない場合や、重なっていたとしてもそれが図３や図４のピーク１のようにピーク幅の算出に影響しない場合には、従来通り、ピークの前半部及び後半部の各々の変曲点における接線とベースラインとの交点間の距離、又はベースラインに平行な直線とピークとの２つの交点間の距離をピーク幅として算出することが望ましい。   If the peaks do not overlap or if they do not affect the calculation of the peak width as in peak 1 of FIG. 3 or FIG. 4, the first half and the second half of the peak as before. It is desirable to calculate the distance between the intersections of the tangent and the base line at each inflection point or the distance between the two intersections of the straight line and the peak parallel to the baseline as the peak width.

本発明に係るクロマトグラフ用データ処理装置によれば、目的とするピークの前半部と後半部の幅のうち一方を求めることさえできれば、そのピーク幅を容易に算出することができる。また、ピーク形状の対称性が保たれている場合には、デコンボリューション演算によりピークを分離するよりも確かなピーク幅を算出することができる。   According to the chromatographic data processing apparatus of the present invention, the peak width can be easily calculated as long as one of the widths of the first half and the second half of the target peak can be obtained. In addition, when the symmetry of the peak shape is maintained, a more reliable peak width can be calculated than when the peaks are separated by deconvolution calculation.

米国薬局方のピーク幅算出法の説明図。Explanatory drawing of the peak width calculation method of the US Pharmacopoeia. 日本薬局方のピーク幅算出法の説明図。Explanatory drawing of the peak width calculation method of the Japanese Pharmacopoeia. 米国薬局方のピーク幅算出法により、実際と異なるピーク幅を算出してしまうピークの例を示す図。The figure which shows the example of the peak which calculates the peak width different from actual by the peak width calculation method of the United States Pharmacopeia. 日本薬局方のピーク幅算出法では、ピーク幅を算出することができないピークの例を示す図。The figure which shows the example of the peak which cannot calculate a peak width by the Japanese Pharmacopoeia peak width calculation method. 本発明に係るデータ処理装置の一実施例を備えたガスクロマトグラフ分析システムの概略構成図。BRIEF DESCRIPTION OF THE DRAWINGS The schematic block diagram of the gas chromatograph analysis system provided with one Example of the data processor concerning this invention. 理論段数算出の処理手順を示すフローチャート。The flowchart which shows the process sequence of theoretical plate number calculation. 本実施例による第１のピーク幅算出の処理手順を示すフローチャート。The flowchart which shows the process sequence of the 1st peak width calculation by a present Example. 第１のピーク幅算出処理の説明図。Explanatory drawing of a 1st peak width calculation process. 本実施例による第２のピーク幅算出の処理手順を示すフローチャート。The flowchart which shows the process sequence of the 2nd peak width calculation by a present Example. 第２のピーク幅算出処理の説明図。Explanatory drawing of a 2nd peak width calculation process.

以下、本発明に係るクロマトグラフ用データ処理装置の一実施例について、図５〜図１０を参照して説明する。   Hereinafter, an embodiment of a chromatographic data processing apparatus according to the present invention will be described with reference to FIGS.

図５は、本実施例によるデータ処理装置を備えるガスクロマトグラフ分析システムの概略構成図である。液体試料はシリンジ１１等により試料気化室１２に注入され、そこで気化されて、キャリアガス流路１３から一定流量で供給されるキャリアガス流に乗ってカラム１４内に送り込まれる。試料に含まれる各種成分は、カラム１４を通過する間に時間的に分離されてカラム１４から出て、順次、検出器１５により検出される。検出器１５による検出信号はデジタルデータに変換された後、逐次、データ処理装置１６に送られ、そこで一旦、ハードディスク等の記憶部に格納される。１つの試料の分析が終了した後（或いは、連続に実行される複数の試料の分析が終了した後）、記憶部に格納されたデータが読み出され、クロマトグラムの作成、ピーク検出等の各種のデータ処理が行われる。   FIG. 5 is a schematic configuration diagram of a gas chromatograph analysis system including a data processing apparatus according to the present embodiment. The liquid sample is injected into the sample vaporizing chamber 12 by a syringe 11 or the like, vaporized there, and sent into the column 14 on the carrier gas flow supplied from the carrier gas channel 13 at a constant flow rate. Various components contained in the sample are temporally separated while passing through the column 14, exit the column 14, and are sequentially detected by the detector 15. The detection signal from the detector 15 is converted into digital data and then sequentially sent to the data processing device 16 where it is temporarily stored in a storage unit such as a hard disk. After the analysis of one sample is completed (or after the analysis of a plurality of samples executed in succession is completed), the data stored in the storage unit is read out, and various data such as creation of chromatogram, peak detection, etc. Data processing is performed.

なお、データ処理装置１６の実体は専用又は汎用のコンピュータであり、所定の処理プログラムを動作させることにより、当該コンピュータをクロマトグラム作成部１６１、ベースライン決定部１６２、ピークトップ検出部１６３、ピーク幅算出部１６４として機能させる他、各種分析のためのデータ処理を実行させる。
The entity of the data processing device 16 is a dedicated or general-purpose computer. By operating a predetermined processing program, the computer is made to have a chromatogram creation unit 161, a baseline determination unit 162, a peak top detection unit 163, a peak width. In addition to functioning as the calculation unit 164, data processing for various types of analysis is executed.

上記データ処理装置１６において理論段数Ｎを算出する際の処理手順を、図６のフローチャートに示す。まずデータ処理装置１６内の記憶部から読み出したクロマトグラムデータに対し、クロマトグラム作成部１６１がクロマトグラムを作成し、ベースライン決定部１６２及びピークトップ検出部１６３が所定のアルゴリズムに基づいて、作成したクロマトグラムのベースライン及びピーク（単数又は複数）を決定・検出する（ステップＳ１、Ｓ２）。そして、各ピーク又は必要なピークのみに関し、ピーク開始時間Ｔｓ、ピークトップ時間Ｔｐ、ピーク終了時間Ｔｅ、ピーク高さＨｐ、ピーク面積Ｓ等、ピークを特徴付けるパラメータを算出する（ステップＳ３）。さらに、上記パラメータとクロマトグラムデータとを利用して、ピーク幅算出部１６４がピーク幅を算出し（ステップＳ４）、ピーク幅の算出方法に対応した式(1)や式(2)等の式により理論段数Ｎを求める（ステップＳ５）。
The processing procedure for calculating the theoretical plate number N in the data processing device 16 is shown in the flowchart of FIG. First, the chromatogram creation unit 161 creates a chromatogram for the chromatogram data read from the storage unit in the data processing device 16, and the baseline determination unit 162 and the peak top detection unit 163 create based on a predetermined algorithm. The baseline and peak (s) of the obtained chromatogram are determined and detected (steps S1, S2). Then, with respect to each peak or only necessary peaks, parameters characterizing the peaks such as peak start time Ts, peak top time Tp, peak end time Te, peak height Hp, peak area S, etc. are calculated (step S3). Further, using the above parameters and chromatogram data, the peak width calculation unit 164 calculates the peak width (step S4), and formulas such as formula (1) and formula (2) corresponding to the peak width calculation method. Thus, the theoretical plate number N is obtained (step S5).

本発明に係るデータ処理の特徴はステップＳ４のピーク幅算出処理にあるので、この処理に関し詳細に説明する。なお、ピーク幅の算出方法は、上記のように、大別して米国薬局方の方法と日本薬局方の方法との２種類があるが、まず米国薬局方の方法に準じてピーク幅Ｗを算出する手順を示す。   Since the feature of the data processing according to the present invention is the peak width calculation processing in step S4, this processing will be described in detail. As described above, there are two types of peak width calculation methods, the US Pharmacopeia method and the Japanese Pharmacopeia method. First, the peak width W is calculated according to the US Pharmacopeia method. Show the procedure.

［第１のピーク幅算出処理］
図７は、本実施例における第１のピーク幅算出処理のフローチャートである。以下、図１及び図８の説明図を参照しつつ、この第１のピーク幅算出処理の手順を説明する。[First peak width calculation process]
FIG. 7 is a flowchart of the first peak width calculation process in the present embodiment. Hereinafter, the procedure of the first peak width calculation process will be described with reference to the explanatory diagrams of FIGS. 1 and 8.

本実施例における第１のピーク幅算出処理では、まず目的とするピークに対し、そのピーク開始時間からピーク終了時間までの間で、ピークトップＰを挟んで前半部と後半部における変曲点を探索する（ステップＳ１１）。なお、変曲点を探索する一般的なアルゴリズムは、次の通りである。   In the first peak width calculation process in the present embodiment, first, inflection points in the first half and the second half of the target peak are sandwiched between the peak start time and the peak end time with the peak top P interposed therebetween. Search is performed (step S11). A general algorithm for searching for an inflection point is as follows.

ピークトップＰから時間的にマイナス側に移動しつつ、ピーク波形上の各点において２次微分値及び１次微分値を算出し、２次微分値が０であって且つ１次微分値が０より大（つまり正値である）であるか否かを判定する。そして、この条件が満たされる際のクロマトグラム上の点を前半部の変曲点Ｃ１として採用する。後半部の変曲点Ｃ２については、ピークトップＰから時間的にプラス側に移動する点と、１次微分値に関する判定条件が０より小（つまり負値である）となる点以外は同様である。なお、ピークトップＰからではなく、ピークの裾野（ピーク開始点及びピーク終了点）から変曲点を探索する方法もある。また、２次微分値及び１次微分値の算出には、例えばサビツキ・ゴーレイ（Savitzky-Golay）の方法など既知のアルゴリズムを用いることができる。   The secondary differential value and the primary differential value are calculated at each point on the peak waveform while moving from the peak top P to the negative side in time. The secondary differential value is 0 and the primary differential value is 0. It is determined whether it is larger (that is, a positive value). And the point on the chromatogram when this condition is satisfied is adopted as the inflection point C1 in the first half. The inflection point C2 in the latter half is the same except that it moves from the peak top P to the plus side in time and the determination condition for the first derivative is smaller than 0 (that is, a negative value). is there. There is also a method of searching for an inflection point not from the peak top P but from the peak base (peak start point and peak end point). For calculating the secondary differential value and the primary differential value, for example, a known algorithm such as a Savitzky-Golay method can be used.

次に、変曲点がピークの前半部と後半部の両方で適切に検出されたか否かを判定する（ステップＳ１２）。なお、ピーク形状が理想的なガウス分布に従う場合、変曲点はピーク高さＨｐの1/e0.5（eは自然対数の底）倍の高さになることが分かっている（特開2004-184148号公報を参照）。従って、例えば検出された変曲点のベースラインからの高さが1/e0.5Ｈｐから所定範囲内に収まっているか否かにより、該変曲点が適切であるか否かを判定することができる。Next, it is determined whether or not the inflection point is properly detected in both the first half and the second half of the peak (step S12). When the peak shape follows an ideal Gaussian distribution, the inflection point is known to be 1 / e 0.5 (e is the base of natural logarithm) times the peak height Hp (Japanese Patent Laid-Open No. 2004-2004). 184148). Therefore, for example, whether or not the inflection point is appropriate can be determined based on whether or not the height from the baseline of the detected inflection point is within a predetermined range from 1 / e 0.5 Hp. it can.

ステップＳ１２において、図１のように変曲点Ｃ１、Ｃ２がピークの前半部と後半部でそれぞれ適切に検出された場合は、各々の変曲点Ｃ１、Ｃ２において１次微分値を用いて接線を引き、該接線とベースラインとの交点Ｂ１、Ｂ２を検出する（ステップＳ１６）。そして、両交点Ｂ１、Ｂ２間の距離をピーク幅Ｗとして採用する（ステップＳ１７）。これらは従来通りの算出処理である。
In step S12, when the inflection points C1 and C2 are appropriately detected in the first half and the latter half of the peak as shown in FIG. 1, the tangent lines are obtained using the primary differential values at the inflection points C1 and C2. And the intersection points B1 and B2 between the tangent line and the base line are detected (step S1 6 ). Then, the distance between the two intersections B1 and B2 is adopted as the peak width W (step S1 7 ). These are conventional calculation processes.

一方、図８のように、ピークの前半部と後半部のいずれかの変曲点が適切に検出されない場合、このままではピーク幅を正しく算出することができない。このような場合、本実施例における第１のピーク幅算出処理では、検出された方の変曲点（図８の例では点Ｃ２）に接線を引き、ベースラインＢＬとの交点Ｂ（又は交点Ｂ２）を検出すると共に、ピークトップＰからベースラインＢＬに下ろした垂線とベースラインＢＬとの交点Ｑを検出する（ステップＳ１３、Ｓ１４）。そして、ステップＳ１３及びＳ１４で検出した２つの交点Ｂ、Ｑ間の距離を求め、この距離の２倍の値をピーク幅Ｗとして採用する（ステップＳ１５）。
これらステップＳ１３〜Ｓ１５の処理により、変曲点の一方が適切に検出されない場合であっても、ピーク幅Ｗを算出することが可能となる。On the other hand, as shown in FIG. 8, when any inflection point in the first half or the second half of the peak is not properly detected, the peak width cannot be calculated correctly as it is. In such a case, in the first peak width calculation process in the present embodiment, a tangent line is drawn to the detected inflection point (point C2 in the example of FIG. 8), and the intersection B (or intersection) with the base line BL. B2) is detected, and an intersection point Q between the perpendicular line dropped from the peak top P to the base line BL and the base line BL is detected (steps S13 and S14). Then, the distance between the two intersections B and Q detected in steps S13 and S14 is obtained, and a value twice the distance is adopted as the peak width W (step S15).
Through the processes in steps S13 to S15, the peak width W can be calculated even when one of the inflection points is not properly detected.

［第２のピーク幅算出処理］
次に、日本薬局方の方法に準じてピーク高さの50%の位置におけるピーク幅Ｗ0.5を算出する処理の手順を、図２及び図１０の説明図を参照しつつ、図９のフローチャートを用いて説明する。[Second peak width calculation process]
Next, the processing procedure for calculating the peak width W 0.5 at the position of 50% of the peak height according to the method of the Japanese Pharmacopoeia is shown in the flowchart of FIG. 9 with reference to the explanatory diagrams of FIGS. It explains using.

本実施例における第２のピーク幅算出処理では、まず目的とするピークのベースラインＢＬからピークトップＰまでの高さの50%の高さにベースラインＢＬに平行な直線ＰＬを引き、該ピークとの交点を検出する（ステップＳ２１）。そして、ステップＳ２１で検出した交点が該ピークの前半部と後半部の両方において得られたか否かを判定する（ステップＳ２２）。   In the second peak width calculation process in this embodiment, first, a straight line PL parallel to the base line BL is drawn to a height of 50% of the height from the base line BL of the target peak to the peak top P, and the peak Is detected (step S21). Then, it is determined whether or not the intersection detected in step S21 is obtained in both the first half and the second half of the peak (step S22).

ステップＳ２２において、ピークの前半部と後半部の両方において直線ＰＬとの交点が存在する場合は、図２に示すように各々の交点Ｄ１、Ｄ２の間の距離を求めることで、ピーク幅Ｗ0.5を算出する（ステップＳ２５）。In step S22, if there is a point of intersection between the straight line PL in both the first half and the second half of the peak, by obtaining the distance between each of the intersection points D1, D2 as shown in FIG. 2, peak width W 0.5 Is calculated (step S25).

一方、図１０のように、ピークの前半部と後半部のいずれかの交点が得られない場合、本実施例における第２のピーク幅算出処理では、ピークトップＰから直線ＰＬに下ろした垂線と該直線Ｓとの交点Ｒを求める（ステップＳ２３）。そして、ステップＳ２１で検出された方の交点（図１０の例では点Ｄ）と交点Ｒの間の距離を算出し、この距離の２倍の値をピーク幅Ｗ0.5として採用する（ステップＳ２４）。
これらステップＳ２３及びＳ２４における処理により、目的とするピークと直線ＰＬとの交点のうちの一方が検出されない場合であっても、ピーク幅Ｗ0.5を算出することが可能となる。On the other hand, as shown in FIG. 10, when the intersection of either the first half or the second half of the peak is not obtained, the second peak width calculation process in the present embodiment uses a perpendicular line dropped from the peak top P to the straight line PL. An intersection R with the straight line S is obtained (step S23). Then, the distance between the intersection (point D in the example of FIG. 10) detected at step S21 and the intersection R is calculated, and a value twice this distance is adopted as the peak width W 0.5 (step S24). .
By the processing in steps S23 and S24, the peak width W0.5 can be calculated even when one of the intersections between the target peak and the straight line PL is not detected.

なお、上記実施例は本発明の一例にすぎないから、上記記載の点以外についても、本発明の趣旨の範囲で適宜に変更や修正を加えることができることは明らかである。例えば、本実施例では第２のピーク幅算出処理を、ベースラインからピークトップまでの高さの50%の高さにおけるピーク幅を求める場合について説明したが、より一般的にＭ%（0＜Ｍ＜100）としても同様の手順により算出することができる。   Since the above embodiment is merely an example of the present invention, it is obvious that changes and modifications can be made as appropriate within the scope of the present invention except for the points described above. For example, in the present embodiment, the second peak width calculation process has been described for the case where the peak width at the height of 50% of the height from the baseline to the peak top is obtained, but more generally M% (0 < M <100) can be calculated by the same procedure.

また、特開2004-184148号公報では、ベースラインからピーク高さの1/e0.5倍の高さにおけるピーク上の点を仮想変曲点として定め、ピークの前半部及び後半部の仮想変曲点における接線とベースラインとの２つの交点間の距離をピーク幅Ｗとして算出する方法を示しているが、このような算出方法に対しては、仮想変曲点の検出及び検出した仮想変曲点の判定については図９のステップＳ２１、Ｓ２２と同様の処理で、判定後のピーク幅の算出については図７のステップＳ１３〜Ｓ１７と同様の処理で、それぞれ行えば良い。Further, in Japanese Patent Laid-Open No. 2004-184148, a point on the peak at a height 1 / e 0.5 times the peak height from the baseline is defined as a virtual inflection point, and virtual inflections in the first half and the latter half of the peak A method of calculating the distance between two intersections of a tangent and a baseline at a point as a peak width W is shown. For such a calculation method, detection of a virtual inflection point and detection of a virtual inflection detected The determination of the point may be performed by the same processing as steps S21 and S22 in FIG. 9, and the calculation of the peak width after the determination may be performed by the processing similar to steps S13 to S17 of FIG.

１１…シリンジ
１２…試料気化室
１３…キャリアガス流路
１４…カラム
１５…検出器
１６…データ処理装置
１６１…クロマトグラム作成部
１６２…ベースライン決定部
１６３…ピークトップ検出部
１６４…ピーク幅算出部DESCRIPTION OF SYMBOLS 11 ... Syringe 12 ... Sample vaporization chamber 13 ... Carrier gas flow path 14 ... Column 15 ... Detector 16 ... Data processing device 161 ... Chromatogram preparation part 162 ... Baseline determination part 163 ... Peak top detection part 164 ... Peak width calculation part

Claims (3)

a) クロマトグラムのベースラインを決定するベースライン決定手段と、
b) 前記クロマトグラムのピークトップを検出するピークトップ検出手段と、
c) 前記ピークトップを挟んでピークの前半部及び後半部において変曲点を検出するとともに該検出した変曲点がピーク本来のものであるかを判定し、ピークの前半部と後半部の両方において前記変曲点がピーク本来のものであると判定したときは該２つの変曲点のそれぞれにおける接線と前記ベースラインとの交点を検出して、これら両交点間の距離をピーク幅として算出し、前記ピークの前半部及び後半部のいずれかにおける変曲点がピーク本来のものであると判定したときは、該ピーク本来のものであると判定した変曲点における接線と前記ベースラインとの交点と、該ピークトップから該ベースラインに下ろした垂線と該ベースラインとの交点と、をそれぞれ検出して、これら両交点間の距離を求め、その距離の２倍の値を該ピークのピーク幅として算出するピーク幅算出手段と、
を備えることを特徴とするクロマトグラフ用データ処理装置。 a) a baseline determination means for determining a baseline of the chromatogram;
b) Peak top detecting means for detecting the peak top of the chromatogram;
c) Inflection points are detected in the first half and the second half of the peak across the peak top, and it is determined whether the detected inflection point is the original peak, and both the first half and the second half of the peak are detected. When it is determined that the inflection point is the original peak, the intersection between the tangent and the base line at each of the two inflection points is detected, and the distance between the two intersections is calculated as the peak width. When the inflection point in either the first half or the second half of the peak is determined to be the original peak , the tangent at the inflection point determined to be the original peak and the baseline And the intersection of the base line and the perpendicular line drawn down from the peak top to the baseline, respectively, to determine the distance between the two intersection points, and to obtain a value twice the distance of the peak And peak width calculating means for calculating as over click width,
A chromatograph data processing apparatus comprising: 前記ピーク幅算出手段が、検出した変曲点がピーク本来のものであるか否かの判定を、前記ベースラインから前記ピークトップまでの高さの1/e0.5（eは自然対数の底）倍の高さから所定範囲内の高さに変曲点が存在するか否かにより行うことを特徴とする請求項１に記載のクロマトグラフ用データ処理装置。 The peak width calculation means determines whether or not the detected inflection point is the original peak, 1 / e 0.5 of the height from the baseline to the peak top (e is the base of natural logarithm) 2. The chromatograph data processing apparatus according to claim 1 , wherein the chromatographic data processing apparatus is performed based on whether or not an inflection point exists from a double height to a height within a predetermined range. a) クロマトグラムのベースラインを決定するベースライン決定手段と、
b) 前記クロマトグラムのピークトップを検出するピークトップ検出手段と、
c) 前記ベースラインからピークトップまでの高さのＭ%の高さに該ベースラインに平行な直線を引き、該ピークトップを挟んでピークの前半部と後半部の両方において該ピークと前記直線との交点を検出可能な場合に、前記ピーク幅算出手段がこれら２つの交点間の距離をピーク幅として算出し、前記ピークの前半部又は後半部のいずれかにおける該直線と該ピークとの交点を検出可能な場合に該交点と、該ピークトップから該直線に下ろした垂線と該直線との交点と、をそれぞれ検出し、これら両交点間の距離を求め、その距離の２倍の値を該ピークのピーク幅として算出するピーク幅算出手段と、
を備えることを特徴とするクロマトグラフ用データ処理装置。 a) a baseline determination means for determining a baseline of the chromatogram;
b) Peak top detecting means for detecting the peak top of the chromatogram;
c) A straight line parallel to the baseline is drawn at a height of M% of the height from the baseline to the peak top , and the peak and the straight line are both in the first half and the second half of the peak across the peak top. When the intersection with the peak can be detected, the peak width calculation means calculates the distance between the two intersections as the peak width, and the intersection of the straight line and the peak in either the first half or the second half of the peak Is detected , the intersection between the intersection and the perpendicular line drawn down from the peak top to the straight line and the straight line is obtained, the distance between the two intersections is obtained, and a value twice the distance is obtained. A peak width calculating means for calculating the peak width of the peak;
A chromatograph data processing apparatus comprising:
JP2013523875A 2011-07-08 2012-06-22 Chromatographic data processor Active JP5737403B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2013523875A JP5737403B2 (en) 2011-07-08 2012-06-22 Chromatographic data processor

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2011151433 2011-07-08
JP2011151433 2011-07-08
JP2013523875A JP5737403B2 (en) 2011-07-08 2012-06-22 Chromatographic data processor
PCT/JP2012/066068 WO2013008611A1 (en) 2011-07-08 2012-06-22 Data processing device for chromatograph

Publications (2)

Publication Number Publication Date
JPWO2013008611A1 JPWO2013008611A1 (en) 2015-02-23
JP5737403B2 true JP5737403B2 (en) 2015-06-17

Family

ID=47505908

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013523875A Active JP5737403B2 (en) 2011-07-08 2012-06-22 Chromatographic data processor

Country Status (3)

Country Link
US (1) US20140149050A1 (en)
JP (1) JP5737403B2 (en)
WO (1) WO2013008611A1 (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103969379B (en) * 2013-01-29 2017-07-18 深圳普门科技有限公司 The measured value adquisitiones that a kind of liquid chromatogram is examined
JP6379463B2 (en) * 2013-09-18 2018-08-29 株式会社島津製作所 Waveform processing support method and waveform processing support device
GB2520684B (en) * 2013-11-27 2020-05-20 Micromass Ltd Method for measurement of ion events
JP2015224895A (en) * 2014-05-26 2015-12-14 有限会社シマムラテック Recycle isolation method in liquid chromatography
JP6540424B2 (en) * 2015-09-24 2019-07-10 富士通株式会社 Estimation device, estimation method, estimation program, engine and mobile device
WO2018025366A1 (en) * 2016-08-04 2018-02-08 株式会社島津製作所 Chromatograph input device and chromatograph
WO2018098502A1 (en) * 2016-11-28 2018-05-31 Board Of Regents, The University Of Texas System Systems, methods and devices for width-based analysis of peak traces
GB201803129D0 (en) 2018-02-27 2018-04-11 Ge Healthcare Bio Sciences Ab Method in liquid chiromatorgraphy
CN113607867B (en) * 2021-07-23 2024-06-11 清华大学合肥公共安全研究院 Double-overlap spectrum peak analysis method based on peak body mapping

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0395481A3 (en) * 1989-04-25 1991-03-20 Spectra-Physics, Inc. Method and apparatus for estimation of parameters describing chromatographic peaks
JPH0798270A (en) * 1993-09-29 1995-04-11 Shimadzu Corp Flow-through-type ultraviolet visible spectrophotometer
US6019896A (en) * 1998-03-06 2000-02-01 Molecular Dynamics, Inc. Method for using a quality metric to assess the quality of biochemical separations
JP3554181B2 (en) * 1998-03-31 2004-08-18 京セラ株式会社 Fluid separation device and method of forming the same
JP3702831B2 (en) * 2001-10-18 2005-10-05 株式会社島津製作所 Chromatographic data processor
JP3985671B2 (en) * 2002-12-02 2007-10-03 株式会社島津製作所 Chromatographic data processor
EP1915717A1 (en) * 2005-08-19 2008-04-30 University Of Tennessee Research Foundation Method and apparatus for allele peak fitting and attribute extraction from dna sample data

Also Published As

Publication number Publication date
JPWO2013008611A1 (en) 2015-02-23
US20140149050A1 (en) 2014-05-29
WO2013008611A1 (en) 2013-01-17

Similar Documents

Publication Publication Date Title
JP5737403B2 (en) Chromatographic data processor
Melamud et al. Metabolomic analysis and visualization engine for LC− MS data
Wei et al. Data preprocessing method for liquid chromatography–mass spectrometry based metabolomics
JP6610678B2 (en) Peak detection method and data processing apparatus
Jiang et al. An automated data analysis pipeline for GC− TOF− MS metabonomics studies
JP6213572B2 (en) Waveform data processing apparatus and waveform data processing program
Vivó-Truyols Bayesian approach for peak detection in two-dimensional chromatography
Tengstrand et al. TracMass 2 A Modular Suite of Tools for Processing Chromatography-Full Scan Mass Spectrometry Data
Zhang Retention time alignment of LC/MS data by a divide-and-conquer algorithm
Jiang et al. Practical and efficient strategy for evaluating oral absolute bioavailability with an intravenous microdose of a stable isotopically-labeled drug using a selected reaction monitoring mass spectrometry assay
WO2015029254A1 (en) Data processing device for chromatograph, and method
Takemura et al. Evaluation of protein-protein docking model structures using all-atom molecular dynamics simulations combined with the solution theory in the energy representation
JP2011122822A (en) Chromatograph mass analyzer
JPWO2017077618A1 (en) Chromatograph mass spectrometry data processing method and processing apparatus
JP6573028B2 (en) Data processing device
JP6036838B2 (en) Comprehensive 2D chromatograph data processor
JP5930066B2 (en) Chromatograph data processing apparatus and data processing method
Mahieu et al. Warpgroup: increased precision of metabolomic data processing by consensus integration bound analysis
JP5211753B2 (en) Chromatographic data processor
Pezzatti et al. Evaluation of different tandem MS acquisition modes to support metabolite annotation in human plasma using ultra high-performance liquid chromatography high-resolution mass spectrometry for untargeted metabolomics
JP3985671B2 (en) Chromatographic data processor
Paritala et al. Quantitative structure retention relationship (QSRR) approach for assessment of chromatographic behavior of antiviral drugs in the development of liquid chromatographic method
JP6493076B2 (en) Peak waveform processing device
Zhou et al. AB-Amy: machine learning aided amyloidogenic risk prediction of therapeutic antibody light chains
Kundrotas et al. Protein–protein alternative binding modes do not overlap

Legal Events

Date Code Title Description
TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20150324

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20150406

R151 Written notification of patent or utility model registration

Ref document number: 5737403

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151