JP5985007B1 - 超音波診断装置 - Google Patents

Abstract

【課題】三次元超音波画像の動的な表示を改善可能な装置を提供する。【解決手段】三次元画像形成部２０は、超音波を立体的に送受して得られた複数時相のボリュームデータに基づいて、それら複数時相に対応した複数フレームの三次元超音波画像を形成する。中間画像形成部５０は、複数フレームの三次元超音波画像の各フレーム間における画像の移動ベクトルに基づいて、各フレーム間に追加される１又は複数フレームの三次元中間画像を形成する。表示画像形成部８０は、複数フレームの三次元超音波画像とそれらの複数フレーム間に追加された複数フレームの三次元中間画像に基づいて三次元動画像を形成し、三次元動画像が表示部８２に表示される。【選択図】図１

Description

本発明は、超音波診断装置に関し、特に、三次元超音波画像を形成する装置に関する。

診断対象を含む空間に超音波を立体的に送受することにより得られるボリュームデータに基づいて、診断対象を立体的に表現した三次元超音波画像を形成する超音波診断装置が知られている。例えば、母体内から得られたボリュームデータに対してレンダリング処理を行うことにより、胎児を立体的に映し出した三次元超音波画像としてボリュームレンダリング画像を得ることができる。

三次元超音波画像を形成するにあたっては、電子的走査または電子的走査と機械的走査の組み合わせにより、超音波ビームが２次元的に走査され、走査空間（ボリューム）内から、ボリュームデータを構成する多数のデータ（ボクセルデータ）が収集される。そのため、超音波ビームを１次元的に走査することにより走査面に対応した断層画像を形成するＢモード画像と比較して、三次元超音波画像の形成には多くの送受信時間を要し、例えば三次元超音波画像を動的に表示する際の一つの課題となる。

ちなみに、テレビジョン映像などを代表例とする一般的な画像に係る画像処理の分野においては、画像を動的に表示する際のフレームレートを改善する技術が提案されている。例えば、特許文献１，２には、画像の動きベクトルを利用して補間フレームを生成することにより画像のフレーム数を増加させてフレームレートを高める技術（ＦＲＣ：フレームレートコントロール）が開示されている。

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

本発明は、上述した背景技術に鑑みて成されたものであり、その目的は、三次元超音波画像の動的な表示を改善することにある。

本発明の具体例として好適な超音波診断装置は、超音波を立体的に送受して得られた複数時相のボリュームデータに基づいて、当該複数時相に対応した複数フレームの三次元超音波画像を形成する超音波画像形成部と、前記複数フレームの三次元超音波画像の各フレーム間における画像の移動情報に基づいて、当該各フレーム間に追加される１又は複数フレームの三次元中間画像を形成する中間画像形成部と、前記複数フレームの三次元超音波画像とそれらの複数フレーム間に追加された前記複数フレームの三次元中間画像に基づいて三次元動画像を表示する表示処理部と、を有することを特徴とする。

上記構成においては、超音波を立体的に送受することにより、例えば複数時相に亘って各時相ごとにボリュームデータが得られる。各時相ごとに得られるボリュームデータは、例えば、超音波が送受される空間（ボリューム）内から収集される複数のデータ（ボクセルデータ）で構成される。

そして、各時相ごとに得られるボリュームデータに基づいて、その時相に対応した各フレームの三次元超音波画像が形成される。三次元超音波画像の好適な具体例は、ボリュームデータに対するレンダリング処理により得られるレンダリング画像であるが、三次元超音波画像の具体例はレンダリング画像に限定されない。

さらに、上記構成においては、複数フレームの三次元超音波画像の各フレーム間に、１又は複数フレームの三次元中間画像が追加される。例えば、三次元超音波画像の互いに隣接する一方フレームと他方フレームの間（各フレーム間）に、１フレーム以上の三次元中間画像が追加される。三次元超音波画像の各フレーム間に追加される三次元中間画像のフレーム数は、例えば、目標とすべきフレームレートに応じて決定される。例えば、複数フレームの三次元超音波画像のフレームレートを２倍にしたいのであれば、三次元超音波画像の各フレーム間に１フレームの三次元中間画像が追加され、フレームレートを４倍にしたいのであれば、三次元超音波画像の各フレーム間に３フレームの三次元中間画像が追加される。つまり、複数フレームの三次元超音波画像のフレームレートをＭ倍（Ｍは２以上の整数）にしたいのであれば、三次元超音波画像の各フレーム間に（Ｍ−１）フレームの三次元中間画像を追加すればよい。

そして、上記構成においては、複数フレームの三次元超音波画像とそれらの複数フレーム間に追加された複数フレームの三次元中間画像に基づいて三次元動画像が表示される。三次元動画像は、三次元超音波画像と三次元中間画像を動的に再生する画像であり、例えば、複数フレームの三次元超音波画像と追加された複数フレームの三次元中間画像を、時相順に次々に表示することにより三次元動画像が実現される。なお、三次元動画像は、実時間に合わせて複数時相に対応した画像を次々に表示させる通常再生（１倍速再生）の他に、スロー再生（低倍速再生）やコマ送り再生により実現されてもよいし、時相順を逆順とする巻き戻し再生等が実現されてもよい。

上記構成の装置によれば、複数フレームの三次元超音波画像とそれらの複数フレーム間に追加された複数フレームの三次元中間画像に基づいて三次元動画像が表示されるため、例えば複数フレームの三次元超音波画像のみに基づいて動画像を表示する場合に比べて、表示される画像のフレームレートが高められる。

ちなみに、特許文献１，２には、テレビジョン映像などを代表例とする一般的な画像のフレームレートを高める技術が開示されているに過ぎず、三次元超音波画像に関する具体的な開示もなければ、三次元超音波画像を対象とすることを示唆する内容すらない。

望ましい具体例において、前記超音波診断装置は、前記移動情報の信頼性の程度に応じて、単位時間あたりのボリュームデータ数であるボリュームレートを制御する制御部をさらに有する、ことを特徴とする。

上記構成において、移動情報の信頼性の程度は、例えば連続的な度合い（連続的に程度を示す評価値など）によって決定されてもよいし、例えば段階的な度合い（「高，中，低」または「高，低」など）によって決定されてもよい。移動情報の信頼性の程度が連続的な度合いによって決定されていれば、その連続的な度合いに応じてボリュームレートが連続的な値となるように制御することができる。移動情報の信頼性の程度が段階的な度合いによって決定されていれば、例えば、その段階的に度合いに応じてボリュームレートも段階的な値となるように制御することができる。

