JP5025400B2 - 超音波診断装置、及び超音波画像処理装置 - Google Patents

Description

この発明は、ドプラスペクトラム画像を取得して、ドプラスペクトラム画像の速度レンジを調整する超音波診断装置、及び超音波画像処理装置に関する。

超音波ドプラ法を用いることで、診断部位の血流情報を得る超音波診断装置が知られている。超音波ドプラ法は、超音波ドプラ法の原理に基づいて被検体内の血流の情報を得る技術である。被検体内の血流などの流れのある診断部位に向けて超音波を送信すると、ドプラ効果により、送信周波数に対して受信周波数が僅かに偏移する。この偏移周波数（ドプラ偏移周波数）は血流速度に比例し、偏移周波数（ドプラ偏移周波数）の周波数解析を行なうことにより血流情報が得られる。

超音波診断装置は、得られたドプラ信号に対して高速フーリエ変換（ＦＦＴ）を行い、その周波数分析の結果を、縦軸に周波数ｆ（速度ｖ）、横軸に時間ｔとしてスペクトラム表示を行い、そのスペクトラム画像を対象として診断で用いる各種項目（パラメータ）の計測が行われる。

ところで、パルスドプラ（ＰＷ）法では、周波数解析のためのサンプリング周波数ｆｓがドプラ偏移周波数よりも低い場合には、エイリアシング現象（折り返り）が起きてしまう。そこで、これを防ぐために、サンプリング周波数に相当するパルス繰り返し周波数（ＰＲＦ：Ｐｕｌｓｅ Ｒｅｐｅｔｉｔｉｏｎ Ｆｒｅｑｕｅｎｃｙ）を上げ、毎回の観測時間間隔を短くする必要がある。この場合、測定したい位置を指定すると、必然的に最大のＰＲＦが決定され、ＰＲＦが決まると、測定できる最高血流速度も決定されてしまう。この計測できる最高血流速度をドプラ速度レンジと称する。

ドプラスペクトラム表示を行ったところ、ドプラ速度レンジが小さすぎると、ドプラ波形の折り返し部分が発生する。このような場合、操作者が手動でドプラ速度レンジを大きく設定することで、折り返し部分がナイキスト周波数（ＰＲＦの半分）内に収まり、表示上、繋がりの良いドプラスペクトラムが得られる。一方、ドプラ速度レンジが大きすぎるとスケールが大きくなるため、ドプラ波形が小さくなってしまい、観測が困難になる。このような場合、操作者が手動でドプラ速度レンジを小さく設定することで、表示画面の上下いっぱいを有効に使ったドプラスペクトラムが得られる。また、速度オフセット（ＢＬＳ：Ｂａｓｅｌｉｎｅ Ｓｈｉｆｔ）を調整することで、ドプラ波形を速度レンジ（縦軸）の所望の位置に表示させることができる。

患者や診断部位によって血流速度が異なり、患者や診断部位によって取得されるドプラ波形が異なるため、そのドプラ波形に応じたドプラ速度レンジを設定する必要がある。そのため、操作者は、患者や診断部位が変わるたびに、ドプラ速度レンジに対応するＰＲＦと速度オフセット（ＢＬＳ）を調整する必要があった。しかしながら、患者や診断部位が変わるたびに、ドプラ速度レンジ（ＰＲＦ）や速度オフセット（ＢＬＳ）を調整することは煩雑である。

そこで、ドプラ速度レンジ（ＰＲＦ）と速度オフセット（ＢＬＳ）を自動的に決定する手法が提案されている（例えば特許文献１）。

従来においては、血流状態を判別して、その判別結果をフィードバックすることでドプラ速度レンジ（ＰＲＦ）と速度オフセット（ＢＬＳ）を求めていた。そして、リアルタイムに血流状態を判別することで、ドプラ速度レンジと速度オフセットをリアルタイムに追跡していた。

しかしながら、過去の血流状態を判別して、ドプラ速度レンジ（ＰＲＦ）と速度オフセット（ＢＬＳ）を推定するため、タイムラグが発生する問題があった。また、過去の血流状態と現在以降の血流状態とが一致するとは限らないため、現在以降に取得されるドプラ波形に適したドプラ速度レンジと速度オフセットを求めることは困難であった。

また、従来においては、血流状態が安定した時点で操作者がトリガーを与えて、それ以降、ドプラ速度レンジ（ＰＲＦ）と速度オフセット（ＢＬＳ）を推定することが行われていた。しかしながら、操作者がトリガーを与えた後に、被検体の呼吸や超音波プローブの位置ずれなどの変動が発生した場合、ドプラ波形に適したドプラ速度レンジと速度オフセットが求められず、再度、推定処理を行う必要があった。

特開２００５−１８５７３１号公報

この発明は上記の問題を解決するものであり、所望のドプラ速度レンジにてドプラスペクトラム画像を表示することが可能な超音波診断装置、及び超音波画像処理装置を提供することを目的とする。

請求項１に記載の発明は、被検体内の運動流体を含む診断部位に向けて、所定のパルス繰り返し周波数（ＰＲＦ）で超音波を送受信するスキャン手段と、前記スキャン手段によって取得された受信信号を記憶する記憶手段と、前記受信信号に対して周波数解析を行なってドプラスペクトラム画像を生成する演算手段と、前記ドプラスペクトラム画像を表示手段に表示させる表示制御手段と、前記表示手段に表示されたドプラスペクトラム画像に対する所望のドプラ速度レンジを指定するための入力手段と、前記入力手段によって入力された前記所望のドプラ速度レンジを受け付けて、前記記憶手段から前記受信信号を読み出し、前記読み出した前記受信信号に対して、前記所望のドプラ速度レンジに応じたサンプリング周波数に従ってリサンプリング処理を行う処理手段と、を有し、前記演算手段は、前記処理手段によってリサンプリング処理が行われた前記受信信号に対して、前記所望のドプラ速度レンジに応じた周波数解析を行なうことで新たなドプラスペクトラム画像を生成し、前記表示制御手段は、前記新たなドプラスペクトラム画像を前記表示手段に表示させることを特徴とする超音波診断装置である。
また、請求項６に記載の発明は、被検体内の運動流体を含む診断部位に向けて、所定のパルス繰り返し周波数（ＰＲＦ）で超音波を送受信することで取得された受信信号を記憶する記憶手段と、前記受信信号に対して周波数解析を行なってドプラスペクトラム画像を生成する演算手段と、前記ドプラスペクトラム画像を表示手段に表示させる表示制御手段と、前記表示手段に表示されたドプラスペクトラム画像に対する所望のドプラ速度レンジを指定するための入力手段と、前記入力手段によって入力された前記所望のドプラ速度レンジを受け付けて、前記記憶手段から前記受信信号を読み出し、前記読み出した前記受信信号に対して、前記所望のドプラ速度レンジに応じたサンプリング周波数に従ってリサンプリング処理を行う処理手段と、を有し、前記演算手段は、前記処理手段によってリサンプリング処理が行われた前記受信信号に対して、前記所望のドプラ測度レンジに応じた周波数解析を行なうことで新たなドプラスペクトラム画像を生成し、前記表示制御手段は、前記新たなドプラスペクトラム画像を前記表示手段に表示させることを特徴とする超音波画像処理装置である。

