(57)【要約】      （修正有）【課題】吸収散乱性物質の分布を精度よく求めることができる光学的測定方法を提供する。【解決手段】光学的測定方法は、吸収散乱性物質を含むサンプルに光を照射または伝搬させたときに生じる反射光と参照光とを干渉させて当該干渉光データを波長の関数として検出し、干渉光データを取得するデータ取得ステップと、複数の波長帯における干渉光データから複数の波長帯における光断層画像の構成要素であるラインデータを作成するラインデータ作成ステップと、複数の波長帯におけるそれぞれのラインデータに対して、光の進行方向において空間分解能よりも狭い範囲の領域で反射光量が小さい向きに凸の溝をノイズとして検出し、当該領域において溝の深さを浅くすることによりノイズを低減するノイズ補正ステップと、ノイズを補正されたラインデータに基づいて吸収散乱性物質の分布を求めるデータ処理ステップと、を備える。【選択図】図６
2015230168.tif 000002
【特許請求の範囲】
【請求項１】
複数の吸収散乱性物質を含むサンプルに光を照射または伝搬させたときに該サンプルにおいて生じる反射光と参照光とを干渉させて当該干渉光データを波長の関数として検出し、分割された複数の波長帯における干渉光データを取得するデータ取得ステップと、
このデータ取得ステップにおいて取得された前記複数の波長帯における干渉光データそれぞれをフーリエ変換し、前記複数の波長帯における光断層画像の構成要素であるラインデータを作成するラインデータ作成ステップと、
このラインデータ作成ステップにおいて作成された前記複数の波長帯におけるそれぞれのラインデータに対して、前記光の進行方向において空間分解能よりも狭い範囲の領域で反射光量が小さい向きに凸の溝をノイズとして検出し、当該領域において前記溝の深さを浅くすることによりノイズを低減するノイズ補正ステップと、
このノイズ補正ステップにおいてノイズを補正されたラインデータに基づいて、前記複数の吸収散乱性物質それぞれについて、入射光量に対する反射光量の減衰率に含まれる吸収散乱成分を求め、この吸収散乱成分に基づいて該吸収散乱性物質の分布を求めるデータ処理ステップと、
を備える、
光学的測定方法。
【請求項２】
前記ノイズ補正ステップが、広い範囲から狭い範囲の順に前記領域の範囲を変えて行われる、
請求項１記載の光学的測定方法。
【請求項３】
前記ノイズ補正ステップが、
前記領域に隣接する前記領域外のデータそれぞれを頂点に含む反射光量が小さい向きに凸の多角形の辺または頂点上に、前記領域内のデータを補正する、請求項１または２記載の光学的測定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、光学的測定方法に関するものである。
【背景技術】
【０００２】
ＯＣＴ（Optical  Coherence  Tomography），ＯＴＤＲ（Optical  Time  Domain  Reflectometer）およびＤＯＴ（Diffuse  Optical  Tomography）等の光学的測定方法では、サンプルに光を照射または伝搬させたときに該サンプルにおいて生じる反射光を検出して、その検出結果に基づいてサンプルの情報を得ることができる。ＯＣＴでは、サンプルにおいて生じる反射光と参照光とを干渉させた光を検出することでサンプルの深さ方向の一次元の光断層画像を取得することができ、また、サンプルへの光照射の位置を走査することで二次元または三次元の光断層画像を取得することができる。
【０００３】
このような光学的測定方法において、サンプルに入射した光量に対してサンプルで生じた反射光の光量は吸収または散乱による減衰が含まれる場合がある。非特許文献１に記載された光学的測定方法は、複数の吸収散乱性物質を含むサンプルに光を照射または伝搬させたときに、入射光量に対する反射光量の減衰率に含まれる吸収散乱成分を複数の吸収散乱性物質それぞれについて互いに分離して求め、サンプル内部の複数の吸収散乱性物質それぞれの分布を求めることができる。
【０００４】
この光学的測定方法は、具体的には、複数の吸収散乱性物質を含むサンプルに光を照射または伝搬させたときに該サンプルにおいて生じる反射光と参照光とを干渉させて当該干渉光データを波長の関数として検出し、分割された複数の波長帯における干渉光データを取得する。そして、この取得した複数の波長帯における干渉光データそれぞれをフーリエ変換し、複数の波長帯における光断層画像を作成する。続いて、この光断層画像に基づいて、複数の吸収散乱性物質それぞれについて、入射光量に対する反射光量の減衰率に含まれる吸収散乱成分を求め、この吸収散乱成分に基づいて該吸収散乱性物質の分布を求める。
【０００５】
ところで、この光学的測定方法では、複数の波長帯における光断層画像それぞれについて、異常値の位置に対応する値がノイズとしてノイズ低減処理される。ここで、異常値とは、入射光量に対する反射光量の減衰率とそのモデル関数との残差平方和の値が閾値を超えるものである。
【先行技術文献】
【非特許文献】
【０００６】
【非特許文献１】Chenyang  Xu,et  al.,"Separation  of  absorption  and  scattering  profiles  in  spectroscopic  optical  coherence  tomography  using  a  least-squares  algorithm",Optics  Express,Vol.12,No.20,pp.4790-4803(2004).
【発明の概要】
【発明が解決しようとする課題】
【０００７】
本発明者は、上記のような光学的測定方法におけるノイズ処理によっては、求められる吸収散乱性物質の分布の精度が向上しない場合があることを見出した。本発明は、上記問題を解消する為になされたものであり、吸収散乱性物質の分布を精度よく求めることができる光学的測定方法を提供することを目的とする。
【課題を解決するための手段】
【０００８】
本発明の光学的測定方法は、複数の吸収散乱性物質を含むサンプルに光を照射または伝搬させたときに該サンプルにおいて生じる反射光と参照光とを干渉させて当該干渉光データを波長の関数として検出し、分割された複数の波長帯における干渉光データを取得するデータ取得ステップと、このデータ取得ステップにおいて取得された複数の波長帯における干渉光データそれぞれをフーリエ変換し、複数の波長帯における光断層画像の構成要素であるラインデータを作成するラインデータ作成ステップと、このラインデータ作成ステップにおいて作成された複数の波長帯におけるそれぞれのラインデータに対して、光の進行方向において空間分解能よりも狭い範囲の領域で反射光量が小さい向きに凸の溝をノイズとして検出し、当該領域において溝の深さを浅くすることによりノイズを低減するノイズ補正ステップと、このノイズ補正ステップにおいてノイズを補正されたラインデータに基づいて、複数の吸収散乱性物質それぞれについて、入射光量に対する反射光量の減衰率に含まれる吸収散乱成分を求め、この吸収散乱成分に基づいて該吸収散乱性物質の分布を求めるデータ処理ステップと、を備える。
【０００９】
本発明の光学的測定方法では、ノイズ補正ステップが、広い範囲から狭い範囲の順に領域の範囲を変えて行われてもよい。また、ノイズ補正ステップが、領域に隣接する領域外のデータそれぞれを頂点に含む反射光量が小さい向きに凸の多角形の辺または頂点上に、領域内のデータを補正してもよい。
【発明の効果】
【００１０】
本発明によれば、吸収散乱性物質の分布を精度よく求めることができる。
【図面の簡単な説明】
【００１１】
【図１】本実施形態の光学的測定方法が適用される光学的測定システム１の概略構成を示す図である。
【図２】本実施形態の光学的測定方法におけるデータ取得ステップおよびラインデータ作成ステップを説明する図である。
【図３】本実施形態の光学的測定方法におけるデータ取得ステップおよびラインデータ作成ステップを説明するグラフである。
【図４】光学的測定システム１により取得される干渉光データ等の例を示す図である。
【図５】ある波長帯域における光断層画像のラインデータの一例を示すグラフである。
【図６】本実施形態の光学的測定方法におけるノイズ補正ステップを説明するグラフである。
【図７】図５のラインデータに対してノイズ補正ステップを行った後のラインデータの例を示すグラフである。
【図８】図７のラインデータに対してノイズ補正ステップを行った後のラインデータの例を示すグラフである。
【図９】行列Ｘの一例を模式的に示す図である。
【図１０】行列Ｌの一例を模式的に示す図である。
【図１１】行列Ｙの一例を模式的に示す図である。
【図１２】行列Ａの一例を模式的に示す図である。
【図１３】本実施形態の光学的測定方法により求められる吸収散乱性物質の濃度積分値の例を示すグラフである。
【図１４】比較例の光学的測定方法により求められる吸収散乱性物質の濃度積分値の例を示すグラフである。
【発明を実施するための形態】
【００１２】
以下、添付図面を参照して、本発明を実施するための形態を詳細に説明する。なお、図面の説明において同一の要素には同一の符号を付し、重複する説明を省略する。本発明は、これらの例示に限定されるものではなく、特許請求の範囲によって示され、特許請求の範囲と均等の意味および範囲内でのすべての変更が含まれることが意図される。
【００１３】
図１は、本実施形態の光学的測定方法が適用される光学的測定システム１の概略構成を示す図である。光学的測定システム１は、ＯＣＴ装置２およびプローブ３を備え、サンプル９の光断層画像を取得する。ＯＣＴ装置２は、光源から出力された光を２分岐して観察光および参照光として出力する。プローブ３は、ＯＣＴ装置２の光分岐部から出力された観察光を導光して先端からサンプル９へ照射するとともに、その光照射によって該サンプル９で生じる反射光を先端に入力してＯＣＴ装置２へ導光する。ＯＣＴ装置２は、プローブ３により導光されて到達した反射光と参照光とを干渉させて当該干渉光を検出する。そして、ＯＣＴ装置２は、この干渉光に基づいて、サンプル９の深さ方向の光断層画像を取得することができる。また、プローブ３によるサンプル９への光照射の位置を走査（スキャン）することで、ＯＣＴ装置２は二次元または三次元の光断層画像を取得することができる。
【００１４】
本実施形態の光学的測定方法は、データ取得ステップ、ラインデータ作成ステップ、ノイズ補正ステップおよびデータ処理ステップを備える。データ取得ステップでは、複数の吸収散乱性物質を含むサンプル９に光を照射または伝搬させたときに該サンプル９において生じる反射光と参照光とを干渉させて当該干渉光データを波長の関数として検出し、分割された複数の波長帯における干渉光データを取得する。ラインデータ作成ステップでは、データ取得ステップにおいて取得された複数の帯域における干渉光データそれぞれをフーリエ変換し、複数の波長帯におけるラインデータを作成する。ノイズ補正ステップでは、ラインデータ作成ステップにおいて作成された複数の波長帯におけるそれぞれのラインデータに対して、光の進行方向において空間分解能よりも狭い範囲の領域で反射光量が小さい向きに凸の溝をノイズとして検出し、当該領域において溝の深さを浅くすることによりノイズを低減する。データ処理ステップでは、ノイズ補正ステップにおいてノイズを補正されたラインデータに基づいて、複数の吸収散乱性物質それぞれについて、入射光量に対する反射光量の減衰率に含まれる吸収散乱成分を求め、この吸収散乱成分に基づいて該吸収散乱性物質の分布を求める。
【００１５】
光学的測定システム１は、図２および図３を用いて説明するように、本実施形態の光学的測定方法による測定を行うこともできる。図２および図３は、本実施形態の光学的測定方法におけるデータ取得ステップおよびラインデータ作成ステップを説明する図である。干渉光のスペクトル（図２（ａ））に対して、帯域を分割する為の複数の窓関数（図２（ｂ））を乗じ、各帯域のスペクトル（図３（ａ））をフーリエ変換する。窓関数は例えばハン関数である。このフーリエ変換の結果、複数の波長帯域におけるラインデータ（図３（ｂ））を作成することができる。ラインデータは、横軸が光の進行方向、即ちサンプル９における深さ方向の位置ｚであり、サンプル９表面でｚ＝０である。また、縦軸が波数ｋにおけるサンプル９内部の位置ｚでのログスケールでの入射光量に対する反射光量の減衰率（反射率（または輝度））Ｙ(ｋ,ｚ)であり（詳しくは後述）、奥行き軸が波数ｋである。奥行き軸の波数ｋは各波数帯域の代表値であり、例えば中心値である。そして、この光断層画像に基づいて、サンプル９内部の吸収性または散乱性を有する吸収散乱性物質（例えば脂質）の分布を求めることができる。
【００１６】
図４は、光学的測定システム１により取得される干渉光データ等の例を示す図である。データ取得ステップにおいて、干渉光を検出することで、波長λに対する干渉光強度（信号強度）を得ることができる。この例では、波長帯域を６分割して波長帯域毎にフーリエ変換することで、波長帯域毎の光断層画像を取得することができる。
【００１７】
図５は、ある波長帯域における光断層画像のラインデータの一例を示すグラフである。すなわち、図５に示されるラインデータは、図３（ｂ）に示される複数の波長帯域における光断層画像の構成要素であるラインデータの一つに相当する。ここで、光の進行方向、即ち、サンプル９の深さ方向の位置ｚにおける光断層画像の空間分解能は、データ取得ステップにおける干渉光データの取得条件（例えば、使用する光源の波長帯域および波長帯域の分割数等）から一義的に決まる。したがって、光の進行方向において当該空間分解能よりも狭い範囲の領域で値が振動するデータについては、高周波ノイズの影響を受けているとみなすことができる。特にＯＣＴでは、光の干渉を利用しているので、取得される光断層画像にスペックルが生じ易いという特性があり、ラインデータにおいて当該スペックルは通常のシグナルよりも急峻な反射光量が小さい向きに凸の溝として現れる。そこで、ノイズ補正ステップでは、以下のようにして光の進行方向において空間分解能よりも狭い範囲の領域で反射光量が小さい向きに凸の溝をノイズとして検出し、当該領域において溝の深さを浅くすることによりノイズを低減する。
【００１８】
図６は、本実施形態の光学的測定方法におけるノイズ補正ステップを説明するグラフである。図６に示される補正前のラインデータは、図５のグラフを一部拡大したものである。ノイズ補正ステップは、領域ａに隣接する領域ａ外のデータそれぞれを頂点に含む反射光量が小さい向きに凸の多角形の辺または頂点上に、領域ａ内のデータを補正する。
【００１９】
一例として、ここでは空間分解能を２０ピクセルとし、これより狭い範囲の領域ａとして１９ピクセルの領域で反射光量が小さい向きに凸の溝を検出する。領域ａは奇数ピクセルの範囲とされる。このとき溝の中央の最下部（極小点）に位置するデータが領域ａの中央のピクセルに対応し、その両側に等しく９ピクセルずつが配置されるようにする。図６中の各データについて、領域ａの中央のデータをデータＤ０とし、データＤ０から深さ方向の位置ｚが大きくなる方へと順番にデータＤ１、データＤ２、データＤ３・・・データＤＮ（Ｎは自然数）とし、データＤ０から深さ方向の位置ｚが小さくなる方へと順番にデータＤ−１、データＤ−２、データＤ−３・・・データＤ−Ｎとする。領域ａには、このうちデータＤ−９からデータＤ９までが含まれる。
【００２０】
溝の検出は、データＤ０の反射率＜データＤ５の反射率＜データＤ９の反射率、かつデータＤ０の反射率＜データＤ−５の反射率＜データＤ−９の反射率という条件で５つのデータを用いて行われる。領域ａ内の一部のデータのみを用いることにより、検出速度が向上するとともに、検出精度が高くなる。なお、領域ａ内の全てのデータを用いて、データＤ０の反射率＜データＤ１の反射率＜データＤ２の反射率＜・・・＜データＤ９の反射率、かつデータＤ０の反射率＜データＤ−１の反射率＜データＤ−２の反射率＜・・・＜データＤ−９の反射率という条件で行ってもよい。
【００２１】
溝の検出に続いて、第１ステップとして、データＤ−９の反射率とデータＤ９の反射率との平均値を求め、当該平均値をデータＤ０の補正後のデータＤ０'の反射率に設定する。次に、第２ステップとして、データＤ−９の反射率を第１ステップで補正されたデータＤ０'の反射率とデータＤ−１０の反射率との間の値に補正するとともに、データＤ９の反射率を上記データＤ０'の反射率とデータＤ１０の反射率との間の値に補正する。
【００２２】
続いて、第３ステップとして、領域ａに隣接する領域ａ外のデータＤ−１０およびデータＤ１０と、第１ステップで補正されたデータＤ０'と、第２ステップで補正されたデータＤ−９'およびデータＤ９'とを頂点とする五角形を作成する。当該五角形は、データＤ−１０とデータＤ−９'との間の辺ｂ１、データＤ−９'とデータＤ０'との間の辺ｂ２、データＤ０'とデータＤ９'との間の辺ｂ３、およびデータＤ９'とデータＤ１０との間の辺ｂ４の四辺を含む。そして、辺ｂ２上にデータＤ−８〜データＤ−１の補正後のデータＤ−８'〜データＤ−１'を補正するとともに、辺ｂ３上にデータＤ１〜データＤ８の補正後のデータＤ１'〜データＤ８'を補正する。同様にして、検出される全ての溝に対し当該ノイズ補正処理を行う。
【００２３】
図７は、図５のラインデータに対してノイズ補正ステップを行った後のラインデータの例を示すグラフである。図７に示されるように、ノイズ補正ステップによって図５のラインデータに存在する高周波ノイズが低減される。続いて、図７に示されるラインデータについて、領域ａの範囲を例えば１１ピクセルに変えて同様に反射光量が小さい向きに凸の溝をノイズとして検出し、上記第１ステップ〜上記第３ステップを行うノイズ補正ステップを繰り返す。このとき、反射光量が小さい向きに凸の溝の検出は、データＤ０の反射率＜データＤ３の反射率＜データＤ５の反射率、かつデータＤ０の反射率＜データＤ−３の反射率＜データＤ−５の反射率という条件で５つのデータを用いて行われる。もちろん全てのデータを用いて行うこともできる。
【００２４】
図８は、図７のラインデータに対してノイズ補正ステップを行った後のラインデータの例を示すグラフである。図８に示されるように、２回のノイズ補正ステップを行うことにより、図５のラインデータに存在した高周波ノイズのうち、１回のノイズ補正ステップで低減されなかったより高周波のノイズについても低減されることとなり、全体として高周波ノイズが一層低減される。図８に示されるラインデータについて、さらに領域ａの範囲を例えば７ピクセルや５ピクセルに変えて同様にノイズ補正ステップを繰り返してもよい。なお、このノイズ補正の手法によれば、領域ａの範囲の最小値は５ピクセルである。
【００２５】
このようなノイズ補正ステップが、広い範囲から狭い範囲の順に領域ａの範囲を変えて、すべての波長帯域における光断層画像のすべてのラインデータについて行われる。続くデータ処理ステップでは、このノイズを補正された光断層画像に基づいて以下のような処理を行う。
【００２６】
干渉光強度Ｉ(ｋ,ｚ)は下記（１）式で表される。上述のようにｋは波数である。また、ｚは、サンプル９における深さ方向の位置であり、サンプル９表面でｚ＝０である。ｍは、サンプル９に含まれる吸収散乱性物質の種類の数である。Ｓ(ｋ)は、光源のスペクトルである。Ｈｓ(ｋ)Ｈｓ(ｚ)は、光学的測定システム１の応答関数である。Ｈｒ(ｋ)Ｈｒ(ｚ)は、サンプル９における後方散乱の応答関数である。εｉ(ｋ)は、サンプル９に含まれるｍ種類の吸収散乱性物質のうちの第ｉの吸収散乱性物質の吸収散乱係数である。ｆｉ(ｚ)は、サンプル９に含まれる第ｉの吸収散乱性物質の空間的濃度分布である。
【００２７】
【数１】
2015230168.tif 000003
【００２８】
このように、干渉光強度Ｉ(ｋ,ｚ)は、光源スペクトルＳ(ｋ)、光学的測定システム１の応答関数Ｈｓ(ｋ)Ｈｓ(ｚ)、サンプル９内部における吸収または散乱による減衰を表す因子（数式中の指数関数の部分）、および、サンプル９における後方散乱の応答関数Ｈｒ(ｋ)Ｈｒ(ｚ)、の積で表される。
【００２９】
波数ｋにおけるサンプル９内部の位置ｚでのログスケールでの入射光量に対する反射光量の減衰率Ｙ(ｋ,ｚ)は下記（２）式で表される。（２）式中のＦｉ(ｚ)は、下記（３）式で表され、サンプル９に含まれる第ｉの吸収散乱性物質による吸収散乱成分であり、第ｉの吸収散乱性物質の濃度のｚ＝０からの積分値である。（２）式中のＣ(ｚ)は、下記（４）式で表され、吸収・散乱が波長に依存しない波長無依存物質による吸収散乱成分であり、この波長無依存物質の濃度のｚ＝０からの積分値である。この波長無依存物質は第(ｍ＋１)の吸収散乱性物質として扱うことができる。
【００３０】
【数２】
2015230168.tif 000004
【００３１】
【数３】
2015230168.tif 000005
【００３２】
【数４】
2015230168.tif 000006
【００３３】
データ取得ステップにおける波長帯の分割数をｎとする。任意の深さｚにおいて、最少二乗法を用いて、下記（５）式で表される目的関数Ｓが最小（目的関数の偏微分がゼロ）となるような解Ｘを求めると、下記（６）式〜（９）式になる。ｋｊは、ｎ個の波長帯のうちの第ｊ波長帯の波数である。（ｍ＋１）行(ｍ＋１)列の行列Ａは、既知の行列要素からなる。(ｍ＋１)行１列の行列Ｘは、データ処理ステップにおいて求められるべき未知の行列要素からなる。（ｍ＋１）行１列の行列Ｙは、データ取得ステップにおいて得られた減衰率と既知の要素からなる。
【００３４】
【数５】
2015230168.tif 000007
【００３５】
【数６】
2015230168.tif 000008
【００３６】
【数７】
2015230168.tif 000009
【００３７】
【数８】
2015230168.tif 000010
【００３８】
【数９】
2015230168.tif 000011
【００３９】
しかし、実際に計算して解Ｘを求めてみると、光断層像の空間的なノイズが重畳された解Ｘが得られる。
【００４０】
このノイズ問題に対処する為、各Ｆｉ(ｚ)がノイズに比べて空間的に滑らかに変化する（空間的な変化が小さい）ことを仮定して、目的関数を修正する。すなわち、各Ｆｉ(ｚ)の二次微分値が小さいと仮定する。修正後の目的関数Ｓ'は下記（１０）式で表される。αは、Ｓに対する二次微分値の重みを決めるパラメータである。α２は、各Ｆｉ(ｚ)の二次微分値を小さくするためのパラメータであるので、平滑化パラメータと呼ぶことにする。目的関数Ｓ'は、減衰率Ｙ，吸収散乱成分Ｆｉおよび吸収散乱成分Ｆｉの空間二次微分を含む関数である。
【００４１】
【数１０】
2015230168.tif 000012
【００４２】
解Ｘは、行列を用いて表現すると下記（１１）式のように表される。Ｌは二次微分演算子（ラプラシアン演算子）である。
【００４３】
【数１１】
2015230168.tif 000013
【００４４】
Ｆｉ(ｚ)は離散的なデータであるので、或る位置での二次微分値は、該位置のデータに加えて前後のデータも用いた計算により求める。計算により求められる二次微分値は、計算手法によって異なり、また、計算に用いるデータの点数によっても異なる。計算手法としては、差分法、Savitzky-Golay  filterを用いた数値微分、前進差分近似、後退差分近似、中心差分近似、Stirlingの公式、などが知られている。
【００４５】
図９は、行列Ｘの一例を模式的に示す図である。図１０は、行列Ｌの一例を模式的に示す図である。図１１は、行列Ｙの一例を模式的に示す図である。図１２は、行列Ａの一例を模式的に示す図である。ここで、深さ方向のデータ数は５１２である。波長帯の分割数はｎである。波長無依存物質を除く吸収散乱性物質の種類の数ｍは１である。Savitzky-Golay  filterを用いた数値微分により３点のデータを用いて二次微分を計算する場合の行列Ｌを示している。二次微分演算子である行列Ｌは、Ｆｉ(ｚ)に作用するが、Ｃ(ｚ)には作用しない。
【００４６】
解Ｘが求められると、すなわち、各吸収散乱性物質についてＦｉ(ｚ)，Ｃ(ｚ)が求められると、これをｚで微分することで各吸収散乱性物質の空間的分布ｆｉ(ｚ)，ｃ(ｚ)が得られる。第ｉの吸収散乱性物質の分布ｆｉ(ｚ)は下記（１２）式で表される。第(ｍ＋１)の吸収散乱性物質（波長無依存物質）の分布ｃ(ｚ)は下記（１３）式で表される。なお、ノイズ低減の為に、各吸収散乱性物質の分布ｆｉ(ｚ)，ｃ(ｚ)について空間上または時間上の平均化処理を行うのが好適である。
【００４７】
【数１２】
2015230168.tif 000014
【００４８】
【数１３】
2015230168.tif 000015
【００４９】
以上説明したように、本実施形態の光学的測定方法は、ラインデータ作成ステップにおいて作成された複数の波長帯における光断層画像それぞれのラインデータに対して、光の進行方向、即ちサンプル９における深さ方向において空間分解能よりも狭い範囲の領域ａで反射光量が小さい向きに凸の溝をノイズとして検出し、当該領域ａにおいて溝の深さを浅くすることによりノイズを低減するノイズ補正ステップを備える。
【００５０】
これによれば、ラインデータにおいて、空間分解能よりも狭い範囲で値が振動するデータについては、高周波ノイズの影響を受けているとみなし、複数の波長帯におけるラインデータすべてに対して各々でノイズ低減処理を直接行うので、分光情報としてのラインデータの精度を向上させることができる。したがって、補正されたラインデータに基づいて求められる吸収散乱成分Ｆｉ（ｚ）の精度を向上させることができるとともに、吸収散乱成分Ｆｉ（ｚ）に基づいて求められる吸収散乱性物質の分布ｆｉ（ｚ）の精度を向上させることができる。
【００５１】
これに対して、比較例としての非特許文献１の光学的測定方法では、分光ＯＣＴ信号としてのＹ（ｋ，ｚ）と、これを吸収散乱成分Ｆｉ（ｚ），Ｃ（ｚ）および吸収散乱係数εｉ(ｋ)を用いて近似したモデル関数との残差平方和を求め、この残差平方和が閾値を超えて異常値となる位置ｚについてのＹ（ｋ，ｚ）を補正する。つまり、当該光学的測定方法では、複数の波長帯におけるラインデータの全てに対して共通の位置ｚについてのＹ（ｋ，ｚ）を補正するにとどまり、波長帯ごとのノイズを直接検出し補正することができない。また、残差平方和はノイズ部分で局所的に大きくなるものの、全空間領域における変動幅もこのピークと同程度の大きさがあるため、閾値でピークを検出することが困難である。さらに、ピークを検出できた場合も、ノイズで残差平方和が大きくなっているのか、モデル関数が合わないのかが区別できない。加えて、当該光学的測定方法では、このようにして検出したノイズに対して、具体的にどのように補正を行うかの詳細が開示されていない。
【００５２】
図１３は、本実施形態の光学的測定方法により求められる吸収散乱性物質の濃度積分値の例を示すグラフであり、図１４は、比較例の光学的測定方法により求められる吸収散乱性物質の濃度積分値の例を示すグラフである。横軸はサンプル９における深さ方向の位置ｚであり、サンプル９表面でｚ＝０である。また、縦軸は、サンプル９に含まれる第ｉの吸収散乱性物質の分布ｆｉ（ｚ）のｚ＝０からの積分値（濃度積分値）であり、第ｉの吸収散乱性物質による吸収散乱成分Ｆｉ（ｚ）である。図１３によれば、深さ方向の位置ｚが約２５０（ピクセル）付近においてグラフの傾きが変化しており、この付近にＦｉ(ｚ)の微分で示される第ｉの吸収散乱性物質の分布ｆｉ(ｚ)に変化があることがわかり易い。一方、図１４によれば、グラフの全体に高周波ノイズが残り、深さ方向の位置ｚが約２５０（ピクセル）付近においてグラフの傾きが変化していることがわかり難くなっている。なお、深さ方向の位置ｚが約７０（ピクセル）付近においてグラフの傾きが変化しているのは、サンプル９とプローブ３との境界の反射による影響である。
【００５３】
このように、本実施形態の光学的測定方法によれば、非特許文献１の光学的測定方法に比べて高周波ノイズが少ない吸収散乱成分Ｆｉ(ｚ)のグラフが得られる。この結果、本実施形態の光学的測定方法によれば、非特許文献１の光学的測定方法に比べてＦｉ(ｚ)の微分で示される吸収散乱性物質の分布ｆｉ(ｚ)を精度よく求めることができる。
【００５４】
また、本実施形態の光学的測定方法では、ノイズ補正ステップが、広い範囲から狭い範囲の順に領域ａの範囲を変えて繰り返し行われる。このように領域ａの範囲を変えてノイズ補正ステップを行うことで、異なる周波数成分のノイズを除去することができる。また、広い範囲から狭い範囲の順に領域ａの範囲を変えてノイズ補正ステップを行うので、広い範囲のノイズを確実に低減することができる。狭い範囲から広い範囲の順に領域ａの範囲を変えてノイズ補正ステップを行った場合は、例えば、狭い範囲の領域ａで反射光量が小さい向きに凸の溝として検出されノイズ低減が行われたために、仮にその溝を含む部分が広い範囲の領域ａでも反射光量が小さい向きに凸の溝として検出されるべき部分であったとしても、検出されなくなる恐れがある。その結果、補正される領域ａの範囲が狭くなり、高周波ノイズの除去が不完全となる恐れがある。
【００５５】
さらに、この光学的測定方法では、ノイズ補正ステップが、領域ａにおいて光の進行方向の中央のデータＤ０の反射率を両端のデータＤ−９の反射率およびデータＤ９の反射率の平均値に補正する第１ステップと、データＤ−９の反射率を第１ステップで補正されたデータＤ０'の反射率とデータＤ−１０の反射率との間の値に補正するとともに、データＤ９の反射率をデータＤ０'の反射率とデータＤ１０の反射率との間の値に補正する第２ステップと、領域ａに隣接する領域ａ外のデータＤ−１０およびデータＤ１０と、第１ステップで補正されたデータＤ０'と、第２ステップで補正されたデータＤ−９'およびデータＤ９'とを頂点とする反射光量が小さい向きに凸の五角形における辺ｂ２上にデータＤ−８〜データＤ−１の補正後のデータＤ−８'〜データＤ−１'を補正するとともに、辺ｂ３上にデータＤ１〜データＤ８の補正後のデータＤ１'〜データＤ８'を補正する第３ステップと、を含む。
【００５６】
これによれば、領域ａにおいて溝の深さを浅くするとともに、補正対象となる領域ａ内のデータと、補正対象外となる領域ａ外のデータとが滑らかにつながるような補正を行うことができる。なお、本実施形態では、領域ａに隣接する領域ａ外のデータＤ−１０およびデータＤ１０に加え、第１ステップで補正されたデータＤ０'と、第２ステップで補正されたデータＤ−９'およびデータＤ９'とを補正に用いる反射光量が小さい向きに凸の多角形の頂点としたが、これに限られない。データＤ−１０とデータＤ１０とを当該頂点として含む反射光量が小さい向きに凸の多角形を補正に用いるのであれば何でもよく、例えば、データＤ−１０とデータＤ１０とデータＤ０'とを頂点とする三角形を補正に用いてもよいし、データＤ−１０とデータＤ１０とデータＤ−５'とデータＤ５'とを頂点とする四角形を補正に用いてもよいし、六角形以上の多角形を補正に用いてもよい。
【符号の説明】
【００５７】
１…光学的測定システム、２…ＯＣＴ装置、３…プローブ、９…サンプル、ａ…領域、ｂ１〜ｂ４…辺、ＤＮ…データ。
【図１】
2015230168.tif 000016
【図２】
2015230168.tif 000017
【図３】
2015230168.tif 000018
【図５】
2015230168.tif 000019
【図６】
2015230168.tif 000020
【図７】
2015230168.tif 000021
【図８】
2015230168.tif 000022
【図９】
2015230168.tif 000023
【図１０】
2015230168.tif 000024
【図１１】
2015230168.tif 000025
【図１２】
2015230168.tif 000026
【図１３】
2015230168.tif 000027
【図１４】
2015230168.tif 000028
【図４】
2015230168.tif 000029