(57)【要約】【課題】レーザドップラ法による非接触式で、精度よく測定対象内の流体の流量を測定できるようにする。【解決手段】測定対象にレーザ光を照射して、前記測定対象内に所定の測定体積となる光収束部分を形成する手段と、前記光収束部分からの後方散乱光を受光する手段と、前記受光した信号を周波数解析する手段と、前記周波数解析の結果から、低周波数帯の信号強度に基づいて前記光収束部分の深さ方向の位置を制御する手段と、前記低周波数帯の信号強度と前記深さ方向の位置の関係から、前記周波数解析の結果の高周波数帯の信号強度に補正を行う手段と、前記補正が行われた高周波数帯の信号強度および周波数に基づいて流量を計算する手段とを備える。【選択図】図３
2018100968.tif 000002
【特許請求の範囲】
【請求項１】
測定対象にレーザ光を照射して、前記測定対象内に所定の測定体積となる光収束部分を形成する手段と、
前記光収束部分からの後方散乱光を受光する手段と、
前記受光した信号を周波数解析する手段と、
前記周波数解析の結果から、低周波数帯の信号強度に基づいて前記光収束部分の深さ方向の位置を制御する手段と、
前記低周波数帯の信号強度と前記深さ方向の位置の関係から、前記周波数解析の結果の高周波数帯の信号強度に補正を行う手段と、
前記補正が行われた高周波数帯の信号強度および周波数に基づいて流量を計算する手段と
を備えたことを特徴とする流量測定装置。
【請求項２】
前記測定対象内における光収束部分の深さ方向の位置を所定の範囲で移動して、低周波数帯の信号強度と深さ方向の位置の関係を示す情報を取得する手段を備え、
前記光収束部分の深さ方向の位置を制御する手段は、ユーザから指定された深さ方向の位置に基づき、前記低周波数帯の信号強度と深さ方向の位置の関係を示す情報から対応する低周波数帯の信号強度を取得し、該信号強度に測定時の低周波数帯の信号強度が一致するように制御を行う
ことを特徴とする請求項１に記載の流量測定装置。
【請求項３】
前記測定対象内における光収束部分の深さ方向の位置を所定の範囲で移動して、低周波数帯の信号強度と深さ方向の位置の関係を示す情報を取得する手段を備え、
前記光収束部分の深さ方向の位置を制御する手段は、前記低周波数帯の信号強度と深さ方向の位置の関係を示す情報から、所定の刻みで決定した深さ方向の位置に対応する低周波数帯の信号強度を取得し、該信号強度に測定時の低周波数帯の信号強度が一致するように制御を行う
ことを特徴とする請求項１に記載の流量測定装置。
【請求項４】
測定対象にレーザ光を照射して、前記測定対象内に所定の測定体積となる光収束部分を形成し、
前記光収束部分からの後方散乱光を受光し、
前記受光した信号を周波数解析し、
前記周波数解析の結果から、低周波数帯の信号強度に基づいて前記光収束部分の深さ方向の位置を制御し、
前記低周波数帯の信号強度と前記深さ方向の位置の関係から、前記周波数解析の結果の高周波数帯の信号強度に補正を行い、
前記補正が行われた高周波数帯の信号強度および周波数に基づいて流量を計算する、
処理を流量測定装置が実行することを特徴とする流量測定方法。
【請求項５】
測定対象にレーザ光を照射して、前記測定対象内に所定の測定体積となる光収束部分を形成し、
前記光収束部分からの後方散乱光を受光し、
前記受光した信号を周波数解析し、
前記周波数解析の結果から、低周波数帯の信号強度に基づいて前記光収束部分の深さ方向の位置を制御し、
前記低周波数帯の信号強度と前記深さ方向の位置の関係から、前記周波数解析の結果の高周波数帯の信号強度に補正を行い、
前記補正が行われた高周波数帯の信号強度および周波数に基づいて流量を計算する、
処理を流量測定装置に実行させることを特徴とする流量測定プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、流量測定装置、流量測定方法および流量測定プログラムに関する。
【背景技術】
【０００２】
臨床医療の分野における新たな治療指標として、末梢血管の血流量値（perfusion）が注目されている。末梢血管の血流量値は、例えば、狭心症・心筋梗塞といった循環器系疾患、精神的ストレス、腫瘍の進行度、火傷の損傷度等と相関があり、定量的な診断指標として応用が期待されている。
【０００３】
患者に負担を与えない観点から非侵襲の測定手法が望まれるところ、非侵襲に生体の血流を測定する手法としてレーザドップラ法（Laser Doppler Flowmetry：ＬＤＦ）が知られている（非特許文献１、２等を参照）。
【０００４】
図１はレーザドップラ法による血流測定の原理を示す図である。図１において、生体の測定対象１にレーザ光を照射すると、生体内で多重反射され、後方散乱光が生じる。後方散乱光のうち、静止組織によって散乱された光は照射光と同じ周波数を保つが、血管１１中の赤血球１２で散乱された光は、赤血球１２の移動速度に応じてドップラシフトした周波数となる。ドップラシフトしていない光とドップラシフトした光は干渉し、様々なビート周波数をもった信号を検知することができる。ここで、ビート周波数は赤血球の速度に対応し、各ビート周波数の信号強度は所定の速度の赤血球の量に対応するため、検知した信号を周波数解析（高速フーリエ変換等）した結果から、式（１）により血流量値Qpを計算することができる。
【０００５】
【数１】
2018100968.tif 000003
【０００６】
式（１）において、ωは信号の角周波数、P(ω)はパワースペクトル、iは低周波数帯（ほぼ直流に近い周波数帯）の信号強度である。i2で割っているのは、積分で求めた値は測定の状態によって左右され、そのままでは比較が困難であるため、値を規格化するためである。
【０００７】
また、レーザドップラ法による測定装置を実装する上では、接触式と非接触式とがある。図２（ａ）は接触式によるプローブ５ａを示しており、プローブ５ａの先端を測定対象１に接触させ、光の照射と受光を行う。図２（ｂ）は非接触式によるプローブ５ｂを示しており、プローブ５ｂの先端を測定対象１に対して非接触に保った状態で光の照射と受光を行う。接触式は、プローブによる接触圧力やプローブの接触による温度変化により血流量に影響を与えてしまうとともに、火傷の損傷部には使用できない。その点で、非接触式が望ましい。
【先行技術文献】
【非特許文献】
【０００８】
【非特許文献１】M. Karlsson, K. Wardell, "Polarized laser Doppler perfusion imaging-reduction of movement-induced artifacts", Journal of Biomedical Optics, 10, 6, (2005), 1-9
【非特許文献２】「体動下で測定可能な小型血流量装置の開発  予備的実験による測定システムの妥当性検証」、橋本将明、鎌田慎、田口良広、第５３回  日本伝熱シンポジウム講演論文集（2016-5）
【発明の概要】
【発明が解決しようとする課題】
【０００９】
上述したように、レーザドップラ法による非接触式の血流測定が望まれるが、体動等に起因するアーチファクト（偽の異常所見）により正確な測定が行えないという問題があった。この点、非特許文献１では、偏光板をディテクタ（検出器）の前に設置することで対処しようとしているが、充分に体動アーチファクトを低減することができないとともに、偏光板による信号光の減衰によりＳＮ比の低下を招くものであった。また、非特許文献２では、体動アーチファクトの低減を目的とした提言がなされているが、精度を高めるための具体的な手法は開示されていない。
【００１０】
なお、血流量測定について説明したが、レーザドップラ法により測定可能な他の流体の流量測定においても同様の問題が生ずる可能性がある。
【００１１】
本発明は上記の従来の問題点に鑑み提案されたものであり、その目的とするところは、レーザドップラ法による非接触式で、精度よく測定対象内の流体の流量を測定できるようにすることにある。
【課題を解決するための手段】
【００１２】
上記の課題を解決するため、本発明にあっては、測定対象にレーザ光を照射して、前記測定対象内に所定の測定体積となる光収束部分を形成する手段と、前記光収束部分からの後方散乱光を受光する手段と、前記受光した信号を周波数解析する手段と、前記周波数解析の結果から、低周波数帯の信号強度に基づいて前記光収束部分の深さ方向の位置を制御する手段と、前記低周波数帯の信号強度と前記深さ方向の位置の関係から、前記周波数解析の結果の高周波数帯の信号強度に補正を行う手段と、前記補正が行われた高周波数帯の信号強度および周波数に基づいて流量を計算する手段とを備える。
【発明の効果】
【００１３】
本発明にあっては、レーザドップラ法による非接触式で、精度よく測定対象内の流体の流量を測定することができる。
【図面の簡単な説明】
【００１４】
【図１】レーザドップラ法による血流測定の原理を示す図である。
【図２】接触式と非接触式の説明図である。
【図３】システム構成例を示す図（その１）である。
【図４】システム構成例を示す図（その２）である。
【図５】システム構成例を示す図（その３）である。
【図６】システム構成例を示す図（その４）である。
【図７】測定体積の例を示す図である。
【図８】空間分解能の説明図である。
【図９】コントローラの情報処理にかかるハードウェア構成例を示す図である。
【図１０】レンズ制御の処理例を示すフローチャート（その１）である。
【図１１】周波数解析結果の例を示す図である。
【図１２】処理過程のデータ等の例を示す図である。
【図１３】血流量測定の処理例を示すフローチャートである。
【図１４】レンズ制御の処理例を示すフローチャート（その２）である。
【図１５】レンズ制御の処理例を示すフローチャート（その３）である。
【図１６】実験構成例を示す図である。
【図１７】実験結果を示す図である。
【発明を実施するための形態】
【００１５】
以下、本発明の好適な実施形態につき説明する。なお、血流量測定について説明するが、レーザドップラ法により測定可能な他の流体の流量測定に適用できることはいうまでもない。
【００１６】
＜構成＞
図３〜図６は実施形態にかかるシステム構成例を示す図であり、図３および図４はセンサヘッド２がチップ型である場合の例、図５および図６はセンサヘッド２が内視鏡型である場合の例を示している。
【００１７】
図３において、センサヘッド２は、いわゆる共焦点光学系を構成するものであり、レーザダイオード等の光源２０１の光路に沿って配置されたピンホール２０２、レンズ２０３、ビームスプリッタ２０４、レンズ２０５を備えている。レンズ２０３からレンズ２０５までの光路は平行光となるように配置されている。また、レンズ２０５には、このレンズ２０５を光軸方向に移動（駆動）可能なＭＥＭＳ（Micro Electro Mechanical Systems）等によるレンズ駆動機構２０６が設けられている。レンズ２０５からは測定対象１に対して収束されたレーザ光が照射されるとともに、レンズ２０５には測定対象１からの後方散乱光が入射される。また、ビームスプリッタ２０４により直角に反射された後方散乱光の光路には、レンズ２０７、ピンホール２０８、検知器（光検知器）２０９が配置されている。
【００１８】
なお、レンズ２０３、２０５、２０７の開口数（Numerical Aperture：ＮＡ）およびピンホール２０２、２０８の径により、測定対象１側で光が収束する部分の大きさ（測定体積）が所定の値となるように設計および調整が行われている。図７（ａ）（ｂ）は、レンズ２０５の光路前方の光が収束する部分の測定体積ＳＶを示している。測定体積ＳＶを小さく、かつ安定した大きさとすることで、図７（ｃ）に示すように、測定対象１の表面からの深さＤ１、Ｄ２、Ｄ３等における分解能が向上するとともに、その後の信号の補正等が容易になる。すなわち、測定対象１内の生体からは、光が収束する測定体積ＳＶ内からの後方散乱光がほとんどを占めることになり、測定体積ＳＶが小さくなることで細かい範囲内の血流量を他と区別して測定することができる。また、測定体積ＳＶの値が一定に保たれることで、測定体積ＳＶの変化に起因する測定値の誤差を考慮しなくて済み、補正が不要となる。なお、測定部分が生体の内部に入るにつれて（深さ位置が測定対象１の内部に入るにつれて）、生体組織による光の吸収量が増え、光が収束する測定体積ＳＶ内からの後方散乱光の減衰が大きくなることから、深さ位置に対する測定値の補正は必要となる。この点については、後述する。
【００１９】
測定体積について、空間分解能との関係で、より詳しく説明する。測定体積は一般的に空間分解能で表現される。図８に示すように、一般に、ｚ軸方向の空間分解能は、平面反射板を光軸方向に移動させた時のデフォーカス量に対するディテクタ検知信号強度分布の半値全幅d(3dB)で定義される。適当なピンホール半径aを選択した時の半値全幅dは、nを屈折率、Mを光学倍率、N.A.をレンズ開口数として、式（２）で表される。
【００２０】
【数２】
2018100968.tif 000004
【００２１】
本実施形態での測定対象は、反射体ではなく、散乱体であるため、上記の式（２）を直接に適用することはできないが、簡易的に式（３）で表すことができる。C0は測定対象の光学特性と装置の仕様に依存する定数である。
【００２２】
【数３】
2018100968.tif 000005
【００２３】
図３に戻り、光源２０１、レンズ駆動機構２０６、検知器２０９はそれぞれ一対の電線３１、３２、３３を束ねたケーブル３を介してコントローラ４に接続される。
【００２４】
コントローラ４は、光源駆動ドライバ４０１とレンズ駆動ドライバ４０２とＡ／Ｄ変換器４０３と周波数解析部４０４とレンズ制御部４０５と入力部４０６とスペクトル補正部４０７と血流量計算部４０８と表示部４０９とを備えている。
【００２５】
光源駆動ドライバ４０１は、電線３１を介して接続されたセンサヘッド２の光源２０１を駆動（駆動電源供給）する機能を有している。レンズ駆動ドライバ４０２は、レンズ制御部４０５による制御のもと、電線３２を介して接続されたセンサヘッド２のレンズ駆動機構２０６を駆動する機能を有している。Ａ／Ｄ変換器４０３は、電線３３を介してセンサヘッド２の検知器２０９から取得した受光信号（アナログ値）をデジタル値に変換する機能を有している。
【００２６】
周波数解析部４０４は、Ａ／Ｄ変換器４０３によりデジタル化された受光信号をＦＦＴ（Fast Fourier Transform）やＬＰＦ（Low-pass filter）等により周波数解析（周波数毎の信号強度であるパワースペクトルの取得）を行う機能を有している。レンズ制御部４０５は、周波数解析部４０４による周波数解析結果のうち低周波数帯の信号強度を用い、レンズ駆動ドライバ４０２およびセンサヘッド２のレンズ駆動機構２０６を介してレンズ２０５を所定の深さ方向の位置に制御する機能を有している。制御に低周波数帯の信号強度を用いるのは、高周波数帯の信号強度が赤血球の速さに応じてドップラシフトした信号に基づくものであるのに対し、低周波数帯の信号強度は測定対象１の表面反射や生体内の静止組織による散乱光に基づくものであり、深さ方向と相関があるためである。入力部４０６は、ユーザから測定開始の操作指示や血流量を測定する測定対象１の表面からの深さ位置等の入力（指定の入力）を行う機能を有している。
【００２７】
スペクトル補正部４０７は、血流量の測定に際して、測定が行われる深さ位置による生体内での減衰等に対するパワースペクトルの補正を行う機能を有している。血流量計算部４０８は、スペクトル補正部４０７により補正が行われたパワースペクトルP(ω)から血流量を計算する機能を有している。表示部４０９は、血流量の測定結果をユーザに対して表示する機能を有している。なお、血流量の測定結果は、測定が行われた深さ位置や測定日時等とともに、データとして蓄積・出力することもできる。
【００２８】
図４は、コントローラ４において、新たにアナログ周波数解析部４１０を設け、レンズ制御部４０５に代えてアナログレンズ制御部４１１を設けることで、デジタル値ではなくアナログ値を用いてレンズ制御を行うようにしたものである。デジタル値の分解能をそれほど高くできない場合やデジタル処理による時間遅れが制御の精度に影響する場合には、アナログ値による方が精度よくレンズ制御を行うことができる。なお、この変形は、後述する図５および図６のコントローラ４にも適用することができる。
【００２９】
図５は、内視鏡型のセンサヘッド２に適用しやすいように、図３のセンサヘッド２から光源２０１と検知器２０９を除去し、コントローラ４内に光源４１２と検知器４１３を配置することで小型化したものである。また、光源４１２および検知器４１３とセンサヘッド２の間は光ファイバ３４、３５により接続され、その端部をレンズ２０３、２０７の光軸に配置することで、ピンホール２０２、２０８も不要にして、いっそうの小型化を図っている。
【００３０】
図６は、図５のセンサヘッド２からビームスプリッタ２０４とレンズ２０７を除去し、光ファイバ３４に往路と復路の光を重畳させたものである。これに対応し、コントローラ４では光ファイバ３４の端部に光サーキュレータ４１４を配置し、往路は光源４１２から光ファイバ３４に光が流れ、復路は光ファイバ３４から検知器４１３に光が流れるようにしている。
【００３１】
図９はコントローラ４の情報処理にかかるハードウェア構成例を示す図である。情報処理としては、主に、周波数解析部４０４、レンズ制御部４０５、スペクトル補正部４０７、血流量計算部４０８における処理と、その他の各部の制御の処理が含まれる。
【００３２】
図９において、コントローラ４は、バス４０を介して相互に接続されたＣＰＵ（Central Processing Unit）４１、ＲＯＭ（Read Only Memory）４２、ＲＡＭ（Random Access Memory）４３、ＨＤＤ（Hard Disk Drive）／ＳＳＤ（Solid State Drive）４４、接続Ｉ／Ｆ（Interface）４５を備えている。ＣＰＵ４１は、ＲＡＭ４３をワークエリアとしてＲＯＭ４２またはＨＤＤ／ＳＳＤ４４等に格納されたプログラムを実行することで、コントローラ４の動作を統括的に制御する。プログラムは、記録媒体を経由して取得されるものでもよいし、ネットワークを経由して取得されるものでもよいし、ＲＯＭ組込でもよい。接続Ｉ／Ｆ４５は、ハードウェア機器とのインタフェースである。
【００３３】
＜動作＞
図１０は、図３〜図６に示したシステム構成によるレンズ制御の処理例を示すフローチャートである。
【００３４】
図１０において、ユーザが測定対象１にセンサヘッド２を近づけると（ステップＳ１０１）、センサヘッド２では光源２０１（４１２）からのレーザ光が光学系を介して測定対象１に照射され、後方散乱光が検知器２０９、４１３により受光される。そして、周波数解析部４０４（アナログ周波数解析部４１０）により周波数解析して得られた低周波数帯の信号強度がレンズ制御部４０５（アナログレンズ制御部４１１）により取得される（ステップＳ１０２）。
【００３５】
図１１は周波数解析部４０４（アナログ周波数解析部４１０）による周波数解析結果の例を示す図であり、横軸は周波数、縦軸は信号強度を示している。図１１では同じ深さ位置における３つの流速の場合の結果を示しており、Ｆ１→Ｆ２→Ｆ３の順に流れが速くなる。すなわち、流れが速いほど高い周波数における信号強度が高くなる。周波数が「０」付近の低周波数帯の信号強度は流速によらずほぼ一定となっている。なお、低周波数帯の信号強度は、表面付近で最大（ピーク）となり、測定箇所が深いほど値が小さくなる。
【００３６】
図１０に戻り、次いで、レンズ制御部４０５（アナログレンズ制御部４１１）は、取得した信号強度の範囲内で目標信号強度を設定し（ステップＳ１０３）、周波数解析して得られた低周波数帯の信号強度が目標信号強度に等しくなるように、レンズ駆動ドライバ４０２およびレンズ駆動機構２０６を介して、レンズ２０５の光軸方向の位置を制御し（ステップＳ１０４）、深さ位置を一定に保つ。なお、レンズ２０５の位置を制御する過程で、図１２（ａ）に示すような、低周波数帯の信号強度のデータを取得して蓄積する。
【００３７】
図１０に戻り、次いで、レンズ制御部４０５（アナログレンズ制御部４１１）は、現在の目標信号強度がピーク強度であるか否か判断する（ステップＳ１０５）。ピーク強度とは、レンズ・測定対象間の距離に対する信号強度の最大値である。そして、現在の目標信号強度がピーク強度でないと判断された場合（ステップＳ１０５のＮｏ）、レンズ制御部４０５（アナログレンズ制御部４１１）は、ピーク強度値が得られるように所定量シフトし（ステップＳ１０６）、レンズ２０５の光軸方向の位置制御（ステップＳ１０４）に戻る。
【００３８】
また、現在の目標信号強度がピーク強度であると判断された場合（ステップＳ１０５のＹｅｓ）、レンズ制御部４０５（アナログレンズ制御部４１１）は、ピーク強度を用いて解析することで、信号強度とレンズ・測定対象間の距離の関係を取得する（ステップＳ１０７）。図１２（ｂ）は信号強度とレンズ・測定対象間の距離の関係を示したものであり、実線の部分が実際に測定されたデータの部分である。この関係から、例えば、ピーク時のレンズ位置を基準として、図１２（ｃ）に示すような、信号強度と深さ位置との関係を取得することができる。
【００３９】
図１０に戻り、次いで、レンズ制御部４０５（アナログレンズ制御部４１１）は、表面からの血流量測定箇所の深さ位置を決定する（ステップＳ１０８）。例えば、入力部４０６からユーザに深さ位置を指定させて決定するか、予め設定された深さ位置に設定する。
【００４０】
次いで、レンズ制御部４０５（アナログレンズ制御部４１１）は、先に取得した信号強度とレンズ・測定対象間の距離の関係（図１２（ｂ）（ｃ））から、目標信号強度を求めて設定する（ステップＳ１０９）。
【００４１】
次いで、レンズ制御部４０５（アナログレンズ制御部４１１）は、周波数解析して得られた低周波数帯の信号強度が目標信号強度に等しくなるように、レンズ駆動ドライバ４０２およびレンズ駆動機構２０６を介して、レンズ２０５の光軸方向の位置を制御し（ステップＳ１１０）、深さ位置を一定に保つ。
【００４２】
次いで、レンズ制御部４０５（アナログレンズ制御部４１１）は、スペクトル補正部４０７および血流量計算部４０８に血流量測定の開始を指示して血流量測定を開始させ（ステップＳ１１１）、その後に血流量測定が終了する（ステップＳ１１２）。
【００４３】
図１３は血流量測定の処理例を示すフローチャートである。図１３において、スペクトル補正部４０７は、周波数解析の結果から生のパワースペクトルP0(ω)を取得し（ステップＳ１２１）、深さ位置を考慮したスペクトル補正により補正後のパワースペクトルP(ω)を取得する（ステップＳ１２２）。例えば、図１２（ｃ）のデータに基づき、ピーク時の信号強度に現在の信号強度が合致するように、ピーク時の信号強度と現在の深さ位置の信号強度の比をパワースペクトルP0(ω)に乗算することでP(ω)を得る。
【００４４】
図１３に戻り、次いで、血流量計算部４０８は、式（４）により血流量値Qpを計算する（ステップＳ１２３）。
【００４５】
【数４】
2018100968.tif 000006
【００４６】
式（４）は、前述した従来の式（１）と比較して、i2で規格化していないが、これは、測定体積が測定対象の光学特性の影響を受けず、上述したスペクトル補正により安定した値として得られるからである。また、式（４）で得られた血流量値Qpは、物理的な流量値（単位時間当たりの質量等）に換算が可能である。
【００４７】
次いで、血流量計算部４０８は、測定した血流量を表示部４０９により表示する（ステップＳ１２４）。なお、血流量の測定結果は、測定が行われた深さ位置や測定日時等とともに、データとして蓄積・出力することもできる。
【００４８】
図１４は、図３〜図６に示したシステム構成によるレンズ制御の他の処理例を示すフローチャートであり、血流測定前にレンズ駆動機構２０６により移動可能な全範囲について事前にスキャンを行って全範囲についてのキャリブレーションを行うようにしたものである。
【００４９】
図１４において、ユーザは測定対象１にセンサヘッド２を近づけた状態で両者を固定する（ステップＳ２０１）。なお、ここでの固定は、体動等による微小な動きについてまで固定するということではなく、スキャンの精度を保つ程度の固定を意味する。
【００５０】
次いで、レンズ制御部４０５（アナログレンズ制御部４１１）は、レンズ駆動ドライバ４０２およびレンズ駆動機構２０６を介してレンズ２０５を駆動し、測定対象をスキャンする（ステップＳ２０２）。そして、レンズ制御部４０５（アナログレンズ制御部４１１）は、スキャンに際しての、信号強度とレンズ・測定対象間の距離の相関関係データを取得（キャリブレーション）する（ステップＳ２０３）。取得されるデータや関係は図１２に示したのと同様であるが、図１２（ｂ）において実際に取得されるデータの範囲（実線部分）はスキャンされる範囲まで拡大される。
【００５１】
図１４に戻り、次いで、レンズ制御部４０５（アナログレンズ制御部４１１）は、表面からの血流量測定箇所の深さ位置を決定する（ステップＳ２０４）。例えば、入力部４０６からユーザに深さ位置を指定させて決定するか、予め設定された深さ位置に設定する。
【００５２】
次いで、レンズ制御部４０５（アナログレンズ制御部４１１）は、先に取得した信号強度とレンズ・測定対象間の距離の関係（図１２（ｂ）（ｃ））から、目標信号強度を求めて設定する（ステップＳ２０５）。
【００５３】
次いで、レンズ制御部４０５（アナログレンズ制御部４１１）は、周波数解析して得られた低周波数帯の信号強度が目標信号強度に等しくなるように、レンズ駆動ドライバ４０２およびレンズ駆動機構２０６を介して、レンズ２０５の光軸方向の位置を制御し（ステップＳ２０６）、深さ位置を一定に保つ。
【００５４】
次いで、レンズ制御部４０５（アナログレンズ制御部４１１）は、スペクトル補正部４０７および血流量計算部４０８に血流測定の開始を指示して血流量測定を開始させ（ステップＳ２０７）、その後に血流量測定が終了する（ステップＳ２０８）。血流量測定の処理は、図１３に示したのと同様である。
【００５５】
図１５は、図３〜図６に示したシステム構成によるレンズ制御の他の処理例を示すフローチャートであり、キャリブレーションを行う点は図１４と同様であるが、測定対象１の固定を不要としたものである。
【００５６】
図１５において、レンズ制御部４０５（アナログレンズ制御部４１１）は、測定対象１にセンサヘッド２が近づいた状態で、体動、揺らぎによるレンズ・測定対象間の距離の時間変化が無視できるよう、高速にレンズ２０５を駆動し、対象をスキャンする（ステップＳ３０１）。
【００５７】
次いで、レンズ制御部４０５（アナログレンズ制御部４１１）は、スキャンに際しての、信号強度とレンズ・測定対象間の距離の相関関係データを取得（キャリブレーション）する（ステップＳ３０２）。取得されるデータや関係は図１２に示したのと同様であるが、図１２（ｂ）において実際に取得されるデータの範囲（実線部分）はスキャンされる範囲まで拡大される。
【００５８】
図１５に戻り、次いで、レンズ制御部４０５（アナログレンズ制御部４１１）は、表面からの血流量測定箇所の深さ位置を決定する（ステップＳ３０３）。例えば、入力部４０６からユーザに深さ位置を指定させて決定するか、予め設定された深さ位置に設定する。
【００５９】
次いで、レンズ制御部４０５（アナログレンズ制御部４１１）は、先に取得した信号強度とレンズ・測定対象間の距離の関係（図１２（ｂ）（ｃ））から、目標信号強度を求めて設定する（ステップＳ３０４）。
【００６０】
次いで、レンズ制御部４０５（アナログレンズ制御部４１１）は、周波数解析して得られた低周波数帯の信号強度が目標信号強度に等しくなるように、レンズ駆動ドライバ４０２およびレンズ駆動機構２０６を介して、レンズ２０５の光軸方向の位置を制御し（ステップＳ３０５）、深さ位置を一定に保つ。
【００６１】
次いで、レンズ制御部４０５（アナログレンズ制御部４１１）は、スペクトル補正部４０７および血流量計算部４０８に血流測定の開始を指示して血流量測定を開始させ（ステップＳ３０６）、その後に血流量測定が終了する（ステップＳ３０７）。血流量測定の処理は、図１３に示したのと同様である。
【００６２】
また、図１４および図１５の処理において、スキャンした範囲内において、所定の刻みで変化させて深さ位置においての血流量測定を自動的に行うこともできる。この場合、キャリブレーション（図１４のステップＳ２０３、図１５のステップＳ３０２）の後に、ピーク位置から下限の深さ位置までを所定の刻みで変化させて行き、その都度に血流量測定を行えばよい。
【００６３】
また、図１０、図１４および図１５の処理において、所定の刻みで測定深さを変化させて、深さ方向の血流量分布測定を自動的に行うこともできる。この場合、信号強度とレンズ・測定対象間の距離の関係を取得（図１０のステップＳ１０７、図１４のステップＳ２０３、図１５のステップＳ３０２）後に、ピーク位置から下限の深さ位置までを所定の刻みで変化させて行き、その都度に血流量測定を行えばよい。
【００６４】
また、上述の方法によって深さ方向の血流量分布から、任意の深さに位置する血管等が判断可能であり、任意の深さ位置に対応する信号強度にレンズを制御することで、選択的に特定箇所（血管等）を測定することが可能である。
【００６５】
＜総括＞
以上説明したように、本実施形態によれば、レーザドップラ法による非接触式で、精度よく測定対象内の流体の流量を測定することができる。すなわち、測定の深さ位置を一定に保った状態で血流量測定を行うため、体動等に影響されることがなくなり、体動等によるアーチファクトを低減できる。また、レーザ光の収束する測定体積を小さく、かつ安定した大きさとできることで、深さ位置に対する分解能を高められるとともに、その後の補正を不要とし、安定した測定値を得ることができる。
【００６６】
＜実験結果＞
図１６は、実験構成例を示す図である。以下、マイクロ流路１００を有する生体を模擬した試料を測定対象１とし、図示するような実験構成で行った実験結果を説明する。
【００６７】
図１７は、実験結果を示す図である。図１６に示すような実験条件で実験した結果を図１７に示す。なお、図１７では、縦軸値「Ｂｌｏｏｄ  ｆｌｏｗ  ｒａｔｅ」を最大値で正規化して示す。
【００６８】
まず、この実験では、図１６に示すように、測定対象１は、深さ方向（図では、「Ｓｃａｎｎｉｎｇ  ｄｉｒｅｃｔｉｏｎ」で示し、左方向である。）に流路幅が２００  μｍ（マイクロメートル）であるマイクロ流路１００を有する。
【００６９】
本実験では、マイクロ流路１００には、赤血球溶液を一定の流速で流した。その上で、本実験では、測定対象１に対して流速測定装置が走査し（図では、左方向となり、「Ｆｏｃａｌ  ｐｏｉｎｔ  ｐｏｓｉｔｉｏｎ」の方向である。）、各走査位置で流速測定装置がスペクトルを取得した。そして、取得されたそれぞれのスペクトルから、血流量値を算出した。
【００７０】
なお、本実験では、スキャンのステップ幅を５０  μｍとした。
【００７１】
血流量値は、測定対象となる体積中において、赤血球速度によって決定する。したがって、ファントムの走査によって深さ方向における分解能が評価できる。
【００７２】
図１７に示すように、本実験では、マイクロ流路１００の流路幅において血流量値が上昇する結果が得られた。具体的には、本実験では、まず、原点（図１６及び図１７における「Ｆｏｃａｌ  ｐｏｉｎｔ  ｐｏｓｉｔｉｏｎ」が「０  μｍ」の位置に相当する。）をファントム表面（図１６における「Ｓｕｒｆａｃｅ」の位置である。）とした。そして、図１６に示すように、測定対象１には、「Ｓｕｒｆａｃｅ」から深さ方向（図では、左方向である。）において「２００  μｍ」乃至「４００  μｍ」の位置にマイクロ流路１００がある。
【００７３】
これに対して、図１７に示す実験結果からわかるように、本実験では、特に「Ｆｏｃａｌ  ｐｏｉｎｔ  ｐｏｓｉｔｉｏｎ」が「２００  μｍ」乃至「４００  μｍ」の範囲、すなわち、マイクロ流路１００がある位置において、血流量値が上昇する結果が測定できた。この結果から、本実施形態における流速測定装置は、深さ方向における分解能を有するといえる。
【００７４】
以上、本発明の好適な実施の形態により本発明を説明した。ここでは特定の具体例を示して本発明を説明したが、特許請求の範囲に定義された本発明の広範な趣旨および範囲から逸脱することなく、これら具体例に様々な修正および変更を加えることができることは明らかである。すなわち、具体例の詳細および添付の図面により本発明が限定されるものと解釈してはならない。
【符号の説明】
【００７５】
１                    測定対象
２                    センサヘッド
２０１                光源
２０２                ピンホール
２０３                レンズ
２０４                ビームスプリッタ
２０５                レンズ
２０６                レンズ駆動機構
２０７                レンズ
２０８                ピンホール
２０９                検知器
３                    ケーブル
３１、３２、３３      電線
３４、３５            光ファイバ
４                    コントローラ
４０１                光源駆動ドライバ
４０２                レンズ駆動ドライバ
４０３                Ａ／Ｄ変換器
４０４                周波数解析部
４０５                レンズ制御部
４０６                入力部
４０７                スペクトル補正部
４０８                血流量計算部
４０９                表示部
４１０                アナログ周波数解析部
４１１                アナログレンズ制御部
４１２                光源
４１３                検知器
４１４                光サーキュレータ
【図１】
2018100968.tif 000007
【図２】
2018100968.tif 000008
【図３】
2018100968.tif 000009
【図４】
2018100968.tif 000010
【図５】
2018100968.tif 000011
【図６】
2018100968.tif 000012
【図７】
2018100968.tif 000013
【図８】
2018100968.tif 000014
【図９】
2018100968.tif 000015
【図１０】
2018100968.tif 000016
【図１１】
2018100968.tif 000017
【図１２】
2018100968.tif 000018
【図１３】
2018100968.tif 000019
【図１４】
2018100968.tif 000020
【図１５】
2018100968.tif 000021
【図１６】
2018100968.tif 000022
【図１７】
2018100968.tif 000023