この発明によると、表示手段に表示されたドプラスペクトラム画像に対して所望のドプラ速度レンジを指定し、その指定されたドプラ速度レンジに応じたサンプリング周波数に従って、記憶手段に記憶されている受信信号に対してリサンプリング処理と周波数解析を行なうことで、時間分解能と周波数分解を劣化させずに、所望のドプラ速度レンジにてドプラスペクトラム画像を表示することが可能となる。

この発明の実施形態に係る超音波診断装置について図１及び図２を参照して説明する。図１は、この発明の実施形態に係る超音波診断装置を示すブロック図である。図２は、この発明の実施形態に係る超音波診断装置に設置された処理部を示すブロック図である。

この実施形態に係る超音波診断装置１は、超音波プローブ２、送受信部３、ドプラ処理部４、Ｂモード処理部１３、ＤＳＣ１４、表示制御部１５、及びユーザインターフェース（ＵＩ）１６を備えている。

超音波プローブ２には、複数の超音波振動子が所定方向（走査方向）に１例に配置された１次元アレイプローブ、又は、複数の超音波振動子が２次元的に配置された２次元アレイプローブが用いられる。

送受信部３は、超音波プローブ２に電気信号を供給して超音波を発生させる送信部（図示しない）と、超音波プローブ２からの信号を受信する受信部（図示しない）とを備えている。

送受信部３の送信部は、図示しないクロック発生回路、送信遅延回路、及びパルサ回路を備えている。クロック発生回路は、超音波信号の送信タイミングや送信周波数を決めるクロック信号を発生する回路である。送信遅延回路は、超音波の送信時に遅延を掛けて送信フォーカスを実施する回路である。パルサ回路は、各超音波振動子に対応した個別経路（チャンネル）の数分のパルサを内蔵し、遅延が掛けられた送信タイミングで駆動パルスを発生し、超音波プローブ２の各超音波振動子に供給するようになっている。送受信部３の送信部は、所定のパルス繰り返し周波数（ＰＲＦ）に従って、超音波プローブ２に超音波を送信させる。

また、送受信部３の受信部は、図示しないプリアンプ回路、Ａ／Ｄ変換回路、及び受信遅延・加算回路を備えている。プリアンプ回路は、超音波プローブ２の各超音波振動子から出力されるエコー信号を受信チャンネルごとに増幅する。Ａ／Ｄ変換回路は、増幅されたエコー信号をＡ／Ｄ変換する。受信遅延・加算回路は、Ａ／Ｄ変換後のエコー信号に対して受信指向性を決定するのに必要な遅延時間を与え、加算する。その加算により、受信指向性に応じた方向からの反射成分が強調される。なお、この送受信部３によって加算処理された信号を「ＲＦ信号」と称する。送受信部３から出力されたＲＦ信号は、ドプラ処理部４又はＢモード処理部１３に出力される。

Ｂモード処理部１３は、エコーの振幅情報の映像化を行い、エコー信号からＢモード超音波ラスタデータを生成する。具体的には、Ｂモード処理部１３は、ＲＦ信号に対してバンドパスフィルタ処理を行い、その後、出力信号の包絡線を検波し、検波されたデータに対して対数変換による圧縮処理を施す。Ｂモード処理部１３により生成されたＢモード超音波ラスタデータはＤＳＣ１４に出力される。

ＤＳＣ１４（Ｄｉｇｉｔａｌ Ｓｃａｎ Ｃｏｎｖｅｒｔｅｒ：デジタルスキャンコンバータ）は、直交座標系で表される画像を得るために、超音波ラスタデータを直交座標で表される画像データに変換する（スキャンコンバージョン処理）。ＤＳＣ１４から表示制御部１５に画像データが出力され、表示制御部１５は画像データに基づく画像を表示部１７に表示させる。例えば、ＤＳＣ１４は、Ｂモード超音波ラスタデータに基づいて２次元情報としての断層像データを生成し、その断層像データを表示制御部１５に出力する。表示制御部１５はその断層像データに基づく断層像を表示部１７に表示させる。

ドプラ処理部４は、直交位相検波部５、レンジゲート（ＲＧ）処理部６、ウォールフィルタ７、第１メモリ８、処理部９、第２メモリ１０、ＦＦＴ演算部１１、及びドプラ波形メモリ１２を備えている。

直交位相検波部５は、実部成分及び虚部成分にそれぞれ対応してデジタル方式のミキサ及びローパスフィルタを備え、ＲＦ信号を直交位相検波する。そして、直交位相検波部５は、直交位相検波されたデジタル量のＩＱ信号を出力する。この検波によりＲＦ信号から、ベースバンドのドプラ信号（実部成分及び虚部成分：ＩＱ信号）が抽出される。

レンジゲート（ＲＧ）処理部６は、直交位相検波部５から出力されたＩＱ信号を受けて、そのＩＱ信号の高周波成分を除去してドプラ偏移周波数成分のみからなるドプラ信号を得て、そのうち、被検体内の所望の深さのドプラ信号を抽出する。

ウォールフィルタ７は、レンジゲート処理部６からのレンジゲートで指定された被検体内の所定位置のドプラ信号（ＩＱ信号）から、比較的動きの遅い血管壁、心臓壁などの不要低周波ドプラ信号を除去して、検出すべき血流のドプラ信号（ＩＱ信号）を抽出する。そして、ウォールフィルタ７から出力されたＩＱ信号は、第１メモリ８に記憶される。

