JP7216225B2

JP7216225B2 - クロマトグラムデータ処理装置、クロマトグラムデータ処理方法、クロマトグラムデータ処理プログラム及び記憶媒体

Info

Publication number: JP7216225B2
Application number: JP2021561190A
Authority: JP
Inventors: 健一郎鮫島
Original assignee: Alps Electric Co Ltd; Alps Alpine Co Ltd
Current assignee: Alps Alpine Co Ltd
Priority date: 2019-11-27
Filing date: 2020-09-29
Publication date: 2023-01-31
Anticipated expiration: 2040-09-29
Also published as: WO2021106356A1; JPWO2021106356A1

Description

本発明は、クロマトグラムデータ処理装置、クロマトグラムデータ処理方法、クロマトグラムデータ処理プログラム及び記憶媒体に関する。

クロマトグラフを使用して試料中の成分を測定する際、クロマトグラムのピーク部分に複数個の成分が重なり合う部分は、ガウシアン(Gaussian)関数や指数関数により修飾されたガウシアン(ＥＭＧ：Exponentially Modified Gaussian)関数を用いてクロマトグラムをフィッティングすることにより、クロマトグラムをそれぞれの成分に分解し、それぞれの試料成分を測定するデータ処理装置がある（例えば、特許文献１参照）。

特許文献１には、クロマトグラムにより試料の成分を測定するクロマトグラフのデータ処理装置において、試料から得られるデータの重なりピークに対してガウシアン関数又はＥＭＧ関数を用いてフィッティングして見積もった非線形パラメータと、標準試料から同様にして得られた非線形パラメータとを比較し、両パラメータが近似している場合にはフィッティングが妥当であることを出力し、試料から得られた非線形パラメータを用いてデータの重なりピークを個々の成分に分解するデータ処理装置が開示されている。

日本国特開２０００－１３１３０５号公報

特許文献１のデータ処理装置では、クロマトグラム内の隣接するピーク同士の重なりが小さいピークについては、ガウシアン関数及び指数関数により修飾されたガウシアン関数のいずれを用いても、精度良くフィッティングすることができる。しかし、クロマトグラム内の隣接するピーク同士の重なりが大きいピークの場合、隣接するピーク同士が互いに影響し合うため、ピークのリーディング又はテーリングが発生する場合がある。その結果、隣接するピーク同士の重なりが大きいピークについては、指数関数により修飾されたガウシアン関数であっても、精度良くフィッティングできない場合がある。

即ち、特許文献１のデータ処理装置を用いても、フィッティング前後の差が大きくなるため、フィッティング後の結果を妥当と見なすことが困難となる可能性が高くなる。

本発明の一態様は、クロマトグラム内のピークを個々の成分ごとにより高精度に分解できるクロマトグラムデータ処理装置を提供することを目的とする。

本発明の一態様に係るクロマトグラムデータ処理装置は、測定対象試料に含まれる複数の成分が分離カラムで分離され、分離された前記成分を検出して得られるクロマトグラムにより前記成分を測定するクロマトグラムデータ処理装置であって、時間と信号強度とをそれぞれ軸として、複数の前記成分に由来するピークを有する前記クロマトグラムを得るクロマトグラム取得部と、前記クロマトグラム内の複数の前記ピークに対して、ガウシアン関数又は指数関数で修飾されたガウシアン関数と、ピアソン関数とを用いて、フィッティング処理を行い、複数の前記ピークを分解するフィッティング部と、前記フィッティング処理を行った後の前記クロマトグラムから前記成分を判定する判定部と、を備える。

本発明の一態様に係るクロマトグラムデータ処理装置は、クロマトグラム内のピークを個々の成分ごとにより高精度に分解できる。

一実施形態に係るクロマトグラムデータ処理装置の機能を示すブロック図である。クロマトグラムの一例を示す図である。微分データの一例を示す図である。ピーク面積の計算の一例を示す図である。ピーク面積の比較の一例を示す図である。ピークの分離度が１．０未満の場合の２つのピークの重なり具合を示す説明図である。ピークの分離度が１．０以上の場合の２つのピークの重なり具合を示す説明図である。クロマトグラムデータ処理装置のハードウェア構成を示すブロック図である。本実施形態に係るクロマトグラムデータ処理方法を説明するフローチャートである。クロマトグラムデータ処理装置を適用した分析装置の一例を示す図である。実施例１で用いたクロマトグラムである。実施例１－１のフィッティング処理した結果を示す図である。比較例１－１のフィッティング処理した結果を示す図である。比較例１－２のフィッティング処理した結果を示す図である。比較例１－３のフィッティング処理した結果を示す図である。比較例１－４のフィッティング処理した結果を示す図である。実施例２で用いたクロマトグラムである。実施例２－１のフィッティング処理した結果を示す図である。比較例２－１のフィッティング処理した結果を示す図である。比較例２－２のフィッティング処理した結果を示す図である。比較例２－３のフィッティング処理した結果を示す図である。比較例２－４のフィッティング処理した結果を示す図である。

以下、本発明の実施の形態について、詳細に説明する。なお、説明の理解を容易にするため、各図面において同一の構成要素に対しては同一の符号を付して、重複する説明は省略する。また、図面における各部材の縮尺は実際とは異なる場合がある。本明細書において数値範囲を示すチルダ「～」は、別段の断わりがない限り、その前後に記載された数値を下限値及び上限値として含むことを意味する。

＜クロマトグラムデータ処理装置＞
一実施形態に係るクロマトグラムデータ処理装置について説明する。図１は、一実施形態に係るクロマトグラムデータ処理装置の機能を示すブロック図である。図１に示すように、クロマトグラムデータ処理装置１は、データ処理部１０と、入力部２０と、出力部３０と、表示部４０とを備える。クロマトグラムデータ処理装置１は、データ処理部１０で測定対象試料に含まれる成分を、クロマトグラフである分離カラムで分離して、分離カラムから溶出した測定対象試料中の分離された成分を検出して得られるクロマトグラムを処理することにより成分を測定する。

測定対象試料(以下、単に試料という)は、液体、気体のいずれでもよい。試料として、例えば、生体由来の物質（血液、汗、唾液、又は尿等）、医薬品、食品添加物、合成された化学物質（農料等）、又は環境負荷物質（工場等から排出される排水、廃液又は地下水等）が挙げられる。本実施形態では、クロマトグラムデータ処理装置１で試料中の成分の分析を行うものとする。