三次元中間画像は移動情報に基づいて得られるため、三次元中間画像の信頼性は移動情報の信頼性に依存する。そこで、移動情報の信頼性の程度に応じて、単位時間あたりのボリュームデータ数であるボリュームレートが制御される。例えば、各時相ごとに比較的高密度にボリュームデータが収集されるとボリュームレートは比較的低くなってしまうものの、移動情報の信頼性が比較的高い場合には三次元中間画像にも信頼性が認められると判断して、比較的低いボリュームレートによる比較的高密度なボリュームデータの収集が維持される。その一方、例えば、移動情報の信頼性が比較的低い場合には三次元中間画像の信頼性が低いことが予想されるため、ボリュームレートを高めて三次元中間画像の信頼性を高める制御が実現される。これにより、例えば、胎児の三次元動画像を得る場合に、移動情報の信頼性が比較的高い場合には、低いボリュームレートによる高密度なボリュームデータの収集を行い、例えば、胎児の動きが激しくなるなどの理由により移動情報の信頼性が低くなった場合に、ボリュームレートを高めることにより、移動情報の信頼性を向上させて例えば胎児の激しい動きに対応することなどが可能になる。

望ましい具体例において、前記制御部は、信頼性の判定条件に基づいて前記移動情報の信頼性が低いと判定された場合に、前記ボリュームレートを増加させて前記三次元超音波画像の単位時間あたりのフレーム数であるフレームレートを増加させる、ことを特徴とする。

望ましい具体例において、前記制御部は、超音波を立体的に送受するボリューム内における空間フレーム数と、各空間フレームを構成するビームライン数と、超音波のパルス繰り返し周期と、を含む複数の送受信パラメータのうちの少なくとも１つを変更することにより、前記ボリュームレートを制御する、ことを特徴とする。

望ましい具体例において、前記中間画像形成部は、前記移動情報として、前記各フレームの三次元超音波画像を構成する複数の画像要素の各画像要素ごとに各フレーム間における相関演算に基づいて移動ベクトルを導出し、複数の画像要素とそれらの移動ベクトルに基づいて当該各フレーム間に対応した複数の中間画像要素を得ることにより、複数の中間画像要素で構成された前記三次元中間画像を形成する、ことを特徴とする。

上記構成において、移動ベクトルの導出には、例えば、各フレーム間における画像同士のパターンマッチング処理を利用することが望ましいものの、勾配法、位相限定相関法などの他の公知の手法が利用されてもよい。

望ましい具体例において、前記超音波診断装置は、前記複数の画像要素の各画像要素ごとに前記相関演算により得られる相関値の大きさに基づいて、前記移動情報として導出される移動ベクトルの信頼性の程度を示す評価値を算出する評価部をさらに有し、前記制御部は、移動ベクトルの信頼性を示す前記評価値に基づいて、単位時間あたりのボリュームデータ数であるボリュームレートを制御する、ことを特徴とする。

本発明により、三次元超音波画像の動的な表示が改善される。例えば、本発明の好適な具体例によれば、複数フレームの三次元超音波画像とそれらの複数フレーム間に追加された複数フレームの三次元中間画像に基づいて三次元動画像が表示されるため、例えば複数フレームの三次元超音波画像のみに基づいて動画像を表示する場合に比べて、表示される画像のフレームレートが高められる。

本発明の実施形態として好適な超音波診断装置の全体構成図である。 三次元超音波画像と三次元中間画像の具体例を示す図である。 移動ベクトルを導出する処理の具体例を示す図である。 移動ベクトルに対する平滑化処理の具体例を示す図である。 三次元中間画像を形成する処理の具体例１を示す図である。 中間ベクトルの具体例１を示す図である。 三次元中間画像を形成する処理の具体例２を示す図である。 中間ベクトルの具体例２を示す図である。 ボリュームレート制御の具体例を示す図である。

図１は、本発明の実施形態として好適な超音波診断装置の全体構成図である。プローブ１０は、診断対象を含む三次元空間内に超音波を送受する超音波探触子である。プローブ１０は、複数の振動素子を備えており、各振動素子が送受信部１２から得られる送信信号に応じて三次元空間に超音波を送波する。また、三次元空間から超音波の反射波（エコー）を受波した各振動素子がその反射波に応じた受波信号を送受信部１２に出力する。

送受信部１２は、プローブ１０が備える複数の振動素子の各々に対応した送信信号を出力してプローブ１０を送信制御する。その送信制御により、超音波の送信ビームが形成され、三次元空間内で送信ビームが走査される。

ビームデータ処理部１４は、プローブ１０が備える複数の振動素子に対応した複数の受波信号を送受信部１２から得て、それら複数の受波信号に対して整相加算処理などのビーム形成処理を施す。これにより、超音波の受信ビームが形成されて三次元空間内で走査される。つまり、受信ビームのビームアドレスを異ならせながら、三次元空間内で受信ビームが走査され、ビームデータ処理部１４は、複数のビームアドレスに対応した複数のラインデータを形成する。各ラインデータは、複数のエコーデータで構成される。

こうして、三次元空間内で超音波ビーム（送信ビームとそれに対応した受信ビーム）が立体的に走査され、三次元空間内から複数のエコーデータが得られる。なお、プローブ１０は、超音波ビームを三次元空間内において走査して立体的にエコーデータを収集する３Ｄプローブである。例えば、一次元的に配列された複数の振動素子（１Ｄアレイ振動子）によって電子的に形成される走査面を機械的に動かすことにより超音波ビームが三次元的に走査される。また、二次元的に配列された複数の振動素子（２Ｄアレイ振動子）を電子的に制御して超音波ビームが三次元的に走査されてもよい。

ボリュームデータ記憶部１６には、三次元空間内から得られた複数のエコーデータに基づくボリュームデータが記憶される。例えば、ビームデータ処理部１４が複数のエコーデータに対して座標変換処理等を施し、後段の処理において好適な座標系にデータを変換してから、変換後のデータがボリュームデータとしてボリュームデータ記憶部１６に記憶される。

例えば、ビームデータ処理部１４は、超音波の走査座標系（例えばｒθφ座標系）で得られた複数のエコーデータに対して、三次元の座標変換処理と補間処理などを施すことにより、三次元の直交座標系（ｘｙｚ座標系）に対応した複数のボクセルデータを形成し、複数のボクセルデータで構成されるボリュームデータがボリュームデータ記憶部１６に記憶される。