ＦＦＴ演算部１１は、ウォールフィルタ７で抽出されたドプラ信号（ＩＱ信号）に対して周波数解析を行なって、その解析結果であるドプラスペクトラム信号を得る。具体的には、ＦＦＴ演算部１１は、短時間フーリエ変換（ＳＴＦＴ）により、時間軸のドプラ信号（ＩＱ信号）に窓関数を掛けて、周波数軸の信号（ドプラスペクトラム信号）に変換する。ＦＦＴ演算部１１から出力されたドプラスペクトラム信号はドプラ波形メモリ１２に記憶される。

表示制御部１５は、ＦＦＴ演算部１１によって求められたドプラスペクトラム画像を表示部１７に表示させる。具体的には、表示制御部１５は、縦軸に周波数ｆ（速度ｖ）、横軸に時間ｔが表されたドプラスペクトラム画像を表示部１７に表示させる。

また、操作者が入力部１８を用いてフリーズ（静止表示）の指示を与えると、表示制御部１５は、ユーザインターフェース（ＵＩ）１６からそのフリーズ指示を受けて、そのフリーズ指示が与えられたタイミングで表示部１７に表示されているドプラスペクトラム画像を静止させて表示部１７に表示させる。すなわち、表示制御部１５は、ある時間の間に取得されたドプラスペクトラム画像を表示部１７に表示させる。その状態で、表示制御部１５は、一定の範囲を囲むことで、所望のドプラ速度レンジと所望の時間とを指定するためのＲＯＩマーカをドプラスペクトラム画像に重ねて表示部１７に表示させる。操作者は、このＲＯＩマーカによって、観察を所望するドプラ速度レンジと時間を指定することができる。

ここで、表示部１７に表示されるドプラスペクトラム画像について、図３を参照して説明する。図３は、ドプラ速度レンジを縮小する場合におけるドプラスペクトラム画像を示す図である。

図３（ａ）に示すように、表示制御部１５は、縦軸が速度ｖ（周波数）で横軸が時間ｔで表されるドプラスペクトラム画像１００を表示部１７に表示させる。さらに、表示制御部１５は、一定の範囲を囲むＲＯＩマーカ２００をドプラスペクトラム画像１００に重ねて表示部１７に表示させる。この実施形態では、表示制御部１５は、矩形状の範囲を囲むＲＯＩマーカ２００を生成して、ドプラスペクトラム画像１００に重ねて表示部１７に表示させる。

操作者は表示部１７に表示されたドプラスペクトラム画像１００を観察しながら、入力部１８を用いて、ＲＯＩマーカ２００をドプラスペクトラム画像１００上で移動させて、所望のドプラ速度レンジと時間とを指定する。表示制御部１５は、入力部１８にて入力された移動指示に従って、ＲＯＩマーカ２００をドプラスペクトラム画像１００上で移動させて表示部１７に表示させる。例えば、操作者は入力部１８を用いて、ＲＯＩマーカ２００の縦軸方向（速度レンジ方向）の幅を変えたり、横軸方向（時間軸方向）の幅を変えたり、ドプラスペクトラム画像１００上での位置を変えたりして、所望のドプラ速度レンジと所望の時間とを指定する。また、操作者は入力部１８を用いて所望の速度オフセット（ＢＬＳ）を指定する。

操作者によって指定された所望のドプラ速度レンジと所望の時間は、ユーザインターフェース１６から表示制御部１５に出力される。表示制御部１５は、指定された時間に取得されたドプラスペクトラム画像をドプラ波形メモリ１２から読み込み、指定されたドプラ速度レンジに従って拡大又は縮小したドプラスペクトラム画像を表示部１７に表示させる。例えば、ドプラ速度レンジを小さくすることで、ドプラスペクトラム画像を拡大表示することができ、ドプラ速度レンジを大きくすることで、ドプラスペクトラム画像を縮小表示することができる。

例えば図３（ａ）に示すように、ドプラ速度レンジが大きすぎるとスケールが大きくなるため、ドプラスペクトラム画像が小さくなってしまう。この場合、図３（ａ）に示すように、ＲＯＩマーカ２００によってドプラスペクトラム画像１００の一部分を指定すると、図３（ｂ）に示すように、表示制御部１５は、ＲＯＩマーカ２００によって指定された範囲に含まれるドプラスペクトラム画像１１０をドプラ波形メモリ１２から読み出して表示部１７に表示させる。これにより、ドプラスペクトラム画像１００の一部分が拡大表示されることになる。つまり、ＲＯＩマーカ２００によってより小さいドプラ速度レンジを指定することで、そのドプラ速度レンジに納まるドプラスペクトラム画像が表示されることになる。なお、この段階においては、ドプラスペクトラム画像１１０は、ドプラスペクトラム画像１１０を単純に拡大させた画像（デジタルズームの画像）であるため、画面上における画素あたりの画質が劣化している。また、操作者が入力部１８を用いて速度オフセット（ＢＬＳ）を調整することで、表示制御部１５は、その速度オフセット（ＢＬＳ）に従って、ドプラスペクトラム画像をドプラ速度レンジ（縦軸）の所望の位置に表示させる。

ドプラスペクトラム画像の別の例について図４を参照して説明する。図４は、ドプラ速度レンジを拡大する場合におけるドプラスペクトラム画像を示す図である。例えばドプラ速度レンジが小さすぎると、図４（ａ）に示すように、ドプラスペクトラム画像３００にエイリアシング現象（折り返り）が発生して表示部１７に表示される。この場合、図４（ｂ）に示すように、表示制御部１５は、ドプラ速度レンジを大きくして、ドプラスペクトラム画像３１０を表示部１７に表示させる。これにより、ドプラスペクトラム画像３００が縮小表示されることになる。さらに、表示制御部１５は、ドプラスペクトラム画像３１０にＲＯＩマーカ４００を重ねて表示部１７に表示させる。そして、ＲＯＩマーカ４００によってドプラスペクトラム画像３１０の一部分を指定すると、表示制御部１５は、ＲＯＩマーカ４００によって指定された範囲に含まれるドプラスペクトラム画像をドプラ波形メモリ１２から読み出して表示部１７に表示させる。これにより、ドプラスペクトラム画像３１０の一部分が拡大表示されることになる。なお、この段階においては、拡大されたドプラスペクトラム画像は、単純に拡大された画像（デジタルズームの画像）であるため、画面上における画素あたりの画質が劣化している。また、操作者が入力部１８を用いて速度オフセット（ＢＬＳ）を調整することで、表示制御部１５は、その速度オフセット（ＢＬＳ）に従って、ドプラスペクトラム画像をドプラ速度レンジ（縦軸）の所望の位置に表示させる。

