特開 2022-62493 マススペクトル処理装置及び方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】公開特許公報(A)

(11)【公開番号】P2022062493

(43)【公開日】2022-04-20

(54)【発明の名称】マススペクトル処理装置及び方法

(51)【国際特許分類】

   G01N 30/86 20060101AFI20220413BHJP

   G01N 27/62 20210101ALI20220413BHJP

   G01N 30/72 20060101ALI20220413BHJP

【ＦＩ】

G01N30/86 L

G01N27/62 C

G01N30/86 D

G01N30/86 E

G01N30/86 R

G01N30/72 A

【審査請求】有

【請求項の数】10

【出願形態】ＯＬ

(21)【出願番号】P 2020170540

(22)【出願日】2020-10-08

(71)【出願人】

【識別番号】000004271

【氏名又は名称】日本電子株式会社

(74)【代理人】

【識別番号】110001210

【氏名又は名称】特許業務法人ＹＫＩ国際特許事務所

(72)【発明者】

【氏名】田中 博一

【テーマコード（参考）】

2G041

【Ｆターム（参考）】

2G041CA01

2G041DA13

2G041EA06

2G041GA03

2G041GA06

2G041HA01

2G041MA05

(57)【要約】

【課題】クロマトグラムに含まれる微小なピークをノイズから区別するための情報を生成する。
【解決手段】クロマトグラム１６８と共に、第１トレンドチャート１７０及び第２トレンドチャート１７２が表示される。第１トレンドチャート１７０は、複数のマススペクトルから求められる複数の第１代表値列に基づいて生成される。第２トレンドチャート１７２は、複数のマススペクトルから求められる複数の第２代表値列に基づいて生成される。第１トレンドチャート１７０及び第２トレンドチャート１７２に基づいて、マススペクトル安定期間が判定される。
【選択図】図１２

【特許請求の範囲】

【請求項1】

時間軸上に並ぶ複数のマススペクトルからなるマススペクトル列に基づいて、クロマトグラムを生成するクロマトグラム生成部と、
前記マススペクトル列を構成するマススペクトルごとに当該マススペクトルに含まれる複数の代表ピークを特定し、これにより複数の代表ピークの時間的変化が反映された複数のグラフからなるトレンドチャートを生成するチャート生成部と、
前記クロマトグラムと共に、前記トレンドチャートを表示し又は前記トレンドチャートから生成されるマススペクトル安定情報を表示する表示処理部と、
を含むことを特徴とするマススペクトル処理装置。

【請求項2】

請求項１記載のマススペクトル処理装置において、
前記チャート生成部は、
前記マススペクトルに含まれる複数の代表ピークを抽出する抽出部と、
前記複数の代表ピークに基づいて、前記マススペクトルの特徴量として、複数の代表値からなる代表値列を特定する特定部と、
前記複数のマススペクトルから求められた複数の代表値列に基づいて、前記複数のグラフとして、代表値列の時間的変化を表す複数のグラフを生成するグラフ生成部と、
を含む、ことを特徴とするマススペクトル処理装置。

【請求項3】

請求項２記載のマススペクトル処理装置において、
前記抽出部は、前記複数の代表ピークとして、前記マススペクトルにおける上位ｋ個（但しｋは２以上の整数）のピークを抽出する、
ことを特徴とするマススペクトル処理装置。

【請求項4】

請求項３記載のマススペクトル処理装置において、
前記特定部は、前記代表値列として、前記上位ｋ個のピークに対応するｋ個のｍ／ｚからなるｍ／ｚ列を求め、
前記グラフ生成部は、保持時間軸とｍ／ｚ軸とを有する座標系に対して、前記複数のマススペクトルから求められた複数のｍ／ｚ列をプロットすることにより、ｋ個のグラフを生成する、
ことを特徴とするマススペクトル処理装置。

【請求項5】

請求項３記載のマススペクトル処理装置において、
前記ｋは３以上であり、
前記上位ｋ個のピークは、最上位ピーク及びｋ－１個の上位ピークにより構成され、
前記特定部は、前記代表値列として、前記最上位ピークの強度に対する前記ｋ－１個の上位ピークの強度の比率からなる比率列を求め、
前記グラフ生成部は、保持時間軸と比率軸とを有する座標系に対して、前記複数のマススペクトルから求められた複数の比率列をプロットすることにより、ｋ－１個のグラフを生成する、
ことを特徴とするマススペクトル処理装置。

【請求項6】

請求項１記載のマススペクトル処理装置において、
前記チャート生成部は、
前記マススペクトルに含まれる複数の第１代表ピークを抽出する第１抽出部と、
前記マススペクトルに含まれる複数の第２代表ピークを抽出する第２抽出部と、
前記複数の第１代表ピークに基づいて、前記マススペクトルの第１特徴量として、複数の第１代表値からなる第１代表値列を求める第１特定部と、
前記複数の第２代表ピークに基づいて、前記第１特徴量とは異なる前記マススペクトルの第２特徴量として、複数の第２代表値からなる第２代表値列を求める第２特定部と、
前記複数のマススペクトルから求められた複数の第１代表値列に基づいて、前記複数のグラフとして、第１代表値列の時間的変化を表す複数の第１グラフを生成する第１グラフ生成部と、
前記複数のマススペクトルから求められた複数の第２代表値列に基づいて、前記複数のグラフとして、第２代表値列の時間的変化を表す複数の第２グラフを生成する第２グラフ生成部と、
を含み、
前記複数の第１グラフにより第１トレンドチャートが生成され、前記複数の第２グラフにより第２トレンドチャートが生成される、
ことを特徴とするマススペクトル処理装置。

【請求項7】

請求項１記載のマススペクトル処理装置において、
前記トレンドチャートの解析によりマススペクトル安定期間を判定するチャート解析部を含み、
前記クロマトグラムと共に、前記マススペクトル安定期間を示す前記マススペクトル安定情報が表示される、
ことを特徴とするマススペクトル処理装置。

【請求項8】

請求項１記載のマススペクトル処理装置において、
前記トレンドチャートに基づいて判定されたマススペクトル安定期間に応じて、前記クロマトグラムに含まれるピークに対する処理を実行する処理部を含む、
ことを特徴とするマススペクトル処理装置。

【請求項9】

時間軸上に並ぶ複数のマススペクトルからなるマススペクトル列に基づいて、マススペクトルごとに当該マススペクトルに含まれる複数の代表ピークからなる代表ピーク列を特定する工程と、
前記複数のマススペクトルから特定された複数の代表ピーク列に基づいて、複数の代表ピークの時間的変化が反映された複数のグラフからなるトレンドチャートを生成する工程と、
を含むことを特徴とするマススペクトル処理方法。