なお、ビームデータ処理部１４が複数のエコーデータに対して二次元の座標変換処理を施すことにより、例えば、走査座標系のうち走査面に対応したｒθ座標系が直交座標系（ｘｙ座標系）に変換されて、変換後のデータがボリュームデータとしてボリュームデータ記憶部１６に記憶されてもよい。また、三次元空間内において超音波を立体的に走査することにより得られた複数のエコーデータが、超音波の立体的な走査に対応した走査座標系（例えばｒθφ座標系）に対応したアドレスを付され、ボリュームデータとしてボリュームデータ記憶部１６に記憶されてもよい。

ボリュームデータは、複数時相に亘って各時相ごとに次々に形成されてボリュームデータ記憶部１６に記憶される。

三次元画像形成部２０は、ボリュームデータ記憶部１６に記憶されたボリュームデータに基づいて三次元超音波画像を形成する。三次元画像形成部２０は、胎児等の診断対象を含む三次元空間に対応したボリュームデータに基づいて、その診断対象を立体的に映し出した三次元超音波画像を形成する。三次元超音波画像の好適な具体例は、公知のボリュームレンダリング処理により得られるレンダリング画像である。

ボリュームレンダリング処理においては、三次元空間に対応したボリュームデータの外側に演算上の仮想的な視点ＶＰが設定され、ボリュームデータを間に挟んで、視点ＶＰと反対側に演算上の二次元平面としてのスクリーンが仮想的に設定される。その視点ＶＰを基準として複数のレイ（透視線）が定義される。各レイは、例えば視点ＶＰを起点としてボリュームデータを貫通してからスクリーン上に達するように設定される。これにより、各レイ上またはそのレイの近傍において、そのレイに対応した複数ボクセルのボクセルデータが対応することになる。そして、各レイごとに、視点ＶＰ側から、そのレイに対応した複数ボクセルに対してレンダリング法に基づくボクセル演算を逐次的に実行すると、最終のボクセル演算の結果としてそのレイに対応した画素値が決定される。そして、複数のレイから得られる複数の画素値をスクリーン上にマッピングすることによりレンダリング画像が得られる。

ボリュームデータ記憶部１６には、複数時相のボリュームデータが記憶されており、三次元画像形成部２０は、各時相ごとに、その時相に対応したボリュームデータに基づいてその時相に対応した三次元超音波画像（レンダリング画像）を形成する。これにより、複数時相に対応した複数フレームの三次元超音波画像が次々に形成される。そして、中間画像形成部５０により、複数フレームの三次元超音波画像の各フレーム間に追加される１又は複数フレームの三次元中間画像が形成される。

図２は、三次元超音波画像と三次元中間画像の具体例を示す図である。図２（Ａ）には複数フレームの三次元超音波画像（３Ｄ超音波画像）の具体例が図示されている。図２（Ａ）には、複数フレームのうちのフレーム１〜５までの３Ｄ超音波画像が代表的に図示されている。そして、複数フレームの３Ｄ超音波画像の各フレーム間に、１又は複数フレームの三次元中間画像（３Ｄ中間画像）が追加されて、三次元動画像（３Ｄ動画像）用のフレーム列が形成される。

図２（Ｂ）には、３Ｄ動画像用のフレーム列の具体例が図示されている。図２（Ｂ）に示す具体例では、３Ｄ超音波画像の各フレーム間に４フレームの３Ｄ中間画像が追加されている。つまり、３Ｄ超音波画像のフレーム１とフレーム２の間に４フレーム（フレーム数４）の３Ｄ中間画像が追加されており、３Ｄ超音波画像のフレーム２とフレーム３の間と、フレーム３とフレーム４の間と、フレーム４とフレーム５の間にも、それぞれ、４フレームの３Ｄ中間画像が追加されている。

３Ｄ中間画像を追加することにより、追加前のフレームレートよりも追加後のフレームレートが高められる。つまり、３Ｄ超音波画像のみの場合における単位時間あたりのフレーム数（フレームレート）よりも、３Ｄ超音波画像に３Ｄ中間画像が追加された３Ｄ動画像用フレーム列における単位時間あたりのフレーム数（フレームレート）が多く（高く）なる。例えば、図２に示す具体例においては、図２（Ａ）に示す３Ｄ超音波画像のみのフレームレートと比較して、図２（Ｂ）に示す３Ｄ動画像用フレーム列のフレームレートは５倍となる。

例えば、複数フレームの３Ｄ超音波画像のフレームレートをＭ倍（Ｍは２以上の整数）にしたいのであれば、３Ｄ超音波画像の各フレーム間に（Ｍ−１）フレームの３Ｄ中間画像を追加すればよい。

図１に戻り、表示画像形成部８０は、三次元画像形成部２０において形成された複数フレームの三次元超音波画像（３Ｄ超音波画像）と、中間画像形成部５０において形成された複数フレームの三次元中間画像（３Ｄ中間画像）に基づいて、三次元動画像用のフレーム列（図２の３Ｄ動画像用フレーム列）を形成する。そして、表示画像形成部８０は、三次元動画像用のフレーム列（３Ｄ動画像用フレーム列）に基づいて、三次元の超音波画像を動的に示した三次元動画像（３Ｄ動画像）を表示部８２に表示させる。

なお、三次元中間画像（３Ｄ中間画像）の形成に利用される移動ベクトルは、移動ベクトル演算部３０において導出される。また、移動ベクトル評価部４０において、移動ベクトルの信頼性の程度が評価される。

制御部１００は、図１の超音波診断装置内を全体的に制御する。制御部１００による全体的な制御には、操作デバイス９０を介して、医師や検査技師などのユーザから受け付けた指示も反映される。

図１に示す構成（符号を付された各部）のうち、送受信部１２，ビームデータ処理部１４，三次元画像形成部２０，移動ベクトル演算部３０，移動ベクトル評価部４０，中間画像形成部５０，表示画像形成部８０の各部は、例えば電気電子回路やプロセッサ等のハードウェアを利用して実現することができ、その実現において必要に応じてメモリ等のデバイスが利用されてもよい。また、上記各部に対応した機能の少なくとも一部がコンピュータにより実現されてもよい。つまり、上記各部に対応した機能の少なくとも一部が、ＣＰＵやプロセッサやメモリ等のハードウェアと、ＣＰＵやプロセッサの動作を規定するソフトウェア（プログラム）との協働により実現されてもよい。