そして、所望の時相間に取得されたドプラスペクトラム画像が所望のドプラ速度レンジ内に納まって表示部１７に表示されるように、操作者は入力部１８を用いて、時間、ドプラ速度レンジ、及び速度オフセット（ＢＬＳ）を調整する。入力部１８を用いて指定された時間、ドプラ速度レンジ、及び速度オフセット（ＢＬＳ）は、ユーザインターフェース（ＵＩ）１６からドプラ処理部４に出力される。

処理部９は、操作者によって指定されたドプラ速度レンジに応じたサンプリング周波数に従って、第１メモリ８から指定された時間に取得されたＩＱ信号を読み出してリサンプリング処理を行う。このとき、指定されたドプラ速度レンジに応じて、パルス繰り返し周波数（ＰＲＦ）に相当するサンプリング周波数を有理数倍に変更したサンプリング周波数でリサンプリング処理を行う。具体的には、インターポレータとデシメータを用いることで、パルス繰り返し周波数（ＰＲＦ）に相当するサンプリング周波数を有理数倍に変更したサンプリング周波数でリサンプリング処理を行う。

ここで、処理部９の構成を図２に示す。図２に示すように、処理部９は、アップサンプリング部９１、複素バンドパスフィルタ（Ｃｏｍｐｌｅｘ ＢＰＦ：ＣＢＰＦ）９２、９３、設定部９４、合成部９５、及びダウンサンプリング部９６を備えている。アップサンプリング部９１は、複素数信号であるドプラ信号（ＩＱ信号）を入力する実部成分用及び虚部成分用の２チャンネル分の回路で構成されている。また、一方のＣＢＰＦ９２は、ドプラ信号（ＩＱ信号）の正側周波数成分（Ｆｏｒｗａｒｄ成分）側に割り当てられている。他方のＣＢＰＦ９３は、ドプラ信号（ＩＱ信号）の負側周波数成分（Ｒｅｖｅｒｓｅ成分）側に割り当てられている。

例えば、操作者によって指定されたドプラ速度レンジに対応するサンプリング周波数ｆｓ’でリサンプリングする場合、以下の式（１）に従って、サンプリング周波数ｆｓ’を規定する。
ｆｓ’＝（Ｎ２／Ｎ１）×ｆｓ・・・式（１）
ここで、ｆｓはサンプリング周波数（パルス繰り返し周波数ＰＲＦに相当する）である。
Ｎ１、Ｎ２は整数である。
すなわち、サンプリング周波数ｆｓを有理数倍に変更することで、サンプリング周波数ｆｓ’を規定する。

アップサンプリング部９１は、入力されたドプラ信号（ＩＱ信号）の信号値の相互間にゼロ（０）値を挿入することで、出力値のサンプリング周波数を増加させる。これは、ゼロ挿入によって、ドプラ信号の時系列方向の長さが整数Ｎ２倍（Ｎ２＝２、３、・・・）になるように、実部チャンネル及び虚部チャンネルごとになされる。具体的には、ドプラ信号の信号値の相互間に（Ｎ２−１）個のゼロ（０）値を挿入することで、出力値のサンプリング周波数をＮ２倍に増加させる。しかしながら、ナイキスト周波数がＮ２倍になるため、アップサンプリングにより新たに生じたスペクトル成分を除去するために、アップサンプリング後の信号に対してデジタルローパスフィルタ処理を施す。このデジタルローパスフィルタ処理は、ＣＢＰＦ９２、９３によって行われる。

ゼロ挿入されたドプラ信号はＣＢＰＦ９２、９３にそれぞれ送られ、バンドパスフィルタリングの処理が施される。このバンドパスフィルタリングを行う場合、ＣＢＰＦ９２、９３それぞれの係数がベースラインのシフト量（ＢＬＳ）に応じて制御される。このシフト量（ＢＬＳ）の情報は、ユーザインターフェース（ＵＩ）１６から処理部９に出力される。

この係数制御されたＣＢＰＦ９２、９３は、ドプラ信号（ＩＱ信号）をフィルタリングする。まず、ドプラスペクトラム画像を表示部１７に表示させた状態で、操作者が入力部１８を操作することで、ドプラスペクトラム画像のベースラインＢＬをシフトさせる。このシフト量（ＢＬＳ）の情報は、ユーザインターフェース（ＵＩ）１６から処理部９に出力される。そして、設定部９４は、ユーザインターフェース（ＵＩ）１６から出力されたベースラインのシフト量（ＢＬＳ）に応じて、例えば特開２００２−３２５７６７号公報に記載の方法によって、ＣＢＰＦ９２、９３のバンド幅及び中心周波数を設定する。この設定は、アップサンプリング部９１によるゼロ挿入の倍数Ｎ２と、ダウンサンプリング部９６による倍数（１／Ｎ１）に対応して行われる。

ここで、ドプラ信号（ＩＱ信号）の正側周波数成分（Ｆｏｒｗａｒｄ成分）側に割り当てられたＣＢＰＦ９２におけるフィルタリング係数Ｚ’とカットオフ周波数ｆｃ（ｆｗｄ）は、以下の式（２）、（３）で与えられる。
フィルタリング係数Ｚ’＝ｅｘｐ（−ｊπ＊（ＢＬＳ／２））・・・式（２）
カットオフ周波数ｆｃ（ｆｗｄ）＝（Ｎ２／Ｎ１）＊ｆｓ＊（０．５＋ＢＬＳ）・・式（３）

また、ドプラ信号（ＩＱ信号）の負側周波数成分（Ｒｅｖｅｒｓｅ成分）側に割り当てられたＣＢＰＦ９３におけるフィルタリング係数Ｚ’’とカットオフ周波数ｆｃ（ｒｅｖ）は、以下の式（４）、（５）で与えられる。
フィルタリング係数Ｚ’’＝ｅｘｐ（−ｊπ＊（０．５−ＢＬＳ）／２）・・・式（４）
カットオフ周波数ｆｃ（ｒｅｖ）＝（Ｎ２／Ｎ１）＊ｆｓ＊（ＢＬＳ−０．５）・・式（５）

ＣＢＰＦ９２、９３は、上記のフィルタリング係数とカットオフ周波数に従って、ドプラ信号（ＩＱ信号）をフィルタリングする。