【請求項10】

情報処理装置において実行されるプログラムであって、
時間軸上に並ぶ複数のマススペクトルからなるマススペクトル列に基づいて、マススペクトルごとに当該マススペクトルに含まれる複数の代表ピークからなる代表ピーク列を特定する機能と、
前記複数のマススペクトルから特定された複数の代表ピーク列に基づいて、複数の代表ピークの時間的変化が反映された複数のグラフからなるトレンドチャートを生成する機能と、
を含むことを特徴とするプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、マススペクトル処理装置及び方法に関し、特に、マススペクトル列に基づいてクロマトグラムを生成するマススペクトル処理装置及び方法に関する。

【背景技術】

【0002】

質量分析システムは、例えば、クロマトグラフ、質量分析装置及びマススペクトル処理装置により構成される。クロマトグラムにおいて、試料から時間的に分離された複数の成分が取り出され、それらが質量分析装置へ順次送られる。質量分析装置においては、導入される複数の成分に対して質量分析が繰り返し実行される。マススペクトル処理装置は、質量分析の結果に基づいて時間軸上に並ぶ複数のマススペクトルからなるマススペクトル列を生成し、マススペクトル列に基づいてクロマトグラムを生成する。クロマトグラムに含まれる複数のピークがクロマトグラフで分離、抽出された複数の成分に対応する。個々のピークの解析により、具体的には、個々のピークに対応するマススペクトルの解析により、各成分が定性分析あるいは定量分析される。

【0003】

なお、特許文献１には、クロマトグラム上にピークトップ位置を示すマークを表示する技術が開示されている。特許文献２には、クロマトグラムに含まれるピークがノイズに相当するか否かを判断する技術が開示されている。特許文献１及び特許文献２には、クロマトグラム評価用のグラフ又はチャートを生成することについては記載されていない。

【先行技術文献】

【特許文献】

【0004】

【特許文献1】特許第６２７３７８５号公報

【特許文献2】特開２００５－２２１２７６号公報

【発明の概要】

【発明が解決しようとする課題】

【0005】

試料から抽出された各成分のマススペクトルを特定するために、クロマトグラムに含まれる個々のピークが自動的に検出されあるいは目視により特定される。ピークの自動検出に際して、微小なピークまで検出しようとすると、ノイズまでが検出されてしまう。ノイズを避けてピーク検出を行おうとすると、微小なピークを見逃してしまう。一方、クロマトグラムにおいて、特にそれ全体が一画面内に表示された場合において、目視によりピーク、特に微小ピークを特定することは困難である。自動的なピーク検出や目視によるピーク特定を支援する技術の実現が求められている。

【0006】

本発明の目的は、クロマトグラムに含まれるピークを検出又は特定し易くすることにある。あるいは、本発明の目的は、クロマトグラムに含まれる微小なピークをノイズから区別するための情報を簡易な構成で生成することにある。

【課題を解決するための手段】

【0007】

本開示に係るマススペクトル処理装置は、時間軸上に並ぶ複数のマススペクトルからなるマススペクトル列に基づいて、クロマトグラムを生成するクロマトグラム生成部と、前記マススペクトル列を構成するマススペクトルごとに当該マススペクトルに含まれる複数の代表ピークを特定し、これにより複数の代表ピークの時間的変化が反映された複数のグラフからなるトレンドチャートを生成するチャート生成部と、前記クロマトグラムと共に、前記トレンドチャートを表示し又は前記トレンドチャートから生成されるマススペクトル安定情報を表示する表示処理部と、を含むことを特徴とする。

【0008】

本発明に係るマススペクトル処理方法は、時間軸上に並ぶ複数のマススペクトルからなるマススペクトル列に基づいて、マススペクトルごとに当該マススペクトルに含まれる複数の代表ピークからなる代表ピーク列を特定する工程と、前記複数のマススペクトルから特定された複数の代表ピーク列に基づいて、複数の代表ピークの時間的変化が反映された複数のグラフからなるトレンドチャートを生成する工程と、を含むことを特徴とする。

【発明の効果】

【0009】

本発明によれば、クロマトグラムに含まれるピークを検出又は特定し易くなる。あるいは、本発明によれば、クロマトグラムに含まれる微小なピークをノイズから区別するための情報を簡易な構成で生成できる。

【図面の簡単な説明】

【0010】

【図1】実施形態に係る質量分析システムを示すブロック図である。

【図2】第１トレンドチャートの生成を説明するための図である。

【図3】マススペクトルにおける複数の代表ピークを示す図である。

【図4】第１トレンドチャートの第１例を示す図である。

【図5】第１トレンドチャートの第２例を示す図である。

【図6】第２トレンドチャートの生成を説明するための図である。

【図7】第２トレンドチャートの第１例を示す図である。

【図8】第２トレンドチャートの第２例を示す図である。

【図9】第２トレンドチャートの第３例を示す図である。

【図10】スペクトル安定期間の判定方法を示す図である。

【図11】第１トレンドチャートの第３例を示す図である。

【図12】第１表示例を示す図である。

【図13】第２表示例を示す図である。

【図14】拡大表示を説明するための図である。

【図15】類似関係が判定された場合における表示例を示す図である。

【図16】類似関係の判定方法を示す図である。

【図17】主抽出位置の補正を説明するための図である。

【図18】２つの副抽出位置の補正を説明するための図である。

【発明を実施するための形態】

【0011】

以下、実施形態を図面に基づいて説明する。

【0012】

（１）実施形態の概要
実施形態に係るマススペクトル処理装置は、クロマトグラム生成部、チャート生成部、及び、表示処理部を有する。クロマトグラム生成部は、時間軸上に並ぶ複数のマススペクトルからなるマススペクトル列に基づいて、クロマトグラムを生成する。チャート生成部は、マススペクトル列を構成するマススペクトルごとに当該マススペクトルに含まれる複数の代表ピークを特定し、これにより複数の代表ピークの時間的変化が反映された複数のグラフからなるトレンドチャートを生成する。表示処理部は、クロマトグラムと共に、トレンドチャートを表示し又はトレンドチャートから生成されるマススペクトル安定情報を表示する。

【0013】

上記構成によれば、マススペクトルの時間的変化が代表ピーク列の時間的変化として捉えられる。代表ピーク列の時間的変化がトレンドチャートとして表示され、あるいは、トレンドチャートから生成されるマススペクトル安定情報が表示される。トレンドチャートの生成に当たって、各マススペクトルそれ全体が参照されるのではなく、各マススペクトルを代表する複数の代表ピークが参照されるので、ノイズの影響を受け難くなり、また計算量を削減できる。