ボリュームデータ記憶部１６は、例えば、半導体メモリやハードディスクドライブなどの記憶デバイスによって実現することができる。表示部８２の好適な具体例は、液晶ディスプレイ等である。操作デバイス９０は、例えばマウス、キーボード、トラックボール、タッチパネル、その他のスイッチ類等のうちの少なくとも一つにより実現できる。そして制御部１００は、例えば、ＣＰＵやプロセッサやメモリ等のハードウェアと、ＣＰＵやプロセッサの動作を規定するソフトウェア（プログラム）との協働により実現することができる。

図１の超音波診断装置の全体構成は以上のとおりである。次に、図１の超音波診断装置により実現される機能の具体例について詳述する。なお、図１に示した構成（符号を付された各部）については、以下の説明において図１の符号を利用する。

図３は、移動ベクトルを導出する処理の具体例を示す図である。移動ベクトル演算部３０は、複数フレームの三次元超音波画像（３Ｄ超音波画像）の互いに隣接する２つのフレームの間において、三次元超音波画像を構成する複数画素の各々についての移動ベクトルを導出（算出）する。移動ベクトルの導出にはパターンマッチング処理が利用される。

図３には、複数フレームのうちの互いに隣接するフレームｎ（ｎは自然数）とフレーム（ｎ＋１）の３Ｄ超音波画像が図示されている。移動ベクトル演算部３０は、フレームｎの３Ｄ超音波画像を構成する複数画素の各画素ごとに、パターンマッチング処理により移動ベクトルを算出し、フレームｎの３Ｄ超音波画像を構成する複数画素、望ましくは全画素について、各画素の移動ベクトルを導出する。

パターンマッチング処理においては、フレームｎの３Ｄ超音波画像の各画素ごとに、例えばその画素を取り囲むようにテンプレートＴが設定される。そして、フレーム（ｎ＋１）の３Ｄ超音波画像内において、例えばテンプレートＴに対応した位置の画像領域を含むように探索領域ＳＡが設定される。探索領域ＳＡの設定には、公知の様々な手法を利用することができる。もちろん、例えばフレーム（ｎ＋１）の３Ｄ超音波画像内の全域が探索領域ＳＡとされてもよい。

テンプレートＴと探索領域ＳＡが設定されると、探索領域ＳＡ内においてテンプレートＴが移動され、各移動位置において、フレームｎの３Ｄ超音波画像のテンプレートＴ内の複数画素と、フレーム（ｎ＋１）の３Ｄ超音波画像のテンプレートＴ内の複数画素とに基づいて、相関演算が行われる。例えば、フレームｎ内のテンプレートＴに対応した位置を初期位置とし、その初期位置からの各変位(dx,dy)ごとに相関値が算出され、探索領域ＳＡ内の全域に亘る複数変位に対応した相関値が算出される。

相関値とは画像データ間（複数画素と複数画素）の相関関係の程度（類似の程度）を示す数値であり、相関値の算出には相関演算の各手法に応じた公知の数式が用いられる。相関値としては、例えば、ＳＳＤ（Sum of Square Difference：差の二乗和）などが好適であり、ＳＳＤは、類似の度合が大きいほど小さな値を示す。なお、例えば、位相限定相関法や相互相関法等により類似の度合が大きいほど大きな値を示す相関値が利用されてもよい。

探索領域ＳＡ内の全域に亘る複数変位に対応した相関値が算出されると、複数変位の中から最も類似の度合が大きい変位が特定され、演算対象となる各画素の移動先が決定される。

図３には、フレームｎの３Ｄ超音波画像を構成する複数画素の代表例として画素Ａが図示されている。そして、画素Ａに関するパターンマッチング処理により、フレーム（ｎ＋１）の３Ｄ超音波画像内において、画素Ａの移動先として画素Ｂが特定される。移動ベクトル演算部３０は、画素Ａの位置（座標）を始点として画素Ｂの位置（座標）を終点とするベクトルを画素Ａの移動ベクトルＡＢとする。移動ベクトル演算部３０は、フレームｎ内の複数画素について、望ましくは全画素について、各画素ごとに、その画素の位置を始点としてその画素の移動先を終点とする移動ベクトルを導出する。

さらに、移動ベクトル演算部３０は、フレームｎの３Ｄ超音波画像を構成する各画素の移動ベクトルに対して、平滑化処理を施すことが望ましい。

図４は、移動ベクトルに対する平滑化処理の具体例を示す図である。移動ベクトル演算部３０は、各注目画素ごとに、その注目画素とその近傍に位置する複数画素に対応した複数の移動ベクトルの平均ベクトルを算出し、算出した平均ベクトルをその注目画素における平滑化後の移動ベクトルとする。

図４には、複数画素の代表例として画素Ａが図示されており、画素Ａとその周辺の８画素からなる９画素に対応した９つの移動ベクトルの平均ベクトルが算出され、算出された平均ベクトルが画素Ａの平滑化後の移動ベクトルとされる。なお、移動ベクトル演算部３０は、平均ベクトルの算出において、複数の移動ベクトルの中から特異な移動ベクトルを除外してから平均ベクトルを算出するようにしてもよい。

各画素の移動ベクトルを平滑化処理することにより、例えば、ある画素の移動ベクトルが誤検出された場合においても、その誤検出に伴う悪影響を低減または除去することができる。

移動ベクトル演算部３０により、複数画素に対応した複数の移動ベクトル（平滑化後の移動ベクトルが望ましい）が算出されると、算出された移動ベクトルに基づいて三次元中間画像（３Ｄ中間画像）が形成される。

図５は、三次元中間画像を形成する処理の具体例１を示す図である。中間画像形成部５０は、複数フレームの三次元超音波画像（３Ｄ超音波画像）の互いに隣接する２つのフレームの間において得られた複数画素の移動ベクトルに基づいて、そのフレーム間に追加される１又は複数フレームの三次元中間画像（３Ｄ中間画像）を形成する。

図５には、３Ｄ超音波画像のフレームｎ（ｎは自然数）とフレーム（ｎ＋１）の間に追加される中間画像を形成する具体例が図示されている。図５（１）に示す移動ベクトルＡＢはフレームｎの３Ｄ超音波画像内における画素Ａに対応した移動ベクトル（図３，図４参照）である。つまり、画素Ａの移動先がフレーム（ｎ＋１）の３Ｄ超音波画像内における画素Ｂの位置（座標）である。中間画像形成部５０は、図５（１）に示す移動ベクトルＡＢから図５（２）に示す中間ベクトルＡＣを得る。