データ処理部１０は、クロマトグラム取得部１１と、微分データ作成部１２と、ピーク解析部１３と、ピーク面積比較部１４と、フィッティング部１５と、判定部１６とを備える。

クロマトグラム取得部１１は、複数の成分に由来するピークを有するクロマトグラムを取得する。図２に示すように、クロマトグラム取得部１１は、溶出時間（保持時間）及び信号強度（電圧）を、それぞれＸ軸及びＹ軸とした時系列データを取得する。取得したクロマトグラムには、試料中の各種成分に由来する複数のピークが現れる。図２に示すクロマトグラムは、目的成分に由来するピーク（以降、単に目的ピークという）の頂点Ｐ１とその周辺部を拡大表示したクロマトグラムである。目的ピークは、頂点Ｐ１と、２つの谷Ｔ１及びＴ２と、２つの変曲点Ｉ１及びＩ２とを有する。谷Ｔ１は、図２においてＰ１の左側に位置する谷であり、谷Ｔ２は、図２においてＰ１の右側に位置する谷である。変曲点Ｉ１は頂点Ｐ１と谷Ｔ１との間に位置し、変曲点Ｉ２は頂点Ｐ１と谷Ｔ２との間に位置する。

クロマトグラム上に出現した、あるピークの信号強度は、あるピークと対応する成分の、試料中における成分濃度に比例する。試料の種類や分析条件等によっては、ある溶出時間の周りに複数の成分が近接して溶出するため、複数のピークが混在した不分離ピークが発生する場合がある。この場合、各ピーク同士が重なり合った状態となっているため、不分離ピークは複数のピークに分解（分割）させて、ピークと対応する成分の検出性能を向上させる。

微分データ作成部１２は、クロマトグラム取得部１１で得たクロマトグラムを平滑化処理して、平滑化処理したクロマトグラムを少なくとも１回微分することで、微分データを取得する。微分データ作成部１２は、クロマトグラムを平滑化処理して１回微分することで、図３に示すような１次微分データを取得することができる。なお、図３の縦軸は信号強度（単位：ｕＶ）を溶出時間（単位：分）で微分した微分値（ｄ（ｕＶ）／ｄｔ）を示す。

図３に示す１次微分データは、クロマトグラムを１回微分して得られたデータである。図３の正の領域は、図２のクロマトグラムの傾きが増大している領域で、図３の負の領域は、図２のクロマトグラムの傾きが減少している領域である。例えば、図２に示す、頂点Ｐ１、谷Ｔ１及びＴ２は、図３の１次微分データで、１次微分値がゼロとなる点である。図３の１次微分データにおいて、１次微分値がゼロとなる点を探索し、その前後において、１次微分値が正から負へ、又は負から正へと推移しているかを判定することで、これらのゼロとなる点が頂点Ｐ１、又は谷Ｔ１及びＴ２であるかを判定することができる。また、１次微分値がゼロとなる点の溶出時間を求めることで、頂点Ｐ１、谷Ｔ１及びＴ２の溶出時間を正確に求めることができる。２つの谷Ｔ１及びＴ２が正確に求められるので、谷Ｔ１と谷Ｔ２との間をピーク面積の境界とした時、ピーク面積が正確に求められる。

また、図２に示すクロマトグラムの変曲点Ｉ１が、図３では検出強度が上昇から下降に変わる点となり、図２に示すクロマトグラムの変曲点Ｉ２が、図３では検出強度が下降から上昇に変わる点となる。変曲点Ｉ１及びＩ２が正確に求められるので、変曲点Ｉ１と変曲点Ｉ２との間をピーク幅とした時、ピーク幅が正確に求められる。

微分データ作成部１２は、クロマトグラムの平滑化処理及び微分処理に、クロマトグラムの平滑化と微分とを同時に行うＳａｖｉｔｚｋｙ－Ｇｏｌａｙ法（ＳＧ法）等を用いることができる。Ｓａｖｉｔｚｋｙ－Ｇｏｌａｙ法は、最小二乗平滑化処理の一例であり、他の方式と比較して比較的容易に複雑な処理を実行することが可能である。なお、クロマトグラムの平滑化及び微分手法としては、Ｓａｖｉｔｚｋｙ－Ｇｏｌａｙ法に限らず、適宜各種手法を使用することができる。

ピーク解析部１３は、図４に示すように、微分データ作成部１２で得られた１次微分データからクロマトグラムの目的成分のピーク面積Ｓ１を計算する。具体的には、ピーク解析部１３は、微分データ作成部１２で得られた１次微分データから目的ピークの谷Ｔ１及びＴ２を導出する。その後、谷Ｔ１から谷Ｔ２までの範囲において、図４に示すように、図２のクロマトグラムにおいて頂点Ｐ１を有する目的ピークを積分することで得られる面積を、ピーク面積Ｓ１として算出する。

ピーク面積比較部１４は、複数のピークが出現する溶出時間領域（ピーク検出領域）において、複数のピークのピーク面積比を比較することで、目的ピークであるかを判断する。具体的には、ピーク面積比較部１４は、予め、目的ピークが出現する溶出時間領域と、その領域内に現れるピーク面積比の範囲とを規定する。そして、ピーク面積比較部１４は、図５に示すようなクロマトグラムでは、ピーク検出領域Ａ１の中にあるすべてのピークに対して、溶出時間及び面積比の何れもが、規定された範囲内にあるか判定を行う。

前記条件を満たすピークが存在した場合は、前記条件を満たすピークが目的ピークとなる。図５では、２つのピークのうち、溶出時間が遅いピークをピーク１とし、溶出時間が早いピークをピーク２とする。ピーク面積比較部１４は、２つのピーク（ピーク１、ピーク２）の溶出時間及び面積の何れもが、所定の範囲内にあるかを判断することで、目的ピークであるかを検出する。図５のピーク１では溶出時間及び面積比の範囲の何れもが所定の範囲内にあり、図５のピーク２では溶出時間のみ所定の範囲内にある。よって、ピーク面積比較部１４は、ピーク１を目的ピークと判断することができ、ピーク１と隣接するピーク２を目的ピークから除外することができる。

フィッティング部１５は、図２に示すクロマトグラム内の複数のピークに対して、下記式（I）に示す放物型分散修正ガウシアン(ＰＶＭＧ：Parabolic-Variance Modified Gaussian)関数又はＥＭＧ関数と、ピアソン（Pearson）関数とを用いて、フィッティング処理を行い、クロマトグラム内の複数のピークを分解する。

（式（I）中、ｘは溶出時間であり、μはあるピーク（頂点）の溶出時間であり、σは標準偏差であり、Ａ及びａ及びｂは任意定数である。）