【0014】

実施形態においては、複数の代表ピークとして、高い強度をもった複数のピークが選ばれる。トレンドチャートは、複数のグラフにより構成され、それはマススペクトルの時間的変化の傾向又はパターンを示すものである。複数のグラフの状態（例えば、密集、偏移、安定、変化等）から、個々の時刻におけるマススペクトルの状態を推認することが可能となる。

【0015】

実施形態において、チャート生成部は、抽出部、特定部、及び、グラフ生成部を有する。抽出部は、マススペクトルに含まれる複数の代表ピークを抽出する。特定部は、複数の代表ピークに基づいて、マススペクトルの特徴量として、複数の代表値からなる代表値列を求める。グラフ生成部は、複数のマススペクトルから求められた複数の代表値列に基づいて、代表値列の時間的変化を表す複数のグラフを生成する。

【0016】

実施形態において、抽出部は、複数の代表ピークとして、マススペクトルにおける上位ｋ個（但しｋは２以上の整数）のピークを抽出する。ｋが小さ過ぎると、マススペクトルの特徴量抽出が不十分となる。ｋが大き過ぎると、演算量が増大してしまう。状況に応じてｋの値を定めるのが望ましい。ここで、ｋは後述するＭ及びＮに相当するパラメータである。なお、ｋ＝１とする変形例も考えられる。

【0017】

実施形態において、特定部は、代表値列として、上位ｋ個のピークに対応するｋ個のｍ／ｚからなるｍ／ｚ列を求める。グラフ生成部は、保持時間軸とｍ／ｚ軸とを有する座標系に対して、複数のマススペクトルから求められた複数のｍ／ｚ列をプロットすることにより、ｋ個のトレンドグラフを生成する。この構成によれば、ｋ個のｍ／ｚの時間変化がｋ個のグラフとして表現される。

【0018】

実施形態において、ｋは３以上である。上位ｋ個のピークは、最上位ピーク及びｋ－１個の上位ピークにより構成される。特定部は、代表値列として、最上位ピークの強度に対するｋ－１個の上位ピークの強度の比率からなる比率列を求める。グラフ生成部は、保持時間軸と比率軸とを有する座標系に対して、複数のマススペクトルから求められた複数の比率列をプロットすることにより、ｋ－１個のグラフを生成する。この構成によれば、最上位ピークと上位ｋ個のピークの強度関係の時間的変化を表示し得る。なお、ｋ＝２とする変形例も考えられる。

【0019】

実施形態において、チャート生成部は、第１抽出部、第２抽出部、第１特定部、第２特定部、第１グラフ生成部、及び、第２グラフ生成部を有する。第１抽出部は、マススペクトルに含まれる複数の第１代表ピークを抽出する。第２抽出部は、マススペクトルに含まれる複数の第２代表ピークを抽出する。第１特定部は、複数の第１代表ピークに基づいて、マススペクトルの第１特徴量として、複数の第１代表値からなる第１代表値列を求める。第２特定部は、複数の第２代表ピークに基づいて、第１特徴量とは異なるマススペクトルの第２特徴量として、複数の第２代表値からなる第２代表値列を求める。第１グラフ生成部は、複数のマススペクトルから求められた複数の第１代表値列に基づいて、第１代表値列の時間的変化を表す複数の第１グラフを生成する。第２グラフ生成部は、複数のマススペクトルから求められた複数の第２代表値列に基づいて、第２代表値列の時間的変化を表す複数の第２グラフを生成する。複数の第１グラフにより第１トレンドチャートが生成される。複数の第２グラフにより第２トレンドチャートが生成される。

【0020】

上記構成によれば、第１トレンドチャート及び第２トレンドチャートの両者を考慮して各時刻におけるマススペクトルの状態を評価し得る。複数の第１代表ピークと複数の第２代表ピークとが同一であってもよい。

【0021】

実施形態に係るマススペクトル処理装置は、トレンドチャートの解析によりマススペクトル安定期間を判定するチャート解析部を含む。クロマトグラムと共に、マススペクトル安定期間を示すマススペクトル安定情報が表示される。マススペクトル安定情報の参照により、解析対象となるマススペクトルを選定し易くなる。

【0022】

実施形態に係るマススペクトル処理装置は、トレンドチャートに基づいて特定されたマススペクトル安定期間に応じて、クロマトグラムに含まれるピークに対する処理を実行する処理部を含む。その処理には、ピーク頂点の特定、ピーク前後のピーク外側位置の特定、ピークに含まれる複数のサブピークの分離、複数のサンプルから得られた複数のクロマトグラムの比較、等が含まれる。

【0023】

実施形態に係るマススペクトル処理方法は、特定工程、及び、生成工程を有する。特定工程では、時間軸上に並ぶ複数のマススペクトルからなるマススペクトル列に基づいて、マススペクトルごとに当該マススペクトルに含まれる複数の代表ピークからなる代表ピーク列が特定される。生成工程では、複数のマススペクトルから特定された複数の代表ピーク列に基づいて、複数の代表ピークの時間的変化が反映された複数のグラフからなるトレンドチャートが生成される。

【0024】

上記マススペクトル処理方法は、ハードウエアの機能として又はソフトウエアの機能として実現され得る。後者の場合、マススペクトル処理方法を実行するプログラムが、ネットワークを介して又は可搬型記憶媒体を介して、情報処理装置にインストールされる。情報処理装置の概念には、コンピュータ、質量分析装置、及び、質量分析システムが含まれる。

【0025】

（２）実施形態の詳細
図１には、実施形態に係る質量分析システムの全体構成がブロック図として示されている。図示された質量分析システム１０は、測定部１２、及び、情報処理装置１４を有する。質量分析システム１０は、いわゆる、ＧＣ－ＭＳ（ガスクロマトグラフ質量分析計）である。具体的には、情報処理装置１４は、コンピュータ等により構成され、それはマススペクトル処理装置として機能する。測定部１２は、ガスクロマトグラフ（ＧＣ）１６、及び、質量分析装置１８により構成される。質量分析装置１８の前段に、他の前処理部が設けられてもよい。

【0026】

ＧＣ１６は、試料に含まれる複数の成分を時間的に分離するカラムを有している。すなわち、カラムにより、試料から複数の成分が取り出される。分離された複数の成分が質量分析装置１８中のイオン源２０に順次送られる。

【0027】

質量分析装置１８は、イオン源２０、質量分析器２２、及び、検出器２４により構成される。イオン源２０は、例えば、ＥＩ法（電子イオン化法）に従うイオン源である。他のイオン源が用いられてもよい。