合成部９５は、ＣＢＰＦ９２から出力されたドプラ信号の正側周波数成分と、ＣＢＰＦ９３から出力されたドプラ信号の負側周波数成分とを合成してダウンサンプリング部９６に出力する。

ダウンサンプリング部９６は、入力されたドプラ信号（ＩＱ信号）のサンプリング周波数を（１／Ｎ１）倍にする。ダウンサンプリング部９６は、入力されたドプラ信号をＮ１個間隔でピックアップすることにより、出力信号のサンプリング周波数を減少させる。しかしながら、ダウンサンプル後はナイキスト周波数が半分になるため、折り返しのエイリアシング成分がドプラ信号成分に重なる可能性がある。そのため、ダウンサンプルする前にデジタルローパスフィルタ処理を施して、ナイキスト周波数以上の高周波成分を除去しておく。このデジタルローパスフィルタ処理は、ＣＢＰＦ９２、９３によって行われる。そして、処理部９から出力されたサンプリング処理後のＩＱ信号は第２メモリ１０に記憶される。

例えば、フリーズ状態のドプラスペクトラム画像において、ドプラスペクトラム画像を拡大表示する場合は、ドプラ速度レンジを縮小する。１例として、ドプラ速度レンジを２／３に縮小する場合、アップサンプリング部９１は、Ｎ２＝２としてアップサンプリングを行い、ダウンサンプリング部９６は、Ｎ１＝３としてダウンサンプリングを行う。一方、ドプラスペクトラム画像を縮小表示する場合は、ドプラ速度レンジを拡大する。１例として、ドプラ速度レンジを３／２に拡大する場合、アップサンプリング部９１は、Ｎ２＝３としてアップサンプリングを行い、ダウンサンプリング部９６は、Ｎ１＝２としてダウンサンプリングを行う。このように、処理部９は、パルス繰り返し周波数ＰＲＦを有理数倍に変更したサンプリング周波数でリサンプリング処理を行う。

なお、ドプラ速度レンジを拡大する場合であって、２＊ｆｓ≦ｆｓ’＜４＊ｆｓとなる場合、処理部９は準同形処理を行う。例えば、処理部９による処理について、ａｎｔｉａｌｉｅｓ（ｆ）という関数を定義し、ｆｓ’’＝ａｎｔｉａｉｅｓ（ｆｓ’）として処理を行う。

ＦＦＴ演算部１１は、第２メモリ１０からリサンプリング後のＩＱ信号を読み出して、短時間フーリエ変換（ＳＴＦＴ）により、時間軸のドプラ信号（ＩＱ信号）に窓関数を掛けて、周波数軸の信号（ドプラスペクトラム信号）に変換する。このとき、ＦＦＴ演算部１１は、操作者によって指定されたドプラ速度レンジ（サンプリング周波数）の変更に応じて、例えば特開平１１−２８５４９５号公報に記載の方法によって、周波数分解能と時間分解能とを調整してドプラスペクトラム信号を生成する。ここで、この処理について簡単に説明する。

時間分解能ΔＸと周波数分解能ΔＹは、以下の式（６）、（７）で求められる。
ΔＸ＝（Ｄｘ×Ｎ×Ｋ×Ｍ）／（ＰＲＦ×Ｓ）・・・式（６）
ΔＹ＝Ｄｙ／（Ｎ×Ｗｗ）・・・式（７）
Ｄｘ：ドプラスペクトラム画像の表示範囲の横軸（時間軸）方向の画素数
Ｄｙ：ドプラスペクトラム画像の表示範囲の縦軸（周波数軸）方向の画素数
Ｎ：ＦＦＴのデータ数
Ｋ：交互段数
Ｍ：間引き率
ＰＲＦ：パルス繰り返し周波数（ドプラ速度レンジに相当する）
Ｓ：スクロール速度
Ｗｗ：窓重み

なお、ＦＦＴのデータ数Ｎは、２のべき乗という制約があるが、目的の周波数分解能を実現するための計算上のデータ数は２のべき乗にならない場合が多い。そこで、ＦＦＴのデータ数Ｎには、計算上のデータ数以上でそれに最も近い２のべき乗を選び、その選んだＦＦＴのデータ数Ｎと、計算上のデータ数との不一致を、ＦＦＴの窓関数を掛けることで埋めるようにしている。この窓の大きさ（窓幅）を決めるパラメータが、窓重みＷｗである。例えば計算上のデータ数が３２０である場合に、ＦＦＴのデータ数Ｎとしては５１２が選択され、この場合、窓重みＷｗは５／８に設定される。例えば、時間分解能に対する周波数分解能の比を一定に維持する場合、ＦＦＴ演算部１１は、パルス繰り返し周波数ＰＲＦ（ドプラ速度レンジ）、スクロール速度Ｓ、交互段数Ｋ、間引き率Ｍ、Ｄｙが与えられると、操作者の所望の比に維持するためのＦＦＴのデータ数Ｎと窓重みＷｗ（窓幅）を計算するようになっている。つまり、パルス繰り返し周波数ＰＲＦ（ドプラ速度レンジ）が変わった場合（時間分解能が変わった場合）、それに応じてＦＦＴのデータ数Ａ４を調整して周波数分解能を変えて、時間分解能に対する周波数分解能の比を一定にすることになる。

このデータ数Ａ４は、上述したようにＦＦＴで扱える２のべき乗には必ずしもならない。そこで、計算上のデータ数Ａ４以上でそれに最も近い２のべき乗を選択し、その選択したＦＦＴのデータ数Ｎと、両者間の不一致を埋めるための窓重みＷｗとを選択する。このような処理によって、ドプラ速度レンジが変わっても、時間分解能ΔＸに対する周波数分解能ΔＹの比を一定に保つことができる。

また、周波数分解能が変動しても、時間分解能を操作者が所望する一定値に固定するような処理を行っても良い。例えば、データ数Ｎが変わったとき、ＦＦＴ演算部１１は、上記式（６）において、データ数Ｎ以外のパラメータのうち少なくとも１つを変えて、時間分解能を一定値に維持する。また、時間分解能が変動しても、周波数分解能を操作者が所望する一定値に固定するような処理を行っても良い。

ドプラ速度レンジと速度オフセットとが変更されたドプラスペクトラム信号は、ＦＦＴ演算部１１からドプラ波形メモリ１２に出力されて、ドプラ波形メモリ１２に記憶される。そして、表示制御部１５は、ドプラ速度レンジなどが変更されたドプラスペクトラム画像をドプラ波形メモリ１２から読み出して、表示部１７に表示させる。

