(57)【要約】【課題】拍動する心筋組織の各心拍時相における組織画像データや心機能パラメータデータ等を正確に収集する。【解決手段】輪郭面検出部６は、被検体の心臓領域に対する３次元走査によって収集された時系列的なボリュームデータをトラッキング処理することにより各心拍時相における心筋の心内膜輪郭面及び心外膜輪郭面を検出する。そして、診断用画像データ生成部１０は、心機能パラメータ計測部９が前記輪郭面に配置された複数の追跡点における心拍時相間の変位を計測して収集した心機能パラメータデータと、ワイヤフレーム形成部７が前記輪郭面に対して形成したワイヤフレームのメッシュパターンと、組織画像データ生成部８が前記輪郭面における前記ボリュームデータのボクセルを切り出して生成した組織画像データとを重畳して診断用画像データを生成し、得られた各心拍時相の診断用画像データを表示部１１において動画像表示する。【選択図】図１
2010115372.tif 000002
【特許請求の範囲】
【請求項１】
被検体の心臓領域に対する３次元走査によって得られた受信信号に基づいて複数の心拍時相におけるボリュームデータを時系列的に生成するボリュームデータ生成手段と、
基準心拍時相のボリュームデータにおける心筋組織の輪郭面を検出し、前記基準心拍時相のボリュームデータと他の心拍時相のボリュームデータとのトラッキング処理により前記他の心拍時相のボリュームデータにおける心筋組織の輪郭面を検出する輪郭面検出手段と、
前記輪郭面における前記ボリュームデータのボクセルあるいは前記輪郭面を含む閉曲面によって囲まれた前記ボリュームデータのボクセルを切り出して各心拍時相における組織画像データを生成する組織画像データ生成手段と、
前記組織画像データを処理して診断用画像データを生成する診断用画像データ生成手段と
前記診断用画像データを表示する表示手段とを
備えたことを特徴とする３次元超音波診断装置。
【請求項２】
前記輪郭面の「変位」に基づいて心機能パラメータを計測し、この計測結果を用いて心機能パラメータデータを生成する心機能パラメータ計測手段を備え、前記診断用画像データ生成手段は、前記組織画像データに前記心機能パラメータデータを重畳して前記診断用画像データを生成することを特徴とする請求項１記載の３次元超音波診断装置。
【請求項３】
前記輪郭面に対してワイヤフレームを形成するワイヤフレーム形成手段を備え、前記診断用画像データ生成手段は、前記ワイヤフレーム形成手段によって形成されたワイヤフレームのメッシュパターンを前記組織画像データに重畳して前記診断用画像データを生成することを特徴とする請求項１又は請求項２に記載した３次元超音波診断装置。
【請求項４】
前記基準心拍時相のボリュームデータを用いて生成された所定断面のＭＰＲ画像データにおける心筋組織に対し輪郭点あるいは輪郭線を設定する輪郭設定手段を備え、前記輪郭面検出手段は、心筋組織に設定された複数からなる前記輪郭点あるいは前記輪郭線を補間処理して前記基準心拍時相のボリュームデータにおける前記心筋組織の輪郭面を検出することを特徴とする請求項１記載の３次元超音波診断装置。
【請求項５】
前記輪郭設定手段は、心筋組織の心内膜及び心外膜の少なくとも何れかに対し前記輪郭線あるいは前記輪郭点を設定することを特徴とする請求項４記載の３次元超音波診断装置。
【請求項６】
前記輪郭設定手段は、前記輪郭面検出手段が検出した前記心筋組織の輪郭面を所望の位置に更新することを特徴とする請求項４記載の３次元超音波診断装置。
【請求項７】
前記輪郭面検出手段は、前記基準心拍時相のボリュームデータにおける前記心筋組織の輪郭面に複数の追跡点を配置し、前記追跡点の各々を中心とした３次元関心領域における前記基準心拍時相のボリュームデータ及び他の心拍時相のボリュームデータに対するトラッキング処理により前記他の心拍時相のボリュームデータにおける心筋組織の輪郭面を検出することを特徴とする請求項１記載の３次元超音波診断装置。
【請求項８】
前記組織画像データ生成手段は、前記輪郭面から切り出した前記ボリュームデータのボクセルに基づいてＭＰＲ（Multi Planar Reconstruction）画像データを生成することを特徴とする請求項１記載の３次元超音波診断装置。
【請求項９】
前記組織画像データ生成手段は、前記閉曲面によって切り出した前記ボリュームデータのボクセルを処理して３次元画像データあるいは投影画像データを生成することを特徴とする請求項１記載の３次元超音波診断装置。
【請求項１０】
前記組織画像データ生成手段は、前記閉曲面によって切り出した前記ボリュームデータのボクセルが有するボクセル値及び前記ボクセルから前記閉曲面までの距離に基づいて前記ボクセルの不透明度を設定し、前記不透明度が設定された前記ボクセルをレンダリング処理して前記３次元画像データを生成することを特徴とする請求項９記載の３次元超音波診断装置。
【請求項１１】
前記心機能パラメータ計測手段は、前記心機能パラメータの計測結果に対し所定の色相と不透明度を設定して前記心機能パラメータデータを生成することを特徴とする請求項２記載の３次元超音波診断装置。
【請求項１２】
前記心機能パラメータ計測手段は、前記輪郭面における心拍時相間の「変位」に基づいて、「歪み（ストレイン）」、「回転」、「捻れ（torsion）」、「速度」、「歪みレート」、「回転レート」、「捻れレート」及び「加速度」の少なくとも何れかを前記心機能パラメータとして計測することを特徴とする請求項２記載の３次元超音波診断装置。
【請求項１３】
前記ワイヤフレーム形成手段は、前記輪郭面に配置された前記複数の追跡点を結線することによって前記ワイヤフレームを形成することを特徴とする請求項７記載の３次元超音波診断装置。
【請求項１４】
前記組織画像データ生成手段は、前記追跡点によって囲まれた前記ワイヤフレームの微小平面に対応する前記ボリュームデータのボクセルを切り出して微小ＭＰＲ画像データを生成し、この微小ＭＰＲ画像データを前記微小平面に貼付してＭＰＲ（Multi Planar Reconstruction）画像データを生成することを特徴とする請求項１３記載の３次元超音波診断装置。
【請求項１５】
前記組織画像データ生成手段は、前記ＭＰＲ画像データに対しスペックルの強調を目的としたフィルタリング処理を行なうことを特徴とする請求項８又は請求項１４に記載した３次元超音波診断装置。
【請求項１６】
前記表示手段は、前記診断用画像データ生成手段が生成した前記診断用画像データを心拍時相に対応させて繰り返し表示することを特徴とする請求項１記載の３次元超音波診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、３次元超音波診断装置に係り、特に、被検体の心臓領域から収集される時系列的なボリュームデータに基づいて組織画像データの生成や心機能パラメータの計測を行なう３次元超音波診断装置に関する。
【背景技術】
【０００２】
超音波診断装置は、超音波プローブに設けられた振動素子から発生する超音波を被検体内に放射し、被検体組織の音響インピーダンスの差異によって生ずる反射波を前記振動素子により受信して生体情報を収集するものであり、超音波プローブを体表に接触させるだけの簡単な操作で超音波画像データのリアルタイム表示が可能となるため、各種臓器の形態診断や機能診断に広く用いられている。
【０００３】
特に、複数の振動素子が一次元配列された超音波プローブを機械的に移動させる方法や複数の振動素子が二次元配列された超音波プローブを用いる方法により被検体の三次元データ（ボリュームデータ）を収集する３次元超音波診断装置が開発され、このボリュームデータに基づいて得られる３次元画像データ、ＭＰＲ（Multi Planar Reconstruction）画像データ、投影画像データ等（以下では、これらを纏めて組織画像データと呼ぶ。）の観察により更に高度な診断や治療が可能となっている。
【０００４】
又、近年の心臓領域における超音波診断では、ボリュームデータのスペックルをトラッキング処理することにより拍動する心筋組織の時間的変化や「歪み（ストレイン）」等によって代表される心機能パラメータを３次元的に捉えることが可能となった。
【０００５】
心機能パラメータの計測では、被検体に対する３次元走査により得られた受信信号に基づいてボリュームデータを時系列的に収集し、基準心拍時相（例えば拡張末期）のボリュームデータにＡＣＴ法等を適用して心筋組織の心内膜輪郭面や心外膜輪郭面を検出する。そして、基準心拍時相のボリュームデータと他の心拍時相のボリュームデータとの３次元トラッキング処理により、心内膜輪郭面や心外膜輪郭面の各部における「変位」等を心機能パラメータとして計測する。
【０００６】
又、上述の方法によって計測された心機能パラメータの分布を表示する際、先ず、各心拍時相の心内膜輪郭面や心外膜輪郭面の位置を示すワイヤフレームを形成し、次いで、心機能パラメータの計測結果に所定の色相を設定して生成した心機能パラメータデータを前記ワイヤフレームに重畳させて表示する方法が一般に行われている（例えば、特許文献１参照。）。この場合、心内膜輪郭面や心外膜輪郭面の運動状態は、ワイヤフレームを示す３次元的なメッシュパターンの動きによって把握することができる。即ち、ワイヤフレームのメッシュパターンに心機能パラメータデータを重畳して表示することにより、心筋組織の運動状態と各種心機能パラメータの同時観測が可能となる。
【０００７】
一方、当該被検体の３次元走査によって収集された心臓領域のボリュームデータの中から、所望関心領域のボリュームデータを切り出すことにより心筋組織等の３次元的な形態観察を容易に行なう方法も開発され、例えば、心臓領域のボリュームデータを所定サイズの直方体表面や任意閉曲面等に基づいて切出す方法が実用化されている。
【特許文献１】特開２００２−３０６４８３号公報
【発明の開示】
【発明が解決しようとする課題】
【０００８】
上述したワイヤフレームのメッシュパターンに心機能パラメータデータを重畳させて表示する従来の方法では組織画像データの表示は行なわれないため、トラッキング状態の把握において現実感が損なわれ、又、前記メッシュパターンにＭＰＲ画像データ等の組織画像データを重畳させて表示する場合には、心筋組織の変位をトラッキングしているにもかかわらずトラッキング状態を示すメッシュパターンと組織画像データを同時に観察することができない。このため、トラッキング状態の全体像を捉えることが困難となり、更に、心筋組織の変位とメッシュパターンの移動に基づくトラッキングとの関連を表示部において直接観察することができないという問題点を有していた。
【０００９】
又、動きの激しい心臓領域から収集されたボリュームデータの中から心筋組織のボリュームデータを切出して３次元画像データ等の組織画像データや心機能パラメータデータ等を生成する場合、心内膜輪郭面及び心外膜輪郭面を含む閉曲面を基準心拍時相にて収集されたボリュームデータの心筋組織に設定しても、前記基準心拍時相に後続する他の心拍時相にて収集されたボリュームデータの心筋組織と前記閉曲面との間に顕著な位置ズレが発生するため、心筋組織のボリュームデータを全ての心拍時相において正確に切出すことが困難となる。このため、所望心拍周期における心筋組織の組織画像データや心機能パラメータデータを正確に収集することができないという問題点を有していた。
【００１０】
本発明は、このような従来の問題点に鑑みてなされたものであり、その目的は、被検体の心臓領域に対する３次元走査によって収集された時系列的なボリュームデータに基づいて心筋組織の組織画像データや心機能パラメータデータを収集する際、心拍動に伴なって変動する前記ボリュームデータをトラッキング処理することにより、心臓の心内膜輪郭面や心外膜輪郭面、あるいは、心内膜輪郭面と心外膜輪郭面に囲まれた心筋組織の組織画像データや心機能パラメータデータを安定かつ正確に収集することが可能な３次元超音波診断装置を提供することにある。
【課題を解決するための手段】
【００１１】
上記課題を解決するために、請求項１に係る本発明の３次元超音波診断装置は、被検体の心臓領域に対する３次元走査によって得られた受信信号に基づいて複数の心拍時相におけるボリュームデータを時系列的に生成するボリュームデータ生成手段と、基準心拍時相のボリュームデータにおける心筋組織の輪郭面を検出し、前記基準心拍時相のボリュームデータと他の心拍時相のボリュームデータとのトラッキング処理により前記他の心拍時相のボリュームデータにおける心筋組織の輪郭面を検出する輪郭面検出手段と、前記輪郭面における前記ボリュームデータのボクセルあるいは前記輪郭面を含む閉曲面によって囲まれた前記ボリュームデータのボクセルを切り出して各心拍時相における組織画像データを生成する組織画像データ生成手段と、前記組織画像データを処理して診断用画像データを生成する診断用画像データ生成手段と前記診断用画像データを表示する表示手段とを備えたことを特徴としている。
【発明の効果】
【００１２】
本発明によれば、被検体の心臓領域に対する３次元走査によって収集された時系列的なボリュームデータに基づいて心筋組織の組織画像データや心機能パラメータデータを収集する際、心拍動に伴なって変動する前記ボリュームデータをトラッキング処理することにより、心臓の心内膜輪郭面や心外膜輪郭面、あるいは、心内膜輪郭面と心外膜輪郭面に囲まれた心筋組織の組織画像データや心機能パラメータデータを安定かつ正確に収集することができる。このため、前記被検体に対する診断精度と診断効率を大幅に向上させることができる。
【発明を実施するための最良の形態】
【００１３】
以下、図面を参照して本発明の実施例を説明する。
【実施例１】
【００１４】
本発明の第１の実施例における３次元超音波診断装置は、被検体の心臓領域に対する３次元走査によって収集された時系列的なボリュームデータをトラッキング処理することにより各心拍時相における心筋組織の心内膜輪郭面及び心外膜輪郭面を検出し、これらの輪郭面に配置した複数の追跡点における心拍時相間の変位を計測することによって得られた心機能パラメータデータと前記輪郭面において形成したワイヤフレームのメッシュパターンと前記輪郭面における前記ボリュームデータのボクセルを切り出して生成した組織画像データとしてのＭＰＲ画像データを重畳して診断用画像データを生成し、得られた各心拍時相の診断用画像データを表示部において動画像表示する。
【００１５】
（装置の構成と機能）
本発明の第１の実施例における３次元超音波診断装置の構成と機能につき図１乃至図６を用いて説明する。尚、図１は、３次元超音波診断装置の全体構成を示すブロック図であり、図２及び図４は、この３次元超音波診断装置が備える送受信部／受信信号処理部及びボリュームデータ生成部の具体的な構成を示すブロック図である。
【００１６】
図１に示す３次元超音波診断装置１００は、被検体の３次元的な心臓領域に対し超音波パルス（送信超音波）を送信し、この送信によって得られた超音波反射波（受信超音波）を電気信号（受信信号）に変換する複数個の振動素子が２次元配列された超音波プローブ３と、前記被検体の所定方向に対し超音波パルスを送信するための駆動信号を前記振動素子に供給し、これらの振動素子から得られた複数チャンネルの受信信号を整相加算する送受信部２と、整相加算後の受信信号を処理してＢモードデータを生成する受信信号処理部４と、前記心臓領域に対する３次元走査によって得られたＢモードデータを超音波送受信方向に対応させて配列することによりボリュームデータを生成し、更に、当該被検体の、例えば、１心拍周期にて生成された心拍時相Ｐ１乃至Ｐｎのボリュームデータを保存するボリュームデータ生成部５を備えている。
【００１７】
又、３次元超音波診断装置１００は、後述の入力部１２から供給される心内膜輪郭線及び心外膜輪郭線の設定情報に基づいて基準心拍時相（例えば、心拍時相Ｐ１）のボリュームデータにおける心内膜輪郭面及び心外膜輪郭面を検出し、更に、基準心拍時相のボリュームデータと心拍時相Ｐ２乃至Ｐｎのボリュームデータとのトラッキング処理によりこれらの心拍時相における心内膜輪郭面及び心外膜輪郭面を検出する輪郭面検出部６と、輪郭面検出部６において検出された心内膜輪郭面及び心外膜輪郭面に対してワイヤフレームを形成するワイヤフレーム形成部７と、輪郭面検出部６において検出された心内膜輪郭面及び心外膜輪郭面の位置情報に基づき、前記ボリュームデータ生成部５が生成した前記心臓領域のボリュームデータの中から心内膜輪郭面及び心外膜輪郭面に対応したボリュームデータのボクセルを切り出して心拍時相Ｐ１乃至ＰｎのＭＰＲ（Multi Planar Reconstruction）画像データを組織画像データとして生成する組織画像データ生成部８と、上述の心内膜輪郭面及び心外膜輪郭面に配置された複数の追跡点（トラッキングポイント）において心機能評価に有効な「歪み（ストレイン）」等の心機能パラメータを計測する心機能パラメータ計測部９を備えている。
【００１８】
更に、３次元超音波診断装置１００は、ワイヤフレーム形成部７にて形成されたワイヤフレームのメッシュパターンと心機能パラメータ計測部９にて得られた心機能パラメータデータを組織画像データ生成部８から供給されるＭＰＲ画像データに重畳して心拍時相Ｐ１乃至Ｐｎの診断用画像データを生成する診断用画像データ生成部１０と、得られたこれらの診断用画像データを心拍時相に対応させて繰り返し（ループ）表示する表示部１１と、被検体情報の入力、ボリュームデータ生成条件の設定、ボリュームデータに対するＭＰＲ断面の設定、前記ＭＰＲ断面にて収集されたＭＰＲ画像データに対する心内膜輪郭線及び心外膜輪郭線の設定、各種コマンド信号の入力等を行なう入力部１２と、心臓領域に対する３次元走査を制御する走査制御部１３と、当該被検体の心拍時相を計測する心拍時相計測ユニット１４と、３次元超音波診断装置１００が有する上述の各ユニットを統括的に制御するシステム制御部１５を備えている。
【００１９】
超音波プローブ３は、２次元配列されたＭ個の図示しない振動素子をその先端部に有し、前記先端部を被検体の体表に接触させて超音波の送受信を行なう。振動素子は電気音響変換素子であり、送信時には電気パルス（駆動信号）を超音波パルス（送信超音波）に変換し、又、受信時には超音波反射波（受信超音波）を電気的な受信信号に変換する機能を有している。そして、これら振動素子の各々は、図示しないＭチャンネルの多芯ケーブルを介して送受信部２に接続されている。尚、本実施例では、Ｍ個の振動素子が２次元配列されたセクタ走査用の超音波プローブ３を用いた３次元超音波診断装置１００について述べるが、リニア走査やコンベックス走査等に対応した超音波プローブを用いても構わない。
【００２０】
次に、図２に示す送受信部２は、超音波プローブ３の振動素子に対して駆動信号を供給する送信部２１と、振動素子から得られた受信信号に対して整相加算（所定方向から得られる受信信号の位相を一致させて加算）を行なう受信部２２を備えている。
【００２１】
送信部２１は、レートパルス発生器２１１、送信遅延回路２１２及び駆動回路２１３を備え、レートパルス発生器２１１は、送信超音波の繰り返し周期を決定するレートパルスを生成して送信遅延回路２１２に供給する。送信遅延回路２１２は、送信に使用されるＭｔ個の振動素子と同数の独立な遅延回路から構成され、送信超音波を撮影領域内の所定の深さに集束するための集束用遅延時間と所定方向（θｐ、φｑ）に送信するための偏向用遅延時間を上記レートパルスに与えて駆動回路２１３へ供給する。駆動回路２１３は、送信遅延回路２１２と同数の独立な駆動回路を有し、超音波プローブ３にて２次元配列されたＭ個の振動素子の中から送信用として選択されたＭｔ（Ｍｔ≦Ｍ）個の振動素子を駆動し、被検体の体内に送信超音波を放射する。
【００２２】
一方、受信部２２は、超音波プローブ３に内蔵されたＭ個の振動素子の中から受信用として選択されたＭｒ（Ｍｒ≦Ｍ）個の振動素子に対応するＭｒチャンネルのＡ／Ｄ変換器２２１及び受信遅延回路２２２と加算器２２３を備えており、受信用の振動素子から供給されたＭｒチャンネルの受信信号は、Ａ／Ｄ変換器２２１にてデジタル信号に変換され、受信遅延回路２２２に送られる。
【００２３】
受信遅延回路２２２は、撮影領域内の所定深さにて反射した受信超音波を集束するための集束用遅延時間と、所定方向（θｐ、φｑ）に対して受信指向性を設定するための偏向用遅延時間をＡ／Ｄ変換器２２１から出力されるＭｒチャンネルの受信信号の各々に与え、加算器２２３は、受信遅延回路２２２からの受信信号を加算する。即ち、受信遅延回路２２２と加算器２２３により、所定方向から得られた受信信号は整相加算される。又、受信部２２の受信遅延回路２２２及び加算器２２３は、その遅延時間の制御によって複数方向に対する受信指向性を同時に形成する所謂並列同時受信を可能とし、並列同時受信の適用により３次元走査に要する時間は大幅に短縮される。尚、上述の送受信部２が備える送信部２１及び受信部２２の一部は超音波プローブ３の内部に設けられていても構わない。
【００２４】
図３は、超音波プローブ３の中心軸をｚ軸とした直交座標（ｘ−ｙ−ｚ）における超音波送受信方向（θｐ、φｑ）を示したものであり、振動素子はｘ軸方向及びｙ軸方向に２次元配列され、θｐ及びφｑは、ｘ−ｚ平面及びｙ−ｚ平面に投影された超音波送受信方向のｚ軸に対する角度を示している。そして、走査制御部１３から供給される走査制御信号に従って送信部２１の送信遅延回路２１２及び受信部２２の受信遅延回路２２２における遅延時間が制御され、被検体内の心臓領域に対する３次元走査が行なわれる。
【００２５】
図２へ戻って、受信信号処理部４は、超音波データとしてのＢモードデータを生成する機能を有し包絡線検波器４１と対数変換器４２を備えている。包絡線検波器４１は、受信部２２の加算器２２３から供給される整相加算後の受信信号を包絡線検波し、対数変換器４２は、包絡線検波された受信信号の振幅を対数変換してＢモードデータを生成する。尚、包絡線検波器４１と対数変換器４２は順序を入れ替えて構成してもよい。
【００２６】
次に、図１に示したボリュームデータ生成部５の具体的な構成につき図４のブロック図を用いて説明する。
【００２７】
ボリュームデータ生成部５は、図４に示すように、超音波データ記憶部５１、補間処理部５２及びボリュームデータ記憶部５３を備え、超音波データ記憶部５１には、当該被検体に対する３次元走査によって得られた受信信号に基づき受信信号処理部４が生成した所定心拍時相のＢモードデータが超音波送受信方向（θｐ、φｑ）を付帯情報として順次保存される。
【００２８】
一方、補間処理部５２は、超音波データ記憶部５１から読み出したＢモードデータを超音波送受信方向（θｐ、φｑ）に対応させて配列することにより３次元Ｂモードデータを生成し、更に、この３次元Ｂモードデータを構成する不等間隔のボクセルを補間処理して、例えば、図３のｘ方向、ｙ方向及びｚ方向に対し等方的なボクセルで構成されるボリュームデータを生成する。そして、心拍時相Ｐ１乃至Ｐｎにおいて得られた超音波ボリュームデータは、心拍時相計測ユニット１４から供給される当該被検体の心拍時相（即ち、心拍時相Ｐ１乃至Ｐｎ）を付帯情報としてボリュームデータ記憶部５３に保存される。
【００２９】
再び図１へ戻って、輪郭面検出部６は、図示しない補間処理部、追跡点設定部及びトラッキング処理部を備え、心拍時相Ｐ１のボリュームデータを用いて組織画像データ生成部８が生成し表示部１１において表示されたＭＰＲ画像データに対し入力部１２が設定した心内膜輪郭線及び心外膜輪郭線の位置情報を受信する。この場合、例えば、心臓の２腔断面に対応したＭＰＲ断面や４腔断面に対応したＭＰＲ断面においてＭＰＲ画像データが生成され、これらＭＰＲ画像データの心内膜及び心外膜に対して入力部１２が設定した心内膜輪郭線及び心外膜輪郭線の位置情報がシステム制御部１５を介して輪郭面検出部６に供給される。そして、輪郭面検出部６の前記補間処理部は、ボリュームデータ生成部５のボリュームデータ記憶部５３から供給される心拍時相Ｐ１（基準心拍時相）のボリュームデータに上述の輪郭線を設定し、これら輪郭線の各々に対してスプライン補間等の補間処理を行なうことにより心拍時相Ｐ１のボリュームデータにおける心内膜輪郭面及び心外膜輪郭面を検出する。
【００３０】
一方、前記追跡点設定部は、上述の心内膜輪郭面及び心外膜輪郭面に複数の追跡点を配置し、前記トラッキング処理部は、これら追跡点の各々を中心とした３次元関心領域における心拍時相Ｐ１のボリュームデータと心拍時相Ｐ２のボリュームデータとをトラッキング処理することによって心拍時相Ｐ２のボリュームデータにおける心内膜輪郭面及び心外膜輪郭面を検出する。更に、前記トラッキング処理部は、同様の手順により心拍時相Ｐ１のボリュームデータと心拍時相Ｐ３乃至Ｐｎのボリュームデータとのトラッキング処理により心拍時相Ｐ３乃至Ｐｎのボリュームデータにおける心内膜輪郭面及び心外膜輪郭面を検出する。
【００３１】
但し、心拍時相Ｐ１のボリュームデータと心拍時相Ｐ３乃至Ｐｎのボリュームデータとのトラッキング処理の替わりに、隣接する心拍時相のボリュームデータに対してトラッキング処理を行なうことによって心拍時相Ｐ３乃至Ｐｎのボリュームデータにおける心内膜輪郭面及び心外膜輪郭面を検出してもよい。尚、上述のトラッキング処理としてボリュームデータ間のパターンマッチングが通常行なわれるが他の方法であってもよく、特に限定されない。
【００３２】
次に、相互相関処理を適用したパターンマッチングの具体例につき図５を用いて説明する。図５において、図示しない心内膜輪郭面あるいは心外膜輪郭面の追跡点Ｐｔを中心に設定された３次元関心領域Ｒｏにおける心拍時相ＰａのボリュームデータＶａが有するボクセル値をＡ（ｆ、ｇ、ｈ）（ｆ＝１〜Ｆ，ｇ＝１〜Ｇ、ｈ＝１〜Ｈ）、心拍時相Ｐａに後続する心拍時相ＰｂのボリュームデータＶｂが有するボクセル値をＢ（ｆ、ｇ、ｈ）（ｆ＝ｊ＋１〜ｊ＋Ｆ，ｇ＝ｋ＋１〜ｋ＋Ｇ、ｈ＝ｓ＋１〜ｓ＋Ｈ）、とすれば、心拍時相Ｐａから心拍時相Ｐｂまでの間に追跡点Ｐｔが変位する方向と距離は次式（１）の変位評価関数γＡＢ（ｊ、ｋ、ｓ）によって計測される。
【数１】
2010115372.tif 000003
【００３３】
そして、３次元関心領域ＲｏのボリュームデータＶａをボリュームデータＶｂに対しｐ方向、ｑ方向及びｈ方向に逐次移動させながら上式（１）の変位評価関数γＡＢ（ｊ、ｋ、ｓ）を算出し、ｊ＝ｊｘ、ｋ＝ｋｘ及びｓ＝ｓｘにおいてγＡＢ（ｊ、ｋ、ｓ）がピーク値を示した場合、心拍時相Ｐａにおける追跡点Ｐｔは心拍時相Ｐｂにおいてｆ方向にｊｘボクセル分、ｇ方向にｋｘボクセル分、ｈ方向にｓｘボクセル分だけ変位したものとして検出される。但し、上式（１）におけるＮは、３次元関心領域Ｒｏにおけるｆ方向のボクセル数Ｆ、ｇ方向のボクセル数Ｇ及びｈ方向のボクセル数Ｈの積を示している。
【００３４】
一方、ワイヤフレーム形成部７は、心拍時相Ｐ１乃至Ｐｎのボリュームデータにおける心内膜輪郭面及び心外膜輪郭面に対して配置された上述の追跡点を所定間隔で抽出して結線することによりワイヤフレームを形成する。図６は、ワイヤフレーム形成部７によって形成されたワイヤフレームの具体例を示したものであり、このワイヤフレームＷｆは、心内膜輪郭面あるいは心外膜輪郭面Ｃｐにおいて配置された図示しない複数の追跡点の中から所定間隔で抽出した追跡点Ｔｐを交点として形成される。例えば、心機能パラメータの計測を目的として上述の輪郭面に配置された複数の追跡点の中から所定間隔で間引かれた数分の１乃至数十分の１の追跡点を結線することによってワイヤフレームＷｆが形成される。
【００３５】
図１へ戻って、組織画像データ生成部８は、ボリュームデータ生成部５のボリュームデータ記憶部５３から供給される基準心拍時相（心拍時相Ｐ１）のボリュームデータに入力部１２から供給される心臓の２腔断面や４腔断面に対応したＭＰＲ断面の位置情報を設定し、これらのＭＰＲ断面に対応する前記ボリュームデータのボクセルを切り出して心内膜輪郭線及び心外膜輪郭線の設定を目的としたＭＰＲ画像データを生成する。
【００３６】
更に、組織画像データ生成部８は、ボリュームデータ記憶部５３から供給された心拍時相Ｐ１乃至Ｐｎのボリュームデータに輪郭面検出部６から供給された心拍時相Ｐ１乃至Ｐｎにおける心内膜輪郭面及び心外膜輪郭面の位置情報を設定し、心内膜輪郭面あるいは心外膜輪郭面に対応する前記ボリュームデータのボクセルを切り出して曲面状のＭＰＲ画像データを生成する。尚、この場合、心内膜輪郭面あるいは心外膜輪郭面に垂直な方向に対しても所定の厚み範囲で複数のボクセルを切り出し、これらのボクセルが有するボクセル値を厚み方向に演算処理（例えば、加算平均処理）して心内膜輪郭面や心外膜輪郭面におけるＭＰＲ画像データを生成してもよい。
【００３７】
次に、心機能パラメータ計測部９は、図示しない演算処理部と不透明度・色相設定部を有し、前記演算処理部は、輪郭面検出部６から供給された心拍時相Ｐ１乃至Ｐｎの心内膜輪郭面及び心外膜輪郭面における追跡点の位置情報に基づいて心内膜輪郭面及び心外膜輪郭面における「変位」を計測し、更に、単位長さ当たりの変位で定義される「歪み（ストレイン）」を心機能パラメータとして計測する。一方、前記不透明度・色相設定部は、前記心機能パラメータの計測結果に対し所定の色相と不透明度を設定して心機能パラメータデータを生成する。尚、心機能パラメータ計測部９の演算処理部が追跡点の位置情報に基づいて計測する心機能パラメータは上述の「歪み（ストレイン）」に限定されるものではなく、例えば、「変位」、「回転」、「捻れ（torsion）」、「速度」等を心機能パラメータとして計測してもよく、これらの時間的変化を示す「歪みレート」、「回転レート」、「捻れレート」、「加速度」等であっても構わない。
【００３８】
一方、診断用画像データ生成部１０は、図示しない画像データ合成部と画像データ記憶部を備え、前記画像データ合成部は、心拍時相Ｐ１乃至Ｐｎにて収集されたボリュームデータに基づきワイヤフレーム形成部７が形成したワイヤフレームのメッシュパターンと心機能パラメータ計測部９にて得られた心機能パラメータデータを組織画像データ生成部８が生成したＭＰＲ画像データに重畳して心拍時相Ｐ１乃至Ｐｎの診断用画像データを生成する。この場合、心機能パラメータの計測結果に所定の色相と透明度を設定して生成された心機能パラメータデータをＭＰＲ画像データに重畳することにより組織画像データ及び心機能パラメータデータの各々を分離して観察することが可能となる。そして、得られたこれらの診断用画像データは、心拍時相を付帯情報として前記画像データ記憶部に保存される。
【００３９】
表示部１１は、図示しない表示制御部、画像データ変換部及びモニタを備え、前記表示制御部は、診断用画像データ生成部１０の画像データ記憶部に保存された診断用画像データを心拍時相に対応させて順次読み出し、前記画像データ変換部は、これらの診断用画像データに対し所定の表示フォーマットに基づいた座標変換やＤ／Ａ変換等の変換処理を行なって前記モニタに表示する。即ち、前記モニタには、変換処理された心拍時相Ｐ１乃至Ｐｎの診断用画像データが繰り返し（ループ）表示される。
【００４０】
入力部１２は、操作パネル上にキーボード、トラックボール、マウス、選択ボタン、入力ボタン等の入力デバイスや表示パネルを備え、心内膜輪郭線及び心外膜輪郭線の設定を目的としたＭＰＲ画像データの収集断面（ＭＰＲ断面）を設定するＭＰＲ断面設定機能１２１や前記ＭＰＲ画像データに対して心内膜輪郭線及び心外膜輪郭線を設定する輪郭設定機能１２２を有している。又、被検体情報の入力、ボリュームデータ生成条件の設定、更には、各種コマンド信号の入力等も上述の表示パネルや入力デバイスを用いて行なわれる。
【００４１】
走査制御部１３は、当該被検体の心臓領域に対し超音波による３次元走査を行なうための遅延時間制御を送信部２１の送信遅延回路２１２及び受信部２２の受信遅延回路２２２に対して行なう。心拍時相計測ユニット１４は、被検体の心電波形に基づいて心拍時相情報を発生する機能を有し、被検体体表面に装着され心電波形を計測する計測用電極と、この計測用電極によって計測された心電波形を所定の振幅に増幅する増幅回路と、増幅された心電波形をデジタル信号に変換するＡ／Ｄ変換器と、デジタル変換された心電波形に基づいて、例えば、心臓の拡張末期を基準とする１心拍周期に対して心拍時相Ｐ１乃至Ｐｎを設定する心拍時相設定部（何れも図示せず）を備えている。
【００４２】
そして、システム制御部１５は、図示しないＣＰＵと記憶回路を備え、入力部１２において入力／設定された各種の情報は前記記憶回路に保存される。そして、前記ＣＰＵは、入力部１２から入力された上述の入力／設定情報や自己の記憶回路に予め保管された各種の入力／設定情報に基づいて３次元超音波診断装置１００の各ユニットを統括的に制御し、心機能パラメータデータやワイヤフレームのメッシュパターンが重畳されたＭＰＲ画像データを診断用画像データとして動画像表示する。
【００４３】
（診断用画像データの生成／表示手順）
次に、本実施例における診断用画像データの生成／表示手順につき図７のフローチャートを用いて説明する。
【００４４】
診断用画像データの生成に先立ち、３次元超音波診断装置１００の操作者は、心拍時相計測ユニット１４が備えた計測用電極を被検体の胸部体表面に装着した後、入力部１２において被検体情報の入力、ボリュームデータ生成条件の設定、診断用画像データ生成条件の設定等を行なう。そして、これらの入力情報や設定情報は、システム制御部１５の記憶回路に保存される（図７のステップＳ１）。
【００４５】
上述の初期設定が終了したならば、操作者は、心拍時相計測ユニット１４を動作状態に切り替えることによって当該被検体に対する心拍時相の計測が開始され（図７のステップＳ２）、更に、入力部１２においてボリュームデータの収集開始コマンドを入力することによって当該被検体の心臓領域に対するボリュームデータの収集が開始される（図７のステップＳ３）。
【００４６】
ボリュームデータの収集に際し、走査制御部１３は、送受信部２の送信遅延回路２１２及び受信遅延回路２２２に対し当該被検体の心臓領域を３次元走査するための走査制御信号を供給する。
【００４７】
一方、被検体から得られた心電波形に基づいて心拍時相計測ユニット１４が発生した最初の心拍時相（心拍時相Ｐ１）の同期信号を、システム制御部１５を介して受信した送信部２１のレートパルス発生器２１１は、予め設定された繰り返し周期（レート周期）を有するレートパルスを前記同期信号に基づいて生成し送信遅延回路２１２へ供給する。
【００４８】
送信遅延回路２１２は、走査制御部１３から供給された走査制御信号に基づいて当該被検体の心臓領域を含む３次元領域の所定の深さに超音波を集束するための集束用遅延時間と、最初の超音波送受信方向（θ１、φ１）に超音波を送信するための偏向用遅延時間を前記レートパルスに与え、このレートパルスをＭｔチャンネルの駆動回路２１３へ供給する。
【００４９】
次いで、駆動回路２１３は、送信遅延回路２１２から供給されたレートパルスに基づいて駆動信号を生成し、この駆動信号を超音波プローブ３に設けられたＭｔ個の送信用振動素子に供給して被検体内に送信超音波を放射する。放射された送信超音波の一部は、音響インピーダンスの異なる臓器境界面や組織にて反射し、超音波プローブ３に設けられたＭｒ個の受信用振動素子によって受信されてＭｒチャンネルの電気的な受信信号に変換される。次いで、この受信信号は、受信部２２のＡ／Ｄ変換器２２１においてデジタル信号に変換され、更に、Ｍｒチャンネルの受信遅延回路２２２において所定の深さからの受信超音波を収束するための集束用遅延時間と超音波送受信方向（θ１、φ１）からの受信超音波に対し強い受信指向性を設定するための偏向用遅延時間が走査制御部１３から供給された上述の走査制御信号に基づいて与えられた後加算器２２３にて整相加算される。そして、整相加算後の受信信号が供給された受信信号処理部４の包絡線検波器４１及び対数変換器４２は、この受信信号に対して包絡線検波と対数変換を行なってＢモードデータを生成し、得られたＢモードデータは、超音波送受信方向（θ１、φ１）を付帯情報としてボリュームデータ生成部５の超音波データ記憶部５１に保存される。
【００５０】
次いで、走査制御部１３は、送受信部２の送信遅延回路２１２及び受信遅延回路２２２における遅延時間を制御してθ方向にΔθ、φ方向にΔφずつ順次更新された超音波送受信方向（θｐ、φｑ）（θｐ＝θ１＋（ｐ−１）Δθ（ｐ＝１〜Ｐ）、φｑ＝φ１＋（ｑ−１）Δφ（ｑ＝１〜Ｑ）、但し、超音波送受信方向（θ１、φ１）を除く）の各々に対し同様の手順で超音波を送受信して３次元走査を行なう。そして、各々の送受信方向にて得られたＢモードデータも上述の超音波送受信方向を付帯情報として超音波データ記憶部５１に保存される。
【００５１】
心拍時相Ｐ１の心臓領域に対するＢモードデータの収集と保存が終了したならば、ボリュームデータ生成部５の補間処理部５２は、超音波データ記憶部５１から読み出した複数のＢモードデータを超音波送受信方向（θｐ、φｑ）に対応させて配列することにより３次元Ｂモードデータを生成し、更に、この３次元Ｂモードデータを構成する不等間隔のボクセルを補間処理して等方的なボクセルで構成されるボリュームデータを生成する。そして、得られたボリュームデータは、心拍時相計測ユニット１４からシステム制御部１５を介して供給される心拍時相Ｐ１を付帯情報としてボリュームデータ生成部５のボリュームデータ記憶部５３に保存される（図７のステップＳ４）。
【００５２】
同様にして、心拍時相Ｐ２乃至Ｐｎの各々に同期した３次元走査によって収集されたボリュームデータも心拍時相Ｐ２乃至Ｐｎを付帯情報としてボリュームデータ記憶部５３に保存される（図７のステップＳ２乃至Ｓ４）。即ち、ボリュームデータ記憶部５３には、当該被検体の、例えば、拡張末期を基準とする１心拍周期において収集された心拍時相Ｐ１乃至Ｐｎのボリュームデータが夫々の心拍時相を付帯情報として保存される。
【００５３】
当該被検体の１心拍周期を構成する心拍時相Ｐ１乃至Ｐｎのボリュームデータが保存されたならば、組織画像データ生成部８は、ボリュームデータ生成部５のボリュームデータ記憶部５３から読み出した基準心拍時相（心拍時相Ｐ１）のボリュームデータに対し入力部１２から供給された心臓の２腔断面や４腔断面に対応するＭＰＲ断面の位置情報を設定し、これらのＭＰＲ断面に対応する前記ボリュームデータのボクセルを切り出してＭＰＲ画像データを生成する。そして、診断用画像データ生成部１０は、これらのＭＰＲ画像データを合成して表示部１１のモニタに表示する（図７のステップＳ５）。
【００５４】
次いで、操作者は、入力部１２に備えられた入力デバイスを用い表示部１１に表示されている複数からなるＭＰＲ画像データの心内膜及び心外膜に対して心内膜輪郭線及び心外膜輪郭線を設定する（図７のステップＳ６）。そして、入力部１２から心内膜輪郭線及び心外膜輪郭線の位置情報を受信した輪郭面検出部６の補間処理部は、ボリュームデータ生成部５のボリュームデータ記憶部５３から供給された心拍時相Ｐ１のボリュームデータに心内膜輪郭線及び心外膜輪郭線を設定し、これらの輪郭線に対してスプライン補間等の補間処理を行なうことにより心内膜輪郭面及び心外膜輪郭面を検出する。
【００５５】
一方、輪郭面検出部６の追跡点設定部は、検出された心内膜輪郭面及び心外膜輪郭面の面上に複数の追跡点を配置し、輪郭面検出部６のトラッキング処理部は、これら追跡点の各々を中心とした３次元関心領域における心拍時相Ｐ１のボリュームデータと心拍時相Ｐ２のボリュームデータとのトラッキング処理によって心拍時相Ｐ２のボリュームデータにおける心内膜輪郭面及び心外膜輪郭面を検出する。更に、前記トラッキング処理部は、同様の手順により前記心拍時相Ｐ１のボリュームデータと心拍時相Ｐ３乃至Ｐｎのボリュームデータとのトラッキング処理により心拍時相Ｐ３乃至Ｐｎのボリュームデータにおける心内膜輪郭面及び心外膜輪郭面を検出する（図７のステップＳ７）。
【００５６】
次に、ワイヤフレーム形成部７は、心拍時相Ｐ１乃至Ｐｎのボリュームデータにおける心内膜輪郭面及び心外膜輪郭面に配置された複数からなる追跡点を所定間隔で抽出し、抽出された追跡点を結線することによりワイヤフレームを形成する（図７のステップＳ８）。一方、組織画像データ生成部８は、ボリュームデータ記憶部５３から供給された心拍時相Ｐ１乃至Ｐｎのボリュームデータに上述の心拍時相Ｐ１乃至Ｐｎの心内膜輪郭面及び心外膜輪郭面を設定し、心内膜輪郭面あるいは心外膜輪郭面に対応する前記ボリュームデータのボクセルを切り出して曲面状のＭＰＲ画像データを組織画像データとして生成する（図７のステップＳ９）。又、心機能パラメータ計測部９は、輪郭面検出部６から供給された心拍時相Ｐ１乃至Ｐｎの心内膜輪郭面及び心外膜輪郭面における追跡点の位置情報に基づいて心内膜輪郭面及び心外膜輪郭面における「変位」を計測し、更に、単位長さ当たりの変位で定義される「歪み（ストレイン）」を心機能パラメータとして計測する（図７のステップＳ１０）。
【００５７】
次に、診断用画像データ生成部１０の画像データ合成部は、心拍時相Ｐ１乃至Ｐｎにて収集されたボリュームデータに基づいてワイヤフレーム形成部７が形成したワイヤフレームのメッシュパターンと心機能パラメータ計測部９にて得られた心機能パラメータデータを組織画像データ生成部８が生成したＭＰＲ画像データに重畳して心拍時相Ｐ１乃至Ｐｎの診断用画像データを生成する。そして、得られたこれらの診断用画像データは心拍時相を付帯情報として自己の画像データ記憶部に保存される（図７のステップＳ１１）。
【００５８】
次いで、表示部１１は、診断用画像データ生成部１０の画像データ記憶部に保存された診断用画像データを心拍時相に対応させて順次読み出し、所定の変換処理を行なって自己のモニタに繰り返し（ループ）表示する（図７のステップＳ１２）。
【００５９】
以上述べた第１の実施例によれば、被検体の心臓領域に対する３次元走査によって収集された時系列的なボリュームデータに基づいて心筋組織のＭＰＲ画像データや心機能パラメータデータを収集する際、心拍動に伴なって変動する前記ボリュームデータをトラッキング処理することにより、心臓の心内膜輪郭面及び心外膜輪郭面におけるＭＰＲ画像データや心機能パラメータデータを安定かつ正確に収集することができる。このため、前記被検体に対する診断精度と診断効率を大幅に向上させることができる。
【００６０】
又、上述の実施例によれば、同一の心拍時相において収集された心内膜輪郭面及び心外膜輪郭面における心機能パラメータデータと組織画像データとしてのＭＰＲ画像データを重畳させて動画像表示しているため、心筋表面の状態と心筋表面における心機能パラメータを同時に観察することができる。特に、心機能パラメータの計測結果に対し好適な色相と透明度を設定して生成された心機能パラメータデータをＭＰＲ画像データに重畳することによりＭＰＲ画像データと心機能パラメータデータの各々を分離させて観察することが可能となる。
【００６１】
更に、心内膜輪郭面及び心外膜輪郭面において形成されたワイヤフレームのメッシュパターンを上述の心機能パラメータやＭＰＲ画像データに重畳させて動画像表示しているため、トラッキング時における心筋組織の運動状態をＭＰＲ画像データにおいて把握することができ、メッシュパターンとＭＰＲ画像データとの比較観察によりボリュームデータ間のトラッキング処理が正確に行なわれているか否かを容易に確認することができる。
【００６２】
尚、上述の第１の実施例における組織画像データ生成部８は、ボリュームデータ記憶部５３から供給されたボリュームデータに輪郭面検出部６から供給された心内膜輪郭面及び心外膜輪郭面の位置情報を設定し、これらの輪郭面に対応する前記ボリュームデータのボクセルを切り出して曲面状のＭＰＲ画像データを生成する場合について述べたが、図８に示すように、ワイヤフレーム形成部７が形成したワイヤフレームＷｆ（図６参照）の追跡点Ｔｐによって囲まれた複数からなる微小平面Ｐｆの位置情報を前記ボリュームデータに設定し、微小平面Ｐｆに対応した前記ボリュームデータのボクセルを切り出して生成した微小ＭＰＲ画像データＤｍをワイヤフレームＷｆの微小平面Ｐｆに貼付してもよい。
【００６３】
又、第１の実施例及び上述の変形例におけるＭＰＲ画像データの生成に際し、ＭＰＲ画像データをフィルタリング処理することによりこの画像データに含まれたスペックルの強調表示を可能とする画像処理部を組織画像データ生成部８に備えてもよい。スペックルが強調されたＭＰＲ画像データに心機能パラメータデータやワイヤフレームのメッシュパターンを重畳することにより心筋組織の運動状態やトラッキング状態の把握が更に容易となる。この場合のフィルタリング処理の選択やフィルタリング特性の設定は入力部１２によって行なうことが可能であるが特に限定されない。
【実施例２】
【００６４】
次に本発明の第２の実施例について説明する。上述の第１の実施例では、組織画像データとして心内膜輪郭面及び心外膜輪郭面に沿ったＭＰＲ画像データを生成する場合について述べたが、本実施例では、心内膜輪郭面及び心外膜輪郭面によって囲まれたボリュームデータをレンダリング処理して３次元画像データを生成する場合について述べる。
【００６５】
即ち、本実施例における３次元超音波診断装置は、被検体の心臓領域に対する３次元走査によって収集された時系列的なボリュームデータをトラッキング処理することにより各心拍時相における心筋の心内膜輪郭面及び心外膜輪郭面を検出し、これらの輪郭面に配置した複数の追跡点における心拍時相間の変位を計測することによって得られる心機能パラメータデータと前記輪郭面において形成されるワイヤフレームのメッシュパターンと前記輪郭面を含む閉曲面を用いて切出したボリュームデータのレンダリング処理によって得られる３次元画像データを重畳して診断用画像データを生成し、得られた各心拍時相の診断用画像データを表示部において動画像表示する。
【００６６】
（装置の構成と機能）
本実施例における３次元超音波診断装置２００の構成と機能につき図９及び図１０を用いて説明する。但し、３次元超音波診断装置２００の全体構成を示す図９のブロック図において、図１に示した３次元超音波診断装置１００が備えるユニットと同様な構成と機能を有するユニットは同一の符号を付加し詳細な説明は省略する。
【００６７】
即ち、図９に示す本実施例の３次元超音波診断装置２００は、３次元超音波診断装置１００と同様な送受信部２、超音波プローブ３、受信信号処理部４、ボリュームデータ生成部５、表示部１１、入力部１２、走査制御部１３、心拍時相計測ユニット１４及びシステム制御部１５を備え、更に、入力部１２から供給される心内膜輪郭線及び心外膜輪郭線の設定情報に基づき、基準心拍時相（心拍時相Ｐ１）のボリュームデータに対し心内膜輪郭面及び心外膜輪郭面を含んだ閉曲面を輪郭面として検出し、更に、前記基準心拍時相に後続する心拍時相Ｐ２乃至Ｐｎのボリュームデータをトラッキング処理することにより前記ボリュームデータにおける閉曲面を検出する輪郭面検出部６ａと、輪郭面検出部６ａにおいて検出された閉曲面に対してワイヤフレームを形成するワイヤフレーム形成部７ａと、前記ボリュームデータ生成部５が生成した心臓領域のボリュームデータの中から上述の閉曲面によって囲まれたボリュームデータを切り出して心拍時相Ｐ１乃至Ｐｎにおける３次元画像データを組織画像データとして生成する組織画像データ生成部８ａと、前記閉曲面の面上に配置された複数の追跡点において心機能評価に有効な「歪み（ストレイン）」等の心機能パラメータを計測する心機能パラメータ計測部９ａを備えている。
【００６８】
又、３次元超音波診断装置２００は、ワイヤフレーム形成部７ａにおいて形成されたワイヤフレームのメッシュパターンと心機能パラメータ計測部９ａにおいて得られた心機能パラメータデータを組織画像データ生成部８ａから供給される３次元画像データに重畳して心拍時相Ｐ１乃至Ｐｎにおける診断用画像データを生成する診断用画像データ生成部１０ａを備えている。
【００６９】
輪郭面検出部６ａは、図示しない補間処理部、追跡点設定部及びトラッキング処理部を備え、前記補間処理部は、ボリュームデータ生成部５のボリュームデータ記憶部５３から供給される心拍時相Ｐ１のボリュームデータに入力部１２から供給される複数の輪郭線情報を設定し、これらの輪郭線を補間処理してボリュームデータにおける心内膜輪郭面及び心外膜輪郭面とその端部によって形成される閉曲面を検出する。
【００７０】
一方、前記追跡点設定部は、上述の閉曲面に複数の追跡点を配置し、前記トラッキング処理部は、これら追跡点の各々を中心とした３次元関心領域における心拍時相Ｐ１のボリュームデータと心拍時相Ｐ２のボリュームデータとのトラッキング処理によって心拍時相Ｐ２のボリュームデータにおける閉曲面を検出する。更に、前記トラッキング処理部は、同様の手順により心拍時相Ｐ１のボリュームデータと心拍時相Ｐ３乃至Ｐｎのボリュームデータとのトラッキング処理により心拍時相Ｐ３乃至Ｐｎのボリュームデータにおける閉曲面を検出する。
【００７１】
ワイヤフレーム形成部７ａは、心拍時相Ｐ１乃至Ｐｎのボリュームデータにおける閉曲面に対して配置された上述の追跡点を所定間隔で抽出して結線することによりワイヤフレームを形成する。一方、組織画像データ生成部８ａは、図１に示した組織画像データ生成部８と同様にボリュームデータ生成部５のボリュームデータ記憶部５３から供給される基準心拍時相（心拍時相Ｐ１）のボリュームデータに入力部１２から供給されるＭＰＲ断面の位置情報を設定し、これらＭＰＲ断面に対応する前記ボリュームデータのボクセルを切り出して心内膜輪郭線及び心外膜輪郭線の設定を目的としたＭＰＲ画像データを生成する。
【００７２】
更に、組織画像データ生成部８ａは、ボリュームデータ記憶部５３から供給される心拍時相Ｐ１乃至Ｐｎのボリュームデータに輪郭面検出部６ａが検出した心拍時相Ｐ１乃至Ｐｎの閉曲面を設定し、閉曲面の内部から切り出した前記ボリュームデータのボクセルをレンダリング処理して３次元画像データを生成する機能を有している。この組織画像データ生成部８ａは、図示しない不透明度・色相設定部とレンダリング処理部を備えている。そして、前記不透明度・色相設定部は、閉曲面の内部から切り出したボリュームデータのボクセルが有するボクセル値及び前記ボクセルから閉曲面までの距離に基づいて各ボクセルの不透明度を設定し、前記レンダリング処理部は、前記不透明度・色相設定部によって設定された不透明度の情報に基づいて閉曲面内のボリュームデータをレンダリング処理し３次元画像データを生成する。
【００７３】
通常のレンダリング処理では、ボクセル値に基づいて３次元画像データの明度（輝度）や不透明度が設定され、奥行き情報や影付け情報等に基づいて３次元画像データの色相が設定されるが、３次元画像データと心機能パラメータを重畳させて表示する本実施例では、ボリュームデータのボクセル値に基づいて３次元画像データの明度や不透明度を設定し、心機能パラメータの計測値に基づいて心機能パラメータデータの色相を設定する方法が好適である。但し、色相情報を有した３次元画像データと心機能パラメータデータとを重畳させても構わない。
【００７４】
図１０は、閉曲面からボクセルまでの距離に基づいて前記不透明度・色相設定部が設定する不透明度特性の具体例を示したものであり、前記不透明度・色相設定部は、心筋組織の表層を示す閉曲面から深部に向って不透明度特性を任意に設定する機能を有し、この不透明度特性の設定は、システム制御部１５の記憶回路に予め保管された設定情報あるいは入力部１２から供給される設定情報に基づいて行なわれる。
【００７５】
例えば、前記不透明度・色相設定部は、図１０（ａ）に示すように心内膜側の閉曲面Ｃｉ及び心外膜側の閉曲面Ｃｅから深部に向って所定間隔で領域Ｒ１乃至Ｒｍを形成し、各々の領域におけるボクセル値に対し図１０（ｂ）に示すような不透明度特性を設定する。具体的には、同一のボクセル値に対し、閉曲面Ｃｉ及びＣｅに近い領域ほど低い不透明度を設定し、深部の領域に対しては高い不透明度を設定する。組織画像データ生成部８ａの前記不透明度・色相設定部にこのような不透明度特性の設定機能をもたせることにより、閉曲面が心筋組織の心内膜あるいは心外膜を正確に捉えていない場合においても閉曲面を多少厚めに設定しておけば閉曲面近傍の不透明度を低く設定することができるため、前記閉曲面の内部から切り出されたボリュームデータに基づいて生成される３次元画像データにおいて心筋組織の表面（即ち、心内膜及び心外膜）を確実に捉えることが可能となる。
【００７６】
図９に戻って、心機能パラメータ計測部９ａは、図示しない演算処理部と不透明度・色相設定部を有し、前記演算処理部は、輪郭面検出部６ａから供給された心拍時相Ｐ１乃至Ｐｎの閉曲面における追跡点の位置情報に基づいて心筋の心内膜及び心外膜における「変位」を計測し、更に、単位長さ当たりの変位で定義される「歪み（ストレイン）」を心機能パラメータとして計測する。一方、前記不透明度・色相設定部は、前記心機能パラメータの計測結果に所定の色相と不透明度を設定して心機能パラメータデータを生成する。尚、上述の不透明度や色相の設定は、システム制御部１５の記憶回路に予め保管された設定情報あるいは入力部１２から供給される設定情報に基づいて行なわれる。
【００７７】
次に、診断用画像データ生成部１０ａは、心拍時相Ｐ１乃至Ｐｎの各々において収集されたボリュームデータに基づいてワイヤフレーム形成部７ａが形成したワイヤフレームと心機能パラメータ計測部９ａが生成した心機能パラメータデータを組織画像データ生成部８ａが生成した３次元画像データに重畳して心拍時相Ｐ１乃至Ｐｎにおける診断用画像データを生成し、前記心拍時相を付帯情報として自己の画像データ記憶部に保存する。
【００７８】
（診断用画像データの生成／表示手順）
次に、本実施例における診断用画像データの生成／表示手順につき図１１のフローチャートを用いて説明する。但し、図１１のフローチャートにおいて、図７に示した第１の実施例と同一のステップは同一の符号を付加し詳細な説明は省略する。
【００７９】
表示部１１に表示されたＭＰＲ画像データの心内膜及び心外膜に対する心内膜輪郭線及び心外膜輪郭線の設定が図１１のステップＳ６において終了したならば、これら輪郭線の位置情報を受信した輪郭面検出部６ａの補間処理部は、ボリュームデータ生成部５のボリュームデータ記憶部５３から供給された心拍時相Ｐ１のボリュームデータに上述の輪郭線を設定し、これらの輪郭線に対しスプライン補間等の補間処理を行なうことにより心内膜輪郭面及び心外膜輪郭面とその端部によって形成される閉曲面を検出する。
【００８０】
一方、輪郭面検出部６ａの追跡点設定部は、上述の閉曲面に複数の追跡点を配置し、輪郭面検出部６ａのトラッキング処理部は、これら追跡点の各々を中心とした３次元関心領域における心拍時相Ｐ１のボリュームデータと心拍時相Ｐ２のボリュームデータとのトラッキング処理により心拍時相Ｐ２のボリュームデータにおける閉曲面を検出する。更に、前記トラッキング処理部は、同様の手順により心拍時相Ｐ１のボリュームデータと心拍時相Ｐ３乃至Ｐｎのボリュームデータとのトラッキング処理により心拍時相Ｐ３乃至Ｐｎのボリュームデータにおける閉曲面を検出する（図１１のステップＳ７ａ）。
【００８１】
次に、ワイヤフレーム形成部７ａは、心拍時相Ｐ１乃至Ｐｎのボリュームデータにおける上述の閉曲面に配置された複数からなる追跡点を所定間隔で抽出し、抽出した追跡点を結線することによりワイヤフレームを形成する（図１１のステップＳ８ａ）。一方、組織画像データ生成部８ａの不透明度・色相設定部は、前記閉曲面の内部から切り出されたボリュームデータのボクセルが有するボクセル値及び前記ボクセルから閉曲面までの距離に基づいて各ボクセルの不透明度や色相を設定し、組織画像データ生成部８ａのレンダリング処理部は、前記不透明度・色相設定部によって設定された不透明度や色相の情報に基づいて閉曲面内のボリュームデータをレンダリング処理し３次元画像データを組織画像データとして生成する（図１１のステップＳ９ａ）。
【００８２】
又、心機能パラメータ計測部９ａは、輪郭面検出部６ａから供給された心拍時相Ｐ１乃至Ｐｎの閉曲面における追跡点の位置情報に基づいて心筋の心内膜及び心外膜における「変位」を計測し、更に、単位長さ当たりの変位で定義される「歪み（ストレイン）」を心機能パラメータとして計測する（図１１のステップＳ１０ａ）。
【００８３】
次に、診断用画像データ生成部１０ａは、心拍時相Ｐ１乃至Ｐｎの各々において収集されたボリュームデータに基づいてワイヤフレーム形成部７ａが形成したワイヤフレームのメッシュパターンと心機能パラメータ計測部９ａが生成した心機能パラメータデータを組織画像データ生成部８ａが生成した３次元画像データに重畳して心拍時相Ｐ１乃至Ｐｎにおける診断用画像データを生成し、得られたこれらの診断用画像データを、心拍時相を付帯情報として自己の画像データ記憶部に保存する（（図１１のステップＳ１１ａ）。そして、前記画像データ記憶部に保存された診断用画像データは、図７のステップＳ１２と同様に心拍時相に対応させて順次読み出され、表示部１１に繰り返し（ループ）表示される（図１１のステップＳ１２）。
【００８４】
以上述べた第２の実施例によれば、被検体の心臓領域に対する３次元走査によって収集された時系列的なボリュームデータに基づいて心筋組織の組織画像データや心機能パラメータデータを収集する際、心拍動に伴なって変動する前記ボリュームデータをトラッキング処理することにより、心臓の心内膜輪郭面及び心外膜輪郭面を含む閉曲面によって囲まれた心筋組織の組織画像データや心機能パラメータデータを安定かつ正確に収集することができる。このため、前記被検体に対する診断精度と診断効率を大幅に向上させることができる。
【００８５】
又、上述の実施例によれば、同一の心拍時相において収集された前記心筋組織の心機能パラメータデータと組織画像データとしての３次元画像データを重畳させて動画像表示しているため、心筋組織の状態と「歪み（ストレイン）」等の心機能パラメータを同時に観察することができる。特に、本実施例では、閉曲面に囲まれた心筋組織のボリュームデータを切出して３次元画像データを生成しているため、心筋組織の動態観測を詳細に観察することが可能となる。又、閉曲面に近い領域ほど低い不透明度を設定し、深部の領域に対しては高い不透明度を設定しているため、閉曲面が心筋組織の心内膜あるいは心外膜を正確に捉えていない場合においてもこの閉曲面の内部から切り出したボリュームデータに基づいて生成される３次元画像データにおいて心筋組織の表面を確実に捉えることが可能となる。
【００８６】
更に、心内膜輪郭面及び心外膜輪郭面において形成されたワイヤフレームのメッシュパターンを上述の心機能パラメータや３次元画像データに重畳させて動画像表示しているため、トラッキング時における心筋組織の運動状態を３次元画像データにおいて詳細に把握することができ、メッシュパターンと３次元画像データとの比較観察によりボリュームデータ間のトラッキング処理が正確に行なわれているか否かを容易に確認することができる。
【００８７】
以上、本発明の実施例（第１の実施例及び第２の実施例）について述べてきたが、本発明は、上述の実施例に限定されるものではなく、変形して実施することが可能である。例えば、上述の実施例では、心拍時相Ｐ１のボリュームデータに基づいて生成されたＭＰＲ画像データの心内膜及び心外膜に対して心内膜輪郭線及び心外膜輪郭線を設定し、これらの輪郭線を補間処理して心内膜輪郭面及び心外膜輪郭面を検出する場合について述べたが、前記ＭＰＲ画像データの心内膜及び心外膜に対して離散的な複数の輪郭点を設定し、これらの輪郭点を補間処理して心内膜輪郭線及び心外膜輪郭線と心内膜輪郭面及び心外膜輪郭面を検出してもよい。又、特開２００４−７３８５０号公報等に記載されたＡＣＴ（Automated Contour Tracking）法等を適用して心内膜輪郭線あるいは心内膜輪郭面を自動検出し、更に、心内膜輪郭面から所定間隔離れた位置に心外膜輪郭面を自動設定してもよい。
【００８８】
又、上述の実施例では、心拍時相Ｐ１のボリュームデータと心拍時相Ｐ２乃至Ｐｎのボリュームデータとのトラッキング処理によって各心拍時相における心内膜輪郭面及び心外膜輪郭面やこれらの輪郭面を含む閉曲面を検出する場合について述べたが、隣接する心拍時相のボリュームデータをトラッキング処理して心拍時相Ｐ２乃至Ｐｎにおける心内膜輪郭面及び心外膜輪郭面や閉曲面を検出してもよい。又、上述のトラッキング処理は、ボリュームデータ間のパターンマッチングによって行なう場合について述べたが他の方法を適用してもよく、特に限定されない。
【００８９】
更に、心機能パラメータ計測部９あるいは心機能パラメータ計測部９ａが追跡点の位置情報に基づいて計測する心機能パラメータは上述の「歪み（ストレイン）」に限定されるものではなく、例えば、「変位」、「回転」、「捻れ（torsion）」、「速度」等を心機能パラメータとして計測してもよく、これらの時間的変化を示す「歪みレート」、「回転レート」、「捻れレート」、「加速度」等であっても構わない。
【００９０】
一方、閉曲面内のボリュームデータに基づいて組織画像データ生成部８あるいは組織画像データ生成部８ａが生成する組織画像データは３次元画像データに限定されるものではなく投影画像データであっても構わない。
【００９１】
又、複数の振動素子が二次元配列された超音波プローブ３を用いてボリュームデータを収集する場合について述べたが、複数の振動素子が一次元配列された超音波プローブを機械的に移動あるいは回動させることによって上述のボリュームデータを収集してもよい。更に、整相加算後の受信信号を処理して得られたＢモードデータを用いてボリュームデータを生成する場合について述べたが、ボリュームデータは、カラードプラデータ等の他の超音波データに基づいて生成しても構わない。
【００９２】
又、上述の実施例における入力部１２に心内膜輪郭面及び心外膜輪郭面の位置や閉曲面の厚み（即ち、心内膜輪郭面と心外膜輪郭面の間隔）を更新する機能をもたせ、操作者は、表示部１１に表示された診断用画像データあるいはＭＰＲ画像データや３次元画像データの観察下にて入力部１２の入力デバイスを操作することにより輪郭面や閉曲面の最適化を行なってもよい。
【００９３】
更に、ワイヤフレームのメッシュパターンの替わりに心内膜輪郭面及び心外膜輪郭面の変位（移動方向と移動量）を示す移動ベクトルや移動成分ベクトルを組織画像データや心機能パラメータデータと共に表示しても構わない。
【００９４】
又、上述の実施例では、受信信号処理部４によって生成された３次元Ｂモードデータを補間処理してボリュームデータを生成する場合について述べたが、前記３次元ボリュームデータそのものをボリュームデータとして用いても構わない。
【図面の簡単な説明】
【００９５】
【図１】本発明の第１の実施例における３次元超音波診断装置の全体構成を示すブロック図。
【図２】同実施例の３次元超音波診断装置が備える送受信部及び受信信号処理部の具体的な構成を示すブロック図。
【図３】同実施例の３次元走査における超音波送受信方向を説明するための図。
【図４】同実施例の３次元超音波診断装置が備えるボリュームデータ生成部の具体的な構成を示すブロック図。
【図５】同実施例のトラッキング処理に用いられるパターンマッチングの具体例を示す図。
【図６】同実施例のワイヤフレーム形成部において形成されるワイヤフレームの具体例を示す図。
【図７】同実施例における診断用画像データの生成／表示手順を示すフローチャート。
【図８】同実施例の組織画像データ生成部が生成するＭＰＲ画像データの変形例を示す図。
【図９】本発明の第２の実施例における３次元超音波診断装置の全体構成を示すブロック図。
【図１０】同実施例の組織画像データ生成部において設定される不透明度特性を示す図。
【図１１】同実施例における診断用画像データの生成／表示手順を示すフローチャート。
【符号の説明】
【００９６】
２…送受信部
２１…送信部
２２…受信部
３…超音波プローブ
４…受信信号処理部
５…ボリュームデータ生成部
６、６ａ…輪郭面検出部
７、７ａ…ワイヤフレーム形成部
８、８ａ…組織画像データ生成部
９、９ａ…心機能パラメータ計測部
１０、１０ａ…診断用画像データ生成部
１１…表示部
１２…入力部
１２１…ＭＰＲ断面設定機能
１２２…輪郭設定機能
１３…走査制御部
１４…心拍時相計測ユニット
１５…システム制御部
１００、２００…３次元超音波診断装置
【図１】
2010115372.tif 000004
【図２】
2010115372.tif 000005
【図３】
2010115372.tif 000006
【図４】
2010115372.tif 000007
【図５】
2010115372.tif 000008
【図６】
2010115372.tif 000009
【図７】
2010115372.tif 000010
【図９】
2010115372.tif 000011
【図１０】
2010115372.tif 000012
【図１１】
2010115372.tif 000013
【図８】
2010115372.tif 000014