(57)【要約】【課題】目的の層以外の層によるノイズの影響を軽減する光吸収係数算出装置、濃度定量装置、光吸収係数算出方法、濃度定量方法、光吸収係数の算出を行うプログラム及び濃度の算出を行うプログラムを提供する。【解決手段】照射手段１０２と、受光手段１０３と、光強度取得手段１０４と、複数の観測対象のモデルにおける短時間パルス光の時間分解波形の所定の時刻における光強度と、複数の観測対象のモデルにおける短時間パルス光の時間分解波形の所定の時刻における無吸収時の光強度と、複数の観測対象のモデルの各々について設定された複数の光散乱媒質の層の各々の層の光吸収係数のモデル値と、に基づいて作成された検量線を取得し、光強度取得手段１０４が取得した光強度を検量線と照合して、任意の層における光吸収係数を算出する光吸収係数算出手段１０６と、を含む。【選択図】図１
2014006063.tif 000002
【特許請求の範囲】
【請求項１】
複数の光散乱媒質の層から形成される観測対象のうち、任意の層における光吸収係数を算出する光吸収係数算出装置であって、
前記観測対象に短時間パルス光を照射する照射手段と、
前記短時間パルス光が前記観測対象によって後方散乱した光を受光する受光手段と、
前記照射手段が短時間パルス光を照射した時刻以降の所定の時刻において前記受光手段が受光した光の強度を取得する光強度取得手段と、
複数の観測対象のモデルにおける前記短時間パルス光の時間分解波形の前記所定の時刻における光強度と、前記複数の観測対象のモデルにおける前記短時間パルス光の時間分解波形の前記所定の時刻における無吸収時の光強度と、前記複数の観測対象のモデルの各々について設定された前記複数の光散乱媒質の層の各々の層の光吸収係数のモデル値と、に基づいて作成された検量線を取得し、前記光強度取得手段が取得した光強度を前記検量線と照合して、前記任意の層における光吸収係数を算出する光吸収係数算出手段と、
を含むことを特徴とする光吸収係数算出装置。
【請求項２】
前記光吸収係数算出手段は、
前記観測対象のモデルの数を示すｋ、前記観測対象から得られる信号を形成する主成分の個数を示すｉ、前記観測対象の層数を示すｊ、ｋ番目の前記観測対象のモデルにおける前記短時間パルス光の入射光強度に対する時間分解波形の時刻ｔにおける光強度を示すＩ’ｓｋ（ｔ）、前記複数の観測対象のモデルに対して共通の、前記短時間パルス光の入射光強度に対する時間分解波形の時刻ｔにおける無吸収時の光強度を示すＮ’ｓ（ｔ）、ｋ番目の前記観測対象のモデルにおける第ｉ主成分の主成分分析スコアを示すαｋｉ、第ｉ主成分の主成分分析ローディングを示すＰｉ（ｔ）、ｋ番目の前記観測対象のモデルにおける第ｊ層の光吸収係数のモデル値を示すμａｋｊ、前記光吸収係数のモデル値μａｋｊと前記主成分分析スコアαｋｉとの対応関係を示す前記検量線としてのＰＩＬＳｉｊ、前記観測対象から得られる信号の第ｉ主成分の主成分スコアを示すαｍｉ、前記受光手段が時刻ｔにおいて受光した入射光強度に対する光強度を示すＩ’ｍ（ｔ）、前記観測対象の第ｊ層の光吸収係数を示すμａｍｊ、を用いて、
【数１】
2014006063.tif 000015
から前記任意の層における光吸収係数を算出する請求項１に記載の光吸収係数算出装置。
【請求項３】
請求項１または２に記載の光吸収係数算出装置と、
前記光吸収係数算出手段で算出した光吸収係数に基づいて、前記任意の層における目的物質の濃度を算出する濃度算出手段と、
を含むことを特徴とする濃度定量装置。
【請求項４】
前記濃度算出手段は、前記観測対象が皮膚であり、前記任意の層が真皮層であるときに、当該真皮層に含まれるグルコースの濃度を定量することを特徴とする請求項３に記載の濃度定量装置。
【請求項５】
複数の光散乱媒質の層から形成される観測対象のうち、任意の層における光吸収係数を算出する光吸収係数算出方法であって、
前記観測対象に短時間パルス光を照射し、
前記短時間パルス光が前記観測対象によって後方散乱した光を受光し、
前記照射手段が短時間パルス光を照射した時刻以降の所定の時刻において前記受光手段が受光した光の強度を取得し、
複数の観測対象のモデルにおける前記短時間パルス光の時間分解波形の前記所定の時刻における光強度と、前記複数の観測対象のモデルにおける前記短時間パルス光の時間分解波形の前記所定の時刻における無吸収時の光強度と、前記複数の観測対象のモデルの各々について設定された前記複数の光散乱媒質の層の各々の層の光吸収係数のモデル値と、に基づいて作成された検量線を取得し、前記光強度取得手段が取得した光強度を前記検量線と照合して、前記任意の層における光吸収係数を算出することを特徴とする光吸収係数算出方法。
【請求項６】
前記観測対象のモデルの数を示すｋ、前記観測対象から得られる信号を形成する主成分の個数を示すｉ、前記観測対象の層数を示すｊ、ｋ番目の前記観測対象のモデルにおける前記短時間パルス光の入射光強度に対する時間分解波形の時刻ｔにおける光強度を示すＩ’ｓｋ（ｔ）、前記複数の観測対象のモデルに対して共通の、前記短時間パルス光の入射光強度に対する時間分解波形の時刻ｔにおける無吸収時の光強度を示すＮ’ｓ（ｔ）、ｋ番目の前記観測対象のモデルにおける第ｉ主成分の主成分分析スコアを示すαｋｉ、第ｉ主成分の主成分分析ローディングを示すＰｉ（ｔ）、ｋ番目の前記観測対象のモデルにおける第ｊ層の光吸収係数のモデル値を示すμａｋｊ、前記光吸収係数のモデル値μａｋｊと前記主成分分析スコアαｋｉとの対応関係を示す前記検量線としてのＰＩＬＳｉｊ、前記観測対象から得られる信号の第ｉ主成分の主成分スコアを示すαｍｉ、前記受光手段が時刻ｔにおいて受光した入射光強度に対する光強度を示すＩ’ｍ（ｔ）、前記観測対象の第ｊ層の光吸収係数を示すμａｍｊ、を用いて、
【数２】
2014006063.tif 000016
から前記任意の層における光吸収係数を算出する請求項５に記載の光吸収係数算出方法。
【請求項７】
請求項５または６に記載の光吸収係数算出方法で算出した光吸収係数に基づいて、前記任意の層における目的物質の濃度を算出することを特徴とする濃度定量方法。
【請求項８】
前記観測対象が皮膚であり、前記任意の層が真皮層であるときに、当該真皮層に含まれるグルコースの濃度を定量することを特徴とする請求項７に記載の濃度定量方法。
【請求項９】
複数の光散乱媒質の層から形成される観測対象のうち、任意の層における光吸収係数の算出を行うプログラムであって、
前記観測対象に短時間パルス光を照射し、
前記短時間パルス光が前記観測対象によって後方散乱した光を受光し、
前記照射手段が短時間パルス光を照射した時刻以降の所定の時刻において前記受光手段が受光した光の強度を取得し、
複数の観測対象のモデルにおける前記短時間パルス光の時間分解波形の前記所定の時刻における光強度と、前記複数の観測対象のモデルにおける前記短時間パルス光の時間分解波形の前記所定の時刻における無吸収時の光強度と、前記複数の観測対象のモデルの各々について設定された前記複数の光散乱媒質の層の各々の層の光吸収係数のモデル値と、に基づいて作成された検量線を取得し、前記光強度取得手段が取得した光強度を前記検量線と照合して、前記任意の層における光吸収係数の算出を行うプログラム。
【請求項１０】
前記観測対象のモデルの数を示すｋ、前記観測対象から得られる信号を形成する主成分の個数を示すｉ、前記観測対象の層数を示すｊ、ｋ番目の前記観測対象のモデルにおける前記短時間パルス光の入射光強度に対する時間分解波形の時刻ｔにおける光強度を示すＩ’ｓｋ（ｔ）、前記複数の観測対象のモデルに対して共通の、前記短時間パルス光の入射光強度に対する時間分解波形のモデルの時刻ｔにおける無吸収時の光強度を示すＮ’ｓ（ｔ）、ｋ番目の前記観測対象のモデルにおける第ｉ主成分の主成分分析スコアを示すαｋｉ、第ｉ主成分の主成分分析ローディングを示すＰｉ（ｔ）、ｋ番目の前記観測対象のモデルにおける第ｊ層の光吸収係数のモデル値を示すμａｋｊ、前記光吸収係数のモデル値μａｋｊと前記主成分分析スコアαｋｉとの対応関係を示す前記検量線としてのＰＩＬＳｉｊ、前記観測対象から得られる信号の第ｉ主成分の主成分スコアを示すαｍｉ、前記受光手段が時刻ｔにおいて受光した入射光強度に対する光強度を示すＩ’ｍ（ｔ）、前記観測対象の第ｊ層の光吸収係数を示すμａｍｊ、を用いて、
【数３】
2014006063.tif 000017
から前記任意の層における光吸収係数を算出する請求項９に記載の光吸収係数の算出を行うプログラム。
【請求項１１】
請求項９または１０に記載の光吸収係数の算出を行うプログラムで算出した光吸収係数に基づいて、前記任意の層における目的物質の濃度の算出を行うプログラム。
【請求項１２】
前記観測対象が皮膚であり、前記任意の層が真皮層であるときに、当該真皮層に含まれるグルコースの濃度を定量することを特徴とする請求項１１に記載の濃度の算出を行うプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明のいくつかの態様は、複数の光散乱媒質の層から形成される観測対象のうち、任意の層における光吸収係数を算出する光吸収係数算出装置、濃度定量装置、光吸収係数算出方法、濃度定量方法、光吸収係数の算出を行うプログラム及び濃度の算出を行うプログラムに関する。
【背景技術】
【０００２】
従来、血糖値の測定は、指先などから採血を行い、血中のグルコースに対する酵素活性を測ることで行っていた。しかし、このような血糖値の測定方法は、指先などから血液を採取して測定しなければならず、採血に手間と痛みを伴うことや、血液を付着させる測定チップが必要なことから、採血を必要としない非侵襲型の血糖値の測定方法が望まれている。
【０００３】
そこで、皮膚に近赤外光を照射し、その光吸収量からグルコースの濃度を求める方法が検討されている（例えば、特許文献１を参照）。具体的には、予めグルコース濃度と照射する光の波長と光の吸収量との関係を示す検量線を作成しておき、モノクロメーター等を用いてある波長域を走査し、その波長域の各波長に対する吸収量を求め、当該波長及び吸収量と検量線とを比較することでグルコース濃度を算出する。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許第３９３１６３８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、従来の非侵襲血糖値測定方法は、光の入出射間距離を定めることによって、真皮層の近赤外スペクトルを測定するため、測定したスペクトルに、真皮層のスペクトルのみならず表皮層や皮下組織層のスペクトルも含まれ、観測される光吸収係数の変化には表皮層や皮下組織層によるノイズが含まれてしまうという問題があった。
【０００６】
本発明は上記の点に鑑みてなされたものであり、その目的は、目的の層以外の層によるノイズの影響を軽減する光吸収係数算出装置、濃度定量装置、光吸収係数算出方法、濃度定量方法、光吸収係数の算出を行うプログラム及び濃度の算出を行うプログラムを提供することにある。
【課題を解決するための手段】
【０００７】
本発明のいくつかの態様は上記の課題を解決するためになされたものであり、複数の光散乱媒質の層から形成される観測対象のうち、任意の層における光吸収係数を算出する光吸収係数算出装置であって、前記観測対象に短時間パルス光を照射する照射手段と、前記短時間パルス光が前記観測対象によって後方散乱した光を受光する受光手段と、前記照射手段が短時間パルス光を照射した時刻以降の所定の時刻において前記受光手段が受光した光の強度を取得する光強度取得手段と、複数の観測対象のモデルにおける前記短時間パルス光の時間分解波形の前記所定の時刻における光強度と、前記複数の観測対象のモデルにおける前記短時間パルス光の時間分解波形の前記所定の時刻における無吸収時の光強度と、前記複数の観測対象のモデルの各々について設定された前記複数の光散乱媒質の層の各々の層の光吸収係数のモデル値と、に基づいて作成された検量線を取得し、前記光強度取得手段が取得した光強度を前記検量線と照合して、前記任意の層における光吸収係数を算出する光吸収係数算出手段と、を含むことを特徴とする。
【０００８】
この構成によれば、任意の層の光吸収係数を選択的に算出することができる。そのため、他の層によるノイズの影響を低減し、精度の高い濃度の定量を行うことができる。
【０００９】
また、本発明のいくつかの態様は、前記光吸収係数算出手段は、前記観測対象のモデルの数を示すｋ、前記観測対象から得られる信号を形成する主成分の個数を示すｉ、前記観測対象の層数を示すｊ、ｋ番目の前記観測対象のモデルにおける前記短時間パルス光の入射光強度に対する時間分解波形の時刻ｔにおける光強度を示すＩ’ｓｋ（ｔ）、前記複数の観測対象のモデルに対して共通の、前記短時間パルス光の入射光強度に対する時間分解波形の時刻ｔにおける無吸収時の光強度を示すＮ’ｓ（ｔ）、ｋ番目の前記観測対象のモデルにおける第ｉ主成分の主成分分析スコアを示すαｋｉ、第ｉ主成分の主成分分析ローディングを示すＰｉ（ｔ）、ｋ番目の前記観測対象のモデルにおける第ｊ層の光吸収係数のモデル値を示すμａｋｊ、前記光吸収係数のモデル値μａｋｊと前記主成分分析スコアαｋｉとの対応関係を示す前記検量線としてのＰＩＬＳｉｊ、前記観測対象から得られる信号の第ｉ主成分の主成分スコアを示すαｍｉ、前記受光手段が時刻ｔにおいて受光した入射光強度に対する光強度を示すＩ’ｍ（ｔ）、前記観測対象の第ｊ層の光吸収係数を示すμａｍｊ、を用いて、下記の（１）式〜（４）式から前記任意の層における光吸収係数を算出することを特徴とする。
【００１０】
【数１】
2014006063.tif 000003
【００１１】
この構成によれば、任意の層の光吸収係数を高精度に算出することができる。
【００１２】
また、本発明のいくつかの態様は、前記光吸収係数算出装置と、前記光吸収係数算出手段で算出した光吸収係数に基づいて、前記任意の層における目的物質の濃度を算出する濃度算出手段と、を含むことを特徴とする。
【００１３】
この構成によれば、他の層によるノイズの影響を低減し、精度の高い濃度の定量を行うことができる。
【００１４】
また、本発明のいくつかの態様は、前記濃度算出手段は、前記観測対象が皮膚であり、前記任意の層が真皮層であるときに、当該真皮層に含まれるグルコースの濃度を定量することを特徴とする。
【００１５】
この構成によれば、他の層によるノイズの影響を低減し、グルコースの濃度の定量を高精度で行うことができる。
【００１６】
また、本発明のいくつかの態様は、複数の光散乱媒質の層から形成される観測対象のうち、任意の層における光吸収係数を算出する光吸収係数算出方法であって、前記観測対象に短時間パルス光を照射し、前記短時間パルス光が前記観測対象によって後方散乱した光を受光し、前記照射手段が短時間パルス光を照射した時刻以降の所定の時刻において前記受光手段が受光した光の強度を取得し、複数の観測対象のモデルにおける前記短時間パルス光の時間分解波形の前記所定の時刻における光強度と、前記複数の観測対象のモデルにおける前記短時間パルス光の時間分解波形の前記所定の時刻における無吸収時の光強度と、前記複数の観測対象のモデルの各々について設定された前記複数の光散乱媒質の層の各々の層の光吸収係数のモデル値と、に基づいて作成された検量線を取得し、前記光強度取得手段が取得した光強度を前記検量線と照合して、前記任意の層における光吸収係数を算出することを特徴とする。
【００１７】
この方法によれば、任意の層の光吸収係数を選択的に算出することができる。そのため、他の層によるノイズの影響を低減し、精度の高い濃度の定量を行うことができる。
【００１８】
また、本発明のいくつかの態様は、前記観測対象のモデルの数を示すｋ、前記観測対象から得られる信号を形成する主成分の個数を示すｉ、前記観測対象の層数を示すｊ、ｋ番目の前記観測対象のモデルにおける前記短時間パルス光の入射光強度に対する時間分解波形の時刻ｔにおける光強度を示すＩ’ｓｋ（ｔ）、前記複数の観測対象のモデルに対して共通の、前記短時間パルス光の入射光強度に対する時間分解波形の時刻ｔにおける無吸収時の光強度を示すＮ’ｓ（ｔ）、ｋ番目の前記観測対象のモデルにおける第ｉ主成分の主成分分析スコアを示すαｋｉ、第ｉ主成分の主成分分析ローディングを示すＰｉ（ｔ）、ｋ番目の前記観測対象のモデルにおける第ｊ層の光吸収係数のモデル値を示すμａｋｊ、前記光吸収係数のモデル値μａｋｊと前記主成分分析スコアαｋｉとの対応関係を示す前記検量線としてのＰＩＬＳｉｊ、前記観測対象から得られる信号の第ｉ主成分の主成分スコアを示すαｍｉ、前記受光手段が時刻ｔにおいて受光した入射光強度に対する光強度を示すＩ’ｍ（ｔ）、前記観測対象の第ｊ層の光吸収係数を示すμａｍｊ、を用いて、下記の（１）式〜（４）式から前記任意の層における光吸収係数を算出することを特徴とする。
【００１９】
【数２】
2014006063.tif 000004
【００２０】
この方法によれば、任意の層の光吸収係数を高精度に算出することができる。
【００２１】
また、本発明のいくつかの態様は、前記光吸収係数算出手段で算出した光吸収係数に基づいて、前記任意の層における目的物質の濃度を算出することを特徴とする。
【００２２】
この方法によれば、他の層によるノイズの影響を低減し、精度の高い濃度の定量を行うことができる。
【００２３】
また、本発明のいくつかの態様は、前記観測対象が皮膚であり、前記任意の層が真皮層であるときに、当該真皮層に含まれるグルコースの濃度を定量することを特徴とする。
【００２４】
この方法によれば、他の層によるノイズの影響を低減し、グルコースの濃度の定量を高精度で行うことができる。
【００２５】
また、本発明のいくつかの態様は、複数の光散乱媒質の層から形成される観測対象のうち、任意の層における光吸収係数の算出を行うプログラムであって、前記観測対象に短時間パルス光を照射し、前記短時間パルス光が前記観測対象によって後方散乱した光を受光し、前記照射手段が短時間パルス光を照射した時刻以降の所定の時刻において前記受光手段が受光した光の強度を取得し、複数の観測対象のモデルにおける前記短時間パルス光の時間分解波形の前記所定の時刻における光強度と、前記複数の観測対象のモデルにおける前記短時間パルス光の時間分解波形の前記所定の時刻における無吸収時の光強度と、前記複数の観測対象のモデルの各々について設定された前記複数の光散乱媒質の層の各々の層の光吸収係数のモデル値と、に基づいて作成された検量線を取得し、前記光強度取得手段が取得した光強度を前記検量線と照合して、前記任意の層における光吸収係数の算出を行うことを特徴とする。
【００２６】
このプログラムによれば、任意の層の光吸収係数を選択的に算出することができる。そのため、他の層によるノイズの影響を低減し、精度の高い濃度の定量を行うことができる。
【００２７】
また、本発明のいくつかの態様は、前記観測対象のモデルの数を示すｋ、前記観測対象から得られる信号を形成する主成分の個数を示すｉ、前記観測対象の層数を示すｊ、ｋ番目の前記観測対象のモデルにおける前記短時間パルス光の入射光強度に対する時間分解波形の時刻ｔにおける光強度を示すＩ’ｓｋ（ｔ）、前記複数の観測対象のモデルに対して共通の、前記短時間パルス光の入射光強度に対する時間分解波形のモデルの時刻ｔにおける無吸収時の光強度を示すＮ’ｓ（ｔ）、ｋ番目の前記観測対象のモデルにおける第ｉ主成分の主成分分析スコアを示すαｋｉ、第ｉ主成分の主成分分析ローディングを示すＰｉ（ｔ）、ｋ番目の前記観測対象のモデルにおける第ｊ層の光吸収係数のモデル値を示すμａｋｊ、前記光吸収係数のモデル値μａｋｊと前記主成分分析スコアαｋｉとの対応関係を示す前記検量線としてのＰＩＬＳｉｊ、前記観測対象から得られる信号の第ｉ主成分の主成分スコアを示すαｍｉ、前記受光手段が時刻ｔにおいて受光した入射光強度に対する光強度を示すＩ’ｍ（ｔ）、前記観測対象の第ｊ層の光吸収係数を示すμａｍｊ、を用いて、下記の（１）式〜（４）式から前記任意の層における光吸収係数を算出することを特徴とする。
【００２８】
【数３】
2014006063.tif 000005
【００２９】
このプログラムによれば、任意の層の光吸収係数を高精度に算出することができる。
【００３０】
また、本発明のいくつかの態様は、前記光吸収係数の算出を行うプログラムで算出した光吸収係数に基づいて、前記任意の層における目的物質の濃度の算出を行うことを特徴とする。
【００３１】
このプログラムによれば、他の層によるノイズの影響を低減し、精度の高い濃度の定量を行うことができる。
【００３２】
また、本発明のいくつかの態様は、前記観測対象が皮膚であり、前記任意の層が真皮層であるときに、当該真皮層に含まれるグルコースの濃度を定量することを特徴とする。
【００３３】
このプログラムによれば、他の層によるノイズの影響を低減し、グルコースの濃度の定量を高精度で行うことができる。
【図面の簡単な説明】
【００３４】
【図１】本発明の一実施形態に係る血糖値測定装置の構成を示す概略ブロック図である。
【図２】シミュレーション部が算出した時間分解波形を示すグラフである。
【図３】皮膚の主要物質の吸収スペクトルを示すグラフである。
【図４】血糖値測定装置が血糖値を測定する動作を示すフローチャートである。
【発明を実施するための形態】
【００３５】
以下、図面を参照して、本発明の実施の形態について説明する。
図１は、本発明の一実施形態に係る血糖値測定装置１００の構成を示す概略ブロック図である。
【００３６】
図１に示すように、血糖値測定装置１００（濃度定量装置）は、光吸収係数算出装置１０１、物質吸収情報記憶部１０７、濃度算出部１０８（濃度算出手段）、濃度単位変換部１０９、及び濃度表示部１１０、を備える。
【００３７】
血糖値測定装置１００は、皮膚（観測対象）の真皮層（任意の層）に含まれるグルコース（目的物質）の濃度を測定する。
【００３８】
光吸収係数算出装置１０１は、照射部１０２（照射手段）、受光部１０３（受光手段）、計測光強度取得部１０４（光強度取得手段）、検量線記憶部１０５、及び光吸収係数算出部１０６（光吸収係数算出手段）、を備える。
【００３９】
シミュレーション部１１１、光強度記憶部１１２、無吸収時光強度記憶部１１３、及び検量線作成部１１４は、血糖値測定装置１００の外部に配置されている。
【００４０】
シミュレーション部１１１は、複数の皮膚モデルの各々について、光を照射するシミュレーションを行う。皮膚モデルとしては、光吸収係数がゼロの皮膚モデル及び所定の光吸収係数を有する皮膚モデルを用いる。なお、所定の光吸収係数を有する皮膚モデルについては、皮膚の各層の光吸収係数の設定値（光吸収係数のモデル値）がばらついた複数の皮膚モデルを準備する。
【００４１】
光強度記憶部１１２は、複数の皮膚モデルの各々について、短時間パルス光の時間分解波形のモデルの所定の時刻における光強度（以下、光強度の時間関数と称することがある）を記憶する。なお、ここで、短時間パルス光とは、パルス幅が１００ｐｓｅｃ程度かそれ以下のパルス光を意味する。また、短時間パルス光として０．１ｐｓｅｃから数ｐｓｅｃの範囲のパルス幅を持つパルス光を用いても良い。
【００４２】
無吸収時光強度記憶部１１３は、複数の皮膚モデルの各々について、短時間パルス光の時間分解波形のモデルの所定の時刻における無吸収時の光強度（以下、無吸収時光強度の時間関数と称することがある）を記憶する。
検量線作成部１１４は、複数の皮膚モデルを用いて検量線を作成する。なお、検量線については後述する。
【００４３】
照射部１０２は、皮膚に対して短時間パルス光を照射する。この照射部１０２が照射する複数の短時間パルス光は、皮膚を構成する主要物質の各々の物質の吸収スペクトル分布の直交性が高くなる波長の光、すなわち、皮膚を構成する主要物質の各々の物質のうち、ある主要物質における特定物質の吸収スペクトルの極大値が他の物質の吸収スペクトルの極大値と大きく異なる波長の光を含んでいる。
【００４４】
受光部１０３は、短時間パルス光が皮膚によって後方散乱した光を受光する。
計測光強度取得部１０４は、受光部１０３が受光した光のある時刻における光強度を取得する。
検量線記憶部１０５は、検量線作成部１１４が作成した検量線を記憶する。
光吸収係数算出部１０６は、短時間パルス光を照射した皮膚の真皮層における光吸収係数を算出する。
【００４５】
物質吸収情報記憶部１０７は、皮膚の主要物質の光吸収係数、またはモル吸光係数を予め記憶する。
濃度算出部１０８は、真皮層に含まれるグルコースの濃度を算出する。
濃度単位変換部１０９は、グルコースの濃度の単位を所望の単位に変換する。
濃度表示部１１０は、グルコースの濃度を表示する。
【００４６】
本実施形態の血糖値測定装置１００においては、照射部１０２は、皮膚に短時間パルス光を照射し、受光部１０３は、短時間パルス光が皮膚によって後方散乱した光を受光し、計測光強度取得部１０４は、時刻ｔにおいて受光部１０３が受光した光の強度を取得する。シミュレーション部１１１は、複数の皮膚モデルの各々について、予め、光強度の時間関数と、無吸収時光強度の時間関数と、を算出する。光強度記憶部１１２は、シミュレーション部１１１が算出した光強度の時間関数を記憶する。無吸収時光強度記憶部１１３は、シミュレーション部１１１が算出した無吸収時光強度の時間関数を記憶する。検量線作成部１１４は、光強度記憶部１１２に記憶された光強度の時間関数と、無吸収時光強度記憶部１１３に記憶された無吸収時光強度の時間関数と、複数の皮膚モデルの各々において設定された皮膚の各々の層の光吸収係数のモデル値と、に基づいて検量線を作成する。検量線記憶部１０５は、検量線作成部１１４が作成した検量線を記憶する。
【００４７】
次に、光吸収係数算出部１０６は、計測光強度取得部１０４が取得した光強度を、検量線記憶部１０５に記憶された検量線と照合して、皮膚の真皮層の光吸収係数を算出する。そして、濃度算出部１０８は、光吸収係数算出部１０６が算出した光吸収係数に基づいて、真皮層におけるグルコースの濃度を算出する。
【００４８】
これにより、真皮層以外の層によるノイズの影響を軽減して、真皮層に含まれるグルコースの濃度を算出することができる。なお、濃度算出部１０８が算出したグルコースの濃度は、濃度単位変換部１０９により所望の単位に変換され、濃度表示部１１０に表示される。
【００４９】
次に、検量線作成部１１４による検量線の作成について説明する。
血糖値測定装置１００は、血糖値を測定する前に、複数の皮膚モデルを用いて検量線を取得しておく必要がある。
【００５０】
検量線の作成は、主成分回帰分析（ＰＣＲ；Principal Component Regression）、部分最小二乗回帰分析（ＰＬＳ；Partial Least Squares）等の多変量解析を用いて作成することが好ましい。本実施形態では、検量線を作成する際に、主成分回帰分析を用いる。なお、これに限らず、部分最小二乗回帰分析を用いてもよい。
【００５１】
初めに、シミュレーション部１１１は、複数の皮膚モデルを生成する。皮膚モデルの生成は、皮膚の各層の光散乱係数、光吸収係数及び厚みを決定することで行う。
【００５２】
シミュレーション部１１１は、皮膚モデルとして、光吸収係数がゼロの皮膚モデル、及び所定の光吸収係数を有する皮膚モデル、を準備する。所定の光吸収係数を有する皮膚モデルとしては、ある光吸収係数の分布を持った複数の皮膚モデルを準備する。複数の皮膚モデルにおいては、皮膚の各々の層（例えば、真皮、表皮、皮下組織）に対して光吸収係数のモデル値が設定され、各層の光吸収係数のモデル値の組み合わせは、皮膚モデルごとに異なるものとされる。
【００５３】
なお、皮膚モデルの数はなるべく多くしたほうが好ましい。皮膚モデルの数を多くしたほうが精度の高い濃度の定量が可能となる。皮膚モデルとして一般的な人の皮膚モデル（代表モデル）のみを考慮した場合、代表モデルに近い皮膚に対しては精度の高い濃度の定量が可能であるが、代表モデルに近くない皮膚に対しては精度の高い濃度の定量が難しくなる。それに対して、ある程度の光吸収係数のばらつきを考慮した複数の皮膚モデルを用いて検量線を作成した場合には、どのような皮膚に対しても比較的高い精度の濃度の定量が可能となる。
【００５４】
本実施形態では、ある光吸収係数の分布を持った複数の皮膚モデルとして、表１に示すように、５つの皮膚モデル（１番目の第１皮膚モデル〜５番目の第５皮膚モデル）を用いて、主成分回帰分析による検量線作成を行う。
【００５５】
【表１】
2014006063.tif 000006
【００５６】
例えば、第１皮膚モデルとしては、表皮層の光吸収係数が０．２０ｍｍ−１、真皮層の光吸収係数が０．５０ｍｍ−１、皮下組織層の光吸収係数が０．２０ｍｍ−１のものを用いる。第２皮膚モデルとしては、表皮層の光吸収数が０．１５ｍｍ−１、真皮層の光吸収係数が０．６５ｍｍ−１、皮下組織層の光吸収係数が０．０５ｍｍ−１のものを用いる。第３皮膚モデルとしては、表皮層の光吸収係数が０．０５ｍｍ−１、真皮層の光吸収係数が０．６０ｍｍ−１、皮下組織層の光吸収係数が０．１０ｍｍ−１のものを用いる。第４皮膚モデルとしては、表皮層の光吸収数が０．０５ｍｍ−１、真皮層の光吸収係数が０．４５ｍｍ−１、皮下組織層の光吸収係数が０．１５ｍｍ−１のものを用いる。第５皮膚モデルとしては、表皮層の光吸収係数が０．１１ｍｍ−１、真皮層の光吸収係数が０．５９７ｍｍ−１、皮下組織層の光吸収係数が０．１０ｍｍ−１のものを用いる。
【００５７】
なお、各層の光散乱係数は生体や皮膚によるばらつきが小さいため、どの皮膚モデルについても同じ値とする。本実施形態において、表皮層の光散乱係数は４．１７ｍｍ−１、真皮層の光散乱係数は１．３９ｍｍ−１、皮下組織層の光散乱係数は０．８８３１ｍｍ−１とする。この場合、短時間パルス光の時間分解波形の所定の時刻における無吸収時の光強度は、どの皮膚モデルについても同じ値（共通）となる。
【００５８】
同様に、各層の厚みは、どの皮膚モデルについても同じ値とする。本実施形態において、表皮層の厚みは略０．３ｍｍ、真皮層の厚みは略１．２ｍｍ、皮下組織層の厚みは略３．０ｍｍとする。
【００５９】
シミュレーション部１１１は、皮膚モデルを生成すると、当該皮膚モデルに光を照射するシミュレーションを行う。このとき、照射部１０２の位置と受光部１０３の位置との間の距離を決定しておく必要がある。シミュレーションは、モンテカルロ法を用いて行うと良い。
【００６０】
シミュレーション部１１１は、照射する光のモデルを光子（光束）とし、当該光子を皮膚モデルに照射する計算を行う。皮膚モデルに照射された光子は、皮膚モデル内を移動する。
【００６１】
図２は、シミュレーション部が算出した時間分解波形の一例を示すグラフである。
図２の横軸は光子の照射からの経過時間を示し、縦軸は受光部１０３が検出した光子数を示している。シミュレーション部１１１は、単位時間毎に受光部１０３に到達した光子の個数を算出することで、複数の皮膚モデルについて、図２に示すような時間分解波形を算出する。
【００６２】
このような処理を複数の皮膚モデルについて行うことにより、シミュレーション部１１１は、複数の波長に対して、複数の皮膚モデルの時間分解波形を算出する。このとき、シミュレーション部１１１は、皮膚の主要物質（水、たんぱく質、脂質、グルコース等）の吸収スペクトルの直交性が高くなる波長について時間分解波形を算出すると良い。
【００６３】
そして、シミュレーション部１１１は、算出した時間分解波形に基づいて、複数の皮膚モデルの各々について、光強度の時間関数と、無吸収時光強度の時間関数と、を算出する。
【００６４】
図３は、皮膚の主要物質の吸収スペクトルを示すグラフである。
図３の横軸は照射する光の波長を示し、縦軸は吸収係数を示している。図３を参照すると、グルコースの吸収係数は、波長が１６００ｎｍのときに極大となり、水の吸収係数は、波長が１４５０ｎｍのときに極大となる。そのため、シミュレーション部１１１は、例えば１４５０ｎｍ、１６００ｎｍといった皮膚の主要物質の吸収スペクトルの直交性が高くなる波長について時間分解波形を算出すると良い。
【００６５】
シミュレーション部１１１は、複数の波長に対する皮膚モデルの時間分解波形を算出し、複数の皮膚モデルにおける光強度の時間関数を算出すると、当該光強度の時間関数の情報を光強度記憶部１１２に記憶させる。また、シミュレーション部１１１は、複数の波長に対する皮膚モデルの無吸収時の時間分解波形を算出し、複数の皮膚モデルにおける無吸収時光強度の時間関数を算出すると、当該無吸収時光強度の時間関数の情報を無吸収時光強度記憶部１１３に記憶させる。
【００６６】
次に、血糖値測定装置１００が血糖値を測定する動作について説明する。
図４は、血糖値測定装置１００が血糖値を測定する動作を示すフローチャートである。
【００６７】
予め、シミュレーション部１１１により、複数の皮膚モデルの各々について、光強度の時間関数と、無吸収時光強度の時間関数と、を算出しておく（ステップＳ１）。
【００６８】
そして、検量線作成部１１４により、光強度記憶部１１２に記憶された複数の皮膚モデルにおける光強度の時間関数と、無吸収時光強度記憶部１１３に記憶された複数の皮膚モデルにおける無吸収時光強度の時間関数と、複数の皮膚モデルの各々において設定された皮膚の各々の層の光吸収係数のモデル値と、に基づいて検量線を作成する。本実施形態において、検量線作成部１１４は、下記の（１）式、（２）式に基づいて、検量線を作成する。
【００６９】
【数４】
2014006063.tif 000007
【００７０】
但し、皮膚モデルの数をｋ（ｋ＝１〜ｘ）、皮膚から得られる信号を形成する主成分の個数をｉ、皮膚の層数をｊ、ｋ番目の皮膚モデルにおける短時間パルス光の入射光強度に対する時間分解波形の時刻ｔにおける光強度（光強度の時間関数）をＩ’ｓｋ（ｔ）、シミュレーション上の入射光強度をＩｓｉｎ、ｋ番目の皮膚モデルにおける時間分解波形の時刻ｔにおける光強度をＩｓｋ（ｔ）、複数の皮膚モデルに対して共通の、短時間パルス光の時間分解波形の入射光強度に対する時刻ｔにおける無吸収時の光強度（無吸収時光強度の時間関数）をＮ’ｓ（ｔ）、時間分解波形の無吸収時の光強度をＮｓ（ｔ）、ｋ番目の皮膚モデルにおける第ｉ主成分の主成分分析スコアをαｋｉ、第ｉ主成分の主成分分析ローディングをＰｉ（ｔ）、ｋ番目の皮膚モデルにおける第ｊ層の光吸収係数のモデル値をμａｋｊ、光吸収係数のモデル値μａｋｊと主成分分析スコアをαｋｉとの対応関係を示す検量線をＰＩＬＳｉｊ、とする。
【００７１】
具体的には、シミュレーション部１１１で算出した、複数の皮膚モデルにおける入射光強度に対する光強度の時間関数Ｉ’ｓｋ（ｔ）及び複数の皮膚モデルにおける入射光強度に対する無吸収時光強度の時間関数Ｎ’ｓ（ｔ）を（１）式の左辺に代入するとともに、第ｉ主成分の主成分分析ローディングＰｉ（ｔ）を（１）式の右辺に代入することにより、ｋ番目の皮膚モデルにおける第ｉ主成分の主成分分析スコアαｋｉを求める。次に、求めた主成分分析スコアαｋｉを（２）式の右辺に代入するとともに、複数の皮膚モデルの各々について設定された光吸収係数のモデル値μａｋｊを（２）式の左辺に代入することにより、検量線ＰＩＬＳｉｊを求める。
【００７２】
なお、上記の（２）式は一般式である。（２）式を、本実施形態における三層構造に適用して行列式に変形すると下記の（５）式となる。
【００７３】
【数５】
2014006063.tif 000008
【００７４】
ここで、光吸収係数のモデル値μａｘ１は、表皮層の光吸収係数のモデル値を示し、光吸収係数のモデル値μａｘ２は、真皮層の光吸収係数のモデル値を示し、光吸収係数のモデル値μａｘ３は、皮下組織層の光吸収係数のモデル値を示す。
【００７５】
図４に戻り、ユーザが血糖値測定装置１００を皮膚にあてがい、測定開始スイッチ（図示せず）の押下等によって血糖値測定装置１００を動作させると、照射部１０２は、皮膚に対して波長λ１の短時間パルス光を照射する（ステップＳ３）。ここで、波長λ１は、シミュレーション部１１１が時間分解波形を算出した複数の波長の中の１つである。
【００７６】
照射部１０２が短時間パルス光を照射すると、受光部１０３は、照射部１０２から照射され、皮膚によって後方散乱した光を受光する（ステップＳ４）。このとき、受光部１０３は、照射開始からの単位時間毎（例えば、１ピコ秒毎）の受光強度を内部メモリに登録しておく。
【００７７】
受光部１０３が受光を完了すると、計測光強度取得部１０４は、受光部１０３の内部メモリに格納されている、異なる時刻ｔにおける受光強度Ｉ’ｍ（ｔ）を取得する（ステップＳ５）。
【００７８】
計測光強度取得部１０４が、受光強度Ｉ’ｍ（ｔ）を取得すると、光吸収係数算出部１０６は、上記の（３）式、（４）式に基づいて、受光強度Ｉ’ｍ（ｔ）を検量線に照合して、皮膚の各層の光吸収係数μａｍｊを算出する（ステップＳ６）。
【００７９】
但し、皮膚から得られる信号の第ｉ主成分の主成分スコアをαｍｉ、受光部１０３が時刻ｔにおいて受光した入射光強度に対する光強度をＩ’ｍ（ｔ）、測定に用いられる入射光強度をＩｍｉｎ、時刻ｔにおいて受光した光強度をＩｍ（ｔ）、皮膚の第ｊ層の光吸収係数をμａｍｊとする。
【００８０】
具体的には、シミュレーション部１１１で算出した、複数の皮膚モデルにおける入射光強度に対する無吸収時光強度の時間関数Ｎ’ｓ（ｔ）、第ｉ主成分の主成分分析ローディングＰｉ（ｔ）、及び受光部１０３が時刻ｔにおいて受光した入射光強度に対する光強度Ｉ’ｍ（ｔ）を（３）式の右辺に代入することにより、実測値としての第ｉ主成分の主成分スコアαｍｉを求める。次に、求めた主成分スコアαｍｉ及び検量線ＰＩＬＳｉｊを（４）式の右辺に代入することにより、実測値としての光吸収係数μａｍｊを求める。
【００８１】
なお、上記の（４）式は一般式である。（４）式を、本実施形態における三層構造に適用して行列式に変形すると下記の（６）式となる。
【００８２】
【数６】
2014006063.tif 000009
【００８３】
ここで、光吸収係数μａｍ１は、実測値としての表皮層の光吸収係数を示し、光吸収係数μａｍ２は、実測値としての真皮層の光吸収係数を示し、光吸収係数μａｍ３は、実測値としての皮下組織層の光吸収係数を示す。
【００８４】
表２は、光吸収係数の真値と、本発明に係る主成分回帰分析により算出した実測値としての光吸収係数と、を比較した表である。
【００８５】
【表２】
2014006063.tif 000010
【００８６】
表２に示すように、表皮層の光吸収係数の真値は０．１１ｍｍ−１、真皮層の光吸収係数の真値は０．５５ｍｍ−１、皮下組織層の光吸収係数は０．１ｍｍ−１である。これに対し、本発明に係る主成分回帰分析により算出した実測値としての表皮層の光吸収係数は０．１０９４ｍｍ−１、真皮層の光吸収係数は０．５５２７ｍｍ−１、皮下組織層の光吸収係数は０．１０５ｍｍ−１である。このようにして算出された実測値としての光吸収係数は、光吸収係数の真値に極めて近い値となることが分かる。
【００８７】
光吸収係数算出部１０６が皮膚の各層の光吸収係数μａｍ１〜μａｍ３を算出すると、光吸収係数算出部１０６は、皮膚の主要物質の種類数と同じ数の波長に対して光吸収係数μａｍ１〜μａｍ３を算出したか否かを判定する（ステップＳ７）。本実施形態では、皮膚の主要物質を水、たんぱく質、脂質、グルコースの４種類として血糖値の測定を行うため、光吸収係数算出部１０６は、４種類の波長λ１〜λ４に対して光吸収係数μａｍ１〜μａｍ３を算出したか否かを判定する。ここで、波長λ１〜λ４は、シミュレーション部１１１が時間分解波形を算出した複数の波長の中から選出する。
【００８８】
光吸収係数算出部１０６が、光吸収係数μａｍ１〜μａｍ３を算出していない波長λ１〜λ４があると判定した場合（ステップＳ７：ＮＯ）、ステップＳ１に戻り、まだ光吸収係数μａｍ１〜μａｍ３を算出していない波長λ１〜λ４の光吸収係数μａｍ１〜μａｍ３の算出を行う。
【００８９】
他方、光吸収係数算出部１０６が、波長λ１〜λ４の光吸収係数μａｍ１〜μａｍ３を算出していると判定した場合（ステップＳ７：ＹＥＳ）、濃度算出部１０８は、式（７）に基づいて真皮質に含まれるグルコースの濃度を算出する（ステップＳ８）。
【００９０】
【数７】
2014006063.tif 000011
【００９１】
但し、第ｍ層における光吸収係数をμａｍ、皮膚を形成する第ｈ物質の光吸収係数をμａｈ、皮膚を形成する第ｈ物質の体積濃度をｃｖｈとする。
なお、上記の（７）式は一般式である。（７）式を、本実施形態における４つの波長に適用するよう変形すると下記の（８）式となる。
【００９２】
【数８】
2014006063.tif 000012
【００９３】
但し、μａ２（λ）は真皮層における波長λ１〜λ４の光吸収係数、μａｗ（λ）は真皮層における波長λ１〜λ４の水の光吸収係数、μａｐ（λ）は真皮層における波長λ１〜λ４のたんぱく質の光吸収係数、μａｌ（λ）は真皮層における波長λ１〜λ４の脂質の光吸収係数、μａｇ（λ）は真皮層における波長λ１〜λ４のグルコースの光吸収係数を示す。また、ｃｖｗは水の体積濃度（体積分率）、ｃｖｐはたんぱく質の体積濃度（体積分率）、ｃｖｌは脂質の体積濃度（体積分率）、ｃｖｇはグルコースの体積濃度（体積分率）を示す。
【００９４】
濃度算出部１０８は、物質吸収情報記憶部１０７に記憶されている測定対象中の主要物質の光吸収係数と、上記の（８）式により算出された真皮層における光吸収係数μａ２とからグルコースの濃度を算出する。
【００９５】
なお、上記の（７）式に替えて下記の（９）式を用いてグルコースの濃度を算出してもよい。
【００９６】
【数９】
2014006063.tif 000013
【００９７】
但し、第ｍ層における光吸収係数をμａｍ、皮膚を形成する第ｈ物質のモル吸光係数をεｈ、皮膚を形成する第ｈ物質のモル濃度をｃｈとする。
なお、上記の（９）式は一般式である。（９）式を、本実施形態における三層構造に適用するよう変形すると下記の（１０）式となる。
【００９８】
【数１０】
2014006063.tif 000014
【００９９】
但し、εｗ（λ）は真皮層における波長λ１〜λ４の水のモル吸光係数、εｐ（λ）は真皮層における波長λ１〜λ４のたんぱく質のモル吸光係数、εｌ（λ）は真皮層における波長λ１〜λ４の脂質のモル吸光係数、εｇ（λ）は真皮層における波長λ１〜λ４のグルコースのモル吸光係数を示す。また、ｃｗは水のモル濃度、ｃｐはたんぱく質のモル濃度、ｃｌは脂質のモル濃度、ｃｇはグルコースのモル濃度を示す。
【０１００】
濃度算出部１０８は、物質吸収情報記憶部１０７に記憶されている測定対象中の主要物質のモル吸光係数と、上記の（１０）式により算出された真皮層におけるモル吸光係数εとからグルコースの濃度を算出する。
【０１０１】
濃度単位変換部１０９は、濃度算出部１０８で算出したグルコースの濃度の単位を所望の単位に変換する。濃度表示部１１０は、グルコースの濃度を表示する。
【０１０２】
このように、本実施形態によれば、複数の皮膚モデルにおける短時間パルス光の時間分解波形の所定の時刻における光強度（光強度の時間関数）と、複数の皮膚モデルにおける短時間パルス光の時間分解波形の所定の時刻における無吸収時の光強度（無吸収時光強度の時間関数）と、複数の皮膚モデルの各々について設定された皮膚の各々の層の光吸収係数のモデル値と、に基づいて作成された検量線を取得し、光強度取得手段が取得した光強度を検量線と照合して、任意の層の光吸収係数を選択的に算出することができる。仮に、光吸収係数を算出する際に、近似式を用いると、得られる結果は近似解となり、特に光吸収係数の大きい層と光吸収係数の小さい層とが隣接している場合に誤差が大きくなる。これに対し、本実施形態においては、光強度取得手段が取得した光強度が検量線と照合されるので、その算出結果は高精度となる。そのため、算出した光吸収係数に基づいて目的物質の濃度を算出することにより、他の層によるノイズの影響を低減し、精度の高い濃度の定量を行うことができる。
【０１０３】
また、本実施形態によれば、上記の（１）式〜（４）式を用いた主成分回帰分析により、光吸収係数が算出される。よって、任意の層の光吸収係数を高精度に算出することができる。
【０１０４】
また、本実施形態によれば、算出した光吸収係数に基づいて真皮層に含まれるグルコースの濃度を算出することにより、他の層によるノイズの影響を低減し、グルコースの濃度の定量を高精度で行うことができる。
【０１０５】
なお、算出した光吸収係数に基づいて真皮層に含まれるグルコースの濃度を算出する方法として、主成分回帰分析や部分最小二乗回帰分析等の多変量解析を用いてもよい。
【０１０６】
以上、図面を参照してこの発明の一実施形態について詳しく説明してきたが、具体的な構成は上述のものに限られることはなく、この発明の要旨を逸脱しない範囲内において様々な設計変更等をすることが可能である。
例えば、本実施形態では、濃度定量方法を血糖値測定装置１００に実装し、皮膚の真皮層に含まれるグルコースの濃度を測定する場合を説明したが、これに限られず、濃度定量方法を、複数の光散乱媒質の層から形成される観測対象の任意の層における目的物質の濃度を定量する他の装置に用いても良い。
【０１０７】
上述の血糖値測定装置１００は内部に、コンピュータシステムを有している。そして、上述した各処理部の動作は、プログラムの形式でコンピュータ読み取り可能な記録媒体に記憶されており、このプログラムをコンピュータが読み出して実行することによって、上記処理が行われる。ここでコンピュータ読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、半導体メモリ等をいう。また、このコンピュータプログラムを通信回線によってコンピュータに配信し、この配信を受けたコンピュータが当該プログラムを実行するようにしても良い。
【０１０８】
また、上記プログラムは、前述した機能の一部を実現するためのものであっても良い。さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であっても良い。
【０１０９】
また、上記実施形態においては、観測対象を層構成で説明したが、格子状（メッシュ構成）で構成される対象にも適用できる。
【符号の説明】
【０１１０】
１００…血糖値測定装置（濃度定量装置）、１０１…光吸収係数算出装置、１０２…照射部（照射手段）、１０３…受光部（受光手段）、１０４…計測光強度取得部（光強度取得手段）、１０６…光吸収係数算出部（光吸収係数算出手段）、１０８…濃度算出部（濃度算出手段）
【図１】
2014006063.tif 000018
【図２】
2014006063.tif 000019
【図３】
2014006063.tif 000020
【図４】
2014006063.tif 000021