図６は、中間ベクトルの具体例１を示す図である。図６には、フレームｎの３Ｄ超音波画像とフレーム（ｎ＋１）の３Ｄ超音波画像の間における中間ベクトルの具体例が図示されている。移動ベクトルＡＢは、フレームｎの３Ｄ超音波画像内における画素Ａに対応した移動ベクトルであり、画素Ａの移動先がフレーム（ｎ＋１）の３Ｄ超音波画像内における画素Ｂの位置（座標）である。

また、図６には、フレームレートをＭ倍（Ｍは２以上の整数）とする具体例が図示されている。つまり、フレームｎの３Ｄ超音波画像とフレーム（ｎ＋１）の３Ｄ超音波画像の間に、（Ｍ−１）フレームの３Ｄ中間画像が追加される。図６に示す具体例では、フレームｎの３Ｄ超音波画像に近い（時相的に近い）順に、フレーム１，２，３，・・・，Ｍ−１の３Ｄ中間画像が配置されている。

中間画像形成部５０は、画素Ａに対応した移動ベクトルＡＢから画素Ａに関する中間ベクトルＡＣ（ＡＣ１，ＡＣ２，ＡＣ３，・・・，ＡＣ（Ｍ−１））を導出する。中間ベクトルＡＣの始点と方向は移動ベクトルＡＢと同じであり、中間ベクトルＡＣの大きさ（ベクトルの長さ）は、３Ｄ中間画像のフレーム番号（１，２，３，・・・，Ｍ−１）に応じて決定される。

例えば、図６において、フレーム１の３Ｄ中間画像に対応した中間ベクトルＡＣ１は、始点と方向は移動ベクトルＡＢと同じであり、大きさ（ベクトルの長さ）が移動ベクトルＡＢの１／Ｍとされる。また。フレーム２の３Ｄ中間画像に対応した中間ベクトルＡＣ２は、始点と方向は移動ベクトルＡＢと同じであり、大きさ（ベクトルの長さ）が移動ベクトルＡＢの２／Ｍとされる。そして、フレーム（Ｍ−１）の３Ｄ中間画像に対応した中間ベクトルＡＣ（Ｍ−１）は、始点と方向は移動ベクトルＡＢと同じであり、大きさ（ベクトルの長さ）が移動ベクトルＡＢの（Ｍ−１）／Ｍとされる。

中間画像形成部５０は、フレームｎの３Ｄ超音波画像を構成する複数画素の各画素ごとに、その画素に対応した移動ベクトルに基づいて、３Ｄ中間画像の各フレームに対応した中間ベクトルを生成する。

図５に戻り、図５（２）に示す中間ベクトルＡＣは、フレームｎの３Ｄ超音波画像の画素Ａに関する中間ベクトルであり、３Ｄ中間画像の注目フレーム（例えば図６におけるフレーム１〜Ｍ−１までの３Ｄ中間画像のいずれか）に対応した中間ベクトルである。

中間画像形成部５０は、フレームｎの３Ｄ超音波画像を構成する複数画素について、望ましくは全画素について、各画素ごとに中間ベクトルを生成する。そして、中間画像形成部５０は、複数画素に対応した複数の中間ベクトルのうち、ベクトルの終点が注目画素Ｚの位置（座標）に最も近い中間ベクトルを選択する。図５（２）に示す具体例では、中間ベクトルＡＣが選択される。

さらに、中間画像形成部５０は、選択した中間ベクトルの終点が注目画素の位置となるように中間ベクトルを平行移動する。図５の具体例では、図５（２）の中間ベクトルＡＣを平行移動した結果として、図５（３）の中間ベクトルＡ´Ｃ´が得られる。

そして、中間画像形成部５０は、平行移動後の中間ベクトルの始点近傍における複数画素の画素値に基づいて注目画素の画素値を決定する。例えば、図５（４）に示すように、フレームｎの３Ｄ超音波画像内における座標Ａ´の近傍４点の画素の画素値から、例えば線形補間処理等により、注目画素Ｚの画素値が算出される。そして、算出された注目画素Ｚの画素値が、３Ｄ中間画像の注目フレームを構成する中間画素とされる。

中間画像形成部５０は、３Ｄ中間画像の注目フレームを構成する全画素の各々を注目画素として、図５を利用して説明した処理を実行し、３Ｄ中間画像の注目フレームを構成する全画素の画素値を得ることにより、その注目フレームに対応した３Ｄ中間画像を形成する。さらに、中間画像形成部５０は、複数フレームの３Ｄ中間画像（例えば図６におけるフレーム１〜Ｍ−１までの３Ｄ中間画像）の各フレームを注目フレームとして、図５を利用して説明した処理を実行することにより、複数フレームの３Ｄ中間画像を形成する。

図７は、三次元中間画像を形成する処理の具体例２を示す図である。図７には、３Ｄ超音波画像のフレームｎ（ｎは自然数）とフレーム（ｎ＋１）の間に追加される中間画像を形成する具体例が図示されている。図７（１）に示す移動ベクトルＡＢはフレームｎの３Ｄ超音波画像内における画素Ａに対応した移動ベクトル（図３，図４参照）である。つまり、画素Ａの移動先がフレーム（ｎ＋１）の３Ｄ超音波画像内における画素Ｂの位置（座標）である。中間画像形成部５０は、図７（１）に示す移動ベクトルＡＢから図７（２）に示す中間ベクトルＣＢを得る。

図８は、中間ベクトルの具体例２を示す図である。図８には、フレームｎの３Ｄ超音波画像とフレーム（ｎ＋１）の３Ｄ超音波画像の間における中間ベクトルの具体例が図示されている。移動ベクトルＡＢは、フレームｎの３Ｄ超音波画像内における画素Ａに対応した移動ベクトルであり、画素Ａの移動先がフレーム（ｎ＋１）の３Ｄ超音波画像内における画素Ｂの位置（座標）である。

また、図８には、フレームレートをＭ倍（Ｍは２以上の整数）とする具体例が図示されている。つまり、フレームｎの３Ｄ超音波画像とフレーム（ｎ＋１）の３Ｄ超音波画像の間に、（Ｍ−１）フレームの３Ｄ中間画像が追加される。図８に示す具体例では、フレームｎの３Ｄ超音波画像に近い（時相的に近い）順に、フレーム１，２，３，・・・，Ｍ−１の３Ｄ中間画像が配置されている。

