(57)【要約】  取得された蛍光指紋情報に対して新規な２次微分処理を施すことにより、精度の高い推定モデル（検量線）を得る。  特定成分の含有量が既知の試験試料を準備し（Ｓ０１），蛍光指紋情報を取得する（Ｓ０２）。この蛍光指紋情報に対し、新規な２次微分処理を含む前処理を行い（Ｓ０３）、蛍光指紋情報から試験試料中の特定成分の含有量を推定する推定モデル（検量線）を作成する（Ｓ０４）。検量線の検証（Ｓ０５）後、未知の試料に適用して未知試料に含有される特定成分の含有量を推定する（Ｓ０６）。
2020031447.tif 000002
【特許請求の範囲】
【請求項１】
試験試料について励起波長・蛍光波長・蛍光強度のデータからなる蛍光指紋情報を取得する蛍光指紋情報取得工程と、
２次微分処理により前記試験試料に特有な前記蛍光強度のピーク値が現れるように軸を設定し、前記軸に沿った前記蛍光指紋情報の２次微分処理を少なくとも含む前処理工程と、
少なくとも前記２次微分処理が実施された蛍光指紋情報を説明変数とし、前記試験試料についての既知の定量値を目標変数として、検量線を取得する推定モデル作成工程と、
を含む、蛍光指紋分析による試験試料の評価方法。
【請求項２】
前記推定モデル作成工程において、多変量解析によって上記検量線を作成することを特徴とする請求項１に記載の蛍光指紋分析による試験試料の評価方法。
【請求項３】
前記多変量解析は、PLS回帰分析であることを特徴とする請求項２に記載の蛍光指紋分析による試験試料の評価方法。
【請求項４】
前記前処理工程において、前記蛍光指紋情報に対して低感度領域の削除処理を行うことを特徴とする請求項１〜３の何れか１項に記載の蛍光指紋分析による試験試料の評価方法。
【請求項５】
前記試験試料が、クロロゲン酸を含む、請求項１〜４の何れか１項に記載の蛍光指紋分析による試験試料の評価方法。
【請求項６】
前記試験試料が、スコポレチンを含む、請求項１〜４の何れか１項に記載の蛍光指紋分析による試験試料の評価方法。
【請求項７】
前記試験試料が、ルチンを含む、請求項１〜４の何れか１項に記載の蛍光指紋分析による試験試料の評価方法。
【請求項８】
前記励起波長が２８５nm近傍値であって前記蛍光波長が４１０，４１５，４２０，４２５，４３０, ４３５，４４０, ４５０，４６０nm近傍値のうちの少なくとも１つである波長のいずれかの励起波長／蛍光波長の組み合わせにより、ルチンの検量線を作成することを特徴とする請求項７に記載の蛍光指紋分析による試験試料の評価方法。
【請求項９】
前記励起波長が２８５nm近傍値であって前記蛍光波長が４１０，４１５，４２０，４２５，４３０, ４３５，４４０, ４５０，４６０nm近傍値の全てである波長の励起波長／蛍光波長の組み合わせの両方により、ルチンの検量線を作成することを特徴とする請求項８に記載の蛍光指紋分析による試験試料の評価方法。
【請求項１０】
前記試験試料は、たばこ製品の原料であることを特徴とする請求項１〜９の何れか１項に記載の蛍光指紋分析による試験試料の評価方法。
【請求項１１】
前記試験試料は、励起光の照射前に粉末状に粉砕・混合されることを特徴とする請求項１〜１０の何れか１項に記載の蛍光指紋分析による試験試料の評価方法。
【請求項１２】
前記試験粉砕によって、試料が１ｍｍ以下の粒径とされることを特徴とする請求項９に記載の蛍光指紋分析による試験試料の評価方法。
【請求項１３】
前記試験材料は、事前に水分量を一定化するために、所定の調和条件で所定時間蔵置されることを特徴とする請求項１〜１２の何れか１項に記載の蛍光指紋分析による試験試料の評価方法。
【請求項１４】
前記試験材料は、たばこ製品の原料であり、前記調和条件は、温度２２℃、湿度６０％の室内という条件であり、前記所定時間は２４時間以上であることを特徴とする請求項１３に記載の蛍光指紋分析による試験試料の評価方法。
【請求項１５】
請求項１〜１４の何れか１項に記載の蛍光指紋分析による試験試料の評価方法により得られた検量線と未知の試料の蛍光指紋情報とに基づき、前記試料に含有される特定成分の含有量を推定することを特徴とする成分量推定方法。
【請求項１６】
前記試料に含有される特定成分の含有量を推定する際に、前記蛍光指紋分析による試験試料の評価方法における前処理と同一の処理を行うことを特徴とする請求項１５に記載の成分量推定方法。
【請求項１７】
前記未知の試料の蛍光指紋情報を取得するために、励起波長が２８５nm近傍値であって蛍光波長が４１０，４１５，４２０，４２５，４３０, ４３５，４４０, ４５０，４６０nm近傍値のうち少なくとも１つである波長のいずれかの励起波長／蛍光波長の組み合わせを用いることを特徴とする請求項１５又は１６に記載の成分量推定方法。
【請求項１８】
前記未知の資料の蛍光指紋情報を取得するために、前記励起波長が２８５nm近傍値であって前記蛍光波長が４１０，４１５，４２０，４２５，４３０, ４３５，４４０，４５０，４６０nm近傍値の全てである波長の励起波長／蛍光波長の組み合わせの両方を用いることを特徴とする請求項１７に記載の成分量推定方法。
【請求項１９】
コンピュータに請求項１〜１８の何れか１項に記載の方法を実行させるためのプログラム。
【請求項２０】
試料についての励起波長・蛍光波長・蛍光強度のデータからなる蛍光指紋情報を入力し、２次微分処理により前記試料に特有な前記蛍光強度のピーク値が現れるように軸を設定し、前記軸に沿った前記蛍光指紋情報の２次微分処理を少なくとも含む前処理手段と、
少なくとも前記２次微分処理が実施された蛍光指紋情報を説明変数とし、前記試験試料についての既知の定量値を目標変数として、検量線を取得する推定モデル作成手段と、
前記推定モデル作成手段により取得された前記検量線と未知の試料の蛍光指紋情報とに基づき、前記未知の試料に含有される特定成分の含有量を推定する成分量推定手段と、
を具備することを特徴とする装置。
【請求項２１】
前記未知の試料の蛍光指紋情報は前記前処理手段で処理され、処理後の蛍光指紋情報が前記成分量推定手段に入力されることを特徴とする請求項２０に記載の装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、蛍光指紋分析を利用した試料の評価・推定方法、プログラム、及び装置に関し、特に、試料中の微量成分の評価・推定に好適に利用し得る評価・推定方法、プログラム、及び装置に関する。
【背景技術】
【０００２】
図２に示すように、蛍光物質を含む試験試料に、段階的に波長を変化させながら励起光を照射し、試験試料から発せられる光（蛍光）を測定すると、励起波長（λＥｘ）、蛍光波長（測定波長）（λＥｍ）、蛍光強度（ＩＥｘ，Ｅｍ）を３直交軸とする3次元空間において、対応するポイントがプロットされる。
これらのポイントの集合を可視化したものを、蛍光指紋、または、励起蛍光マトリクス（Excitation Emission Matrix；ＥＥＭ）と呼ぶ。
蛍光指紋は、各ポイントの蛍光強度を等高線形状や色分布等で表示することにより、３次元グラフとして表すことができ（図３参照）、また、２次元グラフとして表すこともできる（図４参照）。【０００３】
このような蛍光指紋は、３次元の膨大な情報を有する試験試料固有のパターンを示しており、各種の鑑別や定量等に使用できる。例えば、非特許文献１には、蛍光指紋を利用した蕎麦粉混合割合の推定手法等が紹介されている。【０００４】
このような蛍光指紋を用いた分析は、他の分光分析手法と比較して、試験試料に対して蛍光染色などの前処理を行うことなく試験試料のキャラクタリゼーションが可能であること、操作が容易で計測も短時間で行えること、情報量が多く定量も比較的容易に行えること、非破壊・非接触での計測が可能であること、紫外〜可視光を用いていることから装置が比較的安価であること、などの長所を有している。【０００５】
蛍光指紋を計測するには、励起波長をスキャンして、それぞれの蛍光スペクトルを連続的に計測する機能を有する蛍光分光光度計が必要であるが、このような機能を有する蛍光分光光度計も市販されている（日立ハイテクサイエンス社製Ｆ−7000等）。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２０１７−３６９９１号公報
【特許文献２】特開２０１７−５１１６２号公報
【非特許文献】
【０００７】
【非特許文献１】杉山純一「光の指紋による食品の鑑別・定量」；「食品と容器」2013 VOL.54 NO.5 P308-315
【発明の概要】
【発明が解決しようとする課題】
【０００８】
蛍光指紋分析を利用した既知の試料の評価・推定手法においては、前処理を行っても、精度の高い推定モデル（検量線）の作成が困難な場合があった。特に、化学構造が類似した成分を含む試料中から特定の成分の定量を行う場合に困難性が高かった。
【課題を解決するための手段】
【０００９】
本発明は、このような課題を解決するために提案されたものであり、実施の態様を例示すれば、以下のとおりである。
（態様１）
試験試料について励起波長・蛍光波長・蛍光強度のデータからなる蛍光指紋情報を取得する蛍光指紋情報取得工程と、２次微分処理により前記試験試料に特有な前記蛍光強度のピーク値が現れるように軸を設定し、前記軸に沿った前記蛍光指紋情報の２次微分処理を少なくとも含む前処理工程と、少なくとも前記２次微分処理が実施された蛍光指紋情報を説明変数とし、前記試験試料についての既知の定量値を目標変数として、検量線を取得する推定モデル作成工程と、を含む、蛍光指紋分析による試験試料の評価方法。
（態様２）
前記推定モデル作成工程において、多変量解析によって上記検量線を作成することを特徴とする態様１に記載の蛍光指紋分析による試験試料の評価方法。
（態様３）
前記多変量解析は、PLS回帰分析であることを特徴とする態様２に記載の蛍光指紋分析による試験試料の評価方法。
（態様４）
前記前処理工程において、前記蛍光指紋情報に対して低感度領域の削除処理を行うことを特徴とする態様１〜３に記載の蛍光指紋分析による試験試料の評価方法。
（態様５）
前記試験試料が、クロロゲン酸を含む、態様１〜４に記載の蛍光指紋分析による試験試料の評価方法。
（態様６）
前記試験試料が、スコポレチンを含む、態様１〜４に記載の蛍光指紋分析による試験試料の評価方法。
（態様７）
前記試験試料が、ルチンを含む、態様１〜４に記載の蛍光指紋分析による試験試料の評価方法。
（態様８）
前記励起波長が２８５nm近傍値であって前記蛍光波長が４１０，４１５，４２０，４２５，４３０, ４３５，４４０, ４５０，４６０nm近傍値のうちの少なくとも１つである波長のいずれかの励起波長／蛍光波長の組み合わせにより、ルチンの検量線を作成することを特徴とする態様７に記載の蛍光指紋分析による試験試料の評価方法。
（態様９）
前記励起波長が２８５nm近傍値であって前記蛍光波長が４１０，４１５，４２０，４２５，４３０, ４３５，４４０, ４５０，４６０nm近傍値の全てである波長の励起波長／蛍光波長の組み合わせの両方により、ルチンの検量線を作成することを特徴とする態様８に記載の蛍光指紋分析による試験試料の評価方法。
（態様１０）
前記試験試料は、たばこ製品の原料であることを特徴とする態様１〜９に記載の蛍光指紋分析による試験試料の評価方法。
（態様１１）
前記試験試料は、励起光の照射前に粉末状に粉砕・混合されることを特徴とする態様１〜１０に記載の蛍光指紋分析による試験試料の評価方法。
（態様１２）
前記試験粉砕によって、試料が１ｍｍ以下の粒径とされることを特徴とする態様９に記載の蛍光指紋分析による試験試料の評価方法。
（態様１３）
前記試験材料は、事前に水分量を一定化するために、所定の調和条件で所定時間蔵置されることを特徴とする態様１〜１２に記載の蛍光指紋分析による試験試料の評価方法。
（態様１４）
前記試験材料は、たばこ製品の原料であり、前記調和条件は、温度２２℃、湿度６０％の室内という条件であり、前記所定時間は２４時間以上であることを特徴とする態様１３に記載の蛍光指紋分析による試験試料の評価方法。
（態様１５）
態様１〜１４に記載の蛍光指紋分析による試験試料の評価方法により得られた検量線と未知の試料の蛍光指紋情報とに基づき、前記試料に含有される特定成分の含有量を推定することを特徴とする成分量推定方法。
（態様１６）
前記試料に含有される特定成分の含有量を推定する際に、前記蛍光指紋分析による試験試料の評価方法における前処理と同一の処理を行うことを特徴とする態様１５に記載の成分量推定方法。
（態様１７）
前記未知の試料の蛍光指紋情報を取得するために、励起波長が２８５nm近傍値であって蛍光波長が４１０，４１５，４２０，４２５，４３０, ４３５，４４０, ４５０，４６０nm近傍値のうち少なくとも１つである波長のいずれかの励起波長／蛍光波長の組み合わせを用いることを特徴とする態様１５、１６に記載の成分量推定方法。
（態様１８）
前記未知の資料の蛍光指紋情報を取得するために、前記励起波長が２８５nm近傍値であって前記蛍光波長が４１０，４１５，４２０，４２５，４３０, ４３５，４４０，４５０，４６０nm近傍値の全てである波長の励起波長／蛍光波長の組み合わせの両方を用いることを特徴とする態様１７に記載の成分量推定方法。
（態様１９）
コンピュータに態様１〜１８に記載の方法を実行させるためのプログラム。
（態様２０）
試料についての励起波長・蛍光波長・蛍光強度のデータからなる蛍光指紋情報を入力し、２次微分処理により前記試料に特有な前記蛍光強度のピーク値が現れるように軸を設定し、前記軸に沿った前記蛍光指紋情報の２次微分処理を少なくとも含む前処理手段と、少なくとも前記２次微分処理が実施された蛍光指紋情報を説明変数とし、前記試験試料についての既知の定量値を目標変数として、検量線を取得する推定モデル作成手段と、前記推定モデル作成手段により取得された前記検量線と未知の試料の蛍光指紋情報とに基づき、前記未知の試料に含有される特定成分の含有量を推定する成分量推定手段と、を具備することを特徴とする装置。
（態様２１）
前記未知の試料の蛍光指紋情報は前記前処理手段で処理され、処理後の蛍光指紋情報が前記成分量推定手段に入力されることを特徴とする態様２０に記載の装置。【００１０】
なお、上述の態様において、「近傍値」という語句は、励起波長や蛍光（測定）波長が誤差を必然的に伴うことから、励起波長や蛍光（測定）波長の誤差範囲に存在する波長をも包含することを明確化するために付したものである。因みに、誤差範囲は測定機器や測定条件等により変動し得るものである。【００１１】
また、「プログラム」とは、任意の言語や記述方法に基づき記述されたデータ処理方法であり、ソースコードやバイナリコード等の形式を問うものではない。また、「プログラム」は単一の形で構成されてもよいが、複数のモジュールやライブラリとして分散構成されてもよく、また、他の既存のプログラムと協働してその機能を達成するように構成されたものであってもよい。【００１２】
また、「装置」は、ハードウエアとして構成されてもよいが、コンピュータのソフトウエアによって各種機能を実現する機能実現手段の組合せとして構成されてもよい。機能実現手段には、例えば、プログラムモジュールが含まれ得る。
【発明の効果】
【００１３】
本発明によれば、精度の高い推定モデル（検量線）の作成が可能になる。また、特に、化学構造が類似した成分を含む試料中から特定の成分の定量を的確に行うことができる。
【図面の簡単な説明】
【００１４】
【図１】図１は、本発明の実施の一態様の概要を説明するためのフローチャートである。
【図２】図２は、計測対象物に励起光を照射した場合に該計測対象物から発せられた蛍光のスペクトルの概要を示す説明図である。
【図３】図３は、蛍光指紋の一例を３次元的に示す等高線形状のグラフである。
【図４】図４は、蛍光指紋の一例を２次元的に示す等高線形状のグラフである。
【図５】図５は、スコポレチン、クロロゲン酸、ルチンを含む混合サンプルの蛍光指紋の一例を表す図である。
【図６Ａ】図６Ａは、を単体で含むサンプルの蛍光指紋の一例を表す図である。
【図６Ｂ】図６Ｂは、クロロゲン酸を単体で含むサンプルの蛍光指紋の一例を表す図である。
【図６Ｃ】図６Ｃは、ルチンを単体で含むサンプルの蛍光指紋の一例を表す図である。
【図７Ａ】図７Ａは、複数のサンプルについて、スコポレチンの既知の含有量と蛍光指紋情報による推定値によって規定される点をプロットしたグラフである。
【図７Ｂ】図７Ｂは、複数のサンプルについて、クロロゲン酸の既知の含有量と蛍光指紋情報による推定値によって規定される点をプロットしたグラフである。
【図７Ｃ】図７Ｃは、複数のサンプルについて、ルチンの既知の含有量と蛍光指紋情報による推定値によって規定される点をプロットしたグラフである。
【図８Ａ】図８Ａは、複数のサンプルについて、ルチン含有量の実測値（化学分析値）と前処理を行った蛍光指紋情報による推定値によって規定される点をプロットしたグラフである。
【図８Ｂ】図８Ｂは、複数のサンプルについて、クロロゲン酸含有量の実測値（化学分析値）と前処理を行った蛍光指紋情報による推定値によって規定される点をプロットしたグラフである。
【図９】図９は、複数のサンプルについて、ルチン含有量の実測値（化学分析値）と前処理を省いた蛍光指紋情報による推定値によって規定される点をプロットしたグラフである。
【図１０】図１０は、本発明の実施の別の一態様の概要を説明するためのブロック図である。
【発明を実施するための形態】
【００１５】
以下、本発明の実施の一態様を説明するとともに、試料中に微量成分として含まれるクロロゲン酸、スコポレチン、ルチンの含有量の本発明に基づく評価・推定方法について説明する。また、特に、タバコ原料を試料とした場合について、たばこ原料中のルチン含有量の本発明に基づく評価・推定方法について説明する。
なお、この実施の一態様により、本発明が限定されるものではないことに留意されたい。
【００１６】
＜本発明の実施の一態様の概要＞
図１は、本発明の実施の一態様の概要を説明するためのフローチャートである。
【００１７】
最初に、成分の種類や含有量が既知の試験試料を準備し（Ｓ０１）、このような既知の試験試料に対して蛍光指紋を測定し、蛍光指紋情報を取得する（Ｓ０２）。
【００１８】
次に、取得された蛍光指紋情報に対する前処理を行う（Ｓ０３）。本発明の実施の一態様は、特に、この前処理に特徴を有しており、前処理の詳細は後述する。
【００１９】
次に、前処理済みの蛍光指紋情報と特定成分の既知の含有量との関係をモデル化し、推定モデル（検量線）を作成する（Ｓ０４）。このモデル化は、具体的には、前処理済みの蛍光指紋情報を説明変数、既知の含有量を目的変数とする推定式を、様々な多変量解析手法やデータマイニング手法を使用して構築することにより行われ、蛍光指紋情報から試験試料中の特定成分の含有量を推定する検量線（回帰式）が作成される。なお、推定式の構築に用いるアルゴリズムはより汎用的かつ非線形現象にも対応する機械学習アルゴリズム、例えばsupport vector machine (SVM), random forest (RF), neural networkなどでも良い。モデル化の際に用いられる多変量解析手法の例は後述する。
【００２０】
このようにして構築された推定モデル（検量線）を検証し、その有効性を確認する（Ｓ０５）。
【００２１】
有効性が確認された推定モデル（検量線）を用いて、未知の試料の蛍光指紋情報に基づき、前記未知の試料に含有される特定成分の含有量を推定する（Ｓ０６）。
【００２２】
＜本発明の実施の別の一態様の概要＞
図１０は、本発明の実施の別の一態様の概要を説明するためのブロック図である。
【００２３】
蛍光指紋分析による試料の評価・推定装置１００は、試料についての励起波長・蛍光波長・蛍光強度のデータからなる蛍光指紋情報を入力とし、２次微分処理により前記試料に特有な前記蛍光強度のピーク値が現れるように軸を設定し、前記軸に沿った前記蛍光指紋情報の２次微分処理を少なくとも含む前処理手段１１０と、前記前処理手段１１０の出力を入力とし、少なくとも前記２次微分処理が実施された蛍光指紋情報を説明変数とし、前記試験試料についての既知の定量値を目標変数として、検量線を取得する推定モデル作成手段１２０と、前記推定モデル作成手段１２０により取得された前記検量線と未知の試料の蛍光指紋情報とに基づき、前記未知の試料に含有される特定成分の含有量を推定する成分量推定手段１３０とを具備している。
【００２４】
まず、既知の蛍光分光光度計等を用いて、成分の種類や含有量が既知の試験試料の蛍光指紋情報を取得する。
【００２５】
次に、取得された蛍光指紋情報を前記前処理手段１１０に入力し、入力された蛍光指紋情報に対する前処理を行う。本態様も、前述の本発明の実施の一態様と同様、特に、この前処理に特徴を有しており、前処理の詳細は後述する。
【００２６】
次に、推定モデル作成手段１２０により、前処理済みの蛍光指紋情報と特定成分の既知の含有量との関係をモデル化し、推定モデル（検量線）を作成する。このモデル化は、前述の本発明の実施の一態様と同様である。そして、このようにして構築された推定モデル（検量線）を検証し、その有効性を確認し、有効性が確認された推定モデル（検量線）を不図示のメモリ等に記憶しておく。
【００２７】
成分量推定手段１３０は、有効性が確認された推定モデル（検量線）を用いて、未知の試料の蛍光指紋情報に基づき、前記未知の試料に含有される特定成分の含有量を推定する。なお、未知の試料の蛍光指紋情報に対しても、前記前処理手段１１０による前処理を行うことが望ましい（図１０の蛍光指紋分析による試料の評価・推定装置１００は、そのような構成を採用している）が、必要に応じてこのような前処理を省略することもできる。
【００２８】
以下、蛍光指紋情報に対する前処理、モデル化の際に用いられる多変量解析手法、試料中に微量成分として含まれるクロロゲン酸、スコポレチン、ルチンの含有量の本発明に基づく評価・推定方法について説明する。なた、特に、タバコ原料を試料とした場合について、たばこ原料中のルチン含有量の本発明に基づく評価・推定方法について説明する。
【００２９】
＜蛍光指紋情報に対する前処理＞
試験試料の蛍光指紋情報を取得する際には、蛍光指紋の計測値（励起波長毎の蛍光スペクトル）をそのまま用いることもできるが、必要に応じて、各種の前処理を行う必要がある。
【００３０】
本発明の実施の一態様は、特に、前処理において、以下に詳述するような新規な手法（本発明に特有な２次微分処理）を採用することにより、測定精度の向上を図るものである。
【００３１】
３次元空間における、励起波長（λＥｘ）、蛍光波長（測定波長）（λＥｍ）、蛍光強度（ＩＥｘ，Ｅｍ）の３直交軸を、それぞれ、ｙ軸、ｘ軸、ｚ軸とし、試験試料の測定で得られたＥＥＭ（蛍光指紋）の各ポイントＰi（１≦ｉ≦Ｎ）の座標を（ｘi，ｙi，ｚi）で表す。また、ｘ軸及びｙ軸で決定される平面をｘ−ｙ平面、ｘ軸及びｚ軸で決定される平面をｘ−ｚ平面、ｙ軸及びｚ軸で決定される平面をｙ−ｚ平面と呼ぶ。
【００３２】
ｘ−ｙ平面上にｗ軸を取り、ｗ軸と交わりｚ軸と平行な軸をｚ’軸とすると、ｘ−ｙ平面と直交するｗ−ｚ’平面が得られる、このｗ−ｚ’平面を平行移動させ、前記蛍光指紋を切って得られた各断面は、ｚ’＝ｆ（ｗ）という関数で表すことができる。この関数の２次微分値（ｄ2ｚ’/ｄｗ2＝（ｄ2/ｄｗ2）ｆ（ｗ））を求める。求めた２次微分値から、勾配の一番急激な変化点をピークとして求める。
なお、ピークは1つとは限らず、複数存在する可能性があり、この場合には、複数のピークに対応する情報を取得する。
【００３３】
そして、このような処理を行うことにより、１次微分では検出が困難な、埋もれ易い微かなピークであっても有効に抽出することができる。また１次微分処理と異なり、スペクトル取得時に問題となる波長依存のベースラインの補正を容易に行うことができる。
【００３４】
そこで、上述のようなｗ軸の設定をどのように行うべきかが問題となる。試験試料に最適なｗ軸（以下、「最適軸」という）を設定するために、以下に例示するような手法が有効である。
（１）ｗ軸の方向を固定し、ｗ−ｚ’平面を平行移動させると、当該ｗ軸について前記蛍光指紋の断面積（∫ｆ（ｗ）ｄｗ）の最大値が得られる。ｗ軸の方向を変えて、この断面積の最大値が最大となるｗ軸（ｗSmax）を求め、このｗSmaxを最適軸とする。
（２）ｗ軸の方向を固定し、ｗ−ｚ’平面を平行移動させると、当該ｗ軸について前記蛍光指紋の断面の底辺（ｘ−ｙ平面に接する辺）近傍の長さの最大値が得られる。ｗ軸の方向を変えて、この長さの最大値が最大となるｗ軸（ｗLmax）を求め、このｗLmaxを最適軸とする。
（３）蛍光指紋の各ポイントＰi（１≦ｉ≦Ｎ）からの距離の二乗の総和が最小となるような直線（最小二乗直線）を求めると、当該最小二乗直線を含むｗ−ｚ’平面が決定される。そこで、これにより決定されるｗ軸を最適軸とする。
（４）蛍光指紋の各ポイントＰi（１≦ｉ≦Ｎ）の中、前記蛍光指紋の底面近傍のポイント（ｘ−ｙ平面近傍のポイント）からの距離の二乗の総和が最小となるような直線（最小二乗直線）を求めると、当該最小二乗直線を含むｗ−ｚ’平面が決定される。そこで、これにより決定されるｗ軸を最適軸とする。
【００３５】
なお、最適軸の設定は上記手法に限られるものではなく、例えば、上記手法により得られたｗ軸を暫定的な最適軸として、この暫定的な最適軸をｘ−ｙ平面上で適宜回転して、前記試験試料に特有な前記蛍光強度のピーク値が現れるようにｗ軸を設定するようにしてもよい。
また、試験試料によっては、最適軸をｘ軸、又は、ｙ軸に平行な軸として設定する簡素化した手法を採用することもできる。
【００３６】
なお、特許文献２には、２次元展開された蛍光指紋に対する信号処理演算として２次微分が例示されているが（段落〔００２６〕等参照）、本発明の上記最適軸に係る技術思想については、記載も示唆もされていないことに留意されたい。
【００３７】
上記本発明に特有な２次微分処理の外に、計測された蛍光指紋からノイズを除去して有効な蛍光指紋情報を得るために、非蛍光成分の除去処理、散乱光の除去処理、低感度領域の削除処理の中の１つ又は複数の組合せを前処理として採用することができる。また、取得された蛍光指紋情報に対し、中心化、規格化、標準化、ベースライン補正、平滑化、オートスケーリング、対数変換（Log10）等の演算処理の中の１つ又は複数の組合せを前処理として採用することもできる。さらに、従前使用されている２次微分処理を併せて使用することもできる。
【００３８】
なお、上述の本発明に特有な２次微分処理とそれ以外の前処理を併用する場合に、その適用順序は適宜設定可能であるが、処理の効率化の観点から、非蛍光成分の除去処理、散乱光の除去処理、低感度領域の削除処理等の処理を先行させることが望ましい。
【００３９】
＜モデル化の際に用いられる多変量解析手法＞
モデル化の際に用いられる多変量解析手法として、ＰＬＳ（Partial Least Squares）回帰分析、重回帰分析、主成分回帰分析、最小二乗法等の各種の解析手法を用いることができる。
【００４０】
ＰＬＳ回帰分析は、主成分と目的変数との共分散が最大になるように主成分を抽出する手法であり、説明変数の間に強い相関がある場合（多重共線性を生ずる場合）に有効な手法である。
【００４１】
主成分回帰分析は、主成分の分散が最大になるように主成分を抽出する手法であり、説明変数のみを用いて主成分分析を行い、得られた主成分と目的変数との間で最小二乗法による重回帰分析を行うものである。
重回帰分析は、説明変数と目的変数との間で最小二乗法を適用するものであり、主成分回帰分析とは異なる特徴を有するものである。
【００４２】
なお、上述の各解析手法自体は周知であり、本発明のモデル化の際に特殊な処理を要請するものでもないので処理内容の詳細は省略するが、ＰＬＳに関しては、検量線の作成との関連で後述する。
【００４３】
＜試料中のクロロゲン酸、スコポルチン、ルチンの含有量の推定への適用例＞
次に、本発明の手法を試料中の微量成分（クロロゲン酸、スコポルチン、ルチン）の含有量の推定に適用した適用例について詳述する。
【００４４】
これらの微量成分（クロロゲン酸、スコポルチン、ルチン）はポリフェノールの一種であり、蕎麦や果皮に多く含まれる栄養素の１つとして知られており、抗酸化作用等を有することから、食品や化粧品等様々な業界から注目されている。
【００４５】
これらの微量成分の定量方法としては、試料中のこれらの微量成分を抽出液で抽出し、これを高速液体クロマトグラフ（HPLC）で定量する方法が知られているが、この方法では分析結果を得るまでに労力と時間を要するという問題がある。
【００４６】
これらの微量成分は、試料中に化学構造が類似したポリフェノールとして混在して存在する場合が多い（例えば、たばこ原料等）。このような場合、化学構造が類似したポリフェノールは類似した蛍光指紋を有するため、蛍光指紋法を適用しても、特定の微量成分だけを抜き出して正確に定量することには困難性が予想される。
【００４７】
上記のような化学構造が類似した微量成分の定量に関する困難性に係る事情等を、図５、図６Ａ〜図６Ｃに基づき説明する。
【００４８】
図５は、クロロゲン酸（１０ｐｐｍ）、ルチン（１ｐｐｍ）、スコポレチン（０．１ｐｐｍ）を含む混合サンプルの蛍光指紋を表している。他方、図６Ａ、図６Ｂ、図６Ｃは、混合サンプルに含まれる、スコポレチン（１ｐｐｍ：（Ａ））、クロロゲン酸（１ｐｐｍ：（Ｂ））、ルチン（１ｐｐｍ：（Ｃ））のそれぞれを単独で含む試験試料の蛍光指紋を表している。以下の表１は、前記混合サンプル、並びに、スコポレチン（１ｐｐｍ）、クロロゲン酸（１ｐｐｍ）、ルチン（１ｐｐｍ）の蛍光指紋のピーク位置近傍の波長（励起波長／蛍光波長）、及び、対応する蛍光強度を纏めたものである。
【００４９】
【表１】
2020031447.tif 000003
【００５０】
図５及び図６Ａ〜図６Ｃ並びに表１から、概略、以下のような知見が得られる。
・混合物の蛍光指紋において、ルチンのピーク位置近傍の波長はスコポレチンと一致している（図５の符号「Ａ」を付した箇所参照）。
・混合物の蛍光指紋において、クロロゲン酸にみられるピーク位置（図５の符号「Ｂ」を付した箇所参照）が現れている。
・混合物の蛍光指紋において、ルチンに関する特徴は、一見したところでは視認できない。
・試料中の含有量が１ｐｐｍの場合、蛍光強度は、スコポレチンとルチンで３桁相違する。
・但し、蛍光指紋のパターンは、スコポレチン、クロロゲン酸、ルチンでわずかに異なる。
【００５１】
そこで、スコポレチン、クロロゲン酸、ルチン間でわずかに異なる蛍光指紋のパターンに着目して、試料中のそれぞれの含有量を推定することを可能にすることが考えられる。本適用例は、取得された蛍光指紋の前処理において、前述のような新規な手法を採用することにより、簡便かつ迅速に、試料中のスコポレチン、クロロゲン酸、ルチンの含有量を推定することができるものである。
【００５２】
以下、本発明の実施の一態様の各工程について説明する。
〔試験試料の準備〕
＜ルチン・クロロゲン酸・スコポレチンを標品で混合した試料を用いる場合＞
ルチン・クロロゲン酸・スコポレチンの標品を様々な割合で混合した４０サンプル（溶媒は５０％エタノール／水溶液）を準備した。各サンプルについては、粒径１ｍｍ以下に粉砕し十分に混合したものを試験試料として用意した。このような粉砕・混合により、ルチン・クロロゲン酸・スコポレチンが原料中に一様に分布せず局在しているようなケースであっても、測定の精度を担保するようにした。
【００５３】
また、試験試料は事前に水分量を一定化するために蔵置されたものを用いた。例えばたばこ原料の場合、調和条件（22度60%の室内）で２４時間以上蔵置することが好ましい。このように測定前に水分量を一定化させておくことで蛍光のピークシフトが起こりにくくなる。
【００５４】
＜ルチン含有量の測定にたばこ原料を試料として用いる場合＞
ルチン含有量が既知の各サンプルについては、粒径１ｍｍ以下に粉砕し十分に混合したものを試験試料として用意した。
【００５５】
試料がたばこ原料であった場合、ルチンはたばこ原料中に均一に存在するのではなく局在することがわかっていることから、このように、測定前に試料を一定粒径（１ｍｍ径）以下に粉砕し、十分に混合してから蛍光指紋を取得することが好ましいものである。なお、各サンプルのルチン含有量については、事前に高速液体クロマトグラフ（HPLC）により定量しておいた。
【００５６】
また、試験試料は事前に水分量を一定化するために蔵置されたものを用いた。たばこ原料の場合は、調和条件（22度60%の室内）で２４時間以上蔵置することが好ましい。このように測定前に水分量を一定化させておくことで蛍光のピークシフトが起こりにくくなる。
【００５７】
但し、ルチンに関しては溶かせる溶媒が限られているため、固体のまま測定するのが望ましい場合もあり、たばこ原料中のルチンの含有量を測定（推定）するために固体のままで蛍光指紋情報を取得する場合もある。
【００５８】
〔蛍光指紋情報の取得〕
試験試料の蛍光指紋情報を取得するために、蛍光指紋測定装置として、日立ハイテクサイエンス社製F-7000を用い、反射法（FrontFace）により測定を行った。
【００５９】
測定条件は、励起光２００−６００nm, 蛍光２００−７００nm, 分解能５nm, スリット幅５nm，フォトマル感度７００Vであった。なお、分解能５nmを考慮すれば、測定波長は少なくとも５nm程度の誤差を許容するものである。
【００６０】
〔蛍光指紋情報に対する前処理〕
＜ルチン・クロロゲン酸・スコポレチンを標品で混合した試料を用いる場合＞
ルチン・クロロゲン酸・スコポレチンの標品を様々な割合で混合した試験試料から取得された蛍光指紋情報に対して前処理を行う。この前処理には、例えば、MatlabやPLS_toolbox等の専用ソフトウエアが使用される。前処理としては、前述の＜蛍光指紋情報に対する前処理＞において詳述した本発明に特有な２次微分処理の外に、好ましくは、各スペクトルに対する従前使用されている２次微分処理や、成分情報に寄与しない波長を除去する処理が挙げられる。因みに、成分情報に寄与しない波長を除去する処理として、例えば、以下のような手法を採用し得るが、各処理手法自体は既知であり、その詳細についての説明は省く。
（ａ）Variable important projection(VIP)
（ｂ）interval PLS(iPLS)
（ｃ）Genetic algorithms(GA)
（ｄ）Jack-knife分析
（ｅ）Forward interval PLS
（ｆ）Backward interval PLS(biPLS)
（ｇ）Synergy interval PLS(siPLS)
（ｈ）LASSO type method
【００６１】
本発明に特有な２次微分処理に関しては、前述の＜蛍光指紋情報に対する前処理＞において詳述した手法を適用するが、本混合試料のケースでみれば、ルチン、クロロゲン酸、スコポレチンの蛍光指紋情報に対して、ｘ軸（蛍光波長軸）に平行な軸に沿った２次微分処理を行うことにより、有効な結果を得ることができる。
【００６２】
＜ルチン含有量の測定のためにたばこ原料を試料として用いる場合＞
基本的に、前述したルチン・クロロゲン酸・スコポレチンを標品で混合した試料を用いる場合と同様の前処理を行えばよい。
【００６３】
〔検量線の作成・検証〕
検量線は、具体的には、取得された蛍光指紋情報を説明変数、既知の成分（ルチン、クロロゲン酸、スコポレチン）の含有量を目的変数とし、ＰＬＳ回帰分析（以下、単に「ＰＬＳ」ということもある）を使用して作成する。
【００６４】
ＰＬＳ回帰分析に関し、その概略を簡潔に説明しておく。
ＰＬＳでは、説明変数Ｘ（行列）と目的変数ｙ（ベクトル）は、以下の二つの基本式（１）、（２）を満たしている。
Ｘ＝ＴＰT＋Ｅ        （１）
ｙ＝Ｔｑ＋ｆ          （２）
ここで、Ｔは潜在変数（行列）、Ｐはローディング（行列）、Ｅは説明変数Ｘの残差（行列）、ｑは係数（ベクトル）、ｆは目的変数の残差（ベクトル）、ＰTはＰの転置行列である。
【００６５】
因みに、ＰＬＳは、説明変数Ｘの情報を目的変数ｙのモデリングに直接用いるのではなく、説明変数Ｘの情報の一部を潜在変数ｔに変換し、潜在変数ｔを用いて目的変数ｙをモデリングするものである。なお、潜在変数の数は、例えば、クロスバリデーションによる予測的説明分散値を指標として決定することができる。また、潜在変数は、主成分と呼ばれることもある。
【００６６】
特に、１成分モデルの場合には、上記の（１）、（２）は、以下の（３）、（４）で表される。
Ｘ＝ｔ1ｐ1T＋Ｅ        （３）
ｙ＝ｔ1ｑ1＋ｆ         （４）
ここで、ｔ1は潜在変数（ベクトル）、ｐ1はローディング（ベクトル）、ｑ1は係数（スカラー）である。
【００６７】
今、ｔ1がＸの線形結合で表されると仮定すると、以下の（５）が成立する。
ｔ1＝Ｘｗ1          （５）
ここで、ｗ1は規格化された重みベクトルである。
【００６８】
ＰＬＳは、ｙとｔ1との共分散ｙTｔ1を、ｗ1のノルムが１（｜ｗ1｜＝１）という条件下で最大化するようなｔ1を求めるものであり、ｔ1の算出には、所謂ラグランジュの未定乗数法を用いればよい。ラグランジュの未定乗数法を用いた計算手法は周知であるから、計算の詳細は省略し、ｗ1、ｐ1、ｑ1に関する計算結果のみ、以下の（６）〜（８）として示す。  
ｗ1＝ＸTｙ／｜ＸTｙ｜    （６）
ｐ1＝ＸTｔ1／ｔ1Tｔ1      （７）
ｑ1＝ｙＴｔ１／ｔ1Tｔ1    （８）
なお、（７）、（８）式のｔ1は、（６）式で求めたｗ1を（５）式に代入することにより算出されたベクトルである。
【００６９】
多成分モデルについても、同様の手法で計算できるが、計算手法は周知であるから詳細は省略する。
【００７０】
まず、ルチン・クロロゲン酸・スコポレチンを標品で混合した試料を用いる場合について、検量線の作成・検証の詳細を説明する。
【００７１】
上述のような手法で得られた検量線を、ルチン、クロロゲン酸、スコポレチンの各成分について作成・検証するために、ルチン、クロロゲン酸、スコポレチンの各成分の含有量が既知の複数のサンプルを、検量線の作成に使用するためのキャリブレーション用サンプル群と、検量線を検証して有効性を確認するためのバリデーション用サンプル群に分けて用意する。
【００７２】
キャリブレーション用サンプル群に対して上述のＰＬＳ回帰分析を適用し、取得された蛍光指紋情報から各成分について含有量を推定する検量線を作成する。なお、検量線の作成に当たり、取得された蛍光指紋に対する前処理を省くこともできるが、（１）本発明に特有な２次微分処理、（２）成分情報に寄与しない波長に対する除去処理等の前処理を行うことが望ましい。（２）の前処理に関しては、例えば、以下のような処理を採用することができる。
・非蛍光成分の除去，散乱光の除去，低感度領域の除去
・対数変換（Log10）→従前の二次微分→規格化（normalize）→オートスケーリング（autoscale）
・VIPによる波長限定
因みに、上記（１）、（２）の前処理の順序は、適宜決めることができるが、（２）を先行させることが望ましい。
【００７３】
次に、バリデーション用サンプル群について、取得された蛍光指紋情報から前記検量線を用いて各成分の含有量を推定し、検量線の検証を行う。
【００７４】
表２は、ルチン、クロロゲン酸、スコポレチンの各成分について、キャリブレーション用サンプル群の決定係数Ｒ2及び検量線作成時の標準誤差ＳＥＣ、並びに、バリデーション用サンプル群の決定係数Ｒ2及び検量線評価（検証）時の標準誤差ＳＥＰを纏めたものである。
【００７５】
【表２】
2020031447.tif 000004
【００７６】
図７Ａ〜図7Ｃは、ルチン、クロロゲン酸、スコポレチンの各成分について、バリデーション用サンプル群のデータをグラフ化したものであり、図７Ａはスコポレチン、図７Ｂはクロロゲン酸、図７Ｃはルチンをそれぞれ表している。
【００７７】
上述の表２及び図７Ａ〜図7Ｃからも明らかなように、ルチン、クロロゲン酸、スコポレチンの検量線は、良好な推定精度を有しており、当該検量線の有効性が確認されている。
【００７８】
次に、ルチン含有量の測定のためにたばこ原料を試料として用いる場合について、キャリブレーション及びバリデーションの推定精度について詳述する。
【００７９】
図８Ａは、横軸に高速液体クロマトグラフ（HPLC）による実測値（化学分析値）、縦軸に蛍光指紋による推定値を取り、バリデーション用サンプル群に属する各サンプルについて、前記前処理を行った場合の対応する点をプロットしたグラフである。
【００８０】
前述のような、ルチン・クロロゲン酸・スコポレチンを標品で混合した試料を用いる場合と同様の前処理を行った場合のキャリブレーション用サンプル群において、決定係数Ｒ2＝０．９７（ＳＥＣ＝０．０８６％）であり、化学分析値と検量線による推定値との間に高い相関を有し、良好な推定精度であることが確認されている。また、図８Ａによれば、バリデーション用サンプル群における推定精度は、決定係数Ｒ2＝０．９１（ＳＥＰ＝０．１９％）あり、当該検量線の有効性が確認されている。
【００８１】
なお、前述のような前処理を省いた場合には、キャリブレーション用サンプル群における推定精度は、決定係数Ｒ2＝０．９３（ＳＥＣ＝０．１６％）、バリデーション用サンプル群における推定精度は、図９に示されるように、決定係数Ｒ2＝０．８７（ＳＥＰ＝０．２２％）であり、化学分析値と検量線による推定値との間に高い相関を有してはいるものの、前述のような前処理を行った場合と比較すると、推定精度について若干の低下が認められる。
【００８２】
次に、クロロゲン酸含有量の測定のためにたばこ原料を試料として用いる場合について、キャリブレーション及びバリデーションの推定精度について詳述する。
【００８３】
図８Ｂは、横軸に高速液体クロマトグラフ（HPLC）による実測値（化学分析値）、縦軸に蛍光指紋による推定値を取り、バリデーション用サンプル群に属する各サンプルについて、前記前処理を行った場合の対応する点をプロットしたグラフである。
【００８４】
前述のような、ルチン・クロロゲン酸・スコポレチンを標品で混合した試料を用いる場合と同様の前処理を行った場合のキャリブレーション用サンプル群において、決定係数Ｒ2＝０．８７（ＳＥＣ＝０．１９％）であり、化学分析値と検量線による推定値との間に高い相関を有し、良好な推定精度であることが確認されている。また、図８Ｂによれば、バリデーション用サンプル群における推定精度は、決定係数Ｒ2＝０．８８（ＳＥＰ＝０．２０％）あり、当該検量線の有効性が確認されている。
【００８５】
〔未知試料における各成分の含有量の推定〕
有効性が確認された検量線を用いて、ルチン、クロロゲン酸、スコポレチンの各成分の含有量が未知の試料の蛍光指紋情報に基づき、前記試料に含有される各成分の含有量を推定する。
【００８６】
なお、各成分の含有量が未知の試料について、取得された蛍光指紋に対する前処理を省くこともできるが、検量線を取得したときと同様、（１）本発明に特有な２次微分処理、（２）成分情報に寄与しない波長に対する除去処理等の前処理を行い、次に、処理後の蛍光指紋から検量線に基づいて未知試料の各成分の含有量を推定することが望ましい。なお、上記（１）、（２）の前処理の順序は、適宜決めることができるが、（２）を先行させることが望ましい。
【００８７】
各成分以外のポリフェノール成分をノイズとして有効に除去し、各成分に特徴的に強い蛍光を示す特定の励起／蛍光波長（各成分の化学構造に基づいて決まる蛍光強度が最大になる波長条件）を使用することにより、蛍光指紋法による各成分の的確な定量が実現できる。
【００８８】
そして、このように有効性が確認された検量線を用いて、ルチン、クロロゲン酸、スコポレチンの各成分の含有量が未知の試料の蛍光指紋情報に基づき、前記試料に含有される各成分の含有量を有効に推定できる。
なお、たばこ原料等のルチンを含有する未知試料中のルチンの定量に関しては、以下に詳述するような簡便化した手法を有効に利用できる。
【００８９】
＜ルチンに関する簡便化した適用例＞
これまでの説明においては、基本的に、全ての蛍光指紋情報を使用してルチン含有量を推定する態様を説明してきたが、精度は多少犠牲にしても、安価なハンディタイプのデバイスへの適用や簡素化された測定工程による測定の迅速化等に対する要請が存在する。
そこで、このような簡易測定という要請に応えるための簡便化した態様を以下説明する。
【００９０】
この簡素化した態様は、励起波長／蛍光波長（測定波長）として、２８５／４１０，４１５，４２０，４２５，４３０，４３５，４４０，４５０，４６０（nm）の９波長のみを用いて蛍光指紋情報を取得し、ルチンの定量を行うものである。定量のための手法は、基本的に、全ての蛍光指紋情報を使用するケースと同様であり、詳細は省略することとするが、このように波長を限定した（概ね１０波長未満の）場合には、検量線の作成に際し、PLS回帰以外に、重回帰分析（MLR）も利用可能である。なお、波長を限定しない場合にはPLS回帰を使用することが望ましい。
【００９１】
なお、２８５／４１０，４１５，４２０，４２５，４３０，４３５，４４０，４５０，４６０（nm）の９波長は、ルチンの化学構造に基づいて蛍光強度が最大値を取る特定の励起／蛍光波長に相当するものである。
【００９２】
このような特定数の波長を使用する簡便化した態様については、キャリブレーション用サンプル群における推定精度は、決定係数Ｒ2＝０．７３（ＳＥＣ＝０．２５％）であった。また、バリデーション用サンプル群における推定精度は、決定係数Ｒ2＝０．８２（ＳＥＰ＝０．２６％）であり、前述の態様と比較すると、推定精度の低下が認められるものの、簡易測定という要請に十分応え得るレベル範囲にあるものといえる。
【００９３】
また、簡便化した態様はこれに限られるものではなく、励起波長／蛍光波長（測定波長）として、２８５／４１０，４１５，４２０，４２５，４３０，４３５，４４０，４５０，４６０（nm）のうち少なくとも１つを含むような態様を採用することができる。この場合には、９波長のみを用いるケースと比較して、測定精度の点では劣るものの、迅速性の向上を実現できる。
【００９４】
そして、上述のような簡便化した態様を、ハンディタイプのデバイスや、簡素化された測定工程として実現することにより、安価且つ迅速に測定結果を取得することができる。
【００９５】
なお、上述のような簡便化した適用例に関しては、迅速な定量を行うために、取得された蛍光指紋に対する前処理を省くこともできる。
【００９６】
また、上述のような簡便化した適用例は、たばこ原料中のルチンの含有量の推定のみに有効なものではなく、ルチンを含有する他の材料に対しても有効に適用することができる。
【００９７】
本発明は、上述した実施の態様以外にも、特許請求の範囲に記載した技術的思想の範囲内において種々の異なる実施の態様を採用し得ることに留意されたい。
【符号の説明】
【００９８】
１００：蛍光指紋分析による試料の評価・推定装置
１１０：前処理手段
１２０：推定モデル作成手段
１３０：成分量推定手段
【図１】
2020031447.tif 000005
【図２】
2020031447.tif 000006
【図３】
2020031447.tif 000007
【図４】
2020031447.tif 000008
【図５】
2020031447.tif 000009
【図６Ａ】
2020031447.tif 000010
【図６Ｂ】
2020031447.tif 000011
【図６Ｃ】
2020031447.tif 000012
【図７Ａ】
2020031447.tif 000013
【図７Ｂ】
2020031447.tif 000014
【図７Ｃ】
2020031447.tif 000015
【図８Ａ】
2020031447.tif 000016
【図８Ｂ】
2020031447.tif 000017
【図９】
2020031447.tif 000018
【図１０】
2020031447.tif 000019
【国際調査報告】
2020031447.tif 000020
2020031447.tif 000021
2020031447.tif 000022
2020031447.tif 000023