以下、添付図面を参照して、超音波診断装置の実施形態を詳細に説明する。
(実施形態)
まず、本実施形態に係る超音波診断装置の構成について説明する。図1は、本実施形態に係る超音波診断装置の構成例を示すブロック図である。図1に例示するように、本実施形態に係る超音波診断装置は、超音波プローブ1と、モニタ2と、入力装置3と、装置本体10とを有する。また、超音波プローブ1には、位置センサ4が取り付けられ、装置本体10の近傍には、トランスミッター5が設置される。また、装置本体10は、ネットワークを介して外部装置6と接続される。
超音波プローブ1は、例えば、複数の圧電振動子を有し、これら複数の圧電振動子は、後述する装置本体10が有する送受信部11から供給される駆動信号に基づき超音波を発生する。また、超音波プローブ1は、被検体Pからの反射波を受信して電気信号に変換する。また、超音波プローブ1は、圧電振動子に設けられる整合層と、圧電振動子から後方への超音波の伝播を防止するバッキング材等を有する。なお、超音波プローブ1は、装置本体10と着脱自在に接続される。
超音波プローブ1から被検体Pに超音波が送信されると、送信された超音波は、被検体Pの体内組織における音響インピーダンスの不連続面で次々と反射され、反射波信号として超音波プローブ1が有する複数の圧電振動子にて受信される。受信される反射波信号の振幅は、超音波が反射される不連続面における音響インピーダンスの差に依存する。なお、送信された超音波パルスが、移動している血流や心臓壁等の表面で反射された場合の反射波信号は、ドプラ効果により、移動体の超音波送信方向に対する速度成分に依存して、周波数偏移を受ける。
装置本体10に接続される超音波プローブ1は、例えば、複数の圧電振動子が一列で配置された1次元超音波プローブである。1次元超音波プローブとしては、セクタ型、リニア型又はコンベックス型等の超音波プローブが挙げられる。或いは、装置本体10に接続される超音波プローブ1は、例えば、一列に配置された複数の圧電振動子により、被検体Pを2次元で走査するとともに、複数の圧電振動子を所定の角度(揺動角度)で揺動させることで、被検体Pを3次元で走査するメカニカル4Dプローブである。或いは、装置本体10に接続される超音波プローブ1は、例えば、複数の圧電振動子がマトリックス状に配置されることで、被検体Pを3次元で超音波走査することが可能な2Dプローブである。2Dプローブは、超音波を集束して送信することで、被検体Pを2次元で走査することも可能である。
なお、本実施形態に係る超音波診断装置は、後述するように、パルス波(PW:Pulsed Wave)ドプラ法、又は、連続波(CW:Continuous Wave)ドプラ法により、ドプラ波形の収集を行なう。本実施形態において、装置本体10に接続される超音波プローブ1は、Bモード画像データ及びカラードプラ画像データの撮影用の超音波送受信とともに、PWドプラ法によるPWモード、又は、CWドプラ法によるCWモードのドプラ波形の収集用の超音波送受信を実行可能な超音波プローブである。
ここで、上述したように、超音波プローブ1には位置センサ4が取り付けられる。また、上述したように、装置本体10の近傍の任意の位置に、トランスミッター5が配置される。位置センサ4及びトランスミッター5は、超音波プローブ1の位置情報(座標及び角度)を検出するための位置検出システムである。例えば、位置センサ4は、超音波プローブ1に取り付けられる磁気センサである。位置センサ4は、例えば、超音波プローブ1の本体の端部に取り付けられる。また、例えば、トランスミッター5は、自装置を中心として外側に向かって磁場を形成する装置である。
位置センサ4は、トランスミッター5によって形成された3次元の磁場の強度と傾きとを検出する。そして、位置センサ4は、検出した磁場の情報に基づいて、トランスミッター5を原点とする空間における自装置の位置(座標及び角度)を算出し、算出した位置を装置本体10に送信する。ここで、位置センサ4は、自装置が位置する3次元の座標及び角度を、超音波プローブ1の3次元位置情報として、装置本体10に送信する。
なお、本実施形態は、上記の位置検出システム以外のシステムにより、超音波プローブ1の位置情報を取得する場合であっても適用可能である。例えば、本実施形態は、ジャイロセンサや加速度センサ等を用いて、超音波プローブ1の位置情報を取得する場合であっても良い。
入力装置3は、後述するインターフェース部19を介して装置本体10と接続される。入力装置3は、マウス、キーボード、ボタン、パネルスイッチ、タッチコマンドスクリーン、フットスイッチ、トラックボール、ジョイスティック等を有し、超音波診断装置の操作者からの各種設定要求を受け付け、装置本体10に対して受け付けた各種設定要求を転送する。
モニタ2は、超音波診断装置の操作者が入力装置3を用いて各種設定要求を入力するためのGUI(Graphical User Interface)を表示したり、装置本体10において生成された各種画像データ等を表示したりする。
外部装置6は、後述するインターフェース部19及びネットワークを介して装置本体10と接続される装置である。例えば、外部装置6は、各種の医用画像のデータを管理するシステムであるPACS(Picture Archiving and Communication System)のデータベースや、医用画像が添付された電子カルテを管理する電子カルテシステムのデータベースなどである。或いは、外部装置6は、例えば、X線CT(Computed Tomography)装置、MRI(Magnetic Resonance Imaging)装置等、本実施形態に係る超音波診断装置以外の各種医用画像診断装置である。
本実施形態に係る装置本体10は、例えば、DICOM(Digital Imaging and Communications in Medicine)に則った画像フォーマットに統一された各種医用画像のデータを、インターフェース部19を介して外部装置6から取得することができる。例えば、装置本体10は、後述するインターフェース部19を介して、自装置で生成した超音波画像データの比較対象となるボリュームデータ(X線CTボリュームデータや、MRIボリュームデータ等)を、インターフェース部19を介して外部装置6から取得する。
装置本体10は、超音波プローブ1が受信した反射波信号に基づいて超音波画像データを生成する装置である。例えば、本実施形態に係る装置本体10は、2次元の反射波データに基づいて2次元の超音波画像データを生成可能な装置である。また、例えば、本実施形態に係る装置本体10は、3次元の反射波データに基づいて3次元の超音波画像データを生成可能な装置である。以下、3次元の超音波画像データを「超音波ボリュームデータ」と記載する。
装置本体10は、図1に示すように、送受信部11と、Bモード処理部12と、ドプラ処理部13と、画像生成部14と、画像メモリ15と、画像処理部16と、内部記憶部17と、制御部18と、インターフェース部19とを有する。
送受信部11は、パルス発生器、送信遅延部、パルサ等を有し、超音波プローブ1に駆動信号を供給する。パルス発生器は、所定のレート周波数で、送信超音波を形成するためのレートパルスを繰り返し発生する。また、送信遅延部は、超音波プローブ1から発生される超音波をビーム状に集束し、かつ送信指向性を決定するために必要な圧電振動子ごとの遅延時間を、パルス発生器が発生する各レートパルスに対し与える。また、パルサは、レートパルスに基づくタイミングで、超音波プローブ1に駆動信号(駆動パルス)を印加する。すなわち、送信遅延部は、各レートパルスに対し与える遅延時間を変化させることで、圧電振動子面から送信される超音波の送信方向を任意に調整する。
なお、送受信部11は、後述する制御部18の指示に基づいて、所定のスキャンシーケンスを実行するために、送信周波数、送信駆動電圧等を瞬時に変更可能な機能を有している。特に、送信駆動電圧の変更は、瞬間にその値を切り替え可能なリニアアンプ型の発信回路、又は、複数の電源ユニットを電気的に切り替える機構によって実現される。
また、送受信部11は、プリアンプ、A/D(Analog/Digital)変換器、受信遅延部、加算器等を有し、超音波プローブ1が受信した反射波信号に対して各種処理を行って反射波データを生成する。プリアンプは、反射波信号をチャネル毎に増幅する。A/D変換器は、増幅された反射波信号をA/D変換する。受信遅延部は、受信指向性を決定するために必要な遅延時間を与える。加算器は、受信遅延部によって処理された反射波信号の加算処理を行なって反射波データを生成する。加算器の加算処理により、反射波信号の受信指向性に応じた方向からの反射成分が強調され、受信指向性と送信指向性とにより超音波送受信の総合的なビームが形成される。
本実施形態に係る送受信部11は、被検体Pを2次元走査するために、超音波プローブ1から2次元の超音波ビームを送信させる。そして、本実施形態に係る送受信部11は、超音波プローブ1が受信した2次元の反射波信号から2次元の反射波データを生成する。また、本実施形態に係る送受信部11は、被検体Pを3次元走査するために、超音波プローブ1から3次元の超音波ビームを送信させる。そして、本実施形態に係る送受信部11は、超音波プローブ1が受信した3次元の反射波信号から3次元の反射波データを生成する。
ここで、本実施形態に係る超音波診断装置は、超音波プローブ1が1次元超音波プローブであっても、上述した位置検出システムを用いて、3次元の反射波データを生成することが可能である。例えば、操作者は、超音波プローブ1を被検体Pの体表に当接させた状態で、超音波プローブ1の位置及び角度を変更して、複数断面を2次元走査することで、被検体Pを3次元走査する。これにより、送受信部11は、複数断面の2次元の反射波データを生成する。例えば、後述する制御部18は、位置検出システムから得られた超音波プローブ1の3次元位置情報から、複数断面の2次元の反射波データを3次元に配置することで、3次元の反射波データを再構成することができる。
なお、送受信部11からの出力信号の形態は、RF(Radio Frequency)信号と呼ばれる位相情報が含まれる信号である場合や、包絡線検波処理後の振幅情報である場合等、種々の形態が選択可能である。
Bモード処理部12は、送受信部11から反射波データを受信し、対数増幅、包絡線検波処理等を行なって、信号強度が輝度の明るさで表現されるデータ(Bモードデータ)を生成する。
ドプラ処理部13は、送受信部11から受信した反射波データから速度情報を周波数解析し、ドプラ効果による血流や組織、造影剤エコー成分を抽出し、速度、分散、パワー等の移動体情報を多点について抽出したデータ(ドプラデータ)を生成する。本実施形態の移動体は、血管内を流動する血液である。
なお、本実施形態に係るBモード処理部12及びドプラ処理部13は、2次元の反射波データ及び3次元の反射波データの両方について処理可能である。
画像生成部14は、Bモード処理部12及びドプラ処理部13が生成したデータから超音波画像データを生成する。すなわち、画像生成部14は、Bモード処理部12が生成したBモードデータから反射波の強度を輝度にて表したBモード画像データを生成する。また、画像生成部14は、ドプラ処理部13が生成したドプラデータから移動体情報(血流情報や組織の移動情報)を表す平均速度画像データ、分散画像データ、パワー画像データ、又は、これらの組み合わせ画像としてのカラードプラ画像データを生成する。
ここで、画像生成部14は、一般的には、超音波走査の走査線信号列を、テレビ等に代表されるビデオフォーマットの走査線信号列に変換(スキャンコンバート)し、表示用の超音波画像データを生成する。具体的には、画像生成部14は、超音波プローブ1による超音波の走査形態に応じて座標変換を行なうことで、表示用の超音波画像データを生成する。また、画像生成部14は、スキャンコンバート以外に種々の画像処理として、例えば、スキャンコンバート後の複数の画像フレームを用いて、輝度の平均値画像を再生成する画像処理(平滑化処理)や、画像内で微分フィルタを用いる画像処理(エッジ強調処理)等を行なう。また、画像生成部14は、超音波画像データに、種々のパラメータの文字情報、目盛り、ボディーマーク等を合成する。
すなわち、Bモードデータ及びドプラデータは、スキャンコンバート処理前の超音波画像データであり、画像生成部14が生成するデータは、スキャンコンバート処理後の表示用の超音波画像データである。なお、Bモードデータ及びドプラデータは、生データ(Raw Data)とも呼ばれる。
更に、画像生成部14は、Bモード処理部12が生成した1走査線上のBモードデータの時系列データから、Mモード画像データを生成することが可能である。また、画像生成部14は、ドプラ処理部13が生成したドプラデータから、血流の速度情報を時系列に沿ってプロットしたドプラ波形を生成することが可能である。
更に、画像生成部14は、Bモード処理部12が生成した3次元のBモードデータに対して座標変換を行なうことで、3次元のBモード画像データを生成する。また、画像生成部14は、ドプラ処理部13が生成した3次元のドプラデータに対して座標変換を行なうことで、3次元のドプラ画像データを生成する。3次元Bモードデータ及び3次元ドプラデータは、スキャンコンバート処理前のボリュームデータとなる。すなわち、画像生成部14は、「3次元のBモード画像データや3次元のドプラ画像データ」を「超音波ボリュームデータ」として生成する。
更に、画像生成部14は、超音波ボリュームデータをモニタ2にて表示するための各種の2次元画像データを生成するために、超音波ボリュームデータに対してレンダリング処理を行なう。画像生成部14が行なうレンダリング処理としては、断面再構成法(MPR:Multi Planer Reconstruction)を行なって超音波ボリュームデータからMPR画像データを生成する処理がある。また、画像生成部14が行なうレンダリング処理としては、超音波ボリュームデータに対して「Curved MPR」を行なう処理や、超音波ボリュームデータに対して「Maximum Intensity Projection」を行なう処理がある。また、画像生成部14が行なうレンダリング処理としては、3次元の情報を反映した2次元画像データを生成するボリュームレンダリング(VR:Volume Rendering)処理がある。
更に、画像生成部14は、他の医用画像診断装置が収集したボリュームデータに対しても、上記の各種レンダリング処理を行なうことができる。かかるボリュームデータは、X線CT装置により収集された3次元のX線CT画像データ(X線CTボリュームデータ)や、MRI装置により収集された3次元のMRI画像データ(MRIボリュームデータ)である。一例として、画像生成部14は、後述する取得部161が取得した情報に基づいて、現時点で生成した2次元の超音波画像データの走査断面に対応する断面を用いたMPR処理により、X線CTボリュームデータから2次元のX線CT画像データを生成する。
画像メモリ15は、画像生成部14が生成した表示用の画像データを記憶するメモリである。また、画像メモリ15は、Bモード処理部12やドプラ処理部13が生成したデータを記憶することも可能である。画像メモリ15が記憶するBモードデータやドプラデータは、例えば、診断の後に操作者が呼び出すことが可能となっており、画像生成部14を経由して表示用の超音波画像データとなる。
画像処理部16は、本実施形態において、ドプラ波形の収集を支援するために、装置本体10に設置される。画像処理部16は、図1に示すように、取得部161と、抽出部162と、算出部163と、ガイド画像生成部164とを有する。
取得部161は、3次元の医用画像データであるボリュームデータにおいて、超音波プローブ1の走査断面に対応する断面の位置を特定する対応関係を取得する。上記のボリュームデータは、超音波診断装置以外の他種の医用画像診断装置により撮影されたボリュームデータであり、例えば、X線CTボリュームデータである。
取得部161により取得される対応関係は、従来、「同時表示機能」に用いられている。「同時表示機能」は、超音波プローブ1の移動により変更される走査断面の2次元超音波画像データと、当該2次元超音波画像データと位置合わせされたX線CTボリュームデータの2次元X線CT画像データとを、モニタ2の画面にリアルタイムで同時に表示する機能である。図2及び図3は、取得部を説明するための図である。
取得部161は、図2に示すように、超音波プローブ1に取り付けた位置センサ4及びトランスミッター5で構成される位置検出システムから、超音波プローブ1の実空間における3次元位置情報(座標及び角度)を取得する。
「同時表示機能」を行なう場合、超音波プローブ1の3軸(X、Y、Z)と、X線CTボリュームデータ100(図2を参照)の3軸との軸合わせが行われる。例えば、操作者は、位置センサ4が取り付けられた超音波プローブ1を、被検体Pに対して垂直に当接し、セットボタンを押下する。取得部161は、セットボタンが押下された時点で取得した超音波プローブ1の3次元位置情報から、実空間での垂直方向で定まる直行3軸をセットする。
次に、操作者は、図2に示すように、X線CTボリュームデータ100のMPR画像データ101aに描出された特徴部分と同一の特徴部分が描出される超音波画像データ200が表示されるように、超音波プローブ1を移動させて、再度セットボタンを押下する。また、操作者は、MPR画像データ101aの特徴部分と、超音波画像データ200の特徴部分とを、例えば、マウスにより指定する。特徴部分としては、例えば、血管や、剣状突起などが用いられる。
取得部161は、セットボタンが再度押下された時点で取得した超音波プローブ1の3次元位置情報と、X線CTボリュームデータ100における特徴部分の位置情報と、超音波画像データ200における特徴部分の位置情報とから、実空間における「任意の走査断面」の位置と、「任意の走査断面」の超音波走査で生成された超音波画像データにおける位置と、X線CTボリュームデータ100における「任意の走査断面」に対応する断面の位置との対応関係を取得する。
かかる対応関係を用いて、例えば、制御部18は、X線CTボリュームデータ100で、走査断面に対応する断面を特定し、画像生成部14は、制御部18が特定した断面を用いてX線CTボリュームデータ100からMPR画像データを生成することができる。これにより、モニタ2は、移動された走査断面の超音波画像データ200Aと同一断面のMPR画像データ101Aとを同時に表示する。
ここで、上記の対応関係は、「任意の走査断面」の超音波走査で生成された超音波画像データの座標を、X線CTボリュームデータ100の座標に変換することができる。また、上記の対応関係は、X線CTボリュームデータ100の任意の位置の座標を、超音波画像データの座標に変換することができる。
図1に戻って、抽出部162は、ボリュームデータに含まれる血管領域を抽出する。例えば、抽出部162は、血液に対応するCT値を有するボクセルを抽出して、X線CTボリュームデータ100に含まれる血管領域を抽出する。なお、本実施形態は、ボリュームデータに含まれる血管領域を、操作者が血管の位置を手動で設定して抽出する場合であっても良い。
本実施形態において、取得部161が取得した対応関係と、抽出部162が抽出する血管領域とを用いて、算出部163及びガイド画像生成部164が行なう処理については、後に詳述する。
内部記憶部17は、超音波送受信、画像処理及び表示処理を行なうための制御プログラムや、診断情報(例えば、患者ID、医師の所見等)や、診断プロトコルや各種ボディーマーク等の各種データを記憶する。また、内部記憶部17は、必要に応じて、画像メモリ15が記憶する画像データの保管等にも使用される。
更に、内部記憶部17は、外部装置6から転送された各種医用画像の保管にも使用される。具体的には、内部記憶部17は、後述するインターフェース部19を経由して、外部装置6から転送されたDICOM規格のボリュームデータを記憶する。本実施形態では、内部記憶部17は、ドプラ波形が収集される被検体Pの血管を含むボリュームデータ(例えば、X線CTボリュームデータやMRIボリュームデータ)を記憶する。
制御部18は、超音波診断装置の処理全体を制御する。具体的には、制御部18は、入力装置3を介して操作者から入力された各種設定要求や、内部記憶部17から読込んだ各種制御プログラム及び各種データに基づき、送受信部11、Bモード処理部12、ドプラ処理部13、画像生成部14、画像処理部16の処理を制御する。また、制御部18は、画像メモリ15や内部記憶部17が記憶する表示用の超音波画像データをモニタ2にて表示するように制御する。また、制御部18は、画像処理部16の処理結果をモニタ2にて表示するように制御する。
インターフェース部19は、入力装置3、ネットワーク及び外部装置6に対するインターフェースである。入力装置3が受け付けた操作者からの各種設定情報及び各種指示は、インターフェース部19により、制御部18に転送される。また、入力装置3が操作者から受け付けた画像データの転送要求は、インターフェース部19により、ネットワークを介して外部装置6に通知される。また、外部装置6が転送した画像データは、インターフェース部19により、内部記憶部17に格納される。
以上、本実施形態に係る超音波診断装置の全体構成について説明した。かかる構成のもと、本実施形態に係る超音波診断装置は、ドプラ情報に基づく血流速度情報を示すドプラ波形の収集を行なう。
従来、ドプラ波形の収集を行なう場合、操作者は、Bモード、又は、カラードプラモードで、ドプラ波形の収集を行ないたい観察部位の血管が描出される走査断面を探索する。そして、操作者は、2次元超音波画像データ(2次元Bモード画像データ、又は、2次元カラードプラ画像データ)を参照して、ドプラ波形収集用の走査断面(血流情報収集用の走査断面)を決定する。以下では、血流速度情報の収集用に超音波プローブ1で走査される走査断面を「収集断面」と記載する。
そして、操作者は、収集断面において、ドプラ波形の収集を行なう範囲を配置する。以下では、収集断面で、血流情報収集用に配置される範囲を「収集範囲」と記載する。PWモードでは、操作者は、収集断面の2次元超音波画像データに描出された血管の走行に合わせて、血管内の観察部位に収集範囲であるサンプリングゲートを配置する。PWモードでは、サンプリングゲート内の血流速度情報を示すドプラ波形が表示される。
また、CWモードでは、操作者は、収集断面の2次元超音波画像データに描出された観察部位を含む血管を通るように、ライン状のサンプリングマーカを配置する。CWモードでは、サンプリングマーカと同一位置に設定された走査線上の全ての血流速度情報を示すドプラ波形が表示される。図4は、PWモードのドプラ波形収集で、従来行なわれている設定の一例を示す図である。
図4の右図は、モニタ2に表示された収集断面のBモード画像データである。具体的には、図4の右図は、観察対象の血管が長軸断面で描出されている収集断面のBモード画像データである。操作者は、Bモード画像データにて血管の長軸断面内にサンプリングゲートを配置する(図4の右図に示す平行2重線を参照)。そして、操作者は、又は、制御部18は、サンプリングゲートを通過する走査線であり、収集断面内で走査可能な走査線の位置に、ラインマーカを配置する(図4の右図に示す点線を参照)。ここで、ドプラ情報で得られる血流速度情報は、血管の走行方向の血流速度情報ではなく、ラインマーカ上の走査線方向の血流速度情報である。ここで、血流方向(血管の走行方向)と、走査線方向との角度を「θ」とする。ドプラ処理部13は、ドプラ情報から得られた血流速度情報を「1/cosθ」で角度補正することで、血管の走行方向の血流速度情報を求める。
操作者は、サンプリングゲート及びラインマーカを配置するとともに、角度補正を行なうための「θ」をドプラ処理部13に通知するために、血管の走行方向に平行となる角度マーカを、サンプリングゲート内に配置する(図4の右図の平行2重線を斜めに横切る2つの線分を参照)。なお、角度マーカの向きは、操作者により調整可能である。ドプラ処理部13は、ラインマーカと角度マーカとの角度から「θ」を求め、角度補正を行なう。
サンプリングゲート及びラインマーカの位置情報は、ドプラ情報をサンプリングする収集範囲を定義する位置情報となる。また、角度補正用の情報である角度マーカの位置情報も、収集範囲を定義する位置情報として含めることができる。
そして、操作者は、サンプリングゲート及びラインマーカの配置と、角度マーカの調整とが完了した後、PWモードでのドプラ波形の収集開始要求を、入力装置3を介して入力する。これにより、ドプラ波形の収集が開始され、モニタ2は、図4の左図に示すように、ドプラ波形を表示する。
なお、CWモードでは、PWモードのラインマーカに対応するサンプリングマーカ及び角度マーカが配置される。CWモードでは、血管を通るサンプリングマーカの位置情報が、ドプラ情報をサンプリングする収集範囲を定義する位置情報となる。また、CWモードでは、角度補正用の情報である角度マーカの位置情報も、収集範囲を定義する位置情報として含めることができる。
操作者は、ドプラ情報を得るために、2次元超音波画像データを参照しながら、超音波プローブ1の当接位置及び当接角度を最適な位置及び角度に調整して、収集範囲を配置する必要がある。ここで、超音波プローブ1の最適な当接位置及び最適な当接角度とは、例えば、「θ」が最小となる位置及び角度である。血管の走行方向と、走査線方向との角度「θ」が「0度」に近いほど、「1/cosθ」が「1」に近づき、血流速度情報の角度補正による誤差が少なくなる。一方、血管の走行方向と、走査線方向との角度「θ」が大きくなるほど、「1/cosθ」が小さくなり、血流速度情報の角度補正による誤差が大きくなる。
しかし、操作者は、2次元超音波画像データを参照しても3次元空間における血管の走行を把握することが困難である。このため、操作者にとって、血管の走行方向と、走査線方向との角度ができるだけ小さくなるように収集範囲を配置することは、困難な操作となる。また、観察対象の血管が長軸断面で描出されていると操作者が判断した場合であっても、画像データが2次元であることから、参照している断面が実際には長軸断面からずれている場合がある。この場合、角度マーカは、血管の走行方向と平行とならないため、「θ」に誤差が生じ、血流速度情報を正確に求めることが出来ない。
そこで、本実施形態では、最適な血流速度情報を表示するため操作者が行なう操作を容易にするために、算出部163及びガイド画像生成部164は、取得部161が取得した対応関係と、抽出部162が抽出する血管領域とを用いて、以下の処理を行なう。なお、以下では、PWモードで行なわれる処理を中心に説明する。なお、以下の処理が開始される前に、取得部161は、ボリュームデータ(例えば、X線CTボリュームデータ100)の対応関係を取得している。
算出部163は、収集範囲及び収集断面のボリュームデータ(他種ボリュームデータ)における位置情報を、当該ボリュームデータに含まれる血管の走行方向に基づいて算出する。具体的には、算出部163は、収集される血流速度情報の誤差が許容範囲内となる収集範囲及び収集断面のボリュームデータ(他種ボリュームデータ)における位置情報を、当該ボリュームデータに含まれる血管の走行方向に基づいて算出する。すなわち、算出部163は、ドプラ波形を収集するための収集範囲と、収集範囲が配置される収集断面との位置情報を、実空間ではなく、血管領域を抽出可能な他種ボリュームデータの仮想空間において算出する。
そして、ガイド画像生成部164は、対応関係及び位置情報を用いて、収集範囲が配置された収集断面が走査される位置まで超音波プローブ1を操作者が移動するためのガイド画像データを生成する。具体的には、ガイド画像生成部164は、ボリュームデータ(他種ボリュームデータ)から生成された2次元画像データを用いる。この2次元画像データは、例えば、X線CTボリュームデータ100のVR画像データである。VR画像データは、画像生成部14が生成する場合であっても、ガイド画像生成部164が生成する場合であっても良い。そして、ガイド画像生成部164は、この2次元画像データに、超音波プローブ1の現時点での走査断面を示すマーカを対応関係に基づく位置に重畳し、収集範囲及び収集断面を示すマーカを位置情報に基づく位置に重畳する。これにより、ガイド画像生成部164は、ガイド画像データを生成する。ガイド画像データは、操作者が超音波プローブ1の移動操作を行なうための、ナビゲーション用の画像となる。
まず、算出部163は、ボリュームデータにおいて、超音波プローブ1が走査可能な3次元走査範囲に対応する3次元範囲を取得する。ここで、3次元走査範囲は、操作者が被検体Pの体表に当接した超音波プローブ1を可能な範囲で移動させることで得られる。そして、算出部163は、3次元走査範囲と、取得部161が取得済みの対応関係とから、3次元範囲の位置を算出する。図5は、3次元範囲の取得処理の一例を示す図である。
例えば、心臓の超音波走査は、肋間に超音波プローブ1を当接した状態で行なわれる。すなわち、超音波プローブ1は、必ずしも、任意の当接位置及び当接角度で観察部位を含む断面を走査することが出来ない。そこで、例えば、操作者は、被検体Pの体表に超音波プローブ1を当接した状態で、図5の左図に示すように、超音波プローブ1を走査可能な範囲で移動することで、3次元走査範囲を設定する。取得部161は、3次元走査範囲の実空間における位置情報を、位置センサ4から取得し、算出部163に通知する。そして、算出部163は、3次元走査範囲の実空間における位置情報と、対応関係とから、図5の右図に示すように、X線CTボリュームデータ100において、3次元走査範囲に対応する3次元範囲300の位置を算出する。3次元範囲300は、収集範囲及び収集断面をX線CTボリュームデータ100において、探索する範囲となる。3次元範囲300に限定することで、最適なドプラ波形を収集可能な収集範囲及び収集断面の位置情報算出処理に要する負荷は、軽減される。
ただし、3次元範囲300の位置は、対応関係を用いずに算出される場合であっても良い。かかる場合、算出部163は、X線CTボリュームデータ100の被検体Pの体表が位置する情報を用いて、3次元範囲の位置を算出する。図6は、3次元範囲の取得処理の別の一例を示す図である。例えば、算出部163は、空気に対応するCT値を用いて、図6に示すように、X線CTボリュームデータ100における被検体Pの体表400の位置を算出する。体表400は、超音波プローブ1が当接可能な範囲である。例えば、算出部163は、体表400のX線CTボリュームデータ100における位置と、超音波プローブ1が2次元走査する断面の形状及び大きさとに基づいて、3次元範囲300の位置を算出する。
なお、例えば、算出部163は、骨に対応するCT値を用いて、X線CTボリュームデータ100における被検体Pの骨領域の位置を算出し、骨領域の位置を更に用いて、3次元範囲300の位置を算出しても良い。或いは、算出部163は、操作者が指定した体表400の一部範囲の位置を更に用いて、3次元範囲300の位置を算出しても良い。
操作者は、3次元走査範囲の走査を終えた後、超音波プローブ1が3次元走査範囲内の断面を走査している状態で、入力装置3を用いてナビゲーション開始要求を入力する。これにより、算出部163及びガイド画像生成部164は、以下に説明するナビゲーション処理を開始する。図7〜図13は、算出部及びガイド画像生成部が行なうナビゲーション処理を説明するための図である。
図7の右図は、現時点の走査断面のBモード画像データである。また、図7の左図は、上述したX線CTボリュームデータ100から生成されたVR画像データ101である。算出部163は、現時点の走査断面の位置情報を取得部161から取得し、取得した位置情報と対応関係とを用いて、X線CTボリュームデータ100において現時点の走査断面に対応する断面(以下、対応走査断面と記載する)の位置を算出する。ガイド画像生成部164は、図7に示すように、VR画像データ101において、算出部163が算出した位置に対応走査断面に該当する枠500を重畳する。また、画像生成部14は、図7に示すように、現時点の走査断面のBモード画像データに、現時点の走査断面に対応する枠600を重畳させる。例えば、枠500及び枠600は、黄色の実線で描画される。なお、VR画像データ101には、血管が描出されている。
そして、入力装置3は、ボリュームデータにおいて、操作者から観察部位の指定を受け付ける。例えば、操作者は、図8の左図に示すように、VR画像データ101の血管内にマウスのカーソルを移動して、クリックする。VR画像データ101は、3次元情報を有するX線CTボリュームデータ100から生成された画像であるため、算出部163は、操作者が観察部位として指定した位置の3次元位置を算出可能である。
そこで、算出部163は、図9に示すように、VR画像データ101で指定された観察部位のX線CTボリュームデータ100における3次元位置を算出する。そして、抽出部162は、図9に示すように、観察部位の3次元位置に位置する血管領域501を、CT値を用いて抽出する。なお、本実施形態は、抽出部162が血管領域501から血管壁を抽出しても良い。また、本実施形態は、抽出部162が予めX線CTボリュームデータ100の血管領域(又は、血管壁)を全て抽出しておいても良い。
そして、算出部163は、血管領域501を用いて観察部位における血管の走行方向を算出する。例えば、算出部163は、図9に示すように、観察部位における血管走行ベクトル502を算出する。図9では、血管走行ベクトル502を両矢印で示している。
そして、算出部163は、血管の走行方向(血管走行ベクトル502)に沿う走査線の位置及び角度をX線CTボリュームデータ100において算出する。血管の走行方向(血管走行ベクトル502)に沿う走査線は、観察部位で収集される血流速度情報の誤差が最小となる最適走査線となる。観察部位で収集される血流速度情報の誤差が最小となる最適走査線の位置及び角度は、観察部位を通り、血管の走行方向と方向が一致する走査線の位置及び角度である。算出部163は、X線CTボリュームデータ100において最適走査線の位置及び角度を算出する。そして、算出部163は、算出した位置及び角度で定まる線を最適走査線に対応する最適ラインとして設定する。そして、算出部163は、最適ラインを用いて、収集範囲及び収集断面の位置情報を算出する。
具体的には、算出部163は、図10に示すように、観察部位を通り血管走行ベクトル502に一致する線を、最適ライン503として配置する。換言すると、算出部163は、血管走行ベクトル502(血管領域を用いて算出した観察部位における血管の走行方向)に対して平行に、最適ライン503を配置する。最適ライン503と血管走行ベクトル502との角度は「0度」となり、仮に、最適ライン503に対応する走査線で収集される血流速度情報の誤差は、最小となる。換言すると、最適ライン503は、X線CTボリュームデータ100において、最適なラインマーカに対応するマーカとなる。そして、算出部163は、図10に示すように、最適ライン503に垂直に、観察部位を含む2重線504を配置する。2重線504は、X線CTボリュームデータ100において、最適なサンプリングゲートに対応するマーカとなる。最適ライン503及び2重線504は、PWモードにより観察部位のドプラ波形を収集するための最適収集範囲となる。
そして、算出部163は、図10に示すように、観察部位を通り、かつ、観察部位における血管の走行方向に沿った走行ライン505を配置する。換言すると、算出部163は、2重線504に垂直に、走行ライン505を配置する。走行ライン505は、X線CTボリュームデータ100において、高精度に角度補正を行なえる最適な角度マーカに対応するマーカとなる。
そして、算出部163は、収集範囲及び収集断面を探索するために、上述した3次元範囲300内で観察部位の位置を固定した状態で最適ライン503を移動する。なお、PWモードであることから、算出部163は、最適ライン503及び2重線504との相対的位置関係を維持した状態で、最適ライン503を移動する。また、算出部163は、走行ライン505については、移動を行なわない。
そして、算出部163は、移動した最適ライン503を走査可能な走査断面であり、移動した最適ライン503と血管走行ベクトル502との角度が最小となる位置を探索する。算出部163は、上記の探索条件で探索を行なって、位置情報を算出する。これにより、算出部163は、図11に示すように、上記の探索条件を満たす「ライン506及び2重線507」の位置を、収集範囲の位置として算出する。また、算出部163は、図11に示すように、上記の探索条件を満たす「枠508」の位置を、収集断面の位置として算出する。これにより、算出部163は、位置情報を算出する。なお、図11に示すように、走行ライン505の位置は、保持されている。
そして、ガイド画像生成部164は、図8の左図に示す「VR画像データ101及び枠500」に、「ライン506及び2重線507」、「枠508」及び「走行ライン505」それぞれを示すマーカを重畳して、ガイド画像データを生成する。すなわち、ガイド画像生成部164は、図12の左図に示すように、枠508を示す点線枠509を重畳し、ライン506を示す点線510を重畳し、2重線507を示す2重線511を重畳し、走行ライン505を示すライン512を角度補正用の角度マーカとして重畳する。なお、ガイド画像生成部164は、黄色の実線で示す枠500と差別化を行なうため、点線枠509を緑色で描画する。そして、ガイド画像データは、制御部18の制御により、モニタ2に表示される。
すなわち、点線枠509は、収集断面を示すマーカとなる。また、点線510及び2重線511は、収集範囲を示すマーカとなる。操作者は、図12の左図に示すガイド画像データを参照して、超音波プローブ1の移動を開始する。超音波プローブ1が移動されると、図12の右図のBモード画像データが更新されるとともに、対応関係に基づいて、枠500の位置も移動する。操作者は、移動する枠500を参照して、枠500が点線枠509と一致するように、超音波プローブ1を移動する。
更に、制御部18は、ガイド画像データとともに、収集断面が走査される位置まで超音波プローブ1を操作者が移動するために必要となる情報をモニタ2に表示させても良い。例えば、ガイド画像生成部164は、撮影部位に対応する3次元ボディーマークを内部記憶部17から取得する。そして、例えば、ガイド画像生成部164は、図13に示すように、3次元ボディーマーク上において、2つのプローブ模式図を配置する。ガイド画像生成部164は、一方のプローブ模式図を、ナビゲーション開始要求時の超音波プローブ1の位置に対応する位置に配置し、他方のプローブ模式図を、収集断面走査時の超音波プローブ1の位置に対応する位置に配置する。そして、例えば、ガイド画像生成部164は、図13に示すように、収集断面が走査される位置の操作として、平行移動を示す矢印と、回転移動を示す矢印とを配置する。そして、制御部18は、図13に示す画像データをモニタ2に表示させる。
なお、ガイド画像生成部164は、枠500の位置に、Bモード画像データを重畳したガイド画像データを生成表示しても良い。また、本実施形態は、ガイド画像データに用いられるX線CTボリュームデータ100の2次元画像データが、MPR画像データであっても良い。かかる場合、操作者は、複数のMPR画像データを用いた複数のガイド画像データをサムネール表示させたり、動画表示させたりしても良い。また、本実施形態は、ガイド画像データを参照した操作者が、収集範囲及び収集断面の位置を調整しても良い。
そして、制御部18は、収集断面が走査される位置まで操作者が超音波プローブ1を移動した場合、収集断面が走査されている旨を操作者に報知する情報を出力させる。また、収集断面が走査される位置まで操作者が超音波プローブ1を移動した場合、制御部18は、収集断面の超音波画像データに収集範囲を重畳表示させる。図14は、ナビゲーション完了時の表示形式を説明するための図である。
超音波プローブ1による走査断面が、点線枠509に対応する位置となった場合、すなわち、枠500が点線枠509と重なった場合、ガイド画像生成部164は、図14に示すように、枠500を削除し、緑色の点線枠509を、緑色の実線の枠509’に変更する。これにより、操作者は、収集断面が走査されていることを認識できる。なお、同時に、ガイド画像生成部164又は画像生成部14は、図14に示すように、黄色の枠600を、緑色の枠600’に変更しても良い。これによっても、操作者は、収集断面が走査されていることを認識できる。
また、枠500が点線枠509と重なった場合、ガイド画像生成部164又は画像生成部14は、図14に示すように、収集断面のBモード画像データに、「ガイド画像データにおける収集範囲のマーカである点線510及び2重線511」に対応する「ラインマーカ及びサンプリングゲート」を重畳する。更に、ガイド画像生成部164又は画像生成部14は、図14に示すように、「ガイド画像データにおける角度補正用のマーカであるライン512」に対応する「角度マーカ」を重畳する。なお、本実施形態は、収集断面のBモード画像データに重畳された「ラインマーカ及びサンプリングゲート」を、操作者が調整しても良い。かかる場合、ガイド画像データにおける各マーカの位置を連動して移動しても良い。
そして、制御部18は、収集範囲での血流速度情報(ドプラ波形)の収集を開始させる。なお、ドプラ波形収集時には、ドプラ処理部13は、Bモード画像データに重畳された角度マーカとラインマーカとの角度を用いない。ドプラ処理部13は、走行ライン505とライン506との角度を角度補正用の「θ」として用いる。
ドプラ波形の表示形式は、操作者により選択される場合であっても良いし、初期設定されている場合であっても良い。以下、図15〜図18を用いて、本実施形態で実行可能なドプラ波形の表示形式について説明する。図15〜図18は、本実施形態に係るドプラ波形の表示形式を説明するための図である。
図15に示す表示形式は、図14に示すナビゲーション完了時の表示形式に更にドプラ波形を表示させるものである。図15では、モニタ2において、ガイド画像データと、Bモード画像データが並列表示され、ガイド画像データ及びBモード画像データの下に、ドプラ波形が表示される。
図16に示す表示形式は、図15に示す表示形式と比較して、VR画像データ101を用いたガイド画像データの代わりに、収集断面を含む断面でX線CTボリュームデータ100を切断したMPR画像データに、収集断面に対応するマーカ、収集範囲に対応するマーカ及び角度マーカに対応するマーカを重畳した画像データを表示させるものである。
図17に示す表示形式は、ガイド画像データの表示領域に、ドプラ波形を表示させて、ドプラ波形とBモード画像データとを並列表示させるものである。図18に示す表示形式は、図17の表示形式で表示されていたBモード画像データの一部に、カラードプラモードのカラードプラ画像データを重畳させるものである。なお、本実施形態は、更に、3次元ボディーマークを表示する表示形式が行なわれても良い。
なお、上記では、PWモードにおいて行なわれるナビゲーションについて説明した。一方、CWモードでナビゲーションが行なわれる場合、例えば、3次元範囲300内で、最適ライン503を移動して、図11に示すライン506を収集範囲とする点以外、上記の方法が適用可能である。
次に、図19を用いて、本実施形態に係る超音波診断装置の処理について説明する。図19は、本実施形態に係る超音波診断装置の処理の一例を説明するためのフローチャートである。なお、図19では、PWモードを行なうための収集範囲及び収集断面の位置情報算出用の他種ボリュームデータの対応関係を取得した後の処理について説明する。
図19に示すように、本実施形態に係る制御部18は、操作者からナビゲーション開始要求を受け付けたか否かを判定する(ステップS101)。ここで、ナビゲーション開始要求を受け付けない場合(ステップS101否定)、制御部18は、ナビゲーション開始要求を受け付けるまで待機する。
一方、ナビゲーション開始要求を受け付けた場合(ステップS101肯定)、制御部18の指示に基づくガイド画像生成部164の処理により、ボリュームデータのVR画像データに、現時点の走査断面のマーカを重畳表示する(ステップS102、図7を参照)。そして、算出部163は、VR画像データにおいて、観察部位の指定を受け付けたか否かを判定する(ステップS103)。ここで、観察部位の指定を受け付けない場合(ステップS103否定)、算出部163は、観察部位の指定を受け付けるまで待機する。
一方、観察部位の指定を受け付けた場合(ステップS103肯定)、算出部163は、指定された観察部位のボリュームデータにおける位置を算出し、観察部位における血管の走行方向を算出する(ステップS104、図9を参照)。そして、算出部163は、ボリュームデータにおける最適ラインの位置を算出する(ステップS105、図10を参照)。そして、算出部163は、最適ラインに2重線(サンプリングゲートに対応)及び走行ライン(角度マーカに対応)を配置する(ステップS106、図10を参照)。
そして、算出部163は、走行ラインの位置は保持したまま、観察部位を中心に2重線及び最適ラインを3次元範囲内で移動して、移動した最適ラインを走査可能な走査断面であり、移動した最適ラインと血管走行方向との角度が最小となる収集範囲及び収集断面を探索する(ステップS107、図11を参照)。そして、算出部163は、探索した収集範囲及び収集断面の位置情報を算出する(ステップS108)。
そして、ガイド画像生成部164は、対応関係と位置情報とを用いてガイド画像データを生成し、モニタ2は、ガイド画像データを表示する(ステップS109、図12を参照)。操作者は、ガイド画像データを参照して、超音波プローブ1の移動を開始する。
そして、制御部18は、走査断面が収集断面と一致したか否かを判定する(ステップS110)。ここで、走査断面が収集断面と一致していない場合(ステップS110否定)、制御部18の制御により、新たな走査断面のガイド画像データの生成表示がステップS109で行なわれる。
一方、制御部18は、走査断面が収集断面と一致した場合(ステップS110肯定)、制御部18は、収集断面が走査されている旨を通知する(ステップS111、図14を参照)。そして、制御部18は、収集断面の超音波画像データに、ラインマーカ、サンプリングゲート及び角度マーカを重畳表示させる(ステップS112)。
そして、制御部18は、ドプラ波形の収集及び表示を開始させ(ステップS113)、ナビゲーション処理を終了する。
上述してきたように、本実施形態では、血管領域を抽出可能なボリュームデータを用いて算出して観察部位の血管走行方向を用いて、3次元範囲内で最適な収集範囲及び収集断面の位置情報を算出し、ガイド画像データの生成表示を行なう。従って、本実施形態では、最適な血流速度情報を表示するため操作者が行なう操作を容易にすることができる。
なお、上記の実施形態では、走査可能な3次元領域として設定された3次元走査範囲に対応するボリュームデータの3次元範囲内で、最適な探索条件により、1セットの収集範囲及び収集断面を探索する場合について説明した。しかし、収集断面は、あくまでも、ボリュームデータの仮想空間で探索された断面である。このため、操作者が実空間の収集断面を走査しようとしても、収集断面を走査できない場合や、実空間の収集断面を走査されても、Bモード画像データに骨等に起因する陰影が発生する場合がある。
かかる場合に対応するため、本実施形態は、以下に説明する第1変形例を行なっても良い。第1変形例では、算出部163は、3次元範囲内で移動した最適ラインと観察部位における血管の走行方向との角度が所定値以下となる位置を探索する。例えば、所定値は、「60度」である。
そして、算出部163は、候補収集範囲及び候補走査断面の候補セットを複数求めて、複数の候補位置情報を算出する。ここで、制御部18は、例えば、角度が小さい順に複数の候補位置情報の順位を決定する。図20は、第1変形例を説明するための図である。
例えば、制御部18は、図20の(A)に示すように、所定値「THθ=60度」である場合、「角度:45度」の候補位置情報を第1候補とし、「角度:48度」の候補位置情報を第2候補とし、「角度:50度」の候補位置情報を第3候補とする。
そして、制御部18は、複数の候補位置情報の中で所定順位の候補位置情報に基づくガイド画像データの変更要求を操作者から受け付けた場合、所定順位の次の順位の候補位置情報に基づくガイド画像データを表示させる。例えば、第1候補の収集断面で適切なBモード画像が表示されなかった場合、操作者は、入力装置3が有する「Next」ボタンを押下する。かかる場合、制御部18は、第1候補の位置情報に基づくガイド画像データを表示させる。なお、制御部18は、候補の順位が下がったことを通知するため、第2候補のガイド画像データに重畳されるマーカの色を、第1候補のガイド画像データに重畳されるマーカの色から変更させても良い。
ここで、制御部18は、角度の値を用いる方法以外の方法で、順位付けを行なっても良い。例えば、制御部18は、角度が最小となる候補位置情報を第1位とし、当該第1位の候補位置情報からの移動操作が少ない順に、残余の候補位置情報の順位を決定する。或いは、例えば、制御部18は、現時点での超音波プローブ1の位置からの移動操作が少ない順に複数の候補位置情報の順位を決定する。
また、所定値を用いて1つの候補セットしか抽出されていない場合がある。かかる状態で、この候補セットの候補位置情報に基づくガイド画像データの変更要求を操作者から受け付けた場合、算出部163は、所定値を大きくした値以下となる候補セットの位置を再度探索する。例えば、図20の(B)に示すように、「THθ=60度」であり、「角度:45度」の第1候補しか探索されてない場合、算出部は、「THθ’=65度」に変更する。これにより、算出部163は、例えば、図20の(B)に示すように、「角度:61度」の位置情報が得られる候補収集範囲及び候補走査断面の候補セットを第2候補とする。そして、制御部18は、「角度:61度」の位置情報に基づくガイド画像データを表示させる。
なお、本実施形態は、所定値を用いて候補セットが探索されなかった場合に、所定値を大きくした値以下となる候補セットの位置を再度探索しても良い。或いは、本実施形態は、所定値を用いて候補セットが探索されなかった場合は、観察部位を走査可能な新たな3次元走査範囲の設定を要求する旨を操作者に通知しても良い。
また、上記の実施形態では、X線CTボリュームデータやMRIボリュームデータ等、血管領域を抽出可能な他種ボリュームデータを用いる場合について説明した。しかし、本実施形態は、超音波送受信により撮影され血流が存在する領域が抽出可能なボリュームデータが用いられる場合であっても良い。この第2変形例では、かかるボリュームデータは、カラーフローマッピング(CFM:Color Flow Mapping)法で3次元走査して得られるカラードプラボリュームデータや、パワーボリュームデータが用いられる。
例えば、操作者は、図5の左図に示す3次元走査範囲を設定する際に、カラードプラモードで走査を行なう。例えば、制御部18は、取得部161が位置センサ4から取得した超音波プローブ1の3次元位置情報から、複数断面それぞれの2次元反射波データを3次元反射波データに再構成し、ドプラ処理部13に送信する。図21は、第2変形例を説明するための図である。
これにより、画像生成部14は、例えば、図21に示すように、パワーボリュームデータ700を生成する。抽出部162は、パワー値に応じた輝度値が割り当てられたボクセルを抽出することで、図21に示すように、血流領域710を抽出する。血流領域710は、血管の走行方向を算出可能である。なお、第2変形例では、超音波プローブ1による走査断面の位置と、パワーボリュームデータ700の位置との対応関係は、位置センサ4から取得される超音波プローブ1の3次元位置情報のみで取得することができる。
また、画像処理部16が行なうナビゲーション処理は、Mモード撮影に適用可能である。例えば、Mモードは、走査線上の心臓弁や心筋の運動情報を得るために行なわれる。図22は、第3変形例を説明するための図である。例えば、第3変形例は、図22に示すように、Mモード用サンプリングラインを最適な位置及び角度で配置する場合に適用される。この第3変形例では、Mモード用サンプリングラインが最適な位置及び角度となる条件は、例えば、正常な心臓弁が移動する方向や、正常な心筋が移動する方向を、心臓ボリュームデータのVR画像データ上で、操作者が指定することで設定可能である。
なお、上記の実施形態、第1変形例、第2変形例及び第3変形例において説明した各処理のうち、自動的に行なわれるものとして説明した処理の全部または一部を手動的に行なうこともでき、或いは、手動的に行なわれるものとして説明した処理の全部または一部を公知の方法で自動的に行なうこともできる。この他、上記文書中や図面中で示した処理手順、制御手順、具体的名称、各種のデータやパラメータを含む情報については、特記する場合を除いて任意に変更することができる。
また、上記の説明で図示した各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部または一部を、各種の負荷や使用状況等に応じて、任意の単位で機能的または物理的に分散・統合して構成することができる。さらに、各装置にて行なわれる各処理機能は、その全部または任意の一部が、CPU及び当該CPUにて解析実行されるプログラムにて実現され、あるいは、ワイヤードロジックによるハードウェアとして実現され得る。
また、上記の実施形態、第1変形例、第2変形例及び第3変形例で説明した画像処理方法は、あらかじめ用意された画像処理プログラムをパーソナルコンピュータやワークステーション等のコンピュータで実行することによって実現することができる。このプログラムは、インターネット等のネットワークを介して配布することができる。また、このプログラムは、ハードディスク、フレキシブルディスク(FD)、CD−ROM、MO、DVD等のコンピュータで読み取り可能な記録媒体に記録され、コンピュータによって記録媒体から読み出されることによって実行することもできる。
以上、説明したとおり、実施形態、第1変形例、第2変形例及び第3変形例によれば、最適な血流速度情報を表示するため操作者が行なう操作を容易にすることができる。
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。