フィッティング部１５は、フィッティング処理を行うことで、例えば、図２に示すクロマトグラムとほぼ同じ溶出時間及びピーク高さを有しつつ、隣接するピーク同士が完全に分解されたデータを、近似的に得ることができる。目的ピークを検出し、検出された目的ピークから目的ピークに由来する成分濃度を算出するには、ピーク面積比を正確に算出できることが重要である。フィッティング部１５では、クロマトグラム内に分解不十分なピークが複数存在したとしても、分解不十分な複数のピークから目的ピークのみをフィッティング処理により分解することができる。そのため、目的ピークに由来する成分の成分濃度や検出精度等が高まる。

ピアソン（Pearson）関数としては、下記式（II）に示すピアソンVII（Pearson VII）関数を用いることができる。

（式（II）中、ｘは溶出時間であり、μはあるピーク（頂点）の溶出時間であり、σは標準偏差であり、ξは形状パラメータであり、Ａは任意定数である。）

フィッティング部１５は、クロマトグラム取得部１１で得られたクロマトグラムの複数のピークに対して、下記式（１）を用いることが好ましい。

（式（１）中、ｘは溶出時間であり、μはあるピーク（頂点）の溶出時間であり、σは標準偏差であり、Ａ及びａ及びｂは任意定数であり、ξは形状パラメータである。）

フィッティング部１５は、フィッティング処理により分解された、複数のピーク中に含まれる目的ピークと、目的ピークと最隣接するピークとの分離度が１．０以上とすることが好ましい。ピークの分離度が１．０未満では、図６に示すように、隣接する２つのピークの重なりが大きく、これらのピーク間の分解状態が不十分となる。結果、それぞれのピーク同士が影響し合い、ピークのリーディング又はテーリングが発生してしまう。即ち、フィッティング処理により、目的ピークを分解しても、目的ピークに由来する成分濃度を正確に算出できなくなり、検出精度が低下する。一方、ピークの分離度が１．０以上では、図７に示すように、隣接する２つのピークの重なりが小さくなるため、複数ピークの中から、フィッティング処理により目的ピークを、容易に分解することができる。即ち、目的ピークに由来する成分濃度を正確に算出することができ、検出精度を高められる。なお、ピークの分離度が１．５以上であれば、隣接する２つのピークは完全に分解したと見なせるため、ピークの分離度は１．５以上であることが望ましい。

判定部１６は、フィッティング部１５でフィッティング処理を行った後のクロマトグラムから成分を判定する。

入力部２０は、クロマトグラムのデータが入力されると共に、データ処理に必要なパラメータ等が入力される。入力部２０は、例えば、キーボード、マウス、操作ボタン、タッチパネル等が挙げられる。

出力部３０は、判定部１６で判定された成分を出力する。出力部３０としては、例えば、プリンタ等が挙げられる。

表示部４０は、判定部１６で判定された成分やクロマトグラムのピークの分解結果等を表示する。表示部４０としては、例えば、モニタディスプレイ等が挙げられる。

＜クロマトグラムデータ処理装置のハードウェア構成＞
次に、クロマトグラムデータ処理装置のハードウェア構成の一例について説明する。図８は、クロマトグラムデータ処理装置１のハードウェア構成を示すブロック図である。図８に示すように、クロマトグラムデータ処理装置１は、情報処理装置（コンピュータ）で構成され、物理的には、演算処理部であるＣＰＵ（ＣｅｎｔｒａｌＰｒｏｃｅｓｓｉｎｇＵｎｉｔ：プロセッサ）１０１、主記憶装置であるＲＡＭ（ＲａｎｄｏｍＡｃｃｅｓｓＭｅｍｏｒｙ）１０２及びＲＯＭ（ＲｅａｄＯｎｌｙＭｅｍｏｒｙ）１０３、入力デバイスである入力装置１０４、出力装置１０５、通信モジュール１０６並びにハードディスク等の補助記憶装置１０７等を含むコンピュータシステムとして構成することができる。これらは、バス１０８で相互に接続されている。なお、出力装置１０５及び補助記憶装置１０７は、外部に設けられていてもよい。

ＣＰＵ１０１は、クロマトグラムデータ処理装置１の全体の動作を制御し、各種の情報処理を行う。ＣＰＵ１０１は、ＲＯＭ１０３又は補助記憶装置１０７に格納されたクロマトグラムデータ処理プログラムを実行して、測定収録画面と解析画面の表示動作を制御する。

ＲＡＭ１０２は、ＣＰＵ１０１のワークエリアとして用いられ、主要な制御パラメータや情報を記憶する不揮発ＲＡＭを含んでもよい。

ＲＯＭ１０３は、基本入出力プログラム等を記憶する。クロマトグラムデータ処理プログラムはＲＯＭ１０３に保存されてもよい。

入力装置１０４は、キーボード、マウス、操作ボタン、タッチパネル等である。

出力装置１０５は、モニタディスプレイ等である。出力装置１０５では、予測結果等が表示され、入力装置１０４や通信モジュール１０６を介した入出力操作に応じて画面が更新される。

通信モジュール１０６は、ネットワークカード等のデータ送受信デバイスであり、外部のデータ収録サーバ等からの情報を取り込み、他の電子機器に解析情報を出力する通信インタフェースとして機能する。

補助記憶装置１０７は、ＳＳＤ（ＳｏｌｉｄＳｔａｔｅＤｒｉｖｅ）、及びＨＤＤ（ＨａｒｄＤｉｓｋＤｒｉｖｅ）等の記憶装置であり、例えば、クロマトグラムデータ処理プログラムやクロマトグラムデータ処理装置１の動作に必要な各種のデータ、ファイル等を格納する。

図１に示すクロマトグラムデータ処理装置１の各機能は、ＣＰＵ１０１、ＲＡＭ１０２等の主記憶装置又は補助記憶装置１０７に所定のコンピュータソフトウェア（クロマトグラムデータ処理プログラムを含む）を読み込ませ、ＲＡＭ１０２、ＲＯＭ１０３又は補助記憶装置１０７に格納された予測プログラム等をＣＰＵ１０１により実行する。入力装置１０４、出力装置１０５及び通信モジュール１０６を動作させると共に、ＲＡＭ１０２、ＲＯＭ１０３及び補助記憶装置１０７等におけるデータの読み出し及び書き込みを行うことで、クロマトグラムデータ処理装置１の各機能は、実現される。すなわち、本実施形態のクロマトグラムデータ処理プログラムをコンピュータ上で実行させることで、クロマトグラムデータ処理装置１は、図１の、クロマトグラム取得部１１と、微分データ作成部１２と、ピーク解析部１３と、ピーク面積比較部１４と、フィッティング部１５と、判定部１６と、入力部２０と、出力部３０と、表示部４０として、それぞれ機能することができる。