【0028】

質量分析器２２は、イオン源２０から出た各イオンの質量電荷比（ｍ／ｚ）を分析する部分である。質量分析器２２として、飛行時間型質量分析器や四重極型質量分析器等を利用し得る。検出器２４において個々のイオンが検出される。検出器２４から出力された検出信号２６は、図示されていない信号処理回路に送られる。検出信号２６は、ｍ／ｚごとのイオン強度を示す信号である。

【0029】

次に、情報処理装置１４について説明する。情報処理装置１４は、プロセッサ（例えばＣＰＵ）、表示器４２、及び、入力器４４を有する。図１においては、プロセッサが実現する複数の機能が複数のブロックとして表現されている。

【0030】

スペクトル生成部２８は、質量分析装置１８から順次出力される検出信号に基づいてマススペクトルを生成する。スペクトル生成部２８により、保持時間軸上において並ぶ複数のマススペクトルからなるマススペクトル列が生成される。図示の構成例では、マススペクトル列が記憶部３０に記憶されている。記憶部３０は、メモリにより構成される。

【0031】

ＴＩＣＣ（トータルイオンカレントクロマトグラム）生成部３４は、マススペクトルごとに、マススペクトルを積算してＴＩＣ（トータルイオンカレント）を演算し、ＴＩＣの時間的変化を表したグラフとして、ＴＩＣＣを生成する。ＴＩＣＣを表すデータが表示処理部３２に送られている。

【0032】

第１トレンドチャート生成部３６は、マススペクトル列に基づいて第１トレンドチャートを生成するモジュールである。第１トレンドチャートは、横軸としての保持時間軸、及び、縦軸としてのｍ／ｚ軸を有し、マススペクトルの時間的変化を複数の第１代表値（具体的にはマススペクトルを代表する複数のｍ／ｚ）の時間的変化として表現した複数の第１グラフの集合体である。第１トレンドチャートからマススペクトル安定期間を自動的に検出し又は目視で特定することが可能となる。各マススペクトルを代表する複数のｍ／ｚは第１特徴量と言い得る。第１トレンドチャートについては後に詳述する。

【0033】

第２トレンドチャート生成部３８は、マススペクトル列に基づいて第２トレンドチャートを生成するモジュールである。第２トレンドチャートは、横軸としての保持時間軸、及び、縦軸としての比率軸を有し、マススペクトルの時間的変化を第２代表値列（具体的にはマススペクトルの状態を示す複数の比率であり、その詳細については後述する。）の時間的変化として表現した複数の第２グラフの集合体である。第２トレンドチャートからマススペクトル安定期間を自動的に検出し又は目視で特定することが可能となる。個々のマススペクトルごとに特定される複数の比率は、第１特徴量とは異なる第２特徴量と言い得る。第２トレンドチャートについては後に詳述する。

【0034】

表示処理部３２は、表示器４２に表示される画像を生成するものである。その画像には、クロマトグラム、第１トレンドチャート、第２トレンドチャート、スペクトル安定情報、等が含まれる。表示器４２には、必要に応じて、マススペクトルや条件設定欄が表示される。

【0035】

表示処理部３２は、チャート解析部４０を有する。チャート解析部４０は、第１トレンドチャート及び第２トレンドチャートに基づいて、マススペクトル安定期間を判定する。また、チャート解析部４０は、類似ピーク検索機能、等を有する。表示処理部３２には、クロマトグラムを解析するクロマトグラム解析部も含まれる。クロマトグラム解析部は、クロマトグラムにおけるピークをノイズと区別しながら検出する機能を備える。その際には、マススペクトル安定期間が参照される。入力器４４を利用して、マススペクトル安定期間の判定条件を設定し得る。

【0036】

図２には、第１トレンドチャート生成アルゴリズムが示されている。マススペクトル列を構成する個々のマススペクトル５０から、強度順で、上位Ｍ個のピークが特定される（符号５２－１から符号５２－ｍを参照）。上位Ｍ個のピークは、Ｍ個の代表ピークと言い得る。ここで、Ｍは２以上の整数であり、例えば、Ｍは５である。その場合、１番目に大きな強度を有する第１ピークから５番目に大きな強度を有する第５ピークまでの５つのピークが特定される。それらはピーク列を構成する。そのピーク列は、代表ピーク列と言い得る。強度順で上位からＭ個のピークを参照することにより、マススペクトル全体を参照する場合に比べて、ノイズの影響を受け難く、また計算量を削減できる。

【0037】

続いて、第１ピークから第Ｍピークに対応するＭ個のｍ／ｚが特定され、それが第１座標系５４にプロットされる（符号５６－１から符号５６－ｍを参照）。Ｍ個のｍ／ｚは、第１ピークに対応するｍ／ｚ、第２ピークに対応するｍ／ｚ、・・・、第Ｍピークに対応するｍ／ｚである。それらはそれぞれ代表値であり、それら全体はマススペクトルの特徴量又は特徴ベクトルと言い得る。第１座標系５４は、横軸としての保持時間軸、及び、縦軸としてのｍ／ｚ軸を有する。マススペクトル生成タイミングごとに、Ｍ個のｍ／ｚを表すＭ個の点がプロットされる。各点の属性情報として、順位（ピーク強度順位）が管理される。

【0038】

今回のマススペクトル及び前回のマススペクトルに対応する２つの点列の間が、Ｍ個の直線で連結される（符号５８－１から符号５８－ｍを参照）。その場合、同じ順位に対応付けられた２つの点ごとに、２つの点が連結される。

【0039】

以上の処理の繰り返しにより、第１グラフから第ＭグラフまでのＭ個のグラフが生成される（符号６０－１から符号６０－ｍを参照）。個々のグラフは折れ線グラフである。折れ線グラフに代えて曲線グラフを採用してもよい。第１グラフ６０－１は、第１順位に対応する複数の点を通るグラフであり、第２グラフ６０－２は第２順位に対応する複数の点を通るグラフであり、・・・、第Ｍグラフ６０－ｍは第Ｍ順位に対応するグラフである。それらのＭ個のグラフそれら全体として、第１トレンドチャート６１が構成される。第１トレンドチャート６１は、代表ｍ／ｚトレンドチャートあるいは高強度ピーク位置トレンドチャートとも言い得る。Ｍの値は、自動的に設定され又はユーザーにより指定される。

【0040】

