(57)【要約】  自動波形処理部１４３によって、クロマトグラムのピークが自動で検出された後、そのピークに対してユーザが手動操作を行うと、手動波形処理部１４４は、そのピークの情報を変更する。時間差算出部１４７は、２つの解析済みデータにおける同一成分のピーク間の溶出時間差を算出する。ピーク選択部１４６は、対象データについて、解析済みデータのピークと同一成分のピークを、分析日時、すなわち、クロマトグラムデータの生成日時の順序に基づいて選択する。ピーク情報反映部１４８は、溶出時間差及び分析日時に基づいて、解析済みデータにおけるピークのピーク情報１３４の変更を、対象データにおいて選択されたピークのピーク情報１３４に反映させる。そのため、対象データにおけるピークのピーク情報１３４を、手動操作によらずに自動的に変更できる。また、対象データにおけるピークのピーク情報１３４の変更を精度よく行うことができる。
2017179154.tif 000002
【特許請求の範囲】
【請求項１】
同一成分を含むサンプルに対するクロマトグラフィー分析を繰り返し行うことにより得られた、分析日時が異なる３つ以上のクロマトグラムデータに対する処理を行うデータ処理装置であって、
前記３つ以上のクロマトグラムデータのそれぞれについて、ピークを自動で検出する自動波形処理部と、
ユーザによる手動操作に基づいて、前記自動波形処理部で検出されたピークのピーク情報を変更する手動波形処理部と、
前記手動波形処理部により同一成分のピークの前記ピーク情報が変更された２つのクロマトグラムデータを解析済みデータとして、当該２つの解析済みデータにおける同一成分のピーク間の溶出時間差を算出する時間差算出部と、
前記３つ以上のクロマトグラムデータのうち前記２つの解析済みデータを除く対象データについて、前記手動波形処理部により前記ピーク情報が変更されたピークと同一成分のピークを分析日時に基づいて選択するピーク選択部と、
前記時間差算出部により算出された溶出時間差及び分析日時に基づいて、前記解析済みデータにおけるピークの前記ピーク情報の変更を、前記ピーク選択部により選択された対象データにおけるピークの前記ピーク情報に反映させるピーク情報反映部とを備えることを特徴とするデータ処理装置。
【請求項２】
前記解析済みデータ及び前記対象データを分析日時の順序に従って表示させる表示処理部をさらに備えることを特徴とする請求項１に記載のデータ処理装置。
【請求項３】
前記表示処理部は、前記解析済みデータ及び前記対象データにおける同一成分のピークに対応付けて、少なくとも分析日時が直前のデータのピーク位置及び分析日時が直後のデータのピーク位置を表示させることを特徴とする請求項２に記載のデータ処理装置。
【請求項４】
前記解析済みデータにおける時間範囲を指定する時間範囲指定部をさらに備え、
前記ピーク選択部は、前記時間範囲指定部により指定された時間範囲及び分析日時に基づいて、前記手動波形処理部により前記ピーク情報が変更されたピークと同一成分のピークを選択することを特徴とする請求項１に記載のデータ処理装置。
【請求項５】
同一成分を含むサンプルに対するクロマトグラフィー分析を繰り返し行うことにより得られた、分析日時が異なる３つ以上のクロマトグラムデータに対する処理を行うデータ処理方法であって、
前記３つ以上のクロマトグラムデータのそれぞれについて、ピークを自動で検出する自動波形処理ステップと、
ユーザによる手動操作に基づいて、前記自動波形処理ステップで検出されたピークのピーク情報を変更する手動波形処理ステップと、
前記手動波形処理ステップにより同一成分のピークの前記ピーク情報が変更された２つのクロマトグラムデータを解析済みデータとして、当該２つの解析済みデータにおける同一成分のピーク間の溶出時間差を算出する時間差算出ステップと、
前記３つ以上のクロマトグラムデータのうち前記２つの解析済みデータを除く対象データについて、前記手動波形処理ステップにより前記ピーク情報が変更されたピークと同一成分のピークを分析日時に基づいて選択するピーク選択ステップと、
前記時間差算出ステップにより算出された溶出時間差及び分析日時に基づいて、前記解析済みデータにおけるピークの前記ピーク情報の変更を、前記ピーク選択ステップにより選択された対象データにおけるピークの前記ピーク情報に反映させるピーク情報反映ステップとを含むことを特徴とするデータ処理方法。
【請求項６】
前記解析済みデータ及び前記対象データを分析日時の順序に従って表示させる表示処理ステップをさらに含むことを特徴とする請求項５に記載のデータ処理方法。
【請求項７】
前記表示処理ステップでは、前記解析済みデータ及び前記対象データにおける同一成分のピークに対応付けて、少なくとも分析日時が直前のデータのピーク位置及び分析日時が直後のデータのピーク位置を表示させることを特徴とする請求項６に記載のデータ処理方法。
【請求項８】
前記解析済みデータにおける時間範囲を指定する時間範囲指定ステップをさらに含み、
前記ピーク選択ステップでは、前記時間範囲指定ステップにより指定された時間範囲及び分析日時に基づいて、前記手動波形処理ステップにより前記ピーク情報が変更されたピークと同一成分のピークを選択することを特徴とする請求項５に記載のデータ処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、クロマトグラムデータに対して自動で波形処理を行うデータ処理装置及びデータ処理方法に関するものである。
【背景技術】
【０００２】
従来より、クロマトグラフィー分析を行うに際して、自動で波形処理を行うデータ処理装置が利用されている。【０００３】
このようなデータ処理装置は、波形処理として、例えば、クロマトグラムデータにおけるピークを検出する処理などを行う。具体的には、データ処理装置は、クロマトグラムデータにおけるピークの開始点及び終了点を検出することで、その開始点から終了点までの波形をピークとして認識（検出）する。そして、データ処理装置は、ピークの開始点と終了点とを結ぶベースラインを作成し、このベースラインとピークとで囲まれる領域の面積を算出する。データ処理装置では、このような処理が自動で行われている（例えば、下記特許文献１参照）。【０００４】
このようなデータ処理装置を用いる場合において、データ処理装置による処理結果がユーザが所望する内容とは異なる場合がある。例えば、ピークの開始点及び終了点が本来の位置からずれて検出される場合がある。この場合には、ユーザは、表示部などに表示されるクロマトグラムを確認しながら、手動操作でピークの開始点及び終了点を変更（補正）する。そして、ユーザは、変更した後の正確なピークに基づいて分析を行う。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００４−２７１４２２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
上記のような従来のデータ処理装置では、ユーザ作業が煩雑化するおそれがあった。具体的には、分析動作を連続して複数回行うような場合には、クロマトグラムにおけるピークの開始点及び終了点が適切な位置からずれて検出されることが繰り返し発生してしまう。そして、ユーザは、クロマトグラムごとに手動操作でピークの開始点及び終了点を変更する必要がある。【０００７】
このような場合において、あるクロマトグラムに対する手動操作によるピーク（ピーク情報）の変更内容を、他のクロマトグラムに対してそのまま用いることで、他のクロマトグラムに対するユーザの手動操作を省くことが検討される。例えば、複数のクロマトグラムのそれぞれから同一成分のピークを検出する場合において、手動操作によって補正を行った後のクロマトグラムにおけるピークの開始点及び終了点の時間を、そのまま、他のクロマトグラムにおいてピークの開始点及び終了点の時間としてピークを検出すれば、他のクロマトグラムに対するユーザの手動操作を省くことが可能である。【０００８】
しかし、分析動作を繰り返し行うような場合（クロマトグラムデータを連続して生成する場合）には、前回分析時の残留物がカラム内に存在することなどにより、同一成分のピークであっても、クロマトグラムごとでピークの溶出時間に個別のずれが生じてしまう。その結果、同一成分のピークであっても、クロマトグラムごとでピークの開始点及び終了点の時間が異なり、ピークを精度よく検出できないという不具合が生じる。【０００９】
本発明は、上記実情に鑑みてなされたものであり、自動で波形処理されたデータを補正する作業を効率化でき、かつ、当該補正を精度よく行うことができるデータ処理装置及びデータ処理方法を提供することを目的とする。
【課題を解決するための手段】
【００１０】
（１）本発明に係るデータ処理装置は、同一成分を含むサンプルに対するクロマトグラフィー分析を繰り返し行うことにより得られた、分析日時が異なる３つ以上のクロマトグラムデータに対する処理を行う。前記データ処理装置は、自動波形処理部と、手動波形処理部と、時間差算出部と、ピーク選択部と、ピーク情報反映部とを備える。前記自動波形処理部は、前記３つ以上のクロマトグラムデータのそれぞれについて、ピークを自動で検出する。前記手動波形処理部は、ユーザによる手動操作に基づいて、前記自動波形処理部で検出されたピークのピーク情報を変更する。前記時間差算出部は、前記手動波形処理部により同一成分のピークの前記ピーク情報が変更された２つのクロマトグラムデータを解析済みデータとして、当該２つの解析済みデータにおける同一成分のピーク間の溶出時間差を算出する。前記ピーク選択部は、前記３つ以上のクロマトグラムデータのうち前記２つの解析済みデータを除く対象データについて、前記手動波形処理部により前記ピーク情報が変更されたピークと同一成分のピークを分析日時に基づいて選択する。前記ピーク情報反映部は、前記時間差算出部により算出された溶出時間差及び分析日時に基づいて、前記解析済みデータにおけるピークの前記ピーク情報の変更を、前記ピーク選択部により選択された対象データにおけるピークの前記ピーク情報に反映させる。【００１１】
このような構成によれば、自動波形処理部によってクロマトグラムのピークが自動で検出された後、そのピークに対してユーザが手動操作を行うと、手動波形処理部は、ピークのピーク情報を変更する。そして、ピーク情報反映部は、その変更されたピーク情報の内容を、他のクロマトグラムデータ（対象データ）におけるピークのピーク情報に反映させる。
そのため、対象データにおけるピークのピーク情報を、手動操作によらずに自動的に変更できる。
その結果、対象データを補正する作業を効率化できる。【００１２】
また、時間差算出部は、２つの解析済みデータにおける同一成分のピーク間の溶出時間差を算出する。ピーク選択部は、対象データについて、解析済みデータのピークと同一成分のピークを、分析日時（分析日時の順序）、すなわち、クロマトグラムデータの生成日時の順序に基づいて選択する。ピーク情報反映部は、これらの溶出時間差及び分析日時（分析日時の順序）に基づいて、解析済みデータにおけるピークのピーク情報の変更を、対象データにおいて選択されたピークのピーク情報に反映させる。【００１３】
すなわち、データ処理装置では、まず、対象データについて、解析済みデータのピークと同一成分のピークを正確に選択する。そして、データ処理装置では、２つの解析済みデータにおける同一成分のピーク間の溶出時間差、及び、分析日時（分析日時の順序）に基づくことにより、クロマトグラムごとでのピークの溶出時間のずれを考慮して、解析済みデータにおけるピーク情報の変更内容を、その選択されたピークのピーク情報に反映させる。
そのため、対象データにおけるピークのピーク情報の変更（対象データの補正）を精度よく行うことができる。
このように、本発明に係るデータ処理装置によれば、自動で波形処理されたデータを補正する作業を効率化でき、かつ、当該補正を精度よく行うことができる。【００１４】
（２）また、前記データ処理装置は、表示処理部をさらに備えてもよい。前記表示処理部は、前記解析済みデータ及び前記対象データを分析日時の順序に従って表示させる。
このような構成によれば、ユーザは、解析済みデータ及び対象データにおける同一成分のピークのピーク位置を、分析日時の順序と照らし合わせて確認できる。
そのため、分析日時の順序と、ピーク位置のずれとの関係を容易に確認できる。【００１５】
（３）また、前記表示処理部は、前記解析済みデータ及び前記対象データにおける同一成分のピークに対応付けて、少なくとも分析日時が直前のデータのピーク位置及び分析日時が直後のデータのピーク位置を表示させてもよい。【００１６】
このような構成によれば、ユーザは、同一成分のピークに関して、あるデータにおけるピーク位置を、分析日時が直前及び直後のデータにおけるピーク位置と比較して確認できる。【００１７】
（４）また、前記データ処理装置は、時間範囲指定部をさらに備えてもよい。前記時間範囲指定部は、前記解析済みデータにおける時間範囲を指定する。前記ピーク選択部は、前記時間範囲指定部により指定された時間範囲及び分析日時に基づいて、前記手動波形処理部により前記ピーク情報が変更されたピークと同一成分のピークを選択してもよい。
このような構成によれば、時間範囲設定部によって時間範囲を指定することにより、特定の成分に限定してピーク情報を変更できる。
そのため、全体の処理時間を短くできる。【００１８】
（５）本発明に係るデータ処理方法は、同一成分を含むサンプルに対するクロマトグラフィー分析を繰り返し行うことにより得られた、分析日時が異なる３つ以上のクロマトグラムデータに対する処理を行う。前記データ処理方法は、自動波形処理ステップと、手動波形処理ステップと、時間差算出ステップと、ピーク選択ステップと、ピーク情報反映ステップとを含む。前記自動波形処理ステップでは、前記３つ以上のクロマトグラムデータのそれぞれについて、ピークを自動で検出する。前記手動波形処理ステップでは、ユーザによる手動操作に基づいて、前記自動波形処理ステップで検出されたピークのピーク情報を変更する。前記時間差算出ステップでは、前記手動波形処理ステップにより同一成分のピークの前記ピーク情報が変更された２つのクロマトグラムデータを解析済みデータとして、当該２つの解析済みデータにおける同一成分のピーク間の溶出時間差を算出する。前記ピーク選択ステップでは、前記３つ以上のクロマトグラムデータのうち前記２つの解析済みデータを除く対象データについて、前記手動波形処理ステップにより前記ピーク情報が変更されたピークと同一成分のピークを分析日時に基づいて選択する。前記ピーク情報反映ステップでは、前記時間差算出ステップにより算出された溶出時間差及び分析日時に基づいて、前記解析済みデータにおけるピークの前記ピーク情報の変更を、前記ピーク選択ステップにより選択された対象データにおけるピークの前記ピーク情報に反映させる。【００１９】
（６）また、前記データ処理方法は、表示処理ステップをさらに含んでもよい。前記表示処理ステップは、前記解析済みデータ及び前記対象データを分析日時の順序に従って表示させる。【００２０】
（７）また、前記表示処理ステップでは、前記解析済みデータ及び前記対象データにおける同一成分のピークに対応付けて、少なくとも分析日時が直前のデータのピーク位置及び分析日時が直後のデータのピーク位置を表示させてもよい。【００２１】
（８）また、前記データ処理方法は、時間範囲指定ステップをさらに含んでもよい。前記時間範囲指定ステップでは、前記解析済みデータにおける時間範囲を指定する。前記ピーク選択ステップでは、前記時間範囲指定ステップにより指定された時間範囲及び分析日時に基づいて、前記手動波形処理ステップにより前記ピーク情報が変更されたピークと同一成分のピークを選択する。
【発明の効果】
【００２２】
本発明によれば、対象データのピークのピーク情報を、手動操作によらずに自動的に変更できる。また、２つの解析済みデータにおける同一成分のピーク間の溶出時間差、及び、分析日時に基づくことにより、連続して分析を行う場合に生じるクロマトグラムごとでのピークの溶出時間のずれを考慮して、対象データのピークのピーク情報を変更できる。そのため、自動で波形処理されたデータを補正する作業を効率化でき、かつ、当該補正を精度よく行うことができる。
【図面の簡単な説明】
【００２３】
【図１】本発明の一実施形態に係るデータ処理装置、及び、当該データ処理装置とともに用いられるガスクロマトグラフの構成例を示した概略図である。
【図２】図１のデータ処理装置の制御部、及び、その周辺の部材の具体的構成を示したブロック図である。
【図３】制御部による処理の一例を示したフローチャートである。
【図４】表示画面の表示態様の一例であって、自動波形処理が行われたクロマトグラムを表示する際の表示例を示した概略図である。
【図５】表示画面の表示態様の一例であって、クロマトグラムに対して手動波形処理が行われる際の表示例を示した概略図である。
【図６】表示画面の表示態様の一例であって、解析済みデータ及び対象データを表示する際の表示例を示した概略図である。
【図７】表示画面の表示態様の一例であって、解析済みデータ及びピーク情報が反映された対象データを表示する際の表示例を示した概略図である。
【発明を実施するための形態】
【００２４】
１．ガスクロマトグラフの全体構成
図１は、本発明の一実施形態に係るデータ処理装置１、及び、当該データ処理装置１とともに用いられるガスクロマトグラフ２の構成例を示した概略図である。
【００２５】
データ処理装置１は、ガスクロマトグラフ２で得られるデータに対して各種のデータ処理を行うための装置である。データ処理装置１は、例えば、コンピュータにより構成される。データ処理装置１の詳しい構成については、後述する。
【００２６】
ガスクロマトグラフ２は、キャリアガスとともに試料ガスをカラム３内に供給することにより分析を行うためのものであり、上記カラム３以外に、カラムオーブン４、試料導入部５及び検出器６などを備えている。
カラム３は、例えば、キャピラリカラムからなる。カラム３は、ヒータ及びファンなど（いずれも図示せず）とともにカラムオーブン４内に収容されている。
カラムオーブン４は、カラム３を加熱するためのものであり、分析時にはヒータ及びファンを適宜駆動させる。
【００２７】
試料導入部５は、カラム３内にキャリアガス及び試料ガスを導入するためのものであり、その内部に試料気化室（図示せず）が形成されている。この試料気化室には、液体試料が注入され、試料気化室内で気化された試料が、キャリアガスとともにカラム３内に導入される。また、試料気化室には、ガス供給流路７及びスプリット流路８が連通している。
ガス供給流路７は、試料導入部５の試料気化室内にキャリアガスを供給するための流路である。
【００２８】
スプリット流路８は、スプリット導入法によりカラム３内にキャリアガス及び試料ガスを導入する際に、試料気化室内のガス（キャリアガス及び試料ガスの混合ガス）の一部を所定のスプリット比で外部に排出するための流路である。
【００２９】
検出器６は、例えば、水素炎イオン化検出器（ＦＩＤ）や、炎光光度検出器（ＦＰＤ）により構成される。検出器６は、カラム３から導入されるキャリアガスに含まれる各試料成分を順次検出する。
【００３０】
ガスクロマトグラフ２において試料を測定する際は、分析対象となる試料が試料導入部５に注入される。試料は、試料気化室において気化される。また、試料導入部５の試料気化室には、ガス供給流路７を介してキャリアガスが供給される。
【００３１】
試料気化室内で気化された試料は、キャリアガスとともにカラム３内に導入される。試料に含まれる各試料成分は、カラム３内を通過する過程で分離されて、検出器６に順次導入される。
【００３２】
そして、検出器６において、カラム３から導入されるキャリアガスに含まれる各試料成分が順次検出される。また、検出器６からの検出信号が、データ処理装置１に入力される。
データ処理装置１では、入力された信号に基づいてクロマトグラムが生成される。そして、ユーザは、得られたクロマトグラムを確認して、各種分析を行う。
【００３３】
２．制御部及び周辺の部材の具体的構成
図２は、データ処理装置１の制御部１４、及び、その周辺の部材の具体的構成を示したブロック図である。
データ処理装置１は、操作部１１と、表示部１２と、記憶部１３と、制御部１４とを備えている。
【００３４】
操作部１１は、例えば、キーボード及びマウスなどにより構成される。ユーザは、操作部１１を操作することにより、分析条件などの各種情報を制御部１４に入力することができる。
表示部１２は、例えば、液晶表示器などにより構成される。表示部１２には、制御部１４の制御により、分析結果などの各種情報が表示される。
【００３５】
記憶部１３は、例えば、ＲＯＭ（Read Only Memory）、ＲＡＭ（Random Access Memory）及びハードディスクなどにより構成されている。記憶部１３は、複数の検出データ１３１を記憶している。各検出データ１３１は、データ処理装置１で生成されるクロマトグラムに関するデータである。各検出データ１３１には、クロマトグラムデータ１３２と、分析日時情報１３３と、ピーク情報１３４とが含まれている。クロマトグラムデータ１３２と、分析日時情報１３３と、ピーク情報１３４とは、互いに対応付けられている。
【００３６】
クロマトグラムデータ１３２は、後述するデータ処理部１４１によって、検出器６からの検出信号に基づいて生成されるクロマトグラムのデータ（クロマトグラムデータ）である。
分析日時情報１３３は、クロマトグラムデータ１３２が生成された日時、具体的には、検出器６によって検出信号が出力された日時の情報である。
【００３７】
ピーク情報１３４は、クロマトグラムデータ１３２によって表されるクロマトグラムにおけるピークに関する情報である。具体的には、ピーク情報１３４には、クロマトグラムにおけるピークの開始点、ピークの終了点、及び、ピーク位置の情報が含まれる。なお、ピークの開始点の情報とは、ピークの開始点の時間の情報であり、ピークの終了点の情報とは、ピークの終了点の時間の情報であり、ピーク位置の情報とは、ピークが一番高い点（ピークトップ）の時間の情報である。
【００３８】
制御部１４は、例えば、ＣＰＵ（Central Processing Unit）を含む構成である。制御部１４は、操作部１１、表示部１２及び検出器６との間で、電気信号の入力又は出力が可能である。制御部１４は、必要に応じて、記憶部１３に情報を格納し、また、記憶部１３に記憶されている情報を読み出す。制御部１４は、ＣＰＵがプログラムを実行することにより、データ処理部１４１、表示処理部１４２、自動波形処理部１４３、手動波形処理部１４４、時間範囲指定部１４５、ピーク選択部１４６、時間差算出部１４７及びピーク情報反映部１４８などとして機能する。
【００３９】
データ処理部１４１は、検出器６からの検出信号に基づいてクロマトグラムを得る。データ処理部１４１で生成されたクロマトグラムのデータは、クロマトグラムデータ１３２として記憶部１３に格納される。また、検出器６から検出信号が出力された日時の情報は、分析日時情報１３３としてクロマトグラムデータ１３２に対応付けて記憶部１３に格納される。
【００４０】
表示処理部１４２は、生成されたクロマトグラムを表示部１２にリアルタイムで表示させる処理を行う。すなわち、データ処理部１４１でクロマトグラムが生成されると、上記したように、そのデータがクロマトグラムデータ１３２として記憶部１３に格納されるとともに、リアルタイムで表示部１２にクロマトグラムが表示される。
【００４１】
自動波形処理部１４３は、記憶部１３に記憶されているクロマトグラムデータ１３２で表されるクロマトグラムにおけるピークを自動で検出する処理を行う。具体的には、記憶部１３は、図示しないピーク検出用のパラメータを複数記憶している。そして、自動波形処理部１４３は、クロマトグラムにおいて、このパラメータを満たす部分の波形をピークとして自動で検出する処理を行う。パラメータとしては、ピークの開始点を検出するためのパラメータ、ピークの終了点を検出するためのパラメータ、及び、ピーク位置（ピークトップの位置）を検出するためのパラメータが挙げられる。すなわち、自動波形処理部１４３は、クロマトグラムにおける、ピークの開始点、ピークの終了点、及び、ピーク位置を自動で検出する。また、自動波形処理部１４３によって検出されたピークの情報（ピークの開始点、ピークの終了点、及び、ピーク位置の情報）は、ピーク情報１３４として記憶部１３に格納される。
【００４２】
手動波形処理部１４４は、詳しくは後述するが、ユーザによる操作部１１の操作に基づいて、クロマトグラムのピークに関する変更（補正）を受け付ける。また、手動波形処理部１４４は、この受け付けた内容に基づいて、記憶部１３に記憶されているピーク情報１３４を書き換える。
【００４３】
表示処理部１４２は、上記した処理に加えて、自動波形処理部１４３によって検出されたクロマトグラムのピークを、表示部１２に表示させる処理を行う。また、表示処理部１４２は、手動波形処理部１４４が受け付けた内容に基づいて、表示部１２に表示されているクロマトグラムのピークを補正して表示する処理を行う。
【００４４】
時間範囲指定部１４５は、詳しくは後述するが、ユーザによる操作部１１の操作に基づいて、ピーク情報１３４を自動で補正（反映処理）する際の時間範囲の指定を受け付ける。
【００４５】
ピーク選択部１４６は、記憶部１３に記憶されている分析日時情報１３３、及び、時間範囲指定部１４５が受け付けた時間範囲に基づいて、クロマトグラムデータ１３２で表されるクロマトグラムにおいて、自動で補正を行う対象となるピークを選択する。
【００４６】
時間差算出部１４７は、詳しくは後述するが、手動波形処理部１４４によって内容が書き換えられた後の２つのピーク情報１３４に基づいて、ピークの溶出時間差を算出する。
【００４７】
ピーク情報反映部１４８は、詳しくは後述するが、記憶部１３に記憶されている分析日時情報１３３、及び、時間差算出部１４７が算出するピークの溶出時間差に基づいて、ピーク選択部１４６が選択するクロマトグラムのピーク（補正対象となるピーク）のピーク情報１３４に対して、手動波形処理部１４４によって変更されたピーク情報１３４の内容を反映する。
【００４８】
表示処理部１４２は、上記した処理に加えて、ピーク情報反映部１４７によって変更された後のピーク情報１３４に基づいて、クロマトグラムのピークを補正して表示する処理を行う。
【００４９】
３．制御部による制御動作、及び、表示画面の表示態様
（１）自動波形処理
図３は、制御部１４による処理の一例を示したフローチャートである。
データ処理装置１及びガスクロマトグラフ２において分析動作が開始されると、データ処理部１４１（図２参照）は、検出器６からの検出信号に基づいてクロマトグラムを生成する。そして、図示しないが、表示部１２には、クロマトグラムがリアルタイムで表示される。また、生成されたクロマトグラムのデータは、検出データ１３１（クロマトグラムデータ１３２及び分析日時情報１３３）として記憶部１３に格納される。なお、この例では、同一成分を含むサンプルに対して、分析日時が異なる分析動作（クロマトグラフィー分析）が４回行われたとして、説明する。
【００５０】
そして、ユーザによって操作部１１が操作されて、自動波形処理のための操作が行われると、自動波形処理部１４３は、記憶部１３に記憶されている４つの検出データ１３１（クロマトグラムデータ１３２及び分析日時情報１３３）を読み出す。また、自動波形処理部１４３は、記憶部１３から図示しないパラメータを読出し、そのパラメータを満たす波形の部分を、各クロマトグラムデータ１３２で表されるクロマトグラムにおけるピークとして検出する（ステップＳ１０１：自動波形処理ステップ）。そして、表示処理部１４２は、自動波形処理部１４３によって検出されるクロマトグラムのピークを、表示画面１２１に表示させる。
【００５１】
図４は、表示画面１２１の表示態様の一例であって、自動波形処理が行われたクロマトグラムを表示する際の表示例を示した概略図である。
上記した表示部１２（図２参照）には、表示画面１２１が含まれている。
【００５２】
自動波形処理が行われると、表示部１２の表示画面１２１には、複数のクロマトグラム領域１２２が表示される。クロマトグラム領域１２２は、自動波形処理によってピークが検出された後のクロマトグラムを表示するための領域である。この例では、表示画面１２１には、クロマトグラム領域１２２として、第１クロマトグラム領域１２３、第２クロマトグラム領域１２４、第３クロマトグラム領域１２５及び第４クロマトグラム領域１２６が表示されている。なお、第１クロマトグラム領域１２３〜第４クロマトグラム領域１２６では、横軸が時間、縦軸が信号強度を表している。
【００５３】
第１クロマトグラム領域１２３、第２クロマトグラム領域１２４、第３クロマトグラム領域１２５及び第４クロマトグラム領域１２６は、表示部１２において、この順序で、上方から下方に向かって並ぶように配置されている。
【００５４】
この例では、分析日時（データが生成された日時）の順序が最初である分析日時情報１３３に対応づけられるクロマトグラムデータ１３２の内容（クロマトグラムデータ１３２で表されるクロマトグラム）が、第１クロマトグラム領域１２３に表示される。また、分析日時の順序が２番目である分析日時情報１３３に対応づけられるクロマトグラムデータ１３２が、第２クロマトグラム領域１２４に表示される。また、分析日時の順序が３番目である分析日時情報１３３に対応づけられるクロマトグラムデータ１３２が、第３クロマトグラム領域１２５に表示される。また、分析日時の順序が４番目である分析日時情報１３３に対応づけられるクロマトグラムデータ１３２が、第４クロマトグラム領域１２６に表示される。具体的には、第１クロマトグラム領域１２３〜第４クロマトグラム領域１２６のそれぞれでは、クロマトグラムＡ〜Ｄのそれぞれが表示されるとともに、自動波形処理部１４３によって検出されたピークが表示される。
【００５５】
このように、表示画面１２１では、上方側に配置されるクロマトグラム領域１２２に表示されるクロマトグラムほど、その生成日時が前であって、下方側に配置されるクロマトグラム領域１２２に表示されるクロマトグラムほど、その生成日時が後となるように、生成日時の順序に従って（分析日時の順序に従って）クロマトグラムが表示される。
【００５６】
具体的には、第１クロマトグラム領域１２３では、クロマトグラムＡ上において、点Ｅと、点Ｆとを結ぶベースラインＧが表示されることにより、点Ｅがピークの開始点として表され、点Ｆがピークの終了点として表され、ベースラインＧの一端（点Ｅ）から他端（点Ｆ）までの間の波形がピークＰ１として表される。
【００５７】
また、第２クロマトグラム領域１２４では、クロマトグラムＢ上において、点Ｈと、点Ｉとを結ぶベースラインＪが表示されることにより、点Ｈがピークの開始点として表され、点Ｉがピークの終了点として表され、ベースラインＪの一端（点Ｈ）から他端（点Ｉ）までの間の波形がピークＰ２として表される。
【００５８】
また、第３クロマトグラム領域１２５では、クロマトグラムＣ上において、点Ｋと、点Ｌとを結ぶベースラインＭが表示されることにより、点Ｋがピークの開始点として表され、点Ｌがピークの終了点として表され、ベースラインＭの一端（点Ｋ）から他端（点Ｌ）までの間の波形がピークＰ３として表される。
【００５９】
また、第４クロマトグラム領域１２６では、クロマトグラムＤ上において、点Ｎと、点Ｏとを結ぶベースラインＱが表示されることにより、点Ｎがピークの開始点として表され、点Ｏがピークの終了点として表され、ベースラインＱの一端（点Ｎ）から他端（点Ｏ）までの間の波形がピークＰ４として表される。
【００６０】
なお、第１クロマトグラム領域１２３〜第４クロマトグラム領域１２６で表示されるクロマトグラムＡ〜ＤにおけるピークＰ１〜Ｐ４は、同一成分を表している。
【００６１】
また、クロマトグラムＡ〜ＤにおけるピークＰ１〜Ｐ４の情報（ピークの開始点の時間、ピークの終了点の時間、ピーク位置の情報時間）は、ピーク情報１３４として記憶部１３に記憶される。
【００６２】
また、各クロマトグラム領域１２２では、各ピークのピーク位置（ピークトップの位置）が矢印で示されるとともに、表示しているクロマトグラムの生成日時（分析日時）に対して、生成日時が直前となるクロマトグラムにおけるピークのピーク位置、及び、生成日時が直後となるクロマトグラムにおけるピークのピーク位置が矢印で示されている。
【００６３】
具体的には、第１クロマトグラム領域１２３では、クロマトグラムＡのピークＰ１のピーク位置が矢印Ａ１で示され、クロマトグラムＡの生成日時に対して、直後の生成日時であるクロマトグラムＢのピークＰ２のピーク位置が矢印Ａ２で示されている。また、第２クロマトグラム領域１２４では、クロマトグラムＢのピークＰ２のピーク位置が矢印Ｂ２で示され、クロマトグラムＢの生成日時に対して、直前の生成日時であるクロマトグラムＡのピークＰ１のピーク位置が矢印Ｂ１で示され、クロマトグラムＢの生成日時に対して、直後の生成日時であるクロマトグラムＣのピークＰ３のピーク位置が矢印Ｂ３で示されている。また、第３クロマトグラム領域１２５では、クロマトグラムＣのピークＰ３のピーク位置が矢印Ｃ３で示され、クロマトグラムＣの生成日時に対して、直前の生成日時であるクロマトグラムＢのピークＰ２のピーク位置が矢印Ｃ２で示され、クロマトグラムＣの生成日時に対して、直後の生成日時であるクロマトグラムＤのピークＰ４のピーク位置が矢印Ｃ４で示されている。また、第４クロマトグラム領域１２６では、クロマトグラムＤのピークＰ４のピーク位置が矢印Ｄ４で示され、クロマトグラムＤの生成日時に対して、直前の生成日時であるクロマトグラムＣのピークＰ３のピーク位置が矢印Ｄ３で示されている。
【００６４】
クロマトグラム領域１２２では、クロマトグラムＡ〜ＤにおけるピークＰ１〜Ｐ４は、同一成分を表すピークであるが、分析日時が後になるほど（下方側に配置されるクロマトグラム領域１２２に表示されるクロマトグラムほど）、右側、すなわち、遅い時点に表れている。また、クロマトグラム領域１２２では、分析日時が１つ後になるたびに、クロマトグラムにおける同一成分のピークが、ほぼ同一時間だけ遅れるように表れている。
【００６５】
これは、データ処理装置１において、連続してクロマトグラフィー分析を行う場合に、前回分析時の残留物がカラム３内に存在することなどに起因している。そこで、本実施形態では、以下のようにして、各クロマトグラムのピークを自動で補正する。
【００６６】
（２）ピーク情報の反映処理
表示画面１２１において、クロマトグラムＡ〜ＤにおけるピークＰ１〜Ｐ４では、その終了点が、クロマトグラムにおいて傾きが水平になる時点よりも前の時点で検出されている。そのため、ユーザは、自動波形処理が行われたクロマトグラムに対して、手動操作でピークを補正する。このとき、ユーザは、クロマトグラム領域１２２に表示されているクロマトグラムを２つ選択し、その選択した各クロマトグラムに対して手動操作でピークの補正を行う（図３のステップＳ１０２でＹＥＳ）。この例では、ユーザは、第１クロマトグラム領域１２３に表示されるクロマトグラムＡ、及び、第３クロマトグラム領域１２５に表示されるクロマトグラムＣを、手動操作で補正する。
具体的には、ユーザは、手動操作でクロマトグラムのピークの補正を行う場合には、まず、補正する対象となるクロマトグラムが表示されているクロマトグラム領域１２２（図４参照）を選択する。
【００６７】
図５は、表示画面１２１の表示態様の一例であって、クロマトグラムＡに対して手動波形処理が行われる際の表示例を示した概略図である。
ここでは、ユーザが、第１クロマトグラム領域１２３を選択した場合を例に挙げて説明する。
【００６８】
ユーザの操作部１１の操作によって、第１クロマトグラム領域１２３が指定されると、表示処理部１４２は、表示画面１２１において、第１クロマトグラム領域１２３のクロマトグラムＡを拡大して表示する。この状態において、ユーザは、例えば、ドラッグ操作などにより、終了点Ｆを選択した後、グラフＡ上の点Ｒを選択（指定）する。
【００６９】
すると、手動波形処理部１４４は、ユーザの指定を受け付け、クロマトグラムＡにおいて、ピークＰ１の終了点を点Ｆから点Ｒに変更する処理を行う（図３のステップＳ１０３：手動波形処理ステップ）。また、表示処理部１４２は、ピークの開始点Ｅと、変更した後のピークの終了点Ｒとを結ぶベースラインＳを表示する。このとき、変更されたピークＰ１の情報（新たなピークの終了点Ｒの時間の情報）は、新たなピーク情報１３４として、記憶部１３において上書きされる。そして、記憶部１３において、この新たなピーク情報１３４に対応づけられているクロマトグラムデータ１３２は、解析済みデータとして記憶される（図３のステップＳ１０４）。また、記憶部１３では、クロマトグラムデータ１３２のうち、解析済みデータ以外のデータが、対象データとされる。
【００７０】
そして、ユーザによる手動操作が終了すると、表示処理部１４２は、図６に示すように、再度、各クロマトグラム領域１２２を表示画面１２１に表示させる。このとき、表示処理部１４２は、変更した後のピークＰ１を含むクロマトグラムＡを、第１クロマトグラム領域１２３に表示させる（表示処理ステップ）。
図６は、表示画面１２１の表示態様の一例であって、解析済みデータ及び対象データを表示する際の表示例を示した概略図である。
【００７１】
この例では、ユーザは、上記の第１クロマトグラム領域１２３に対する手動操作と同様にして、第３クロマトグラム領域１２５を指定し、クロマトグラムＣにおいて、ピークＰ３を補正する。これにより、第３クロマトグラム領域１２５では、クロマトグラムＣにおいて、ピークの開始点Ｋと、新たなピークの終了点Ｔとを結ぶベースラインＵが表示されている。また、ピークに関する情報は、新たなピーク情報１３４として、記憶部１３において上書きされる。そして、記憶部１３において、この新たなピーク情報１３４に対応づけられているクロマトグラムデータ１３２は、解析済みデータとして記憶される。また、この状態において、図４に示す場合と同様に、各クロマトグラム領域１２２では、分析対象であるクロマトグラムのピークのピーク位置が矢印で示される。さらに、各クロマトグラム領域１２２では、このピークと同一成分のピークであって、分析日時が直前であるクロマトグラムのピーク位置、及び、分析日時が直後であるクロマトグラムのピーク位置のそれぞれが、分析対象のピークに対応付けるように矢印で示されている。
【００７２】
このように、表示画面１２１のクロマトグラム領域１２２には、ユーザの手動操作によってピークに関する情報が変更されたクロマトグラム（解析済みデータに対応するクロマトグラム）、及び、自動波形処理のみが行われたクロマトグラム（対象データに対応するクロマトグラム）が表示される。具体的には、図６では、第１クロマトグラム領域１２３及び第３クロマトグラム領域１２５に、解析済みデータに対応するクロマトグラムＡ及びＣが表示され、第２クロマトグラム領域１２４及び第４クロマトグラム領域１２６に、対象データに対応するクロマトグラムＢ及びＤが表示される。また、これらは、分析日時の順序に従って、上下方向に沿って並ぶように配置されている。
【００７３】
この状態から、ユーザは、操作部１１を操作して、解析済みデータにおけるピーク情報１３４の内容を、対象データのピーク情報１３４に反映させる。具体的には、ユーザは、操作部１１を操作して、解析済みデータにおける時間範囲を設定する。この例では、ユーザは、第１クロマトグラム領域１２３において、時間ｔ１及び時間ｔ２を設定する。そして、時間範囲指定部１４５は、設定された時間範囲、すなわち、時間ｔ１から時間ｔ２までを指定する（図３のステップＳ１０５でＹＥＳ：時間範囲指定ステップ）。
【００７４】
また、ピーク選択部１４６は、時間範囲指定部１４５が指定した時間範囲、及び、各クロマトグラム領域１２２に表示されるクロマトグラムの分析日時の順序（クロマトグラムの生成日時の順序）に基づいて、対象データにおけるピークを選択する（図３のステップＳ１０６：ピーク選択ステップ）。
【００７５】
具体的には、ピーク選択部１４６は、対象データのクロマトグラムのうち、分析日時がクロマトグラムＡの直後であるクロマトグラムＢを選択し、そのクロマトグラムＢにおいて、クロマトグラムＡのピークＰ１（手動波形処理部１４４によって内容が変更されたピークであって、指定された時間範囲に含まれるピーク）に最も近い時間に位置するピークＰ２を、ピークＰ１と同一成分のピークとして選択する。また、ピーク選択部１４６は、対象データのクロマトグラムのうち、分析日時がクロマトグラムＣの直後であるクロマトグラムＤを選択し、そのクロマトグラムＤにおいて、クロマトグラムＡのピークＰ１に最も近い時間に位置するピークＰ４を、ピークＰ１（手動波形処理部１４４によって内容が変更されたピークであって、指定された時間範囲に含まれるピーク）と同一成分のピークとして選択する。
そして、時間差算出部１４７は、２つの解析済みデータにおける同一成分のピーク間の溶出時間差を算出する（図３のステップＳ１０７：時間差算出ステップ）。
【００７６】
具体的には、時間差算出部１４７は、第１クロマトグラム領域１２３に表示されるクロマトグラムＡのピーク位置の時間（矢印Ａ１で示される時間）ｔ３と、第３クロマトグラム領域１２５に表示されるクロマトグラムＣのピーク位置の時間（矢印Ｃ３で示される時間）ｔ４との間の溶出時間差Δｔを算出する。
【００７７】
また、ピーク情報反映部１４８は、時間差算出部１４７が算出した溶出時間差Δｔ、及び、各クロマトグラム領域１２２に表示されるクロマトグラムの分析日時の順序（クロマトグラムの生成日時の順序）に基づいて、解析済みデータにおけるピーク情報の内容を、ピーク選択部１４６が選択した対象データのピーク情報に反映させる（図３のステップＳ１０８：ピーク情報反映ステップ）。
図７は、表示画面１２１の表示態様の一例であって、解析済みデータのピーク情報が反映された対象データを表示する際の表示例を示した概略図である。
【００７８】
クロマトグラム領域１２２では、クロマトグラムの分析日時が後になるたびに、同一成分のピークが同じ時間だけ遅れることから、ピーク情報反映部１４８は、以下のように、ピーク情報を反映させる。
【００７９】
具体的には、ピーク情報反映部１４８は、第２クロマトグラム領域１２４において、クロマトグラムＡのピークの終了点Ｒの時間であるｔ５に、溶出時間差Δｔの半分を加えた時間をｔ６とし、クロマトグラムＢにおいて、この時間ｔ６に該当する点を新たな終了点Ｖとする。そして、表示処理部１４２は、第２クロマトグラム領域１２４において、ピークの開始点Ｈと、新たなピークの終了点Ｖとを結ぶベースラインＷを表示させる。
これは、一方の解析済みデータであるクロマトグラムＡの分析日時が１番目であって、かつ、他方の解析済みデータであるクロマトグラムＣの分析日時が３番目であること、及び、ピーク選択部１４６によって選択されたピークＰ２を含むクロマトグラムＢの分析日時が２番目であることに基づいている。したがって、例えば、一方の解析済みデータの分析日時が１番目であり、他方の解析済みデータの分析日時が２番目である場合には、対象データのピークの終了点に加えられる時間は、溶出時間差Δｔの半分ではなく、溶出時間差Δｔそのものとなる。
【００８０】
また、ピーク情報反映部１４８は、第４クロマトグラム領域１２６において、クロマトグラムＣのピークの終了点Ｔの時間であるｔ７に、溶出時間差Δｔの半分を加えた時間をｔ８とし、クロマトグラムＤにおいて、この時間ｔ８に該当する点を新たな終了点Ｘとする。そして、表示処理部１４２は、第４クロマトグラム領域１２６において、ピークの開始点Ｎと、新たなピークの終了点Ｘとを結ぶベースラインＹを表示させる。
これは、ピーク選択部１４６によって選択されたピークＰ４を含むクロマトグラムＤの分析日時が、解析済みデータであるクロマトグラムＣの分析日時の直後であることに基づいている。
【００８１】
このように、データ処理装置１では、２つの解析済みデータの溶出時間差Δｔを算出し、クロマトグラムの分析日時の順序に基づいて、解析済みのデータのピークに関する時間を、溶出時間差Δｔに基づく時間だけずらす処理を行うことにより、解析済みデータにおけるピーク情報１３４の変更内容を、対象データのピーク情報１３４に反映させる。
【００８２】
４．作用効果
（１）本実施形態では、自動波形処理部１４３によって、各クロマトグラムデータ１３２で表されるクロマトグラムのピークが自動で検出された後（自動波形処理ステップ）、そのピークに対してユーザが手動操作を行うと、手動波形処理部１４４は、そのピークの情報（ピーク情報）を変更する（手動波形処理ステップ）。また、変更された情報は、新たなピーク情報１３４として、記憶部１３において上書きされる。そして、ピーク情報反映部１４８は、その変更されたピーク情報１３４の内容を、他のクロマトグラムデータ（対象データ）におけるピークのピーク情報１３４に反映させる（ピーク情報反映ステップ）。
【００８３】
そのため、対象データにおけるピークのピーク情報１３４を、手動操作によらずに自動的に変更できる。
その結果、対象データを補正する作業を効率化できる。
【００８４】
また、時間差算出部１４７は、２つの解析済みデータにおける同一成分のピーク間の溶出時間差Δｔを算出する（時間差算出ステップ）。ピーク選択部１４６は、対象データについて、解析済みデータのピーク（クロマトグラムＡのピークＰ１）と同一成分のピークを、分析日時（分析日時の順序）、すなわち、クロマトグラムデータの生成日時の順序に基づいて選択する（ピーク選択ステップ）。ピーク情報反映部１４８は、溶出時間差Δｔ及び分析日時（分析日時の順序）に基づいて、解析済みデータにおけるピークのピーク情報１３４の変更を、対象データにおいて選択されたピークのピーク情報１３４に反映させる（ピーク情報反映ステップ）。
【００８５】
すなわち、データ処理装置１では、まず、ピーク選択部１４６によって、対象データについて、解析済みデータのピーク（クロマトグラムＡのピークＰ１）と同一成分のピークを正確に選択する。そして、データ処理装置１では、２つの解析済みデータにおける同一成分のピーク間の溶出時間差Δｔ、及び、分析日時（分析日時の順序）に基づくことにより、クロマトグラムごとでのピークの溶出時間のずれを考慮して、解析済みデータにおけるピーク情報１３４の変更内容（新たなピークの時間に関する情報）を、その選択されたピークのピーク情報１３４に反映させる。
【００８６】
そのため、対象データにおけるピークのピーク情報１３４の変更（対象データの補正）を精度よく行うことができる。
このように、本実施形態のデータ処理装置１及びデータ処理方法によれば、自動で波形処理されたデータ（ピーク）を補正する作業を効率化でき、かつ、当該補正を精度よく行うことができる。
【００８７】
（２）また、本実施形態では、表示処理部１４２は、ユーザの手動操作によってピークに関する情報が変更されたクロマトグラム（解析済みデータに対応するクロマトグラム）、及び、自動波形処理のみが行われたクロマトグラム（対象データに対応するクロマトグラム）を、表示画面１２１において、分析日時の順序に従って、上下方向に沿って並ぶように表示させる（表示処理ステップ）。
【００８８】
そのため、ユーザは、各クロマトグラムにおける同一成分のピークのピーク位置を、分析日時の順序と照らし合わせて確認できる。
その結果、分析日時の順序と、ピーク位置のずれとの関係を容易に確認できる。
【００８９】
（３）また、本実施形態では、表示処理部１４２は、各クロマトグラム領域１２２において、分析対象であるクロマトグラムのピークのピーク位置を矢印で示す。さらに、表示処理部１４２は、各クロマトグラム領域１２２において、このピークと同一成分のピークであって、分析日時が直前であるクロマトグラムのピーク位置、及び、分析日時が直後であるクロマトグラムのピーク位置のそれぞれを、分析対象のピークに対応付けるように矢印で示す。
【００９０】
そのため、ユーザは、同一成分のピークに関して、あるデータにおけるピーク位置を、分析日時が直前及び直後のデータにおけるピーク位置と比較して確認できる。
【００９１】
（４）また、本実施形態では、時間範囲指定部１４５は、解析済みデータにおける時間範囲を指定する。ピーク選択部１４６は、時間範囲指定部１４５により指定された時間範囲及び分析日時に基づいて、ピーク情報１３４が変更されたピークと同一成分のピークを選択する。
【００９２】
そのため、時間範囲指定部１４５によって時間範囲を指定することにより、特定の成分に限定してピーク情報１３４を変更できる。
そのため、全体の処理時間を短くできる。
【００９３】
５．変形例
以上の説明では、データ処理装置１は、ガスクロマトグラフ２とともに用いられるとして説明した。しかし、データ処理装置１は、ガスクロマトグラフ（ＧＣ）以外の装置とともに用いることが可能である。例えば、データ処理装置１は、液体クロマトグラフ（ＬＣ）などとともに用いることが可能である。
【００９４】
また、以上の説明では、ユーザは、手動操作によって、クロマトグラムにおけるピークの終了点を変更するとして説明した。しかし、ユーザは、クロマトグラムにおけるピークに関して、その他の情報（ピーク情報）を変更することも可能である。そして、ピーク情報反映部１４８は、ユーザによって変更されたその他のピーク情報の内容を、他のクロマトグラムのピーク情報に反映させることも可能である。例えば、ユーザは、手動操作によって、クロマトグラムにおけるピークの開始点、及び、ピーク位置を変更することが可能である。そして、ピーク情報反映部１４８は、変更されたピークの開始点、及び、ピーク位置の情報を、他のクロマトグラムのピーク情報に反映させることも可能である。
【００９５】
また、以上の説明では、対象データは、２つであるとして説明した。しかし、対象データは、３つ以上であってもよい。
【符号の説明】
【００９６】
１      データ処理装置
１４      制御部
１３１      検出データ
１３２      クロマトグラムデータ
１３４      ピーク情報
１４１      データ処理部
１４２      表示処理部
１４３      自動波形処理部
１４４      手動波形処理部
１４５      時間範囲指定部
１４６      ピーク選択部
１４７      時間差算出部
１４８      ピーク情報反映部
【図１】
2017179154.tif 000003
【図２】
2017179154.tif 000004
【図３】
2017179154.tif 000005
【図４】
2017179154.tif 000006
【図５】
2017179154.tif 000007
【図６】
2017179154.tif 000008
【図７】
2017179154.tif 000009
【国際調査報告】
2017179154.tif 000010
2017179154.tif 000011