中間画像形成部５０は、画素Ａに対応した移動ベクトルＡＢから、移動ベクトルＡＢの終点である画素Ｂを基準として、中間ベクトルＣＢ（Ｃ１Ｂ，Ｃ２Ｂ，Ｃ３Ｂ，・・・，Ｃ（Ｍ−１）Ｂ）を導出する。中間ベクトルＣＢの終点と方向は移動ベクトルＡＢと同じであり、中間ベクトルＣＢの大きさ（ベクトルの長さ）は、３Ｄ中間画像のフレーム番号（１，２，３，・・・，Ｍ−１）に応じて決定される。

例えば、図８において、フレーム１の３Ｄ中間画像に対応した中間ベクトルＣ１Ｂは、終点と方向は移動ベクトルＡＢと同じであり、大きさ（ベクトルの長さ）が移動ベクトルＡＢの１／Ｍとされる。また。フレーム２の３Ｄ中間画像に対応した中間ベクトルＣ２Ｂは、終点と方向は移動ベクトルＡＢと同じであり、大きさ（ベクトルの長さ）が移動ベクトルＡＢの２／Ｍとされる。そして、フレーム（Ｍ−１）の３Ｄ中間画像に対応した中間ベクトルＣ（Ｍ−１）Ｂは、終点と方向は移動ベクトルＡＢと同じであり、大きさ（ベクトルの長さ）が移動ベクトルＡＢの（Ｍ−１）／Ｍとされる。

中間画像形成部５０は、フレームｎの３Ｄ超音波画像を構成する複数画素の各画素ごとに、その画素に対応した移動ベクトルに基づいて、３Ｄ中間画像の各フレームに対応した中間ベクトルを生成する。

図７に戻り、図７（２）に示す中間ベクトルＣＢは、フレームｎの３Ｄ超音波画像の画素Ａに関する中間ベクトルであり、３Ｄ中間画像の注目フレーム（例えば図８におけるフレーム１〜Ｍ−１までの３Ｄ中間画像のいずれか）に対応した中間ベクトルである。

中間画像形成部５０は、フレームｎの３Ｄ超音波画像を構成する複数画素について、望ましくは全画素について、各画素ごとに中間ベクトルを生成する。そして、中間画像形成部５０は、複数画素に対応した複数の中間ベクトルのうち、ベクトルの始点が注目画素Ｚの位置（座標）に最も近い中間ベクトルを選択する。図７（２）に示す具体例では、中間ベクトルＣＢが選択される。

さらに、中間画像形成部５０は、選択した中間ベクトルの始点が注目画素の位置となるように中間ベクトルを平行移動する。図７の具体例では、図７（２）の中間ベクトルＣＢを平行移動した結果として、図７（３）の中間ベクトルＣ´Ｂ´が得られる。

そして、中間画像形成部５０は、平行移動後の中間ベクトルの終点近傍における複数画素の画素値に基づいて注目画素の画素値を決定する。例えば図７（４）に示すように、フレーム（ｎ＋１）の３Ｄ超音波画像内における座標Ｂ´の近傍４点の画素の画素値から、例えば線形補間処理等により、注目画素Ｚの画素値が算出される。そして、算出された注目画素Ｚの画素値が、３Ｄ中間画像の注目フレームを構成する中間画素とされる。

中間画像形成部５０は、３Ｄ中間画像の注目フレームを構成する全画素の各々を注目画素として、図７を利用して説明した処理を実行し、３Ｄ中間画像の注目フレームを構成する全画素の画素値を得ることにより、その注目フレームに対応した３Ｄ中間画像を形成する。さらに、中間画像形成部５０は、複数フレームの３Ｄ中間画像（例えば図８におけるフレーム１〜Ｍ−１までの３Ｄ中間画像）の各フレームを注目フレームとして、図７を利用して説明した処理を実行することにより、複数フレームの３Ｄ中間画像を形成する。

例えば、図５と図６を利用して説明した具体例１または図７と図８を利用して説明した具体例２により、複数フレームの３Ｄ超音波画像の複数フレーム間に追加される複数フレームの３Ｄ中間画像が形成される。

上述したように、三次元中間画像（３Ｄ中間画像）の形成には、移動ベクトル演算部３０において導出される移動ベクトルが利用される。移動ベクトル演算部３０は、三次元超音波画像（３Ｄ超音波画像）のフレームｎ内の複数画素について、望ましくは全画素について、各画素ごとに、その画素の位置を始点としてその画素の移動先を終点とする移動ベクトルを導出する（図３参照）。

移動ベクトル評価部４０は、３Ｄ超音波画像の各フレーム内の複数画素に対応した複数の移動ベクトルに関する信頼性を評価する。例えば、移動ベクトル演算部３０によるパターンマッチング処理において、最も相関が強い（最も類似の度合が大きい）として選択された各画素の移動先における相関値に基づいて移動ベクトルの信頼性が評価される。

具体的には、相関が強い（類似の程度が大きい）ほど小さな値となる相関値が利用されており、算出された相関値が閾値を超えた場合（または閾値以上の場合）には、その相関値が信頼できる値ではなく、その相関値に基づいて得られた移動ベクトルの信頼性が低いと判定される。そして、例えば、３Ｄ超音波画像の各フレームごとに、そのフレームを構成する複数画素のうちの、移動ベクトルの信頼性が低いと判定された画素数の割合が、そのフレームにおける移動ベクトルの信頼性の評価値とされる。この具体例の場合、各フレームごとに得られる信頼性の評価値が大きいほど、そのフレームの移動ベクトルの信頼性が低いことを示している。

そして、例えば、各フレームの移動ベクトルの信頼性が低い場合には、例えば、信頼性の評価値が閾値よりも大きい場合には、例えば、パターンマッチングにおける探索領域ＳＡ（図３）が広げられ、各画素の移動先の選択範囲を広げて移動先の特定精度を高めることが望ましい。なお、探索領域ＳＡを広げるとパターンマッチングの処理時間が増大するため、例えば、その処理時間の増大に応じて、３Ｄ超音波画像の各フレーム間に追加される３Ｄ中間画像のフレーム数（フレームの枚数）を少なくするようにしてもよい。また、各フレームの移動ベクトルの信頼性が低い場合に、テンプレートＴ（図３）を大きくすることにより、各画素の移動先の特定精度を高めるようにしてもよい。