クロマトグラムデータ処理プログラムは、以下の構成のプログラムを用いることができる。

すなわち、クロマトグラムデータ処理プログラムは、
測定対象試料に含まれる複数の成分が分離カラムで分離され、分離された前記成分を検出して得られるクロマトグラムにより前記成分を測定する際にコンピュータに実行させるであって、
時間と信号強度とをそれぞれ軸として、複数の前記成分に由来するピークを有する前記クロマトグラムを得るクロマトグラム取得部と、
前記クロマトグラム内の複数の前記ピークに対して、ガウシアン関数又は指数関数で修飾されたガウシアン関数と、ピアソン関数とを用いて、フィッティング処理を行い、複数の前記ピークを分解するフィッティング部と、
前記フィッティング処理を行った後の前記クロマトグラムから前記成分を判定する判定部と、
をコンピュータに実行させる。

クロマトグラムデータ処理プログラムは、例えば、ＲＡＭ２２やＲＯＭ２３の主記憶装置又は補助記憶装置２４等のコンピュータが備える記憶装置内に格納することができる。また、クロマトグラムデータ処理プログラムは、その一部又は全部が、通信回線等の伝送媒体を介して伝送され、コンピュータが備える通信モジュール等により受信されて記録（インストールを含む）される構成としてもよい。

また、クロマトグラムデータ処理プログラムは、その一部又は全部が、ＣＤ－ＲＯＭ、ＤＶＤ－ＲＯＭ、フラッシュメモリ等の携帯可能な記憶媒体に格納し、コンピュータで読み取り可能としてもよい。これにより、クロマトグラムデータ処理プログラムは記憶媒体に格納された状態から、コンピュータで読み取らせてコンピュータ内に記録（インストールを含む）させることができる。また、クロマトグラムデータ処理プログラムは記憶媒体に格納された状態で移動させることができる。

＜クロマトグラムデータ処理方法＞
次に、本実施形態に係るクロマトグラムデータ処理装置を用いて、本実施形態に係るクロマトグラムデータ処理方法について説明する。本実施形態に係るクロマトグラムデータ処理方法は、図１に示すような構成を有するクロマトグラムデータ処理装置１において、測定対象試料に含まれる成分を、予め分離カラムで分離して、分離カラムから溶出した測定対象試料中の分離された成分を検出して得られるクロマトグラムにより成分を測定する。

図９は、本実施形態に係るクロマトグラムデータ処理方法を説明するフローチャートである。図９に示すように、本実施形態に係るクロマトグラムデータ処理方法は、クロマトグラム取得工程（ステップＳ１１）、微分データ作成工程（ステップＳ１２）、ピーク解析工程（ステップＳ１３）、ピーク面積比較工程（ステップＳ１４）、フィッティング工程（ステップＳ１５）及び判定工程（ステップＳ１６）を含む。

クロマトグラムデータ処理装置１は、クロマトグラム取得部１１により、図２に示すような溶出時間及び信号強度をそれぞれＸ軸及びＹ軸とした時系列データ、即ち、試料中の各種成分に由来する複数のピークが出現するクロマトグラム（生データ）を入力情報として得る（クロマトグラム取得工程：ステップＳ１１）。

取得したクロマトグラムには、試料中の各種成分に由来する複数のピークが現れる。図２に示すように、クロマトグラムには、頂点Ｐ１と、２つの谷Ｔ１及びＴ２と、２つの変曲点Ｉ１及びＩ２とを有する目的ピークが出現する。

次に、クロマトグラムデータ処理装置１は、微分データ作成部１２により、クロマトグラム取得部１１で得たクロマトグラムを平滑化処理して１回微分することで、図３に示すような１次微分データを得る（微分データ作成工程：ステップＳ１２）。

図２に示す、目的ピークの頂点Ｐ１、谷Ｔ１及びＴ２は、図３の微分データでは、１次微分値がゼロとなる点であるため、１次微分値がゼロとなる点の溶出時間を求めることで、頂点Ｐ１、谷Ｔ１及びＴ２の溶出時間を正確に求めることができる。また、図３の微分データより、変曲点Ｉ１及びＩ２が正確に求められるので、変曲点Ｉ１と変曲点Ｉ２との間をピーク幅とした時、ピーク幅も正確に求められる。

次に、クロマトグラムデータ処理装置１は、ピーク解析部１３により、図４に示すように、微分データ作成工程（ステップＳ１２）で得られた１次微分データから、目的ピークの谷Ｔ１及びＴ２を導出する。その後、谷Ｔ１から谷Ｔ２までの間で目的ピークを積分計算することで、目的ピークのピーク面積Ｓ１を算出する（ピーク解析工程：ステップＳ１３）。

次に、クロマトグラムデータ処理装置１は、ピーク面積比較部１４により、複数のピークを含む溶出時間領域において、複数のピークのピーク面積比を比較して、目的ピークを判断する（ピーク面積比較工程：ステップＳ１４）。

具体的には、予め、目的ピークが出現するピーク検出領域と、その領域内に現れるピークの面積比の範囲とを規定する。そして、図５に示すようなクロマトグラムでは、予め規定した溶出時間のピーク検出領域Ａ１の中で、溶出時間及び面積比の何れもが所定の範囲内にあるピークが存在するかを判定する。図５に示す２つのピーク（ピーク１、ピーク２）のうち、溶出時間が遅く且つ面積比が範囲内にあるピークをピーク１、溶出時間が早く且つ面積比が範囲外にあるピークをピーク２とする。このとき、クロマトグラムデータ処理装置１は、ピーク１を目的ピークと判断することができ、ピーク１と隣接するピーク２を目的ピークから除外することができる。

次に、クロマトグラムデータ処理装置１は、フィッティング部１５により、クロマトグラム内の複数のピークに対して、上記式（I）に示すＰＶＭＧ関数又はＥＭＧ関数と、ピアソン関数とを用いて、フィッティング処理を行い、クロマトグラム内の複数のピークを分解する（フィッティング工程：ステップＳ１５）。