ここで、新たに生成されたドプラスペクトラム画像の１例を図３（ｃ）と図４（ｃ）にそれぞれ示す。例えば図３（ｃ）に示すドプラスペクトラム画像１２０は、ドプラ速度レンジと速度オフセットとが変更されて、処理部９にてリサンプリング処理が施され、さらに、ＦＦＴ演算部１１にてＳＴＦＴ処理が施されて得られた画像である。一方、図３（ｂ）に示すドプラスペクトラム画像１１０は、ドプラスペクトラム画像１１０を単純に拡大して表示した画像である。ドプラスペクトラム画像１１０は、単純なデジタルズームの画像であるため、画面上における画素あたりの画質が劣化している。一方、再度サンプリング処理とＳＴＦＴ処理が施されて得られたドプラスペクトラム画像１２０は、時間分解能と周波数分解能が劣化していない。このように、この実施形態に係る超音波診断装置１によると、時間分解能と周波数分解能を劣化させずに、所望のドプラ速度レンジと速度オフセットで表されるドプラスペクトラム画像１２０を生成することが可能となる。

また、図４（ｃ）に示すドプラスペクトラム画像３２０も同様に、ドプラ測度レンジと速度オフセットとが変更されて、処理部９にてリサンプリング処理が施され、さらに、ＦＦＴ演算部１１にてＳＴＦＴ処理が施されて得られた画像である。このドプラスペクトラム画像３２０も、時間分解能と周波数分解能が劣化せずに生成されている。

また、操作者によって所望のドプラ速度レンジ（ＰＲＦ）と速度オフセット（ＢＬＳ）が決定されると、ドプラ速度レンジ（ＰＲＦ）がユーザインターフェース（ＵＩ）１６から送受信部３に出力され、速度オフセット（ＢＬＳ）がユーザインターフェース（ＵＩ）１６からＦＦＴ演算部１１に出力される。送受信部３は、ユーザインターフェース（ＵＩ）１６から出力されたドプラ速度レンジ（ＰＲＦ）を受けると、そのドプラ速度レンジ（ＰＲＦ）に従って超音波プローブ２に超音波を送受信させる。また、ＦＦＴ演算部１１は、ユーザインターフェース（ＵＩ）１６から速度オフセット（ＢＬＳ）を受けると、そのシフト量だけＦＦＴ処理の読み出しアドレスを変化させて速度のオフセットを調整する。このように、ドプラ速度レンジと速度オフセットが変更されると、変更されたドプラ速度レンジに従ってスキャンを行ない、変更された速度オフセットによってＦＦＴ処理を行う。

なお、ユーザインターフェース（ＵＩ）１６は、表示部１７と入力部１８とを備えている。表示部１７は、ＣＲＴや液晶ディスプレイなどのモニタで構成されており、画面上に断層像、３次元画像又はドプラスペクトラム画像などが表示される。入力部１８は、キーボード、マウス、トラックボール、又はＴＣＳ（Ｔｏｕｃｈ Ｃｏｍｍａｎｄ Ｓｃｒｅｅｎ）などで構成されている。操作者は、入力部１８を用いることで、表示部１７の画面上でＲＯＩマーカを所望の位置に移動させることができる。

また、入力部１８に、ＲＯＩマーカの大きさと位置を変えるスイッチを設けても良い。例えば、ＲＯＩマーカの縦方向（ドプラ速度レンジ）の拡大又は縮小を指示するためのスイッチＡや、ＲＯＩマーカの横方向（時間軸）の拡大又は縮小を指示するためのスイッチＢを入力部１８に設けても良い。また、入力部１８に設置されたトラックボールによって、ＲＯＩマーカの位置を変えて、観察対象となる時相を選択するようにしても良い。そして、ドプラスペクトラム画像を静止させた状態で、操作者がスイッチＡ、スイッチＢ、及びトラックボールを用いることで、ＲＯＩマーカを所望の位置に移動させ、また、ＲＯＩマーカの大きさを所望の大きさに変える。これにより、所望のドプラ速度レンジと所望の時間を指定することができる。

また、超音波診断装置１には、図示しない制御部が設置されている。制御部は、送受信部３、ドプラ処理部４、Ｂモード処理部１３、ＤＳＣ１４、表示制御部１５、及びユーザインターフェース（ＵＩ）１６に接続され、各部の動作を制御する。

（超音波画像処理装置）
また、超音波の送受信によって取得された受信信号に基づいてドプラスペクトラム画像を生成する超音波画像処理装置を、超音波診断装置の外部に設けても良い。この超音波画像処理装置は、上述したドプラ処理部４、表示制御部１５、及びユーザインターフェース（ＵＩ）１６を備えている。そして、超音波画像処理装置は、操作者によって指定された所望のドプラ速度レンジに応じたサンプリング周波数に従って、受信信号に対してリサンプリング処理を行い、更に、ＳＴＦＴ処理を行うことで、そのドプラ速度レンジで表されたドプラスペクトラム画像を生成する。

（動作）
次に、この発明の実施形態に係る超音波診断装置による一連の動作について図５を参照して説明する。図５は、この発明の実施形態に係る超音波診断装置による一連の動作を示すフローチャートである。

（ステップＳ０１）
まず、超音波プローブ２と送受信部３とによって、初期設定された所定のパルス繰り返し周波数（ＰＲＦ）に従って、被検体に対して超音波を送信し、被検体からの反射波を受信する。送受信部３によって受信された受信信号は、直交位相検波部５、レンジゲート処理部６、及びウォールフィルタ７によって処理が施される。そして、ドプラ信号（ＩＱ信号）が第１メモリ８に順次、記憶されていく。ＦＦＴ演算部１１は、ドプラ信号に対して周波数解析を行なうことでドプラスペクトラム画像を順次、生成していく。ＦＦＴ演算部１１によって生成されたドプラスペクトラム画像は、ドプラ波形メモリ１２に順次、記憶される。

（ステップＳ０２）
表示制御部１５は、ドプラ波形メモリ１２からドプラスペクトラム画像を読み出して、縦軸に周波数ｆ（速度ｖ）、横軸に時間ｔが表されたドプラスペクトラム画像を表示部１７に表示させる。

（ステップＳ０３）
そして、操作者は表示部１７に表示されたドプラスペクトラム画像を観察しながら、入力部１８を用いて、所望のタイミングでフリーズ（静止表示）の指示を与える。