図３には、マススペクトルの一例が示されている。横軸はｍ／ｚ軸であり、縦軸は強度軸である。マススペクトル６２には、多数のピークが含まれる。図示の例では、多数のピークの内で、上位５個のピークが代表ピーク列として特定される（符号６４－１から符号６４－５を参照）。続いて、上位５個のピークに対応する５個のｍ／ｚが特定される（ｍ／ｚ－１からｍ／ｚ－５を参照）。それらは第１代表値列であり、それら全体がマススペクトルの第１特徴量に相当する。

【0041】

図４の上段には、第１例に係るクロマトグラム６６が示され、図４の下段には、第１例に係る第１トレンドチャート６８が示されている。クロマトグラム６６が有する保持時間軸のスケールと第１トレンドチャート６８が有する保持時間軸のスケールは同一である。このことは、以下に説明する、複数の保持軸を示した図５、図７－１５、図１７－１８においても同様である。

【0042】

図４に示す例において、第１トレンドチャート６８は、５つのグラフ６８ａ,６８ｂ，６８ｃ，６８ｄ，６８ｅにより構成される。保持時間軸に沿って、５つのグラフ６８ａ,６８ｂ，６８ｃ，６８ｄ，６８ｅは、交差したり平行に走行したりしている。符号７２Ａ，７２Ｂ，７２Ｃは、５つのグラフ６８ａ,６８ｂ，６８ｃ，６８ｄ，６８ｅにおいて、概ね平行に走行している３つの部分を示しており、それらはクロマトグラム６６に含まれる３つのピーク７０Ａ，７０Ｂ，７０Ｃに対応している。

【0043】

すなわち、第１トレンドチャート６８において、交差がない又は交差があまりない単調な部分は、マススペクトルが安定している期間であり、つまり成分が抽出されている期間であると考えられる。例えば、第１トレンドチャートにおいて、符号８６で示す部分においては、５つのグラフ６８ａ,６８ｂ，６８ｃ，６８ｄ，６８ｅが不安定となっており、少なからずの交差箇所つまり頻繁な交差が生じている。そのような期間においては、マススペクトルの形態又は内容が不安定であることが推認され、そのような期間において、クロマトグラムにピークがあったとしても、それはノイズである可能性が高いと評価し得る。符号７８，８０，８４で示す部分についても同様のことが言える。

【0044】

ちなみに、第１トレンドチャート６８において、交差がまったく生じていない安定した期間７３には、クロマトグラム６６内のピーク７０Ａが含まれ、更にピーク７０Ａの後側（図４において右側）の裾部分も含まれる。その裾部分もピーク７０Ａの一部であると評価し得る。

【0045】

以上のように、第１トレンドチャートによれば、マススペクトルの時間的変化又は各時点での様子を分かり易く表現することが可能であり、その第１トレンドチャートを通じてマススペクトルが安定している期間を比較的に容易に特定することが可能である。これにより、化合物由来のピークであるか否かを判定することが可能となる。

【0046】

なお、第１トレンドチャート６８において、部分７２Ｂの両端部には、孤立した交差７４，７６が認められる。ユーザーにより又は自動的に設定された判定条件に従って、交差７４，７６を避けて安定期間が判定され、あるいは、それも含めて安定期間が判定される。

【0047】

図５には、第２例に係るクロマトグラム８８及び第２例に係る第１トレンドチャート９０が示されている。第１トレンドチャート９０において、符号９２，９６で示す部分は安定しており、一方、符号１００で示す部分は不安定である。クロマトグラム８８においては、緩やかなピーク９４，９８が認められ、それらは符号９２，９６で示す部分に対応している。すなわち、第１トレンドチャート９０に基づいて、クロマトグラム８８に含まれる個々の山状部分について真のピークであるかノイズであるかを判別し得る。例えば、符号１０２で示す部分には複数の微小な山が認められるが、それに対応する部分１００を参照すると、その部分１００は不安定である。このことから、符号１０２で示す部分に含まれる複数の微小な山はノイズであるとみなせる。

【0048】

安定している部分９２の長さとして、安定期間９３を画定でき、安定期間９３からピーク９４の両端を特定し得る。例えば、ピークトップに相当する主抽出位置、及び、ピーク外側の２つのノイズ参照位置に相当する２つの副抽出位置を安定期間９３から容易に特定し得る。

【0049】

図６には、第２トレンドチャート生成アルゴリズムが示されている。マススペクトル列を構成する個々のマススペクトル５０から、強度順で、上位Ｎ個のピークが特定される（符号１０４－１から符号１０４－ｎを参照）。上位Ｎ個のピークは、Ｎ個の代表ピークと言い得る。ここで、Ｎは３以上の整数であり、例えば、Ｎは６である。その場合、１番目に大きな強度を有する第１ピークから６番目に大きな強度を有する第６ピークまでが特定される。そのピーク列は、代表ピーク列と言い得る。強度順で上位Ｎ個のピークを参照することにより、マススペクトル全体を参照する場合に比べて、ノイズの影響を受け難く、また計算量を削減できる。

【0050】

具体的には、Ｎ個のピークは、１番目に大きな強度を有する第１ピーク（最上位ピーク）、及び、２番目からＮ番目に大きなピークを有するＮ－１個の上位ピークからなる。最上位ピークの強度が基準強度とされる。Ｎ－１個の上位ピークの強度を基準強度で除することにより、Ｎ－１個の比率が演算され、Ｎ－１個の比率が第２座標系１０６にプロットされる（符号１０８－１から符号１０８－ｎ－１を参照）。プロットされる各点に対しては、その属性情報として、比率の大きさの順位（比率順位）が対応付けられる。Ｎ－１個の比率はそれぞれ代表値であり、それら全部はマススペクトルの特徴量又は特徴ベクトルと言い得る。

【0051】

第２座標系１０６は、横軸としての保持時間軸、及び、縦軸としての比率軸を有する。マススペクトルごとにＮ－１個の点（点列）をプロットしていく過程の中で、隣り合う２つの点列がＮ－１個の線で連結される。その場合、比率の大きさ順で、同じ順位を有する２つの点ごとに、２つの点が相互に連結される（符号１１０－１から符号１１０－ｎ－１を参照）。

【0052】

以上の処理の繰り返しにより、第１グラフから第Ｎ－１グラフまでのＮ－１個のグラフが生成される（符号１１２－１から符号１１２－ｎ－１を参照）。個々のグラフは折れ線グラフである。折れ線グラフに代えて曲線グラフを採用してもよい。

【0053】

