JP4230736B2 - Diagnostic imaging equipment - Google Patents
Diagnostic imaging equipment Download PDFInfo
- Publication number
- JP4230736B2 JP4230736B2 JP2002239220A JP2002239220A JP4230736B2 JP 4230736 B2 JP4230736 B2 JP 4230736B2 JP 2002239220 A JP2002239220 A JP 2002239220A JP 2002239220 A JP2002239220 A JP 2002239220A JP 4230736 B2 JP4230736 B2 JP 4230736B2
- Authority
- JP
- Japan
- Prior art keywords
- light
- unit
- living body
- ultrasonic
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Description
【0001】
【発明が属する技術分野】
この発明は画像診断装置に関し、特に生体光計測装置と超音波診断装置とを組み合わせ、散乱体である生体内の吸収物質の変化量を画像化するための画像診断装置に関する。
【0002】
【従来の技術】
可視から近赤外領域の波長の光を用いて、生体の酸素代謝機能を簡便かつ非侵襲的に計測する診断装置として、生体光計測装置が開発、実用化されている。生体光計測装置では、生体内に存在するヘモグロビン等の特定の吸収物質が、酸素状態と脱酸素状態とで特定波長の光の吸収性が変化することを利用して、生体の酸素代謝機能を画像として表示するようにしたものである。さらに生体光計測装置を発展させたものとして、複数の光照射部と複数の光計測部とをマトリックス状に配置して、比較的広い範囲の生体機能情報を画像として表示できる光トポグラフィ装置がある。光トポグラフィ装置は、例えば脳の活性部位の特定やてんかん患者の発作時における焦点位置などの生体情報を得ることができ、診断装置としての重要性が増している。
【0003】
しかし、このような生体光計測装置は、光照射部から光測定部までの光路における光の減衰量(その変化)を測定するものであるため、吸収物質の絶対値を直接計測することはできない。即ち、光照射部と光測定部との間で照射された光が通る経路は、光照射部−光測定部間の距離のみならず生体の構造にも依存し、光路長が異なれば光の減衰量も異なることとなり、生体光計測で計測された値には、このような構造に依存する光路長のファクタも含まれることになる。
【0004】
従って、吸収物質の絶対値を求めるためには、生体表面から生体内に照射された光が光計測部に達するまでの生体内での光路長と、その光路長における散乱特性や吸収特性を知る必要がある。
【0005】
光路長を解析する手法としては、モンテカルロシミュレーションや光拡散方程式の数値解析などが知られており、またそれらを応用して複雑な境界形状からなる領域について吸収物質の定量を可能にした手法やアルゴリズムが提案されている(例えば、「光強度に依存しない吸収成分定量法−光路長平均と分散を用いるMVS法−」、土屋ら、pp.118〜121、OSJ-BIOP 2000-34、「有限差分法を用いた新生児頭部モデルにおける光伝播解析」、福井ら、pp.20〜23、OSJ-BIOP 2001-6など)。また特開平8-29329号公報には、散乱体である生体内の吸収物質の吸光係数と、モンテカルロ法により求めた模擬的な散乱体(吸収物質を含まない散乱体)における光路長とを用いて吸収物質濃度の空間分布を求め表示する手法が提案されている。
【0006】
これらの技術を実際の患者に適用する場合には、その患者について計測部位の構造を知る必要があるが、上記技術はファントムや頭部モデルを用いた光路長解析手法に限定されており、そのまま実用することはできない。一般に、生体の形態を計測する装置として、MRI、X線CT、超音波診断装置等が実用化されているが、他のモダリティから得た形態情報は、同一患者の情報であっても測定箇所によって異なるため、そのまま用いることはできず、形態情報と生体光計測装置の計測部位との位置合わせのための処理が必要となる。また生体光計測とは別にこのような形態の計測を行っておく必要があるため、リアルタイム性を充足できない。
【0007】
ところで超音波診断装置はMRIやX線CTに比べ簡便に形態情報を得ることができる装置であり、この超音波診断装置を生体光計測装置と組み合わせて、生体光計測装置で得られる生体機能情報を補い、形態情報をも表示できるようにした装置が提案されている(特許第3251417号)。しかしこの装置では、生体機能情報と併せて形態情報を観察することはできるが、吸収物質の絶対値表示を行えない点は、従来の生体光計測装置と同様である。
【0008】
また従来の生体光計測装置では、光照射部と光計測部の間隔は、頭部の計測の場合、深度約30mmを計測するのに適した間隔に固定されており、被検体や目的部位に応じた調整を行うことはできなかった。このため、生体光計測で得られる情報の精度に限界があった。
【0009】
【発明が解決しようとする課題】
そこで本発明は、生体光計測においてリアルタイムで吸収物質やその変化量の絶対値を画像化することができる画像診断装置を提供することを目的とする。また本発明は、被検体に応じて光路長の調整、即ち光照射部と光計測部との間隔の最適化を行うことができ、これにより精度の高い計測を行うことが可能な画像診断装置を提供することを目的とする。
【0010】
【課題を解決するための手段】
上記目的を達成するため、本発明の画像診断装置は、光照射部と、光照射部から生体内に照射され、生体内又は生体表層で散乱・透過した光を検出する光計測部と、光計測部が検出した光の強度に対応する信号を処理し、生体内の吸収物質の変化量を求め、画像化する信号処理部と、光照射部と光計測部との間に配置され、生体に超音波を送波すると共に生体内部で反射された超音波を受波する超音波探触子と、超音波探触子で受波した超音波に対応する信号を処理し、生体の形態画像を形成する超音波画像形成部とを備え、前記信号処理部は、前記超音波画像形成部で形成された形態画像をもとに生体内に照射された光の光路長を推定する手段と、推定された光路長及び前記吸収物質の変化量をもとに前記変化量の絶対値を算出する手段とを備える。
【0011】
この画像診断装置によれば、超音波画像をもとに生体光計測における光路長を推定し、これによって生のデータである吸収物質の変化量を補正できるので、測定部位の光路長に依存しない正確な情報を得ることができる。
【0012】
また本発明の画像診断装置は、光照射部と、光照射部から生体内に照射され、生体内又は生体表層で散乱・透過した光を検出する光計測部と、光計測部が検出した光の強度に対応する信号を処理し、生体内の吸収物質の変化量を求め、画像化する信号処理部と、光照射部と光計測部との間に配置され、生体に超音波を送波すると共に生体内部で反射された超音波を受波する超音波探触子と、超音波探触子で受波した超音波に対応する信号を処理し、生体の形態画像を形成する超音波画像形成部とを備え、前記信号処理部は、前記超音波画像形成部で形成された形態画像をもとに生体内に照射された光の光路長及び到達深度を推定する手段と、推定された光路長及び到達深度をもとに光照射部及び光計測部の最適な位置を算出する手段とを備える。
【0013】
この画像診断装置によれば、超音波画像をもとに生体光計測における光路長を推定し、これによって光照射部及び光計測部の位置を最適化できるので、目的とする深度の情報を得ることができる。
【0014】
さらに最適化した光照射部及び光計測部の位置における光計測結果について、光路長を用いた補正を行なうことにより、目的の測定部位について光路長に依存しない正確な情報を得ることができる。
【0015】
本発明の画像診断装置は、光照射部及び光計測部がそれぞれ複数備えられた装置に好適に適用することができる。
【0016】
また本発明の画像診断装置は、表示の機能として、超音波画像形成部で形成された超音波形態画像と生体光計測で計測された吸収物質の変化量或いはその絶対値を同時に表示する手段を備える。
この画像診断装置によれば、各測定位置における形態画像と生体光計測による測定結果とを併せて観察できるほか、光路長やパラメータの設定を確実に行なうことができる。
【0017】
さらに、本発明の画像診断装置は、光照射部及び光計測部がそれぞれ複数備えられている場合には、表示の機能として、信号処理部は、全ての測定点について、変化量の絶対値を表示する手段を備える。
これにより複数の光照射部及び光計測部がカバーする広い領域について、吸収物質の変化を観察することができる。
【0018】
【発明の実施の形態】
以下、本発明の画像診断装置の実施形態を説明する。
【0019】
図1は、本発明の画像診断装置の一実施形態の概要を示す図である。この画像診断装置は、光トポグラフィユニット10と、超音波ユニット20と、プローブ30と、制御ユニット40とから構成される。プローブ30は、被検体50の光計測と超音波計測とを同時に行なうために、複数の光照射部と複数の光計測部とを所定の配列で配置すると共に、光照射部と光計測部との間に超音波探触子を配置したものであり、その構成の詳細は後述するが、複数の光照射部及び光計測部は光ファイバー60a、60bを介して光トポグラフィユニット10に接続され、超音波探触子はケーブルを介して超音波ユニット20に接続されている。
【0020】
光トポグラフィユニット10は、光計測を行なうためのユニットで、既存の生体光計測装置と同様の構成を有している。簡単に説明すると、生体内の特定の吸収物質が吸収する可視から近赤外の波長の光を発生するとともに、この光に複数の光照射部に対応して異なる変調を与えるレーザー発生部11と、この特定波長の光を生体に照射し、生体の表層を透過した光或いは表層で散乱した光(以下、まとめて透過光ともいう)を、複数の光計測部毎に検出するための光電変換器、例えばアバランシェダイオード12と、複数の光計測部からの信号をもとに、光照射部と光計測部との位置で決まる測定チャンネル毎の信号に変換するロックインアンプ13と、連続可変アンプ14とを備えている。連続可変アンプ14の信号は制御ユニット40に送られ、ここで吸収物質の変化量を示す信号(例えばヘモグロビン信号)に変換される。
【0021】
超音波ユニット20は、プローブ30に備えられた超音波探触子からの信号を受けて、被検体の超音波画像を形成するもので、通常の超音波診断装置と同様の構成を有している。この場合、複数の超音波探触子からの信号をそれぞれ処理するために、複数の超音波探触子と同数の接続を備え、複数の超音波探触子からの信号を同時処理することが可能なように構成されている。或いは、複数の超音波ユニット20を択一的に接続するための切り替え器が設けられ、この切り替え器を介して連続的に切り替えを行い、データを取得し、擬似的に同時計測を実現する。
【0022】
超音波ユニット20は、簡単に説明すると、超音波探触子に備えられた複数の超音波振動子に特定周波数のパルス信号を送出するとともに各超音波振動子からの信号を受信するアナログ送受信回路21、アナログ送受信回路21による信号の送受信を制御する制御回路22、アナログ送受信回路21からの信号を整相して加算し、受信ビーム信号を形成するデジタル整相回路23、受信ビーム信号に対しフィルタリング処理、圧縮処理、検波処理、エンハンス等の信号処理を施す信号処理回路24及び信号処理後の受信信号を画像データに変換するデジタルスキャンコンバータ(DSC)25を備えている。また図示していないが、通常の超音波診断装置が備えるドップラ演算部を備えていても良い。
【0023】
制御回路22、信号処理回路24及びDSC25はデータバス26に接続されており、互いに必要な信号の受け渡しを行なう。
【0024】
制御ユニット40は、例えば、中央演算装置(CPU)、記憶装置(HDD)、キイボード46やモニター44等の入出力装置等を備えたパーソナルコンピューター(PC)上に構築することができ、光トポグラフィユニットからの信号をA/D変換してPCに入力するための光トポグラフィインターフェイス部41、超音波ユニット20のDSC25からの画像信号をPCに入力するための画像インターフェイス部42、超音波ユニット20の制御回路22との間で制御信号の受け渡しを行なう制御インターフェイス部43を備えている。
【0025】
また図示していないが、CPUは光トポグラフィインターフェイス部41を介して入力された信号を元に、生体内の吸収物質の変化信号、例えばヘモグロビン信号を作成し、これをタイムコースとして或いは2次元トポグラフィ画像としてモニター44に表示させる。このようなCPUの機能は、通常の生体光計測装置の機能と同様である。更に本発明においてCPUは、画像インターフェイス部42或いはキイボード46を介して入力された超音波計測から得た情報やパラメータを利用して、生のデータ(吸収物質の変化信号)に対する絶対値演算や計測点の最適化処理などを行なう。
【0026】
光トポグラフィインターフェイス部41、画像インターフェイス部42及び制御インターフェイス部43は、PCバス45に接続されており、互いに必要な信号の受け渡しを行なう。また画像インターフェイス部42は、モニター44に接続されており、超音波ユニット20で形成された画像をモニター44に表示させる。
【0027】
また制御ユニット40は、好適には、本発明の画像診断装置の動作モードを切り替えるための手段を備えている。動作モードは、光トポグラフィユニット20のみを動作させる光トポモード、超音波ユニット30のみを動作させる超音波モード、光トポグラフィと超音波計測の複合計測を行なう複合計測モードである。動作モードの切り替えためは、例えば、動作切り替えボタンによって実現できる。或いは、PCのモニターにモード選択メニューを表示し、ユーザーが所望のモードを選択するようにしてもよい。
【0028】
プローブ30は、被検体への光の照射及び被検体からの透過光の検出を行なうとともに、超音波探触子として機能するものであり、図2に示すように、フレキシブルな薄板状の装着部材(頭部の場合にはヘッドギア)31に、光照射部32a及び光計測部32bから成る光トポプローブ32と、超音波探触子35とが固定されている。装着部材31には、図示しない固定用の紐或いはベルトが備えられ、これによって被検体の所定部位、例えば頭部にプローブ30を固定することができる。
【0029】
光トポプローブ32は、それぞれ着脱可能な一対のソケット33、34から構成され、ソケット33の一方は光ファイバ60(a,b)の先端に固定され、他方は装着部材31に固定されている。光ファイバ先端に固定された一方のソケットを装着部材31に固定された他方に嵌着することにより、光ファイバ先端を被検体の装着部位に当接することができる。光トポプローブ32は光照射部32aと光計測部32bとが交互に配列し、これらの中間点が測定点となる。なお、図では一列の配列しか示していないが、光照射部32aと光計測部32bは、図示する配列方向と直交する方向(紙面に垂直な方向)にも配列し、広がりのある領域について計測するようになっている。
【0030】
超音波探触子35は、上述のように配列した光照射部32aと光計測部32bとの間に、固定治具36によって装着部材31に固定される。従って光照射部32aと光計測部32bとで決まる測定チャンネル毎に超音波探触子35からの情報が得られることになる。
【0031】
超音波探触子35の構造は、通常の超音波探触子と同様であり、例えば図3に示すように、被検体50に接触する側から順に、音響レンズ351、マッチング層352、積層振動子353、ソフトバッキング354、バッキング材355から成り、積層振動子353とソフトバッキング354との間に積層振動子353に信号を送るためのフレキシブル基板356が挿入されている。フレキシブル基板356の図示しない端部は、超音波ユニット20のアナログ送受信回路21に接続されている。
【0032】
次に、以上のような構成における画像診断装置の動作を説明する。本発明の画像診断装置は、上述のように光トポグラフィユニット10と超音波ユニット20とを備えているものであるから、これらを別個に動作させて通常の生体光計測装置や超音波診断装置として使用することも可能であり、またこれら各ユニットで得られた画像を別個に或いは同時にモニターに表示させることも可能である。しかし、本発明の特徴の一つは、超音波診断装置で得られた画像情報を元に光計測における光経路(光路長や到達深度)を導出し、この情報を生体光計測に反映させることである。具体的には、推定された光路長をもとに生体内吸収物質の変化量の絶対値を求めること、推定された光路長及び到達深度をもとに光照射部及び光計測部の最適な位置を算出し、光計測の高精度化を図ることである。以下、そのような動作について説明する。
【0033】
最初に、本発明の画像診断装置を用いた生体内吸収物質の変化量の絶対値を求める手順を図4を参照して説明する。図4は、絶対値算出のための手順の一例を示すフロー図である。この実施形態では、生体光計測の代表的な計測部位である頭部を対象として、その血流内ヘモグロビン値を計測する場合を説明する。
【0034】
プローブ30に設けられた複数の超音波探触子35の順次を駆動し、それぞれ脳表画像を作成する(ステップ401)。DSC25で作成された画像信号は、画像インターフェイス42を介して、制御ユニット40の記憶装置に記憶されるとともにモニター44に表示される。次にこれら画像を用いて頭部の構造をモデル化する(ステップ402)。頭部構造のモデル化は、例えば図5(a)、(b)に示すように、各超音波探触子35に対応するチャンネル毎に、モニター44にステップ401で作成した超音波画像を表示させ、各層(頭皮、頭蓋骨、脳髄液)の厚さd1、d2、d3を求めることにより行なう。次いで生体光計測で用いる光に対する各層の吸収係数、散乱係数等のパラメータを設定する(ステップ403)。光トポグラフィユニット10で発生する光は近赤外光であり、例えば780nm、830nmの二波長の光が用いられる。これら波長の光について各層の吸収係数及び散乱係数は、文献で知られた値であり、予め設定しておいたものを選択するようにしてもよい。
【0035】
各チャンネルについて、頭部モデル化、パラメータの設定がなされると、それをもとにチャンネル毎に当該超音波探触子を挟んで配置された光照射部と光計測部との間の光経路、光路長、光最深部到達位置を推定する(ステップ404)。既に述べたように、頭部モデルについて光路長等を求める手法は公知であり、例えばモンテカルロ法、光拡散方程式の数値解析、これらの改良法(例えば、「有限差分法を用いた新生児頭部モデルにおける光伝播解析」、福井ら、OSJ-BIOP 6001-6)を採用することができる。
【0036】
図6に、一つのチャンネル(第kチャンネル)について求めた光経路、平均光路長Lk、最深部到達距離Mkを模式的に示す。これら光経路、平均光路長及び最深部到達距離を、例えば図7に示すように、超音波画像に重ねてプローブ30の配列に対応させて表示する(ステップ405)。この表示を利用して、超音波画像に推定処理結果を順次重ねていくことにより、チャンネル毎のモデル化及びパラメータ設定(ステップ402、403)と(ステップ404)を確認的に実行することができる。
【0037】
各チャンネルについて光経路、平均光路長Lk、最深部到達距離Mkを推定した後、光トポグラフィユニット10を駆動し、プローブ30の各光照射部32aを介して被検体50の頭部に所定の波長の光を照射すると共に透過光を光測定部32bを介して光トポグラフィユニット10で検出する(ステップ406)。ここで計測される値は、各測定点(チャンネル)毎のヘモグロビンの変化量であり、各チャンネル毎の光路長のバラツキを含む値である。そこで、各チャンネル毎の計測値C1、C2・・・Cnを、そのチャンネルについて推定された推定光路長Lで除し、バラツキを補正した絶対値としての変化量C’1、C’2・・・C’nを得る(ステップ407)。
【0038】
【数1】
【0039】
この変化量の絶対値を用いて、光計測信号のタイムコース及び2次元トポグラフィを作成し、表示する(ステップ408)。表示の一例を図8に示す。この表示では、表示部(a)、(b)にそれぞれ所定のチャンネルの超音波断層画像と光計測信号を表示し、表示部(c)に各チャンネルの光計測信号値を等高線状に表示した画像(トポグラフィ画像)を、(d)に頭部全体における測定位置を示す画像を表示している。表示部(a)、(b)については、例えば、表示部(d)に表示された測定位置のチャンネルをマウス等の入力装置で選択することにより、所望のチャンネルの情報を表示するようにしてもよい。
【0040】
こうして表示されたトポグラフィは、チャンネルの光路長のバラツキを含まない画像であるので、測定点のヘモグロビン変化を正確に反映したものであり、診断の信頼性を向上することができる。
【0041】
尚、以上の実施形態では、ヘモグロビンの変化量について、それを推定された光路長で除した絶対値を求めた場合を示したが、ヘモグロビン濃度について、このような絶対値を求めることも可能である。即ち、照射光強度をI0、検出光強度をId、光路長Lkに存在する吸収物質(例えばヘモグロビン)の吸収係数をε、その濃度をCとすると、次式が成り立つ。
【0042】
【数2】
【0043】
ここでDsは減衰定数であり、直接求めることは困難であるため、照射光として波長の異なる2種類の光を用い、各波長毎に計測を行いDsを消去する。こうして求めた吸収物質の濃度を光路長で除することにより、光路長のバラツキを補正したヘモグロビン濃度(光路長における平均値)を求めることができる。
【0044】
次に本発明の画像診断装置を用いて、光照射部及び光計測部の位置の最適化を行なう実施形態を説明する。光照射部及び光計測部の位置の最適化を実現するのに好適なプローブ30’の一実施形態を図9に示す。
【0045】
このプローブ30’は、光トポプローブ32として光照射部Sと光計測部Dとが交互にマトリックス状に配列し、これらの間に超音波探触子35が設けられる点は、図2のプローブ30と同様である。また光ファイバ先端を装着部材に固定するための機構(ソケット)及び超音波探触子35の固定方法も図2のプローブ30と同様である。
【0046】
但し、このプローブ30’は、装着部材31の光トポプローブ固定位置にソケット34を移動するための溝37が形成されており、ソケット34の固定位置を調整できるように構成されている。溝37は、例えば、光トポプローブ32の配列方向に対し45°の角度を持つ十字形であり、ソケット34は、装着部材31の面に垂直な方向への移動が規制された状態で溝37に嵌合している。このような構成において、ソケット34を十字の溝37に沿って移動することにより、隣接するソケットとの相対位置を変化させることができる。
【0047】
次に、このような装置の機構を前提として、光照射部及び光計測部の位置の最適化の手順を図10を参照して説明する。図10は、最適化の手順の一例を示すフロー図である。
【0048】
まずプローブ30’に設けられた複数の超音波探触子35の順次を駆動し、それぞれ脳表画像を作成する(ステップ101)。これら画像を用いて各層(頭皮、頭蓋骨、脳髄液)の厚さd1、d2、d3を求め、頭部の構造をモデル化する(ステップ102)。次いで各層の近赤外光の吸収係数、散乱係数等のパラメータを設定する(ステップ103)。これら頭部モデル及びパラメータを用いて、各チャンネル毎にモンテカルロ法等により光経路、光路長Lk及び光最深部到達位置Mkを求め(ステップ104)、これを超音波画像に重ねて表示する(ステップ105)。以上のステップ101〜105は、図4に示したステップ401〜405と同様である。
【0049】
次に各チャンネル毎に表示された画像を確認しながら、光照射部と光計測部間の理想距離r’1、r’2、・・・r’nを設定する(ステップ106)。即ち、計測しようとする部位が表示された光最深部到達位置Mkがよりも深い場合には、距離を大きく設定し、計測しようとする部位がMkよりも浅い場合には距離を小さく設定する。
【0050】
こうして各チャンネルについて理想距離が設定されると、次式で表される実際の距離r1、r2、・・・rnと理想距離との差の評価量eを最小とするプローブ位置を算出する(ステップ107)。
【0051】
【数3】
【0052】
プローブ位置は、例えば図9に示す十字溝37の交差部を原点として、右上、右下、左上、左下の4つの方向のいずれかについて何ミリメートルというように求める。これは図9に示す複数の光トポプローブ32の一つを原点とする座標を設定し、その座標上で上述の計算により、各光トポプローブの最適位置を算出し、算出された位置に最も近い溝37上の位置を指定することにより実現できる。
【0053】
こうして実質的に理想距離となるプローブ位置が求められたならば、それに従い光トポプローブ32を移動した後(ステップ108)、光トポグラフィユニットを作動させて光計測を行なう(ステップ109)。これにより測定すべき部位からの情報を確実に測定することができる。
【0054】
この場合にも、新たに設定された光照射部と光測定部の位置について、ステップ104に戻り光路長を求めることができ(ステップ110)、この光路長を用いることにより、計測したヘモグロビン変化量の絶対値を算出し(ステップ111)、それを用いた光トポグラフィを表示することができる(ステップ112)。またステップ104に戻るのではなく、光路長と光照射部及び光測定部間の距離との相関の概算を予め求めておき、ステップ106で設定した光照射部及び光測定部間の理想距離に対応する光路長を用いることも可能である。この場合には、光路長自体の正確性は劣るが、時間のかかるモンテカルロ法等の演算を省くことができ、測定のリアルタイム性を向上することができる。
【0055】
尚、以上の実施形態では、光照射部と光測定部とをマトリックス状に配置し、2次元の生体光計測を行なう場合を説明したが、例えば、図11に示すように、一対の光照射部と光測定部との間に超音波探触子を配置した構造についても適用することができる。
【0056】
【発明の効果】
本発明によれば、生体光計測装置と超音波診断装置とを組み合わせ、超音波画像から得られる情報を利用して生体光計測の高精度化を図ることができる。
【図面の簡単な説明】
【図1】本発明の画像診断装置の概要を示す図
【図2】本発明の画像診断装置用のプローブの一実施形態を示す図
【図3】本発明の画像診断装置用の超音波探触子の一実施形態を示す図
【図4】本発明の画像診断装置の動作の一実施形態を示すフロー図
【図5】本発明による頭部構造のモデル化を説明する図
【図6】本発明による光経路の推定結果を示す図
【図7】本発明の画像診断装置の表示例を示す図
【図8】本発明の画像診断装置の表示例を示す図
【図9】本発明の画像診断装置用のプローブの他の実施形態を示す図
【図10】本発明の画像診断装置の動作の他の実施形態を示すフロー図
【図11】本発明の画像診断装置用のプローブの他の実施形態を示す図
【符号の説明】
10・・・光トポグラフィユニット
20・・・超音波ユニット
30・・・プローブ
32・・・光トポプローブ
35・・・超音波探触子
40・・・制御ユニット
60・・・光ファイバ[0001]
[Technical field to which the invention belongs]
The present invention relates to an image diagnostic apparatus, and more particularly, to an image diagnostic apparatus for combining a living body optical measurement apparatus and an ultrasonic diagnostic apparatus and imaging a change amount of an absorbing substance in a living body that is a scatterer.
[0002]
[Prior art]
A biological light measurement device has been developed and put to practical use as a diagnostic device that simply and noninvasively measures the oxygen metabolism function of a living body using light having a wavelength in the visible to near-infrared region. In the biological light measurement device, a specific absorbing substance such as hemoglobin that exists in the living body changes the absorbability of light of a specific wavelength between the oxygen state and the deoxygenated state, thereby improving the oxygen metabolism function of the living body. It is displayed as an image. As a further development of the biological light measurement device, there is an optical topography device that can display a relatively wide range of biological function information as an image by arranging a plurality of light irradiation units and a plurality of light measurement units in a matrix. . The optical topography device can obtain biological information such as the identification of the active site of the brain and the focal position at the time of seizure of an epileptic patient, and is becoming increasingly important as a diagnostic device.
[0003]
However, since such a biological light measurement device measures the attenuation amount (change) of light in the optical path from the light irradiation unit to the light measurement unit, it cannot directly measure the absolute value of the absorbing substance. . In other words, the path through which the light irradiated between the light irradiation unit and the light measurement unit passes depends not only on the distance between the light irradiation unit and the light measurement unit, but also on the structure of the living body. The attenuation amount is also different, and the value measured by the biological light measurement includes an optical path length factor depending on such a structure.
[0004]
Therefore, in order to obtain the absolute value of the absorbing substance, the optical path length in the living body until the light irradiated into the living body from the living body surface reaches the optical measurement unit, and the scattering characteristics and the absorption characteristics in the optical path length are known. There is a need.
[0005]
Known methods for analyzing the optical path length include Monte Carlo simulation and numerical analysis of the light diffusion equation, and by applying these methods, algorithms and algorithms that enable the quantification of absorbing substances in regions consisting of complex boundary shapes. (For example, “quantitative absorption component determination method independent of light intensity—MVS method using optical path length average and dispersion”), Tsuchiya et al., Pp. 118-121, OSJ-BIOP 2000-34, “finite difference Analysis of light propagation in a neonatal head model using the method ", Fukui et al., Pp.20-23, OSJ-BIOP 2001-6). Japanese Patent Application Laid-Open No. 8-29329 uses the extinction coefficient of the absorbing material in the living body, which is a scatterer, and the optical path length of the simulated scatterer (scattering material not including the absorbing material) obtained by the Monte Carlo method. Thus, a method for obtaining and displaying the spatial distribution of the concentration of absorbing substances has been proposed.
[0006]
When applying these techniques to an actual patient, it is necessary to know the structure of the measurement site for that patient, but the above technique is limited to optical path length analysis methods using phantoms and head models. It cannot be put to practical use. In general, MRI, X-ray CT, ultrasonic diagnostic equipment, etc. have been put to practical use as devices for measuring the form of a living body, but the morphological information obtained from other modalities can be measured even if it is information on the same patient. Therefore, it cannot be used as it is, and processing for positioning the form information and the measurement site of the biological light measurement device is required. Moreover, since it is necessary to perform such a form of measurement separately from biological light measurement, the real-time property cannot be satisfied.
[0007]
By the way, the ultrasonic diagnostic apparatus is an apparatus that can easily obtain morphological information as compared with MRI and X-ray CT, and the biological function information obtained by the biological optical measurement apparatus by combining the ultrasonic diagnostic apparatus with the biological optical measurement apparatus. Has been proposed (Japanese Patent No. 3251417). However, this apparatus can observe the morphological information together with the biological function information, but it is the same as the conventional biological light measuring apparatus in that the absolute value of the absorbing substance cannot be displayed.
[0008]
Also, in the conventional biological light measurement device, the distance between the light irradiation unit and the light measurement unit is fixed at an interval suitable for measuring a depth of about 30 mm in the case of head measurement. It was not possible to make adjustments accordingly. For this reason, there is a limit to the accuracy of information obtained by biological light measurement.
[0009]
[Problems to be solved by the invention]
Therefore, an object of the present invention is to provide an image diagnostic apparatus capable of imaging an absorbing substance and an absolute value of a change amount thereof in real time in biological light measurement. In addition, the present invention can adjust the optical path length according to the subject, that is, can optimize the interval between the light irradiation unit and the optical measurement unit, and thereby can perform high-accuracy measurement. The purpose is to provide.
[0010]
[Means for Solving the Problems]
In order to achieve the above object, an image diagnostic apparatus according to the present invention includes a light irradiator, an optical measurement unit that detects light that is irradiated into the living body from the light irradiator, and scattered / transmitted in the living body or the surface of the living body, A signal corresponding to the intensity of light detected by the measurement unit is processed to determine the amount of change in the absorption material in the living body, and is arranged between the signal processing unit to be imaged, the light irradiation unit, and the optical measurement unit. An ultrasonic probe for transmitting ultrasonic waves to the body and receiving ultrasonic waves reflected inside the living body, and processing signals corresponding to the ultrasonic waves received by the ultrasonic probe to form a morphological image of the living body Means for estimating the optical path length of the light irradiated into the living body based on the morphological image formed by the ultrasonic image forming unit; Means for calculating an absolute value of the change amount based on the estimated optical path length and the change amount of the absorbing material; Provided.
[0011]
According to this diagnostic imaging apparatus, it is possible to estimate the optical path length in living body light measurement based on the ultrasound image, and thereby correct the amount of change in the absorbing material, which is raw data, so it does not depend on the optical path length of the measurement site. Accurate information can be obtained.
[0012]
Further, the diagnostic imaging apparatus of the present invention includes a light irradiation unit, a light measurement unit that detects light that is irradiated into the living body from the light irradiation unit, and is scattered or transmitted through the living body or the surface of the living body, and light detected by the light measurement unit. It is placed between the signal processing unit, the light irradiation unit, and the optical measurement unit that processes the signal corresponding to the intensity of the light, determines the amount of change in the absorption material in the living body, and transmits the ultrasonic wave to the living body. An ultrasonic probe that receives ultrasonic waves reflected inside the living body, and an ultrasonic image that forms a morphological image of the living body by processing signals corresponding to the ultrasonic waves received by the ultrasonic probe Forming means, and the signal processing unit is estimated with means for estimating the optical path length and the reaching depth of the light irradiated into the living body based on the morphological image formed by the ultrasonic image forming unit Means for calculating the optimum positions of the light irradiation unit and the optical measurement unit based on the optical path length and the reaching depth .
[0013]
According to this diagnostic imaging apparatus, it is possible to estimate the optical path length in living body light measurement based on the ultrasound image, and thereby optimize the positions of the light irradiation unit and the light measurement unit, so that information on the target depth is obtained. be able to.
[0014]
Further, by correcting the optical measurement results at the positions of the optimized light irradiation unit and optical measurement unit using the optical path length, accurate information independent of the optical path length can be obtained for the target measurement site.
[0015]
The diagnostic imaging apparatus of the present invention can be suitably applied to an apparatus provided with a plurality of light irradiation sections and light measurement sections.
[0016]
Further, the diagnostic imaging apparatus of the present invention has means for simultaneously displaying, as a display function, the ultrasonic morphological image formed by the ultrasonic image forming unit and the amount of change or the absolute value of the absorbent measured by the biological light measurement. Prepare.
According to this diagnostic imaging apparatus, the morphological image at each measurement position and the measurement result by the biological light measurement can be observed together, and the optical path length and parameters can be set reliably.
[0017]
Furthermore, when the image diagnostic apparatus of the present invention includes a plurality of light irradiation units and light measurement units, as a display function, the signal processing unit calculates the absolute value of the change amount for all measurement points. Means for displaying are provided.
Thereby, it is possible to observe the change of the absorbing material in a wide area covered by the plurality of light irradiation units and the optical measurement unit.
[0018]
DETAILED DESCRIPTION OF THE INVENTION
Hereinafter, embodiments of the diagnostic imaging apparatus of the present invention will be described.
[0019]
FIG. 1 is a diagram showing an outline of an embodiment of an image diagnostic apparatus of the present invention. The diagnostic imaging apparatus includes an optical topography unit 10, an
[0020]
The optical topography unit 10 is a unit for performing optical measurement, and has the same configuration as that of an existing biological optical measurement device. Briefly, a
[0021]
The
[0022]
The
[0023]
The
[0024]
The control unit 40 can be constructed on, for example, a personal computer (PC) equipped with a central processing unit (CPU), a storage device (HDD), input / output devices such as a
[0025]
Although not shown in the figure, the CPU creates a change signal of a substance absorbed in the living body, for example, a hemoglobin signal, based on a signal input through the optical topography interface unit 41, and uses this as a time course or a two-dimensional topography. The image is displayed on the
[0026]
The optical topography interface unit 41, the
[0027]
The control unit 40 preferably includes means for switching the operation mode of the diagnostic imaging apparatus of the present invention. The operation mode is an optical topography mode in which only the
[0028]
The
[0029]
The optical topo probe 32 includes a pair of
[0030]
The
[0031]
The structure of the
[0032]
Next, the operation of the diagnostic imaging apparatus having the above configuration will be described. Since the image diagnostic apparatus of the present invention includes the optical topography unit 10 and the
[0033]
First, the procedure for obtaining the absolute value of the amount of change in the bioabsorbable substance using the diagnostic imaging apparatus of the present invention will be described with reference to FIG. FIG. 4 is a flowchart showing an example of a procedure for calculating an absolute value. This embodiment demonstrates the case where the hemoglobin value in the bloodstream is measured for the head which is a typical measurement part of biological light measurement.
[0034]
The plurality of
[0035]
When head modeling and parameter setting are made for each channel, the optical path between the light irradiation unit and the optical measurement unit arranged with the ultrasonic probe sandwiched for each channel based on that. Then, the optical path length and the position reaching the deepest part of the light are estimated (step 404). As already mentioned, methods for obtaining the optical path length and the like for the head model are well known. For example, the Monte Carlo method, numerical analysis of the light diffusion equation, these improved methods (for example, “newborn head model using finite difference method” Light propagation analysis in Fukui et al., OSJ-BIOP 6001-6).
[0036]
FIG. 6 schematically shows the optical path, average optical path length Lk, and deepest reach distance Mk obtained for one channel (k-th channel). These optical paths, average optical path lengths and deepest reach distances are displayed in correspondence with the arrangement of the
[0037]
After estimating the optical path, average optical path length Lk, and deepest reach distance Mk for each channel, the optical topography unit 10 is driven, and a predetermined wavelength is applied to the head of the subject 50 via each
[0038]
[Expression 1]
[0039]
Using the absolute value of the change amount, the time course and two-dimensional topography of the optical measurement signal are created and displayed (step 408). An example of the display is shown in FIG. In this display, ultrasonic tomographic images and optical measurement signals of predetermined channels are displayed on the display units (a) and (b), respectively, and optical measurement signal values of each channel are displayed in contour lines on the display unit (c). An image (topography image) is displayed, and an image showing the measurement position in the entire head is displayed in (d). For the display units (a) and (b), for example, by selecting the channel at the measurement position displayed on the display unit (d) with an input device such as a mouse, information on the desired channel is displayed. Also good.
[0040]
Since the topography displayed in this way is an image that does not include variations in the optical path length of the channel, it accurately reflects the change in hemoglobin at the measurement point, and the diagnostic reliability can be improved.
[0041]
In the above embodiment, the case where the absolute value obtained by dividing the change amount of hemoglobin by the estimated optical path length is obtained, but such an absolute value can be obtained for the hemoglobin concentration. is there. That is, if the irradiation light intensity is I0, the detection light intensity is Id, the absorption coefficient of the absorbing substance (for example, hemoglobin) existing in the optical path length Lk is ε, and its concentration is C, the following equation is established.
[0042]
[Expression 2]
[0043]
Here, Ds is an attenuation constant and is difficult to obtain directly. Therefore, two types of light having different wavelengths are used as irradiation light, and measurement is performed for each wavelength to erase Ds. By dividing the thus determined concentration of the absorbing substance by the optical path length, it is possible to determine the hemoglobin concentration (average value in the optical path length) corrected for variations in the optical path length.
[0044]
Next, an embodiment for optimizing the positions of the light irradiation unit and the light measurement unit using the image diagnostic apparatus of the present invention will be described. FIG. 9 shows an embodiment of a
[0045]
In this probe 30 ', the light irradiation section S and the light measurement section D are alternately arranged in a matrix form as the optical topo probe 32, and the
[0046]
However, the
[0047]
Next, a procedure for optimizing the positions of the light irradiation unit and the light measurement unit will be described with reference to FIG. FIG. 10 is a flowchart showing an example of the optimization procedure.
[0048]
First, the plurality of
[0049]
Next, the ideal distances r′1, r′2,..., R′n between the light irradiation unit and the light measurement unit are set while confirming the image displayed for each channel (step 106). That is, the distance is set to be larger when the deepest light reaching position Mk on which the part to be measured is displayed is deeper, and the distance is set to be smaller when the part to be measured is shallower than Mk.
[0050]
When the ideal distance is thus set for each channel, the probe position that minimizes the evaluation amount e of the difference between the actual distances r1, r2,. 107).
[0051]
[Equation 3]
[0052]
For example, the probe position is obtained in millimeters in any of the four directions of upper right, lower right, upper left, and lower left with the intersection of the
[0053]
When the probe position that is substantially the ideal distance is found in this way, the optical topography probe 32 is moved accordingly (step 108), and then the optical topography unit is operated to perform optical measurement (step 109). As a result, information from the site to be measured can be reliably measured.
[0054]
In this case as well, the optical path length can be obtained for the newly set positions of the light irradiation unit and the light measurement unit by returning to step 104 (step 110), and by using this optical path length, the measured hemoglobin change amount Can be calculated (step 111), and an optical topography using the absolute value can be displayed (step 112). Instead of returning to step 104, an approximate correlation between the optical path length and the distance between the light irradiation unit and the light measurement unit is obtained in advance, and the ideal distance between the light irradiation unit and the light measurement unit set in
[0055]
In the above embodiment, the case where the light irradiation unit and the light measurement unit are arranged in a matrix and two-dimensional biological light measurement is performed has been described. For example, as shown in FIG. The present invention can also be applied to a structure in which an ultrasonic probe is disposed between the unit and the light measurement unit.
[0056]
【The invention's effect】
ADVANTAGE OF THE INVENTION According to this invention, the biological optical measurement apparatus and an ultrasonic diagnostic apparatus can be combined, and the high precision of biological optical measurement can be achieved using the information obtained from an ultrasonic image.
[Brief description of the drawings]
FIG. 1 is a diagram showing an outline of an image diagnostic apparatus according to the present invention.
FIG. 2 is a diagram showing an embodiment of a probe for a diagnostic imaging apparatus according to the present invention.
FIG. 3 is a diagram showing an embodiment of an ultrasonic probe for an image diagnostic apparatus according to the present invention.
FIG. 4 is a flowchart showing an embodiment of the operation of the diagnostic imaging apparatus of the present invention.
FIG. 5 is a diagram illustrating modeling of a head structure according to the present invention.
FIG. 6 is a diagram showing an optical path estimation result according to the present invention.
FIG. 7 is a view showing a display example of the diagnostic imaging apparatus of the present invention.
FIG. 8 is a view showing a display example of the diagnostic imaging apparatus of the present invention.
FIG. 9 is a view showing another embodiment of the probe for the diagnostic imaging apparatus of the present invention.
FIG. 10 is a flowchart showing another embodiment of the operation of the diagnostic imaging apparatus of the present invention.
FIG. 11 is a view showing another embodiment of the probe for the diagnostic imaging apparatus of the present invention.
[Explanation of symbols]
10 ... Optical topography unit
20 ... Ultrasonic unit
30 ... Probe
32 ... Optical topo probe
35 ・ ・ ・ Ultrasonic probe
40 ... Control unit
60 ・ ・ ・ Optical fiber
Claims (8)
前記信号処理部は、前記超音波画像形成部で形成された形態画像をもとに生体内に照射された光の光路長を推定する手段と、推定された光路長及び前記吸収物質の変化量をもとに前記変化量の絶対値を算出する手段とを備えたことを特徴とする画像診断装置。Processes signals corresponding to the intensity of light detected by the light measurement unit, the light measurement unit that detects the light that is irradiated into the living body from the light irradiation unit, and is scattered or transmitted through the living body or the surface of the living body. The signal processing unit that obtains and changes the amount of absorption material in the living body, and is arranged between the light irradiating unit and the optical measuring unit, and transmits an ultrasonic wave to the living body and reflects the ultrasonic wave inside the living body. Image diagnostic apparatus comprising an ultrasonic probe that receives a sound wave, and an ultrasonic image forming unit that processes a signal corresponding to the ultrasonic wave received by the ultrasonic probe and forms a morphological image of a living body Because
The signal processing unit includes means for estimating an optical path length of light irradiated into the living body based on the morphological image formed by the ultrasonic image forming unit, and the estimated optical path length and the amount of change in the absorbing substance And a means for calculating an absolute value of the amount of change based on the image diagnostic apparatus.
前記信号処理部は、前記超音波画像形成部で形成された形態画像をもとに生体内に照射された光の光路長及び到達深度を推定する手段と、推定された光路長及び到達深度をもとに光照射部及び光計測部の最適な位置を算出する手段とを備えたことを特徴とする画像診断装置。Processes signals corresponding to the intensity of light detected by the light measurement unit, the light measurement unit that detects the light that is irradiated into the living body from the light irradiation unit, and is scattered or transmitted through the living body or the surface of the living body. The signal processing unit that obtains and changes the amount of absorption material in the living body, and is arranged between the light irradiating unit and the optical measuring unit, and transmits an ultrasonic wave to the living body and reflects the ultrasonic wave inside the living body. Image diagnostic apparatus comprising an ultrasonic probe that receives a sound wave, and an ultrasonic image forming unit that processes a signal corresponding to the ultrasonic wave received by the ultrasonic probe and forms a morphological image of a living body Because
The signal processing unit is configured to estimate an optical path length and a reaching depth of light irradiated into the living body based on a morphological image formed by the ultrasonic image forming unit, and to estimate the estimated optical path length and the reaching depth. An image diagnostic apparatus comprising: means for calculating optimal positions of the light irradiation unit and the light measurement unit.
前記信号処理部は、前記超音波画像形成部で形成された形態画像をもとに生体内に照射された光の光路長及び到達深度を推定する手段と、推定された光路長及び到達深度をもとに光照射部及び光計測部の最適な位置を算出する手段と、光照射部及び光計測部の最適位置における、推定された光路長及び吸収物質の変化量をもとに前記変化量の絶対値を算出する手段とを備えたことを特徴とする画像診断装置。Corresponds to a plurality of light irradiation units, a plurality of light measurement units that detect light scattered in and transmitted through the living body or on the surface of the living body, and light intensity detected by the light measurement unit. It is arranged between the signal processing unit that processes the signal, determines the amount of change in the absorbed substance in the living body, and images it, and each light irradiation unit and the optical measurement unit adjacent thereto, An ultrasonic probe that transmits ultrasonic waves and receives ultrasonic waves reflected inside the living body, and a signal corresponding to the ultrasonic waves received by the ultrasonic probe, and processes the morphological image of the living body. An image diagnostic apparatus comprising an ultrasonic image forming unit to be formed,
The signal processing unit is configured to estimate an optical path length and a reaching depth of light irradiated into the living body based on a morphological image formed by the ultrasonic image forming unit, and to estimate the estimated optical path length and the reaching depth. Based on the estimated light path length and the amount of change in the absorbing material at the optimal positions of the light irradiation unit and the light measurement unit based on the means for calculating the optimal position of the light irradiation unit and the light measurement unit An image diagnostic apparatus comprising: means for calculating an absolute value of
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002239220A JP4230736B2 (en) | 2002-08-20 | 2002-08-20 | Diagnostic imaging equipment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002239220A JP4230736B2 (en) | 2002-08-20 | 2002-08-20 | Diagnostic imaging equipment |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2004073559A JP2004073559A (en) | 2004-03-11 |
JP4230736B2 true JP4230736B2 (en) | 2009-02-25 |
Family
ID=32022382
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002239220A Expired - Fee Related JP4230736B2 (en) | 2002-08-20 | 2002-08-20 | Diagnostic imaging equipment |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4230736B2 (en) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4517111B2 (en) * | 2004-06-07 | 2010-08-04 | 独立行政法人情報通信研究機構 | Brain function measuring device, brain function measuring method, and brain function measuring program |
JP4739878B2 (en) * | 2005-09-21 | 2011-08-03 | シチズンホールディングス株式会社 | Cerebral blood flow measuring device |
JP5216387B2 (en) * | 2008-03-31 | 2013-06-19 | 株式会社日立製作所 | Probe device |
JPWO2011027548A1 (en) * | 2009-09-04 | 2013-01-31 | パナソニック株式会社 | Probe and image reconstruction method using the probe |
JP5527658B2 (en) * | 2010-04-30 | 2014-06-18 | 浜松ホトニクス株式会社 | Scattering absorber measurement method and apparatus |
JP6144131B2 (en) * | 2012-10-31 | 2017-06-07 | 株式会社東芝 | Ultrasonic diagnostic equipment |
JP6154613B2 (en) * | 2013-01-18 | 2017-06-28 | 浜松ホトニクス株式会社 | Cross-sectional image measuring apparatus and measuring method |
GB2519335A (en) | 2013-10-17 | 2015-04-22 | Univ Loughborough | Opto-physiological sensor and method of design |
EP3136967B1 (en) * | 2014-04-29 | 2021-10-20 | The Board of Regents of The University of Texas System | Systems for detecting sub-tissue anomalies |
JP6486733B2 (en) | 2015-03-17 | 2019-03-20 | 株式会社東芝 | Ultrasonic diagnostic apparatus and biopsy apparatus |
US10499836B2 (en) | 2016-03-11 | 2019-12-10 | Fujita Medical Instruments Co., Ltd. | Oxygen saturation measuring sensor, and oxygen saturation measuring apparatus |
JP6060321B1 (en) * | 2016-04-11 | 2017-01-11 | 株式会社フジタ医科器械 | Oxygen saturation measuring sensor and oxygen saturation measuring device |
CN110946615B (en) * | 2019-11-19 | 2023-04-25 | 苏州佳世达电通有限公司 | Ultrasonic diagnostic apparatus and operation method using the same |
-
2002
- 2002-08-20 JP JP2002239220A patent/JP4230736B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2004073559A (en) | 2004-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2036487B1 (en) | Measurement method and measurement apparatus | |
JP5541662B2 (en) | Subject information acquisition apparatus and control method thereof | |
JP5460000B2 (en) | Imaging apparatus and imaging method | |
US8364414B2 (en) | Apparatus and method for processing biological information | |
US9339254B2 (en) | Object information acquiring apparatus | |
US20170071476A1 (en) | Image generating apparatus, image generating method, and program | |
JP4230736B2 (en) | Diagnostic imaging equipment | |
JP6440140B2 (en) | Subject information acquisition apparatus, processing apparatus, and signal processing method | |
KR101984824B1 (en) | Method and apparatus for analyzing elastography of tissue using ultrasound | |
WO2011052061A1 (en) | Photo-acoustic device | |
KR101273585B1 (en) | Ultrasound imaging apparatus and display method of ultrasound image | |
US20130116536A1 (en) | Acoustic wave acquiring apparatus and acoustic wave acquiring method | |
JP5773578B2 (en) | SUBJECT INFORMATION ACQUISITION DEVICE, CONTROL METHOD AND PROGRAM FOR SUBJECT INFORMATION ACQUISITION DEVICE | |
WO2012011242A1 (en) | Image information acquiring apparatus, image information acquiring method and image information acquiring program | |
JP6659766B2 (en) | Optical measuring device and ultrasonic diagnostic device | |
Liu et al. | Toward wearable healthcare: a miniaturized 3d imager with coherent frequency-domain photoacoustics | |
JP6049780B2 (en) | Photoacoustic device | |
JP2017038917A (en) | Subject information acquisition device | |
WO2015174273A1 (en) | Breast measurement method and measurement device | |
JP4739878B2 (en) | Cerebral blood flow measuring device | |
US20170265749A1 (en) | Processing apparatus and processing method | |
JP6000609B2 (en) | Subject information acquisition apparatus and control method thereof | |
JP2016152879A (en) | Subject information acquisition apparatus | |
JP6686066B2 (en) | Photoacoustic device | |
JP6336013B2 (en) | Photoacoustic device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050701 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080729 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080908 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20081202 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20081204 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111212 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111212 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121212 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131212 Year of fee payment: 5 |
|
LAPS | Cancellation because of no payment of annual fees |