（ステップＳ０４）
表示制御部１５は、ユーザインターフェース（ＵＩ）１６からそのフリーズ指示を受けて、そのフリーズ指示が与えられたタイミングで表示部１７に表示されているドプラスペクトラム画像を静止させて表示部１７に表示させる。

（ステップＳ０５）
さらに、表示制御部１５は、一定の範囲を囲むＲＯＩマーカをドプラスペクトラム画像に重ねて表示部１７に表示させる。例えば図３（ａ）に示すように、表示制御部１５は、ドプラスペクトラム画像１００を固定させて表示部１７に表示させ、そのドプラスペクトラム画像１００上にＲＯＩマーカ２００を重ねて表示させる。そして、操作者は表示部１７に表示されているドプラスペクトラム画像１００を観察しながら、入力部１８を用いて、ＲＯＩマーカ２００の縦軸方向（速度レンジ方向）の幅を変えたり、横軸方向（時間軸方向）の幅を変えたり、位置を変えたりして、所望のドプラ速度レンジと所望の時間とを指定する。また、操作者は入力部１８を用いて所望の速度オフセット（ＢＬＳ）を指定する。このように、操作者はＲＯＩマーカを用いて、観察を所望するドプラ速度レンジと時間とを指定する。

（ステップＳ０６）
表示制御部１５は、ＲＯＩマーカ２００によって指定された時間に取得されたドプラスペクトラム画像をドプラ波形メモリ１２から読み込み、指定されたドプラ速度レンジに従って拡大又は縮小したドプラスペクトラム画像を表示部１７に表示させる。例えば図３（ａ）に示すように、ＲＯＩマーカ２００によってドプラスペクトラム画像１００の一部分を指定すると、図３（ｂ）に示すように、表示制御部１５は、ＲＯＩマーカ２００によって指定された範囲に含まれるドプラスペクトラム画像１１０をドプラ波形メモリ１２から読み出して表示部１７に表示させる。これにより、ドプラスペクトラム画像１００の一部分が拡大表示されることになる。なお、この段階においては、ドプラスペクトラム画像１１０は、ドプラスペクトラム画像１００を単純に拡大させた画像（デジタルズームの画像）であるため、画面上における画素あたりの画質が劣化している。

（ステップＳ０７）
そして、操作者は拡大表示又は縮小表示されたドプラスペクトラム画像を観察し、そのドプラスペクトラム画像がドプラ速度レンジに収まって表示されているか確認する。また、表示されているドプラスペクトラム画像の速度オフセットが適した値になっているか確認する。すなわち、観察を所望するドプラスペクトラム画像に対して、ドプラ測度レンジが小さすぎたり大きすぎたりしないか確認する。拡大表示又は縮小表示されたドプラスペクトラム画像がドプラ速度レンジに納まって表示されていない場合や、速度オフセットが適した値になっていない場合は（ステップＳ０７、Ｎｏ）、操作者はＲＯＩマーカによって範囲を指定し、その範囲内のドプラスペクトラム画像を観察することで、ドプラスペクトラム画像を適したドプラ速度レンジ内に表示させる（ステップＳ０５、Ｓ０６）。すなわち、観察を所望するドプラスペクトラム画像に対してドプラ速度レンジが小さすぎたり大きすぎたりする場合は、再度、ＲＯＩマーカによって所望の範囲を指定する。

一方、拡大表示又は縮小表示されたドプラスペクトラム画像がドプラ速度レンジに納まって表示され、また、速度オフセットが適した値になっている場合（ステップＳ０７、Ｙｅｓ）、操作者は入力部１８を用いて再処理の指示を与える。これにより、操作者によって与えられた所望のドプラ速度レンジと速度オフセットとが、ユーザインターフェース（ＵＩ）１６から処理部９に出力される。

（ステップＳ０８）
処理部９は、第１メモリ８からドプラ信号（ＩＱ信号）を読み出し、指定されたドプラ速度レンジに応じて、パルス繰り返し周波数（ＰＲＦ）に相当するサンプリング周波数を整数倍に変えたサンプリング周波数で、そのドプラ信号に対してアップサンプリングと、複素バンドパスフィルタによるフィルタリングと、ダウンサンプリングとを行うことで、リサンプリング処理を行う。

（ステップＳ０９）
ＦＦＴ演算部１１は、新たなサンプリング周波数に従ってリサンプリング処理が施されたドプラ信号に対してＳＴＦＴ処理を施すことで、新たなドプラスペクトラム信号を生成する。このドプラスペクトラム信号はドプラ波形メモリ１２に記憶される。

（ステップＳ１０）
そして、表示制御部１５は、ドプラ波形メモリ１２から新たなドプラスペクトラム画像を読み込んで、表示部１７に表示させる。これにより、例えば図３（ｃ）に示すように、新たなドプラスペクトラム画像１２０は、時間分解能と周波数分解能が劣化せずに生成されて表示されることになる。

以上のように、この実施形態に係る超音波診断装置１によると、フリーズ状態のドプラスペクトラム画像を参照しながらドプラ速度レンジと速度オフセットを変更することで、ドプラスペクトラム画像に適したドプラ速度レンジと速度オフセットを簡便に設定することができる。すなわち、静止しているドプラスペクトラム画像を拡大表示又は縮小表示することで、所望のドプラ速度レンジと速度オフセットを確認することができるため、所望の値を簡便に設定することが可能となる。

そして、設定されたドプラ速度レンジと速度オフセットに基づいて、ＦＦＴ処理前のドプラ信号（ＩＱ信号）に対してリサンプリング処理を行い、更に、ＳＴＦＴ処理を施すことで、時間分解能と周波数分解能とを劣化させずに、所望のドプラ速度レンジで表されたドプラスペクトラム画像を生成して表示することが可能となる。すなわち、単純にドプラスペクトラム画像をドプラ速度レンジに応じて拡大又は縮小させるだけでなく、設定されたドプラ速度レンジに基づいて再度、リサンプリング処理とＳＴＦＴ処理を行うことで、時間分解能と周波数分解能を劣化させずに、所望のドプラ速度レンジで表されたドプラスペクトラム画像が得られる。このように、画質を劣化させずに、所望のドプラ測度レンジで表されたドプラスペクトラム画像が得られるため、そのドプラスペクトラム画像に基づいて所望の計測を行うことが可能となる。

