(57)【特許請求の範囲】
【請求項１】
それぞれが６００ｎｍ乃至１０００ｎｍの波長領域のうちの異なる波長を有する複数の種類の光を被検体に照射することによって、前記被検体の内部の吸収特性と散乱特性のうち少なくとも一方を測定する測定部と、
コラーゲンと脂肪を含む複数の構成要素の全体に対するコラーゲンと脂肪のそれぞれの割合を、各構成要素の吸収係数と前記測定部の測定結果からフィッティング手法を用いて算出し、脂肪及びコラーゲンのそれぞれのフィッティング係数と生体組織の状態の関係と、算出されたコラーゲンの割合と脂肪の割合から前記被検体の生体組織の状態を判定する演算処理部と、
前記演算処理部による処理の結果を表示する表示部と、
を有することを特徴とする測定装置。
【請求項２】
前記測定部は、前記複数の種類の光に、６００ｎｍ乃至７００ｎｍの波長領域のうちの所定の波長を有する光と、７３０ｎｍ乃至７６０ｎｍの波長領域のうちの異なる波長を有する少なくとも２種類の光とを使用することを特徴とする請求項１記載の測定装置。
【請求項３】
前記複数の構成要素は水を更に含み、
前記測定部は、前記複数の種類の光に、９００ｎｍ乃至１０００ｎｍの波長領域のうちの異なる波長を有する２種類の光を使用し、
前記演算処理部は、脂肪と水のそれぞれの量、前記複数の構成要素の全体に対する割合、並びに、脂肪と水の合計の量の分布のうちの少なくとも一つを算出することを特徴とする請求項１又は２記載の測定装置。
【請求項４】
前記演算処理部は、コラーゲンと脂肪のそれぞれの量、前記複数の構成要素の全体に対する割合、並びに、コラーゲンと脂肪の合計の量の分布のうちの少なくとも一つを算出することを特徴とする請求項１乃至３のうちいずれか一項記載の測定装置。
【請求項５】
前記複数の構成要素は、還元ヘモグロビンと酸化ヘモグロビンを更に含み、
前記測定部は、前記複数の種類の光に、７６０ｎｍ乃至８５０ｎｍの波長領域のうちの異なる波長を有する少なくとも２種類の光を使用し、
前記演算処理部は、還元ヘモグロビンと酸化ヘモグロビンのそれぞれの量と、前記複数の構成要素の全体に対する割合、並びに、還元ヘモグロビンと酸化ヘモグロビンの合計の量の少なくとも一つを前記測定部の測定結果から算出することを特徴とする請求項１乃至４のうちいずれか一項記載の測定装置。
【請求項６】
前記演算処理部は、還元ヘモグロビンのモル濃度と酸化ヘモグロビンのモル濃度の合計に対する酸化ヘモグロビンのモル濃度の割合で与えられるヘモグロビンの酸素飽和度を場所ごとに算出することを特徴とする請求項５記載の測定装置。
【請求項７】
前記表示部は、前記状態の分布と前記ヘモグロビンの酸素飽和度の分布を重ねて表示する請求項６記載の測定装置。
【請求項８】
前記測定部は、前記被検体から射出される拡散光を検出する光検出器を有し、
前記演算処理部は、前記被検体の吸収係数を算出し、前記被検体の吸収係数から前記被検体の各構成要素の割合を算出することを特徴とする請求項１乃至７のうちいずれか一項に記載の測定装置。
【請求項９】
前記測定部は、
超音波を発生して被検部位に集束する超音波トランスデューサアレイと、
被検部位における入射光と集束された超音波の相互作用によって発生する信号を検出する光検出器と、
を有し、
前記演算処理部は、前記信号から前記被検体の吸収係数を算出し、前記被検体の吸収係数から前記被検体の各構成要素の割合を算出することを特徴とする請求項１乃至７のうちいずれか一項に記載の測定装置。
【請求項１０】
前記測定部は、前記被検体の内部の被検部位からの弾性波を検出する超音波トランスデューサアレイを有し、
前記演算処理部は、前記被検部位の位置と吸収係数を算出し、前記被検体の吸収係数から前記被検体の各構成要素の割合を算出することを特徴とする請求項１乃至７のうちいずれか一項に記載の測定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体内部の分光特性を測定する測定装置に関する。
【背景技術】
【０００２】
光を用いて測定するマンモグラフィなど、被検体内部の代謝と関係する分光特性を測定して分光特性の分布を画像化する測定装置が知られている。ここで、分光特性は吸収（分光）特性と散乱（分光）特性を含む。医療診断からは、分光特性から生体組織の状態を簡単に判定する技術を確立することが好ましい。ここでいう「生体組織の状態」とは、正常組織、良性腫瘍、悪性腫瘍などを意味する。従来、腫瘍の特定やその種別を行う方法が幾つか提案されている。
【０００３】
特許文献１及び非特許文献１は、近赤外光を利用し、還元ヘモグロビン、酸化ヘモグロビン及びその合計の量に対する酸化ヘモグロビンの割合（酸素飽和度）を求め、腫瘍を特定する方法を開示している。この方法（以下、「ヘモグロビン法」と呼ぶ）は、腫瘍が正常組織よりも還元ヘモグロビンと酸化ヘモグロビンの合計の量は多く、酸素飽和度は少ないことを利用している。
【０００４】
非特許文献２は、被検体から微小組織を抽出し、ラマン分光法を用いてコラーゲンを検出し、検出されたコラーゲンを利用して腫瘍を検出する方法を開示している。正常組織に対する腫瘍のコラーゲンの割合はヘモグロビンの酸素飽和度よりも高いため、コラーゲンを利用する方法は、腫瘍の検出精度が高いという特長を有する。更に、非特許文献２は、脂肪とコラーゲンの全体に対する割合を調べることで腫瘍の状態の判断ができると報告している。
【０００５】
非特許文献３では、幅広い波長領域の近赤外光を利用して被検体の吸収係数を求め、コラーゲンを含む構成要素の割合をフィッティングにより求めている。
【０００６】
非特許文献４では、６３７ｎｍ、６８０ｎｍ、７８５ｎｍ、９０５ｎｍ、９３３ｎｍ、１０６０ｎｍの波長を用いて被検体の構成要素の割合をフィッティングにより求めている。非特許文献４は、コラーゲンをパラメータとして使用することによって精度良くヘモグロビンを見積もり、腫瘍の検出精度を高めている。更に、非特許文献４によれば、コラーゲン量が多いとされる高いマンモグラフィ濃度の領域は腫瘍になる可能性が高い。また、間質の主成分がコラーゲンで、間質の構造と構成の変化が良性・悪性両方の病変の状況と関連していることから、コラーゲンは乳癌発癌の初期の役割を担っている。
【特許文献１】特許第３１０７９１４号明細書
【非特許文献１】Ｎｅａｒ−Ｉｎｆｒａｒｅｄ  Ｃｈａｒａｃｔｅｒｉｚａｔｉｏｎ  ｏｆ  Ｂｒｅａｓｔ  Ｔｕｍｏｒｓ  Ｉｎ  Ｖｉｖｏ  ｕｓｉｎｇ  Ｓｐｅｃｔｒａｌｌｙ−Ｃｏｎｓｔｒａｉｎｅｄ  Ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ，  Ｄａｒｔｍｏｕｔｈ，  Ｐｏｇｕｅ，  ２００５
【非特許文献２】Ｄｉａｇｎｏｓｉｎｇ  ｂｒｅａｓｔ  ｃａｎｃｅｒ  ｂｙ  ｕｓｉｎｇ  Ｒａｍａｎ  ｓｐｅｃｔｒｏｓｃｏｐｙ，  ＭＩＴ，  Ｈａｋａ，  Ｐｒｏｃ  Ｎａｔｌ  Ａｃａｄ  Ｓｃｉ  ＵＳＡ，  ２００５
【非特許文献３】Ａｂｓｏｒｐｔｉｏｎ  ｐｒｏｐｅｒｔｉｅｓ  ｏｆ  ｂｒｅａｓｔ：  ｔｈｅ  ｃｏｎｔｒｉｂｕｔｉｏｎ  ｏｆ  ｃｏｌｌａｇｅｎ，  ＵＬＴＲＡＳ−ＣＮＲ−ＩＮＦＭ  ａｎｄ  ＩＦＮ−ＣＮＲ，  Ｐｏｌｉｔｅｃｎｉｃｏ  ｄｉ  Ｍｉｌａｎｏ，  ２００６
【非特許文献４】Ａｓｓｅｓｓｍｅｎｔ  ｏｆ  ｃｏｌｌａｇｅｎ  ａｂｓｏｒｐｔｉｏｎ  ａｎｄ  ｒｅｌａｔｅｄ  ｐｏｔｅｎｔｉａｌ  ｄｉａｇｎｏｓｔｉｃ ａｐｐｌｉｃａｔｉｏｎｓ，  ＳＰＩＥ−ＯＳＡ  Ｖｏｌ．  ６６２９  ６６２９０Ｄ−１
【発明の開示】
【発明が解決しようとする課題】
【０００７】
しかし、ヘモグロビン法は腫瘍特性精度が低く、ラマン分光法は被検体を切開して組織を抽出するために被検体に負担をかける。
【０００８】
非特許文献３に開示された手法は、幅広い波長領域の光を使用するために測定時間が長く実用的ではない。また、非特許文献４はコラーゲンを用いて見積もられたヘモグロビンにより生体組織の状態を判定しているだけであり、直接コラーゲンを用いて生体組織の状態を判定しているわけではない。
【０００９】
このように、従来の測定方法は、いずれも被検体の生体組織の状態を高精度かつ簡単に（即ち、測定に長時間かけたり、被検体に切開などの負担をかけたりすることなく）判定することができなかった。
【００１０】
本発明は、被検体の生体組織の状態を高精度かつ簡単に判定する測定装置に関する。
【課題を解決するための手段】
【００１１】
本発明の一側面としての測定装置は、それぞれが６００ｎｍ乃至１０００ｎｍの波長領域のうちの異なる波長を有する複数の種類の光を被検体に照射することによって、前記被検体の内部の吸収特性と散乱特性のうち少なくとも一方を測定する測定部と、コラーゲンと脂肪を含む複数の構成要素の全体に対するコラーゲンと脂肪のそれぞれの割合を、各構成要素の吸収係数と前記測定部の測定結果からフィッティング手法を用いて算出し、脂肪及びコラーゲンのそれぞれのフィッティング係数と生体組織の状態の関係と、算出されたコラーゲンの割合と脂肪の割合から前記被検体の生体組織の状態を判定する演算処理部と、前記演算処理部による処理の結果を表示する表示部と、を有することを特徴とする。
【００１２】
本発明の更なる目的又はその他の特徴は、以下、添付図面を参照して説明される好ましい実施例によって明らかにされるであろう。
【発明の効果】
【００１３】
本発明によれば、被検体の生体組織の状態を高精度かつ簡単に判定する測定装置を提供することができる。
【発明を実施するための最良の形態】
【００１４】
本発明の実施例の測定装置は、測定部と、演算処理部と、表示部と、を有する。
【００１５】
測定部は、それぞれが６００ｎｍ乃至１０００ｎｍの波長領域のうちの異なる波長を有する複数の種類の光を吸収散乱体としての被検体に照射することによって、被検体の内部の分光特性を測定する。このように測定部は、近赤外光の生体の窓と呼ばれる波長領域を利用し、特許文献２のような被検体の切開を必要としない。
【００１６】
測定部は、後述するように、拡散光トモグラフィ（ＤＯＴ）、音響光学トモグラフィ（ＡＯＴ）、光音響トモグラフィ（ＰＡＴ）のいずれの測定方法も使用することができる。ＤＯＴ（Ｄｉｆｆｕｓｅｄ Ｏｐｔｉｃａｌ Ｔｏｍｏｇｒａｐｈｙ）は、近赤外光を被検体に導入してその拡散光を検出する。入射光は光源の光の強度を変調してもよいし、パルス光を利用してもよい。ＡＯＴ（Ａｃｏｕｓｔｏ−Ｏｐｔｉｃａｌ Ｔｏｍｏｇｒａｐｈｙ）は、被検部位にコヒーレント光及び集束超音波を照射し、超音波集束領域で光が変調される効果（音響光学効果）を利用し、変調光を光検出器で検出する。ＰＡＴ（Ｐｈｏｔｏ−Ａｃｏｕｓｔｉｃ Ｔｏｍｏｇｒａｐｈｙ）は、被検部位とそれ以外の組織との光エネルギーの吸収率差を利用し、被検部位が照射された光エネルギーを吸収して瞬間的に膨張する際の超音波（光音響信号）をトランスデューサで受信する。
【００１７】
演算処理部は、コラーゲンと脂肪を含む複数の構成要素の全体に対するコラーゲンと脂肪のそれぞれの割合を、各構成要素の吸収係数と測定部の測定結果とから算出する。次に、演算処理部は、脂肪及びコラーゲンのそれぞれのフィッティング係数と生体組織の状態の関係と、算出されたコラーゲンの割合と脂肪の割合から被検体の生体組織の状態を判定する。以下の実施例は、被検体の生体組織を、正常組織、繊維腫瘍、線維嚢胞性変化及び浸潤癌のいずれの状態であるかを判定するが、本発明では被検体の生体組織の状態の種類と数は限定されない。
【００１８】
表示部は、演算処理部による処理の結果を表示する。
【００１９】
以下、測定装置の実施例について説明する。
【実施例１】
【００２０】
図１は、本発明の実施例１の測定装置１００のブロック図である。被検体Ｅは、吸収散乱体としての測定対象物であり、具体的には、乳房などの生体組織である。図２は、測定装置１００を被験者Ｂの被検体（乳房）Ｅの乳癌の検出に適用した構成の概略断面図である。
【００２１】
測定装置１００の測定部は、信号発生器１、光源２、光ファイバ３及び１１、測定容器４、光検出器１２、信号抽出部１３を有する。
【００２２】
被検体Ｅは測定容器４に収められる。被検体Ｅと測定容器４の間は、光の屈折率や散乱係数及び、超音波の音響特性が被検体Ｅにほぼ同等とみなせる、特性が既知で均一の媒質（マッチング材５）で満たされる。半導体レーザーなどを利用した光源２からは、信号発生器（正弦波発信器）１を用いて周波数ｆ１で強度変調された光が射出される。一般に生体計測においては、数十〜数百ＭＨｚの正弦波で変調することが望ましい。光源２からの変調光は、光ファイバ３を通して測定容器４の側面から入射される。容器内に入射された変調光は、光拡散理論から導かれるように、変調周波数ｆ１のエネルギー密度波（拡散光子密度波）Ｗとして被検体Ｅの中を伝播する。この変調周波数ｆ１で伝わる拡散光子密度波Ｗを信号光として、光ファイバ１１を通してＰＭＴ（Ｐｈｏｔｏ  ＭｕｌｔｉＰｌｉｅｒ）やＡＰＤ（Ａｖａｌａｎｃｈｅ  Ｐｈｏｔｏ  Ｄｉｏｄｅ）などの光検出器１２で検出する。拡散光子密度波Ｗから必要な情報を信号抽出部１３で抽出する。
【００２３】
信号抽出部１３が抽出した分光特性の情報に基づいて、演算処理部１４は、被検体内部の吸収係数と散乱係数を算出し、吸収係数と散乱係数から被検体Ｅの構成要素の割合を算出する。演算処理部１４は、算出した吸収係数と散乱係数及び構成要素の割合をメモリ１５に保存する。画像生成部１６は、保存した値をマップ化する。表示部１７は、３次元的な吸収係数と散乱係数および構成要素の割合の分布を表示する。
【００２４】
図３は、被検体Ｅの一断面の断層像を得るための測定装置１００の動作を説明するためのフローチャートである。まず、ステップ１０１では、信号発生器１が、光源２を周波数ｆ１（数十〜数百ＭＨｚ）で強度変調させて、ある波長の光源として駆動させる。ステップ１０２で、光源２は、光ファイバ３を介して、ある位置から光を被検体Ｅに入射させる。次にステップ１０３で、光検出器１２は、光ファイバ１１を介して、変調光の振幅ＩＡＣ（ｒ，ｔ）と位相Φ（ｒ，ｔ）を計測する。一般的に拡散光子密度波を伝播する拡散光子密度波Ｗは紡錘形状となり、得られる吸収係数と散乱係数はこの光路上の平均の値となる。
【００２５】
得られた振幅と位相からステップ１０４で吸収係数μａと等価散乱係数μｓ’を求める。なお、本実施例は、光拡散方程式から導かれる近似解を用いるが、本発明は、より厳密な解にも適用することができる。
【００２６】
均一な吸収散乱体内の点光源から距離ｒの位置で時刻ｔにおける光強度Ｉ（ｒ，ｔ）［光子／ｓｅｃ・ｍｍ２］は次式で与えられる。
【００２７】
【数１】
0005219440.tif 000002
【００２８】
但し、
【００２９】
【数２】
0005219440.tif 000003
【００３０】
【数３】
0005219440.tif 000004
【００３１】
【数４】
0005219440.tif 000005
【００３２】
ここで、ＩＤＣは検出光強度のバイアス成分［光子／ｓｅｃ・ｍｍ２］、εは任意の位相項、Ａは第二の光源における入射光子数［光子／ｓｅｃ］、Ｄは拡散係数［ｍｍ］、νは吸収散乱体内部の光速度［ｍｍ／ｓｅｃ］である。また、μａは吸収係数［ｍｍ−１］、μＳ’は等価散乱係数［ｍｍ−１］、ωは超音波で変調された変調波の角周波数［ｒａｄ／ｓｅｃ］（変調周波数をｆとするとω＝２πｆ）である。
【００３３】
振幅ＩＡＣ（ｒ，ｔ）と位相Φ（ｒ，ｔ）を測定すれば、数式２及び３より、吸収係数μａと等価散乱係数μｓ’を算出することができる。
【００３４】
次にステップ１０５で、光ファイバ３の被検体Ｅに対する位置をΔｄだけずらし、同様に拡散光子密度波Ｗを光検出器１２で検出する。あるいは、予めΔｄだけずらした位置に光ファイバ３を複数取り付けておき、拡散光子密度波Ｗを光検出器１２で検出してもよい。このようにΔｄだけずらすことにより被検体Ｅに対する拡散光子密度波の光路がずれ、ずれた位置での吸収係数及び散乱係数を算出することが可能となる。以上の手順を繰り返し、被検体Ｅの吸収・散乱特性をステップ１０７でマッピングすれば、一断面の断層像を得ることができる。これは、光ファイバ１１と被検体Ｅの一方又は両方の位置をずらして測定しても同様に被検体Ｅの一断面の断層像を得ることができる。この断面を紙面鉛直方向に走査することで、最終的に被検体Ｅの３次元的な吸収・散乱情報を取得可能となる。
【００３５】
次に、ステップ１０８では、光源の波長を変える。この方法については、波長の異なる複数の光源を別々の場所から入射させても、ファイバで束ね一箇所から入射させてもよく、あるいは白色の光源を用意し、回折格子を用いて分光させたり、波長フィルターを用いて特定の波長だけを選択したりしてもよい。これらのいずれかの方法によって別の波長で被検体Ｅに変調光を入射させ、上記のステップ１０１から１０７の操作を繰り返し、測定された振幅ＩＡＣ（ｒｉ，θｊ）と位相Φ（ｒｉ，θｊ）、及び算出されたμａとμｓ’は、ステップ１０９で各波長ごとにメモリ１５に随時記憶される。用いる波長数を被検体Ｅの構成要素の数と同じにすれば一意に求まる。しかし、測定したμａとμｓ’に何らかの誤差（光源２の出力のぶれ、光源２と光検出器１２の位置のぶれ、外乱光によるぶれなど）で測定値が真値とずれてしまうと構成割合の見積もり誤差となる。そこである程度構成要素の数よりも多い波長を入射させる。
【００３６】
次にステップ１１０で、演算処理部１４は、後述する方法でフィッティングを行い、測定対象物の構成割合、特にコラーゲン、脂肪、水の割合を求める。
【００３７】
測定する波長は、求めたい構成要素に応じて変えるのが一般的である。それは各々の構成要素の吸収係数が、波長によって独特の吸収係数のスペクトル（以下、「吸収スペクトル」と呼ぶ）を示すため、構成要素に応じた特徴的な波長を用いて測定することで見積もり精度が向上するためである。
【００３８】
図４は、生体組織の主要な構成要素であるコラーゲン、脂肪、水、還元ヘモグロビン、酸化ヘモグロビンの吸収スペクトルである。例えば、還元ヘモグロビン及び酸化ヘモグロビンの吸収スペクトルに着目する。両者は約８００ｎｍでお互いが交差し、還元ヘモグロビンは約７６０ｎｍで鋭いピーク（周囲の波長よりも吸収係数が高い箇所）が存在するが、酸化ヘモグロビンは約９２０ｎｍで緩やかなピークがあるのみである。また、両者は短波長になるほど大きくなり、逆に長波長にいくほど小さくなる傾向にある。
【００３９】
従来のヘモグロビン法は、この特徴に鑑みて両者の特定を試みており、６００ｎｍから７００ｎｍと８００ｎｍ前後の波長を主として使用する。また、還元ヘモグロビン及び酸化ヘモグロビンのみならず、水、脂肪を特定する際には９００ｎｍ以上の波長が用いる。これは、前述の特徴的なピークが９００ｎｍ以上に存在するためである。非特許文献１が、６６１ｎｍ、７６１ｎｍ、７８５ｎｍ、８０８ｎｍ、８２６ｎｍ、８４９ｎｍの波長を用いるのもかかる理由によるものと考えられる。
【００４０】
一方、コラーゲンを特定するためには６００ｎｍから７００ｎｍと８００ｎｍ前後及び９００ｎｍ以上の波長では精度良く見積もれない。これは、図３に示すコラーゲンの吸収スペクトルがこの波長範囲でピークがないからである。
【００４１】
但し、他の構成要素の吸収スペクトルとの比較により、以下の特徴的な波長を見出すことができる。第１に、６００ｎｍから７００ｎｍの範囲で他の吸収スペクトルと比較してコラーゲンの吸収係数が大きい。第２に、７３０ｎｍから７６０ｎｍの範囲で他の吸収スペクトルと比較して、コラーゲンの吸収係数のみが負の傾き（波長が大きくなると吸収係数が小さくなる）となっている。
【００４２】
図５は、コラーゲンの吸収スペクトルに関する実測値と非特許文献３に記載された値との比較を示すグラフである。同図を参照するに、絶対値の違いが大きいものの、上記２点において、両者に違いはなくコラーゲンが特定できることが理解される。更に、図４から分かるように１０００ｎｍ付近で水の吸収係数が大きくなるため、被検体Ｅ中に水が多く含まれる場合には水による吸収が顕著になり、生体組織深部の測定が困難になる。
【００４３】
コラーゲンのみの構成割合を誤差に埋もれにくく見積もるために、６００ｎｍ乃至７００ｎｍの波長を有する光と、７３０ｎｍ乃至７６０ｎｍの波長領域のうちの少なくとも２つの異なる波長を有する光を使用することが好ましい。一般に、求めたい構成要素の吸収係数が大きい波長領域で測定すると、その構成要素の敏感度が高く見積もり精度が向上し易いために６００ｎｍ乃至７００ｎｍの波長領域の波長を選択する。なお、吸収係数が大きくなると被検体から射出される光量も小さくなるため、吸収係数が大きくなりすぎても測定精度が低下することに注意する必要がある。また、上記第２の特徴と傾きを見るためには最低２波長が必要であることから、７３０ｎｍ乃至７６０ｎｍの波長領域のうちの少なくとも２つの異なる波長を有する光を使用する。それらの波長で測定すると、マッピングしたある部分の吸収係数μａの波長ごとのデータをμａ＿λと置き、上記３波長で測定すると、３つのμａ、即ち、
【００４４】
【数５】
0005219440.tif 000006
【００４５】
がメモリ１５に記憶される。構成割合をＣｃｏｎｔｅｎｔとすると、
【００４６】
【数６】
0005219440.tif 000007
【００４７】
と表すことができる。これに、図４に示す既知の構成要素の吸収係数μａ＿λ＿ｃｏｎｔｅｎｔの分布を乗ずることで吸収係数を得ることができる。
【００４８】
【数７】
0005219440.tif 000008
【００４９】
数式７を数式５に最小二乗法などのフィッティング手法でＣｃｏｎｔｅｎｔをフィッティングすることで、測定としては全ての構成割合を求めることが可能である。ただし上記波長の測定範囲においては、コラーゲン以外の求めた構成割合に誤差が多い。それは前述した通りコラーゲンのみに特徴的な波長を用いたことによるものである。一方コラーゲンのみ、誤差に埋もれにくい構成割合を算出することができ、コラーゲンの量あるいは全体の量に対する割合を求めることが可能となる。
【００５０】
以上のように、被検体Ｅの一断面の一部分についてコラーゲンの量あるいは全体の量に対する割合を求めた後、マッピングされた吸収係数μａを用いて同様の手順をとる。これによって、被検体Ｅの一断面のコラーゲンの量あるいは全体の量に対する割合をマッピングすることができる。更に、この断面を走査することで、最終的に被検体Ｅの３次元の吸収・散乱情報およびコラーゲンの量あるいは全体の量に対する割合の分布を取得できる。
【００５１】
このコラーゲンの割合をマッピングした例を図６に示す。被検体内部のそれぞれの場所で算出したＣｃｏｌｌａｇｅｎの値を対応する位置に表示し、色の濃さによって値の大小を定めた。このように表示することで、コラーゲンがどの位置に多く、どの位置に少ないかを一目で確認することができる。
【００５２】
以上の波長に加えて９００ｎｍ乃至１０００ｎｍの波長領域で２波長を用いることで、コラーゲン、脂肪及び水の量あるいは全体の量に対する割合を誤差に埋もれにくく求めることができる。これは、図４に示すように、脂肪・水の吸収係数は９００ｎｍ以上で高くなっているからである。例えば、９１０ｎｍと９７０ｎｍを付加して、６４０ｎｍ、７３０ｎｍ、７６０ｎｍ、９１０ｎｍ及び９７０ｎｍを用いて測定する。マッピングしたある部分の吸収係数μａの波長ごとのデータをμａ＿λとおき、上の５波長で測定すると、以下の５つのμａがメモリ１５に保存される。
【００５３】
【数８】
0005219440.tif 000009
【００５４】
構成割合をＣｃｏｎｔｅｎｔとすると次式で表される。
【００５５】
【数９】
0005219440.tif 000010
【００５６】
数式９を図３に示す既知の構成要素の吸収係数μａ＿λ＿ｃｏｎｔｅｎｔの分布を乗ずることで吸収係数の値となる。
【００５７】
【数１０】
0005219440.tif 000011
【００５８】
数式１０を数式８に最小二乗法などのフィッティング手法でＣｃｏｎｔｅｎｔをフィッティングすることで、コラーゲン、脂肪、水のみ誤差に埋もれにくい構成割合を算出することができ、コラーゲン、脂肪、水の量あるいは全体の量に対するそれらの割合を求めることが可能である。
【００５９】
以上のように、被検体Ｅの一断面の一部分についてコラーゲン、脂肪、水の量あるいは全体の量に対するそれらの割合を求めた後、マッピングされた吸収係数μａを用いて同様の手順をとる。これによって、被検体Ｅの一断面のコラーゲン、脂肪、水の量あるいは全体の量に対するそれらの割合をマッピングすることができる。更に、この断面を走査することで、最終的に被検体Ｅの３次元の吸収・散乱情報およびコラーゲン、脂肪、水の量あるいは全体の量に対するそれらの割合の分布を取得できる。
【００６０】
図７は、脂肪の割合と水の割合をマッピングした例を示す。被検体Ｅ内部のそれぞれの場所で算出したＣｌｉｐｉｄとＣｗａｔｅｒの値を対応する位置に表示し、色の濃さによって値の大小を定めた。このように表示することで、脂肪と水がどの位置に多く、どの位置に少ないかを一目で確認することができる。
【００６１】
更に、測定装置１００は、コラーゲンと脂肪の割合を評価することで病変部位の確認も行うことが可能となる。図８は、脂肪のフィッティング係数（Ｆｉｔ  Ｃｏｅｆｆｉｃｉｅｎｔ：ＦＣ）を横軸にし、コラーゲンのフィッティング係数を縦軸にしたときの生体組織の状態の関係を示すグラフである。ＦＣはフィッティングしたときの係数で、全体の量を１としたときの割合で値を決めている。図８により、脂肪のＦＣとコラーゲンのＦＣの組み合わせによって被検体Ｅの生体組織の状態を判定する。被検部位の脂肪とコラーゲンのＦＣが見積もられたとき、図８にその値をプロットすることによって生体組織の状態を判定する。生体組織は、Ｎｏｒｍａｌ（正常組織）、Ｆｉｂｒｏｃｙｓｔｉｃ  ｃｈａｎｇｅ（線維腺腫）
、Ｆｉｂｒｏａｄｅｎｏｍａ（線維嚢胞性変化）及びＩｎｆｉｌｔｒａｔｉｎｇ  Ｃａｒｃｉｎｏｍａ（浸潤癌）のいずれか状態であるかを判定する。図９は、状態ごとに色分け及びマッピングした結果を示す。表示部１７が図９に示す被検体Ｅの一部又は全部の状態の分布を表示することによって、被検体Ｅの腫瘍の位置及び種類を確認することができる。
【００６２】
非特許文献２は、腫瘍は正常組織と比較するとコラーゲンが約４倍、脂肪が１／４になると記載している。比較対象として周りの組織との比較や、過去の同じ分布との比較、肺や乳房など左右で同じ構成を取る部位の場合は片方の部位と比較して、コラーゲンが約４倍かつ脂肪が１／４倍になっていれば腫瘍と疑うことができる。腫瘍と正常組織で４倍変化するという量が、ヘモグロビン法の場合の約１．７倍よりもはるかに大きく誤差に埋もれにくく測定精度を高めるという利点も有する。
【００６３】
コラーゲンのみを見積もる波長に、７６０ｎｍ乃至８５０ｎｍの波長領域のうちの２波長を用いることで、コラーゲン、還元ヘモグロビン、酸化ヘモグロビンの量あるいは全体の量に対する割合を誤差に埋もれにくく求めることができる。これは、図４に示すように、還元ヘモグロビンと酸化ヘモグロビンの交差点が８００ｎｍ近傍に存在し、その前後の波長で測定することにより、両者の変化を検出しやすいからである。
【００６４】
例えば、８００ｎｍと８５０ｎｍの２波長を加えて、６４０ｎｍ、７３０ｎｍ、７６０ｎｍ、８００ｎｍ、８５０ｎｍを用いて測定する。マッピングしたある部分の吸収係数μａの波長ごとのデータをμａ＿λとおいて上記５波長で測定すると、数式１１に示す５つのμａがメモリ１５に保存される。また、構成割合をＣｃｏｎｔｅｎｔは数式１２のように表すことができる。
【００６５】
【数１１】
0005219440.tif 000012
【００６６】
【数１２】
0005219440.tif 000013
【００６７】
数式１２に、図３に示す既知の構成要素の吸収係数μａ＿λ＿ｃｏｎｔｅｎｔの分布を乗ずることで吸収係数の値となる。
【００６８】
【数１３】
0005219440.tif 000014
【００６９】
数式１３を数式１１に最小二乗法などのフィッティング手法でＣｃｏｎｔｅｎｔをフィッティングすることで、コラーゲン、還元ヘモグロビン、酸化ヘモグロビンのみ誤差に埋もれにくい構成割合を算出することができ、コラーゲン、還元ヘモグロビン、酸化ヘモグロビンの量あるいは全体の量に対するそれらの割合を求めることができる。
【００７０】
以上のように、被検体Ｅの一断面の一部分についてコラーゲン、還元ヘモグロビン、酸化ヘモグロビンの量あるいは全体の量に対するそれらの割合を求めたのち、マッピングされた吸収係数μａを用いて同様の手順をとる。これによって、被検体Ｅの一断面のコラーゲン、還元ヘモグロビン、酸化ヘモグロビンの量あるいは全体の量に対するそれらの割合をマッピングすることができる。更に、この断面を走査することで、最終的に被検体Ｅの３次元の吸収・散乱情報及びコラーゲン、還元ヘモグロビン、酸化ヘモグロビンの量あるいは全体の量に対するそれらの割合の分布を取得できる。
【００７１】
図１０は、還元ヘモグロビン、酸化ヘモグロビンの割合をマッピングした結果である。被検体Ｅの内部のそれぞれの場所で算出したＣＨｂＯ２とＣＨｂの値を対応する位置に表示し、色の濃さによって値の大小を定めた。表示部１７がこのように表示することで、還元ヘモグロビン、酸化ヘモグロビンがどの位置に多く、どの位置に少ないかを一目で確認することができる。
【００７２】
このように、コラーゲン、還元ヘモグロビン及び酸化ヘモグロビンの量あるいは全体の量に対するそれらの割合の分布を得ると、ヘモグロビン法のみ又はコラーゲンで癌を検出する方法のみを使用するよりも腫瘍の検出精度が高まる。
【００７３】
本実施例では、例えば６４０ｎｍ、７３０ｎｍ、７６０ｎｍ、８００ｎｍ、８５０ｎｍ、９１０ｎｍ及び９７０ｎｍの波長を有する複数の種類の光を用いて被検体Ｅの分光特性を測定する。これにより、コラーゲン、脂肪、水、還元ヘモグロビン及び酸化ヘモグロビンの量あるいは全体の量に対する割合を誤差に埋もれにくく求めることができる。マッピングしたある部分の吸収係数μａの波長ごとのデータをμａ＿λと置き、上記７波長で測定すると、数式１４で示す７つのμａがメモリ１５に保存される。
【００７４】
【数１４】
0005219440.tif 000015
【００７５】
また、求めたい構成割合をＣｃｏｎｔｅｎｔとすると、数式１５の表すことができる。
【００７６】
【数１５】
0005219440.tif 000016
【００７７】
これに、図３に示す既知の構成要素の吸収係数μａ＿λ＿ｃｏｎｔｅｎｔの分布を乗ずることで、数式１６に示すように、吸収係数の値となる。
【００７８】
【数１６】
0005219440.tif 000017
【００７９】
数式１６を数式１４に最小二乗法などのフィッティング手法でＣｃｏｎｔｅｎｔをフィッティングすると、高精度にコラーゲン、還元ヘモグロビン、酸化ヘモグロビン、脂肪、水の量又は全体の量に対するそれらの割合を求めることができる。
【００８０】
以上のように、被検体Ｅの一断面の一部分についてコラーゲン、還元ヘモグロビン、酸化ヘモグロビン、脂肪、水の量あるいは全体の量に対するそれらの割合を求める。その後、マッピングされた吸収係数μａを用いて同様の手順をとることによって、被検体Ｅの一断面のコラーゲン、還元ヘモグロビン、酸化ヘモグロビン、脂肪、水の量あるいは全体の量に対するそれらの割合をマッピングする。更に、この断面を走査することで、最終的に被検体Ｅの３次元の吸収・散乱情報およびコラーゲン、還元ヘモグロビン、酸化ヘモグロビン、脂肪、水の量あるいは全体の量に対するそれらの割合の分布を取得できる。
【００８１】
コラーゲン、還元ヘモグロビン、酸化ヘモグロビン、脂肪、水の量又は全体の量に対するそれらの割合の分布を得ることによって、腫瘍の位置と量を検出でき、腫瘍の良性・悪性を判別し、更には腫瘍自体の検出精度を増すことができる。
【００８２】
還元ヘモグロビン・酸化ヘモグロビンを用いた正常か腫瘍かの判定は、数式１７で表される酸素飽和度（ＳＴＯ２）を用いて行われる。なお、［Ｘ］はＸの１リットル辺りのモル濃度を表す。
【００８３】
【数１７】
0005219440.tif 000018
【００８４】
従来は、酸素飽和度で正常か腫瘍かの判定が行われてきたが、その変化量は小さく誤差に埋もれ易く、腫瘍の有無を判定する精度が低い。
【００８５】
そこで、本実施例は、コラーゲンと脂肪のＦＣによる判定を主要な判定として使用し、さらに酸素飽和度を用いた判定を補助的判定として使用することで、正常と腫瘍の判定の精度を上げることができる。場所ごとのコラーゲンの量と脂肪の量を求め、生体組織の状態を判定したのと同様に、場所ごと還元ヘモグロビンと酸化ヘモグロビンの量から数式１７に従って酸素飽和度を求め、その酸素飽和度によっても判定を行う。
【００８６】
図１１は、酸素飽和度の値ごとに色分け及びマッピングした結果であり、表示部１７が表示する。
【００８７】
図１２は、生体組織の状態の分布と酸素飽和度の分布を同時に重ねて表示してマッピングした結果であり、表示部１７が表示する。この２種類の分布の重なり領域がより腫瘍の可能性が高い領域となる。この表示によって、医師や患者に腫瘍の可能性が高い領域を高精度にかつ簡単に提供することができる。その際、この２種類の分布のうちの最初の分布は信頼性が高いので重みを与えて表示してもよい。
【００８８】
その他、水の量や還元ヘモグロビンと酸化ヘモグロビンの量の合計や、求めた散乱係数を表示部１７に表示させて、腫瘍の可能性の判断の補助に用いてもよい。例えば、散乱係数が高い領域は、腫瘍になる可能性が高いマンモグラフィ濃度が高い領域と比較的一致する。そのため、診断において表示部１７に表示された散乱係数の分布の見ることで、付加的な情報を得ることができる。
【００８９】
一般に、波長数をｎとすると、数式１８に示すように、フィッティングを行う。
【００９０】
【数１８】
0005219440.tif 000019
【００９１】
以上のように、被検体Ｅの一断面の一部分について構成割合を求めたのち、マッピングされた吸収係数μａを用いて同様の手順をとることによって、被検体Ｅ一断面の構成要素をマッピングすることができる。更に、この断面を走査することで、最終的に被検体Ｅの３次元の吸収・散乱情報および構成要素の分布を取得できる。
【００９２】
その他、光源２を変調する代わりに、数ピコ秒の短パルスを被検体Ｅに入射させ、出力される時間波形からμａとμｓ’を見積もることもできる。取得したμａとμｓ’から構成要素を見積もる方法は前述した通りである。この手法を用いてもなお被検体Ｅの３次元の吸収・散乱情報および構成要素の分布を取得できる。また、被検体Ｅに一定の強度の光を照射し、その透過光強度からμａを求め、取得したμａから構成要素を見積もることも可能である。
【実施例２】
【００９３】
図１３は、本発明の実施例２の測定装置１００Ａのブロック図である。測定装置１００ＡはＡＯＴを利用して被検体Ｅの分光特性を測定する。図１３において、図１と同様の部材には同一の参照符号を付している。測定装置１００Ａも図２に示す測定装置１００の代わりに使用することができる。
【００９４】
レーザーなどのコヒーレントな光が光源２から連続的に射出される。射出光は、光ファイバ３を通って測定容器４の側面から入射される。容器内に入射された光は、媒質内部を吸収と散乱を繰り返しながら伝播する。
【００９５】
また、信号発生器１でｆの正弦波信号を用いて、測定容器４の底面に配置された超音波トランスデューサアレイ７を駆動する。超音波トランスデューサアレイ７を制御して、測定容器４の局所的な領域である被検部位に音圧が集束するように超音波を照射する。超音波集束領域Ｐでは、音圧による媒質の密度変化が生じ、媒質の屈折率や散乱係数が変化する。
【００９６】
超音波集束領域Ｐを光が通過すると、媒質の屈折率や散乱係数の変化により、光の位相が変調される。超音波集束領域Ｐで変調を受けると、超音波の駆動周波数ｆで変調された光として以後媒質中を伝播する。
【００９７】
光検出器１２には、超音波によって変調された光と変調されていない光が入射し、変調光と非変調光の両方からの信号が検出される。
【００９８】
この信号を、信号抽出部１３が検出信号をフーリエ変換し、非変調信号Ｉ1と超音波の周波数ｆによる変調信号I2とに分離する。分離した信号と参照信号を利用して、測定対象の内部の吸収係数と散乱係数を算出する。同じく演算処理部１４で吸収係数と散乱係数から被検体Ｅの構成要素の割合を算出し、算出した吸収係数と散乱係数および構成要素の割合はメモリ１５に保存される。その保存された値を画像生成部１６でマップ化され、表示部１７で３次元的な吸収係数と散乱係数および構成要素の割合の分布を表示する。
【００９９】
図１４は、被検体Ｅの一断面の断層像を得るための測定装置１００Ａの動作を説明するためのフローチャートである。まず、ステップ２０１で、信号発生器１が光源２を駆動する。ステップ２０２で、光源２は光ファイバ３を介して、光を被検体Ｅに入射させる。ステップ２０３で、超音波トランスデューサアレイ７は、被検体Ｅに超音波を照射及び集束させる。ステップ２０４で光検出器１２は光を検出する。ステップ２０５で、光検出器１２が検出した光強度から、数式２及び３を用いて吸収係数μａと等価散乱係数μｓ’を算出する。
【０１００】
次に、ステップ２０６で、超音波トランスデューサアレイ７を制御して音圧の集束位置をずらす。次に、ステップ２０７で、光と超音波の相互作用によってタグ付けされた局所領域（超音波集束領域Ｐ）を測定容器４の内部の吸収散乱体全域にくまなく配置する。被検体Ｅの既知の領域と未知の領域の差分をとり帰納的に求めることで、光と超音波の相互作用によってタグ付けされた局所的な領域の吸収係数・散乱係数を取得することが可能となる。ステップ２０８が、それらの値をマッピングすれば、被検体Ｅの一断面の断層像が容易に得ることができる。
【０１０１】
次に、ステップ２０９が駆動させる光源の波長を変える。この方法については、波長の異なる複数の光源を別々の場所から入射させても、ファイバで束ね一箇所から入射させてもよい。あるいは白色の光源を用意し、回折格子を用いて分光させたり、波長フィルターを用いて特定の波長だけを選択したりしてもよい。これらのいずれかの方法によって別の波長の光を被検体Ｅに入射させ、ステップ２０１から２０７の操作を繰り返し、、ステップ２１０で各波長ごとに算出されたμａとμｓ’メモリ１５に記憶する。用いる波長数を被検体Ｅの構成要素の数よりも多くする点は実施例１と同様である。
【０１０２】
入射させる波長は実施例１と同様に、検出する構成要素によって決定される。測定する波長及びそこから構成要素を求める手段は実施例１に記載したとおりである。本実施例では、光と超音波の相互作用によってタグ付けされた局所的な領域の水、脂肪、酸化ヘモグロビン、還元ヘモグロビン、コラーゲンの構成割合を取得することが可能となる。この断面を走査することで、最終的に被検体Ｅの３次元の吸収・散乱情報（ステップ２１０）及び構成要素の分布（ステップ２１１）を取得することができる。
【実施例３】
【０１０３】
図１５は、本発明の実施例３の測定装置１００Ｂのブロック図である。測定装置１００ＢはＰＡＴを利用して被検体Ｅの分光特性を測定する。図１６において、図１と同様の部材には同一の参照符号を付している。測定装置１００Ｂも図２に示す測定装置１００の代わりに使用することができる。
【０１０４】
半導体レーザーなどを利用した光源２からは、信号発信器１を用いてナノ秒オーダーのパルス光が射出される。光源２からのパルス光は、光ファイバ３を通して測定容器４の側面から入射される。容器内に入射された光は被検体Ｅの中を伝播し、伝播した光が吸収体（被検部位）Ｒに当たると媒質の膨張収縮で弾性波が発生し、吸収体Ｒからの弾性波は被検体Ｅを伝播し、超音波トランスデューサアレイ７で検出される。
【０１０５】
信号抽出部１３において信号発生器１からの同期を取ることにより、アレイごとに検出する時間の差により吸収体Ｒの位置を求めることができる。弾性波の強度から吸収体Ｒの吸収係数μａを算出することができる。演算処理部１４が吸収係数から被検体Ｅの構成要素の割合を算出し、算出した吸収係数と散乱係数、並びに、構成要素の割合はメモリ１５に保存される。その保存された値を画像生成部１６でマップ化され、表示部１７で３次元的な吸収係数と散乱係数および構成要素の割合の分布を表示する。
【０１０６】
図１６は、被検体Ｅの一断面の断層像を得るための測定装置１００Ｂの動作を説明するためのフローチャートである。まず、ステップ３０１が光源２を駆動し、ステップ３０２がナノ秒オーダーの光を被検体Ｅに入射させる。ステップ３０３が吸収体Ｒからの超音波を検出する。ステップ３０４は、その強度からμａを算出し、ステップ３０５が時間差から位置を算出してマッピングを行う。ステップ３０６は、ステップ２０９と同様に、波長を変化させ、ステップ３０７はフィッティングにより被検体Ｅの構成割合を求める。また、算出されたμａとμｓ’と構成割合をメモリ１５に記憶する。用いる波長数を被検体Ｅの構成要素の数よりも多くする点は実施例１と同様である。
【０１０７】
測定する波長及びそこから構成要素を求める方法は実施例１と同様である。本実施例では、局所的な領域の水、脂肪、酸化ヘモグロビン、還元ヘモグロビン、コラーゲンの構成割合を取得することが可能となる。この断面を走査することで、最終的に被検体Ｅの３次元の吸収・散乱情報および構成要素の分布を取得できる。
【０１０８】
以上、本発明の好ましい実施例を説明したが、本発明はこれらに限定されずその要旨の範囲内で種々の変形及び変更が可能である。
【図面の簡単な説明】
【０１０９】
【図１】本発明の実施例１の測定装置のブロック図である。
【図２】図１に示す測定装置を乳癌の検出に適用した構成の概略断面図である。
【図３】図１に示す測定装置の動作を説明するためのフローチャートである。
【図４】生体組織の主要な構成要素であるコラーゲン、脂肪、水、還元ヘモグロビン、酸化ヘモグロビンの吸収スペクトルである。
【図５】コラーゲンの吸収スペクトルに関する実測値と非特許文献３に記載された値との比較を示すグラフである。
【図６】コラーゲンの割合をマッピングした結果である。
【図７】脂肪の割合と水の割合をマッピングした結果である。
【図８】脂肪のフィッティング係数とコラーゲンのフィッティング係数との関係を示すグラフである。
【図９】図８に示す状態を色分け及びマッピングした結果である。
【図１０】還元ヘモグロビン、酸化ヘモグロビンの割合をマッピングした結果である。
【図１１】酸素飽和度の値ごとに色分け及びマッピングした結果である。
【図１２】生体組織の状態の分布と酸素飽和度の分布を同時に表示してマッピングした結果である。
【図１３】本発明の実施例２の測定装置のブロック図である。
【図１４】図１３に示す測定装置の動作を説明するためのフローチャートである。
【図１５】本発明の実施例３の測定装置のブロック図である。
【図１６】図１５に示す測定装置の動作を説明するためのフローチャートである。
【符号の説明】
【０１１０】
２                      光源
１４                    演算処理部
１７                    表示部
１００−１００Ｂ        測定装置
【図１】
0005219440.tif 000020
【図２】
0005219440.tif 000021
【図３】
0005219440.tif 000022
【図４】
0005219440.tif 000023
【図５】
0005219440.tif 000024
【図６】
0005219440.tif 000025
【図７】
0005219440.tif 000026
【図８】
0005219440.tif 000027
【図９】
0005219440.tif 000028
【図１０】
0005219440.tif 000029
【図１１】
0005219440.tif 000030
【図１２】
0005219440.tif 000031
【図１３】
0005219440.tif 000032
【図１４】
0005219440.tif 000033
【図１５】
0005219440.tif 000034
【図１６】
0005219440.tif 000035