WO2007114305A1

WO2007114305A1 - 画像処理装置、およびそれを備えた超音波撮像装置、並びに画像処理方法

Info

Publication number: WO2007114305A1
Application number: PCT/JP2007/056999
Authority: WO
Inventors: Masayuki Fukuzawa; Masayoshi Yamada
Original assignee: National University Corporation Kyoto Institute Of Technology
Priority date: 2006-03-31
Filing date: 2007-03-29
Publication date: 2007-10-11
Also published as: JPWO2007114305A1; US20100234730A1; EP2005890A2; JP5034054B2; EP2005890A9

Abstract

　画像処理部（８）は、被検体断面を超音波ビームにて走査して得られる被検体断面の各点におけるエコー強度フレームｐｘ，ｙの時系列変化ｐｘ，ｙ（ｔｉ）またはドップラ速度フレームｑｘ，ｙの時系列変化ｑｘ，ｙ（ｔｉ）を周波数解析して特定周波数成分の振幅フレームＰＡおよび特定周波数成分の位相フレームＰＰを検出する周波数解析部（１２）と、検出した特定周波数成分の振幅および特定周波数成分の位相に基づく二次元画像を生成する画像生成部（１３）と、を有する。これにより、被検体の二次元形状、被検体の周期的な動きの強弱、および被検体の周期的な動きの位相の空間的連続性を観察できる画像処理装置を提供する。

Description

画像処理装置、およびそれを備えた超音波撮像装置、並びに画像処理方法

技術分野

[0001] 本発明は、エコー強度を画像化する画像処理装置、およびそれを備えた超音波撮像装置、並びに画像処理方法に関するものである。

背景技術

[0002] 従来、生体内の組織の運動や機能を超音波ビームを用いて画像化して観察する超音波撮像装置 (超音波診断装置)が用いられて！/ヽる。

[0003] 画像ィ匕は、超音波ビームを被検体である生体組織に照射し、反射された超音波ビ一ムカも得たエコー信号を用いて行っていた。このような超音波撮像装置におけるェコー信号の画像化方法としては、超音波ビームを電子走査してエコー強度の二次元分布を静止画像化あるいは動画像化する B (Brightness)モード、および超音波ビームを固定して所定期間中のエコー強度の時系列変化を静止画像ィ匕する M (Motion) モードが用いられてきた。

[0004] 例えば心臓、肝臓、または脳などの主要組織では、医師が着目する組織動きの多くは、心拍や血流の拍動に関係しているので、周期的な組織動きの強弱分布や位相の空間的連続性を観察することが有用である。

[0005] 上記の Bモード静止画像および Bモード動画像は、組織の二次元形状、および組織動きの直接観察に用いられている。しかし、 Bモード動画像は、時々刻々と表示が更新されるため、組織動きの周期性が観察できない。また、 Bモード静止画像および Bモード動画像では、組織の二次元形状は観察できるが、組織動きの周期の強弱や位相を観察することができな、。

[0006] 一方、 Mモード静止画像は、超音波プローブに接近離散する組織動きの履歴が得られるので、例えば心臓弁の動きのパターンなど、周期的な組織動きの強弱や位相の観察に用いられている。しかし、 Mモード静止画像は、超音波ビームを走査しないため、組織の二次元形状が観察できない、という問題がある。 Mモード静止画像では組織動きの周期性は観察できるが、周期的な組織動きの強弱や位相が被検体内でどのように分布し、空間的連続性を持つのかを観察することはできな、。

[0007] つまり、 Bモードを用いると、二次元形状は観察できるが、動きの周期性を観察することができず、一方、 Mモードを用いると、動きの周期性は観察できるが、二次元形状は観察できないという問題があった。また、特許文献 1に記載の超音波診断装置は、超音波ドッブラ法という方法を用いて組織の瞬時速度を検出している。これにより、関心領域における組織の瞬時速度を表示している。しかし、特許文献 1に記載の超音波診断装置では、組織の瞬時速度に基づいているため、動きの周期性を観察することができない。

[0008] 従って、医師は Bモード画像と Mモード画像を頭の中に多数記憶し、解剖学的知識を用いて組織形状と組織動きの関係を組み立てなければならな力つた。そこで、組織の二次元形状を観察でき、かつ、動きの周期性を観察することができる画像処理が望まれていた。

[0009] これに対して、本発明者らは、周期的な組織動きの強弱を二次元画像を用いて観察するため、 Bモード動画像の複数フレーム間画像演算によって、エコー強度の特定周波数成分の強弱を抽出して画像ィヒする方法を非特許文献 1および非特許文献 2において提案している。これにより、動きの周期性および二次元形状を観察することができ、上記した Bモードおよび Mモードの問題点を解決して!/、る。

[0010] また、特許文献 2に記載の超音波診断装置は、被検体内の運動体の速度を順次求め、このそれぞれの速度を基準速度を用いて補正し、この補正された速度データを表示する装置である。これにより、血流の補正速度を 2次元的または 3次元的に表示することができ、拍動性を簡便且つ効果的に表示することができる。

[0011] また、特許文献 3に記載の超音波画像処理装置は、組織追跡イメージング法において、体動等を原因とする並進速度成分や回転速度成分を取り除いた速度分布画像を、時相毎に生成し、この速度分布画像により、組織の所定位置に関する追跡処理を行い、運動情報画像を生成することで、より信憑性の高い診断画像を提供するものである。

[0012] また、特許文献 4に記載の超音波診断装置は、 2D超音波探触子を用いた 3D画像の生成に関するものであり、心臓のような動きの速い人体臓器の正常な拍動および

Zまたは不整脈の拍動を 3次元動画像で観察できるようにすることを目的とするものである。特許文献 4に記載の超音波診断装置は、具体的には、 2D探触子で心臓を走査すると、走査時間内に心臓の収縮期と拡張期とが複数含まれることの解決にあたり、心電信号を使って断層面の位相を特定する手段を有して、る。

[0013] また、特許文献 5に記載の超音波診断画像の記録表示装置は、病状変化の様子を容易に把握することが可能な装置であり、具体的には、現在と過去とにおいて得られた超音波診断画像が切替表示され、それらの画像の動画が非同期の状態で観察される。

[0014] さらに、特許文献 6に記載の発明は、周期的動きの分布を表す画像を生成する画像生成方法および周期的動きの分布を表す画像を生成し表示することが出来る超音波診断装置を提供することを課題としており、空間に分布した複数のサンプリング点から時系列データを取得し、その時系列データから周期性を求め、サンプリング点の分布と周期性の関係を表現した画像を生成している。