また、計測したい時間に取得されたドプラスペクトラム画像を、所望のドプラ速度レンジと速度オフセットに基づいて再度処理して新たなドプラスペクトラム画像を生成しているため、計測したい時間の新たなドプラスペクトラム画像に基づいて所望の計測を行うことができる。一方、従来技術においては、過去に取得されたドプラスペクトラム画像に基づいて現在以降のドプラ速度レンジを推定するため、推定した後に変動が発生した場合、再度、推定処理を行う必要があるため、安定して計測を行うことが困難であった。これに対して、この実施形態に係る超音波診断装置１によると、過去のドプラスペクトラム画像に基づく推定処理を行う必要がないため、安定してドプラスペクトラム画像による計測を行うことが可能となる。

また、所望のドプラ速度レンジと速度オフセットが決定されると、送受信部３は、そのドプラ速度レンジ（ＰＲＦ）に従って超音波プローブ２に超音波を送受信させる。また、ＦＦＴ演算部１１は、その速度オフセット（ＢＬＳ）に従い、シフト量だけＦＦＴ処理の読み出しアドレスを変化させて速度のオフセットを調整する。これにより、ドプラ速度レンジ（ＰＲＦ）と速度オフセット（ＢＬＳ）が更新される。

この発明の実施形態に係る超音波診断装置を示すブロック図である。 この発明の実施形態に係る超音波診断装置に設置された処理部を示すブロック図である。 ドプラ速度レンジを縮小する場合におけるドプラスペクトラム画像を示す図である。 ドプラ速度レンジを拡大する場合におけるドプラスペクトラム画像を示す図である。 この発明の実施形態に係る超音波診断装置による一連の動作を説明するためのフローチャートである。

符号の説明

１ 超音波診断装置
２ 超音波プローブ
３ 送受信部
４ ドプラ処理部
５ 直交位相検波部
６ レンジゲート処理部
７ ウォールフィルタ
８ 第１メモリ
９ 処理部
１０ 第２メモリ
１１ ＦＦＴ演算部
１２ ドプラ波形メモリ
１５ 表示制御部
１６ ユーザインターフェース（ＵＩ）
１７ 表示部
１８ 入力部
９１ アップサンプリング部
９２、９３ ＣＢＰＦ（Ｃｏｍｐｌｅｘ ＢＰＦ）
９４ 設定部
９５ 合成部
９６ ダウンサンプリング部

Claims (6)

1. 被検体内の運動流体を含む診断部位に向けて、所定のパルス繰り返し周波数（ＰＲＦ）で超音波を送受信するスキャン手段と、
   前記スキャン手段によって取得された受信信号を記憶する記憶手段と、
   前記受信信号に対して周波数解析を行なってドプラスペクトラム画像を生成する演算手段と、
   前記ドプラスペクトラム画像を表示手段に表示させる表示制御手段と、
   前記表示手段に表示されたドプラスペクトラム画像に対する所望のドプラ速度レンジを指定するための入力手段と、
   前記入力手段によって入力された前記所望のドプラ速度レンジを受け付けて、前記記憶手段から前記受信信号を読み出し、前記読み出した前記受信信号に対して、前記所望のドプラ速度レンジに応じたサンプリング周波数に従ってリサンプリング処理を行う処理手段と、
   を有し、
   前記演算手段は、前記処理手段によってリサンプリング処理が行われた前記受信信号に対して、前記所望のドプラ速度レンジに応じた周波数解析を行なうことで新たなドプラスペクトラム画像を生成し、
   前記表示制御手段は、前記新たなドプラスペクトラム画像を前記表示手段に表示させることを特徴とする超音波診断装置。
2. 前記表示制御手段は、前記入力手段を介して前記ドプラスペクトラム画像の静止指示を受け付けて、前記ドプラスペクトラム画像を静止させた状態で前記表示手段に表示させ、前記静止させられたドプラスペクトラム画像上に、所望の範囲を囲むことで前記所望のドプラ速度レンジと所望の時間幅とを指定するためのＲＯＩマーカを重ねて前記表示手段に表示させ、
   前記処理手段は、前記記憶手段から前記所望の時間幅に取得された受信信号を読み出し、前記読み出した前記受信信号に対して、前記ＲＯＩマーカによって指定された前記所望のドプラ速度レンジに応じたサンプリング周波数に従ってリサンプリング処理を行うことを特徴とする請求項１に記載の超音波診断装置。
3. 前記処理手段は、前記所望のドプラ速度レンジに応じて、前記所定のパルス繰り返し周波数を有理数倍に変えたサンプリング周波数に従ってリサンプリング処理を行うことを特徴とする請求項１又は請求項２のいずれかに記載の超音波診断装置。
4. 前記処理手段は、前記所望のドプラ速度レンジに応じて、前記所定のパルス繰り返し周波数を整数倍に変えたサンプリング周波数でアップサンプリングし、その後、前記所定のパルス繰り返し周波数を整数倍に変えたサンプリング周波数でダウンサンプリングすることで、前記リサンプリング処理を行うことを特徴とする請求項３に記載の超音波診断装置。
5. 前記演算手段は、前記処理手段によってリサンプリング処理が行われた前記受信信号に対して、前記所望のドプラ速度レンジに応じて、データ数と窓関数の幅とを変えて前記周波数解析を行なうことで新たなドプラスペクトラム画像を生成することを特徴とする請求項１から請求項４のいずれかに記載の超音波診断装置。
6. 被検体内の運動流体を含む診断部位に向けて、所定のパルス繰り返し周波数（ＰＲＦ）で超音波を送受信することで取得された受信信号を記憶する記憶手段と、
   前記受信信号に対して周波数解析を行なってドプラスペクトラム画像を生成する演算手段と、
   前記ドプラスペクトラム画像を表示手段に表示させる表示制御手段と、
   前記表示手段に表示されたドプラスペクトラム画像に対する所望のドプラ速度レンジを指定するための入力手段と、
   前記入力手段によって入力された前記所望のドプラ速度レンジを受け付けて、前記記憶手段から前記受信信号を読み出し、前記読み出した前記受信信号に対して、前記所望のドプラ速度レンジに応じたサンプリング周波数に従ってリサンプリング処理を行う処理手段と、
   を有し、
   前記演算手段は、前記処理手段によってリサンプリング処理が行われた前記受信信号に対して、前記所望のドプラ測度レンジに応じた周波数解析を行なうことで新たなドプラスペクトラム画像を生成し、
   前記表示制御手段は、前記新たなドプラスペクトラム画像を前記表示手段に表示させることを特徴とする超音波画像処理装置。

Images