(57)【要約】【課題】例えばピークを構成するデータ点数が少なくピーク波形の曲線形状からは変曲点が求まらない場合でも、高い精度でピーク幅を求め、正しい理論段数を算出する。【解決手段】ピークトップからの前方サーチ・後方サーチにより変曲点が求まらなかった場合、ベースラインＢＬからピークトップまでの高さｈｐの１／ｅ１／２倍（ｅは自然対数の底）の高さにベースラインＢＬと平行な線を引き、その線とクロマトグラムとの２つの交点を仮想的な変曲点Ｃ１、Ｃ２として定める。そして、その２つの変曲点Ｃ１、Ｃ２のそれぞれにおける接線とベースラインＢＬとの交点Ｂ１、Ｂ２を検出し、両交点Ｂ１、Ｂ２間の距離をピーク幅ｗとする。【選択図】  図８
2004184148.tif 000002
【特許請求の範囲】
【請求項１】
ａ）クロマトグラムのデータに基づき、ベースラインを決定するベースライン決定手段と、
ｂ）前記クロマトグラムのデータ基づき、ピークトップを検出するピークトップ決定手段と、
ｃ）前記ベースラインからピークトップまでの高さの１／ｅ１／２倍（ｅは自然対数の底）の高さに前記ベースラインと平行な線を引き、その線とクロマトグラムとの２つの交点を仮想的な変曲点として定める仮想変曲点決定手段と、
ｄ）前記２つの仮想的な変曲点のそれぞれにおける接線又はそれに相当する直線と前記ベースラインとの交点を検出し、両交点間の距離をピーク幅として求めるピーク幅算出手段と、
を備えることを特徴とするクロマトグラフ用データ処理装置。
【請求項２】
前記ピークトップを挟んでピークの前半部と後半部とで、そのクロマトグラムの曲線の形状に基づいて変曲点を検出する変曲点検出手段をさらに備え、該変曲点検出手段による変曲点が求まらない場合に、前記仮想変曲点決定手段により算出される仮想変曲点を採用することを特徴とする請求項１に記載のクロマトグラフ用データ処理装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、ガスクロマトグラフ、液体クロマトグラフ等のクロマトグラフ用のデータ処理装置に関する。
【０００２】
【従来の技術】
一般に、クロマトグラフのカラムの分離性能を表す指標として、理論段数Ｎが用いられる。或るカラムに関して理論段数Ｎを算出することにより、そのカラムの劣化度合等を判断することができる。理論段数Ｎは、クロマトグラムにおける或る成分の保持時間ｔｒとそのピーク幅ｗとから、次の（１）式により求めることができる。
Ｎ＝１６（ｔｒ／ｗ）２      …（１）
【０００３】
米国食品医薬品局（ＦＤＡ）が管轄する米国薬局方（ＵＳＰ）では、ピーク幅ｗを次のように規定している（例えば、非特許文献１参照）。すなわち、図１０に示すように、ピークトップＰを挟んでピークの前半部及び後半部でそれぞれの変曲点Ｃ１、Ｃ２を求め、この変曲点Ｃ１、Ｃ２においてピーク波形に接線を引く。この両接線とベースラインＢＬとの交点Ｂ１、Ｂ２を求め、両交点Ｂ１、Ｂ２間の距離をピーク幅ｗとする。
【０００４】
上記変曲点Ｃ１、Ｃ２を求める一般的なアルゴリズムは、次の通りである。まず、クロマトグラム上でピークトップＰを検出し、前半部の変曲点Ｃ１については、ピークトップＰから前の方向（時間的にマイナス側）に移動しつつ、ピーク波形上の各点において２次微分値を算出する。そして、この２次微分値がゼロとなる点を変曲点Ｃ１として採用する。後半部の変曲点Ｃ２については、ピークトップＰから後の方向（時間的にプラス側）に移動してゆく以外は、同様である。なお、ピークトップＰからではなく、ピークの裾野（ピーク開始点及びピーク終了点）から変曲点をサーチしてゆく方法もある。
【０００５】
【非特許文献１】
“レビュワー・ガイダンス（Ｒｅｖｉｅｗｅｒ  Ｇｕｉｄａｎｃｅ）  バリデイション・オブ・クロマトグラフィック・メソッズ（Ｖａｌｉｄａｔｉｏｎ  ｏｆ  Ｃｈｒｏｍａｔｏｇｒａｐｈｉｃ  Ｍｅｔｈｏｄｓ）”、［Ｏｎｌｉｎｅ］、センター・フォー・ドラッグ・エバリュエーション・アンド・リサーチ（Ｃｅｎｔｅｒ  ｆｏｒ  Ｄｒｕｇ  Ｅｖａｌｕａｔｉｏｎ  ａｎｄ  Ｒｅｓｅａｒｃｈ  （ＣＤＥＲ））、［平成１４年１１月２９日検索］、インターネット〈ＵＲＬ：  ｈｔｔｐ：／／ｗｗｗ．ｆｄａ．ｇｏｖ／ｃｄｅｒ／ｇｕｉｄａｎｃｅ／ｃｍｃ３．ｐｄｆ〉
【０００６】
【発明が解決しようとする課題】
クロマトグラムにおいてピークの幅が小さかったり、或いはピーク幅が大きくてもデータのサンプリング時間間隔が長かったりする場合には、ピークの開始点から終了点までのクロマトグラムデータ点数が相対的に少ない。こうした場合、例えば図１１に示すようにピークの前半部及び後半部が折れ線状になり、上記アルゴリズムに拠っては変曲点が求まらない場合があった。
【０００７】
本発明はかかる課題を解決するために成されたものであり、その目的とするところは、例えばピークの開始点から終了点までのクロマトグラムのデータ点数が少なく、そのピークの立ち上がり・立ち下がりの形状からは変曲点が求まらない場合でも、ピーク幅を正確に求めることができ、それによって正しい理論段数を計算することのできるクロマトグラフ用データ処理装置を提供することにある。
【０００８】
【課題を解決するための手段、及び発明の効果】
クロマトグラムにおけるピークは、理論的には、次の（２）式で表されるガウス分布（正規分布）に従うとされている。
ｆ（ｘ）＝［１／（２π）１／２］・ｅｘｐ［−（ｘ２／２）］      …（２）
この場合、ピークの前半部及び後半部の変曲点の高さは、１／ｅ１／２である。もちろん、実際のクロマトグラムは様々な要因により上記理論通りのピーク形状とはならないから、本来は、ピークトップを挟んでピークの前半部と後半部とで、そのクロマトグラムの曲線の形状に基づいてそれぞれ変曲点を検出することが望ましい。しかしながら、こうした方法で変曲点が適当に求まらない場合、そのピークが上記ガウス分布に従うという仮定の下に仮想的な変曲点を求めても、大きな誤差は生じない筈である。
【０００９】
こうしたことから、上記課題を解決するために成された本発明に係るクロマトグラフ用データ処理装置は、
ａ）クロマトグラムのデータに基づき、ベースラインを決定するベースライン決定手段と、
ｂ）前記クロマトグラムのデータ基づき、ピークトップを検出するピークトップ決定手段と、
ｃ）前記ベースラインからピークトップまでの高さの１／ｅ１／２倍（ｅは自然対数の底）の高さに前記ベースラインと平行な線を引き、その線とクロマトグラムとの２つの交点を仮想的な変曲点として定める仮想変曲点決定手段と、
ｄ）前記２つの仮想的な変曲点のそれぞれにおける接線又はそれに相当する直線と前記ベースラインとの交点を検出し、両交点間の距離をピーク幅として求めるピーク幅算出手段と、
を備えることを特徴としている。
【００１０】
なお、１／ｅ１／２は約０．６０６５３である。本発明の構成によれば、ピークを構成するデータ点数が少なく、ピークトップを挟んだピークの前半部及び後半部で変曲点が適当に求まらないような場合であっても、上述したように、本来の変曲点（ピーク幅に対して充分に小さな時間幅で以てサンプリングできたとした場合のピークの変曲点）とほぼ同じか又はかなり近いような仮想的な変曲点が求まる。従って、ピーク幅算出手段は、例えば米国薬局方（ＵＳＰ）によるピーク幅の算出方法に上記仮想的な変曲点を採用することによって、高い精度でピーク幅を算出することができる。このようにして求めたピーク幅を用い、（１）式のような理論段数計算を行うことにより、正確な理論段数を算出することができる。
【００１１】
また、ピークを構成するデータ点数が充分に得られる場合には、上記のようにピークの前半部及び後半部のピーク曲線の形状に基づいて変曲点を検出することが望ましい。そこで、本発明に係るクロマトグラフ用データ処理装置において、好ましくは、前記ピークトップを挟んでピークの前半部と後半部とで、そのクロマトグラムの曲線の形状に基づいて変曲点を検出する変曲点検出手段を備え、該変曲点検出手段による変曲点が求まらない場合に、前記仮想変曲点決定手段により算出される仮想変曲点を採用する構成とするとよい。
【００１２】
この構成によれば、ピーク幅に対してデータのサンプリング時間間隔が充分に小さくデータ点数が充分にあるような場合には、そのクロマトグラムの曲線形状を反映した正確な変曲点を得ることができ、そうした変曲点が求まらない場合でも、仮想的な変曲点を定めてピーク幅を高い精度で算出することが可能となる。
【００１３】
【実施例】
以下、本発明に係るクロマトグラフ用データ処理装置の一実施例について、図１〜図５を参照して説明する。
【００１４】
図１は、本実施例によるデータ処理装置を備えるガスクロマトグラフ分析システムの概略構成図である。液体試料はシリンジ１１等により試料気化室１２に注入され、そこで気化されて、キャリアガス流路１３から一定流量で供給されるキャリアガス流に乗ってカラム１４内に送り込まれる。試料に含まれる各種成分は、カラム１４を通過する間に時間的に分離されてカラム１４から出て、順次、検出器１５により検出される。検出器１５による検出信号はデジタルデータに変換された後、逐次、データ処理装置１６に送られ、そこで一旦、ハードディスク等の記憶部に格納される。１つの試料の分析が終了した後（或いは、連続に実行される複数の試料の分析が終了した後）、記憶部に格納されたデータが読み出され、クロマトグラムの作成、ピーク検出等の各種のデータ処理が行われる。
【００１５】
なお、データ処理装置１６の実体は専用又は汎用のコンピュータであり、当該コンピュータで所定の処理プログラムを動作させることにより、クロマトグラムの作成のほか、後述するような各種の処理が実行される。
【００１６】
上記データ処理装置１６において理論段数Ｎを算出する際の処理手順を、図２のフローチャートに示す。すなわち、記憶部から読み出したクロマトグラムデータに対し、所定のアルゴリズムに基いてベースラインＢＬ及びピーク（単数又は複数）を検出する（ステップＳ１、Ｓ２）。そして、各ピーク又は必要なピークのみに関し、ピークの開始時間ｔｓ、ピークトップ時間ｔｐ、ピーク終了時間ｔｅ、ピーク高さｈｐ、ピーク面積Ｓ等、ピークを特徴付けるパラメータを算出する（ステップＳ３）。さらに、上記パラメータとクロマトグラムデータとを利用して、ピーク幅ｗを算出し（ステップＳ４）、算出されたピーク幅ｗを用い、（１）式より理論段数Ｎを求める（ステップＳ５）。
【００１７】
本発明に係るデータ処理の特徴はステップＳ４のピーク幅算出処理にあるので、この処理に関し詳細に説明する。図３はピーク幅算出処理のフローチャートである。
【００１８】
まず、図１０に示すように、目的とするピークのピークトップＰよりも時間的にマイナス側のピーク前半部における変曲点Ｃ１を検出する（ステップＳ１０）。この処理の詳細は後に述べる。次に、検出された変曲点Ｃ１におけるピーク曲線の接線を求め、その接線とベースラインＢＬとの交点Ｂ１を検出する（ステップＳ１１）。同様にして、ピークトップＰよりも時間的にプラス側のピーク後半部における変曲点Ｃ２を検出し（ステップＳ１２）、その変曲点Ｃ２における接線とベースラインＢＬとの交点Ｂ２を検出する（ステップＳ１３）。最後に、これら両交点Ｂ１、Ｂ２間の距離を算出し、これをピーク幅ｗとする（ステップＳ１４）。
【００１９】
上記ステップＳ１０のピーク前半部の変曲点Ｃ１の検出処理、及び上記ステップＳ１２のピーク後半部の変曲点Ｃ２の検出処理について、２つの方法が考え得る。そこで、まず第１の方法について、図４及び図５を参照して詳しく説明する。
【００２０】
〔第１の変曲点検出方法〕
図４は、ピーク前半部の変曲点Ｃ１の検出処理を示すフローチャートである。まず初期化処理として、時間変数ｔに上記ピークトップ時間ｔｐをセットする（ステップＳ２０）。次に、ピークトップＰから前方向（時間的にマイナス側）にサーチするべく、時間変数ｔをサンプリング時間の分だけ差し引き、時間的に前方向に移動させる（ステップＳ２１）。そして、その時間ｔ及びその前後の所定個数（例えば、前後各３個の合計７個）のデータより、そこでの２次微分値ｄｄｙ及び１次微分値ｄｙを算出する（ステップＳ２２）。２次微分値ｄｄｙ及び１次微分値ｄｙの算出には、既知のアルゴリズム（例えばサビツキ・ゴーレイ（Ｓａｖｉｔｚｋｙ−Ｇｏｌａｙ）の方法）を用いることができる。
【００２１】
次いで、２次微分値ｄｄｙが０であって且つ１次微分値ｄｙが０より大（つまり正値である）であるか否かを判定し（ステップＳ２３）、その条件が満たされる場合には、時間ｔにおけるクロマトグラム上の点を変曲点Ｃ１と定め（ステップＳ２４）、処理を終了する。但し、採取されたデータは離散的であるため、常にこのような変曲点が現れるとは限らない。そこで、ステップＳ２３の判定に際しては、２次微分値ｄｄｙが負から正に切り替わる場合にも、ｄｄｙ＝０と等価であると看做す処理を行う。
【００２２】
ステップＳ２３で判定条件が満たされない場合には、ステップＳ２５へ進み、現在の時間ｔがピーク開始時間ｔｓ以上である（つまり後方である）か否かを判定する。前方サーチが未だピーク開始時間ｔｓに達していない場合、ステップＳ２１に戻り、時間変数ｔを更に時間的に前方向へと移動させて、ステップＳ２３により同様の判定を行う。その判定条件を満たす変曲点が見つからない限り、上記処理を時間変数ｔがピーク開始時間ｔｓに到達するまで繰り返す。そして、ステップＳ２５で時間変数ｔがピーク開始時間ｔｓに到達したと判定された時点で、前方サーチを終了する。前方サーチの途中で変曲点Ｃ１が見い出されていればステップＳ２３→Ｓ２４→〔終了〕と進んでいる筈であるから、ステップＳ２５でＮｏと判定されるのは、前方サーチによって変曲点Ｃ１が見つからなかったときである。
【００２３】
そこでステップＳ２６へと進み、次のようにして求まる点Ｔｈ１を変曲点Ｃ１として定める。すなわち、図８に示すようにベースラインＢＬとピークトップＰとの間の距離であるピーク高さｈｐのｅ−１／２（約０．６０６５３）倍の位置に、ベースラインＢＬと平行に線を引く。そして、その線とクロマトグラムのピーク前半部とが交差する点を検出し、その点の時間をＴｈ１とする。図８に示すように、ピーク形状をガウス分布とした場合、上記交差点は変曲点にあたる。つまり、理論上の変曲点である。従って、ステップＳ２６の処理により求まる仮想的な変曲点は、仮にそのピークが充分に細かい時間間隔でサンプリングされたとしたときにその結果から求まる本来の変曲点に、非常に近いものとなる。
【００２４】
図５はピーク後半部の変曲点Ｃ２の検出処理を示すフローチャートである。全体の流れとしては図４と同様であり、次の４点において異なるだけである。（１）後方サーチのため、時間変数ｔをサンプリング時間間隔ずつ増加させる（ステップＳ３１）。（２）変曲点の判定条件において、ｄｙ＜０とする（ステップＳ３３）。（３）サーチの終点をピーク終了時間ｔｅとする（ステップＳ３５）。（４）変曲点が求まらなかった場合のピーク後半の変曲点Ｃ２をＴｈ２と定める（ステップＳ３６）。このＴｈ２は、ピーク高さｈｐのｅ−１／２倍の位置にベースラインＢＬと平行に引かれた線とクロマトグラムのピーク後半部とが交差する点の時間である。これにより、ピーク後半部においても、変曲点が見つからなかった場合に、理論に基づく確からしい仮想的な変曲点が求まる。
【００２５】
すなわち、変曲点Ｃ１、Ｃ２を上述のようにして検出することにより、たとえ前方サーチ・後方サーチによっては変曲点が見つからなかった場合でも、理論に則った確度の高い仮想的な変曲点Ｃ１、Ｃ２を採用して、図９に示すようにピーク幅ｗを計算することができる。なお、このときの接線は、変曲点Ｃ１、Ｃ２における１次微分値ｄｙから引くことができるが、それ以外に、Ｃ１、Ｃ２を変曲点とするガウス分布曲線（図８参照）を想定し、そのガウス分布曲線での接線を引くようにしてもよい。
【００２６】
次に、ピーク前半部の変曲点Ｃ１の検出（ステップＳ１０）及びピーク後半部の変曲点Ｃ２の検出（ステップＳ１２）のアルゴリズムが上記第１の方法とは異なる第２の方法について、図６、図７を参照して説明する。この第２の検出方法における基本的な変曲点のサーチは、本発明者が特願２００１−３２０２０２号で提案しているものである。
【００２７】
〔第２の変曲点検出方法〕
ピーク前半部の変曲点Ｃ１を検出する場合、まず初期化処理として、時間変数ｔにピークトップ時間ｔｐをセットするとともに、後述する最近判定基準値ｍｉｎｄｔに計算機が扱い得る最大値をセットする（ステップＳ４０）。次に、ピークトップＰから前方にサーチしつつ変曲点Ｃ１を見つけるため、上記ステップＳ２１〜Ｓ２３と同様のステップＳ４１〜Ｓ４３の処理を実行する。ステップＳ４１で設定された点が変曲点の判定条件を満たさない場合（ステップＳ４３でＮｏ）、上記ステップＳ２５と同様のステップＳ４７に進み、現在の時間ｔがピーク開始時間ｔｓ以上であるか否かを判定する。そして、前方サーチが未だピーク開始時間ｔｓに達していない場合、ステップＳ４１へと戻り、処理を繰り返す。
【００２８】
上記判定条件によりその点が変曲点であると判定されると（ステップＳ４３でＹｅｓ）、その点と所定の点Ｔｈ１との距離ｄｔ（＝｜ｔ−Ｔｈ１｜）を算出し（ステップＳ４４）、その距離ｄｔが最近判定基準値ｍｉｎｄｔよりも小さいか否かを判定する（ステップＳ４５）。初期化処理により最近判定基準値ｍｉｎｄｔには充分大きな値がセットされているため、当初はステップＳ４５の判定は必ずＹｅｓとなる。従って、ステップＳ４６へ進み、距離ｄｔを最近判定基準値ｍｉｎｄｔに代入し、その点ｔを暫定的に変曲点Ｃ１とする。なお、点Ｔｈ１は第１の検出方法で出てきたものと同じ定義である。
【００２９】
上記のような処理を、時間変数ｔがピーク開始時間ｔｓに到達するまで繰り返す。その間、ステップＳ４３において新たな変曲点が検出された場合には、その点と所定点Ｔｈ１との距離ｄｔがそれ以前に検出された変曲点とＴｈ１との距離の最小値ｍｉｎｄｔよりも小さいか否かを判定し（ステップＳ４４、Ｓ４５）、小さい場合にはそれを新たに暫定的な変曲点Ｃ１として置換する（ステップＳ４６）。
【００３０】
ステップＳ４７で時間変数ｔがピーク開始時間ｔｓに到達したと判定された時点で、前方サーチを終了する。その時点で残っている変曲点Ｃ１が最終的に確定された変曲点である。次いで、前方サーチによって変曲点Ｃ１が求まったか否かを判定し（ステップＳ４８）、求まっている場合には処理を終了する。一方、変曲点Ｃ１が求まらなかった場合には、上記ステップＳ２５でＮｏと判定された場合と状況は同じであるから、ステップＳ２６と同様のステップＳ４９の処理を実行し、仮想的な変曲点Ｃ１を求める。
【００３１】
また、図７に示す変曲点Ｃ２の検出処理については、全体の流れとしては図６と同様であり、図４と図５との相違と同じように、次の４点において異なるだけである。すなわち、（１）後方サーチのため、時間変数ｔをサンプリング時間間隔ずつ増加させる（ステップＳ５１）。（２）変曲点の判定条件において、ｄｙ＜０とする（ステップＳ５３）。（３）サーチの終点をピーク終了時間ｔｅとする（ステップＳ５７）。（４）距離ｄｔを算出する際の所定点をピーク後半部における理論上の値であるＴｈ２とするとともに、変曲点が求まらなかった場合のピーク後半部の変曲点Ｃ２をＴｈ２と定める（ステップＳ５４、Ｓ５９）。
【００３２】
なお、この第２の検出方法では第１の検出方法と異なり、ピークトップＰ付近にノイズが存在し、そこに擬似的な変曲点が現れている場合でも、それがピークの変曲点として採用されることはなく、ピーク高さｈｐのｅ−１／２倍に最も近い変曲点がピークの変曲点として正しく採用されるという利点がある。また同様の利点を得るために、所定点Ｔｈ１、Ｔｈ２に最も近い変曲点をピークの変曲点として採用する以外に、例えばピーク高さｈｐの８０〜２０％（或いは７０〜３０％）等、所定の範囲内に入る変曲点をピークの変曲点として採用する等、適宜の変更を行ってもよい。
【００３３】
また、上記実施例は本発明の一例にすぎないから、上記記載の点以外についても、本発明の趣旨の範囲で適宜に変更や修正を加えることができることは明らかである。
【図面の簡単な説明】
【図１】本発明に係るデータ処理装置を備えたガスクロマトグラフ分析システムの概略構成図。
【図２】理論段数算出の処理手順を示すフローチャート。
【図３】ピーク幅算出の処理手順を示すフローチャート。
【図４】第１の方法によるピーク前半の変曲点検出処理のフローチャート。
【図５】第１の方法によるピーク後半の変曲点検出処理のフローチャート。
【図６】第２の方法によるピーク前半の変曲点検出処理のフローチャート。
【図７】第２の方法によるピーク後半の変曲点検出処理のフローチャート。
【図８】前方サーチ・後方サーチで変曲点が求まらない場合の変曲点の算出方法の説明図。
【図９】図８の算出方法によって求まる変曲点を採用したときのピーク幅決定方法の説明図。
【図１０】米国薬局方（ＵＳＰ）によるピーク幅決定法の説明図。
【図１１】前方サーチ・後方サーチで変曲点が求まらない場合のピーク波形の一例を示す図。
【符号の説明】
１１…シリンジ
１２…試料気化室
１３…キャリアガス流路
１４…カラム
１５…検出器
１６…データ処理装置
【図１】
2004184148.tif 000003
【図２】
2004184148.tif 000004
【図３】
2004184148.tif 000005
【図４】
2004184148.tif 000006
【図５】
2004184148.tif 000007
【図６】
2004184148.tif 000008
【図７】
2004184148.tif 000009
【図８】
2004184148.tif 000010
【図９】
2004184148.tif 000011
【図１０】
2004184148.tif 000012
【図１１】
2004184148.tif 000013