具体的には、第１グラフ１１２－１は、第１順位に対応する複数の点を通るグラフであり、第２グラフ１１２－２は第２順位に対応する複数の点を通るグラフであり、・・・、第Ｎ－１グラフ１１２－ｎ－１は第Ｎ－１順位に対応するグラフである。それらのＮ－１個のグラフにより第２トレンドチャート１１４が構成される。第２トレンドチャート１１４は、比率トレンドチャート、強度関係トレンドチャート、又は、Ｉ（Identification ion）／Ｑ（Quantification ion）チャート、とも言い得る。Ｎの値は、自動的に設定され又はユーザーにより指定される。

【0054】

図７の上段には、図４に示したクロマトグラム６６が示され、図４の下段には、第１例に係る第２トレンドチャート１１６が示されている。

【0055】

第２トレンドチャート１１６は、図示の例において、５つのグラフ１１６ａ,１１６ｂ，１１６ｃ，１１６ｄ，１１６ｅにより構成される。保持時間軸に沿って、５つのグラフ１１６ａ,１１６ｂ，１１６ｃ，１１６ｄ，１１６ｅが上下に並んで走行している。符号１１８Ａ，１１８Ｂ，１１８Ｃは、５つのグラフ１１６ａ,１１６ｂ，１１６ｃ，１１６ｄ，１１６ｅにおいて概ね平行に走行している安定な部分を示している。特に、それらの部分１１８Ａ，１１８Ｂ，１１８Ｃにおいては、５つのグラフ１１６ａ,１１６ｂ，１１６ｃ，１１６ｄ，１１６ｅのそれぞれが単調であり、それらの中心又は重心が下方にシフトしている。それらの部分１１８Ａ，１１８Ｂ，１１８Ｃは、クロマトグラム６６に含まれる３つのピーク７０Ａ，７０Ｂ，７０Ｃに対応している。

【0056】

すなわち、第２トレンドチャート１１６において、単調な部分つまり安定している部分は、マススペクトルが安定していることを示しており、つまり成分が抽出されている期間に相当すると考えられる。例えば、第２トレンドチャート１１６において、符号１２２、１２４で示す部分は不安定であり、それは、マススペクトルの形態が安定していないことを示している。そのような期間において、クロマトグラムにピークが含まれていたとしても、それはノイズである可能性が高い。

【0057】

第２トレンドチャート１１６においては、複数のグラフ１１６ａ,１１６ｂ，１１６ｃ，１１６ｄ，１１６ｅの相互間の接触の可能性はあるものの、それらの交差は生じないので、その内容を観察し易いという利点を得られる。

【0058】

図８の上段には、図５に示したクロマトグラム１０２が示されており、図８の下段には、第２例に係る第２トレンドチャート１２４が示されている。第２トレンドチャート１２４には、下側に偏移した単調な部分１２６，１２８が含まれる。それらの部分１２６，１２８は、クロマトグラム１０２に含まれるピーク９４，９８に対応している。例えば、５つのグラフ１１６ａ,１１６ｂ，１１６ｃ，１１６ｄ，１１６ｅのばらつきや平均レベル等から、単調な部分１２６，１２８が決定されてもよい。

【0059】

図９の上段には、保持時間（測定時間）の全体にわたるクロマトグラム１３２が示されている。図９の下段には、第３例に係る第２トレンドチャート１３４が示されている。それも保持時間の全体に対応している。

【0060】

第２トレンドチャート１３４は、１０本のグラフにより構成されており、つまり上記Ｎは１１である。第２トレンドチャート１３４において、Ｎの値を大きくしていくと、グラフ集合体の下縁が下方に下がっていく。具体的には、Ｎの値を大きくしていくと、Ｎ－１番目のピークの強度が基準強度に対してより小さくなり、Ｎ－１番目のピークから演算された比率が０により近付いていく。

【0061】

第２トレンドチャート１３４において、複数の谷１３８が存在しており、それらはクロマトグラム１３２に含まれる複数のピーク１３６に対応している。例えば、クロマトグラム１３２において、第２トレンドチャート１３４における谷１３９に対応する位置には明確なピークが認められないが、クロマトグラム１３２に対する横軸方向及び縦軸方向の拡大により、谷１３９に対応するピークを発見できる可能性が高い。

【0062】

以上のように、第２トレンドチャートによれば、マススペクトルの時間的変化又は各時点での様子を分かり易く表現することが可能であり、第２トレンドチャートを通じてマススペクトルが安定している期間を比較的に容易に特定することが可能である。これにより、クロマトグラムに含まれる各ピークが化合物に由来するものであるか否かを容易に判定することが可能である。

【0063】

図１０には、マススペクトルの安定を判定するためのアルゴリズムが示されている。スペクトル列１４０に基づいて第１トレンドチャート１４２が生成され、また、スペクトル列１４０に基づいて第２トレンドチャート１４４が生成される。第１トレンドチャート１４２の内で、第１安定判定条件１４６を満たす部分が第１安定期間として判定される。同様に、第２トレンドチャート１４４の内で、第２安定判定条件１４８を満たす部分が第２安定期間として判定される。総合判定１５０では、第１安定期間及び第２安定期間に基づいて、最終的な安定期間が判断される。その場合には、ＡＮＤ条件が適用されてもよいし、ＯＲ条件が適用されてもよい。なお、一方の安定期間のみを採用してもよい。

【0064】

第１安定判定条件に関し、例えば、保持時間軸上で並ぶｋ個のマススペクトルにわたって上位Ｍ個のｍ／ｚの並びが同一である場合に、安定を判定してもよい。マススペクトルが安定していても、強度がほぼ等しい２つの代表ピーク間で、それらの順位が入れ替わる可能性もある。それを考慮し、上位Ｍ個のｍ／ｚ中の所定割合数で一致が認められた場合に、安定を判定してもよい。あるいは、並び（順列）の一致ではなく、組み合わせの一致に基づいて、安定を判定してもよい。上記ｋについて上限を設定してもよい。一致回数が上限に達した場合には安定判定を終了させることにより、同じ成分が出続けているような場合にそれを判定対象から除外し得る。

【0065】

図１１の上段には第３例に係るクロマトグラム１５２が示され、図１１の下段には第３例に係る第１トレンドチャート１５４が示されている。クロマトグラム１５２にはピーク１５６が含まれる。第１トレンドチャート１５４において、ピーク１５６に対応する部分は区間１６４である。区間１６４内の中央において、孤立した交差１６２が生じている。第１判定条件を厳しく設定した場合、交差１６２が生じた時点において安定が判定されず、分離した２つの区間１５８，１６０がそれぞれ安定期間であると判定される。その場合、ピーク１５６の全体を１つのピーク範囲内であると判断できなくなる。そこで、例えば、第１判定条件として、組み合わせの一致を求めることが考えられる。その場合、交差１６２があっても、安定を判定し続けることが可能となる。結果として、区間１６４それ全体を安定期間として定めることが可能となる。