フィッティング処理を行うことで、例えば、図２に示すクロマトグラムとほぼ同じ溶出時間及びピーク高さを有しつつ、隣接するピーク同士が完全に分解されたデータを、近似的に得ることができる。フィッティング処理により、クロマトグラム内の分解不十分な複数ピークから、目的ピークのみを分解することができるので、目的ピーク由来成分の成分濃度や検出精度等を高めることができる。

ピアソン（Pearson）関数としては、上記式（II）に示すピアソン VII（Pearson VII）関数を用いることが好ましい。

クロマトグラムデータ処理装置１は、クロマトグラムの複数のピークに対して、上記式（１）を用いて、フィッティング処理を行うことが好ましい。

次に、クロマトグラムデータ処理装置１は、判定部１６により、フィッティング工程(ステップＳ１５）でフィッティング処理を行った後のクロマトグラムから、試料中の目的成分を判定する（判定工程：ステップＳ１６）。

クロマトグラムデータ処理装置１は、判定部１６による判定結果を、出力部３０により出力してもよいし、表示部４０に表示してもよい。

以上の通り、クロマトグラムデータ処理装置１は、クロマトグラム取得部１１、フィッティング部１５及び判定部１６を備える。クロマトグラムデータ処理装置１は、クロマトグラム取得部１１で、試料中の各種成分に由来する複数のピークを有するクロマトグラムを得る。そして、クロマトグラムデータ処理装置１は、フィッティング部１５で、クロマトグラム内の複数のピークに対して、ＰＶＭＧ又はＥＭＧ関数と、ピアソン関数とを用いてフィッティング処理を行い、複数のピークから目的ピークを分解する。クロマトグラムデータ処理装置１は、フィッティング部１５でフィッティング処理したクロマトグラムから、判定部１６で分解した目的ピークに由来する成分を判定する。

分離カラムで試料中の各種成分を分離して得られるクロマトグラムには、ピーク形状が急峻で左右対称なピーク、急峻で裾野の長い左右非対称なピーク、起伏に乏しくブロードなピーク等、様々な形状のピークが現れる。ピーク形状は、主に、分離カラムの大きさや材質（カラムの基材や修飾基）等で決定される分離性能によって、大きく左右され易い。例えば、試料体積が極めて微量で、且つ、一般的な粒子カラムを用いて成分を分離する場合、カラム体積を可能な限り抑えないと、クロマトグラムのピーク形状がブロードになってしまう。このような場合、ＥＭＧ関数、又はガウシアン関数をローレンツ関数で重み付けした関数等を用い、クロマトグラムのフィッティング処理を行うことで、目的ピークに由来する成分を分解することができる。なお、ローレンツ関数とは、下記式（III）のように示される。

（式（III）中、ｘは溶出時間であり、μはあるピーク（頂点）の溶出時間であり、σは標準偏差であり、Ａは任意定数である。）

一方、無機材料や有機材料等からなるモノリス構造体で構成されるカラム（以下、単にモノリスカラムという）の場合、モノリスカラムは、モノリス構造体の骨格構造（マクロ孔）と、骨格構造中に存在する微小な空隙（メソ孔）によって、粒子カラムと同じカラム体積であっても、粒子カラムよりも低い送液圧力で、高い分離性能を発揮することができる。よって、試料体積が極めて微量で、且つ、モノリスカラムを用いて成分の分離を行う場合、精度良く分離することが可能となる。モノリスカラムを用いて得られるクロマトグラムのピーク形状は、一般的な粒子カラムと比較して、ピーク頂点近傍が急峻で且つピーク裾野が広がり易い傾向にあるため、指数関数により修飾されたガウシアン関数（ＥＭＧ関数）等では、精度良くフィッティングすることが困難となる。

クロマトグラムデータ処理装置１は、フィッティング部１５でフィッティング処理を行う際に、ガウシアン関数又はＥＭＧ関数と、ピアソン関数を合成した関数により、モノリスカラムを用いて得られたクロマトグラムのピークを、精度良くフィッティングすることができる。この結果、クロマトグラムデータ処理装置１は、判定部１６において、フィッティング処理後の曲線データ（分解後の目的ピーク）から、試料中に含まれる目的成分を高精度に判定できる。

よって、クロマトグラムデータ処理装置１は、試料が少量であって、試料中の各種成分が十分に分離されていない場合であっても、クロマトグラム上の各種成分由来のピークをフィッティングし分解することができるので、クロマトグラム上の各種成分濃度や検出精度等を高精度で算出することができる。そのため、クロマトグラムデータ処理装置１は、特に、小型で送液圧力が低く、高い分離性能を有する分離カラム（例えば、モノリスカラム）を用いた場合に有効に用いることができる。

クロマトグラムデータ処理装置１は、フィッティング部１５で、フィッティング処理を行う際、クロマトグラムの複数のピークに対して、上記式（１）を用いることができる。上記式（１）は、ＰＶＭＧ関数とピアソンVII関数を合成（乗算）した関数である。そのため、クロマトグラム上のピークが、モノリスカラム等の高い分離性能を有する分離カラムを用いて分離した場合に得られるような、ピーク頂点近傍が急峻で且つピーク裾野が広がり易い場合でも、フィッティング処理を精度よく行うことができる。この結果、クロマトグラムデータ処理装置１は、判定部１６において、フィッティング処理後の曲線データから試料中に含まれる目的成分を高精度に判定することができる。

クロマトグラムデータ処理装置１は、微分データ作成部１２で、クロマトグラムの平滑化処理を行い、平滑化処理したクロマトグラムを少なくとも１回微分することで、微分データを得ることができる。クロマトグラムの平滑化及び微分を行うことで、得られる微分データのノイズ成分を低減することができる。特に、１次微分値がゼロとなる点を探索する場合、微分データ作成部１２がクロマトグラムの平滑化処理を行うことによって、頂点Ｐ１、谷Ｔ１及びＴ２の溶出時間をより正確に求めることができる。谷Ｔ１及びＴ２の溶出時間をより正確に求めることで、各ピークのピーク面積をより正確に求めることができる。

クロマトグラムデータ処理装置１は、ピーク解析部１３で、微分データ作成部１２で得られた１次微分データから目的ピークの積分範囲を求めて目的ピークのピーク面積を計行し、ピーク面積比較部１４で、複数のピークを含む溶出時間領域において、複数のピークのピーク面積を比較して、目的ピークを判定することができる。これにより、クロマトグラムデータ処理装置１は、目的ピークをより正確に判定できるので、試料中に含まれる目的成分をより高精度に判定することができる。

このように、クロマトグラムデータ処理装置１は、分離する試料が少量であっても、フィッティング処理を行うことで、クロマトグラムの各種成分を分解できるため、試料中の各種成分濃度や検出精度等を高精度で算出することができる。特に、小型で送液圧力が低く、高い分離性能を有する分離カラム（例えば、モノリスカラム）を通過させた試料に対して、有効に用いることができる。そのため、クロマトグラムデータ処理装置１は、例えば、血液中に含まれるタンパク質や核酸等の血液成分、工場等から排出される排水中に含まれる化学物質、地下水に含まれる微量な物質等の分析に有効に用いることができる。特に、血液中のタンパク質濃度の分析に有効に用いることができる。そのため、流路プレート１０Ａは、臨床検査、食物検査、環境検査、又は診療や看護等の医療現場等、様々な用途において好適に用いることができる。

本実施形態に係るクロマトグラムデータ処理装置１を用いて分離カラムを含む流路プレートから排出される試料を分析する分析装置について説明する。なお、ここでは、試料として血液を用い、血液成分として、赤血球中に含まれるタンパク質の一種であるヘモグロビンに対しグルコースが結合したグリコヘモグロビン（ＨｂＡ１ｃ）を測定する場合について説明する。

図１０は、クロマトグラムデータ処理装置１を適用した分析装置の一例を示す図である。図１０に示すように、分析装置６０は、クロマトグラムデータ処理装置１と、流路プレート６１と、分離カラム６２と、容器６３と、送液ポンプ６４と、流路切替弁（切替バルブ）６５と、供給管６６Ａと、排出管６６Ｂと、検出器６７とを備える。

流路プレート６１は、２枚のプレートで構成されており、各プレート同士が接する面上に、液体が通過する流路部及び分離カラム収容部を有している。流路プレート６１は、前記流路及び収容部が正対する状態で、板厚方向に前記プレートを積層することで構成されている。流路プレート６１には、例えば、オレフィン系樹脂、アクリル系樹脂、スチレン系樹脂、ビニル系樹脂、フッ素系樹脂、エンジニアリングプラスチック、スーパーエンジニアリングプラスチック、熱硬化性樹脂、ガラス等を用いることができる。流路プレート６１は、その内部に、液体が通過する流路を有して、主面に流入口６１ａ及び流出口６１ｂを有する。流路プレート６１は、流路の途中に分離カラム６２を収容している。

分離カラム６２は、流路プレート６１内の流路の途中に配置され、２枚のプレートの間に封止された状態で収容されている。分離カラム６２は、血液中の成分を分離するものであり、例えば、液体クロマトグラフィー用の分離カラムが用いられる。

分離カラム６２は、円柱状に加工された多孔質体や微粒子の集合体で構成されている。分離カラム６２は、通過する液体に含まれる各種成分と分離カラム６２との相互作用（例えば、疎水性相互作用、イオン交換作用等）によって、血液中に含まれる各種成分を分離させる。具体的には、相互作用によって生じる、成分同士の移動速度の差や分子の大きさ及び荷電状態等を利用して、血液中に含まれる各種成分を分離する。

固定相の材料は、無機材料や有機材料等から選択される。固定相は、無機材料や有機材料等からなるモノリス構造体を含むことができる。モノリス構造体は、骨格構造（マクロ孔）と、骨格構造中に存在する微小な空隙（メソ孔）の大きさを、目的に合わせて適宜設計することができる。モノリス構造体としては、シリカモノリス等が用いられる。

容器６３は、外部から注入された、移動相を一時的に収容する容器６３Ａと、試料である血液を一時的に収容する容器６３Ｂとを有している。

送液ポンプ６４は、容器６３Ａ内の移動相を送液する送液ポンプ６４Ａと、容器６３Ｂ内の血液を送液する送液ポンプ６４Ｂとを有する。

切替バルブ６５は、送液ポンプ６４Ａ又は送液ポンプ６４Ｂと流路プレート６１とを連結する流路を切り替える。切替バルブ６５で流路を切り替えることで、移動相又は血液が流路プレート６１に送液される。

供給管６６Ａ及び排出管６６Ｂは、流路プレート６１の主面に設けた流入口６１ａ及び流出口６１ｂに挿脱可能に配置されている。供給管６６Ａは、流入口６１ａに挿入され、血液を流路プレート６１内の流路に供給する。排出管６６Ｂは、流出口６１ｂに挿入され、血液を流路プレート６１内の流路から排出する。

検出器６７は、血液中の分離された成分を検出する。検出器６７としては、ＵＶ検出器、三次元検出器、示差屈折率検出器、質量検出器等が用いられる。血液中の分離された成分を、検出器６７で検出することにより得られるクロマトグラムは、クロマトグラムデータ処理装置１にクロマトグラムデータとして入力される。

分析装置６０では、容器６３中に注入された血液を送液ポンプ６４が吸引して一定流量で供給管６６Ａに供給する。血液は、供給管６６Ａを介して流入口６１ａから流路プレート６１内の流路に送られる。血液は、流路内を通って、分離カラム６２に送られ、分離カラム６２を通過することで、血液中の各種成分が分離されて、分離カラム６２の出口から溶出する。分離カラム６２を通過した血液は、流路内を通り、流出口６１ｂから排出管６６Ｂに流れて外部に排出される。排出管６６Ｂから外部に排出された血液は、検出器６７で検出することにより、溶出時間に対する信号強度（クロマトグラム）として変換される。測定に用いられた血液は、不図示の回収容器等で回収されてもよいし、そのまま廃棄してもよい。検出器６７で得られた血液のクロマトグラムは、クロマトグラムデータとしてクロマトグラムデータ処理装置１のデータ処理部１０に入力される。

クロマトグラムデータ処理装置１は、データ処理部１０でクロマトグラムデータを処理して、血液中に含まれる各種タンパク質の中から、特定のタンパク質であるＨｂＡ１ｃを検出する。クロマトグラムデータ処理装置１は、データ処理部１０で処理されたデータを出力部３０で出力すると共に、表示部４０に測定結果を表示する。

このように、分析装置６０は、血液が少量であっても、クロマトグラムデータ処理装置１で、クロマトグラム上の各種タンパク質由来のピークをフィッティング処理により分解できるので、特定のタンパク質であるＨｂＡ１ｃを高精度に検出することができる。そのため、分析装置６０は、分離カラム６２として、小型で送液圧力が低く高い分離性能を有するシリカモノリス等を用いて、血液中のタンパク質を分離することで、少量の血液で、血液中のＨｂＡ１ｃを高精度に検出することができる。

そのため、分析装置６０は、例えば、糖尿病の診断基準として用いられるＨｂＡ１ｃ値を高精度に算出することができるため、糖尿病の診断等に好適に用いることができる。

また、分析装置６０は、分離カラム６２が小型であり、血液が通る流路を小さくできるので、送液ポンプ６４の送液圧力を低くすることができる。

さらに、分析装置６０は、分離カラム６２を流路プレート６１内に収容することができるので、分析装置６０の小型化を図ることができると共に、取り扱いを簡単にすることができる。

なお、本実施形態では、試料が液体である場合について説明したが、試料は気体でもよい。

本実施形態においては、データ処理部１０は、ピーク解析部１３及びピーク面積比較部１４を備えなくてもよいし、微分データ作成部１２、ピーク解析部１３及びピーク面積比較部１４を備えなくてもよい。

本実施形態においては、微分データ作成部１２は、クロマトグラム取得部１１で得たクロマトグラムを２回以上微分した微分データを取得してもよい。

また、本実施形態においては、微分データ作成部１２は、クロマトグラムの平滑化処理が不要な場合等には、クロマトグラムの平滑化処理を行わず、微分のみを行うようにしてもよい。

本実施形態においては、微分データ作成部１２がクロマトグラム取得部１１で得たクロマトグラムの平滑化及び微分を行っているが、データ処理部１０は、クロマトグラムの平滑化処理を行う平滑化処理部と、クロマトグラムの微分のみを行う微分データ作成部１２とを別々に備えてもよい。

以下、実施例及び比較例を示して実施形態を更に具体的に説明するが、実施形態はこれらの実施例及び比較例により限定されるものではない。

＜実施例１＞
［実施例１－１］
（血液中の成分の分析）
図１０に示す分析装置６０を準備し、測定対象液体（試料）として血液を用いた。血液を流路プレート６１の流入口６１ａから供給して、分離カラム６２で血液中の成分を分離させた後、流路プレート６１の流出口６１ｂから血液を排出させた。排出された血液を検出器６７に供給し、血液中の各種成分を検出することで、図１１に示すような血液成分のクロマトグラムを得た。図１１中の破線領域内に、赤血球中のグリコヘモグロビン（ＨｂＡ１ｃ）に由来するピークが２つ検出された。図１１に示すクロマトグラムの破線領域内では、隣接する２つのピークの分離度が低かった。図１１中の破線領域内の２つのピークのうち、溶出時間が約１．２分程度で検出されたピークが不安定型Ａ１ｃ（ＬａｂｉｌｅＡ１ｃ：Ｌ－Ａ１ｃ）に由来するピークであり、溶出時間が約１．６分程度で検出されたピークが安定型Ａ１ｃ（ＳｔａｂｌｅＡ１ｃ:Ｓｔ－Ａ１ｃ）に由来するピークである。

（フィッティング処理）
図１１に示すクロマトグラムの破線領域を、上記式（１）を用いてフィッティング処理した。フィッティング処理した結果を図１２に示す。図１２中、細線が、図１１に示すクロマトグラムのデータ（生データ）を示し、太線が、フィッティング処理後のデータを示し、破線が、フィッティング処理により分解した成分のピークを示す。

［比較例１－１］
実施例１－１において、フィッティング処理を上記式（１）に代えてＥＭＧ関数を用いたこと以外は、実施例１－１と同様にして行った。フィッティング処理した結果を図１３に示す。図１３中の細線、太線及び破線は、図１２と同様である。

［比較例１－２］
実施例１－１において、フィッティング処理を上記式（１）に代えてＰＶＭＧ関数を用いたこと以外は、実施例１－１と同様にして行った。ＰＶＭＧ関数として下記式（I）を用いた。フィッティング処理した結果を図１４に示す。図１４中の細線、太線及び破線は、図１２と同様である。

［比較例１－３］
実施例１－１において、フィッティング処理を上記式（１）に代えてローレンツ関数を用いたこと以外は、実施例１－１と同様にして行った。ローレンツ関数とは、下記式（III）を用いた。フィッティング処理した結果を図１５に示す。図１５中の細線、太線及び破線は、図１２と同様である。

［比較例１－４］
実施例１－１において、フィッティング処理を上記式（１）に代えてピアソンVII関数を用いたこと以外は、実施例１－１と同様にして行った。ピアソンVII関数として下記式（II）を用いた。フィッティング処理した結果を図１６に示す。図１６中の細線、太線及び破線は、図１２と同様である。

＜実施例２＞
［実施例２－１］
（クロマトグラムの準備）
実施例１－１で得た図１１に示すクロマトグラムに代えて、図１７に示すような、隣接する２つのピークの分離度が高いクロマトグラムを用いたこと以外は、実施例１－１と同様にして行った。図１７中の破線領域内の２つのピークのうち、溶出時間が約０．９分程度で検出されたピークがＬ－Ａ１ｃに由来するピークであり、溶出時間が約１．３分程度で検出されたピークがＳｔ－Ａ１ｃに由来するピークである。フィッティング処理した結果を図１８に示す。図１８中の細線、太線及び破線は、図１２と同様である。

［比較例２－１］
実施例２－１において、フィッティング処理を上記式（１）に代えてＥＭＧ関数を用いたこと以外は、実施例２－１と同様にして行った。フィッティング処理した結果を図１９に示す。図１９中の細線、太線及び破線は、図１２と同様である。

［比較例２－２］
実施例２－１において、フィッティング処理を上記式（１）に代えて上記の比較例１－２で用いたＰＶＭＧ関数を用いたこと以外は、実施例２－１と同様にして行った。フィッティング処理した結果を図２０に示す。図２０中の細線、太線及び破線は、図１２と同様である。

［比較例２－３］
実施例２－１において、フィッティング処理を上記式（１）に代えて上記の比較例１－３で用いたローレンツ関数を用いたこと以外は、実施例２－１と同様にして行った。フィッティング処理した結果を図２１に示す。図２１中の細線、太線及び破線は、図１２と同様である。

［比較例２－４］
実施例２－１において、フィッティング処理を上記式（１）に代えて上記の比較例１－４で用いたピアソンVII関数を用いたこと以外は、実施例２－１と同様にして行った。フィッティング処理した結果を図２２に示す。図２２中の細線、太線及び破線は、図１２と同様である。

図１２に示すように、実施例１―１では、隣接する２つのピークの分離度が低いクロマトグラムにおいて、フィッティング処理後のデータの溶出時間及びピーク高さは、生データのものと略一致した。また、フィッティング処理前後で、隣接する２つのピーク形状に大きな差がなかったことから、フィッティング処理により分解した成分のピークは、近似的には、完全分離されたＬ－Ａ１ｃ及びＳｔ－Ａ１ｃと見なせることが確認された。一方、図１３～図１６に示すように、比較例１―１～１－４では、いずれも、フィッティング処理後の、データの溶出時間、ピーク高さ及びデータの隣接する２つのピーク形状は、生データのものから乖離した。よって、フィッティング処理により分解した成分のピークは、完全分離されたＬ－Ａ１ｃ及びＳｔ－Ａ１ｃと見なすのは困難であることが確認された。

図１８に示すように、実施例２―１では、隣接する２つのピークの分離度が高いクロマトグラムにおいて、フィッティング処理後のデータの溶出時間及びピーク高さは、生データのものとほぼ一致した。また、フィッティング処理前後で、隣接する２つのピーク形状に大きな差がないことから、フィッティング処理により分解した成分のピークは、近似的には、完全分離されたＬ－Ａ１ｃ及びＳｔ－Ａ１ｃと見做せることが確認された。一方、図１９～図２２に示すように、比較例２―１～２－４では、いずれも、フィッティング処理後の、データの溶出時間、ピーク高さ及びデータの隣接する２つのピーク形状は、生データのものから乖離した。よって、フィッティング処理により分解した成分のピークは、完全分離されたＬ－Ａ１ｃ及びＳｔ－Ａ１ｃと見なすのは困難であることが確認された。

よって、本実施形態に係るクロマトグラムデータ処理装置を用いれば、クロマトグラム内のＨｂＡ１ｃに由来する２つのピークを、Ｌ－Ａ１ｃ及びＳｔ－Ａ１ｃに分解することができるため、分解して得られたＳｔ－Ａ１ｃの面積比から、ＨｂＡ１ｃ値を高精度に算出することができる。また、隣接する２つのピークの分離度が低いクロマトグラムであっても、ＨｂＡ１ｃに由来する２つのピークを、Ｌ－Ａ１ｃ及びＳｔ－Ａ１ｃに分解することができる。したがって、本実施形態に係るクロマトグラムデータ処理装置は、夾雑物等の影響を受けたクロマトグラムから、Ｓｔ－Ａ１ｃを分解することができるので、可搬性のある小型な液体クロマトグラフィー装置として有効に用いることができるといえる。

以上の通り、実施形態を説明したが、上記実施形態は、例として提示したものであり、上記実施形態により本発明が限定されるものではない。上記実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の組み合わせ、省略、置き換え、変更等を行うことが可能である。これら実施形態やその変形は、発明の範囲や要旨に含まれると共に、特許請求の範囲に記載された発明とその均等の範囲に含まれる。

本出願は、２０１９年１１月２７日に日本国特許庁に出願した特願２０１９－２１４４５０号に基づく優先権を主張するものであり、特願２０１９－２１４４５０号の全内容を本出願に援用する。

１クロマトグラムデータ処理装置
１０データ処理部
１１クロマトグラム取得部
１２微分データ作成部
１３ピーク解析部
１４ピーク面積比較部
１５フィッティング部
１６判定部
２０入力部
３０出力部
４０表示部

Claims

測定対象試料に含まれる複数の成分が分離カラムで分離され、分離された前記成分を検出して得られるクロマトグラムにより前記成分を測定するクロマトグラムデータ処理装置であって、
時間と信号強度とをそれぞれ軸として、複数の前記成分に由来するピークを有する前記クロマトグラムを得るクロマトグラム取得部と、
前記クロマトグラム内の複数の前記ピークに対して、ガウシアン関数又は指数関数で修飾されたガウシアン関数と、ピアソン関数とを用いて、フィッティング処理を行い、複数の前記ピークを分離するフィッティング部と、
前記フィッティング処理を行った後の前記クロマトグラムから前記成分を判定する判定部と、
を備えることを特徴とするクロマトグラムデータ処理装置。
測定対象試料に含まれる複数の成分が分離カラムで分離され、分離された前記成分を検出して得られるクロマトグラムにより前記成分を測定するクロマトグラムデータ処理方法であって、
時間と信号強度とをそれぞれ軸として、複数の前記成分に由来するピークを有する前記クロマトグラムを得るクロマトグラム取得工程と、
前記クロマトグラム内の複数の前記ピークに対して、ガウシアン関数又は指数関数で修飾されたガウシアン関数と、ピアソン関数とを用いて、フィッティング処理を行い、複数の前記ピークを分解するフィッティング工程と、
前記フィッティング処理を行った後の前記クロマトグラムから前記成分を判定する判定工程と、
を含むことを特徴とするクロマトグラムデータ処理方法。
前記フィッティング工程は、前記クロマトグラムの複数の前記ピークに対して、
下記式（１）:

（式（１）中、ｘは溶出時間であり、μはあるピーク（頂点）の溶出時間であり、σは標準偏差であり、Ａ及びａ及びｂは任意定数であり、ξは形状パラメータである。）
を用いることを特徴とする請求項２に記載のクロマトグラムデータ処理方法。
前記クロマトグラムを少なくとも１回微分して、微分データを得る微分データ作成工程を含むことを特徴とする請求項２又は３に記載のクロマトグラムデータ処理方法。
前記微分データ作成工程で得られた前記微分データから前記クロマトグラムの前記ピークのピーク面積を計算するピーク解析工程と、
複数の前記ピークを含む時間において、複数の前記ピークのピーク面積を比較して、目的とするピークを判断するピーク面積比較工程と、
を含むことを特徴とする請求項４に記載のクロマトグラムデータ処理方法。
測定対象試料に含まれる複数の成分が分離カラムで分離され、分離された前記成分を検出して得られるクロマトグラムにより前記成分を測定する際にコンピュータに実行させるクロマトグラムデータ処理プログラムであって、
時間と信号強度とをそれぞれ軸として、複数の前記成分に由来するピークを有する前記クロマトグラムを得るクロマトグラム取得部と、
前記クロマトグラム内の複数の前記ピークに対して、ガウシアン関数又は指数関数で修飾されたガウシアン関数と、ピアソン関数とを用いて、フィッティング処理を行い、複数の前記ピークを分解するフィッティング部と、
前記フィッティング処理を行った後の前記クロマトグラムから前記成分を判定する判定部と、
をコンピュータに実行させることを特徴とするクロマトグラムデータ処理プログラム。
請求項６に記載のクロマトグラムデータ処理プログラムを格納したコンピュータで読み取り可能な記憶媒体。