また、移動ベクトル評価部４０において得られた移動ベクトルの信頼性の評価値に基づいて、制御部１００が超音波の送受信を制御するようにしてもよい。例えば、制御部１００は、移動ベクトル評価部４０から得られる移動ベクトルの信頼性の評価値に基づいて、ボリュームレートを制御してもよい。

図９は、ボリュームレート制御の具体例を示す図である。図９（１）には、複数時相に亘って各時相ごとに得られるボリュームデータの具体例が図示されている。つまり、時相１〜時相６に対した複数のボリュームデータが代表的に図示されている。なお、図示省略した時相７以降においても各時相ごとにボリュームデータが得られてもよい。

図９（２）には、複数フレームの三次元超音波画像（３Ｄ超音波画像）の具体例が図示されている。各フレームの３Ｄ超音波画像は、そのフレームに対応した時相のボリュームデータに基づいて形成される。図９に示す具体例において、時相１のボリュームデータに基づいてフレーム１の３Ｄ超音波画像が形成され、時相２のボリュームデータに基づいてフレーム２の３Ｄ超音波画像が形成される。つまり、ボリュームデータの時相番号と３Ｄ超音波画像のフレーム番号が互いに対応している。

図９（３）には、三次元中間画像（３Ｄ中間画像）の具体例が図示されており、図９（４）には、三次元動画像（３Ｄ動画像）の具体例が図示されている。図９に示す具体例では、例えば、３Ｄ超音波画像のフレーム１とフレーム２の間に３フレーム（フレーム数３）の３Ｄ中間画像が追加され、３Ｄ超音波画像のフレーム２とフレーム３の間にも３フレーム（フレーム数３）の３Ｄ中間画像が追加され、さらに、３Ｄ超音波画像のフレーム３以降にも各フレーム間に３Ｄ中間画像が追加されて、３Ｄ動画像のフレーム列が構成されている。

図９に示す具体例では、３Ｄ超音波画像のフレーム１とフレーム２の間において得られた移動ベクトルの信頼性が低い。移動ベクトルの信頼性が低い場合に、制御部１００は、ボリュームレートを増加させる。但し、図９に示す具体例では、時相２のボリュームデータが得られた直後から、標準レート（標準のボリュームレート）で時相３のボリュームデータの取得が開始されているため、その標準レートを維持したまま時相３のボリュームデータが取得される。そして、時相３のボリュームデータが取得された直後からボリュームレートが増大され、高レート（標準よりも高いボリュームレート）で時相４と時相５のボリュームデータが取得される。なお、時相２のボリュームデータが得られた直後から高レートに切り換えて時相３のボリュームデータを得るようにしてもよい。

また、図９に示す具体例では、３Ｄ超音波画像のフレーム３とフレーム４の間において得られた移動ベクトルの信頼性が高い。移動ベクトルの信頼性が高い場合に、制御部１００は、ボリュームレートを減少させてもよい。但し、図９に示す具体例では、時相４のボリュームデータが得られた直後から、高レートで時相５のボリュームデータの取得が開始されているため、その高レートを維持したまま時相５のボリュームデータが取得される。そして、時相５のボリュームデータが取得された直後からボリュームレートが減少され、つまり標準レートに戻され、標準レートで時相６のボリュームデータが取得される。なお時相４のボリュームデータが得られた直後から標準レートに切り換えて時相５のボリュームデータを得るようにしてもよい。

制御部１００は、超音波の送受に係る複数の送受信パラメータに基づいてボリュームレートを制御する。複数の送受信パラメータには、例えば、空間フレーム数（空間フレーム密度）、ビームライン数（ライン密度）、ＰＲＴ（pulse repetition time）などが含まれる。

空間フレーム数は、超音波が立体的に送受されるボリューム（三次元領域）を構成する空間的なフレーム（空間フレーム）の枚数（フレーム数）である。例えば、ボリュームの大きさが一定で空間フレーム数が変化すると空間フレーム密度も変化する。ビームライン数は、各空間フレームを構成する受信ビームラインの本数である。例えば、各空間フレームの大きさが一定でビームライン数が変化するとライン密度も変化する。そして、ＰＲＴは、超音波のパルス繰り返し周期である。

これら複数の送受信パラメータによりボリュームレートが決定される。標準レート（標準のボリュームレート）に対応した複数の送受信パラメータは、例えば、ボリュームデータのデータ密度を重視して設定される。一方、高レートに対応した複数の送受信パラメータは、標準レート（標準のボリュームレート）よりもボリュームレートが高くなるように設定される。

例えば、高レートの空間フレーム数は標準レートの空間フレーム数よりも小さく（少なく）設定され、高レートのビームライン数は標準レートのビームライン数よりも小さく（少なく）設定され、高レートのＰＲＴは標準レートのＰＲＴよりも小さく（短く）設定される。

制御部１００は、標準レートに対応した複数の送受信パラメータで送受信部１２を制御することにより、標準レートによるボリュームデータの取得を実現し、高レートに対応した複数の送受信パラメータで送受信部１２を制御することにより、高レートによるボリュームデータの取得を実現する。

なお、送受信パラメータとして、パラレル受信（パラレル受信の本数（１本以上））やＴＨＩ（tissue harmonic imagingを行うか否か）が変更されて、ボリュームレートが変更されてもよい。

図９に示す具体例では、ボリュームレートの変更に係わらず、３Ｄ動画像のフレームレートを一定に維持することができる。

例えば、高レート（高いボリュームレート）が標準レート（標準のボリュームレート）の２倍である場合に、標準レートに対応した３Ｄ超音波画像の各フレーム間に３フレームの３Ｄ中間画像を追加し、高レートに対応した３Ｄ超音波画像の各フレーム間に１フレームの３Ｄ中間画像を追加することにより、３Ｄ動画像のフレーム間隔（時相間隔）が一定に維持され、３Ｄ動画像のフレームレートを一定に維持することができる。なお、図９に示す具体例はあくまでも一例であり、３Ｄ動画像のフレームレートを必ずしも一定にする必要はない。

また、３Ｄ動画像のフレームレートは、表示部８２において対応可能な表示フレームレートに応じて決定されることが望ましい。例えば、表示部８２の表示フレームレートが６０Ｈｚ（ヘルツ）であるならば、３Ｄ動画像のフレームレートが６０Ｈｚとなるように、ボリュームレートと３Ｄ中間画像のフレーム数（３Ｄ超音波画像の各フレーム間に追加される３Ｄ中間画像のフレーム数）を決定すればよい。