【0066】

図１０に示した第２安定判定条件に関し、比率列についての評価値、例えば、相対標準偏差に基づいて、判定を判定してもよい。その場合、注目しているマススペクトルから得られた比率列に加えて、その前後に存在する所定個のマススペクトルから得られた所定個の比率列を参照し、それらの比率列から評価値を求めてもよい。順位ごとに評価値を演算し、複数の評価値を総合考慮して、安定を判定してもよい。第１判定条件及び第２判定条件として様々な条件を採用し得る。

【0067】

第１トレンドチャート及び第２トレンドチャートの一方又は両方に基づいて、自動ピーク検出条件を自動的に設定し、あるいは、自動ピーク検出結果を事後的に評価してもよい。例えば、マススペクトルが安定していない期間をピークサーチ条件から除外してもよい。換言すれば、マススペクトルが安定している期間内において微小ピークの探索を行ってもよい。また、自動ピーク検出によって検出されたピークが安定期間内に属するか否かによって、検出されたピークを、化合物に由来するピークかノイズに相当するピークかを識別してもよい。

【0068】

図１２には、第１表示例が示されている。画面１６６上には、クロマトグラム１６８、第１トレンドチャート１７０及び第２トレンドチャート１７２が表示されている。それらに含まれる３つの保持時間軸は互いに平行であり、３つの保持時間軸上のスケールは同一である。クロマトグラム１６８上には３つのマーカー１７４Ａ，１７４Ｂ，１７４Ｃが表示されている。３つのマーカー１７４Ａ，１７４Ｂ，１７４Ｃは、総合判定の結果、導き出された３つの安定期間を示すものである。各マーカー１７４Ａ，１７４Ｂ，１７４Ｃは、着色エリアであり、あるいは、ボックス状の図形である。

【0069】

第１トレンドチャート１７０上には、３つのマーカー１７６Ａ，１７６Ｂ，１７６Ｃが表示されている。３つのマーカー１７６Ａ，１７６Ｂ，１７６Ｃは、第１判定条件により導き出された３つの安定期間を示すものである。各マーカー１７６Ａ，１７６Ｂ，１７６Ｃは、着色エリアであり、あるいは、ボックス状の図形である。

【0070】

第２トレンドチャート１７２上にも、３つのマーカー１７８Ａ，１７８Ｂ，１７８Ｃが表示されている。３つのマーカー１７８Ａ，１７８Ｂ，１７８Ｃは、第２判定条件により導き出された３つの安定期間を示すものである。各マーカー１７８Ａ，１７８Ｂ，１７８Ｃは、着色エリアであり、あるいは、ボックス状の図形である。個々のマーカーの表示態様としては様々なものを採用し得る。

【0071】

クロマトグラム１６８の解析に際して、３つのマーカー１７４Ａ，１７４Ｂ，１７４Ｃを参照することにより、各ピークが化合物に由来するピークであるのか否かを容易に評価することが可能である。その際、第１トレンドチャート１７０の内容及び第２トレンドチャート１７２の内容を参照してもよいし、また、それらに重畳表示されたマーカー１７６Ａ，１７６Ｂ，１７６Ｃ，１７８Ａ，１７８Ｂ，１７８Ｃを参照してもよい。

【0072】

なお、クロマトグラム１６８、及び、安定情報としてのマーカー１７４Ａ，１７４Ｂ，１７４Ｃだけを表示してもよい。

【0073】

図１３には、第２表示例が示されている。画面１８０内には、クロマトグラム１８２、第１トレンドチャート１８４及び第２トレンドチャート１８６に加えて、第１判定条件設定欄１９０及び第２判定条件設定欄１９２が含まれる。第１判定条件設定欄１９０を利用して第１判定条件がユーザーにより設定される。第２判定条件設定欄１９２を利用して第２判定条件がユーザーにより設定される。

【0074】

更に、画面１８０内に、マススペクトル１８８が表示されてもよい。その場合に、指定された時相のマススペクトルが表示されてもよいし、ノイズ成分除去後のマススペクトルが表示されてもよい。例えば、ピーク内の主抽出位置に対応するマススペクトルをＡと表現し、ピーク両外側の２つの副抽出位置（２つのノイズ成分位置）に対応する２つのマススペクトルをＢ及びＣと表現した場合、（Ａ－（Ｂ＋Ｃ）／２）で算出されるマススペクトルがノイズ成分除去後のマススペクトルとされる。

【0075】

図１４の上段には、保持時間全体に対応するクロマトグラム１９４が示されている。図１４の下段には、クロマトグラム１９４の一部１９６の拡大表示が示されている。拡大された一部１９８において、個々のピーク２００－２０６ごとに安定期間が特定されている。例えば、ピーク２００についての安定期間の始期が符号２１０で示され、その終期が符号２１２で示されている。各安定期間内に対して、縦軸方向に個別的に拡大率２０８が定められている。例えば、符号２０６で示す部分の強度は、クロマトグラム１９４を参照すればわかるように、非常に小さいが、縦軸方向の拡大率を大きくすることにより、その形態を認識することが可能となっている。安定期間内に属するピークは仮にその強度が小さいものであっても化合物由来のピークである可能性が高い。

【0076】

なお、複数の安定期間の間に非安定期間（ノイズ期間）が存在する場合、非安定期間を除外して複数の安定期間を連結して表示してもよい。その場合においては、保持時間軸上において省略した期間を特定する情報を付加した方がよい。

【0077】

実施形態に係るトレンドチャートの利用形態として様々なものが考えられる。その内で幾つかを以下に説明する。

【0078】

例えば、図１５及び図１６に示されるように、トレンドチャートに基づいて、クロマトグラムに含まれる複数のピーク間での類似関係（又は同一関係）を判定してもよい。図１５の上段には、クロマトグラム２１４が示され、図１５の下段には、第１トレンドチャート２１６が示されている。図示の例では、第１トレンドチャート２１６において、安定期間２１８を代表するｍ／ｚ列と、安定期間２２０を代表するｍ／ｚ列とが一致している（符号２２２を参照）。クロマトグラム２１４においては、２つの安定期間を示すマーカー２１８Ａ，２２０Ａが表示されている。マーカー２１８Ａ、２２０Ａに対応付けて、それらの近傍に、同一内容を有する２つのラベル２２６，２２８が表示される。

【0079】

２つのｍ／ｚ列の一致（符号２２２を参照）に従って、同一ラベルを付与する処理が実施される（符号２２４を参照）。これにより、ユーザーにおいて、２つのピーク間でのマススペクトルの類似性を認識することが可能となる。例えば、上記表示内容から異性体関係を特定し得る。

