JP3583771B2 - Chromatogram analysis method - Google Patents
Chromatogram analysis method Download PDFInfo
- Publication number
- JP3583771B2 JP3583771B2 JP2002297130A JP2002297130A JP3583771B2 JP 3583771 B2 JP3583771 B2 JP 3583771B2 JP 2002297130 A JP2002297130 A JP 2002297130A JP 2002297130 A JP2002297130 A JP 2002297130A JP 3583771 B2 JP3583771 B2 JP 3583771B2
- Authority
- JP
- Japan
- Prior art keywords
- chromatogram
- peak
- analysis method
- function
- waveform
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Treatment Of Liquids With Adsorbents In General (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Description
【0001】
【産業上の利用分野】
本発明は、液体クロマトグラフィ、ガスクロマトグラフィ等のクロマトグラフィ技術に係り、特にクロマトグラム解析方法(以下、解析方法という)およびそれに使用されるクロマトグラフ装置、前記解析方法およびクロマトグラフ装置に用いられるデ−タ処理装置に関するものである。
【0002】
【従来の技術】
従来の技術の説明に先立ち、技術用語の定義をする。
本出願においては、コンボリューションとは、二つの関数x(t)、y(t)に対するたたみ込み積分として定義され、一般に計測機器からの出力波形は、入力波形と装置関数(電気回路ではインパルス応答、分光器ではスリット関数)とのコンボリューションとなる。
ディコンボリューションとは、コンボリューションの逆演算として定義される。具体的にはコンボリューションであらわされている積分方程式を解くことに相当する。
リコンボリューションとは、装置関数、すなわち拡がりを表わす関数h(t)を用いてディコンボリューションした波形に対し、ふたたびその関数h(t)を用いてコンボリューションすることをいう。
【0003】
さて、従来から一般に計測機器からの出力波形デ−タは、その計測機器の動特性、すなわち装置関数の影響によって歪をうける。
クロマトグラフ装置においても同様であり、装置関数、すなわち拡がりを表わす関数h(t)によりクロマトグラム上に歪を生じたり、クロマトグラム上にピ−クの重なりが生じる。
このような場合、それぞれクロマトグラムのピ−ク面積を正確に求めることが困難であった。
【0004】
従来、上記歪をうけた場合は、ディコンボリュ−ション処理により、上記ピ−クの重なりに対しては、通常曲線適合法を基本とした合成分離法を用い、デ−タ処理的に分解し、定量することが提案されている。
これに関する一般的な解決手法については、南茂夫編著「科学計測のための波形デ−タ処理」、CQ出版株式会社刊、(1986)の記載技術がある。
また、多くの改善技術が、例えば特開昭62−17465号公報記載の技術、特開昭63−151851号公報記載の技術等が提案されている。
【0005】
前記特開昭62−17465号公報記載の技術は、立上り側の半値幅が、立上り地点から終了点までの距離の所定倍になるガウシアンにて近似するものであるが、クロマトグラムが重なりピ−ク波形の場合、必ずしもこのような近似で置換しうる場合は少ない。
この技術は、記憶容量を減らすことには有効であるけれども、ピーク波形が重なった場合、ピーク面積を正確に求めることについては、必ずしも有効ではないという欠点があった。
【0006】
前記特開昭63−151851号公報記載の技術においては、従来、重なりピークの分割をいわゆる垂直分割、スキミング分割処理をしていたものを二つのガウシアンに分割して近似し、最小二乗法と呼ばれる手法にてディコンボリュ−ション処理を行っているが、物理的に意味のある解が得られないという欠点があった。
【0007】
また、米国特許4807148号、米国特許4941101号明細書等記載の技術においては、ディコンボリュ−ション処理において数値解析の手法が用いられ、Janssonの方法、ヤコビ法、ガウス−ザイデル法等があり、演算時間、記憶語数の点で有利となっていたが、解の不連続点の拡大等に欠点があった。
【0008】
【特許文献1】
特開昭62−17465号公報
【特許文献2】
特開昭63−151851号公報
【特許文献3】
米国特許4807148号明細書
【特許文献4】
米国特許4941101号明細書
【非特許文献1】
南茂夫編著「科学計測のための波形デ−タ処理」、CQ出版株式会社刊、(1986)
【0009】
【発明が解決しようとする課題】
上記従来の種々のディコンボリュ−ション処理は、クロマトグラムが歪をうけた場合のディコンボリュ−ション処理、クロマトグラムの重なりピ−ク波形を生じた場合に対する合成分離については、それぞれ個々に解決法を提供するが、歪をうけることによる重なりピ−ク波形を生じた場合には、歪をうけたことを考慮し、重なり波形を分離するところまでは詳しく検討され、対処されていなかった。
【0010】
また、上記従来におけるクロマトグラムピ−ク波形の重なりに対する合成分離法は、次のような問題がある。
1.ガウシアン等の非線形関数を仮定した場合、その関数を選択する妥当性が疑問である。
2.UVスペクトル等のマルチチャンネル情報を導入する場合、フォトダイオ−ドアレイ検出器等の高価な検出器が必要となる。
3.ピ−ク分解においては最小自乗法を用いるため、解が得られるまでに長い計算時間を要し、また物理的に意味のある解が得られない場合もある。
【0011】
本発明は、上記従来技術の問題点を解決するためになされたもので、クロマトグラムの波形の歪を補正し、重なりのピ−ク波形をシャ−プにすることにより分離を向上させ、ピ−ク面積をより正確に求め、定量分析の精度が改良されたクロマトグラム解析方法を提供することを目的とする。
【0012】
【課題を解決するための手段】
上記目的を達成するために、クロマトグラム解析方法に係る本発明の構成は、分離カラムを通過した分析対象試料を検出器によって検出し、検出によって得られたクロマトグラムの解析を行うクロマトグラム解析方法において、
前記分析対象試料を前記分離カラムへ導入と別に、前記分離カラムによって保持されない標準試料を前記分離カラムへ導入するステップと、
前記クロマトグラム上から、前記標準試料のピークを同定するステップと、
前記クロマトグラムから前記分析対象試料の任意のピークのある時間の始点と終点を設定するステップと、
前記設定された時間区間内のデータを、前記同定された標準試料のピーク波形を用いてディコンボリューション処理を行うステップと、を有することを特徴とする。
前記に記載のクロマトグラム解析方法において、
前記ディコンボリューション処理を行った後のピークが複数のピークで、且つ重なりあった状態のままであった場合に、当該重なりピークに対して、垂直分割、あるいは非線形関数を用いた面積分割を行うことを特徴とする。
【0013】
上記目的を達成するために、クロマトグラム解析方法に係る本発明の他の構成は、分離カラムを通過した分析対象試料を検出器によって検出し、検出によって得られたクロマトグラム内のピークを標準偏差σを係数として有するガウシアン関数を用いて表し解析を行うクロマトグラム解析方法において、
前記ガウシアン関数の標準偏差σに任意の値を設定するステップと、
前記クロマトグラムに対してディコンボリューション処理を行うステップと、
当該ディコンボリューション処理後のクロマトグラムについて、クロマトグラムが発散していないか、負の下限値を下回っていないか、不分離ピークがピーク間の谷を形成しているか、のいずれかの基準を満たしているかどうかの判定を行うステップと、
当該判定で不合格の場合、前記ガウシアン関数の標準偏差値を調整し、再度ディコンボリューション処理を行うステップと、を有することを特徴とする。
前記に記載のクロマトグラム解析方法において、
前記ディコンボリューション処理に使用した拡がり関数を用いて、前記ディコンボリューション処理によって分離されたピークに対してリコンボリューション処理後の波形を前記ディコンボリューション処理前のピーク波形と重ねて表示することを特徴とする。
【0014】
より詳しく説明すると、クロマトグラムのピ−クは、種々の原因により拡がっている。例えば配管や検出器のフロ−セル等の流路を試料が流れることにより、試料ゾ−ンが拡がる。これはカラム外の拡がりと呼ばれている。
このカラム外の拡がりに代表される拡がりによる寄与は、オリジナルのピ−ク波形にコンボリュ−ション演算され、検出波形に現われる。
つまり、拡がり関数が既知であれば、検出波形をディコンボリュ−ション演算することによりオリジナルのピ−ク波形を計算することができる。
【0015】
拡がり関数h(t)を決定する手法はいろいろなものが考えられる。
手法を大きく2つに分類すると、実験的に検出デ−タから求める方法と波形を特徴づけるいくつかのパラメ−タを持つ解析数学的なh(t)を仮定する方法とがある。
前者はh(t)が既知の方法、後者はh(t)が未知の方法と考えることができる。
【0016】
前記h(t)を測定から求める方法として、例えばカラム外の拡がりを得るために、カラムをジョイントと交換し、試料を注入し、検出器で測定された波形をh(t)とみなす方法がある。
試料内の成分の種類に依存しないようなカラム内の拡がりも付加するなら、保持しない試料を注入し、そのピ−ク波形をh(t)とみなすことができる。
【0017】
保持しないピ−クを得る方法には、他に溶離液に溶出能力の極めて強い溶媒を使用すること、カラム充填剤に官能基を修飾しないものを作成し保持しないようにする等が考えられる。
また、標準試料の中に保持しない成分を添加し、非保持ピ−クを標準試料のクロマトグラムから得る方法も有効である。
【0018】
h(t)が未知の場合、次のような点を前提条件し、ピ−ク幅の一部分に相当する拡がり関数h(t)をディコンボリュ−ション処理し、よりシャ−プなピ−クを持つ処理クロマトグラムからピ−ク面積をもとに正確に定量することができる。
【0019】
その前提条件とは、検出されるピ−ク波形はオリジナルにはデルタ関数的あるいは非常にシャ−プな長方形的な濃度分布関数に対して、拡がり関数をコンボリュ−ション演算したものであると考えられる。
ここで拡がり関数は、ほぼガウシアンであり、それが若干テ−リング的あるいはリ−ディング的に修正された波形をしている。
つまり、拡がり関数自体もガウシアンに対してテ−リングかリ−ディングの修正関数がコンボリュ−ト演算されていると考えられる。
【0020】
また、ガウシアンには2つのガウシアンをコンボリュ−ション演算すると、それもまたガウシアンであるという加成性がある。
以上により拡がり関数の一部分を構成しているガウシアンをディコンボリュ−ション演算することにより、よりシャ−プなオリジナルな波形に近づいたピ−クを計算することができることになる。
【0021】
【作用】
上記各技術的手段の働きは次のとおりである。
本発明の構成によれば、検出クロマトグラムから既知のh(t)またはガウシアンを合理的に選択し、ディコンボリュ−ション演算することにより、検出クロマトグラムの重なりピ−クよりシャ−プなピ−クを持つクロマトグラムに分離・独立させ、分離・独立した各ピ−クを積分した時、正確にピ−ク面積が求められ、正確な定量をすることができる。
上記ディコンボリュ−ション演算は、例えば非線形関数の最小自乗法と異なり、短い計算時間により物理的に意味のある解を必ず得ることができる。
【0022】
また、上記ディコンボリュ−ション演算は、重なりピ−クが分離・独立されたときに、そのクロマトグラムを観察し、その演算の妥当性を判断しているが、熟練が必要となる。そこで、ディコンボリュ−ション演算に用いた拡がり関数h(t)を用い、独立した各ピークにリコンボリュ−ション演算を施し、最初の重なりピ−クを構成している各ピークを再現することができる。このリコンボリュ−ション演算後の各ピークと最初の重なりピ−クとを重ね描きすることにより上記ディコンボリュ−ション演算の妥当性を正確に判断することができる。
【0023】
また、本発明の構成によれば、例えばフォトダイオ−ドアレイ検出器等の高価な検出器を省略し、簡単な構成によりクロマトグラム解析方法を実施することができる。
【0024】
【実施例】
〔実施例 1〕
以下、本発明に基づく一実施例を図1ないし図11を参照して説明する。
図1は、本発明における一実施例に係るクロマトグラム解析方法(以下、解析方法という)に用いるマイクロボアHPLCシステムの構成を示すブロック図、図2は、図1の解析方法におけるクロマトグラムを示す線図、図3は、図1の解析方法におけるクロマトグラムの拡がりメカニズムの略示説明図、図4は、図1の解析方法におけるディコンボリュ−ション演算説明図、図5は、図1の解析方法に係る処理工程のフロ−チヤ−ト、図6は、図1の解析方法におけるクロマトグラムの重ね描きプロット図である。
【0025】
また、図7は、図1の解析方法によるアミノ酸分析のクロマトグラムを示す線図、図8は、図1の解析方法によるグリコヘモグロビン分析のクロマトグラムを示す線図、図9は、図1の解析方法において拡がり関数が未知の場合の処理フロ−チャ−ト、図10は、図1の解析方法におけるカテコ−ルアミン分析のクロマトグラムを示す線図、図11は、図1の解析方法におけるリコンボリュ−ション演算による重ね描き説明図である。
【0026】
本発明に係るクロマトグラム解析方法を実施するため、マイクロボアHPLCシステムを使用する。図1は、マイクロボアHPLCシステムの全体構成を示している。
図1において、40はポンプ、41はカラム、42は検出器、43はサンプラ、44は制御部、45はデ−タ解析部、46はCRT表示部、47はプリンタ、49は溶離液、51は未知試料、52は標準試料、54は入力部である。
【0027】
本システムを構成する各部は、デ−タ解析部45(〔実施例 2〕において後述)を除き、通常のクロマトグラム解析方法を用いられるものである。したがって、詳細な説明は煩瑣となるので省略し、簡単に説明する。
【0028】
前記マイクロボアHPLCシステムでは、カラム外の拡がりが分離に大きく影響し、無視できない。制御部44の命令によりポンプ40は溶離液49を送液する。可動ニ−ドル弁を有するサンプラ43は、未知試料51あるいは標準試料52をカラム41への流路に注入する。
試料は、溶離液49とともにカラム41に送りこまれ、含有成分が分離展開され、検出器42で検出される。
【0029】
前記検出されたクロマトグラムは、デ−タ解析部45に記憶される。
図2は、検出されたクロマトグラムである。前記標準試料52には非保持成分を添加しているため、非保持ピ−ク61としてそれが検出される。次に溶出されるピ−ク62は標準試料の一成分であり、ピ−ク61と異なりカラム内に保持され、ピ−ク幅も若干広くなっている。
【0030】
図3を参照し、ピ−ク波形が形成されるメカニズムについて説明する。
図3は、図1の一実施例に係る解析方法におけるクロマトグラムの拡がりのメカニズムの説明図である。
図3において、図中、図1の符号と同一符号は同等部分であるので詳細な説明は省略する。
上段(a)図は、試料流路を示すブロック図であり、下段(b)図は、上段(a)図の試料流路における試料の濃度を示している。
【0031】
単一成分を含む試料は、サンプラ43により流路内に注入される。この時の流路内での試料濃度Cの空間分布図は、流れの方向をxとすれば、長方形状の分布線図71である。
この長方形状の分布線図71がカラム41を通過すると拡大する。前記カラム41を通過した段階は、分布線図72のようなガウシアンの形状である。
【0032】
さらに、検出器42において、前記ガウシアンの濃度分布を検出する時には流路配管や検出器のフローセル等により濃度分布が拡げられ、分布線図73のように分布線図72よりも広い分布形状になっている。
【0033】
これを数式で記述すると、(2)式のようなモデルで表現できる。
【数2】
上式において、時刻tにおける検出波形をy(t),サンプル注入時のオリジナルな長方形波形をr(t),拡がり関数をh(t)とする。上式においては、検出波形y(t)は、オリジナルな長方形波形r(t)と、拡がり関数h(t)とのコンボリュ−ションであることを示している。
【0034】
ピ−クの形状の拡がりについては、上記説明でよいが、リテンションタイムについても考慮すると、もう一つDiracのデルタ関数がコンボリュ−ション演算されていると考えられる。これを(3)式に示す。
【数3】
【0035】
上式において、EMG(t)は、Exponentially Modified Gaussinであり、σ、τ、tRは、EMG(t)のパラメータ、δ(t)はデルタ関数、m(t)は修正関数である。
拡がりのメカニズムについては、“DYNAMICS OF CHROMATOGRAPY”,J.C.Giddings著,Marcel Dekker刊,New York,1965に詳しく説明されているのでそれに譲る。
【0036】
拡がり関数h(t)は、カラム内での拡がりとカラム外での拡がりの両方からの影響により形成される。つまりカラム内での拡がり関数とカラム外での拡がり関数のコンボリュ−ションである。
【0037】
図2における非保持ピ−ク61に着目し、図3を参照し説明すると、注入時は長方形波形71だったものが、カラム41を通過し、ガウシアン波形72に拡げられ、カラム外の拡がりの影響を受け、波形73となり、ピ−ク61として検出されていることになる。
【0038】
一方、保持ピ−ク62も同様なメカニズムで拡げられ、ピ−クが形成されているのであるが、前記カラム41内に保持されているために、若干広めのピ−ク幅に拡げられている。
すなわち、ピ−ク62に影響する拡がり関数は、ピ−ク61と共通の拡がり関数とそれ以外の余剰の拡がり関数とのコンボリュ−ションになっている。
【0039】
この場合、ピ−ク62の拡がり関数からピ−ク61の拡がり関数の影響を取り除き、余剰の拡がり関数だけを得ることができる。この処理をディコンボリュ−ション処理と呼ぶ。
【0040】
前記ディコンボリュ−ション演算は、図2に示すピ−ク64とピ−ク65とのような重なりピ−クに対して効果がある。
図4を参照して説明する。図4は、図2に示す重なりピ−クに対しディコンボリュ−ション演算を説明している。
【0041】
前記ピ−ク64波形と前記ピ−ク65波形とに対して前記ピ−ク61波形の拡がり関数をディコンボリュ−ション演算すると、図4に示すような分離された余剰分のピ−クが得られ、より正確にそれぞれのピ−クの面積が計算できる。
実際は、ピ−ク波形61を面積1に規格化した規格化関数をディコンボリュ−ション演算に利用する。規格化はピ−ク面積を保存するために行われる。
【0042】
厳密に言えば、この規格化関数は、(式)2の注入時の長方形関数r(t)とピ−ク波形61に影響する拡がり関数h(t)の規格化されたコンボリュ−ションである。
実際には、ピ−ク波形61以外のピ−クにも長方形関数の影響があるため、ピ−ク61波形の規格化関数でディコンボリュ−ション演算しても問題にはならないわけである。
【0043】
以上の図4のディコンボリュ−ション手順を図1のマイクロボアHPLCシステムのブロック図と、図5に示す検出波形処理のフロ−チャ−トとを参照して説明する。
まず、非保持成分を添加した標準試料52が注入され、工程31にて波形処理が開始される。
【0044】
次に、工程32で非保持ピ−ク61を同定する。これはクロマトグラムの最初の時刻に現われるある大きさ以上のピ−クとして自動的に同定するか、またはオペレ−タが入力部54から時刻指定を行うか、CRT46に表示しているクロマトグラムからマウスを使用し、ピ−ク61をピックアップすることで指定することができる。
【0045】
工程33において、非保持ピ−ク61の面積を1に規格化する。この時、CRT46に規格化関数を表示することもできる。
工程34において、オペレ−タが正確に定量したい重なりピ−クのある時間の始点と終点を設定する。前記CRT46には、図4の左図のように、この部分だけを拡大したクロマトグラムを表示することもできる。ここでオペレ−タはクロマトグラムの全ての時間を設定することも可能である。
【0046】
工程35において規格化関数を用いて、設定された時間区間内でディコンボリュ−ション処理を行う。
工程36において処理前後のクロマトグラムをCRT46に図4の右側と左側とのように表示する。また、図6に示す如く、プリンタ47に重ね描きしたクロマトグラムをプロットすることもできる。
【0047】
工程37において、シャ−プになった図4の右側のピ−クを対象に積分を実行し、定量計算を行う。
この処理後のクロマトグラムでもまだピ−クが重なっている場合には、従来通りの垂直分割等で面積計算する。
ここでガウシアン、EMG等の非線形関数で最小自乗法により面積分割することもできる。ここでCRT46かプリンタ47に定量値を出力する。
工程38において処理を終了する。
【0048】
ディコンボリュ−ション処理をする拡がり関数を標準試料52中の非保持成分のピ−ク波形の規格化関数で代用したが、この方法の他にも拡がり関数を得ることができる。
標準試料に添加する代わりに、予め拡がり関数を求めるために非保持成分だけ注入し、ピ−ク波形を得ることができる。
【0049】
また、逆相クロマトグラフィ−ならメタノ−ル,アセトニトリル,THF,クロロホルム等有機溶媒100%の溶離液を流し、保持の弱い成分を注入し、非保持ピ−クを得ることができる。
イオン交換クロマトグラフィ−では、例えばアミノ酸分析計の場合、pHの高い再生液を流し、保持の弱いホスホセリン等を注入し、そのピ−ク波形を得ることができる。
【0050】
また、カラム充填剤を操作し、非保持ピ−クを得ることもできる。イオン交換クロマトグラフィ−なら、イオン交換基の修飾しない坦体だけの充填剤を作成し、非保持ピ−クを得ることもできる。
【0051】
ポストカラム反応LCの実施例としてアミノ酸分析計をあげる。このシステムもカラム外の拡がりを無視できない。カラム外の拡がりは、カラムをジョイントに交換し、非保持ピ−クの波形を拡がり関数として見積ることができる。
【0052】
図7に示す如く検出されたクロマトグラムにおいては、スレオニン(Thr)とセリン(Ser)の重なり部分が14%であるのに対し、ディコンボリュ−ション処理されたクロマトグラムでは4%と1/3以下に向上した。
この処理クロマトグラムは、ポストカラム反応LCにおいて、カラムの直後で検出できたとした場合のクロマトグラムに相当している。
【0053】
拡がり関数の共通部分の除去例として、3.3分間処理のグリコヘモグロビン分析計をあげる。
拡がり関数を図8の左側に示すHbFのピ−ク波形から見積ることができる。これはHbA1cのピ−ク幅が、リテンションタイムの順序を考慮し、HbFの幅よりも若干広いと考えられるためである。
【0054】
図8の左側の検出クロマトグラムより孤立したピ−クであるHbFの波形を規格化し、拡がり関数を得る。その規格化関数を用いて、不安定型と安定型とのHbA1cの重なりピ−クをディコンボリュ−ション処理すれば、図8の右側に示すような分離の向上したクロマトグラムが求められるわけである。
【0055】
ディコンボリュ−ション処理は、この場合処理クロマトグラムを直接表示する必要はないが、定量精度を向上するためには、一般的なHPLCであるグリコヘモグロビン分析計においても非常に有効な手法である。
【0056】
次にディコンボリュ−ション演算を正しく行うために演算前に利用すべき処理について説明する。図2に示されているノイズ67やドリフト68がディコンボリュ−ション演算を妨害する。
これらをそのままディコンボリュ−ション処理すると、擬似的なピ−クとして増幅されたり、ディコンボリュ−ション処理されるべきピ−クが正しく処理されないといった弊害が起こる。
【0057】
まず、ノイズ除去処理としては一般的に検出クロマトグラムのデ−タ点をいくつか足し算して低減する方法が採られる。
この他に平滑化の方法や公知の高周波成分を除去する方法も利用できる。
【0058】
ドリフト除去処理としては一般的にピ−ク間の谷から谷を直線で結だり、強制的に水平に直線を描き、検出クロマトグラムからベ−スラインを引き算する方法が採られる。
この他に公知の低周波成分を除去する方法,微積分を利用し低次の多項式成分を除去する方法、微分の特徴点からベ−スラインを探索する方法,ニュ−ラルネットワ−クを利用する方法等がある。
【0059】
ここまでは拡がり関数を測定デ−タより求めてきたが、拡がり関数が未知である場合でも、ガウシアンをディコンボリュ−ト演算することができる。
前述したように検出クロマトグラムのピ−ク波形は、(2)式で示されるように、拡がり関数h(t)と長方形波形のr(t)とのコンボリュ−ションで表わせる。
【0060】
この拡がり関数h(t)は、前述の“DYNAMICS OF CHROMATOGRAPY”J.C.Giddings著により、カラム内の拡がりがどのような物理学的,化学的相互作用に起因していても、カラム外の拡がりが無視できないくらい大きくとも、ガウシアン様の形状をしている。
【0061】
h(t)がリ−ディンク的かまたはテ−リング的にガウシアンから多少歪んでいたとしても、下記の(4)式のようにガウシアンG(t)と修正関数m(t)とのコンボリュ−ションとして記述できる。
【数4】
【0062】
上式において、ガウシアンをG(t)、修正関数をm(t)、σは標準偏差である。
例として、h(t)がEMG(t)により表わせるテ−リングピ−クであるなら、前述の(3)式のように解析数学的に修正関数が求められる。
【0063】
また、重なりピ−クは、(5)式のように表わせる。
【数5】
ここで、σ1、σ2はそれぞれ標準偏差である。
【0064】
ガウシアンは、上記の如く、加成性があるため、2つのガウシアンのコンボリュ−ションに分解できるため、下記(6)式のような関係を持っている。
【数6】
【0065】
つまり、(5)式で表わされている2つのピ−ク波形はガウシアンを用いて、ディコンボリュ−ション演算ができる。
この時そのガウシアンは、2つのピ−クの持つ標準偏差σ1とσ2よりも小さな標準偏差σ0を持つ幅の狭い波形でなければならない。
【0066】
こうして、検出された重なりピークよりも幅の狭いピークを得ることができる。
最終的には、y(t)をガウシアンG(t,σ0)によりディコンボリューション演算すると、クロマトグラムは、幅のより狭いガウシアンともとのままの修正関数とのコンボリューションにより表わすことができるわけである。
【0067】
【数7】
上式の〔 〕内のように、σ1,σ2よりピーク幅の狭いピークが得られる。より詳しくいえば、注入時の長方形波形はこの修正関数によりコンボリュ−ション演算されている。
【0068】
この拡がり関数h(t)が未知の場合の処理を図9を参照して説明する。
工程71で処理を開始する。
工程72で検出クロマトグラムからノイズ・ドリフトを除去し、ディコンボリュ−ション処理の準備をする。
工程73でディコンボリュ−ション処理を手動で行なうか、自動で行なうのかを入力する。
【0069】
自動の場合、工程81で、検出クロマトグラムのピ−ク幅を全体的に見て、最小の幅よりもやや小さな標準偏差を持つガウシアンを拡がり関数として設定する。これは面積1に規格化されている。
ここで、リテンションタイムの短いピ−クと長いピ−クでピ−ク幅が2倍以上異なる場合は、時間を分割してディコンボリュ−ション演算することもできる。
【0070】
ここで自動的にピ−ク幅を見つけ出すためにフ−リエ変換を利用することもできる。工程76でディコンボリュ−ション処理を実行する。工程77で処理クロマトグラムが発散していないか、負の下限値を下回っていないか判定する。また不分離ピ−クがピ−ク間の谷を形成しているか等を判定する。
【0071】
ここで判定が合格となれば、工程79へ進む。判定が不合格の場合、もう一度工程81へ戻り、判定を参照し標準偏差を加減し、ディコンボリュ−ション処理を行なう。
判定合格の場合、工程79で処理クロマトグラムの表示,定量結果の表示等を行なう。工程80で終了する。
【0072】
手動の場合、工程74で入力操作を行う。一般的にはディコンボリュ−ション処理するガウシアンの標準偏差を入力する。他に処理区間の始・終点の入力もできる。また特定の重なりピ−クの分離を向上するために、ガウシアンの最適な標準偏差を探索することもできる。
【0073】
この場合、標準偏差の初期値,推定されるピ−クの個数,大きさ,幅,保持時間等を参考値として入力することもできる。また制限値として、ピ−ク個数は3個までとか、重なりピ−クの幅が0.2分以上0.4分以下のように入力することもできる。
【0074】
工程75で工程74にて入力された情報を演算し、ディコンボリュ−ション処理の設定値を求め、判定の基準値を算出する。工程76でディコンボリュ−ション処理を実行する。工程77で判定を行う。
不合格の場合、その理由を工程78で表示し、工程74で再入力を促す。合格の場合、工程79で処理クロマトグラムと定量結果を表示する。
工程80で終了する。
【0075】
カテコ−ルアミン分析計により得られた生の10分間のクロマトグラムを図10の上段(a)に示す。
これは、血漿や尿のような生体試料の分析例である。生体試料を分析すると数多くのピ−クがしばしば現われる。これを適切なガウシアンでディコンボリュ−ション処理すると、中段(b)のようなクロマトグラムが得られる。
各ピ−クがベ−スラインまで完全に分離している。
【0076】
しかし、これをそのまま表示しても、処理クロマトグラムの意味するところは理解困難である。
むしろ、これを下段の(c)に示した公知の表示法にしても差し支えない。
また、定量値のみを出力するだけで、ディコンボリュ−ション処理をオペレ−タ−に意識させないのがよいと考えられる。
【0077】
以上、本発明の一実施例について説明したが、本発明は、上記実施例に限定されるものでなく、種々の変形例が考えられる。
本発明に係るクロマトグラム解析方法の他の一方法を説明する。
図11は、本発明に係る解析方法におけるディコンボリュ−ション、リィコンボリュ−ションによる重ね描きの説明したものである。
【0078】
図11(1)は、検出されたオリジナルの波形である。このオリジナルの波形にディコンボリュ−ション処理し、各波形ピークを分離・独立させる。この状態を図11(2)に示す。
分離・独立させた波形ピークをそれぞれディコンボリュ−ションしたときと同じ拡がり関数h(t)によりリコンボリュ−ションする。
リコンボリュ−ション演算された波形を検出されたオリジナルの波形に重ね描きをする。図11(3)は、この状態を示したものである。
【0079】
〔実施例 2〕
本発明に係るクロマトグラフ解析装置の他の一実施例を説明する。図12は本発明に係る他の一実施例であるクロマトグラフ解析装置の構成を示すブロック図である。
本発明に係るクロマトグラフ解析装置の全体構成については、〔実施例 1〕において、図1のマイクロボアHPLCシステムについて説明したので、煩瑣となるので再度の説明を省略し、データ解析部45のみを説明する。
【0080】
前記データ解析部45は、ノイズ・ドリフト除去機構と、拡がり関数解析機構と、ディコンボリューション演算機構と、定量計算機構と、判定機構とから構成されている。図中、図示されてないが、A/D変換器、記憶装置、インタフェース部等も備えていることはいうまでもない。
前記構成の内、ノイズ・ドリフト除去機構、定量計算機構、判定機構は、公知の技術であるが順序としてあわせて説明をする。なお、前記各機構部は、バスによって接続されている。
【0081】
前記ノイズ・ドリフト除去機構は、平滑化機構と、ベースライン算出機構と、高速フーリエ変換機構と、ニューラル・ネットワーク機構と、微分・積分機構と、ファジイ推論機構とを備えており、前記平滑化機構は、検出器43からのクロマトグラムの測定データ群を移動平均法、Savitzky−Golay法、Kawata−Minami法、周波数領域法のいずれかにより滑らかな曲線にフイットし、前記ベースライン算出機構は、適切なアルゴリズムによりベースラインを決定する。
【0082】
前記高速フーリエ変換機構は、高周波成分であるノイズおよび低周波成分であるベースライン変動をそれぞれ周波数分布に基づき除去する。
また、前記ニューラル・ネットワーク機構は、ドリフト形状を予め学習しておき、測定データから所定の結合定数に従い、ドリフトを除去する。
【0083】
前記微分・積分機構は、微分・積分により低次の多項式で表現されるベースライン成分を除去する。この整理された測定データ点群をファジイ推論機構が所定の結合定数に従い平滑化する。このようにして検出器43からのクロマトグラムのノイズおよびドリフトは除去される。
【0084】
前記拡がり関数解析機構は、Gaussian発生機構と、面積1規格化機構と、検出ピーク波形設定機構とを備え、解析数学的拡がり関数の場合は、前記Gaussian発生機構によりGaussianを発生する。実測的拡がり関数の場合は、前記面積1規格化機構により実測されたピーク波形を規格化し、前記検出ピーク波形設定機構により前記規格化された実測波形を拡がり関数としている。このようにして実測的ないし解析数学的拡がり関数を発生する。
【0085】
前記ディコンボリューション演算機構は、時間区間設定機構と、拡がり関数設定機構と、ディコンボリューション処理機構とを備え、前記時間区間設定機構は、ディコンボリューション処理する区間を設定し、前記拡がり関数設定機構は、上記の発生させた実測的ないし解析数学的拡がり関数から拡がり関数を設定する。
前記ディコンボリューション処理機構は、設定されたディコンボリューション処理する区間の重なり波形を前記拡がり関数によりディコンボリューション処理し、重なり波形を分離・独立させる。
【0086】
前記定量計算機構は、積分機構を備え、前記積分機構は、前記分離・独立させたピーク波形の面積をそれぞれ算出する。このようにしてピーク面積ないし高さに基づき定量を計算する。前記定量結果は、CRT46、プリンタ部47に出力される。
【0087】
各種の判定を行う判定機構は、ディコンボリューション処理の成功を判定するために利用する判別式の値を計算する判別式計算機構、ディコンボリューション処理するために利用するGaussinを探索するときに制限となる値を算出する制限値算出機構と、所定の結合係数に従いその最適なGaussinを探索する拡がり関数探索機構を具備し、上記各処理段階において、各機構が所定の働きが可能となるようサポートしている。
【0088】
図12における制御部44、データ解析部45、CRT部46、プリンタ部47等の機能を一体化して、A/D変換部と、メモリ部と、制御部と、インタフェース部と、出力機器とからなり、クマトグラムの検出ピ−ク波形から拡がりを表わす関数h(t)を算出・出力し、検出ピ−ク波形および前記関数h(t)を用いて前記検出ピ−ク波形をディコンボリューション処理し、ディコンボリューション処理後のクロマトグラムを出力することを特徴とするクロマトグラム解析方法とクロマトグラフ装置とに用いられるデ−タ処理装置を構成しても差し支ええない。
【0089】
【発明の効果】
以上詳細に説明したように、本発明の構成によれば、クロマトグラムの波形の歪を補正し、重なりのピ−ク波形をシャ−プにすることにより分離を向上させ、ピ−ク面積をより正確に求められ、定量分析の精度が改良されたクロマトグラム解析方法を提供することができる。
【図面の簡単な説明】
【図1】本発明における一実施例に係る解析方法に用いるマイクロボアHPLCシステムの構成を示すブロック図である。
【図2】図1の解析方法におけるクロマトグラムを示す線図である。
【図3】図1の解析方法におけるクロマトグラムの拡がりメカニズムの略示説明図である。
【図4】図1の解析方法におけるクロマトグラムの重なりピークに対するディコンボリュ−ション演算説明図である。
【図5】図1の解析方法に係る処理工程のフロ−チヤ−トである。
【図6】図1の解析方法におけるクロマトグラムの重ね描きプロット図である。
【図7】図1の解析方法によるアミノ酸分析のクロマトグラムを示す線図である。
【図8】図1の解析方法によるグリコヘモグロビン分析のクロマトグラムを示す線図である。
【図9】図1の解析方法において拡がり関数が未知の場合の処理フロ−チャ−トである。
【図10】図1の解析方法におけるカテコ−ルアミン分析のクロマトグラムを示す線図である。
【図11】図1の解析方法におけるリィコンボリューション演算による重ね描き説明図である。
【図12】本発明に係る他の一実施例に係るクロマトグラフ解析装置の構成を示すブロック図である。
【符号の説明】
40…ポンプ
41…カラム
42…検出器
43…サンプラ
44…制御部
45…デ−タ解析部
46…CRT表示部
47…プリンタ
49…溶離液
51…未知試料
52…標準試料
54…入力部[0001]
[Industrial applications]
The present invention relates to chromatography techniques such as liquid chromatography and gas chromatography, and more particularly to a chromatogram analysis method (hereinafter referred to as an analysis method) and a chromatographic apparatus used therein, and data used in the analytical method and the chromatographic apparatus. The present invention relates to a processing device.
[0002]
[Prior art]
Prior to the description of the conventional technology, technical terms are defined.
In the present application, the convolution is defined as a convolution integral of two functions x (t) and y (t). In general, an output waveform from a measuring instrument is an input waveform and an apparatus function (an impulse response in an electric circuit). , And a convolution with a spectroscope).
Deconvolution is defined as the inverse operation of convolution. Specifically, this corresponds to solving an integral equation represented by convolution.
The reconvolution means that a waveform deconvolved using a device function, that is, a function h (t) representing a spread is convoluted again using the function h (t).
[0003]
Conventionally, generally, output waveform data from a measuring instrument is distorted by the dynamic characteristics of the measuring instrument, that is, the influence of the device function.
The same applies to a chromatographic apparatus, in which distortion occurs on the chromatogram and peaks overlap on the chromatogram due to the apparatus function, that is, the function h (t) representing the spread.
In such a case, it has been difficult to accurately determine the peak area of each chromatogram.
[0004]
Conventionally, when the above-mentioned distortion is applied, the overlap of the peaks is decomposed by a deconvolution process using a synthesis separation method based on a normal curve fitting method, and is subjected to data processing. It has been proposed to quantify.
A general solution to this problem is described in "Waveform Data Processing for Scientific Measurement", edited by Shigeo Minami, published by CQ Publishing Co., Ltd. (1986).
Also, many improvement techniques have been proposed, for example, the technique described in JP-A-62-17465 and the technique described in JP-A-63-151851.
[0005]
In the technique described in Japanese Patent Application Laid-Open No. Sho 62-17465, the half-width on the rising side is approximated by Gaussian, which is a predetermined multiple of the distance from the rising point to the end point. In the case of peak waveforms, there are few cases where such approximation can be used for replacement.
Although this technique is effective in reducing the storage capacity, it has a drawback that it is not necessarily effective in accurately obtaining the peak area when the peak waveforms overlap.
[0006]
In the technique described in Japanese Patent Application Laid-Open No. 63-151851, the overlapped peak is divided into two Gaussians, which has been so-called vertical division and skimming division, and is approximated by a method of least squares. Although the deconvolution process is performed by the technique, there is a disadvantage that a physically meaningful solution cannot be obtained.
[0007]
Further, in the techniques described in U.S. Pat. Nos. 4,807,148 and 4,941,101, a numerical analysis technique is used in the deconvolution processing, and there are Jansson's method, Jacobi method, Gauss-Seidel method, and the like. Although the number of stored words is advantageous, there is a disadvantage in that the number of discontinuous points in the solution is increased.
[0008]
[Patent Document 1]
JP-A-62-17465
[Patent Document 2]
JP-A-63-151851
[Patent Document 3]
U.S. Pat. No. 4,807,148
[Patent Document 4]
U.S. Pat. No. 4,941,101
[Non-patent document 1]
Minami Shigeo, "Waveform Data Processing for Scientific Measurement", CQ Publishing Co., Ltd., (1986)
[0009]
[Problems to be solved by the invention]
The various conventional deconvolution processes described above provide respective solutions for the deconvolution process when the chromatogram is distorted and the synthesis and separation when the chromatogram has an overlapping peak waveform. However, in the case where an overlapping peak waveform is generated due to the distortion, it has been studied in detail up to the point of separating the overlapping waveform in consideration of the distortion and has not been dealt with.
[0010]
Further, the above-mentioned conventional synthesis separation method for the overlap of chromatogram peak waveforms has the following problem.
1. Assuming a non-linear function such as Gaussian, the validity of selecting that function is questionable.
2. Introducing multi-channel information such as the UV spectrum requires an expensive detector such as a photodiode array detector.
3. In the peak decomposition, since the least squares method is used, a long calculation time is required until a solution is obtained, and a physically meaningful solution may not be obtained in some cases.
[0011]
SUMMARY OF THE INVENTION The present invention has been made to solve the above-mentioned problems of the prior art, and improves the separation by correcting waveform distortion of a chromatogram and sharpening overlapping peak waveforms. The present invention aims to provide a chromatogram analysis method in which the work area is more accurately determined and the accuracy of quantitative analysis is improved.
[0012]
[Means for Solving the Problems]
In order to achieve the above object, a configuration of the present invention relating to a chromatogram analysis method includes a chromatogram analysis method in which a sample to be analyzed having passed through a separation column is detected by a detector, and a chromatogram obtained by the detection is analyzed. At
A step of introducing a standard sample not retained by the separation column into the separation column, separately from introducing the sample to be analyzed into the separation column,
Identifying a peak of the standard sample from the chromatogram;
Setting a start point and an end point of a certain time of the peak of the sample to be analyzed from the chromatogram,
Performing a deconvolution process on the data within the set time interval using a peak waveform of the identified standard sample.
In the chromatogram analysis method described above,
When the peaks after performing the deconvolution processing are a plurality of peaks and remain in an overlapping state, perform vertical division or area division using a nonlinear function on the overlapping peaks. It is characterized by.
[0013]
In order to achieve the above object, another configuration of the present invention according to the chromatogram analysis method includes detecting a sample to be analyzed that has passed through a separation column with a detector, and standard deviation of a peak in the chromatogram obtained by the detection. In a chromatogram analysis method of performing analysis using a Gaussian function having σ as a coefficient,
Setting an arbitrary value to the standard deviation σ of the Gaussian function;
Performing a deconvolution process on the chromatogram;
Regarding the chromatogram after the deconvolution process, the chromatogram does not diverge, does not fall below the negative lower limit, or an unseparated peak forms a valley between peaks, and satisfies one of the following criteria. Determining whether or not
Adjusting the standard deviation value of the Gaussian function when the determination is unsuccessful, and performing deconvolution processing again.
In the chromatogram analysis method described above,
Using the spreading function used in the deconvolution processing, a waveform after the reconvolution processing is displayed on the peak separated by the deconvolution processing so as to overlap the peak waveform before the deconvolution processing. .
[0014]
More specifically, the chromatogram peaks are spread for various reasons. For example, when a sample flows through a flow path such as a pipe or a flow cell of a detector, the sample zone is expanded. This is called extra-column spreading.
The contribution due to the spread typified by the spread outside the column is convoluted to the original peak waveform and appears in the detected waveform.
That is, if the spread function is known, the original peak waveform can be calculated by performing a deconvolution operation on the detected waveform.
[0015]
There are various methods for determining the spread function h (t).
When the methods are roughly classified into two, there are a method of experimentally obtaining from detected data and a method of assuming an analytical mathematical h (t) having some parameters characterizing a waveform.
The former can be considered as a method in which h (t) is known, and the latter can be considered as a method in which h (t) is unknown.
[0016]
As a method for obtaining the h (t) from the measurement, for example, a method of exchanging the column with a joint, injecting a sample, and regarding the waveform measured by the detector as h (t) in order to obtain the spread outside the column is available. is there.
If the extension in the column that does not depend on the type of components in the sample is also added, a sample not retained can be injected, and the peak waveform can be regarded as h (t).
[0017]
Other methods for obtaining peaks that are not retained include the use of a solvent having an extremely high elution ability as an eluent, and the preparation of a column packing material that does not modify a functional group so as not to retain it.
It is also effective to add a component not retained in the standard sample to obtain a non-retained peak from the chromatogram of the standard sample.
[0018]
If h (t) is unknown, the spread function h (t) corresponding to a part of the peak width is deconvoluted under the following conditions, and a sharper peak is obtained. It can be accurately quantified based on the peak area from the processed chromatogram.
[0019]
The precondition is that the detected peak waveform is originally a convolution operation of a spread function with respect to a delta function or a very sharp rectangular density distribution function. Can be
Here, the spread function is almost Gaussian, which has a waveform slightly modified in a tailing or reading manner.
In other words, it can be considered that the spreading function itself is a convolution operation of the tailing or reading correction function for Gaussian.
[0020]
Further, Gaussian has an additive property that when two Gaussians are subjected to a convolution operation, they are also Gaussians.
As described above, by performing the deconvolution operation on the Gaussian constituting a part of the spread function, it is possible to calculate a peak approaching a sharper original waveform.
[0021]
[Action]
The function of each of the above technical means is as follows.
According to the configuration of the present invention, a known h (t) or Gaussian is rationally selected from the detected chromatogram, and deconvolution operation is performed to obtain a peak sharper than the overlapping peak of the detected chromatogram. When the peaks are separated and separated into chromatograms having peaks, and the separated and independent peaks are integrated, the peak area is accurately obtained, and accurate quantification can be performed.
In the deconvolution operation, unlike the least square method of a nonlinear function, for example, a physically meaningful solution can always be obtained in a short calculation time.
[0022]
In the deconvolution operation, when the overlapping peaks are separated and separated, the chromatogram is observed to judge the validity of the operation, but it requires skill. Therefore, the independent peaks can be subjected to the reconvolution operation using the spread function h (t) used for the deconvolution operation, and each peak constituting the first overlapping peak can be reproduced. By overlaying each peak after the reconvolution operation and the first overlapping peak, the validity of the deconvolution operation can be accurately determined.
[0023]
Further, according to the configuration of the present invention, an expensive detector such as a photodiode array detector is omitted, and the chromatogram analysis method can be performed with a simple configuration.
[0024]
【Example】
[Example 1]
Hereinafter, an embodiment according to the present invention will be described with reference to FIGS.
FIG. 1 is a block diagram showing a configuration of a microbore HPLC system used in a chromatogram analysis method (hereinafter, referred to as an analysis method) according to one embodiment of the present invention, and FIG. 2 shows a chromatogram in the analysis method of FIG. FIG. 3 is a schematic explanatory view of a spread mechanism of a chromatogram in the analysis method of FIG. 1, FIG. 4 is an explanatory view of deconvolution calculation in the analysis method of FIG. 1, and FIG. 5 is an analysis method of FIG. FIG. 6 is a plot plot of the chromatogram obtained by the analysis method shown in FIG.
[0025]
FIG. 7 is a diagram showing a chromatogram of amino acid analysis by the analysis method of FIG. 1, FIG. 8 is a diagram showing a chromatogram of glycohemoglobin analysis by the analysis method of FIG. 1, and FIG. 9 is a diagram of FIG. FIG. 10 is a processing flowchart in the case where the spread function is unknown in the analysis method, FIG. 10 is a diagram showing a chromatogram of catecholamine analysis in the analysis method of FIG. 1, and FIG. 11 is a reconvolution in the analysis method of FIG. FIG. 18 is an explanatory diagram of superposition by a calculation operation.
[0026]
In order to carry out the chromatogram analysis method according to the present invention, a microbore HPLC system is used. FIG. 1 shows the overall configuration of the microbore HPLC system.
In FIG. 1, 40 is a pump, 41 is a column, 42 is a detector, 43 is a sampler, 44 is a controller, 45 is a data analyzer, 46 is a CRT display, 47 is a printer, 49 is an eluent, and 51 is an eluent. Is an unknown sample, 52 is a standard sample, and 54 is an input unit.
[0027]
Except for the data analysis unit 45 (described later in [Embodiment 2]), each unit constituting the present system uses a normal chromatogram analysis method. Therefore, detailed description is omitted because it is complicated, and will be described briefly.
[0028]
In the microbore HPLC system, the spread outside the column greatly affects the separation and cannot be ignored. The
The sample is sent to the
[0029]
The detected chromatogram is stored in the
FIG. 2 is the detected chromatogram. Since a non-retained component is added to the
[0030]
With reference to FIG. 3, a mechanism for forming a peak waveform will be described.
FIG. 3 is an explanatory diagram of a mechanism for expanding a chromatogram in the analysis method according to the embodiment of FIG.
In FIG. 3, the same reference numerals as those in FIG. 1 denote the same parts, and a detailed description thereof will be omitted.
The upper diagram (a) is a block diagram showing a sample flow channel, and the lower (b) diagram shows the sample concentration in the sample flow channel of the upper (a) diagram.
[0031]
A sample containing a single component is injected into the channel by the
When this rectangular distribution diagram 71 passes through the
[0032]
Further, when the concentration distribution of Gaussian is detected by the
[0033]
If this is described by a mathematical expression, it can be represented by a model as shown in Expression (2).
(Equation 2)
In the above equation, the detected waveform at time t is y (t), the original rectangular waveform at the time of sample injection is r (t), and the spread function is h (t). In the above equation, the detected waveform y (t) is a convolution of the original rectangular waveform r (t) and the spread function h (t).
[0034]
The spread of the peak shape may be described above. However, in consideration of the retention time, it is considered that another Dirac delta function is subjected to the convolution operation. This is shown in equation (3).
(Equation 3)
[0035]
In the above equation, EMG (t) is Exponentially Modified Gaussin, and σ, τ, t R Is a parameter of EMG (t), δ (t) is a delta function, and m (t) is a correction function.
The mechanism of spreading is described in "DYNAMICS OF CHROMATOGRAPHY", J. Amer. C. Giddings, published by Marcel Dekker, New York, 1965;
[0036]
The spreading function h (t) is formed by the influence from both the spreading inside the column and the spreading outside the column. That is, a convolution of the spread function inside the column and the spread function outside the column.
[0037]
Focusing on the non-holding peak 61 in FIG. 2 and referring to FIG. 3, a
[0038]
On the other hand, the holding peak 62 is also expanded by a similar mechanism to form a peak. However, since the holding peak 62 is held in the
That is, the spread function affecting the peak 62 is a convolution of the spread function common to the peak 61 and the other surplus spread functions.
[0039]
In this case, it is possible to remove the influence of the spread function of the peak 61 from the spread function of the peak 62 and obtain only the surplus spread function. This process is called a deconvolution process.
[0040]
The deconvolution operation is effective for overlapping peaks such as
This will be described with reference to FIG. FIG. 4 illustrates the deconvolution operation for the overlapping peaks shown in FIG.
[0041]
When the spread function of the peak 61 waveform is deconvoluted with respect to the
Actually, a normalization function obtained by normalizing the peak waveform 61 to the
[0042]
Strictly speaking, this normalization function is a standardized convolution of the rectangular function r (t) at the time of injection of (Equation 2) and the spread function h (t) affecting the peak waveform 61. .
Actually, since the peak function other than the peak waveform 61 is also affected by the rectangular function, it does not matter if the deconvolution operation is performed by the normalization function of the peak 61 waveform.
[0043]
The above-described deconvolution procedure of FIG. 4 will be described with reference to the block diagram of the micro-bore HPLC system of FIG. 1 and the flowchart of the detection waveform processing shown in FIG.
First, the
[0044]
Next, in
[0045]
In
In
[0046]
In
In
[0047]
In
If peaks still overlap in the chromatogram after this processing, the area is calculated by a conventional vertical division or the like.
Here, it is also possible to divide the area by the least square method using a nonlinear function such as Gaussian or EMG. Here, the quantitative value is output to the
In
[0048]
Although the spread function for the deconvolution processing is substituted by the normalization function of the peak waveform of the non-retained component in the
Instead of adding to the standard sample, only the non-retained component is injected in order to obtain a spread function in advance, and a peak waveform can be obtained.
[0049]
In the case of reversed-phase chromatography, an eluent containing 100% of an organic solvent such as methanol, acetonitrile, THF, chloroform or the like is flowed, and a component having a low retention is injected to obtain a non-retention peak.
In ion exchange chromatography, for example, in the case of an amino acid analyzer, a regenerated solution having a high pH is allowed to flow, and phosphoserine or the like having a low retention is injected to obtain a peak waveform.
[0050]
In addition, a non-retained peak can be obtained by manipulating the column packing material. In the case of ion-exchange chromatography, it is also possible to prepare a packing consisting only of a carrier without modification of the ion-exchange group to obtain a non-retained peak.
[0051]
An example of the post-column reaction LC is an amino acid analyzer. This system cannot ignore the spread outside the column. Spread outside the column can be estimated by exchanging the column with a joint and estimating the waveform of the unretained peak as a spread function.
[0052]
In the chromatogram detected as shown in FIG. 7, the overlapping portion of threonine (Thr) and serine (Ser) is 14%, whereas in the chromatogram subjected to the deconvolution treatment, it is 4% or less than 1/3. Improved.
This processed chromatogram corresponds to a chromatogram in the case where detection was possible immediately after the column in the post-column reaction LC.
[0053]
As an example of removing the common part of the spread function, a glycated hemoglobin analyzer processed for 3.3 minutes will be described.
The spread function can be estimated from the peak waveform of HbF shown on the left side of FIG. This is because the peak width of HbA1c is considered to be slightly wider than the width of HbF in consideration of the order of the retention times.
[0054]
The waveform of HbF, which is an isolated peak, is normalized from the detection chromatogram on the left side of FIG. 8 to obtain a spread function. If the overlapping peaks of the unstable type and the stable type HbA1c are subjected to deconvolution processing using the standardized function, a chromatogram with improved separation as shown on the right side of FIG. 8 is obtained.
[0055]
In this case, it is not necessary to directly display the processed chromatogram, but the deconvolution treatment is a very effective technique in a glycohemoglobin analyzer which is a general HPLC in order to improve the quantitative accuracy.
[0056]
Next, a description will be given of processing to be used before the operation in order to correctly perform the deconvolution operation. Noise 67 and drift 68 shown in FIG. 2 interfere with the deconvolution operation.
If these are deconvoluted as they are, there will be problems such as amplification as pseudo peaks and improper processing of peaks to be deconvoluted.
[0057]
First, as the noise removal processing, a method of reducing the number of data points of a detected chromatogram by adding some data points is generally adopted.
In addition, a smoothing method and a known method for removing high-frequency components can be used.
[0058]
As the drift removal processing, generally, a method of connecting a valley between peaks with a straight line or forcibly drawing a straight line horizontally and subtracting a base line from a detected chromatogram is employed.
Other known methods for removing low-frequency components, methods for removing low-order polynomial components using calculus, methods for searching for a base line from differential feature points, methods using a neural network, etc. There is.
[0059]
Up to this point, the spread function has been obtained from the measured data. However, even when the spread function is unknown, the Gaussian can be deconvoluted.
As described above, the peak waveform of the detected chromatogram can be represented by the convolution of the spread function h (t) and the rectangular waveform r (t), as shown in equation (2).
[0060]
This spreading function h (t) is described in “DYNAMICS OF CHROMATOGRAPHY” J. C. According to Giddings, a Gaussian-like shape is obtained regardless of the extent of physical or chemical interaction within the column, even if the extent outside the column is not negligible.
[0061]
Even if h (t) is slightly distorted from Gaussian in reading or tailing, the convolution of Gaussian G (t) and correction function m (t) as shown in the following equation (4). Can be described as
(Equation 4)
[0062]
In the above equation, G (t) is a Gaussian, m (t) is a correction function, and σ is a standard deviation.
As an example, if h (t) is a tailing peak that can be represented by EMG (t), a correction function is obtained analytically and mathematically as shown in the above equation (3).
[0063]
The overlapping peak can be expressed as in equation (5).
(Equation 5)
Where σ 1 , Σ 2 Is the standard deviation.
[0064]
Since Gaussian is additive as described above, it can be decomposed into two Gaussian convolutions, and therefore has the relationship as shown in the following equation (6).
(Equation 6)
[0065]
That is, the two peak waveforms expressed by the equation (5) can be subjected to deconvolution operation using Gaussian.
At this time, the Gaussian is the standard deviation σ of the two peaks. 1 And σ 2 It must be a narrow waveform with a smaller standard deviation σ0.
[0066]
Thus, a peak narrower than the detected overlapping peak can be obtained.
Finally, y (t) is converted to Gaussian G (t, σ 0 ), The chromatogram is converted to a narrower Gaussian with the original width. Osamu Mano It can be expressed by convolution with a positive function.
[0067]
(Equation 7)
As shown in [] above, σ 1 , Σ 2 A narrower peak is obtained. More specifically, a rectangular waveform at the time of injection is subjected to a convolution operation by this correction function.
[0068]
The processing when the spread function h (t) is unknown will be described with reference to FIG.
The process starts at
In
[0069]
In the automatic case, in step 81, a Gaussian having a standard deviation slightly smaller than the minimum width is set as a spread function by looking at the peak width of the detected chromatogram as a whole. This is standardized to
Here, when the peak width differs between the peak having a short retention time and the peak having a long retention time more than twice, the deconvolution operation can be performed by dividing the time.
[0070]
Here, Fourier transform can be used to automatically find the peak width. In
[0071]
If the determination is successful, the process proceeds to step 79. If the determination is unsuccessful, the process returns to step 81 again, and the deconvolution processing is performed by referring to the determination and adjusting the standard deviation.
If the judgment is successful, display of the processing chromatogram, display of the quantitative result, and the like are performed in
[0072]
In the case of manual operation, an input operation is performed in step 74. Generally, the standard deviation of the Gaussian to be deconvolved is input. In addition, the start and end points of the processing section can be input. It is also possible to search for an optimal Gaussian standard deviation in order to improve the separation of specific overlapping peaks.
[0073]
In this case, the initial value of the standard deviation, the number, the size, the width, and the holding time of the estimated peak can be input as reference values. As the limit value, the number of peaks can be input up to three, or the width of the overlapping peaks can be set to 0.2 to 0.4 minutes.
[0074]
In a
If not, the reason is displayed at
The process ends at
[0075]
The raw 10-minute chromatogram obtained by the catecholamine analyzer is shown in the upper part (a) of FIG.
This is an example of analyzing a biological sample such as plasma or urine. Many peaks often appear when analyzing biological samples. When this is deconvoluted with an appropriate Gaussian, a chromatogram as shown in the middle part (b) is obtained.
Each peak is completely separated up to the base line.
[0076]
However, even if this is displayed as it is, it is difficult to understand the meaning of the processed chromatogram.
Rather, this may be changed to a known display method shown in (c) below.
It is also preferable that only the quantitative value is output and the operator is not conscious of the deconvolution processing.
[0077]
As mentioned above, although one Example of this invention was described, this invention is not limited to the said Example, Various modifications are possible.
Another method for analyzing a chromatogram according to the present invention will be described.
FIG. 11 is a view for explaining the overdrawing by the deconvolution and the reconvolution in the analysis method according to the present invention.
[0078]
FIG. 11A shows the detected original waveform. The original waveform is subjected to deconvolution processing to separate and independent each waveform peak. This state is shown in FIG.
The separated and independent waveform peaks are reconvolved by the same spread function h (t) as when deconvolved.
The reconvolution-computed waveform is overlaid on the detected original waveform. FIG. 11 (3) shows this state.
[0079]
[Example 2]
Another embodiment of the chromatographic analyzer according to the present invention will be described. FIG. 12 is a block diagram showing a configuration of a chromatographic analyzer according to another embodiment of the present invention.
The overall configuration of the chromatographic analyzer according to the present invention has been described in [Example 1] with respect to the micro-bore HPLC system of FIG. explain.
[0080]
The
The noise drift elimination mechanism, the quantitative calculation mechanism, and the determination mechanism among the above-described configurations are well-known technologies, but will be described together in order. In addition, each said mechanism part is connected by the bus.
[0081]
The noise drift elimination mechanism includes a smoothing mechanism, a baseline calculation mechanism, a fast Fourier transform mechanism, a neural network mechanism, a differentiation / integration mechanism, and a fuzzy inference mechanism. Fits the measured data group of the chromatogram from the
[0082]
The fast Fourier transform mechanism removes noise that is a high-frequency component and baseline fluctuation that is a low-frequency component, respectively, based on the frequency distribution.
The neural network mechanism learns the drift shape in advance and removes the drift from the measured data according to a predetermined coupling constant.
[0083]
The differentiation / integration mechanism removes a baseline component represented by a low-order polynomial by differentiation / integration. The arranged measurement data points are smoothed by the fuzzy inference mechanism according to a predetermined coupling constant. Thus, the noise and drift of the chromatogram from the
[0084]
The spread function analysis mechanism includes a Gaussian generation mechanism, an
[0085]
The deconvolution operation mechanism includes a time section setting mechanism, a spreading function setting mechanism, and a deconvolution processing mechanism.The time section setting mechanism sets a section for deconvolution processing, and the spreading function setting mechanism includes: A spread function is set from the generated actual or analytical mathematical spread function.
The deconvolution processing mechanism performs deconvolution processing on the set overlapping waveform in the section to be subjected to the deconvolution processing using the spreading function, and separates and separates the overlapping waveform.
[0086]
The quantitative calculation mechanism includes an integration mechanism, and the integration mechanism calculates an area of each of the separated and independent peak waveforms. In this way, the quantification is calculated based on the peak area or height. The quantitative result is output to the
[0087]
The determination mechanism that performs various determinations is a discriminant calculation mechanism that calculates the value of a discriminant used to determine the success of deconvolution processing, and is limited when searching for a Gaussin used for deconvolution processing. A limit value calculation mechanism for calculating a value, and a spread function search mechanism for searching for an optimum Gaussin according to a predetermined coupling coefficient, and in each of the processing steps, support for enabling each mechanism to perform a predetermined function. I have.
[0088]
The functions of the
[0089]
【The invention's effect】
As described above in detail, according to the configuration of the present invention, the distortion of the chromatogram waveform is corrected, the overlapping peak waveform is sharpened, and the separation is improved, and the peak area is reduced. It is possible to provide a chromatogram analysis method that is more accurately obtained and has improved quantitative analysis accuracy.
[Brief description of the drawings]
FIG. 1 is a block diagram showing a configuration of a microbore HPLC system used for an analysis method according to one embodiment of the present invention.
FIG. 2 is a diagram showing a chromatogram in the analysis method of FIG. 1;
FIG. 3 is a schematic explanatory view of a spreading mechanism of a chromatogram in the analysis method of FIG. 1;
FIG. 4 is an explanatory diagram of a deconvolution operation for an overlapping peak of a chromatogram in the analysis method of FIG. 1;
FIG. 5 is a flowchart of a processing step according to the analysis method of FIG. 1;
FIG. 6 is an overlaid plot of a chromatogram in the analysis method of FIG. 1;
FIG. 7 is a diagram showing a chromatogram of amino acid analysis by the analysis method of FIG. 1;
8 is a diagram showing a chromatogram of glycohemoglobin analysis by the analysis method of FIG.
FIG. 9 is a processing flowchart in the case where a spread function is unknown in the analysis method of FIG. 1;
10 is a diagram showing a chromatogram of catecholamine analysis in the analysis method of FIG.
FIG. 11 is an explanatory diagram of superposition by a reconvolution operation in the analysis method of FIG. 1;
FIG. 12 is a block diagram showing a configuration of a chromatographic analyzer according to another embodiment of the present invention.
[Explanation of symbols]
40 ... pump
41 ... column
42 ... Detector
43 ... Sampler
44 ... Control unit
45 Data analysis unit
46 CRT display
47… Printer
49 ... Eluent
51 ... unknown sample
52: Standard sample
54 ... input section
Claims (6)
前記ガウシアンの標準偏差σに任意の値を設定するステップと、
前記クロマトグラムに対してディコンボリューション処理を行うステップと、
当該ディコンボリューション処理後のクロマトグラムについて、
クロマトグラムが発散していないか、
負の下限値を下回っていないか、
不分離ピークがピーク間の谷を形成しているか、
のいずれかの基準を満たしているかどうかの判定を行うステップと、
当該判定で不合格の場合、前記ガウシアンの標準偏差値を調整し、再度ディコンボリューション処理を行うステップと、を有することを特徴とするクロマトグラム解析方法。A chromatogram analysis method in which a sample to be analyzed that has passed through a separation column is detected by a detector, and a peak in a chromatogram obtained by the detection is represented and analyzed using a Gaussian having a standard deviation σ as a coefficient as a spreading function. ,
Setting an arbitrary value to the standard deviation σ of the Gaussian;
Performing a deconvolution process on the chromatogram;
For the chromatogram after the deconvolution processing,
Is the chromatogram diverging?
Is not below the negative lower bound,
Whether the unseparated peaks form valleys between the peaks,
Determining whether any of the criteria is satisfied; and
And a step of adjusting the standard deviation value of the Gaussian and performing the deconvolution process again if the judgment is unsuccessful.
前記ディコンボリューション処理によって分離されたピークに対して、前記ディコンボリューション処理に用いた拡がり関数を用いてリコンボリューション処理を行い、リコンボリューション処理後の波形を前記ディコンボリューション処理前のピーク波形と重ねて表示することを特徴とするクロマトグラム解析方法。In the chromatogram analysis method according to claim 1,
The peak separated by the deconvolution process is subjected to a reconvolution process using the spread function used for the deconvolution process, and a waveform after the reconvolution process is displayed by being superimposed on a peak waveform before the deconvolution process. A chromatogram analysis method.
ガウシアンの標準偏差σを設定する前に、前記クロマトグラムからノイズ及びドリフトを除去するステップを含むことを特徴とするクロマトグラム解析方法。In the chromatogram analysis method according to claim 1,
A step of removing noise and drift from the chromatogram before setting the Gaussian standard deviation σ.
前記拡がり関数h(t)は、
h(t)=G(t)* m(t)
(ここで、G(t)は上記ガウシアン、m(t)は修正関数)
前記ガウシアンは、(1)式で表されることを特徴とするクロマトグラム解析方法。
The spreading function h (t) is
h (t) = G (t) * m (t)
(Where G (t) is the Gaussian and m (t) is a correction function)
The said Gaussian is represented by Formula (1), The chromatogram analysis method characterized by the above-mentioned.
前記標準偏差σは、検出されたクロマトグラムのピーク幅よりも小さな標準偏差を設定することを特徴とするクロマトグラム解析方法。In the chromatogram analysis method according to claim 1,
A method for analyzing a chromatogram, wherein the standard deviation σ is set to a standard deviation smaller than a peak width of a detected chromatogram.
前記クロマトグラムが拡がり関数が未知のものであることを特徴とするクロマトグラム解析方法。In the chromatogram analysis method according to claim 1,
A chromatogram analysis method, wherein the chromatogram has a spread function whose function is unknown.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002297130A JP3583771B2 (en) | 2002-10-10 | 2002-10-10 | Chromatogram analysis method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002297130A JP3583771B2 (en) | 2002-10-10 | 2002-10-10 | Chromatogram analysis method |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP24933593A Division JP3391862B2 (en) | 1993-10-05 | 1993-10-05 | Chromatogram analysis method |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2003161725A JP2003161725A (en) | 2003-06-06 |
JP3583771B2 true JP3583771B2 (en) | 2004-11-04 |
Family
ID=19197260
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002297130A Expired - Fee Related JP3583771B2 (en) | 2002-10-10 | 2002-10-10 | Chromatogram analysis method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3583771B2 (en) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7759130B2 (en) * | 2003-11-21 | 2010-07-20 | Eisai R & D Management Co., Ltd. | Mass spectrometric quantitation method for biomolecules based on metabolically labeled internal standards |
JP5068498B2 (en) | 2006-08-31 | 2012-11-07 | 株式会社日立ハイテクノロジーズ | Chromatographic analyzer |
GB2563748B (en) * | 2016-01-28 | 2022-01-19 | Hitachi High Tech Corp | Base sequence determination apparatus, capillary array electrophoresis apparatus, and method |
JP6631373B2 (en) * | 2016-04-13 | 2020-01-15 | 株式会社島津製作所 | Data processing device |
JP6911543B2 (en) * | 2017-06-01 | 2021-07-28 | 東ソー株式会社 | Baseline determination method |
JP7119602B2 (en) * | 2018-06-07 | 2022-08-17 | 東ソー株式会社 | How to calculate chromatogram similarity |
EP3865872A4 (en) * | 2018-11-19 | 2022-08-10 | Canon Kabushiki Kaisha | Information processing device, information processing device control method, program, calculation device, and calculation method |
JP7282801B2 (en) * | 2018-11-20 | 2023-05-29 | 株式会社日立ハイテク | Analysis apparatus having multiple liquid chromatographs and analysis method thereof |
CN112415078B (en) * | 2020-11-18 | 2023-11-14 | 深圳市步锐生物科技有限公司 | Mass spectrum data spectrogram signal calibration method and device |
CN113607867B (en) * | 2021-07-23 | 2024-06-11 | 清华大学合肥公共安全研究院 | Double-overlap spectrum peak analysis method based on peak body mapping |
CN114444542A (en) * | 2022-01-10 | 2022-05-06 | 中国科学院苏州生物医学工程技术研究所 | Liquid chromatography peak noise estimation method, device, storage medium and system |
CN117907511B (en) * | 2024-03-20 | 2024-06-14 | 浙江灵析精仪科技发展有限公司 | Automatic analysis method and device for multi-component overlapping peaks and electronic equipment |
-
2002
- 2002-10-10 JP JP2002297130A patent/JP3583771B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2003161725A (en) | 2003-06-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3391862B2 (en) | Chromatogram analysis method | |
JP3583771B2 (en) | Chromatogram analysis method | |
US7110886B2 (en) | Method and apparatus for analyzing multi-channel chromatogram | |
JP4973628B2 (en) | Chromatograph mass spectrometry data analysis method and apparatus | |
US6787761B2 (en) | Median filter for liquid chromatography-mass spectrometry data | |
US5969228A (en) | Method and devices for chromatographic pattern analysis employing chromatographic variability characterization | |
JP3130139B2 (en) | Chromatogram analysis method and chromatographic device | |
JP5757264B2 (en) | Chromatographic mass spectrometry data processor | |
JP2012177568A (en) | Data processing device, data processing method, and data processing program | |
Allen et al. | Semi-automated alignment and quantification of peaks using parallel factor analysis for comprehensive two-dimensional liquid chromatography–diode array detector data sets | |
JP7310839B2 (en) | PEAK ANALYSIS METHOD AND WAVEFORM PROCESSING DEVICE | |
US9347921B2 (en) | Method for filtering a chromatogram | |
JP3550479B2 (en) | Chromatography analysis method and chromatography analyzer | |
JP2000131284A (en) | Chromatograph mass spectrometer | |
Zhang et al. | Resolution of overlapping capillary electrophoresis peaks by using chemometric analysis: improved quantification by using internal standard | |
JP3025145B2 (en) | Multichannel chromatogram analysis method and apparatus | |
Wenig | Post-optimization of Py-GC/MS data: A case study using a new digital chemical noise reduction filter (NOISERA) to enhance the data quality utilizing OpenChrom mass spectrometric software | |
JP2000088833A (en) | Analyzing device for chromatogram | |
JP7334788B2 (en) | WAVEFORM ANALYSIS METHOD AND WAVEFORM ANALYSIS DEVICE | |
US20230304975A1 (en) | Peak shape estimation device and peak shape estimation method | |
JP2002286704A (en) | Data processing device for chromatograph | |
JP7375928B2 (en) | Peak tracking device, peak tracking method and peak tracking program | |
Solis et al. | Accelerated multiple‐pass moving average: A novel algorithm for baseline estimation in CE and its application to baseline correction on real‐time bases | |
CN116029965A (en) | Data processing method, data processing system and storage medium | |
Chilcote et al. | Use of automatic digital data acquisition and on-line computer analysis in high-resolution liquid chromatography |
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: 20040727 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20040729 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20070806 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080806 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080806 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090806 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100806 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100806 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110806 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120806 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130806 Year of fee payment: 9 |
|
LAPS | Cancellation because of no payment of annual fees |