特許文献 1 :日本国公開特許公報「特開平 8— 84729号公報 (公開日：平成 8年 4月 2日）」

特許文献 2 :日本国公開特許公報「特開 2003— 61958号公報 (公開日：平成 15年 3 月 4日）」

特許文献 3 :日本国公開特許公報「特開 2005— 124636号公報 (公開日：平成 17年 5月 19日）」

特許文献 4:日本国公開特許公報「特開 2002— 336255号公報 (公開日：平成 14年 11月 26日）」

特許文献 5 :日本国公開特許公報「特開平 3— 97451号公報 (公開日：平成 3年 4月 23日）」

特許文献 6 :日本国公開特許公報「特開平 8— 173417号公報 (公開日：平成 8年 7 月 9曰）

非特許文献 1： Jpn.J.Appl.Phys.Vol.34,pp.2854-2856 (平成 7年 3月 18日公表）非特許文献 2 :Jpn.J.Appl.Phys.Vol.38,pp.3385- 3387 (平成 11年 2月 5日公表）発明の開示

[0015] しかしながら、上記特許文献 6、非特許文献 1、および非特許文献 2の提案手法を用いた場合には、エコー強度の特定周波数成分の強弱を求めることができるが、特定周波数成分の強弱を求めるだけだと、それらのつながり（連続性)がわ力もない、という問題があった。

[0016] また、特許文献 4、 5には、ある断層画像全体の撮影時刻が、どの心拍位相かを同定する手段は記載されているが、周期的な動きが断層面内のどこにあるかを抽出する手段がなぐこれを抽出できないという問題があった。さらに、特許文献 2、 3にも同じく周期的な動きが断層面内のどこにあるかを抽出する手段がなぐこれを抽出できないという問題があった。

[0017] 本発明は、上記の問題点に鑑みてなされたものであり、その目的は、被検体の二次元形状、被検体の周期的な動きの強弱、および被検体の周期的な動きの位相の空間的連続性を観察できる画像処理装置、および画像処理方法を提供することにある

[0018] 本発明の画像処理装置は、上記課題を解決するために、被検体断面を超音波ビームにて走査して得られる被検体断面の各点におけるエコー強度またはドッブラ速度の時系列変化を周波数解析して特定周波数成分の振幅および特定周波数成分の位相を検出する周波数解析手段と、検出した特定周波数成分の振幅および特定周波数成分の位相に基づく二次元画像を生成する第 1の画像生成手段と、を有することを特徴としている。

[0019] また、本発明の画像処理方法は、超音波撮像装置に用いる画像処理方法であって、被検体断面を超音波ビームにて走査して得られる被検体断面の各点におけるェコ一強度またはドッブラ速度の時系列変化を周波数解析して特定周波数成分の振幅および特定周波数成分の位相を検出し、検出した特定周波数成分の振幅および特定周波数成分の位相に基づく二次元画像を生成することを特徴としてヽる。

[0020] 上記画像処理装置および画像処理方法は、超音波撮像装置に用いることができる。そして、エコー強度またはドッブラ速度の時系列変化を、自装置 (画像処理装置)の外部から得ている。このエコー強度は、被検体断面を超音波ビームにて走査することにより得ることができる。特に、本発明では、このエコー強度またはドッブラ速度の時系列変化を周波数解析して、特定周波数成分の振幅だけでなぐ特定周波数成分の位相を検出している。特定周波数成分の振幅を検出して、この特定周波数成分の振幅に基づく二次元画像を生成することにより、被検体の周期的な糸且織動きの強弱分布を求めることができる。さら〖こ、本発明では、特に、特定周波数成分の位相を検出しているので、この位相に基づいて二次元画像を生成することにより、特定周波数成分の強弱のつながりのある画像を生成することができる。従って、位相の空間的連続性のある二次元画像を得ることができる。

[0021] それゆえ、組織の動きの周期性および組織の二次元形状、および位相の空間的連続性を観察することができる。

[0022] 本発明の他の目的、特徴、および優れた点は、以下に示す記載によって十分分かるであろう。また、本発明の利点は、添付図面を参照した次の説明で明白になるであろう。

図面の簡単な説明

[0023] [図 1]本発明の実施形態を示すものであり、画像処理部の要部構成を示すブロック図である。

[図 2(a)]Bモード画像の模式図およびこの画像における各点の輝度値と時刻との関係を示すグラフである。

[図 2(b)]図 2 (a)に示すグラフに振幅を求める周波数解析を行ったグラフである。

[図 2(c)]図 2 (a)に示すグラフに位相を求める周波数解析を行ったグラフである。

[図 3]本発明の実施形態を示すものであり、超音波撮像装置の概略構成を示すプロック図である。

[図 4]本発明の実施形態を示すものであり、図 3に示す超音波撮像装置の変形例を示すブロック図である。

[図 5]本発明の他の実施形態を示すものであり、画像処理部の要部構成を示すプロック図である。

符号の説明

[0024] 5 超音波送受信部 (超音波送受信手段） 8 画像処理部（画像処理装置）

12 周波数解析部 (周波数解析手段）

13 画像生成部 (第 1の画像生成手段）

23 フーリエ変換部 (フーリエ変換手段）

25 振幅 LUT (振幅ルックアップテーブル）

26 位相 LUT (位相ルックアップテーブル）

36 画像生成部 (第 1の画像生成手段）

37 二次元画像 LUT (二次元画像ルックアップテーブル）

41 プロファイル画像生成部（第 2の画像生成手段）

43 プロファイル LUT (第 2の二次元画像ルックアップテーブル）

45 位相連続曲線検出部 (検出手段）

P エコー強度フレーム

，

q ドッブラ速度フレーム

，

p (_t ) エコー強度フレームの時間履歴 (エコー強度フレームの時系列変化），

V (t ) 瞬時オプティカルフロー速度 (エコー強度フレームの時系列変化），

PA 特定周波数成分の振幅フレーム

PP 特定周波数成分の位相フレーム

t

HB 心拍周期

発明を実施するための最良の形態

[0025] 〔実施の形態 1〕

本発明の一実施の形態について図面を用いて説明する。

[0026] 図 3は、本発明の超音波撮像装置の概略構成を示すブロック図である。超音波撮像装置は、例えば、母体内の胎児の様子などを超音波のエコーを用いて詳細に観察することができ、主に医療の現場で用いられている。

[0027] 本発明の超音波撮像装置は、同図に示すように、心電波形収集部 1、制御部 2、操作部 3、超音波探触子 4、超音波送受信部 (超音波送受信手段) 5、 Bモード処理部（ Bモード処理手段） 6、フレームメモリ 7、および画像処理部（画像処理装置） 8を有している。さら〖こ、 Bモード処理部 6は、同図に示すように、フィルタ 9、対数圧縮器 10、および検波器 11を有している。また、画像処理部 8は、同図に示すように、周波数解析部 (周波数解析手段） 12および画像生成部（二次元画像生成部;第 1の画像生成手段） 13を有している。

[0028] 操作部 3は、医師などのユーザが、画像のモードを切り替えたりするインターフエ一スである。なお、図示しないが、画像処理部 8の後段には、液晶ディスプレイなどの表示部が設けられている。

[0029] 心電波形収集部 1は、外部から心電波形を収集し、該心電波形から心拍周期および時相を割り出して（求めて）、制御部 2に送る。制御部 2は、主に、この心拍周期および時相、または、操作部 3から手動で指定された解析周波数および初期位相を周波数解析部 12に送る。画像処理部 8は、既存の超音波撮像装置における画像処理部と取り替えることができると共に、既存の超音波撮像装置の内部に新機能として組み込むこともできる。なお、本発明の超音波撮像装置には、画像処理部 8をもともと備えた超音波撮像装置も含まれることは、うまでもな、。

[0030] さらに、制御部 2は、超音波探触子 4に対して、走査方向を示す信号を送り、超音波送受信部 5に対してパルス変調信号を送り、対数圧縮器 10に対してゲインを送り、画像生成部 13に対して、画像生成モードを指定する信号を送る。

[0031] なお、制御部 2は、ゲインを、操作部 3におけるユーザの指定値を基準に決定し、画像生成モードを指定する信号を、操作部 3におけるユーザの選択により生成する。

[0032] 超音波探触子 4は、それぞれが別々の方向、例えば扇形の半径方向に放射状に超音波 (超音波ビーム)を照射するように配された図示しな、複数のトランスデューサ一アレイを有している。各トランスデューサーアレイは、超音波を被検体断面（の組織 )に対して送信し、送信した走査方向からの超音波を受信する。これらのトランスデューサーアレイは、制御部 2からの走査方向を示す信号にて電子的に ON'OFFの選択がされるようになっており、これにより、超音波探触子 4が超音波を送受信する走査方向を順次切り替えることができる。なお、それぞれのトランスデューサーアレイは、超音波探触子 4に固定されており、それぞれのトランスデューサーアレイの超音波送受信方向は、予め決まっている、それゆえ、トランスデューサーアレイが決まれば超音波送受信の走査方向が一義的に決まるようになつている。 [0033] 超音波送受信部 5は、制御部 2から受け取ったパルス変調信号に基づき、超音波探触子 4のトランスデューサーアレイが被検体断面に向けて送信する超音波 (送信超音波）をパルス変調して生成すると共に、トランスデューサーアレイが被検体力も受け取った超音波 (受信超音波)を増幅する役割を有して!/、る。

[0034] Bモード処理部 6は、図 3に示すように、超音波送受信部 5の後段に配されており、フィルタ 9、対数圧縮器 10、および検波器 11を超音波送受信部 5側からこの順に備えている。この Bモード処理部 6は、超音波探触子 4のそれぞれの走査方向における被検体断面の各点のエコー強度を超音波探触子 4から被検体断面までの距離 dの関数にて示される Bモード信号 1 (d)として表す役割を有している。なお、 Bモード信号 1 (d)は、エコー強度を距離 dの関数にて示したものであるため、エコー強度信号ともいう。

[0035] フィルタ 9は、超音波送受信部 5から受け取った超音波のうち、エコー信号のみを選択して透過させる。対数圧縮器 10は、制御部 2から入力されるゲインを用いてフィルタ 9を透過したエコー信号に対して対数圧縮を行なう。検波器 11は、対数圧縮された超音波のパルス変調成分を検出して、 Bモード信号 1 (d)を生成し、出力する役割を有している。

[0036] フレームメモリ 7は、 Bモード処理部 6から受け取った超音波の異なる走査方向毎の複数の Bモード信号 1 (d)を用いて、被検体の断面内の各点 (x、 y)におけるエコー強度フレーム p を決定し、このエコー強度フレーム p の二次元分布画像である Bモード画像 Fi (i= l, 2, · ··, n)を次々と生成し、更新する（図 2 (a)参照)。ここで、 iは、二次元分布画像の生成順序を示す変数である。なお、上記したように、それぞれのトランスデューサーアレイの走査方向は予め決まっているため、被検体断面の各点 (X 、 y)は、トランスデューサーアレイと、超音波探触子 4から被検体までの距離 dが決まれば一義的に決まる。また、フレームメモリ 7には、制御部 2から画像化範囲信号が入力される。ここで、画像ィ匕範囲信号とは、エコー強度信号 (Bモード信号 1 (d) )のうち、 Bモード画像として扇形に含める信号範囲を表している。

[0037] すなわち、超音波探触子 4が超音波を送受信可能な範囲、すなわち、 I (d)における dの範囲と、 Bモード画像として生成する扇形の半径 rは、必ずしも一致せず、 r< d と成り得る。

[0038] 次に、本発明の最重要部分である、画像処理部 8について説明する。画像処理部 8は、図 1に示すように、周波数解析部 12および画像生成部 13を有している。周波数解析部 12は、同図に示すように、複数のメモリ（複数のメモリをマルチフレームメモリともいう） 20、オプティカルフロー抽出部 21、第 1の選択機構 (外部から入力される判定信号に基づいて 2種類の信号のいずれ力を後段に送る機構) 22、およびフーリェ変換部 23を有している。この周波数解析部 12は、フレームメモリ 7にて次々と生成される Bモード画像 Fl〜Fiの全部または一部を受け取り、これら力求めたエコー強度フレーム p の時系列変化 (p (t )または瞬時オプティカルフロー速度 V (t ) ) ，を周波数解析し、制御部 2から指定された特定周波数成分の振幅フレーム PA (Perio dical Amplitude) および位相フレーム PP (Periodical Phase) を求めて、画像生成部 13へ出力する。

[0039] マルチフレームメモリ 20は、フレームメモリ 7から受け取った Bモード画像 Fl〜Fiの全部または一部を時系列に蓄積して、一定時間内でのエコー強度フレーム p の時，間履歴 P (t )を生成する。

，

[0040] オプティカルフロー抽出部 21は、マルチフレームメモリ 20からエコー強度フレーム p の時間履歴 p (t )を受け取り、例えばグラディエント法を用いてエコー強度フレーム p の時間履歴 p (t )から瞬時オプティカルフロー速度 V (t )を抽出する。

，

[0041] 第 1の選択機構 22は、マルチフレームメモリ 20とフーリエ変換部 23との接続と、ォプティカルフロー抽出部 21とフーリエ変換部 23との接続とを切り替える役割を有している。

[0042] 第 1の選択機構 22は、エコー強度フレーム p のコントラストの被検体毎の変化が，

小さい場合には、マルチフレームメモリ 20とフーリエ変換部 23とを接続する。つまり、マルチフレームメモリ 20が生成したエコー強度の時間履歴 p (t )がフーリエ変換部，

23に入力される。一方、第 1の選択機構 22は、エコー強度フレーム p のコントラスト，

が被検体によって大きく変化する場合には、オプティカルフロー抽出部 21とフーリエ変換部 23とを接続する。つまり、オプティカルフロー抽出部 21が生成したォプティカルフロー速度 V (t )がフーリエ変換部 23に入力される。

， [0043] フーリエ変換部 23は、エコー強度の時間履歴 p (t )、または、オプティカルフロー

，

速度 V (t )を受け取り、これらをフーリエ変換して、特定周波数成分の振幅および，

特定周波数成分の位相を検出する。また、オプティカルフロー速度 V (t )は、点 (X

，

、 y)におけるエコー強度の時系列変化、および、点 (x、 y)の近傍におけるエコー強度の時系列変化を用いて、例えばグラディエント法によって求めることができる。

[0044] また、特定周波数成分の振幅 PAおよび特定周波数成分の位相 PPの検出は、例えば、図 1に示すように、エコー強度の時間履歴のフーリエ変換から、制御部 2から周波数解析部 12のフーリエ変換部 23へ送られる解析周波数 ω に相当する周波数

ΗΒ

成分を抽出することによって行われる。

[0045] 次に、周波数解析部 12の動作について、図 2 (a)〜図 2 (c)を用いて説明する。なお、以下では、説明の便宜上、エコー強度の時間履歴 p (t )を用いた場合につい

，

て説明するが、オプティカルフロー速度 V (t )を用いた場合も同様である。

，

[0046] 図 2 (a)〖こ示すように、マルチフレームメモリ 20は、 n個（n枚）の Bモード画像 F1〜F nを蓄積する。さらに、マルチフレームメモリ 20は、これらの Bモード画像の座標（x、 y )におけるエコー強度の時間履歴 p (t)を求める。ここで、座標 (x、 y)は、特許請求

，

の範囲に記載の「被検体断面の各点」の座標に対応して、る。

[0047] ここでは、一例として、図 2 (a)に示すように、点 a (20、 10)、点 b (21、 10)、点 c (22 、 10)、点 d (23、 10)における、輝度値 p (t)と時間 tとの関係を示す。点 aに対応する波形を波形 (e)、点 bに対応する波形を波形 (f)、点 cに対応する波形を波形 (g)、点 dに対応する波形を波形 (h)とする。

[0048] ここで、図 2 (a)に示すように、波形） (h)は、ノイズのような波形である。つまり、点

(a) (d)では、ランダムに動いていることが分かる。これに対して、波形 (f) (g)は、一定の周期性を有している。この 1周期は、図 2 (b)に示すように、 t (ハートビート)で

HB

示され、心臓の 1柏の周期とだいたい同じ周期である。波形 (f)と波形 (g)の周期は、共に t で同じであるが、 p (t)が 0になる時刻が δ t分ずれている。

HB

[0049] 次に、フーリエ変換部 23は、波形 (e)〜(h)に対して、特定周波数成分の振幅を検出する周波数解析を行なう。この結果を図 2 (b)に示す。同図では、波形 (e)〜(h)をそれぞれ周波数解析した結果を波形 (i)〜 (1)にて示してヽる。図 2 (b)では、縦軸が周波数成分 I Ρ ( ω ) Iを示しており、横軸が角周波数 ωを示している。

[0050] 同図に示すように、周期性を持たないノイズのような波形である波形 (e) (h)を周波数解析した、波形 (i) (1)は、周波数成分が突出したピークはなぐ周波数成分の大きさがほぼ同じ波形となっている。これに対して、一定の周期性を有している波形 (f) (g )を周波数解析した、波形 (j) (k)は、角周波数が ω のときにピーク (周波数成分が

ΗΒ

突出した部分)を有している。ここで、ピークの周波数成分の値をそれぞれ、振幅フレーム ΡΑ 、振幅フレーム ΡΑ としている。また、 ω は、超音波撮像装置を用

21、 10 22、 10 ΗΒ

いる医師が適当に決めることができる。また、 ω は、心電図から同期をとる心電同

ΗΒ

期にて決めてもよい。

[0051] さらに、フーリエ変換部 23は、波形 (e)〜(h)に対して、特定周波数成分の位相を検出する周波数解析を行なう。この結果を図 2 (c)に示す。同図では、波形 (e)〜(h) をそれぞれ周波数解析し、波形 (f)を位相基準とする位相を検出した結果を波形 (m ；)〜 (p)にて示している。図 2 (c)では、縦軸が位相 argPを示しており、横軸が角周波数 ωを示している。同図に示すように、周期性を持たないノイズのような波形である波形 (e) (h)に対応する、波形 (m) (p)は、 ω での位相の値は不定となるが、波形 (f)

HB

に対応する波形 (n)では、 ω での位相の値が 0となる。一方、波形 (g)に対応する

HB

波形（o)は、 ω での位相の値が非ゼロ値 ΡΡ を有している。

ΗΒ 22、 10

[0052] このような非ゼロ値 ΡΡ は、図 2 (a) (b)に示すように、波形 (f) (g)が輝度値 P (t)

22、 10

が一定の周期性を有し、かつ、波形 (f)と波形 (g)との位相のずれ δ tが発生しているために生じている。なお、上記の非ゼロ値 PP は、 ± 180° の範囲内で周期的に

22、 10

変化するが、 S t^t の 1Z2以下（以内）であれば、 δ tの値に比例して大きくなる。

HB

[0053] 画像生成部 13は、図 1に示すように、振幅 LUT (ルックアップテーブル） 25、位相 L UT26、振幅判定部 27、第 2の選択機構 28、第 3の選択機構 29、第 4の選択機構 3 0、および D/A変翻 31を有している。

[0054] そして、フレームメモリ 7からエコー強度フレーム pを受け取ると共に、周波数解析部 12から特定周波数の振幅フレーム PAおよび位相フレーム PPを受け取り、これらの V、ずれかの値またはこれらの間の演算結果を評価値として、グレースケール階調画像またはカラーグラデーション画像を生成する。生成されたこれらの画像は、後段のディスプレイに送られる。より具体的には、振幅フレーム PAは振幅 LUT25および振幅判定部 27に入力され、位相フレーム PPは位相 LUT26に入力され、エコー強度フレーム pは、第 2の選択機構 28および第 3の選択機構 29に入力される。

[0055] 振幅 LUT25は、振幅フレーム PAの値に応じたカラーグラデーション値を格納している。制御部 2は、このカラーグラデーション値に応じて、振幅フレーム PAの値をカラ一グラデーション値に変換する。一方、位相 LUT26は、位相フレーム PPの値に応じたカラーグラデーション値を格納している。制御部 2は、このカラーグラデーション値に応じて、位相フレーム PPの値をカラーグラデーション値に変換する。なお、振幅 L UT25、および位相 LUT26としては、 RAM、 ROM, EEPROMなど振幅フレーム P Aまたは位相フレーム PPの値に応じたカラーグラデーション値を格納できるものであれば、特に限定されない。

[0056] 振幅判定部 27には、制御部 2から閾値 Tが送られ、振幅判定部 27はこの閾値 T

H H

と振幅 (振幅フレーム) PAとを比較し、これらの大小に応じた切替信号を第 2の選択機構 28および第 3の選択機構 29に送り、これらの選択機構 28 · 29の後段へ送る信号をフレーム上の各画素毎に切り替える。より具体的には、振幅フレーム ΡΑが閾値 Τ未満の場合には、エコー強度フレーム ρが第 2の選択機構 28および第 3の選択機

Η

構 29からディスプレイ (後段)側へ送られ、振幅フレーム ΡΑおよび位相フレーム ΡΡ は遮断される。一方、振幅フレーム ΡΑが閾値 Τ以上の場合には、第 2の選択機構 2

Η

8および第 3の選択機構 29から振幅フレーム ΡΑおよび位相フレーム ΡΡがディスプレィ (後段)側に配された第 4の選択機構 30へ送られ、エコー強度フレーム ρは遮断される。なお、エコー強度フレーム ρのフレームレートが最も高ぐ振幅フレーム ΡΑおよび位相フレーム ΡΡのフレームレートがエコー強度フレーム ρのフレームレートよりも低くなることがある。つまり、エコー強度フレーム ρ、振幅フレーム ΡΑ、および位相フレーム ΡΡのフレームレートは必ずしも一致するとは限らず、異なる場合もある。

[0057] 第 2の選択機構 28は、振幅 LUT25の後段に配されており、振幅判定部 27からの切替信号によって、この第 2の選択機構 28からディスプレイ側へ送られる信号を、振幅フレーム ΡΑとエコー強度フレーム ρとで切り替えられるようになつている。一方、第 3の選択機構 29は、位相 LUT26の後段に配されており、振幅判定部 27からの切替信号によって、この第 3の選択機構 29からディスプレイ側へ送られる信号を位相フレーム PPとエコー強度フレーム pとで切り替えられるようになって、る。

[0058] 第 4の選択機構 30は、第 2の選択機構 28および第 3の選択機構 29の後段に配されている。そして、第 2の選択機構 28からの出力（各点 (x、 y)におけるエコー強度フレーム pまたは振幅フレーム PAの値）、または、第 3の選択機構 29からの出力（各点（ x、 y)におけるエコー強度フレーム pまたは位相フレーム PPの値）のいずれをデイスプレイ側へ出力するかを切り替える。この切り替えは、制御部 2から第 4の選択機構 3 0に入力される画像生成モード信号に基づいて行われる。

さらに、第 4の選択機構 30の後段に配された DZA変にて、デジタルデータ力もアナログデータへと変換され、その後、ディスプレイ側へ出力される。

[0059] ところで、図 3に示す上記の Bモード処理部 6は、 Bモード信号 I (d)を次々に生成して、フレームメモリ 7を常に最新のエコー強度フレーム pに更新する。マルチフレームメモリ 20は、フレームメモリ 7が更新されるたびに最も古いエコー強度フレーム pを捨てて、新しいエコー強度フレーム pのデータを蓄積することにより、最も新しい n枚のェコー強度フレーム Pを常に保持する。

[0060] 周波数解析部 12は、マルチフレームメモリ 20中の解析開始時点における最新の m 枚 (m≤n)のエコー強度フレーム pから特定周波数成分の振幅フレーム PAと位相フレーム PPとを生成する。ここで、 mを周波数解析部 12の解析フレーム数という。

[0061] なお、マルチフレームメモリ 20におけるエコー強度フレーム pの保持枚数 nを、周波数解析部 12における解析フレーム数 mよりも十分大きくすれば、振幅フレーム PAおよび位相フレーム PPを 1組生成するのに要する時間は、フレームメモリ 7が更新される時間と必ずしも一致する必要はな、。

[0062] 画像生成部 13は、エコー強度フレーム p、特定周波数成分の振幅フレーム PA、または特定周波数成分の位相フレーム PPのいずれかが更新されるたびに表示フレームを再生成し、繰り返し表示する。これにより、超音波撮像装置を用いる医師に動画像を提供することができる。

[0063] 以上により、例えば液晶ディスプレイなどのディスプレイに表示される画像は、特定周波数成分の振幅 (周期的な組織動きの強弱 ·強度) PAが閾値 T未満の画素は、エコー強度 Pのモノクログラデーション値となる、つまり、 Bモード画像となる。一方、特定周波数成分の振幅 (周期的な組織動きの強弱 ·強度) PAが閾値 T

H以上である画素は、位相フレーム PPのカラーグラデーション値となる。従って、一定以上の周期的な組織動き強度を示す領域のみが色づけされ、組織間の位相関係が色のグラデーシヨンとして表現される。それゆえ、周期的な組織動きの強弱だけでなぐ位相の空間的連続性を二次元画像として表現することができる。

[0064] 特に、本発明では、特定周波数成分の位相フレーム PPを求めて、色分けを行っている。そのため、例えば拍動組織などの血流のある組織の連続性 (つながり）を観ることができると!、う効果を奏する。

[0065] なお、上記は、特定周波数成分の振幅、および特定周波数成分の位相をフーリエ変換にて行っている力これは単なる一例にすぎず、例えば自己回帰モデルにて求めてもよい。

[0066] また、特定周波数成分としては、例えば、心拍周期に相当する周波数成分が挙げられるがこれに限定されず、心拍周期の 2倍波や 3倍波などの整数倍の第 2次高調波や第 3次高調波を用いた周波数成分も含まれる。つまり、特定周波数成分には、心拍周期以外の周波数成分も含まれる。

[0067] また、本発明の画像処理装置および超音波撮像装置は、新生児の頭部などの拍動組織 (周期的に動く組織）に特に好適に利用することができる。

[0068] また、上記図 1に示した画像生成部 13の代わりに、図 5に示すような画像生成部（第 1の画像生成手段） 36を用いてもよい。この画像生成部 36は、図 5に示すように、二次元画像 LUT (第 1の二次元画像ルックアップテーブル） 37、振幅判定部 27、および二次元画像合成部 38を備えてヽる。

[0069] 二次元画像 LUT37は、周波数解析部 12から振幅フレーム PAおよび位相フレーム PPを受け取ると共に、制御部 2 (図 3)から二次元画像生成モードを示す信号を受け取る。

[0070] 二次元画像 LUT37には、振幅フレーム PAの値、位相フレーム PPの値、またはこれらの値の組み合わせに応じたカラーグラデーション値がそれぞれ格納されており、これらのカラーグラデーション値は、二次元画像生成モードに応じて書き換えられる。ここで、振幅フレーム PAの値および位相フレーム PPの値の組み合わせによるカラーグラデーション値の例としては、振幅を色の明度で表し、位相を色の色相で表した場合のカラーグラデーション値が考えられる。

[0071] 二次元画像合成部 38は、周波数解析部 12からエコー強度フレーム pを受け取ると共に、振幅判定部 27から切替信号を受け取る。さらに、この切替信号に基づいて、二次元画像 LUT37から上記のグラデーション値を取り出す。

[0072] そして、切替信号に基づ!/、て、二次元画像の各画素毎に、エコー強度フレーム p、または、カラーグラデーション値のいずれかを選択し、二次元画像を合成し、出力する。この二次元画像合成部 38では、振幅フレーム PA、位相フレーム PP、またはェコ一強度フレーム Pのヽずれかが更新されるたびに二次元画像を再合成する。なお、振幅フレーム PA、位相フレーム PP、またはエコー強度フレーム pの更新周期は必ずしも同一である必要はない。

[0073] 〔実施の形態 2〕

さらに、本発明の他の実施の形態の超音波撮像装置について、図 4を用いて説明する。なお、上記の実施の形態との相違点のみ説明し、同一の構成および機能については、その説明を省略する。

[0074] 本実施の形態の超音波撮像装置は、上記の図 3に示す実施の形態 1の構成に加えて、ドッブラ処理部 35、ドッブラ速度フレームメモリ 49、および選択機構 50を備えている。具体的には、図 4に示すように、超音波送受信部 5からの出力が Bモード処理部 6だけでなぐドッブラ処理部 35にも入力されており、このドッブラ処理部 35の後段にドッブラ速度フレームメモリ 49が配されている。つまり、 Bモード処理部 6およびフレームメモリ 7と並列にドッブラ処理部 35およびドッブラ速度フレームメモリ 49が設けられている。

[0075] 選択機構 50は、フレームメモリ 7およびドッブラ速度フレームメモリ 49の後段、かつ、周波数解析部 12の前段に配されており、フレームメモリ 7の出力信号およびドッブラ速度フレームメモリ 49の出力信号のうち、いずれか一方を周波数解析部 12に入力するようになつている。つまり、選択機構 50は、フレームメモリ 7と周波数解析部 12との接続と、ドッブラ速度フレームメモリ 49と周波数解析部 12との接続とを切り替え可能としている。

[0076] ドッブラ処理部 35は、超音波送受信部 5から受け取った超音波エコー信号からドッブラ偏移成分を検出することによって、被検体内の組織の移動速度に対応したドッブラ速度信号 vd (d)を生成する。

[0077] ドッブラ速度フレームメモリ 49は、ドッブラ処理部力も受け取ったドッブラ速度信号 V d (d)を用いて、被検体の断面内の各点 (x、 y)におけるドッブラ速度フレーム q を，決定する。第 5の選択機構 50は、制御部 2からの制御に基づき、被検体組織が不均一で、糸且織動きがエコー強度のコントラストに反映されにくい場合には、ドッブラ速度フレームメモリ 49と周波数解析部 12とを接続する。なお、第 5の選択機構 50がフレームメモリ 7と周波数解析部 12とを接続した場合には、実施の形態 1と同じであるため、その説明を省略する。

[0078] 画像処理部 8の構成は、実施の形態 1と同様であり、異なるのは、周波数解析部 12 の入力として、新たにドッブラ速度フレーム q が加わっている点である。

，

[0079] エコー強度フレーム p は被検体の超音波の反射率を反映するのに対し、ドッブラ，

速度フレーム q は被検体組織の移動速度を反映する。それゆえ、ドッブラ速度フレ，

ームメモリ 49と周波数解析部 12とが接続された場合の周波数解析部 12の出力は、被検体組織の移動速度の周期的変化に基づく。従って、本実施の形態においても、被検体の二次元形状、被検体の周期的な動きの強弱、および被検体の動きの位相の空間的連続性を観察することができる。

[0080] なお、本実施の形態にぉ、て、第 5の選択機構 50は、必ずしも必須の構成ではない。つまり、第 5の選択機構 50を取り除き、フレームメモリ 7と周波数解析部 12とを接続せずに、ドッブラ速度フレームメモリ 49と周波数解析部 12とを接続してもよい。つまり、周波数解析部 12にドッブラ速度フレーム q のみを入力する構成としてもよい。

，

[0081] 〔実施の形態 3〕

さらに、本発明の他の実施の形態の超音波撮像装置について、図 5を用いて説明する。なお、上記の実施の形態と同一の構成'機能を有する部材についてはその説明を省略する。

[0082] 本実施の形態の超音波撮像装置は、図 5に示すように、実施の形態 1 · 2の構成に加えて、周波数解析部 12の後段に画像生成部 36と並列に配されたプロファイル画像生成部（第 2の画像生成手段) 41と、これら画像生成部 36およびプロファイル画像生成部 41の後段に配された表示画像合成部 42と、を備えている。

[0083] なお、画像生成部 36の代わりに画像生成部 13を用いてもよい。なお、周波数解析部 12および画像生成部 36については、既に説明した内容と同一であるため、その説明を省略する。

[0084] プロファイル画像生成部 41は、プロファイル LUT (第 2の二次元画像ルックアップテ一ブル) 43、プロファイル入力切替部 (入力切替部) 44、位相連続曲線検出部（曲線検出部;検出手段) 45、マルチプロファイルメモリ（プロファイルメモリ) 46、プロフアイル曲線保持部（曲線保持部) 47、およびプロファイル履歴画像合成部 (履歴画像合成部) 48を備えている。

[0085] プロファイル LUT43は、周波数解析部 12から振幅フレーム PAおよび位相フレーム PPを受け取ると共に、制御部 2 (図 3)力もプロファイル履歴画像生成モード信号を受け取る。

[0086] このプロファイル LUT43には、上記の二次元画像 LUT37と同様に、振幅フレーム PAの値、位相フレーム PPの値、またはこれらの値の組み合わせに応じたカラーダラデーシヨン値がそれぞれ格納されており、二次元画像 LUT37と同様の機能を有する。但し、プロファイル LUT43は、図 5に示すように、二次元画像 LUT37とは独立して配されており、カラーグラデーション値の書き換えは、二次元画像生成モード信号とは異なる信号である、プロファイル履歴画像生成モード信号に応じて行なわれる。

[0087] 入力切替部 44は、フレームメモリ 7からエコー強度フレーム pを受け取り、制御部 2 力プロファイル履歴画像生成モード信号を受け取り、プロファイル LUT43からカラ一グラデーション値を受け取る。そして、この入力切替部 44は、プロファイル履歴画像生成モード信号に基づいて、エコー強度フレーム pとカラーグラデーション値によるフレームとのうちのいずれかのフレームを後段のマルチプロファイルメモリ 46へ送る。

[0088] 曲線検出部 45は、振幅フレーム PAおよび位相フレーム PPを受け取り、これら振幅フレーム PAおよび位相フレーム PPから、位相が連続する曲線（の座標群）を自動的に抽出する。なお、位相の連続性に加えて、振幅が一定以上の領域を検出することが有効な場合もあり、この場合には、図 5に示すように、振幅フレーム PAの代わりに、振幅判定部 27からの切替信号を入力してもよい。

[0089] 曲線保持部 47は、曲線検出部 45から受け取った位相が連続する曲線 (の座標群）

、または、操作者が手動で指定した曲線 (の座標群)を保持し、これらの曲線をプロフアイルメモリ 46に通知する。

[0090] なお、操作者による曲線の指定は、診察上の関心領域に基づいて行なわれる。例えば、脳内出血児で出血塊付近に動脈がある場合、 Bモード画像上で出血塊の周囲に曲線を描き、その線上に沿って動きを観察するなどにより行なうことができる。なお、実際の曲線の指定は、操作者が Bモード画像上にタツチペンで描くなどの操作にて実現することができる。

[0091] プロファイルメモリ 46は、プロファイル入力切替部 44から受け取ったフレームを現在力も遡って所定のフレーム数蓄積すると共に、これらのフレームについて曲線保持部

47から通知された曲線 (の座標群）にあたる画素値群を保持する。

[0092] 履歴画像合成部 48は、プロファイルメモリ 46に格納された曲線における画素値群の時間履歴を二次元画像 (プロファイル履歴画像）に合成し、出力する。例えば、プ口ファイル曲線に沿って画素値を取り出し、それらを画像の縦 1ラインとする。

[0093] 履歴の各時刻における上記「縦 1ライン」を現在力過去に遡って右力左に並べると、縦がプロファイル曲線に沿った空間座標、横が経過時間となるような画像が生成できる。

[0094] このようにして生成した画像（出力画像)の特別な場合として、曲線として超音波ビームの方向に沿った直線を選び、プロファイル履歴画像生成モードとしてエコー強度 pのモノクログラデーションを選んだ場合、履歴画像合成部 48の出力は、いわゆる M モード画像と等価となる。

[0095] 履歴画像合成部 48は、振幅フレーム PA、位相フレーム PP、またはエコー強度フレーム pのいずれかが更新されるたびに出力画像の縦 1ラインのみを再合成する。

[0096] 表示画像合成部 42は、画像生成部 36から二次元画像を受け取ると共に、プロファィル画像生成部 41からプロファイル履歴画像を受け取る。そして、制御部 2から受け取った表示合成モード信号に基づき、 V、ずれか一方の画像をディスプレイへ出力する。

[0097] すなわち、上記構成では、 Mモード画像 (Mモードエコー画像)から動きを抽出するのではなぐ Bモード動画像力抽出した動きの周期性（二次元分布）の中から 1次元プロファイルを抜き出し、これを縦 1ラインとして描き、該縦 1ラインを時間の経過と共に、右側へ描き、これを繰り返すことにより、動きの周期性と位相の履歴を画像化することができる。

[0098] また、 Mモードにおける 1次元プロファイルは超音波の照射方向に限られている力本実施の形態では、 1次元プロファイルの選び方は、超音波ビームの選び方は超音波ビームの照射方向に限定されない。また、画像化する曲線は、上記したように、曲線検出部 45にて位相が連続する曲線 (の座標群）を自動的に抽出することによって求めるか、または、操作者が外部力も手動で指定する。

[0099] また、上記したように、入力切替部 44にて、エコー強度フレーム p、もしくは、特定周波数成分の振幅フレーム PAの値、位相フレーム PPの値、またはこれらの値の組み合わせに応じたカラーグラデーション値によるフレームのいずれかを後段のプロファィノレメモリ 46【こ送って!/ヽる。

[0100] まず、特定周波数成分の振幅フレーム PAの値、位相フレーム PPの値、またはこれらの値の組み合わせに応じたカラーグラデーション値によるフレームをプロファイルメモリ 46へ送る場合にっヽて説明する。

[0101] この場合、例えば、脳室内出血において出血塊周囲の血流阻害を観察する際に、出血塊周囲の関心領域に医師が曲線を描、て、この曲線上での組織動きの履歴を表示すれば、動画像の観察に比べて、より直感的に組織動きを把握することができる。従って、出血塊周囲などの局所的関心領域の診断に役立つ。さらに、位相の連続した曲線を抽出して曲線を設定する場合には、特に、 Bモードの直接観察では関心領域の同定が困難な際でも、動脈周囲組織の探索を支援することができる。従って、医師による局所関心領域の探索を支援することができる。

[0102] 本発明は上述した各実施形態に限定されるものではなぐ請求項に示した範囲で種々の変更が可能である。すなわち、請求項に示した範囲で適宜変更した技術的手段を組み合わせて得られる実施形態についても本発明の技術的範囲に含まれる。 [0103] また、本発明の画像処理装置では、第 1の画像生成手段は、特定周波数成分の振幅の値と特定周波数成分の位相の値との組み合わせに応じて二次元画像の色のグラデーシヨン値が割り当てられた第 1の二次元画像ルックアップテーブルを有してヽることが好ましい。上記構成によれば、二次元画像を生成する際に、視覚的に周期的な組織動きの強弱および位相の空間的連続性を観察することができる。

[0104] また、本発明の画像処理装置では、第 1の二次元画像ルックアップテーブルには、特定周波数成分の振幅の値に応じた二次元画像の色のグラデーション値および特定周波数成分の位相の値に応じた二次元画像の色のグラデーション値がさらに割り当てられていることが好ましい。上記構成によれば、二次元画像を生成する際に、視覚的に周期的な組織動きの強弱および位相の空間的連続性を観察することができる

[0105] また、本発明の画像処理装置では、画像生成手段は特定周波数成分の振幅の値に応じて二次元画像の色のグラデーション値が割り当てられた振幅ルックアップテーブルと、位相の値に応じて二次元画像の色のグラデーション値が割り当てられた位相ルックアップテーブルと、を有して!/、ることが好まし！/、。

[0106] 上記構成によれば、画像解析手段にて求めた振幅および位相の値に応じて、これらの値を二次元画像の色のグラデーション値に変換することができる。それゆえ、二次元画像を生成する際に、視覚的に周期的な組織動きの強弱および位相の空間的連続性を観察することができる。

[0107] また、本発明の画像処理装置では、特定周波数成分の振幅に基づき生成した二次元画像と特定周波数成分の位相に基づき生成した二次元画像とを切り替えて出力することが好ましい。上記構成によれば、特定周波数成分の振幅を用いた二次元画像と特定周波数成分の位相を用いた二次元画像とを切り替えて出力することができるので、ユーザが所望の映像を得ることができる。

[0108] また、本発明の画像処理装置では、第 1の画像生成手段は、生成する二次元画像を繰り返し更新することが好ましい。上記構成によれば、更新のたびに新しい二次元画像を表示することができる。

[0109] また、本発明の画像処理装置では、特定周波数成分の振幅および特定周波数成分の位相に基づいて、被検体断面上の任意の曲線における画素値の時間履歴を画像ィ匕することにより、二次元画像を求める第 2の画像生成手段と、をさらに備えていることが好ましい。

[0110] また、本発明の画像処理装置では、第 2の画像生成手段は、特定周波数成分の振幅および特定周波数成分の位相に基づ!/ヽて、被検体断面上の任意の曲線における画素値群を決定し、この画素値群の時間履歴を画像化することが好ましヽ。

[0111] 上記構成では、いわゆる Mモードエコー画像力動きを抽出するのではなくて、 B モード画像力も動きを抽出している。また、上記構成によれば、任意の曲線における、動きを画像ィ匕することができる。

[0112] また、本発明の画像処理装置では、第 2の画像生成手段は、特定周波数成分の振幅の値と特定周波数成分の位相の値との組み合わせに応じて二次元画像の色のグラデーシヨン値が割り当てられた第 2の二次元画像ルックアップテーブルを有してヽることが好ましい。上記構成によれば、二次元画像を生成する際に、視覚的に周期的な組織動きの強弱および位相の空間的連続性を観察することができる。

[0113] また、本発明の画像処理装置では、特定周波数成分の振幅および特定周波数成分の位相に基づいて位相が連続している座標群を検出することにより任意の曲線を求める検出手段を備えていることが好ましい。上記構成によれば、 Bモードの直接観察では、関心領域の同定が困難な場合でも、位相の連続した曲線を検出することができるので、例えば新生児頭部における動脈周囲組織の検索をより容易に行なうことができる。

[0114] また、本発明の画像処理装置では、任意の曲線は外部から入力可能となっていることが好ましい。上記構成によれば、操作者が曲線を入力することができるので、例えば新生児頭部における出血塊周囲などの局所的関心領域の診断を行なうことができる。

[0115] また、本発明の画像処理装置では、第 2の画像生成手段は、生成する二次元画像を繰り返し更新することが好ましい。上記構成によれば、更新のたびに新しい二次元画像を表示することができる。

[0116] また、本発明の画像処理装置では、周波数解析手段は、特定周波数成分の振幅および特定周波数成分の位相の検出をフーリエ変換にて行なうフーリエ変換手段を有していることが好ましい。

[0117] また、本発明の画像処理装置では、上記特定周波数成分は、心拍周期に相当する周波数成分であることが好ましい。上記構成によれば、特定周波数成分は心拍周期に相当する周波数成分であるので、動脈血流に伴う周期的な動きを示す組織を見つけることができる。

[0118] また、本発明の超音波撮像装置では、上記いずれの画像処理装置と、被検体の断面を超音波ビームで走査してエコー信号を繰り返し受信する超音波送受信手段とを有することが好ましい。上記構成によれば、被検体の二次元形状、被検体の周期的な動きの強弱、被検体の周期的な動きの位相の空間的連続性を観察できる超音波撮像装置を実現することができる。

[0119] 本発明の画像処理装置は、以上のように、被検体断面を超音波ビームにて走査して得られる被検体断面の各点におけるエコー強度またはドッブラ速度の時系列変化を周波数解析して特定周波数成分の振幅および特定周波数成分の位相を検出する周波数解析手段と、検出した特定周波数成分の振幅および特定周波数成分の位相に基づく二次元画像を生成する第 1の画像生成手段と、を有している。

[0120] 本発明の画像処理方法は、以上のように、超音波撮像装置に用いる画像処理方法であって、被検体断面を超音波ビームにて走査して得られる被検体断面の各点におけるエコー強度またはドッブラ速度の時系列変化を周波数解析して特定周波数成分の振幅および特定周波数成分の位相を検出し、検出した特定周波数成分の振幅および特定周波数成分の位相に基づく二次元画像を生成している。

[0121] 従って、被検体の二次元形状、被検体の周期的な動きの強弱、および被検体の周期的な動きの位相の空間的連続性を観察できる。

[0122] 発明の詳細な説明の項においてなされた具体的な実施形態または実施例は、あくまでも、本発明の技術内容を明らかにするものであって、そのような具体例にのみ限定して狭義に解釈されるべきものではなぐ本発明の精神と次に記載する請求の範囲内で、、ろ、ろと変更して実施することができるものである。

産業上の利用可能性本発明の画像処理装置は、医療現場にて用いることができ、特に、新生児の頭部の観察に好適に用いることができる。

Claims

請求の範囲

[1] 被検体断面を超音波ビームにて走査して得られる被検体断面の各点におけるェコ一強度またはドッブラ速度の時系列変化を周波数解析して特定周波数成分の振幅および特定周波数成分の位相を検出する周波数解析手段と、

検出した特定周波数成分の振幅および特定周波数成分の位相に基づく二次元画像を生成する第 1の画像生成手段と、を有することを特徴とする画像処理装置。

[2] 第 1の画像生成手段は、特定周波数成分の振幅の値と特定周波数成分の位相の値との組み合わせに応じて二次元画像の色のグラデーション値が割り当てられた第 1 の二次元画像ルックアップテーブルを有していることを特徴とする請求の範囲第 1項に記載の画像処理装置。

[3] 第 1の二次元画像ルックアップテーブルには、特定周波数成分の振幅の値に応じた二次元画像の色のグラデーション値および特定周波数成分の位相の値に応じた二次元画像の色のグラデーション値がさらに割り当てられていることを特徴とする請求の範囲第 2項に記載の画像処理装置。

[4] 第 1の画像生成手段は、特定周波数成分の振幅の値に応じて二次元画像の色のグラデーション値が割り当てられた振幅ルックアップテーブルと、特定周波数成分の位相の値に応じて二次元画像の色のグラデーション値が割り当てられた位相ルックアップテーブルと、を有して、ることを特徴とする請求の範囲第 1項に記載の画像処理装置。

[5] 特定周波数成分の振幅に基づき生成した二次元画像と特定周波数成分の位相に基づき生成した二次元画像とを切り替えて出力することを特徴とする請求の範囲第 4 項に記載の画像処理装置。

[6] 第 1の画像生成手段は、生成する二次元画像を繰り返し更新することを特徴とする請求の範囲第 1項力第 5項のいずれか 1項に記載の画像処理装置。

[7] 特定周波数成分の振幅および特定周波数成分の位相に基づいて、

被検体断面上の任意の曲線における画素値の時間履歴を画像ィ匕することにより、二次元画像を求める第 2の画像生成手段と、をさらに備えて、ることを特徴とする請求の範囲第 1項に記載の画像処理装置。

[8] 特定周波数成分の振幅および特定周波数成分の位相に基づいて、被検体断面上の任意の曲線における画素値群を決定し、この画素値群の時間履歴を画像化することにより二次元画像を求める第 2の画像生成手段と、をさらに備えていることを特徴とする請求の範囲第 1項に記載の画像処理装置。

[9] 第 2の画像生成手段は、特定周波数成分の振幅の値と特定周波数成分の位相の値との組み合わせに応じて二次元画像の色のグラデーション値が割り当てられた第 2 の二次元画像ルックアップテーブルを有していることを特徴とする請求の範囲第 7項または第 8項に記載の画像処理装置。

[10] 特定周波数成分の振幅および特定周波数成分の位相に基づ!/、て位相が連続している座標群を検出することにより任意の曲線を求める検出手段を備えていることを特徴とする請求の範囲第 7項力第 9項のいずれ力 1項に記載の画像処理装置。

[11] 任意の曲線は外部力入力可能となっていることを特徴とする請求の範囲第 7項から第 10項のいずれか 1項に記載の画像処理装置。

[12] 第 2の画像生成手段は、生成する二次元画像を繰り返し更新することを特徴とする請求の範囲第 7項力第 11項のいずれか 1項に記載の画像処理装置。

[13] 周波数解析手段は、特定周波数成分の振幅および特定周波数成分の位相の検出をフーリエ変換にて行なうフーリエ変換手段を有していることを特徴とする請求の範囲第 1項から第 12項のいずれか 1項に記載の画像処理装置。

[14] 特定周波数成分は、心拍周期に相当する周波数成分であることを特徴とする請求の範囲第 1項力も第 13項のいずれか 1項に記載の画像処理装置。

[15] 請求の範囲第 1項力第 14項のいずれか 1項に記載の画像処理装置と、

被検体の断面を超音波ビームで走査してエコー信号を繰り返し受信する超音波送受信手段とを有することを特徴とする超音波撮像装置。

[16] 超音波撮像装置に用いる画像処理方法であって、

被検体断面を超音波ビームにて走査して得られる被検体断面の各点におけるェコ一強度またはドッブラ速度の時系列変化を周波数解析して特定周波数成分の振幅および特定周波数成分の位相を検出し、検出した特定周波数成分の振幅および特定周波数成分の位相に基づく二次元画像を生成することを特徴とする画像処理方

91

6669S0/.00Zdf/X3d SOCMT/.OOZ OAV