【0080】

ラベルとして、番号、文字、記号、色相等の様々な表示要素を用い得る。最初からラベルセットを表示し続けると、ラベルセットが目障りとなるのであれば、所定条件が満たされた時点でラベルセットを一時的に表示してもよい。例えば、安定期間やピークがポインタによって指定され又は選択された場合に、ラベルセットが一時的に表示されるようにしてもよい。なお、２つのｍ／ｚ列の間の類似関係又は同一関係の判定に際しては様々な条件を設定し得る。順序までを含む完全一致を求めてもよいし、順序を問わず組み合わせの一致を求めてもよい。ｍ／ｚ列の比較に加えて、比率列の比較を行ってもよい。

【0081】

図１６を用いて具体的に説明する。第１トレンドチャート及び第２トレンドチャートに基づいて、符号２３０で示すように、複数の安定期間が判定される。それらの安定期間については、符号２３２で示すように、複数の安定期間を代表する複数のｍ／ｚ列が特定される。符号２３４で示すように、複数のｍ／ｚ列の中に類似関係を有するｍ／ｚ列群が存在するか否かが調査される。類似関係を有するｍ／ｚ列群が存在する場合、符号２３６で示すように、それらに対応する複数の安定期間又は複数のピークに対して同一内容を有する複数のラベルが付与される。

【0082】

類似関係の調査に際しては、複数の安定期間を代表する複数の比率列を考慮してもよい。例えば、複数のｍ／ｚ列に基づく第１段階の調査で特定された複数の安定期間が、更に比率列の類似性の観点から、絞り込まれてもよい。図示の例では、符号２３８で示すように、第１段階の調査で特定された複数の安定期間を代表する複数の比率列が特定され、符号２４０で示すように、それらについて個別的に評価値が演算される。つまり評価値群が演算される。評価値群に基づいて、第１段階の調査で特定された複数の安定期間の内で、類似関係を認める複数の安定期間が判定される。つまり、第１段階の調査結果に対して第２段階の調査が適用される。類似関係が認められた複数の安定期間に対して同一内容を有する複数のラベルが付与される。それを参照することにより、クロマトグラム中に含まれる特定の複数のピークの相互関係を認識することが可能となる。

【0083】

上記評価値として、比較対象となった２つの比率列の間で相関係数を演算してもよい。その場合、相関係数が優良条件を満たす場合に類似を判定してもよい。また、上記評価値として、比較対象となった２つの比率列の間で順位ごとに差分を演算してもよい。その場合、各差分が閾値以下である場合に類似を判定してもよい。

【0084】

第１トレンドチャート（又は第２トレンドチャート）に基づいて、マススペクトル抽出位置が補正されてもよい。これについて以下に説明する。

【0085】

なお、ピークを代表するマススペクトルを抽出する位置を主抽出位置と表現し、ピークの両側に定められるノイズ成分を抽出する２つの抽出位置をそれぞれ副抽出位置と表現する。主抽出において上記マススペクトルＡが抽出され、２つの副抽出位置において上記マススペクトルＢ，Ｃが抽出される。

【0086】

図１７の上段には、クロマトグラム２４４が示されており、図１７の下段には、第１トレンドチャート２４６が示されている。一般に、クロマトグラム２４４に含まれるピーク２４８のピークトップの位置２５２が主抽出位置として特定され、その位置２５２に基づいてマススペクトル列からマススペクトルＡが抽出される。また、ピーク２４８についてはピーク期間２５０が特定され、その両端の位置２５４，２５６がそれぞれ副抽出位置とされる。２つの副抽出位置に基づいてマススペクトル列からマススペクトルＢ，Ｃが抽出される。以下のように、それらの位置２５２，２５４，２５６の内の一部又は全部が安定期間に基づいて補正され得る。

【0087】

具体的に説明すると、第１トレンドチャートに１つの安定期間のみが含まれる場合、安定期間内に主抽出位置が定められる。同様に、安定期間を避けて２つの副抽出位置が定められる。第１トレンドチャートに複数の安定期間が含まれる場合も基本的に同様である。

【0088】

例えば、図１７中の第１トレンドチャート２４６においては、２つの安定期間２５８，２６０が判定されている。ピーク期間２５０と２つの安定期間２５８，２６０の重複期間が符号２６２，２６４で示されている。２つの安定期間２５８，２６０の内で、大きい方の安定期間２６０が優先され、その安定期間２６０の中で、最も大きなＴＩＣとなる位置２６５が特定され、その位置２６５が、補正された主抽出位置とされる。ピークトップの位置２５２ではなく、そこからシフトした位置２６５が主抽出位置とされる。

【0089】

２つの副抽出位置は、通常、位置２５４及び位置２５６であるが、それらが補正されてもよい。例えば、安定期間２６２の手前側（左側）の位置が第１の副抽出位置とされてもよいし、安定期間２６４の奥側（右側）の位置が第２の副抽出位置とされてもよい。

【0090】

図１８に示す例では、クロマトグラム２６６において、ピーク期間が符号２７０で示されている。その両端の位置が符号２７４，２７６で示されている。第１トレンドチャート２６８において、安定期間が符号２７２で示されている。図示の例においては、安定期間２７２の手前側の位置２７４Ａが第１の副抽出位置とされ、安定期間２７２の奥側の位置２７６Ａが第２の副抽出位置とされる。

【0091】

クロマトグラムにおいて、複数のサブピークが重なって１つのピークとして見えている場合に、代表値列又は安定期間を利用して、複数のサブピークを分離してもよい。複数のサンプルから得られた複数のクロマトグラムを相互に比較する場合において、複数のトレンドチャートから求められる複数の代表値列を比較してもよい。その比較結果から、複数のクロマトグラム間において、同一性又は類似性を有する複数のピークが特定されてもよい。

【0092】

上記実施形態においては、代表ピークから求められるｍ／ｚを第１代表値とし、また、代表ピーク間の強度比を第２代表値としたが、代表ピークから求められる他の情報を代表値としてもよい。例えば、代表ピーク強度の時間的変化、代表ピーク間の面積比、等を代表値としてもよい。

【符号の説明】

【0093】

１０ 質量分析システム、１２ 測定部、１４ 情報処理装置、１６ ガスクロマトグラフ、１８ 質量分析装置、２８ スペクトル生成部、３４ ＴＩＣＣ生成部、３６ 第１トレンドチャート生成部、３８ 第２トレンドチャート生成部、４０ チャート解析部。

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】

【図15】

【図16】

【図17】

【図18】