(57)【要約】【課題】光音響装置により算出される酸素飽和度の値の精度を向上させる。【解決手段】被検体に複数の波長の光が照射されて発生した音響波に由来する信号を用いて、被検体の内部の酸素飽和度分布を示す画像データを取得する処理部を備える画像処理装置であって、処理部は、信号を用いた画像再構成により生成された、被検体の血管の特定位置における酸素飽和度を取得し、血管の特定位置における酸素飽和度を補正するための第１の補正量を取得し、第１の補正量に基づいて被検体の特定位置ではない部分の酸素飽和度を補正するための第２の補正量を取得する。【選択図】図１
2018102359.tif 000002
【特許請求の範囲】
【請求項１】
被検体に複数の波長の光が照射されて発生した音響波に由来する信号を用いて、前記被検体の内部の酸素飽和度分布を示す画像データを取得する処理部を備える画像処理装置であって、
前記処理部は、
前記信号を用いた画像再構成により生成された、前記被検体の血管の特定位置における酸素飽和度に関する情報を取得し、
前記血管の前記特定位置における酸素飽和度に関する情報を補正するための第１の補正量を取得し、
前記第１の補正量に基づいて前記被検体の前記特定位置ではない部分の酸素飽和度に関する情報を補正するための第２の補正量を取得する
ことを特徴とする処理装置。
【請求項２】
前記血管は動脈であることを特徴とする請求項１に記載の処理装置。
【請求項３】
前記処理部は、前記第２の補正量を用いて、前記血管の前記特定位置の周辺部の酸素飽和度に関する情報を補正する
ことを特徴とする請求項１または２に記載の処理装置。
【請求項４】
前記処理部は、前記第２の補正量を用いて、前記血管の前記特定位置から所定の距離に位置する領域の酸素飽和度に関する情報を補正する
ことを特徴とする請求項３に記載の処理装置。
【請求項５】
前記処理部は、入力装置を用いてなされた術者からの入力に基づいて前記特定位置を設定する
ことを特徴とする請求項１ないし４のいずれか１項に記載の処理装置。
【請求項６】
前記処理部は、前記画像再構成により得られた画像を解析することにより前記特定位置を設定する
ことを特徴とする請求項１ないし４のいずれか１項に記載の処理装置。
【請求項７】
前記処理部は、前記被検体の血管において、複数の前記特定位置における酸素飽和度に関する情報を取得する
ことを特徴とする請求項１ないし６のいずれか１項に記載の処理装置。
【請求項８】
前記処理部は、複数の前記特定位置のそれぞれについて取得された前記第２の補正量に基づく補間処理を行う
ことを特徴とする請求項７に記載の処理装置。
【請求項９】
前記処理部は、前記特定位置における、画像再構成により生成された酸素飽和度に関する情報と、光学測定により取得された酸素飽和度に基づいて、前記第１の補正量を取得する
ことを特徴とする請求項１ないし８のいずれか１項に記載の処理装置。
【請求項１０】
前記光学測定は、パルスオキシメータにより行われる
ことを特徴とする請求項９に記載の処理装置。
【請求項１１】
前記入力装置は、前記被検体の表面から所定の深さまでの酸素飽和度に関する情報の投影画像に基づく、前記術者による入力を受け付ける
ことを特徴とする請求項５に記載の処理装置。
【請求項１２】
前記被検体に超音波を送受信する超音波プローブをさらに有する
ことを特徴とする請求項１ないし１１のいずれか１項に記載の処理装置。
【請求項１３】
前記処理部は、前記超音波プローブによる前記超音波の送受信により取得された前記被検体の形状を用いて画像再構成を行う
ことを特徴とする請求項１２に記載の処理装置。
【請求項１４】
前記処理部は、前記超音波プローブを用いて取得されたドップラ信号により前記血管の動脈を検出する
ことを特徴とする請求項１２または１３に記載の処理装置。
【請求項１５】
被検体に複数の波長の光が照射されて発生した音響波に由来する信号を用いて、前記被検体の内部の酸素飽和度に関する情報分布を示す画像データを取得する処理方法であって、
前記信号を用いた画像再構成により生成された、前記被検体の血管の特定位置における酸素飽和度を取得する工程と、
前記血管の前記特定位置における酸素飽和度に関する情報を補正するための第１の補正量を取得する工程と、
前記第１の補正量に基づいて前記被検体の前記特定位置ではない部分の酸素飽和度に関する情報を補正するための第２の補正量を取得する工程と、
を有することを特徴とする処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、処理装置および処理方法に関する。
【背景技術】
【０００２】
近年、光イメージング技術の一つとして、光音響トモグラフィー（Ｐｈｏｔｏａｃｏｕｓｔｉｃ  Ｔｏｍｏｇｒａｐｈｙ）が提案されている。この原理は以下のようなものである。まず、パルスレーザ光などの光を被検体に照射して被検体内を伝搬・拡散させる。すると、その光エネルギーを吸収した被検体内の吸収体が熱膨張して光音響波が発生する。音響波探触子でその光音響波を検出し、検出した信号を再構成することで、光音響画像が得られる。
【０００３】
特許文献１は、近赤外光を用いた光音響装置が開示している。近赤外光は血液によく吸収されるため、光音響装置では血管の分布に関する情報を取得できる。また、複数の波長の光を用いることによって血管の酸素飽和度に関する情報を取得できる。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１４−０９４２２５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
このように、光音響装置を用いることで、被検体の光学特性に関する情報、例えば血管の分布や酸素飽和度の情報を得られる。また、特許文献１のような光音響装置においては、背景の光学係数が分かれば正確に酸素飽和度を算出できるとされている。
【０００６】
ここで通常、同じ血管から分岐する血管のそれぞれにおける酸素飽和度は、互いに近い値になると想定される。しかし、光音響装置が実際に算出した酸素飽和度の値を見ると、同じ血管から分岐する血管であっても、場所に依存して値が変動することがある。さらに、光音響装置が実際に算出した酸素飽和度の値を見ると、正常であれば９５〜９９％の範囲内で一定とされる動脈の酸素飽和度が、その範囲を超えた値として算出されることがある。このように、光音響装置の算出値が必ずしも実際の酸素飽和度値を反映できていない可能性があり、算出精度の向上が求められている。
【０００７】
本発明は上記課題に鑑みてなされたものである。本発明は、光音響装置により算出される酸素飽和度の値の精度を向上させることを目的とする。
【課題を解決するための手段】
【０００８】
本発明は、以下の構成を提供する。すなわち、
被検体に複数の波長の光が照射されて発生した音響波に由来する信号を用いて、前記被検体の内部の酸素飽和度分布を示す画像データを取得する処理部を備える画像処理装置であって、
前記処理部は、
前記信号を用いた画像再構成により生成された、前記被検体の血管の特定位置における酸素飽和度に関する情報を取得し、
前記血管の前記特定位置における酸素飽和度に関する情報を補正するための第１の補正量を取得し、
前記第１の補正量に基づいて前記被検体の前記特定位置ではない部分の酸素飽和度に
関する情報を補正するための第２の補正量を取得する
ことを特徴とする処理装置である。
【０００９】
本発明はまた、以下の構成を提供する。すなわち、
被検体に複数の波長の光が照射されて発生した音響波に由来する信号を用いて、前記被検体の内部の酸素飽和度に関する情報分布を示す画像データを取得する処理方法であって、
前記信号を用いた画像再構成により生成された、前記被検体の血管の特定位置における酸素飽和度を取得する工程と、
前記血管の前記特定位置における酸素飽和度に関する情報を補正するための第１の補正量を取得する工程と、
前記第１の補正量に基づいて前記被検体の前記特定位置ではない部分の酸素飽和度に関する情報を補正するための第２の補正量を取得する工程と、
を有することを特徴とする処理方法である。
【発明の効果】
【００１０】
本発明によれば、光音響装置により算出される酸素飽和度の値の精度を向上させることができる。
【図面の簡単な説明】
【００１１】
【図１】本発明の光音響装置の構成を示す図
【図２】本発明の処理を示すフローチャート
【図３】本発明が提供する解析画面を示す図
【図４】解析画面における血管の位置を説明する図
【図５】本発明における周辺部の補正量を決める方法を説明する図
【発明を実施するための形態】
【００１２】
以下に図面を参照しつつ、本発明の好適な実施の形態について説明する。ただし、以下に記載されている構成部品の寸法、材質、形状およびそれらの相対配置などは、発明が適用される装置の構成や各種条件により適宜変更されるべきものである。よって、この発明の範囲を以下の記載に限定する趣旨のものではない。
【００１３】
本発明は、被検体から伝搬する音響波を検出し、被検体内部の特性情報を生成し、取得する技術に関する。よって本発明は、被検体情報取得装置またはその制御方法、音響波受信装置またはその制御方法、あるいは被検体情報取得方法や信号処理方法として捉えられる。本発明はまた、画像処理装置または画像処理方法としても捉えられる。本発明はまた、これらの方法をＣＰＵやメモリ等のハードウェア資源を備える情報処理装置に実行させるプログラムや、そのプログラムを格納した、コンピュータにより読み取り可能な非一時的な記憶媒体としても捉えられる。
【００１４】
本発明の被検体情報取得装置には、被検体に光（電磁波）を照射すると光音響効果により被検体内で発生する音響波を受信して、被検体の特性情報を取得する、光音響イメージング装置を含む。この場合、特性情報とは、受信された光音響波に由来する受信信号を用いて生成される、被検体内の複数位置のそれぞれに対応する特性値の情報である。
【００１５】
光音響測定により取得される特性情報は、光エネルギーの吸収量や吸収率を反映した値である。例えば、単一の波長の光照射によって生じた音響波の発生源、被検体内の初期音圧、あるいは初期音圧から導かれる光エネルギー吸収密度や吸収係数を含む。また、互いに異なる複数の波長により得られる特性情報から、組織を構成する物質の濃度を取得できる。物質濃度として酸素化ヘモグロビン濃度と脱酸素化ヘモグロビン濃度を求めることに
より、酸素飽和度分布を算出できる。また、物質濃度としては、グルコース濃度、コラーゲン濃度、メラニン濃度、脂肪や水の体積分率なども求められる。
【００１６】
被検体内の各位置の特性情報に基づいて、二次元または三次元の特性情報分布が得られる。分布データは画像データとして生成され得る。特性情報は、数値データとしてではなく、被検体内の各位置の分布情報として求めてもよい。すなわち、初期音圧分布、エネルギー吸収密度分布、吸収係数分布や酸素飽和度分布などの分布情報である。
【００１７】
本発明でいう音響波とは、典型的には超音波であり、音波、音響波と呼ばれる弾性波を含む。トランスデューサや音響変換素子等により音響波から変換された電気信号を音響信号とも呼ぶ。ただし、本明細書における超音波または音響波という記載は、それらの弾性波の波長を限定する意図ではない。光音響効果により発生した音響波は、光音響波または光超音波と呼ばれる。光音響波に由来する電気信号を光音響信号とも呼ぶ。
【００１８】
以下の実施形態では、被検体にパルス光を照射し、光音響効果により被検体から音響波を受信し解析する事により被検体内の光吸収体の分布を取得する光音響装置を取り上げる。このような光音響装置は、人や動物の血管疾患や悪性腫瘍などの診断や化学治療の経過観察に好適である。被検体の例として、被検者の乳房や手のような生体の一部、マウスなどヒト以外の動物、無生物、ファントムなどを挙げられる。
【００１９】
＜実施例１＞
本実施例における光音響装置を図１に示す。図１（ａ）は、光音響装置の構造を示す模式図であり、探触子ユニット１００の断面図と、装置制御部のブロック図とを示す。図１（ｂ）は、探触子ユニット１００の上面図である。
【００２０】
（全体的な構成）
探触子ユニット１００は、半球容器状の支持体１０１の内面に沿って、音響波探触子１０２がスパイラル状に５１２個配置された構成である。さらに支持体１０１の底部には、光照射部１０３から照射された、光音響波を発生させるための計測光が通過する空間１０５が設けられている。被検体は保持部材１０６に配置され、計測光はｚ軸の負の方向から照射される。なお、支持体１０１の内部や保持部材１０６の内部には、光音響波の伝達効率を向上させるために、水、油またはジェルなどの音響整合材を満たすことが好ましい。支持体１０１と被検体の相対的な位置関係は、ＸＹステージのような走査機構（不図示）によって変えることができる。走査機構が支持体１０１を走査することで、被検体の広い範囲を測定できる。
【００２１】
被検体から伝搬し、音響波探触子１０２に受信された光音響波は、電気信号に変換されてデータ取得部１０７に出力される。データ取得部１０７により増幅やＡＤ変換を施された電気信号をデータ処理部１０９が再構成することにより、被検体内部の光学特性値分布を反映した三次元の光音響画像が得られる。
【００２２】
また本実施例の光音響装置は、音響波探触子１０２とは別に、被検体の形態測定に用いられるリニア型の超音波プローブ１０４を備えている。リニア型の超音波プローブ１０４は支持体１０１と共に走査可能な構成となっている。
【００２３】
装置制御部１０８は、光照射部１０３の光照射、データ取得部１０７の受信制御、ＸＹステージの移動、超音波プローブ１０４の送受信など、装置全体の制御に関する指令を行う。また、装置制御部１０８はユーザインターフェースを備えており、術者からの指示を元に、測定パラメータの変更、測定の開始・終了、画像の処理方法の選択、患者情報や画像の保存、データの解析などを実行できる。さらに、詳細なデータ処理はデータ処理部１
０９で行われ、解析画面などが表示装置１１０にされる。
【００２４】
（光照射部）
光照射部１０３は、光音響波を励起する計測光の光源や、計測光を被検体に伝達する光伝送部を備える。光源としては、大出力を得るためにレーザー光源が望ましい。ただし、発光ダイオードやフラッシュランプ等でもよい。レーザーを用いる場合、固体レーザー、ガスレーザー、色素レーザー、半導体レーザーなど様々なものが使用できる。光音響波を効果的に発生させるためには、被検体の熱特性に応じて十分短い時間に光を照射させなければならない。そのため、被検体が生体である場合、パルス幅が１０〜５０ナノ秒程度のパルス光が好適である。また、パルス光の波長は、被検体内部まで光が伝搬する波長であることが望ましい。具体的には、生体の場合、７００ｎｍ以上１１００ｎｍ以下である。
【００２５】
本実施例では、固体レーザーであるチタンサファイアレーザーを用い、波長は７６０ｎｍおよび８００ｎｍとする。複数波長の光を用いることで波長ごとの吸収係数の違いを利用して、酸素飽和度に関する情報の算出が可能となる。光の照射のタイミング、波形、強度等は装置制御部１０８によって制御される。レーザーは１０Ｈｚで２波長を交互に照射することができる。装置制御部１０８は、探触子ユニット１００の走査と光照射および音響波受信のタイミングを連動させて適切に受信データを取得し、不図示のメモリに保存する。
【００２６】
光学系としては、バンドルファイバ、ミラー、プリズム、レンズなどが好適である。その他、光を所望の形状および強度で被検体に照射できるものであれば、どのような光学部品を用いてもよい。
【００２７】
（光音響用の音響波探触子）
音響波探触子１０２は、光音響波を受信する素子である。代表的にはＰＺＴ（圧電セラミックス）やＣＭＵＴ（容量性マイクロマシン探触子）があり、本実施例では、素子としてＣＭＵＴを用いる。音響波探触子１０２は単素子で、φ３ｍｍの開口を持ち、帯域は０．５ＭＨｚ〜４ＭＨｚであり、装置の分解能は０．５ｍｍである。低周波数帯域を含むことによって、太さ３ｍｍ程度の血管であっても良好な画像が取得できる。すなわち、血管の中が抜けてリング状に見えるような状況が発生し難くなる。サンプリング周波数は４０ＭＨｚで、２０４８サンプリングを行う。また、取得するデータは符号付きの１２ビットとする。
【００２８】
（データ取得部）
音響波探触子１０２によって電気信号に変換された信号は、データ取得部１０７に伝送され、増幅器により増幅され、Ａ／Ｄ変換器でデジタル信号に変換され、不図示の記憶手段（例えばフラッシュメモリ、ＨＤＤなど）にデータとして保存される。このとき、データ取得部１０７または装置制御部１０８が備える記憶手段を用いても良いし、外部の記憶装置を用いてもよい。すなわち本発明は、光音響装置に関する発明としても捉えられるし、光音響装置が取得して記憶手段に保存したデータを処理する情報処理装置または情報処理方法に関する発明としても捉えられる。データ取得部１０７は、典型的には増幅器、Ａ／Ｄ変換器、ＦＰＧＡチップなどで構成される。
【００２９】
（装置制御部）
装置制御部１０８は、ＣＰＵやＧＰＵなどのプロセッサー、ＲＯＭやＲＡＭなどのメモリ、通信装置、ユーザインターフェースなどの資源を備え、メモリに展開されたプログラムの各ステップに従って動作する情報処理装置である。装置制御部１０８の各機能は、複数の装置を組み合わせて実現してもよい。ユーザインターフェースとしては、マウス、キーボード、タッチパネルなどの、ユーザーが情報を入力する入力装置や、音声や画像によ
る提示装置などがある。
【００３０】
（リニア型の超音波プローブ）
リニア型の超音波プローブ１０４は、被検体に超音波を送受信して、形態画像やドップラ画像を得る。本実施例ではリニア型の超音波プローブとして、２５６個のＰＺＴ素子を配置されたプローブを用いる。素子の帯域は５ＭＨｚ〜１０ＭＨｚである。また、サンプリング周波数は４０ＭＨｚで２０４８サンプリングを行う。リニア型の超音波プローブ１０４を設ける代わりに、探触子ユニット１００を用いて超音波の送信とエコー波の受信を行ってもよい。
【００３１】
（支持体）
支持体１０１としては、音響整合材の重量や、走査機構による移動に対する耐久性を持つように、ある程度の強度を持つ部材が好ましい。典型的には金属や樹脂が用いられる。また本実施例では、半球容器状の支持体１０１を用いる。これにより、複数の音響波探触子１０２それぞれの感度の高い方向（指向軸）が集まる高感度領域が形成されるので、被検体の特性情報分布の精度が向上する。図１（ｂ）に示すように、本実施例では複数の音響波探触子１０２を３次元スパイラル状に配置しているが、配置方法はこれに限られない。支持体の形状は半球に限られず、リニア状、平面状、お椀状などでも良い。
【００３２】
（保持部材）
保持部材１０６としては、被検体を支持する強度を備えるとともに、光と音響波を透過させる特性を有するものを用いる。好ましい材質として、ポリエチレンテレフタラートやアクリルなどがある。また、薄いゴムやフィルムなども利用できる。保持部材１０６を被検体に沿った形状にすることが好ましい。例えば被検体が乳房であれば、カップ状の保持部材１０６を用いる。
【００３３】
（データ処理部）
データ処理部１０９は、記憶手段に保存された光音響信号を用いた画像再構成により、被検体内部の光音響画像および機能画像である酸素飽和度画像を生成する。その他、光量分布の計算や背景の光学係数取得に必要な情報処理、信号補正など所望の処理を実施する。データ処理部１０９はプロセッサーやメモリなどを備える情報処理装置により構成できる。プロセッサーで動作するプログラムの各モジュールによりデータ処理部の各機能を実現できる。なお、装置制御部１０８とデータ処理部１０９を共通の情報処理装置により構成しても良い。データ処理部１０９は、本発明の処理部に相当する。
【００３４】
光音響画像や酸素飽和度画像は、表示装置１１０であるディスプレイに表示される。ディスプレイとしては、高解像度でカラー表示のできる３０型、コントラスト比１０００：１以上が望ましい。液晶ディスプレイ、プラズマディスプレイ、有機ＥＬディスプレイなど任意のディスプレイを、表示装置１１０として利用できる。
【００３５】
（画像再構成）
詳細な画像再構成はデータ処理部１０９にて行われる。画像再構成は、ユニバーサルバックプロジェクション法や整相加算法など既知の再構成手法を用いる。ここでは、ユニバーサルバックプロジェクション法を用いる。光音響測定で発生する初期音圧分布Ｐ（ｒ）は数式（１）で表わされる。
【数１】
2018102359.tif 000003
【００３６】
このとき投影データに相当する項ｂ（ｒ０，ｔ）を、数式（２）に示す。ここで、ｐｄ（ｒ０）は音響波探触子１０２で検出される光音響信号、ｒ０はそれぞれの音響波探触子１０２の位置、ｔは時間、Ω０は音響波探触子１０２の立体角である。データ取得部１０７で取得したデータを数式（１）に基づいて処理をすることにより初期音圧分布Ｐ（ｒ）を得ることができる。
【数２】
2018102359.tif 000004
【００３７】
次に、吸収係数分布は、初期音圧分布Ｐ（ｒ）から算出できる。吸収体に光が照射された時に発生する音圧Ｐ（ｒ）は、数式（３）で表される。
【数３】
2018102359.tif 000005
Γは弾性特性値であるグリューナイセン（Ｇｒｕｎｅｉｓｅｎ）係数であり、体積膨張係数（β）と音速（ｃ）の二乗の積を比熱（Ｃｐ）で割ったものである。μａは吸収体での吸収係数である。また、Φ（ｒ）は局所的な領域で吸収体に照射された光量である。数式（３）を吸収係数について解くことによって、吸収係数分布μａ（ｒ）を得ることができる。なお、背景の光学係数は、吸収体の吸収係数より十分小さいため吸収係数分布には表れない。
【００３８】
光量分布Φ（ｒ）は、深さ方向に一様に減衰するような場合、例えば数式（４）のように変数ｚを用いて表すことができる。
【数４】
2018102359.tif 000006
Φ０は、表面での入射光の光量である。μｅｆｆは、被検体内での平均的な等価減衰係数で、被検体内の背景の散乱係数μｂｓや吸収係数μｂａを反映したものであり、例えば数式（５）のように表わされる。
【数５】
2018102359.tif 000007
【００３９】
なお、被検体内の背景の散乱係数や吸収係数は、例えば近赤外分光装置などによって測定することができる。ただし、均質媒体を仮定したときの数値を算出するため皮膚、脂肪層、乳腺層、血管など様々な構造を有する乳房において正確に算出することは難しい。本発明においては、あらかじめ取得した背景の光学係数分布から性別、年齢、体重などに応じた値を選択しても良い。動脈と静脈が判別できる程度の画像が得られればよく、その後に補正をするためである。
【００４０】
次に、吸収係数分布からヘモグロビン分布を求める方法を示す。吸収体の吸収係数μａ（λ）は、単位体積当たりの酸化ヘモグロビンＣＨｂＯと還元ヘモグロビンＣＨｂＲの吸
収によって決まる。酸化ヘモグロビンと還元ヘモグロビンの吸収係数をそれぞれの波長でＥＨｂＯ（λ１），ＥＨｂＲ（λ１），ＥＨｂＯ（λ２），ＥＨｂＲ（λ２）とすると、数式（６）のように表わされる。なお、左辺は数式（３）によって得た吸収係数分布である。
【数６】
2018102359.tif 000008
【００４１】
数式（６）から酸化ヘモグロビンと還元ヘモグロビンについて変形するとそれぞれ数式（７）のようになる。
【数７】
2018102359.tif 000009
【００４２】
なお、総ヘモグロビン（ｔＨｂ）は酸化ヘモグロビンと還元ヘモグロビンを合わせた量であるため数式（８）のように表わされる。
【数８】
2018102359.tif 000010
【００４３】
また、酸素飽和度ＳｔＯは、総ヘモグロビンに示す酸化ヘモグロビンの割合であるため数式（９）のように表わされる。
【数９】
2018102359.tif 000011
【００４４】
（処理のフロー）
本発明の処理のフローチャートを図２に示す。
Ｓ１工程で、測定が開始される。この状態では、被検体が保持部材１０６に接触するように挿入されている。また、保持部材１０６と被検体の間には空気が入らないように密着し、被検体周辺には音響整合材である水が充填されている。なお、光音響画像を取得する前に超音波プローブ１０４で超音波測定を行ってもよい。超音波画像により、被検体の病変部などの座標を特定するためである。また、パルスオキシメータなどの光学測定によって酸素飽和度を測定しても良い。一般的な透過型のパルスオキシメータを用いて指における酸素飽和度を測定する。ただし、これに限られたものでなく反射型のパルスオキシメータによりその他の部位を測定しても良い。
【００４５】
Ｓ２工程で、光音響測定が行われる。まず、術者が装置制御部１０８のユーザインター
フェースから測定の指示を出す。それにより、走査機構であるＸＹステージが支持体１０１を所望の位置に移動させる。そして、光照射部１０３からパルス光を照射し、これと同期して音響波探触子１０２が光音響波を受信する。パルス光の照射の際には、支持体１０１をスパイラル状に移動させながら、７６０ｎｍと８００ｎｍの２波長の光を交互に照射する。最終的に、それぞれの波長で１０２４の位置で光音響波を取得する。撮像範囲は選択により、直径８０、１２０、１６０ｍｍとすることができる。なお、交互に異なる波長を照射して２波長分の光音響画像を得る方が、それぞれの波長で２回測定して得る光音響画像に比べて、体動などによる波長間の位置ずれの影響が抑制される。光音響波から変換され、さらにデジタル化された電気信号は、記憶装置に記憶される。
【００４６】
Ｓ３工程で、光音響画像から機能画像が生成される。データ処理部１０９は、まず、記憶装置から受信データを読み出して画像再構成を行い、波長ごとの初期音圧分布を得る。続いて、初期音圧分布と光量分布を用いて波長ごとの吸収係数分布を得る。そして、機能画像である酸素飽和度画像を得る。なお、機能画像は、初期音圧分布などにフィルタを用いてノイズ成分を除去した後の信号を用いて算出することができる。フィルタとしては一般的な平滑化フィルタ、ガウシアンフィルタ、メディアンフィルタ、バイラテラルフィルタ等を用いる。局所的なノイズに左右されることを防ぐためである。また、酸素飽和度は、吸収係数分布において、閾値以上の部分を選択的に表示することが望ましい。これにより血管の位置のみを表示することができる。
【００４７】
なお、光量分布は、被検体の背景の光学係数、および、被検体の形態情報を用いた計算により取得できる。形態情報については、本実施例の場合は超音波プローブ１０４が取得した超音波エコーデータに基づいて算出できる。また、再構成により得られた初期音圧分布画像から、被検体の境界を抽出しても良い。さらにまた、保持部材１０６に被検体を収める場合、簡易的に保持部材１０６の形状を被検体形状としても良い。
【００４８】
Ｓ４工程で、被検体において動脈だと判明している位置（特定位置）の設定が行われる。これについて図３に示される解析画面３００を用いて説明する。解析画面３００は表示装置１１０に表示され、機能画像表示部分３０１、ポインタ３０５、ポインタの中心のｘ座標、ｙ座標を示す補助線３０３および３０４、ポインタ部分のヒストグラム３０２、酸素飽和度に関する表示部３０６，３０７等から構成される。ヒストグラム３０２は、ポインタに含まれる１２×１２画素それぞれにおける酸素飽和度値（％）の数値範囲ごとに、出現頻度を積み上げたものである。
【００４９】
術者は解析画面３００において、入力装置を用いて、最大値投影に用いる深さ方向の画像の枚数、最大値投影のための深さ方向の開始位置、画像の拡大、縮小機能などを設定できる。ここでの最大値投影とは、ｘｙ平面の各点において、ｚ方向における最大の吸収係数値を取得し、その最大の吸収係数値を示す位置における酸素飽和度値を取得し、取得した酸素飽和度値をｘｙ平面上に並べることを言う。このようにすることで、吸収係数の高い血管位置における酸素飽和度が選択的に表示される。以後、このような画像を「酸素飽和度の投影画像」とも呼ぶ。
【００５０】
本実施例では画像の解像度に関し、一画素（ピクセルまたはボクセル）を幅０．０６２５ｍｍとしている。また画像化領域に関しては、５０×５０ｍｍの範囲を表示させている。また、ポインタ部分は１２×１２画素で構成され、術者が入力装置によって移動させることができる。さらに、ポインタ部分の酸素飽和度の平均値および分散値をそれぞれ表示部３０６，３０７に表示させることができる。
【００５１】
図４（ａ）は、図３の機能画像表示部分３０１に表示される、酸素飽和度の投影画像の一例である。ここでは、被検体の皮膚の所定の領域において、皮膚表面から所定の深さ（
ここでは１０ｍｍ）までの範囲を対象として、各位置での深さ方向における吸収係数の最大値を求め、吸収係数が最大値を示す深さにおける酸素飽和度値を求めた。図４（ａ）は、乳房を正面から見たＣｏｒｏｎａｌ方向の模式図である。
【００５２】
図４（ｂ）は、動脈において選択した複数の特定位置（マーキング４０４〜４０６）を示す図である。図４（ａ）の中心から領域４０１の方向に内胸動脈があり、領域４０３の方に乳頭がある。ここでは酸素飽和度の投影画像を白黒画像で示しているが、酸素飽和度の０％を青色、１００％を赤色に対応させて、酸素飽和度に応じてグラデーション表示を行うカラーマップを用いることも好ましい。
【００５３】
光音響画像から動脈を特定するときは、次のような情報（１）〜（５）を用いる。
（１）動脈と静脈は対になっている。
動脈から出た血液は体内の組織を通って静脈に戻るためである。２本の血管が近接する場合（例えば両者の距離が１０ｍｍ以内）、一方が動脈で他方が静脈と推定できる。当然、一本の血管から分岐する場合は同種の血管であると判断できる。被験体が乳房の場合、内胸動脈から血管が直接分岐している部分において、動脈と静脈が並走した伴行部分が観察されることがある。これは動静脈が対となっている典型例で、どちらかが確実に動脈であるため特に注目する場所である。図４（ａ）においては、領域４０１が伴行部分に相当する。領域４０２、領域４０３においても、動脈と静脈が対になっていると判断できる。ただし、実際にこれらの動脈と静脈の間で血液の流れがあることを保証するものではない。
【００５４】
（２）動脈と静脈では、動脈の酸素飽和度が高い。
酸素飽和度値は、動脈において９５％〜９９％程度、静脈において６０〜８０％程度とされている。これらの範囲内の数値を示す血管はそれぞれ動脈および静脈と推定できる。ただし、実際の測定では形状、背景の光学係数、ノイズなどの影響により、算出値がこの範囲に収まらないことがある。しかし、このような場合でも、浅い位置にある対となった動脈と静脈の酸素飽和度の相対的な大きさは逆転しにくい。数式（４）で表わされるように光量分布の誤差は、深くなるほど大きくなるためである。例えば、対になった血管の一方の酸素飽和度が１０５％で、他方が９０％程度になっている場合、相対的に酸素飽和度値の高い血管が動脈と推定できる。なお、酸素飽和度値の算出値の実際の値からの誤差は、被検体の表面に近い部位のほうが小さい。そこで、対になった血管のうち、表面に近い側の血管の酸素飽和度値を基準として比較を行うと良い。
【００５５】
（３）動脈と静脈を比べると一般的に動脈の方が細い。
例えば距離１０ｍｍ以内に２本の血管が近くにある場合、細い方が動脈と推定できる。ただし、酸素飽和度の値を優先的に用いて判断をする。
【００５６】
（４）連続する部分は同種である。
例えば、画像中に管状の吸収体があり、血管であるということは推定できているものの、部分的にノイズが多くて動脈か静脈かの判断ができない場合がある。あるいは、対になった血管が伴行している部分においては、２本の血管の酸素飽和度値が近い値（例えば９５％）として算出されてしまう場合がある。この原因は不明であるが、可能性として血管径が１ｍｍ程度の場合で近接することによる信号の混入が考えられる。また、実際の光量分布が仮定した物理モデルと異なっている可能性、血管が近接しているため動脈の酸素を消費しないまま静脈に戻るため正しい値である可能性も考えられる。これらの場合でも、連続する１本の血管は動脈か静脈のいずれかであり、途中で種類が切り替わることがないという事実を利用すれば、動脈を特定できる。すなわち、管状の吸収体の一部に動脈か静脈かを判定できない部分があったとしても、同じ管状の吸収体の他の部分で動脈であることを特定できていれば、その管状吸収体の全体が動脈であると判断できる。乳房の場合、
内胸動脈、外側胸動脈など特定の動脈から分岐しているので、連続している部分の一部が動脈であると判断できれば、その他の部分も比較的判断しやすくなる。
【００５７】
（５）静脈と比べると動脈の方が深い位置にある。
画像上深い位置から浅い位置に出た血管が先細りして無くなっているような場合は動脈と推定できる。一方、２〜３ｍｍ程度の太さのまま表面に沿うように走っている血管は、静脈と推定できる。なお、深い位置にあるかどうかの判断には、乳房を頭側の平面に投影するＡｘｉａｌ画像、それに垂直な頭尾方向と平行な平面に投影するＳａｇｉｔｔａｌ画像を用いても良い。ただし、酸素飽和度の値を優先的に用いて判断をする。
【００５８】
ところで、本実施例のように装置がリニア型の超音波プローブ１０４を備えており、所望の位置に超音波プローブ１０４を移動できる場合、パルスドップラ法によって拍動を観察することによって動脈を特定しても良い。細い血管においてドップラ信号を検出することは難しいが、内胸動脈に近い位置の２ｍｍ程度の太さの動脈であればドップラ信号を検出できる場合がある。
【００５９】
血管における特定位置の指定は、術者が上記情報に基づき行い、入力装置を用いて設定しても良いし、データ処理部が上記情報に基づく画像解析を行って設定しても良い。前者の場合、ユーザインターフェースが入力を受け付ける入力装置となる。なお、入力する特定位置が多いほど細かく補正の設定ができる。また、選択する血管は太い方が望ましく、例えば１ｍｍ以上とする。細い血管は体動などによる誤差の影響を受けやすく、またノイズの影響を受けやすいからである。このため、所望の太さ以上の血管のみを表示させても良い。太さの判断は、略等しい酸素飽和度の画素が連続しているときに、ｘｙｚ方向それぞれにおける画素数に基づいて判断できる。
【００６０】
さらに、データ処理部における画像処理によって動脈の位置を抜き出しても良い。動脈の抜き出しは、酸素飽和度画像において動脈位置の開始点を与えて、酸素飽和度の近い点を繋げていくことによって行う。血管が交差している部分、伴行している部分などは、術者が適宜判断して補正しても良い。また、血管像のテンプレート処理を行ってもよい。
【００６１】
動脈である位置の入力は、図４（ｂ）のように画像上にマーキング４０４〜４０６を配置することにより行う。マーキングは、図３のように画像上にポインタを表示してマウス等でクリックして配置する。この際、術者は酸素飽和度のヒストグラムや平均値を見ながら適切な位置を設定する。ヒストグラムにおいて、９０〜１０５％の酸素飽和度値を示す画素が、狭い範囲に集中している場合（例えば、酸素飽和度で±１％以内に多数の頻度が積み上がっている場合）は、ポインタのある位置が動脈位置であると判断可能である。逆に出現頻度がばらついているような場合は、ノイズを多く含んでいる、血管以外の部分を含んでいる、動脈と静脈が近接しているなどの理由により、動脈であるとは判断できない。なお、このような場所を設定した場合は、警告表示を出しても良い。
【００６２】
ところで、選定する動脈は正常部位であることが望ましい。正常部であれば、ほぼ均一な数値と考えられるからである。一方、病変部位などでは、酸素分圧の低下により、酸素飽和度が低下することがある。病変部位と正常部位を識別するためには、リニア型の超音波プローブ１０４で得た３次元の超音波画像や、ＭＲＩなどで別途取得した画像に基づいて、腫瘍部位と思われる範囲をマスクし、それ以外の正常部位のみ動脈であることを示すマーキングを設置できるようにしても良い。
【００６３】
Ｓ５工程で、動脈位置における補正量を計算する。一般的に動脈の酸素飽和度は体内で変わらないとされている。このため、パルスオキシメータなどの酸素飽和度測定装置を用いた光学測定で測定した値を、動脈の酸素飽和度値とすることができる。動脈位置におけ
る酸素飽和度が１０２％で、パルスオキシメータの測定値が９８％であれば、補正量は−４％となる。この補正量は、動脈の特定位置における補正量であり、本発明の第１の補正量に相当する。なお、局所的なノイズに補正量が左右されないようにするために、補正できる幅に制限をつけても良い。
【００６４】
Ｓ６工程で、動脈位置周辺における補正量を計算する。この補正量は、特定位置ではない部分の酸素飽和度を補正するためのものであり、本発明の第２の補正量に相当する。図５を用いて周辺部の補正方法を説明する。この工程では、周辺部に含まれる動脈および静脈の酸素飽和度が補正される。図５（ａ）は、動脈位置が一箇所の場合である。本実施例では、内胸動脈に近い図４（ｂ）のマーキング４０４において補正をする。内胸動脈に近い動脈は太いため、酸素飽和度値がより適切であると考えられる。この理由として、Ｓ３工程における算出値がパルスオキシメータでされる値に近いと考えられること、光音響装置から算出される酸素飽和度のサンプリング数を増やしてノイズの低減ができることが挙げられる。
【００６５】
図４（ｂ）におけるマーキング４０４に対応する座標は、座標５０１（ｘ１，ｙ１）で、補正量はｄ１（％）である。ここでは、座標（ｘ１，ｙ１）を中心に周辺部の補正量が可及的に少なくなるように、補正を行う。例えば、補正対象となる範囲を動脈の特定位置からの距離が３０ｍｍ以内の領域とする。これにより、補正の影響を動脈位置の周辺部に留めることができるので、局所的な補正の影響が被検体全体に反映されることがなくなる。なお、酸素飽和度の数値が全体的に一様にシフトしている場合は、補正を被検体の全体に反映させても良い。また、図４（ｂ）のように動脈を抜き出した場合は、動脈位置の補正量を算出した後、最短距離にある動脈の補正量を適用して周辺部を補正しても良い。
【００６６】
図５（ｂ）は、特定された動脈位置が二箇所ある場合である。図４（ｂ）におけるマーキング４０４と４０５に対応した座標は、それぞれ座標５０２（ｘ２，ｙ２，ｚ２）、５０３（ｘ３，ｙ３，ｚ３）で、補正量はそれぞれｄ２、ｄ３（％）である。なお、ｚ位置は、ｘｙ座標それぞれの位置においてに投影した最大値のある深さである。内胸動脈に近い座標５０２の方が深く（表面から遠く）、動脈から分岐した座標５０３の方が浅い（表面に近い）。周辺部の補正量は、ｚ１、ｚ２の位置の補正量に基づく補間処理や補外処理により、求めることができる。補間処理の方法は限定されず、直線で補間しても良いし、指数関数を用いても良い。この場合、同じ深さでは同じ補正量となる。
【００６７】
図５（ｃ）は、特定された動脈位置が三箇所の場合である。図４（ｂ）におけるマーキング４０４，４０５，４０６に対応した座標が、それぞれ座標５０４（ｘ１，ｙ１）、５０５（ｘ２，ｙ２）、５０６（ｘ３，ｙ３）で、補正量がそれぞれｄ１、ｄ２、ｄ３（％）である。この場合、平面の方程式によってｘｙ平面上の周辺部の補正量を求めることができる。なお、図４（ｂ）のマーキングを増やせばより細かく周辺部の補正量を求めることができる。
【００６８】
ところで、装置の特性上、動脈または静脈で酸素飽和度の値が実際より高めまたは低めに算出されるときは、動脈または静脈の一方のみをさらに補正しても良い。
当然上述した処理は適宜組み合わせても良い。このようにして酸素飽和度の補正を行い、酸素飽和度画像を得ることができる。
最後に、Ｓ７工程で、処理を終了する。
【００６９】
以上のように本発明によれば、光音響装置において、場所によらず略一定であるとされる動脈の酸素飽和度を基準にして、酸素飽和度の算出値を妥当な範囲に補正できる。その結果、酸素飽和度分布を示す表示画像の精度を向上させて良好な診断が可能になる。
【符号の説明】
【００７０】
１０３：光照射部、１０８：装置制御部、１０９：データ処理部
【図１】
2018102359.tif 000012
【図２】
2018102359.tif 000013
【図３】
2018102359.tif 000014
【図４】
2018102359.tif 000015
【図５】
2018102359.tif 000016