(57)【要約】【課題】画像データ生成に要する時間を抑制しつつ、画像中のアーチファクトを低減する方法を提供する。【解決手段】被検体から伝搬する音響波を、有限の面積を持つ平面上に所定の間隔で配置された複数の検出位置において受信し、電気信号に変換する素子と、電気信号を用いて被検体内部の特性情報を示す画像を構成する構成手段と、画像を周波数空間の信号に変換したのち、複数の検出位置間の距離ｄと、検出位置と画像上の各位置を結ぶ線が平面から垂直な方向となす角度θを用いて、周波数空間の信号から画像におけるストリークアーチファクトに相当する空間周波数成分を低減するアーチファクト低減手段を有する被検体情報取得装置を用いる。【選択図】図１
2015062551.tif 000002
【特許請求の範囲】
【請求項１】
被検体から伝搬する音響波を、有限の面積を持つ平面上に所定の間隔で配置された複数の検出位置において受信し、電気信号に変換する素子と、
前記電気信号を用いて前記被検体内部の特性情報を示す画像を構成する構成手段と、
前記画像を周波数空間の信号に変換したのち、前記複数の検出位置間の距離ｄと、前記検出位置と前記画像上の各位置を結ぶ線が前記平面から垂直な方向となす角度θを用いて、前記周波数空間の信号から前記画像におけるストリークアーチファクトに相当する空間周波数成分を低減するアーチファクト低減手段と、
を有することを特徴とする被検体情報取得装置。
【請求項２】
前記アーチファクト低減手段は、前記距離ｄと前記角度θに基づき算出した空間周波数フィルタを用いて、前記ストリークアーチファクトに相当する空間周波数成分を低減することを特徴とする請求項１に記載の被検体情報取得装置。
【請求項３】
前記構成手段は、前記被検体内部の吸収体に対応する特性情報を構成し、
前記アーチファクト低減手段は、前記平面に対して平行な方向および垂直な方向のそれぞれにおける、前記検出位置から前記吸収体の予想位置までの距離を用いて、前記空間周波数フィルタを算出する
ことを特徴とする請求項２に記載の被検体情報取得装置。
【請求項４】
前記アーチファクト低減手段は、
前記平面に対して平行な方向における前記検出位置から前記ストリークアーチファクトまでの距離をｘとしたときに、下式を用いて、当該平行な方向における前記ストリークアーチファクトの周期を算出し、
【数１】
2015062551.tif 000009
前記平面に対して垂直な方向における前記検出位置から前記ストリークアーチファクトまでの距離をｚとしたときに、下式を用いて、当該垂直な方向における前記ストリークアーチファクトの周期を算出する
【数２】
2015062551.tif 000010
ことを特徴とする請求項１に記載の被検体情報取得装置。
【請求項５】
前記アーチファクト低減手段により前記ストリークアーチファクトに相当する空間周波数成分を低減されたのち、実空間に変換された実空間画像を表示する表示手段をさらに有する
ことを特徴とする請求項１ないし４のいずれか１項に記載の被検体情報取得装置。
【請求項６】
前記表示手段は、前記実空間画像をＭＩＰ画像、スライス画像または３Ｄ画像により表示する
ことを特徴とする請求項５に記載の被検体情報取得装置。
【請求項７】
前記アーチファクト低減手段は、前記平面に対して垂直な方向における距離に対応した
空間周波数フィルタを算出する
ことを特徴とする請求項２または３に記載の被検体情報取得装置。
【請求項８】
前記アーチファクト低減手段は、前記空間周波数フィルタの窓関数として、ガウス窓、ハミング窓またはハン窓を使用する
ことを特徴とする請求項２または３に記載の被検体情報取得装置。
【請求項９】
操作者が前記空間周波数フィルタを変更するためのユーザーインタフェースを有する
ことを特徴とする請求項２または３に記載の被検体情報取得装置。
【請求項１０】
前記ユーザーインタフェースは、前記空間周波数フィルタに用いる周波数、ならびに、窓関数の種類および大きさの変更を可能とする
ことを特徴とする請求項９に記載の被検体情報取得装置。
【請求項１１】
前記平面上に複数の前記素子が配列された探触子をさらに有する
ことを特徴とする請求項１ないし１０のいずれか１項に記載の被検体情報取得装置。
【請求項１２】
前記探触子を前記被検体上で走査させる駆動手段をさらに有する
ことを特徴とする請求項１１に記載の被検体情報取得装置。
【請求項１３】
素子と、構成手段と、アーチファクト低減手段と、を有する被検体情報取得装置の制御方法であって、
前記素子が、被検体から伝搬する音響波を、有限の面積を持つ平面上に所定の間隔で配置された複数の検出位置において受信し、電気信号に変換する変換ステップと、
前記構成手段が、前記電気信号を用いて前記被検体内部の特性情報を示す画像を構成する構成ステップと、
前記アーチファクト低減手段が、前記画像を周波数空間の信号に変換したのち、前記複数の検出位置間の距離ｄと、前記検出位置と前記画像上の各位置を結ぶ線が前記平面から垂直な方向となす角度θを用いて、前記周波数空間の信号から前記画像におけるストリークアーチファクトに相当する空間周波数成分を低減する低減ステップと、
を有することを特徴とする被検体情報取得装置の制御方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体情報取得装置およびその制御方法に関する。
【背景技術】
【０００２】
音響波を利用して生体等の被検体の内部を画像化する装置として、例えば医療診断に用いられる光音響トモグラフィー（ＰＡＴ：Ｐｈｏｔｏａｃｏｕｓｔｉｃ  ｔｏｍｏｇｒａｐｈｙ）装置が提案されている。光音響トモグラフィー装置は、レーザパルス光を被検体に照射し、被検体内の組織が照射光のエネルギーを吸収した結果生じる音響波を探触子で受信する。この受信した信号を用いて再構成処理を行うことで、被検体内部の光学特性値に関連した画像データが得られる。
光音響トモグラフィー装置の例として、特許文献１には、被検体を音響整合材によって挟み画像化する構成が開示されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１０−１６７００１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
光音響トモグラフィー装置は、被検体からの信号を広い開口上で取得することにより、ノイズに強く、被検体内の吸収体形状に依存しない、良好な画像を取得することができる。
しかし、特許文献１に記載の光音響トモグラフィー装置は音響整合材を用いて被検体を保持しているために、被検体から発生する信号のうち、ある一定角度からの信号のみしか受信できない。そのため、開口が有限な平面となっており、得られる信号が開口制限を受けるような構成となっている。開口制限を受けた信号を用いて再構成された画像には、再構成アーチファクトが多く含まれ、被検体内部の特性情報の診断学的価値を損なう可能性がある。
【０００５】
上記のような装置構成を変更せず再構成アーチファクトを低減させるために、モデルベース法に代表される制約条件付きの再構成方法が有効であり、研究されている。しかし、制約条件付きの再構成方法は多回数の繰り返し演算が必要となり、画像データ生成に要する時間が長くなる。
【０００６】
本発明は上記課題に鑑みてなされたものであり、画像データ生成に要する時間を抑制しつつ、画像中のアーチファクトを低減する方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
本発明は、以下の構成を採用する。すなわち、
被検体から伝搬する音響波を、有限の面積を持つ平面上に所定の間隔で配置された複数の検出位置において受信し、電気信号に変換する素子と、
前記電気信号を用いて前記被検体内部の特性情報を示す画像を構成する構成手段と、
前記画像を周波数空間の信号に変換したのち、前記複数の検出位置間の距離ｄと、前記検出位置と前記画像上の各位置を結ぶ線が前記平面から垂直な方向となす角度θを用いて、前記周波数空間の信号から前記画像におけるストリークアーチファクトに相当する空間周波数成分を低減するアーチファクト低減手段と、
を有することを特徴とする被検体情報取得装置である。
【０００８】
本発明はまた、以下の構成を採用する。すなわち、
素子と、構成手段と、アーチファクト低減手段と、を有する被検体情報取得装置の制御方法であって、
前記素子が、被検体から伝搬する音響波を、有限の面積を持つ平面上に所定の間隔で配置された複数の検出位置において受信し、電気信号に変換する変換ステップと、
前記構成手段が、前記電気信号を用いて前記被検体内部の特性情報を示す画像を構成する構成ステップと、
前記アーチファクト低減手段が、前記画像を周波数空間の信号に変換したのち、前記複数の検出位置間の距離ｄと、前記検出位置と前記画像上の各位置を結ぶ線が前記平面から垂直な方向となす角度θを用いて、前記周波数空間の信号から前記画像におけるストリークアーチファクトに相当する空間周波数成分を低減する低減ステップと、
を有することを特徴とする被検体情報取得装置の制御方法である。
【発明の効果】
【０００９】
本発明によれば、画像データ生成に要する時間を抑制しつつ、画像中のアーチファクトを低減する方法を提供できる。
【図面の簡単な説明】
【００１０】
【図１】被検体情報取得装置の構成図
【図２】吸収体と受信素子の関係図および信号波形図
【図３】ストリークアーチファクトの概念図
【図４】ストリークアーチファクト周期算出式の記号説明図
【図５Ａ】ストリークアーチファクト周期結果を示すグラフ
【図５Ｂ】ストリークアーチファクト周期結果を示すグラフ
【図６】再構成画像、実施例の結果画像および空間周波数フィルタの例を示す図
【図７】実施例３に記載の発明を説明する作業フロー図
【図８】実施例４に記載の表示部を示す図
【図９】実施例４に記載の空間周波数フィルタ手動決定画面を示す図
【発明を実施するための形態】
【００１１】
以下に図面を参照しつつ、本発明の好適な実施の形態について説明する。ただし、以下に記載されている構成部品の寸法、材質、形状およびそれらの相対配置などは、発明が適用される装置の構成や各種条件により適宜変更されるべきものであり、この発明の範囲を以下の記載に限定する趣旨のものではない。
【００１２】
本発明において、音響波とは、音波、超音波、光音響波、光超音波と呼ばれる弾性波を含む。つまり本発明の被検体情報取得装置とは、被検体に光（電磁波）を照射することにより光音響効果に従って被検体内で発生した音響波を受信して、被検体内部の特性情報を取得する装置である。したがって本発明の被検体情報取得装置は、光音響トモグラフィー装置と呼ぶこともできる。
【００１３】
このとき取得される被検体内部の特性情報とは、光照射によって生じた音響波の初期音圧、初期音圧から導かれる光エネルギー吸収密度や吸収係数、組織を構成する物質の濃度等を反映した被検体情報である。物質の濃度とは例えば、酸素飽和度またはオキシヘモグロビン濃度もしくはデオキシヘモグロビン濃度である。また、特性情報としては、数値データではなく、被検体内の各位置の分布情報を取得しても良い。つまり、吸収係数分布や酸素飽和度分布等の分布情報を画像データとして取得しても良い。
【００１４】
以下、図面を参照しつつ、本発明を詳細に説明する。なお、同一の構成要素には原則と
して同一の符号を付して、説明を省略する。本発明は被検体情報取得装置やその作動方法、制御方法としても捉えられる。本発明はまた、制御方法を情報処理装置等に実施させるプログラムとしても捉えられる。
【００１５】
＜基本的な実施形態＞
まず、図１を参照しつつ、被検体情報取得装置の構成と、平面開口を持つ光音響トモグラフィーについて説明する。
【００１６】
（装置構成）
装置は、光源１０１、光照射部１０２、保持板１０３、音響整合材１０５、探触子１０６、探触子スキャン駆動部１１１、スキャン・信号取得制御部１１２、信号処理部１０７、再構成部１０８、アーチファクト低減部１０９、表示部１１０を備える。
【００１７】
被検体情報取得装置は、光源１０１からパルス光を発生させ、そのパルス光をバンドルファイバにより光照射部１０２まで導き、保持板１０３越しに被検体に照射する。被検体内に照射されたパルス光によって被検体１０４内の吸収体は光音響波を発生する。
この光音響波は音響整合材１０５を通り、探触子１０６に受信され、電気信号に変換される。探触子１０６は、探触子の音響波検出面と平行な面に沿って、探触子スキャン駆動部１１１によってスキャンされる。探触子スキャンと信号取得のタイミングはスキャン・信号取得制御部１１２によって制御されている。
【００１８】
信号処理部１０７は、電気信号のデジタル変換や増幅処理を行う。信号処理部１０７には不図示のメモリが含まれており、同位置で取得された信号等を記録する。信号処理部１０７は、スキャン終了時に被検体上の同位置で取得された信号を平均し出力する。再構成部１０８は、信号処理部１０７からの出力信号に遅延処理、再構成などを行い被検体内の吸収係数分布を作成する。アーチファクト低減部１０９は、再構成アーチファクトが低減された吸収係数分布を作成する。表示部１１０は吸収係数分布を表示する。
【００１９】
光源１０１は、７００ｎｍ以上の波長をもつナノ秒オーダーのパルス光を発生するものが好ましい。７００ｎｍより短い波長では、ヘモグロビンやコラーゲン等に多く吸収されるために、光が被検体内深部まで十分に到達できない可能性がある。大出力を得るためにはレーザーが好ましいが、レーザーのかわりに発光ダイオードなどを用いることも可能である。レーザーとしては、固体レーザー、ガスレーザー、色素レーザー、半導体レーザーなど様々なレーザーを使用できる。照射のタイミング、波形、強度などは不図示の光源制御部によって制御される。
【００２０】
また、光を光源から被検体に導くために、光を反射するミラーや、光を集光したり拡大したり形状を変化させるレンズ、光を分散・屈折・反射するプリズム、光を伝搬させる光ファイバ、拡散板等の光学部材を用いても良い。
光は、被検体に対して探触子と同じ側から照射してもよく、反対側から照射してもよい。さらに被検対の両側から照射してもよい。図１では、光照射部１０２を３箇所に設けている。
【００２１】
探触子１０６は、音響波を受信する素子を有する音響検出器である。素子を探触子の面内に複数並べられているタイプであれば、一度に複数の位置の信号を取得可能である。これにより受信時間を短縮できると共に、被検体の振動などの影響を低減できる。
ただし、素子が１つだけ、または少数の場合でも、素子を検出位置の間で走査させることにより、多数の素子を有する探触子と同様の効果を得られる。このような走査は、被検体上の仮想的な平面に所定の間隔で配置された複数の検出位置を想定し、当該検出位置間で素子を移動させることにより行われる。この場合、以下の記載における素子間距離は、
検出位置間の距離と読み替えられる。
【００２２】
探触子１０６に用いられる素子は、圧電現象を用いた変換素子、光の共振を用いた変換素子、容量の変化を用いた変換素子などがある。音響波を受信して電気信号に変換できるものであれば、どのようなものでもよい。
【００２３】
探触子スキャン駆動部１１１は探触子を被検体上で走査させるための駆動手段である。スキャンを行うことにより、被検体の広範囲からの光音響信号を取得可能となる。スキャン・信号取得制御部１１２は、スキャンするタイミングやスキャンピッチ、スキャン中の信号取得タイミングを制御する。信号取得タイミングは、光源制御部による発光タイミングと同期して制御される。探触子スキャン駆動部によって探触子が動くスキャン範囲や信号平均化回数などは不図示のメモリに記録されており、また、外部からの操作によって変更することもできる。探触子は、探触子の検出面と平行な平面上をスキャンされる。
【００２４】
信号処理部１０７は、探触子から入力された電気信号の増幅やアナログデジタル変換を行う。また、増幅やアナログデジタル変換の前後において、探触子スキャンによって同じ検出位置で撮られた電気信号を足し合わせ、平均化することも好ましい。これにより体動の影響除去や、素子特性のばらつきの緩和、ＳＮ比向上が期待できる。
【００２５】
画像再構成部１０８は、信号処理部１０７から入力された電気信号を用いて画像再構成を行い、被検体内の吸収係数の分布を表す画像データである、吸収係数分布画像情報を生成する。吸収係数だけでなく、初期音圧や物質濃度、酸素飽和度を生成しても良い。
本実施形態の再構成の方法として、トモグラフィー技術で用いられるユニバーサルバックプロジェクションを採用した。他の再構成方法としてフーリエドメインでの逆投影、合成開口法、タイムリバーサル法などが挙げられる。
【００２６】
上記に挙げた再構成方法によれば、受信領域が球、円柱、無限平面である場合、光学特性分布を一意に決定できる。しかし、本発明における開口は、有限の面積を持つ平面開口であるため、吸収係数分布が一意に決まらない。再構成アーチファクトのストリークアーチファクトは、光学特性係数が一意に決定しないために現れるアーチファクトである。
【００２７】
アーチファクト低減部１０９は、再構成部１０８から出力された画像に含まれるストリークアーチファクトを低減する。
【００２８】
（アーチファクトの発生原理）
ストリークアーチファクトの発生原理を以下に説明する。
図２Ａに示すような、探触子中の素子Ａ（２０２）と、素子Ａから距離ｄ離れた素子Ｂ（２０３）、素子Ａからｚ方向に距離ｒＡ離れた吸収体２０１を考える。なお、縦軸は素子の配列方向を示し、ｘ方向とも呼ぶ。横軸は被検体深部に向かうｚ方向を示す。吸収体は被検体内に存在し、被検体の音速をｖとする。
なお、探触子の各素子は所定の素子幅を持っている。そこで、これ以降、素子間の距離や素子と吸収体の距離を考える場合、素子の中心点を基準として考察を行う。
【００２９】
吸収体から発生した光音響波は、図２Ｂに示すような波形となる。このとき素子Ａが観測する波形Ａが出現する時間ｔＡ、素子Ｂが観測する波形Ｂが出現する時間ｔＢは、式（１）で表される。
ｒＡ＝ｖ・ｔＡ
ｒＢ＝ｖ・ｔＢ      …（１）
【００３０】
再構成では、式（１）を用いて時間を距離に変換し、素子から波形のある距離ｒＡおよ
びｒＢのどこかに吸収体２０１が存在すると考え、被検体内での吸収体があるだろう予想位置を算出する。
【００３１】
図３に、素子信号から吸収画像を算出する概念図を示す。素子Ａ（２０２）から式（１）で算出された距離ｒＡ離れた位置２０４が、吸収体の予想位置である。より詳細に言うと、位置２０４は、素子Ａから距離ｒＡ離れた点を結んだ線の上の位置である。また、素子Ｂ（２０３）から距離ｒＢ離れた位置２０５も、吸収体の予想位置である。
再構成では、予想位置２０４と２０５が交差する位置２０６に吸収体２０１があるとして画像を作成する。よって、素子の数が多く、素子が吸収体からみて広い開口領域にわたって配置されていると、吸収体の存在位置をより正確に特定しやすい。一方、本発明では素子が有限な面積を持つ平面上にあるため、吸収体が存在する確からしい位置に揺らぎが出て、吸収体位置２０６がぼやけてしまうことがある。
【００３２】
ストリークアーチファクトはこの揺らぎやぼやけが画像データに反映されたものであり、吸収体位置２０６以外の予想位置２０４，２０５を示す線上に出現する。ストリークアーチファクトは、吸収体から尾を引く線として出現する。
【００３３】
図４Ａに素子Ａ（２０２）および素子Ｂ（２０３）と、素子Ａと素子Ｂがそれぞれ作成するストリークアーチファクト２０４および２０５の関係を示す。素子Ａの作成するスリークアーチファクト２０４と素子Ｂの作成するストリークアーチファクト２０５を、探触子の検出面上の距離ｘ、探触子の検出面と垂直な方向の距離ｚを用いて表現すると、式（２）、式（３）となる。
ｘ２＋ｚ２＝ｒＡ２      …（２）
（ｘ−ｄ）２＋ｚ２＝ｒＢ２      …（３）
【００３４】
ここで、図４Ａに示すように、距離ｒＡは素子Ａと吸収体間の距離、距離ｒＢは素子Ｂと吸収体間の距離である。ｄは、素子間距離である。よって、式（４）となる。
ｒＢ２＝ｒＡ２＋ｄ２      …（４）
【００３５】
ここで、探触子面からｚ方向に、ある距離ｚ０だけ離れた位置での、素子Ａと素子Ｂがそれぞれ作成するストリークアーチファクトについて検討する。それぞれのストリークアーチファクトのｘ軸方向位置をｘＡ、ｘＢとすると、ｘＡとｘＢの距離差Δｘは、式（２）、式（３）に基づき、式（５）、式（６）のように求められる。
ｘＡ２＋ｚ０２＝ｒＡ２      …（５）
（ｘＢ−ｄ）２＋ｚ０２＝ｒＢ２      …（６）
ここに式（４）を代入し整理すると、式（７）となる。
【数１】
2015062551.tif 000003
【００３６】
また、「ｘＡ＝ｒＡ・ｓｉｎθ」であるので、式（７）は式（８）と変形される。角度θは、前記検出位置と前記画像上の各位置を結ぶ線が前記平面から垂直な方向となす角度である。
【数２】
2015062551.tif 000004
よって、ストリークアーチファクトは、素子間距離ｄと吸収体までの距離ｒＡ、素子か
らストリークアーチファクトまでの角度θであらわされる周期性を持つことがわかる。
【００３７】
ここで、「ｄ２／ｘＡ２≪１」の場合、「Δｘ＝ｄ＋ｄ２／２ｘＡ」なので、式（９）となる。
Δｘ＝ｄ＋ｄ２／２ｒＡｓｉｎθ      …（９）
またここから、「ｄ２≪２ｒＡｓｉｎθ」のとき、式（１０）となる。
Δｘ＝ｄ      …（１０）
これは吸収体が素子から十分遠いときには式（９）の間隔でｘ方向にストリークアーチファクトが現れ、距離が遠くなるほど素子間距離ｄに近づいていくことを示す。
【００３８】
同様に、図４Ｂに示すような、ｘ方向のある距離ｘ０で、素子Ａと素子Ｂがそれぞれ作成するストリークアーチファクトについて検討する。それぞれのストリークアーチファクトのｚ軸方向位置をｚＡ、ｚＢとすると、ｚＡとｚＢの距離差Δｚを算出すると、式（１１）となる。
【数３】
2015062551.tif 000005
ここで、「２ｘ０ｄ／ｚＡ２≪１」の場合、式（１２）となる。
Δｚ＝ｘ０ｄ／ｚＡ      …（１２）
【００３９】
また、「ｘ０＝ｒＡ・ｃｏｓθ」、および、「ｚＡ＝ｒＡ・ｓｉｎθ」となる。
ここで、θは素子Ａからの角度であり、素子の指向角をθｌｉｍとすると、「０＜θ＜θｌｉｍ」である。すると式（１１）は、式（１１’）と表される。
【数４】
2015062551.tif 000006
また式（１２）は、式（１３）となる。
Δｚ＝ｄ・ｔａｎθ      …（１３）
【００４０】
これらの式から、ストリークアーチファクトはｚ方向に式（１１）で示される間隔で現れ、吸収体が遠くなると式（１３）の間隔で現れることがわかる。
【００４１】
なお、検出位置が並んだ平面と平行な面におけるストリークアーチファクトを表す式（７）は、一般的には、ｘ方向におけるストリークアーチファクトの出現位置として、下式のように記述できる。
【数５】
2015062551.tif 000007
また、検出位置が並んだ平面と垂直な面におけるストリークアーチファクトを表す式（１１）は、一般的には、ｚ方向におけるストリークアーチファクトの出現位置として、下式のように記述できる。
【数６】
2015062551.tif 000008
【００４２】
図５Ａと図５Ｂは、音速１５００ｍ／ｓの被検体に現れるストリークアーチファクトの周期性に関するグラフである。図５Ａは、式（７）に対応するｚ方向の周期性を示す。図５Ｂは、式（１３）に対応するx方向の周期性を示す。それぞれ縦軸は周期性（１／ｍｍ
）を示し、横軸は時間を示す。
このように、探触子の検出平面方向のストリークアーチファクトは主に素子幅ｄの周期を持ち、検出面に垂直方向には角度θに依存した周期を持つ。
【００４３】
（アーチファクトの低減）
アーチファクト低減部は、探触子検出面と平行な方向と垂直な方向に沿って画像を周波数空間に変換し、空間周波数フィルタをかけてアーチファクトを低減する。空間周波数フィルタは、式（７）および式（１３）で示されたストリークアーチファクトの周期を低減するフィルタである。
【００４４】
ここで、ストリークアーチファクトに相当する空間周波数成分の求め方について述べる。複数の検出位置の配列方向（ｘ方向）と平行方向に周期的に表れるストリークアーチファクトに相当する空間周波数ｆｘは、１／Δｘとなる。なお、Δｘは式（７）〜（１１）に示すように求められる。素子配列が２次元アレイ状の場合、y方向に周期的に表れるス
トリークアーチファクトに相当する空間周波数についても、同様の方法で求められる。
また、ｘ方向と垂直なｚ方向に周期的に表れるストリークアーチファクトに相当する空間周波数ｆｚは、１／Δｚとなる。なお、Δｚは式（１１）〜式（１３）に示すように求められる。
【００４５】
したがって本発明で用いるフィルタとしては、周波数空間に変換された画像から、上記のように求められた空間周波数ｆｘ、ｆｚ成分を低減させるものが用いられる。
アーチファクト低減部は、このような条件で選ばれた空間周波数フィルタを周波数空間の画像に作用させて、ストリークアーチファクトの成分を低減させる。その後、画像を実空間に逆変換することで、ストリークアーチファクトが低減された画像が得られる。
【００４６】
ストリークアーチファクトの周波数を低減できる方法であれば、前述した以外の方法をとってもかまわない。たとえばアーチファクト低減部は、空間周波数フィルタを入力画像の空間周波数に掛けたのち実空間画像に戻すことでストリークアーチファクト画像を作成し、この作成した画像を実空間の入力画像から差し引いてもよい。この方法でも、ストリークアーチファクトに対応数空間周波数成分が低減できる。
このストリークアーチファクト画像に先鋭化などの画像処理を行うことも可能である。フィルタの窓関数にはガウス窓、ハミング窓、ハン窓などが使用できる。
【００４７】
アーチファクト低減部から出力された画像は表示部に表示される。表示方法はＭＩＰ画像（ＭＩＰ：Ｍａｘｉｍｕｍ  Ｉｎｔｅｎｓｉｔｙ  Ｐｒｏｊｅｃｔｉｏｎ）や、スライス画像が好適である。ただし、他の表示方法も適用可能である。例えば、３Ｄ画像を異なる複数の方向から表示する方法がある。また、表示画像の傾きや表示領域、ウインドウレベルやウインドウ幅を、利用者が表示を確認しながら変更する方法もある。
【００４８】
また、表示部がアーチファクト低減処理を行う前後の画像を表示するようにすれば、低減効果の把握が容易となる。
装置は、アーチファクト低減部で使用されるフィルタの窓関数や幅を変更するユーザー
インタフェースを備えても良い。また、フィルタが変更された画像や低減処理前後の画像を比較することができるように画像を小ウインドウで表示することが望ましい。
【００４９】
以上述べたように、本発明の実施の形態によれば、画像中のストリークアーチファクトの低減処理を効果的に行い、良好な画像を取得可能となる。
【００５０】
＜実施例１＞
本発明における光音響診断装置の一例について説明する。
光源にはＴｉ：Ｓ（チタンサファイアレーザー）を用い、約８００ｎｍの波長をもつ光を照射した。これらの波長の光密度は１５［ｍＪ／ｃｍ２］であった。Ｔｉ：Ｓの励起にＮｄ：ＹＡＧレーザー光（波長１０６４ｎｍのナノ秒オーダーのパルス光）を用いた。被検体の深部まで届く７００ｎｍ以上の波長を用いて測定を行った。
被検体を保持板と音響整合材間に挟んで保持し、パルス光を照射した。被検体から放射された光音響波を、音響整合材を通して、探触子で受信した。
【００５１】
素子幅１ｍｍ、素子ピッチ１ｍｍ、大きさ２０ｍｍ×３０ｍｍの探触子を使用した。探触子の中心周波数は２ＭＨｚであり、指向性は４０度であった。光音響波を受信したのち、探触子スキャン駆動部によって探触子が動かされ、１５０ｍｍ×９０ｍｍの測定範囲を測定した。被検体厚は５３ｍｍであった。信号は信号処理部で増幅された後にデジタル信号に変換され、再構成処理部に入力された。
【００５２】
再構成部ではユニバーサルバックプロジェクション法を用いて、吸収係数分布画像を再構成した。本実施例では探触子と被検体間に音響整合材が挿入されており、音速やインピーダンスが異なるため、超音波信号は被検体、音響整合材、探触子間で屈折する。再構成時には屈折角を計算し画像再構成を行っている。再構成後の画像を図６Ａに示す。画像の軸はｘ、ｙ、ｚが探触子の走査方向ｘ、ｙと被検者厚み（信号の時間）方向に対応している。図の軸は縦軸がｘ方向、横軸はｚ方向を示している。図に示した画像はｙ方向において中心部２．５ｍｍのＭＩＰ値を表示している。
【００５３】
画像の探触子から遠い範囲、つまり式（９）および式（１０）が成り立つ範囲において、ストリークアーチファクト５０１の強度が大きい。これは、吸収体の収集係数が大きいためと、ｚ軸距離が大きい範囲では、より吸収体から見た信号を受ける開口が狭くなり、吸収体の存在位置が確定できにくくなるためである。
【００５４】
アーチファクト低減部に再構成後の画像を入力する。入力された画像はｚ−ｘ面上で２次元フーリエ変換され空間周波数空間に変換される。空間周波数空間に変換すれば、フーリエ変換以外の手段でもよく、線形フィルタなどを用いてもよい。また、実空間に戻す際も逆フーリエ変換、線形フィルタなどの手段を用いてもよい。
式（１１）から、ｘ方向のストリークアーチファクトの周期は１ｍｍと算出できる。式（１３）と指向角４０度から、ｚ方向のストリークアーチファクトの周期は０〜０．８３ｍｍと算出できる。
【００５５】
ｚ距離が遠い範囲の画像では３０度〜４０度のストリークアーチファクトが多くみられた。これは再構成画像のボクセルが０．２５ｍｍであるために０．２５ｍｍ以下の周期をもつストリークアーチファクトが画像として現れにくく、逆に０．５ｍｍ以上の周期のアーチファクトは現れやすいからである。
【００５６】
よって、空間周波数フィルタをｘ方向１ｍｍ、ｚ方向に０．５７ｍｍ〜０．８３ｍｍの周期に対応する空間周波数のｘ方向１／ｍｍと１．４９／ｍｍ位置に窓関数の中心をおいた。窓関数はハン窓関数を使用し大きさは３２／ｍｍ×１９／ｍｍを用いた。ハン窓の大
きさは低減周波数の範囲により変更される。再構成時には被検体、音響整合材、探触子間の屈折や音速を考慮したが、アーチファクト低減部では屈折による誤差が少ないとして計算した。しかし、屈折を考慮し、ストリークアーチファクトの周期性を算出してもかまわない。
【００５７】
作成した空間周波数フィルタを図６Ｈに示す。図６Ｈの右上が周波数空間の中心である。空間周波数フィルタを空間周波数に変換された画像に掛け、空間周波数の元の画像から差し引いてストリークアーチファクト低減処理をおこない。空間周波数画像を逆フーリエ変換することにより実空間に戻し、ストリークアーチファクト低減画像を作成した。
ストリークアーチファクト低減処理を同じ手順でｙ−ｚ面でも実行した。
ストリークアーチファクトを低減した画像と元画像との差分を図６ＢとＣに示す。ｚ方向が遠い範囲のストリークアーチファクトが大幅に減少しており、効果が確認された。
【００５８】
＜実施例２＞
実施例１においてｚ方向が遠い範囲のみを考えた空間周波数フィルタを用いてストリークアーチファクトを低減した。本実施例ではｚ方向が浅い範囲の空間周波数も考えて空間周波数フィルタを作成し、ストリークアーチファクトを低減した。ここでは、実施例１と異なる部分を特に詳しく説明する。
【００５９】
被検体を保持板と音響整合材間に挟んで保持し、パルス光を照射した。被検体から放射された光音響波を、音響整合材を通して、探触子で受信した。１５０ｍｍ×９０ｍｍの測定範囲を測定した。被検体厚は５３ｍｍであった。
再構成部から出力された画像はアーチファクト低減部に入力される。入力画像をｘ−ｚ平面で２次元フーリエ変換した。
【００６０】
空間周波数フィルタの空間周波数は、ｘ方向は式（７）、ｚ方向は式（１１）から算出した。探触子からｚ方向が１０ｍｍ以降の領域を画像化したので、式（７）におけるｒＡは１０ｍｍ以上５３ｍｍ以下の範囲となる。Δｘの範囲は２ｍｍ〜１．０ｍｍとなる。また、画像のボクセルが０．２５ｍｍであるので、Δｚが０．５０ｍｍ以上のストリークアーチファクトがより多く出現する。指向性４０度と合わせると、θの範囲は約２７度以上４０度以下である。よって、Δｚの範囲は０．５７ｍｍ〜０．８３ｍｍと算出された。
【００６１】
上記の算出した空間周波数をもとに空間周波数フィルタを作成した。空間周波数フィルタを空間周波数に変換された画像に掛け、元の画像から差し引いてストリークアーチファクト低減画像を作成した。
空間周波数フィルタを空間周波数に変換された画像に掛け、空間周波数の元の画像から差し引いてストリークアーチファクト低減処理をおこない。空間周波数画像を逆フーリエ変換することにより実空間に戻し、ストリークアーチファクト低減画像を作成した。
ストリークアーチファクト低減処理を同じ手順でｙ−ｚ面でも実行した。
【００６２】
ストリークアーチファクトを低減した画像と元画像との差分を図６ＤとＥに示す。実施例１の結果図６Ｂと比較すると画像の探触子に近い範囲でのストリークアーチファクトが減少しており、効果が確認された。
【００６３】
＜実施例３＞
実施例１、２において、空間周波数フィルタは１つの画像に１つの空間周波数フィルタが使用された。ストリークアーチファクトはｚ軸方向の距離によって周期性を変化させる。深さに対応した空間周波数フィルタを作成し、ストリークアーチファクトを低減することにより、画像の情報を不必要に削除することなくストリークアーチファクトを低減することができる。本実施例では、深さによって空間周波数フィルタを変更してストリークア
ーチファクトを低減する。ここでは、実施例１、２と異なる部分を特に詳しく説明する。
【００６４】
被検体を保持板と音響整合材間に挟んで保持し、パルス光を照射した。被検体から放射された光音響波を、音響整合材を通して、探触子で受信した。１５０ｍｍ×９０ｍｍの測定範囲を測定した。被検体厚は５３ｍｍであった。
再構成部から出力された画像はアーチファクト低減部に入力される。
【００６５】
図７を用いて、本実施例でアーチファクト低減部が行う処理フローを説明する。
ステップＳ１：入力画像をあらかじめ設定された分割数に分割する。
ステップＳ２：分割した各画像をｘ−ｚ平面で空間周波数画像に変換する。
ステップＳ３：各画像の深さごとに異なった周波数フィルタを作用させ、ストリークアーチファクトを低減する。
ステップＳ４：アーチファクトが低減された空間周波数画像を実空間画像に変換する。
ステップＳ５：ｙ−ｚ平面に対して、ステップＳ２からステップＳ４を行う。
【００６６】
本実施例では、入力画像は７．５ｍｍ分重なるように１０分割された。各画像の深さは分割画像の中心深さを周波数フィルタ作成時に使用する深さとした。各画像をｘ−ｚ平面に対し空間周波数画像に変換した。
【００６７】
空間周波数フィルタの空間周波数は、ｘ方向は式（７）、ｚ方向は式（１１）から算出された。再構成部は探触子からｚ方向が１０ｍｍ以降の領域を画像化しているため、式（７）におけるｒＡは１０ｍｍ以上５３ｍｍ以下の範囲となる。Δｘの範囲は２ｍｍ〜１．０ｍｍ、Δｚの範囲は０．５７ｍｍ〜０．８３ｍｍと算出された。Δｘを初めの画像で２〜１．２ｍｍとし、最後の画像で１．２〜１．０ｍｍとなった。Δｚはほとんど角度によって決まるため、本実施例では深さごとに変化をつけなかった。しかし、角度を考慮してΔｚを変化させてもよい。算出した周波数ごとに空間周波数フィルタを作成し、各画像の空間周波数画像に掛けた。
【００６８】
低減処理後の空間周波数画像を実空間画像へ変換し、画像を足し合わせてアーチファクト低減画像を作成した。
ストリークアーチファクト低減処理を同じ手順でｙ−ｚ面でも実行した。
【００６９】
ストリークアーチファクトを低減した画像と元画像との差分を図６ＦとＧに示す。実施例１と比較すると画像の探触子に近い範囲でのストリークアーチファクトが減少しており、効果が確認された。また、実施例２と比較すると画像の深い範囲でｘ方向の低周波数成分が削除されておらず、余分な周波数を削除しない効果が確認できた。
【００７０】
＜実施例４＞
本実施例の装置は、操作者が、表示部１１０を用いて空間周波数フィルタの周波数を変更するためのユーザーインタフェースを有する。表示部１１０は、図８のようにＯＳ上で動作するアプリケーションに対応するウインドウとして提供されても良いし、専用の表示装置に表示されても良い。
【００７１】
図８において、画像表示窓８０１には、アーチファクト低減を行う前の画像が表示される。窓関数種類選択ボタン８０２には、アーチファクト低減処理に使用する空間周波数フィルタを選択するためのラジオボタンが表示される。テキストボックス８０３には、操作者により、窓関数の大きさが記入される。テキストボックス８０４には、探触子に平行方向の面（ＰＲＯＢＥ  ＤＩＲ）の周波数と、垂直な方向（ＤＥＰＴＨ  ＤＩＲ）の周波数が記入される。
【００７２】
また、ボタン８０５は、手動による空間周波数フィルタの作成を可能とする。すなわち、窓関数種類選択ボタン８０２において手動でのフィルタ作成を選択し、ボタン８０５を押した場合、図９のような周波数空間画像が別ウインドウで開く。
【００７３】
図９の画面には、実施例１で算出された周波数に対応する位置９０１が強調されて表示される。なお、表示する周波数は実施例１のものに限らず、実施例２や３の周波数のものでもよい。図９では空間周波数位置９０１を星型ポインタで示しているが、算出された空間周波数が複数ある場合はポインタの種類や数を増やしても良い。空間周波数を領域として指定したい場合は、領域を任意に塗りつぶしても良い。
【００７４】
操作者が周波数画像上をクリックしたり、周波数領域９０２を選択したりしたのち、ボタン９０３を押すと、選択した窓関数を使用し周囲をぼかした空間周波数フィルタが作成される。作成した空間周波数フィルタは、図８のボタン８０６を押すと小ウインドウが開き、確認できる。
【００７５】
図８のボタン８０７を押すと、今まで選択および指定した情報にて周波数領域が決定し、アーチファクト低減処理が行われる。低減処理後の画像は小ウインドウが開き表示される。表示した画像はスライスやＭＩＰの有無、ＭＩＰスライス厚や表示平面が変更できるように作成した。表示はＭＩＰに限られず、平均を用いてもよい。また画像表示する平面を選択するのではなく、３断面を同時に表示してもよい。また、空間周波数フィルタを作成する方法を選択するボタンをつけてもよい。
【００７６】
以上に例を挙げて述べたように、本発明によれば、被検体の画像に生じるアーチファクトに適応した空間周波数フィルタを算出することにより、画像データ生成に要する時間を抑制しつつ、画像中のアーチファクトを低減できる。とりわけ、平面開口を持つ光音響トモグラフィー装置においても、被検体内情報を損なうことなくストリークアーチファクトを低減した再構成画像を提供可能となる。
【符号の説明】
【００７７】
１０１：光源，１０２：光照射部，１０６：探触子，１０７：信号処理部，１０８：再構成部，１０９：アーチファクト低減部
【図１】
2015062551.tif 000011
【図２】
2015062551.tif 000012
【図３】
2015062551.tif 000013
【図４】
2015062551.tif 000014
【図５Ａ】
2015062551.tif 000015
【図５Ｂ】
2015062551.tif 000016
【図６】
2015062551.tif 000017
【図７】
2015062551.tif 000018
【図８】
2015062551.tif 000019
【図９】
2015062551.tif 000020