さらに、表示画像形成部８０は、３Ｄ動画像をスロー再生するようにしてもよい。そして、スロー再生時には、再生速度に応じて３Ｄ中間画像のフレーム数（３Ｄ超音波画像の各フレーム間に追加される３Ｄ中間画像のフレーム数）を決定することが望ましい。

例えば、通常再生（１倍速再生）時における各フレーム間の３Ｄ中間画像の枚数（フレーム数）がＫ枚（フレーム数Ｋ：但しＫは自然数）の場合に、通常再生の１／２倍速のスロー再生時における各フレーム間の３Ｄ中間画像の枚数（フレーム数）を２Ｋとし、通常再生の１／３倍速のスロー再生時における各フレーム間の３Ｄ中間画像の枚数（フレーム数）を３Ｋとし、通常再生の１／４倍速のスロー再生時における各フレーム間の３Ｄ中間画像の枚数（フレーム数）を４Ｋとする。つまり、通常再生の１／Ｓ倍速（Ｓは自然数）のスロー再生時における各フレーム間の３Ｄ中間画像の枚数（フレーム数）をＳ×Ｋとする。

一般に従来のスロー再生では、通常再生で利用される複数の画像フレームのうちの各画像フレームが複数の表示フレームに対応付けられて繰り返し再生される。つまり、同一の画像フレームが複数の表示フレームで利用される。

これに対し、スロー再生時に再生速度に応じて３Ｄ中間画像のフレーム数を決定する上記処理により、同一の画像フレーム（３Ｄ動画像の各フレーム）が利用される表示フレーム数を少なくすることができ、望ましくは、各画像フレームと各表示フレームとを１対１に対応付けたスロー再生を実現することができる。これにより、例えば、画像内容が滑らかに変化するスロー再生が可能になる。

以上、本発明の好適な実施形態を説明したが、上述した実施形態は、あらゆる点で単なる例示にすぎず、本発明の範囲を限定するものではない。本発明は、その本質を逸脱しない範囲で各種の変形形態を包含する。

１０ プローブ、１２ 送受信部、１４ ビームデータ処理部、１６ ボリュームデータ記憶部、２０ 三次元画像形成部、３０ 移動ベクトル演算部、４０ 移動ベクトル評価部、５０ 中間画像形成部、８０ 表示画像形成部、８２ 表示部、９０ 操作デバイス、１００ 制御部。

Claims (6)

1. 超音波を立体的に送受して得られた複数時相のボリュームデータに基づいて、当該複数時相に対応した複数フレームの三次元超音波画像を形成する超音波画像形成部と、
   前記複数フレームの三次元超音波画像の各フレーム間における画像の移動情報に基づいて、当該各フレーム間に追加される１又は複数フレームの三次元中間画像を形成する中間画像形成部と、
   前記複数フレームの三次元超音波画像とそれらの複数フレーム間に追加された前記複数フレームの三次元中間画像に基づいて三次元動画像を表示する表示処理部と、
   前記移動情報の信頼性の程度に応じて、単位時間あたりのボリュームデータ数であるボリュームレートを制御する制御部と、
   を有する、
   ことを特徴とする超音波診断装置。
2. 請求項１に記載の超音波診断装置において、
   前記制御部は、信頼性の判定条件に基づいて前記移動情報の信頼性が低いと判定された場合に、前記ボリュームレートを増加させて前記三次元超音波画像の単位時間あたりのフレーム数であるフレームレートを増加させる、
   ことを特徴とする超音波診断装置。
3. 請求項１または２に記載の超音波診断装置において、
   前記制御部は、超音波を立体的に送受するボリューム内における空間フレーム数と、各空間フレームを構成するビームライン数と、超音波のパルス繰り返し周期と、を含む複数の送受信パラメータのうちの少なくとも１つを変更することにより、前記ボリュームレートを制御する、
   ことを特徴とする超音波診断装置。
4. 請求項１から３のいずれか１項に記載の超音波診断装置において、
   前記中間画像形成部は、前記移動情報として、前記各フレームの三次元超音波画像を構成する複数の画像要素の各画像要素ごとに各フレーム間における相関演算に基づいて移動ベクトルを導出し、複数の画像要素とそれらの移動ベクトルに基づいて当該各フレーム間に対応した複数の中間画像要素を得ることにより、複数の中間画像要素で構成された前記三次元中間画像を形成する、
   ことを特徴とする超音波診断装置。
5. 請求項４に記載の超音波診断装置において、
   前記複数の画像要素の各画像要素ごとに前記相関演算により得られる相関値の大きさに基づいて、前記移動情報として導出される移動ベクトルの信頼性の程度を示す評価値を算出する評価部をさらに有し、
   前記制御部は、移動ベクトルの信頼性を示す前記評価値に基づいて、単位時間あたりのボリュームデータ数であるボリュームレートを制御する、
   ことを特徴とする超音波診断装置。
6. 超音波を立体的に送受して得られた複数時相のボリュームデータに基づいて、当該複数時相に対応した複数フレームの三次元超音波画像を形成する超音波画像形成部と、
   前記複数フレームの三次元超音波画像の各フレーム間における画像の移動情報に基づいて、当該各フレーム間に追加される１又は複数フレームの三次元中間画像を形成する中間画像形成部であって、前記移動情報として、前記各フレームの三次元超音波画像を構成する複数の画像要素の各画像要素ごとに各フレーム間における相関演算に基づいて移動ベクトルを導出し、複数の画像要素とそれらの移動ベクトルに基づいて当該各フレーム間に対応した複数の中間画像要素を得ることにより、複数の中間画像要素で構成された前記三次元中間画像を形成する中間画像形成部と、
   前記複数の画像要素の各画像要素ごとに前記相関演算により得られる相関値の大きさに基づいて、前記移動情報として導出される移動ベクトルの信頼性の程度を示す評価値を算出する評価部と、
   移動ベクトルの信頼性を示す前記評価値に基づいて、単位時間あたりのボリュームデータ数であるボリュームレートを制御する制御部と、
   前記複数フレームの三次元超音波画像とそれらの複数フレーム間に追加された前記複数フレームの三次元中間画像に基づいて三次元動画像を表示する表示処理部と、
   を有する、
   ことを特徴とする超音波診断装置。

Images