(57)【要約】  蛍光画像取得方法は以下のステップ（Ａ）から（Ｄ）を含む。（Ａ）試料から離れた位置から、広がる光束の励起光で試料を照射して試料中の蛍光物質を励起して蛍光を発生させ、発生した蛍光を二次元検出器で検出して蛍光画像を得るステップ、（Ｂ）前記二次元検出器により蛍光を検出しようとする時点において最も強い蛍光を発する部分からの蛍光が前記二次元検出器に入射するのを実質的に遮断した状態で、前記ステップ（Ａ）と同様に試料に励起光を照射して前記二次元検出器により蛍光画像を得る少なくとも１つのステップ、（Ｃ）ステップ（Ａ）で得た蛍光画像とステップ（Ｂ）で得た少なくとも１つの蛍光画像からなる少なくとも２つの蛍光画像を１つの画像に統合するステップ、及び（Ｄ）ステップ（Ｃ）で統合された画像を表示するステップ。
2015037055.tif 000002
【特許請求の範囲】
【請求項１】
試料台と、
前記試料台に置かれた試料を、試料から離れた位置から、広がる光束の励起光で照射し試料中の蛍光物質を励起して蛍光を発生させる励起光源と、
試料からの蛍光を検出する二次元検出器を備えた撮像部と、
試料と前記二次元検出器との間に配置され、試料からの蛍光の一部が前記二次元検出器に入射するのを遮断する遮断機構と、
前記撮像部により撮像しようとする時点において試料からの最も強い蛍光が前記二次元検出器に入射するのを実質的に遮断するように前記遮断機構を制御する遮断制御部と、
前記遮断機構によって試料の一部が遮断された状態で前記撮像部が得た少なくとも１つの蛍光画像を含む少なくとも２つの蛍光画像を１つの画像に統合するように計算処理する画像統合部と、
前記画像統合部が得た統合画像を表示する表示装置と、
を備えた蛍光画像取得装置。
【請求項２】
前記撮像部が得た蛍光画像データから試料内部を計算して蛍光画像を再構成する三次元分布計算処理部をさらに備え、
前記画像統合部は前記三次元分布計算処理部により再構成された蛍光画像を統合するように構成されている請求項１に記載の蛍光画像取得装置。
【請求項３】
前記遮断制御部は、前記三次元分布計算処理部により再構成された蛍光画像に基づいて、最も強い蛍光が前記二次元検出器に入射するのを実質的に遮断するように前記遮断機構を制御するように構成されている請求項２に記載の蛍光画像取得装置。
【請求項４】
試料の内部情報を入力する入力部をさらに備え、
前記遮断制御部は、前記入力部から入力された試料内部情報に基づいて、最も強い蛍光が前記二次元検出器に入射するのを遮断するように前記遮断機構を制御するように構成されている請求項１又は２に記載の蛍光画像取得装置。
【請求項５】
前記遮断制御部は、蛍光を受光する方向から見たときの蛍光物質の試料内部での深さに応じて遮断範囲を変化させるように前記遮断機構を制御するように構成されている請求項１から４のいずれか一項に記載の蛍光画像取得装置。
【請求項６】
前記遮断機構は、前記遮断制御部からの電気制御によって自動動作する機構を備えている請求項１から５のいずれか一項に記載の蛍光画像取得装置。
【請求項７】
前記表示装置は、統合画像を強度範囲毎に異なる色彩で表示するように構成されている請求項１から６のいずれか一項に記載の蛍光画像取得装置。
【請求項８】
以下のステップ（Ａ）から（Ｄ）を含む蛍光画像取得方法。
（Ａ）試料から離れた位置から、広がる光束の励起光で試料を照射して試料中の蛍光物質を励起して蛍光を発生させ、発生した蛍光を二次元検出器で検出して蛍光画像を得るステップ、
（Ｂ）前記二次元検出器により蛍光を検出しようとする時点において最も強い蛍光を発する部分からの蛍光が前記二次元検出器に入射するのを実質的に遮断した状態で、前記ステップ（Ａ）と同様に試料に励起光を照射して前記二次元検出器により蛍光画像を得る少なくとも１つのステップ、
（Ｃ）ステップ（Ａ）で得た蛍光画像とステップ（Ｂ）で得た少なくとも１つの蛍光画像からなる少なくとも２つの蛍光画像を１つの画像に統合するステップ、及び
（Ｄ）ステップ（Ｃ）で統合された画像を表示するステップ。
【請求項９】
ステップ（Ａ）と（Ｂ）で得た蛍光画像から試料内部を計算してそれぞれの蛍光画像を再構成する再構成ステップをさらに備え、
前記画像統合ステップ（Ｃ）は前記再構成ステップにより再構成された蛍光画像を統合する請求項８に記載の蛍光画像取得方法。
【請求項１０】
前記ステップ（Ｂ）は、前記再構成ステップで得た試料内部情報に基づいて、前記二次元検出器により蛍光を検出しようとする時点における最も強い蛍光を発する部分を決定する請求項９に記載の蛍光画像取得方法。
【請求項１１】
前記ステップ（Ｂ）は、外部から入力された試料内部情報に基づいて、前記二次元検出器により蛍光を検出しようとする時点における最も強い蛍光を発する部分を決定する請求項８又は９に記載の蛍光画像取得方法。
【請求項１２】
前記ステップ（Ｂ）は、蛍光を受光する方向から見たときの蛍光物質の試料内部での深さに応じて遮断範囲を変化させる請求項８から１１のいずれか一項に記載の蛍光画像取得方法。
【請求項１３】
前記ステップ（Ｄ）は、統合された画像を強度範囲毎に異なる色彩で表示する請求項８から１２のいずれか一項に記載の蛍光画像取得方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は試料中の蛍光物質を画像として取得する装置と方法に関し、例えば、小動物などの生体試料を対象とする生体内の蛍光画像を取得する装置と方法に関する。
【背景技術】
【０００２】
蛍光画像取得装置には大きく分けて２つのシステムが存在する。【０００３】
１つは、（Ａ）走査型蛍光画像取得装置である。この装置では、十分に集光された励起光を観察対象である試料の１点から照射して試料内に存在する蛍光物質を励起し、その蛍光物質から発生する蛍光を受光することにより蛍光画像を取得する。その照射点が試料全域に及ぶように走査することで試料全体の蛍光画像を得る（特許文献１，２参照。）。また、得られた各蛍光画像データから蛍光物質の位置を推定する断層画像化の機能をもつ装置も存在する。【０００４】
走査型蛍光画像取得装置の特徴として、観察領域を指定できることが挙げられる。例えば、生体試料の脳を観察したい場合、脳領域のみを走査点として設定すれば、脳以外からのバックグランド蛍光を除去することができ、十分なダイナミックレンジを得ることができる。しかしながら、in vivo観察においてこのようなバックグランド蛍光も重要な情報をもつ場合がある。例えば薬物代謝評価を行う場合、対象領域以外で薬剤が集積してしまうとそれが薬の副作用として影響するおそれがある。このため、まずは生体全体を走査することが必要となるが、走査型装置では観察対象全域の走査に膨大な時間がかかる。【０００５】
もう１つは、（Ｂ）全照射型蛍光取得装置である。この装置では、励起光で試料の広域を同時に照射して広域の蛍光画像を一度に得る（特許文献３参照。）。また、データ数は少ないものの、得られた各蛍光画像データから蛍光物質の位置を推定する断層画像化の機能をもつ装置も存在する。【０００６】
全照射型蛍光取得装置の特徴として、走査型の欠点であった観察対象全域の蛍光画像を短時間で取得できる点が挙げられる。例えば、マウスを利用した薬物スクリーニングや手術中のモニタリング等を目的とした使用方法においては短時間に蛍光画像を得ることが期待できる。しかしながら、全域観察においてはバックグランド蛍光の影響によって十分なダイナミックレンジを得ることができない場合がある。例えば抗体医薬品の評価を行う場合、目的とする対象以外への集積（特に肝臓集積）によって対象領域の微弱な蛍光が検出できない場合がある。また、例えばＩＣＧ（Indo Cyanine Green）を利用したリンパ節観察の場合、ＩＣＧを注入した箇所で蛍光強度が極めて強くなるため、リンパ節末端の微弱な蛍光が検出できない場合がある。【０００７】
このため、（Ｃ）全照射型に励起光の強度分布を与えてダイナミックレンジを広くする方法が考えられている（特許文献４参照。）。この方法は蛍光物質の深さ方向を考慮しない場合（ＤＮＡ検出など）においては有効であるが、マウスやヒト等の生体計測においては、媒体の深さ方向によって励起条件を変えていく必要があるので、（Ｃ）の方法を単純に転用できない。また、励起される箇所によって励起光を弱めたとしてもその箇所の蛍光物質の量が多ければ蛍光強度が大きくなってしまうため、蛍光検出側にも何らかの処理機構が必要となる。【０００８】
なお、蛍光二次元検出器にハレーションを起こさない機構を取り付けることも考えられるが、この場合、蛍光強度の線形性が担保できないために定量性解析等に問題が生じる。
【先行技術文献】
【特許文献】
【０００９】
【特許文献１】米国特許第６６１５０６３号
【特許文献２】米国特許出願公開第２０１２／００１７８４２号
【特許文献３】米国特許第８３４５９４１号
【特許文献４】米国特許第７３０７２６１号
【特許文献５】米国特許出願公開第２０１１／０１９９５００号
【発明の開示】
【発明が解決しようとする課題】
【００１０】
前述のように、走査型は広いダイナミックレンジを確保することができるが観察時間が極めて長くなり、一方、全照射型は短時間・連続時間観察が可能であるがダイナミックレンジが狭く、また相対的に微弱な蛍光が検出できない、といった問題点が挙げられる。【００１１】
本発明は生体の蛍光画像を取得する場合に限らず、生体以外の試料を対象として蛍光画像を取得する場合も含み、観察対象全域の蛍光物質の画像を広いダイナミックレンジで効率よく取得することを目的とする。
【課題を解決するための手段】
【００１２】
本発明の蛍光画像取得装置は、試料への励起光による照射を全照射型とするため、試料台と、試料台に置かれた試料を、試料から離れた位置から、広がる光束の励起光で照射し試料中の蛍光物質を励起して蛍光を発生させる励起光源と、試料からの蛍光を検出する二次元検出器を備えた撮像部と、撮像部が取得した蛍光画像に基づく画像を表示する表示装置とを備えている。【００１３】
本発明の蛍光画像取得装置は、さらに、試料と前記二次元検出器との間に配置され、試料からの蛍光の一部が前記二次元検出器に入射するのを遮断する遮断機構と、前記撮像部により撮像しようとする時点において試料からの最も強い蛍光が前記二次元検出器に入射するのを実質的に遮断するように前記遮断機構を制御する遮断制御部と、前記遮断機構によって試料の一部が遮断された状態で前記撮像部が得た少なくとも１つの蛍光画像を含む少なくとも２つの蛍光画像を１つの画像に統合するように計算処理する画像統合部とを備えている。そして、前記表示装置は、前記画像統合部が得た統合画像を表示するように構成されている。【００１４】
ここで、遮断制御部が二次元検出器に入射する最も強い蛍光を実質的に遮断するときの「実質的に遮断する」とは、最も強い蛍光の次に強い蛍光が検出できるようになる程度にまで、最も強い蛍光を遮断することを意味する。実質的に遮断することの一例として、後で図７を用いて説明する例では、遮断部分以外の領域のダイナミックレンジを５倍に高めて微弱な蛍光の検出を可能にするために、最も強い蛍光の強度を約１／５に減衰させている。その場合の約１／５に減衰させることは、実質的に遮断することに相当する。減衰率は特に限定されるものではなく、１／５よりも大きくても小さくてよい。もちろん、１００％遮断する場合も含む。【００１５】
本発明の蛍光画像取得方法は以下のステップ（Ａ）から（Ｄ）を含む。
（Ａ）試料から離れた位置から、広がる光束の励起光で試料を照射して試料中の蛍光物質を励起して蛍光を発生させ、発生した蛍光を二次元検出器で検出して蛍光画像を得るステップ、【００１６】
（Ｂ）前記二次元検出器により蛍光を検出しようとする時点において最も強い蛍光を発する部分からの蛍光が前記二次元検出器に入射するのを実質的に遮断した状態で、前記ステップ（Ａ）と同様に試料に励起光を照射して前記二次元検出器により蛍光画像を得る少なくとも１つのステップ、【００１７】
（Ｃ）ステップ（Ａ）で得た蛍光画像とステップ（Ｂ）で得た少なくとも１つの蛍光画像からなる少なくとも２つの蛍光画像を１つの画像に統合するステップ、及び【００１８】
（Ｄ）ステップ（Ｃ）で統合された画像を表示するステップ。【００１９】
ステップ（Ａ）では試料の観察対象全域での蛍光画像を取得する。観察対象全域は試料全体のこともあるし、試料の一部のこともある。この場合、試料中の蛍光物質からの蛍光は二次元検出器のダイナミックレンジに依存して検出される。【００２０】
ステップ（Ｂ）では、試料内部での蛍光物質の位置情報に基づいて、蛍光を物理的に遮断する遮断機構を作動させて試料からの最も強い蛍光を実質的に遮断する。試料内部での蛍光物質の位置情報は、例えば蛍光画像から生体内部情報を求める計算処理部によって、又は入力された試料内部情報に基づいて得ることができる。最も強い蛍光を遮断するには、最も強い蛍光を発生する蛍光物質に励起光が照射されないように励起光を遮断してもよく、又は最も強い蛍光を発生する蛍光物質に励起光が照射されて蛍光が発生してもその蛍光が二次元検出器に入射しないようにその蛍光を遮断してもよい。そして、このように撮像時点の試料からの最も強い蛍光を実質的に遮断した状態で蛍光画像を取得する。このとき、最も強い蛍光物質からの蛍光は実質的に遮断されているため、その次に蛍光強度の強い蛍光物質からの蛍光が検出されることになる。この検出は、ステップ（Ａ）での撮像条件に比べて、二次元検出器の露光時間を長くするなどの方法により、感度を上げて行う。【００２１】
さらにステップ（Ｂ）の処理を必要に応じて繰り返す。２回目のステップ（Ｂ）では、ステップ（Ａ）で検出される蛍光はすでに実質的に遮断されているので、１回目のステップ（Ｂ）で検出される蛍光が２回目のステップ（Ｂ）での最も強い蛍光となる。２回目のステップ（Ｂ）での検出は、１回目のステップ（Ｂ）での撮像条件に比べて、二次元検出器の露光時間を長くするなどの方法により、感度を上げて行う。３回目以上のステップ（Ｂ）を行う場合も同様である。【００２２】
このように、二次元検出器の撮像条件（露光時間等）を各ステップ毎に変えていくことで最終的にダイナミックレンジを拡げることができ、強い蛍光から微弱な蛍光まで線形性を担保したまま蛍光検出できるようになる。ステップ（Ａ）とステップ（Ｂ）を含む処理はせいぜい２〜５回程度の繰り返しであるため、走査型で全領域を撮像するよりも短時間ですみ、格段に効率的なプロセスである。
【発明の効果】
【００２３】
本発明は、撮像時点の試料からの最も強い蛍光が二次元検出器に入射するのを実質的に遮断した状態で撮像部が得た蛍光画像を含む少なくとも２つの蛍光画像を１つの画像に統合する。これにより、全照射型に比べて観察対象全域の蛍光を広いダイナミックレンジで検出することができ、かつ走査型に比べて観察対象全域の蛍光を短時間で検出することができるので、効率がよくなる。例えば、生体内の複数個所に蛍光物質が集積する場合や、強い蛍光に隠れて見逃すことの多い微弱な蛍光を検出する際に特に有用である。
【図面の簡単な説明】
【００２４】
【図１】蛍光画像取得装置の一実施形態を示すブロック図である。
【図２】一実施形態の動作の一例を示すフローチャートである。
【図３】一実施形態の動作の一例を示す概要図であり、（Ａ）は観察対象全域での蛍光画像を取得する状態を示し、（Ｂ）は特定の場所からの蛍光が二次元検出器に実質的に入射しないように遮断した状態で蛍光画像を取得する状態を示している。
【図４】遮断機構の第１の具体例を示す図である。
【図５】遮断機構の第２の具体例を示す図である。
【図６】遮断機構によって遮断する範囲を決定する方法を説明する図であり、（Ａ）は生体内に設置された蛍光物質と蛍光を示す断面図、（Ｂ）は生体試料表面で検出される蛍光強度のプロファイルを示す図である。
【図７】図６（Ｂ）のプロファイルを距離２０ｍｍでの蛍光強度が１となるように規格化した蛍光強度のプロファイルを示す図である。
【図８】図４に示される遮断機構の動作の一例を示すフローチャートである。
【図９】図５に示される遮断機構の動作の一例を示すフローチャートである。
【図１０】円筒型ファントムの観察対象全域の蛍光検出を示す図であり、（Ａ）はファントムとそれに埋め込まれた蛍光物質入りチューブを示す図、（Ｂ）は観察対象全域の蛍光画像を示す画像である。
【図１１】図１０（Ｂ）の蛍光画像データに基づいて再構成された蛍光画像を示す縦断面図（左）と横断面図（右）である。
【図１２】同円筒型ファントムの一部を遮断した状態での蛍光検出を示す図であり、（Ａ）はファントムとその一部を覆う遮断機構を示す斜視図、（Ｂ）は検出された蛍光画像を示す画像である。
【図１３】図１２（Ｂ）の蛍光画像データに基づいて再構成された蛍光画像を示す縦断面図（左）と横断面図（右）である。
【図１４】図１０、１２の蛍光画像を異なる２色を用いて同じ画面上に表示する方法を示す図である。
【図１５】図１１、１３の蛍光画像を異なる２色を用いて同じ画面上に表示する方法を示す図である。
【発明を実施するための形態】
【００２５】
蛍光画像取得装置の一実施形態を図１に示す。
試料台１０上に測定対象の試料１２が載せられる。試料１２は特に限定されるものではないが、例えばマウス等の小動物のような生体試料である。励起光源１４は試料台１０に置かれた試料１２を、試料１２から離れた位置から、広がる光束の励起光１６で照射し試料１２中の蛍光物質を励起して蛍光を発生させる励起光源を備えている。励起光源としては、必要な励起波長の励起光を発光するＬＤ（レーザダイオード）やＬＥＤ（発光ダイオード）のほか、ハロゲンランプなどのランプ類を用いることができる。励起光源１４には、検出しようとする蛍光波長の光を除去する励起側干渉フィルタを励起光源に組み合わせて使用するのが好ましい。
【００２６】
励起光源１４は、さらに、好ましい実施形態として三次元分布計算処理部を備えている場合には、撮像された外観画像から順問題解析処理と逆問題解析処理で使用する三次元表面形状データを得るために、試料を照明して試料の外観画像を撮像するための白色ＬＥＤなどの照明用光源も備えている。
【００２７】
撮像部１８は試料１２からの蛍光を検出する二次元検出器２０を備えている。二次元検出器２０は例えばＣＣＤ（電荷結合素子）である。撮像部１８は試料１２からの蛍光を二次元検出器２０に導くための集光レンズ等の光学系と、励起光成分が二次元検出器２０に入射するのを防ぐために励起光成分を除去し検出しようとする蛍光を透過させる蛍光側干渉フィルタも含んでいる。
【００２８】
試料１２の所定の箇所からの蛍光が二次元検出器２０に入射するのを遮断するために、遮断機構２２が設けられている。遮断機構２２の具体的な例は後述するが、遮断制御部２４からの電気制御によって自動動作する機構を備えているのが好ましい。遮断制御部２４は遮断機構２２が撮像時点において試料１２からの最も強い蛍光が二次元検出器２０に入射するのを実質的に遮断するように遮断機構２２を制御する。
【００２９】
画像統合部２６は遮断機構２２によって試料１２の一部が遮断された状態で撮像部１８が得た蛍光画像を含む少なくとも２つの蛍光画像を１つの画像に統合するように計算処理する。表示装置２８は画像統合部２６が得た統合画像を表示する。表示装置２８は、特に限定するものではないが、例えば液晶表示装置である。表示装置２８は、統合画像を強度範囲毎に異なる色彩で表示するように構成されているのが好ましい。
【００３０】
遮断制御部２４は、試料内部情報に基づいて、撮像時点の試料１２からの最も強い蛍光が二次元検出器２０に入射するのを実質的に遮断するように遮断機構２２を制御する。そのための試料内部情報を得るために、一実施形態では、撮像部１８が得た蛍光画像データから試料内部情報を得る三次元分布計算処理部３０を備えている。その場合は、遮断制御部２４は、三次元分布計算処理部３０が得た試料内部情報に基づいて遮断機構２２を制御する。
【００３１】
遮断制御部２４が遮断機構２２を制御するための試料内部情報を得るために、他の実施形態では、検出対象である試料の内部情報を入力する試料内部情報入力部３２をさらに備えいる。その場合は、遮断制御部２４は、入力部３２から入力された試料内部情報に基づいて遮断機構２２を制御する。
【００３２】
三次元分布計算処理部３０を備えている場合、画像統合部２６は三次元分布計算処理部３０により再構成された蛍光画像を統合するように構成されていてもよい。
【００３３】
遮断制御部２４、画像統合部２６及び三次元分布計算処理部３０はコンピュータ３４により実現される。コンピュータ３４はこの蛍光画像取得装置の専用コンピュータ又はパーソナルコンピュータである。入力部３２はコンピュータ３４にデータを入力するキーボードやマウスを含む入力装置である。
【００３４】
図２と図３を参照して、一実施形態の動作の一例を説明する。
【００３５】
（手順１）（図２（１））
まず始めに、励起光源１４からの励起光１６で試料１２の観察対象全域を照射し、撮像部１８の二次元検出器２０により観察対象全域での蛍光画像を取得する（図３（Ａ））。観察対象全域は必ずし試料全体でなくてもよく、試料の観察しようとする一部であってもよい。試料は特に限定されるものではないが、例えば生体試料とする。この処理では、最も蛍光強度の強い場所が感度よく検出される。図３（Ａ）の例では、記号ａで示される場所が蛍光強度の最も強い場所であり、その場所からの蛍光が感度よく検出される。一方、その場所ａからの蛍光を検出する条件で別の箇所ｂの微弱な蛍光まで感度よく撮像しようとすると、二次元検出器の検出感度を上げなければならないため、二次元検出器のダイナミックレンジが足りずに場所ａからの強い蛍光によってハレーションを起こす状況が発生する。
【００３６】
（手順２）（図２（２））
次に、前述の処理で得た蛍光画像データを元に、三次元分布計算処理部３０により生体内部情報を求める。生体内部情報を求める処理手順の具体例として、順問題解析と逆問題解析から断層画像を計算する方法や、異なる波長から深部情報を導き出す方法（特許文献５参照。）等がある。また、皮下の血管やリンパ節など集積箇所が予め判明している場合については、三次元分布計算処理部３０により生体内部情報を求めるのに代えて、試料内部情報入力部３２から生体内部情報を入力するようにしてもよい。
【００３７】
ここで、生体内部情報を求める処理手順の一例として、順問題解析と逆問題解析から断層画像を計算する方法について、その処理手順を説明する。順問題解析は、任意の方向から励起光を生体試料へ照射した場合に励起光が生体内部を伝搬する過程、及び生体内部の任意の位置に蛍光物質が存在した場合に蛍光が生体内部を伝搬する過程から成り立つ。それぞれの過程は次に示す光拡散方程式（１）、（２）によって計算することができる。
2015037055.tif 000003
ここで、
Ｄ：拡散定数、
μａ：吸収係数、
φex(j)：ｊ番目の方向から励起光を与えた場合の位置ｒにおける励起光フルエンスレート、
φem(k)：ｋ番目の位置に蛍光物質が存在した場合の位置ｒにおける蛍光フルエンスレート、
ε：モル吸光係数、
γ：量子収率、
Ｍ：モル濃度、
Ｓ（ｊ）：j番目の位置における励起光強度、
を表す。
【００３８】
（１）式と（２）式から励起・蛍光伝搬の理論値を求めることができ、この計算結果を元にシステム行列Ａを作成する。システム行列Ａとは、蛍光物質の空間分布ｆを与えた場合に生体表面で検出される蛍光分布の理論値を得るための行列である。具体的には、システム行列の列ベクトルｖは、任意の位置における生体表面で検出される蛍光分布の理論値であり、この列ベクトルをｋ番目の位置までそれぞれ求めたのがシステム行列である。
ｇ＝Ａｆ          （３）
Ａ＝［ｖ（１），ｖ（２），…，ｖ（ｋ）］        （４）
ここで、
Ａ：システム行列、
ｆ：蛍光物質の空間分布ベクトル、
ｇ：生体表面での蛍光分布ベクトル、
を表す。例えば、蛍光物質の空間分布ｆ＝［０  ０  ０.５  １  ０.５  ０  …  ０］Tが与えられた場合、このベクトルにシステム行列を掛けると（Ａｆ）、生体表面での蛍光分布ベクトルが得られる。［Ａ  Ｂ  Ｃ］Tは転置行列である。
【００３９】
逆問題解析での演算は、装置で得られた蛍光検出データと順問題解析で計算した理論値から蛍光物質の空間分布を求める過程から成り立つ。まず、順問題解析によって（３）式におけるシステム行列Ａを求めることができる。一方、蛍光生体画像取得装置によって（４）式における生体表面での蛍光分布ベクトルgを得ることができる。逆問題解析では、順問題解析で求められたシステム行列Ａと生体表面での蛍光分布ベクトルgから、未知の蛍光物質の空間分布ベクトルfを求める。ベクトルfの求め方の一般的な方法は最小二乗法であり、以下の（５）式の評価関数を最小化することで求めることができる。
2015037055.tif 000004
【００４０】
一方、今回対象とする問題では、光の吸収・散乱によって測定ノイズの影響が強くなるため、（５）式の最小化問題では解が発散する可能性が高い。そこで、ノイズによる発散を防ぐ目的で以下の(６)式の最小化問題を考える。この最小化問題（チホノフ正則化）は逆問題の研究分野で一般的に利用されており、ベクトルfのノルムを評価関数に組み入れることで解の発散を防ぐことができる。
2015037055.tif 000005
ここで、
λ：正則化パラメータ、
を表す。このようにして、生体内部情報を求めることができ、生体内部での蛍光物質の画像を再構成することができる。
【００４１】
（手順３）（図２（３））
手順２によって蛍光物質のある場所ａの生体内部での位置情報が得られるので、遮断制御部２４はこの情報に基づいて場所ａに対して蛍光を物理的に遮断する遮断機構２２を作動させて場所ａからの蛍光を除去する（図３Ｂ）。場所ａからの蛍光が二次元検出器２０に入射しないように遮断するには、励起光が場所ａの蛍光物質に入射しないように遮断機構２２で励起光の一部を遮断するか、又は仮に励起光が場所ａの蛍光物質に入射しても場所ａの蛍光物質から発生する蛍光が二次元検出器２０に入射しないように遮断機構２２で蛍光を遮断する。
【００４２】
（手順４）（図２（４））
手順３によって、場所ａからの蛍光が二次元検出器２０に実質的に入射しないように限定された観測範囲に対して、再度蛍光画像を取得する（図３Ｂ）。このとき、場所ｂからの蛍光が十分検出できるよう、二次元検出器２０の取得条件を変えて撮像を行う。
【００４３】
さらに、図３Ｂに記号ｃで示されるようなさらに微弱な蛍光が検出される場合に、その蛍光画像も検出するときも手順２〜４を繰り返す。さらに微弱な蛍光を検出するときも同様である。
【００４４】
（手順５）（図２（５））
得られた複数の蛍光画像データに対して画像データを統合する。画像データの統合の際、二次元検出器２０の取得条件を揃えた上で画像データを統合する。統合する画像データは、二次元検出器２０が取得した蛍光画像データであってもよく、手順２で再構成された蛍光画像データであってもよい。
【００４５】
（手順６）（図２（６））
手順５で得たデータを、表示部２８に表示する。観察時にダイナミックレンジを拡げて測定することができたとしても、表示の際に表示部２８のダイナミックレンジが不足する可能性もある。そのため、表示については統合画像に対して、（ａ）表示範囲を任意に設定した上で表示する、（ｂ）対数表示する、（ｃ）蛍光強度をいくつかの範囲に分割して異なる色彩で表示する、といった方法で実施する。
【００４６】
次に、蛍光の遮断機構２２の具体例を説明する。
【００４７】
図４は、マウスのような小さい生体からなる試料の周囲に設置する遮断機構２２Ａの例である。直径が７０ｍｍ程度の円筒型遮断板４２Ａ，４２Ｂが生体試料１２の長手方向の前後の両側２ヶ所に設置されている。生体試料１２の長手方向とは、生体の頭部から尾部を結ぶ直線の方向である。遮断板４２Ａ，４２Ｂは生体試料１２の長手方向に沿って移動できるように、それぞれの電動スライダ４４Ａ，４４Ｂに支持されている。遮断板４２Ａ，４２Ｂの材質はアルミニウムやプラスチック等の軽量なものが好ましく、表面が黒色に塗装されていて光反射をなるべく抑えたものであることが好ましい。遮断板４２Ａ，４２Ｂは、遮断制御部２４からの指示により電動スライダ４４Ａ，４４Ｂによって駆動され、任意の領域まで移動して光遮断を行う。
【００４８】
図５は、励起光源１４及び二次元蛍光検出器２０と、試料１２の間に配置される遮断機構２２Ｂの例である。この例では、二次元蛍光検出器２０の手前に２つの励起光源１４が配置されている。励起光源１４の数は１つでもよく、３つ以上でもよい。励起光源１４は試料（図示されていない）に対して広い角度にわたって励起光１６を照射する。遮断機構２２Ｂは励起光源１４及び二次元蛍光検出器２０の設置位置と試料との間に配置されている。遮断機構２２Ｂは励起光源１４からの励起光１６の一部が試料に照射されるのを遮断するとともに、試料からの蛍光の一部が二次元蛍光検出器２０に入射するのを遮断する。
【００４９】
遮断機構２２Ｂは同一平面内に配置された４枚の遮断板４６Ａ〜４６Ｄを備え、それらの遮断板４６Ａ〜４６Ｄはそれぞれの電動スライダ４８Ａ〜４８Ｄによって互いに独立して移動するように支持されている。遮断板４６Ａ〜４６Ｄの材質はアルミニウムやプラスチック等の軽量なものが好ましく、表面が黒色に塗装されていて光反射をなるべく抑えたものであることが好ましい。遮断板４６Ａ〜４６Ｄは、遮断制御部２４からの指示により電動スライダ４８Ａ〜４８Ｄによって駆動され、任意の領域まで移動して光遮断を行う。
【００５０】
蛍光物質の生体内部での位置情報に対応して、励起・蛍光を物理的に遮断する遮断機構を動作させることが好ましい。その際、例えば、二次元検出器２０からみて生体内部の浅い位置に最も強い蛍光を発する蛍光物質がある場合は、蛍光が生体内部を透過する際の散乱による影響が小さいため、その蛍光を除去するための遮断範囲を狭めることができる。一方、二次元検出器２０からみて生体内部の深い位置にその蛍光物質がある場合は、蛍光が生体内部を透過する際の散乱による影響が大きいため、その蛍光を除去するための遮断範囲を拡げる必要がある。このようにして、生体計測において問題となる散乱問題を考慮しつつ、最適な撮像条件を与えることができる。
【００５１】
遮断機構２２によって遮断する範囲は、例えば以下のように決定する。
【００５２】
まず、生体内部の蛍光物質から発せられる蛍光が生体試料表面でどの程度広がって検出されるかのシミュレーションを実施する。例を図６に示す。図６（Ａ）に示すように、直径１ｍｍの光源を生体内部（吸収係数０．０２［／ｍｍ］、散乱係数１０［／ｍｍ］）に設置するものとする。図６（Ｂ）は、その設置する深さを異ならせた場合に、生体試料表面で検出される蛍光強度のプロファイルを示している。このプロファイルは光拡散方程式（１）、（２）によって計算することができる。
【００５３】
図６（Ｂ）から、距離２０ｍｍでの蛍光強度が１となるように規格化した蛍光強度のプロファイルを図７に示す。図７より、生体試料内部の蛍光がどの程度散乱されて表面で拡がって検出されるかについての指標を得ることができる。図７から得た情報に基づき、蛍光を遮断すべき範囲を決定する。例えば、蛍光物質が二次元検出器２０からみて生体試料の深さ２ｍｍの場所に存在した場合、蛍光強度のピーク地点から３ｍｍ程度の範囲を遮断すれば、この蛍光物質由来の蛍光を約１／５に減衰できる。この遮断によって、遮断部分以外の領域のダイナミックレンジを約５倍に高めることができるようになり、微弱光であっても検出できるようになる。１／５の減衰は実質的な遮断の一例を示す。
【００５４】
図７から得た知見に基づき、位置情報と遮断範囲を対応付けした表を予め作成してコンピュータ３４に保持しておけば、遮断制御部２４は遮断範囲を効率よく指定できる。また、肝臓などの各臓器や皮下の血管やリンパ節など、蛍光物質の集積箇所が予め判明している場合については、生体試料内での蛍光物質の位置情報が予めわかっているので、図７のようなシミュレーションを一度実施しておけば、予め遮断する範囲を入力部３２からコンピュータ３４に入力して保持することができる。
【００５５】
以上の遮断機構の動作の一例を示す。遮断機構２２が図４に示される遮断機構２２Ａである場合について図８により説明する。この遮断機構２２Ａは試料がマウスのような小さい生体が測定対象の場合に適する。コンピュータ３４は深さに対する遮断範囲の対応表を保持しておく。その対応表の一例は表１に示すものである。
【００５６】
【表１】
2015037055.tif 000006
【００５７】
この対応表は、遮断により二次元検出器２０に入射する蛍光の光減衰が１／５になるように予め計算によって求められている。二次元検出器２０に入射する蛍光の光減衰の度合いは、その蛍光の次に強い強度をもつ蛍光が検出できるように実質的に遮断するための度合いであり、必要に応じて変えることができる。
【００５８】
遮断機構２２Ａを使用しないで、観察対象の全領域からの蛍光画像を検出する１回目の撮像はすでに行われているものとする。
【００５９】
２回目の蛍光画像撮像にあたり、遮断制御部２４は、図８に示されるように、三次元分布計算処理部３０からの構成処理（図２（２））による蛍光物質の深さ情報を入力する。表１のデータを読み出し、深さ情報を基にしてから遮断範囲を決定（入力）する。遮断制御部２４はその遮断範囲に基づいて電動スライダ２３Ａ，２３Ｂの動作を制御して遮断板４２Ａ，４２Ｂを所定の位置まで移動させる。その結果、全画像範囲で最も蛍光信号が強かった箇所箇所からの蛍光の約８０％が遮断されて、２回目の蛍光画像撮像時（図２（４））では微弱な蛍光を検出することができるようになる。３回目以降の蛍光画像撮像は、それまでの蛍光画像撮像のそれぞれで最も強度の強かった蛍光画像の全てを実質的に遮断した状態で行う。
【００６０】
遮断機構２２が図５に示される遮断機構２２Ｂの動作も同様に制御することができる。
【００６１】
リンパ節、血管又は心臓というような生体臓器が検出対象の場合は図５に示される遮断機構２２Ｂが適する。コンピュータ３４は深さに対する遮断範囲の対応表を保持しておく。対応表は例えば表２で示されるものであり、リンパ節、血管及び心臓のそれぞれの位置も、二次元検出器からみて試料表面から何ｍｍの深さにあるかが予めわかっているので、それぞれの深さに対して遮断範囲をいくらにするかが予め求められている。この対応表も、遮断により二次元検出器２０に入射する蛍光の光減衰が１／５になるように予め計算によって求められたものである。
【００６２】
【表２】
2015037055.tif 000007
【００６３】
遮断機構２２Ｂを使用しないで、観察対象の全領域からの蛍光画像を検出する１回目の撮像はすでに行われているものとする。
【００６４】
２回目の蛍光画像撮像にあたり、１回目の撮像で最も蛍光強度の強かった臓器が何であるかがわかるので、それを遮断対象とする。遮断制御部２４は、図９に示されるように、入力部３２から遮断対象が何であるかを入力する。遮断制御部２４は表２のデータを読み出し、遮断対象に対する遮断範囲を決定（入力）する。遮断制御部２４はその遮断範囲に基づいて電動スライダ４８Ａ〜４８Ｄの動作を制御して遮断板４６Ａ〜４６Ｄを所定の位置まで移動させる。その結果、全画像範囲で最も蛍光信号が強かった臓器からの蛍光が実質的に遮断されて、２回目の蛍光画像撮像時（図２（４））では遮断対象の臓器以外の臓器からの微弱な蛍光を検出することができるようになる。３回目の蛍光画像撮像は、それまでの蛍光画像撮像のそれぞれで最も強度の強かった蛍光画像の全てを実質的に遮断した状態で行う。
【００６５】
（検証結果）
本発明を円筒型ファントムで検証した結果を以下に示す。
【００６６】
図１０（Ａ）に示されるように、直径２５ｍｍの円筒型ファントム５０（吸収係数０．０２［／ｍｍ］、散乱係数１０［／ｍｍ］）に、濃度が１桁異なるＩＣＧ（インドシアニン・グリーン）入りの２つのチューブ５２Ａと５２Ｂを設置したものを試料として用意した。両チューブ５２Ａ，５２Ｂはファントム５０の円柱の中心部ではなく、ある周面から６ｍｍの深さの場所とし、両チューブ５２Ａ，５２Ｂ間の距離を１０ｍｍ程度とした。ＩＣＧ濃度はチューブ５２Ａが５ｐｐｍ、チューブ５２Ｂが０．５ｐｐｍである。
【００６７】
この条件で観察対象全域の蛍光画像を取得すると、図１０（Ｂ）に示されるように、チューブ５２Ａ由来の蛍光が観察された。チューブ５２Ｂ由来の蛍光も微弱ではあるが観察された。図１０（Ｂ）の数値は蛍光検出強度を表わし、例えば「５ｅ＋２」は「５×１０2」を意味する。図１２、図１４においても同様である。
【００６８】
次に、図１０（Ｂ）の蛍光画像データを利用して、内部情報を計算すると、図１１のように二次元検出器側の周面から深さ６〜７ｍｍの場所にチューブ５２Ａの蛍光画像が再構成された。なお、今回の計算では、１桁程度しか濃度差がないということもあって、チューブ５２Ｂの蛍光画像ついても１／１０程度の濃度で再構成された。
【００６９】
次に、図１１の生体内部情報を元に、予め作成しておいた位置情報と遮断範囲を対応付けした表から蛍光の遮断範囲を算出した。この対応表として表１を使用した。チューブ５２Ａからの蛍光画像は二次元検出器側のファントム周面から深さ６〜７ｍｍの場所に再構成されたので、深さを７ｍｍに該当する遮断範囲として、表１のデータから５．８ｍｍを採用した。
【００７０】
次に、表１から得た情報を基に、図１２（Ａ）のように遮断機構２２Ｃを動作させてチューブ５２Ａからの蛍光が二次元検出器に入射するのを実質的に遮断し、その状態で再度蛍光画像を取得したのが図１２（Ｂ）である。図１０（Ｂ）の結果結果ではチューブ５２Ａからの強い蛍光に阻害されてチューブ５２Ｂからの蛍光画像を鮮明に検出することが難しかった。しかし、チューブ５２Ａからの蛍光を実質的に遮断したことによって二次元検出器の検出感度を上げることができ、図１２（Ｂ）に示されるように、チューブ５２Ｂからの微弱蛍光がコントラストよく検出された（図１２（Ｂ）中の矢印）。
【００７１】
次に、図１２（Ｂ）の蛍光画像データを利用して、内部情報を計算すると、図１３のように二次元検出器側のファントム周面から深さ６〜７ｍｍの場所にチューブ５２Ｂからの蛍光画像が再構成された。
【００７２】
最後に、図１０〜１３のように得た蛍光観察結果および生体内部情報の計算結果を、一つの画像上に統合して表示した。
【００７３】
図１４は、図１０、１２のように得られた蛍光画像を異なる２色を用いて同じ画面上に表示したものである。表示にあたり、ファントム表面での蛍光画像は散乱の影響で拡がっているため、コントラスト強調（ある閾値以下の強度はカット）を実施した上で画像を表示した。
【００７４】
図１５は、図１１、１３のように得られた再構成された蛍光結果を異なる２色を用いて同じ画面上に表示したものである。再構成された蛍光結果を用いることにより、励起光などによるノイズを除去することができる。
【符号の説明】
【００７５】
１０      試料台
１２      試料
１４      励起光源
１６      励起光
１８      撮像部
２０      二次元検出器
２２，２２Ａ，２２Ｂ，２２Ｃ      遮断機構
２４      遮断制御部
２６      画像統合部
２８      表示装置
３０      三次元分布計算処理部
３４      コンピュータ
【図１】
2015037055.tif 000008
【図２】
2015037055.tif 000009
【図７】
2015037055.tif 000010
【図８】
2015037055.tif 000011
【図９】
2015037055.tif 000012
【図３】
2015037055.tif 000013
【図４】
2015037055.tif 000014
【図５】
2015037055.tif 000015
【図６】
2015037055.tif 000016
【図１０】
2015037055.tif 000017
【図１１】
2015037055.tif 000018
【図１２】
2015037055.tif 000019
【図１３】
2015037055.tif 000020
【図１４】
2015037055.tif 000021
【図１５】
2015037055.tif 000022
【国際調査報告】
2015037055.tif 000023
2015037055.tif 000024
2015037055.tif 000025
2015037055.tif 000026