(57)【要約】【課題】  心臓壁筋の全体に関する運動情報を一枚の画像によって観察することができる超音波診断装置等を提供すること提供すること。【解決手段】  心臓に関するボリュームデータにおいて、視点を心腔内の中心付近に設定される軸の各点におき、心腔の基部から中央部付近においては、軸に対して垂直な平面上で視点を中心として視線を３６０°回転させるようにし、心尖付近においては、軸に対して垂直な平面上で視点を中心として視線を３６０°回転させると共に、所定の立体角だけ回転させるようにする。各視点、各視線を用いて、心臓内膜の全体を一覧可能な展開図上の対応する位置に、心臓組織の視点から視線方向に沿った輝度情報分布を一つの値に集約させる所定のマッピング処理を実行し、展開画像を生成し、そのまま表示あるいは展開画像を用いて生成された運動情報のマッピング画像を表示する。【選択図】      図２
2010094287.tif 000002
【特許請求の範囲】
【請求項１】
被検体内の心臓を含む領域に超音波を送信し、前記送信した超音波に基づく反射波を受信し、前記受信した反射波に用いて前記心臓を含む領域に関するボリュームデータを生成するデータ生成手段と、
前記ボリュームデータを用いて、心臓内膜を一覧可能な展開図上に、前記心臓の組織に関する輝度情報分布をマッピングする所定のマッピング処理を実行することで、前記心臓の形態に関する三次元情報を示す第１の展開画像を生成する画像生成手段と、
生成された前記第１の展開画像を表示する表示手段と、
を具備することを特徴とする超音波診断装置。
【請求項２】
前記画像生成手段は、
前記ボリュームデータの心腔に対応する領域に中心軸を設定し、
前記中心軸上に複数の視点を設定し、
前記ボリュームデータの心尖に対応する領域おいては、前記中心軸に対して垂直な平面上で視点を中心として視線を３６０°回転させるようにし、
前記ボリュームデータの心尖に対応しない領域おいては、前記中心軸に対して垂直な平面上で視点を中心として視線を３６０°回転させると共に、所定の立体角だけ回転させるようにして、前記所定のマッピング処理を実行すること、
を特徴とする請求項１記載の超音波診断装置。
【請求項３】
前記第１の展開図は、心尖を中心とし、長軸方向を半径方向とし、短軸方向を回転角方向とする極座標系展開図であることを特徴とする請求項１又は２記載の超音波診断装置。
【請求項４】
時相の異なる複数の前記第１の展開画像を用いて、前記心臓の組織に関する運動情報を算出する算出手段をさらに具備し、
前記画像生成手段は、前記運動情報を前記展開図上の対応する位置にマッピングするマッピング処理を実行することで、前記心臓の前記運動情報に関する三次元情報を示す第２の展開画像を生成し、
前記表示手段は、前記第２の展開画像を所定の形態で表示すること、
を特徴とする請求項１乃至３のうちいずれか一項記載の超音波診断装置。
【請求項５】
前記表示手段は、前記第１の展開画像と前記第２の展開画像とを同時に又は交互に表示することを特徴とする請求項４記載の超音波診断装置。
【請求項６】
前記画像生成手段は、前記第１の展開画像と前記第２の展開画像との間の差異を示す第３の展開画像を生成し、
前記表示手段は、前記第１の展開画像、前記第２の展開画像、前記第３の画像を同時に又は交互に表示することを特徴とする請求項４記載の超音波診断装置。
【請求項７】
前記表示手段は、収集された時刻の異なる複数の前記第１の展開画像、又は収集された時刻の異なる複数の前記第２の展開画像を同時に又は交互に表示することを特徴とする請求項４記載の超音波診断装置。
【請求項８】
前記画像生成手段は、収集された時刻の異なる複数の前記第１の展開画像を同時に又は交互に表示することを特徴とする請求項４記載の超音波診断装置。
【請求項９】
前記画像生成手段は、収集された時刻の異なる二つの前記第２の展開画像の間の差異を示す第４の展開画像を生成し、
前記表示手段は、前記二つの第２の展開画像と前記第４の展開画像とを同時に又は交互に表示することを特徴とする請求項４記載の超音波診断装置。
【請求項１０】
前記画像生成手段は、
前記ボリュームデータにおいて、前記心臓の組織境界に対応する領域を検出し、
前記所定のマッピング処理において用いられる投影面と前記組織境界に対応する領域との間に存在するボクセルを除去するノイズ除去処理を実行した後、前記所定のマッピング処理を実行すること、
を特徴とする請求項１乃至９のうちいずれか一項記載の超音波診断装置。
【請求項１１】
前記画像生成手段は、
前記投影面上の各ピクセルから当該投影面に垂直な方向に探索方向を設定し、
前記探索方向上に重み付き加算区間を設定し、
前記探索方向に沿って前記加算区間を逐次移動させ、各移動位置における前記加算区間についての重み付き加算値を取得し、
前記重み付き加算値を所定の閾値とに基づいて前記組織境界に対応する領域を検出し、
前記投影面と前記組織境界に対応する領域との間に存在するボクセルの値をゼロにする、或いは前記ボクセルの値に減弱係数を乗ずることで、前記ノイズ除去処理を実行すること、
を特徴とする請求項１０記載の超音波診断装置。
【請求項１２】
前記画像生成手段は、
前記投影面上の各ピクセルから当該投影面に垂直な方向に探索方向を設定し、
前記探索方向上のボクセルに対して、投影面に近い側のボクセルから前記探索方向に沿って累積加算値を計算し、
前記累積加算値と所定の閾値とに基づいて前記組織境界に対応する領域を検出し、
前記投影面と前記組織境界に対応する領域との間に存在するボクセルの値をゼロにする、或いは前記ボクセルの値に減弱係数を乗ずることで、前記ノイズ除去処理を実行すること、
を特徴とする請求項１０記載の超音波診断装置。
【請求項１３】
被検体の心臓を含む領域に超音波を送信して得られる反射波を用いて生成された、前記心臓を含む領域に関するボリュームデータを記憶する記憶手段と、
前記ボリュームデータを用いて、心臓内膜を一覧可能な展開図上に、前記心臓の組織に関する輝度情報分布をマッピングする所定のマッピング処理を実行することで、前記心臓の形態に関する三次元情報を示す展開画像を生成する画像生成手段と、
生成された前記展開画像を表示する表示手段と、
を具備することを特徴とする超音波信号処理装置。
【請求項１４】
コンピュータに、
被検体の心臓を含む領域に超音波を送信して得られる反射波に基づいて生成された、前記心臓を含む領域に関するボリュームデータを用いて、心臓内膜を一覧可能な展開図上に、前記心臓の組織に関する輝度情報分布をマッピングする所定のマッピング処理を実行することで、前記心臓の形態に関する三次元情報を示す展開画像を生成させる画像生成機能と、
生成された前記展開画像を表示させる表示機能と、
を実現させることを超音波信号処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、組織の運動に関する時間変化を映像化する機能を有する超音波診断装置等に関する。
【背景技術】
【０００２】
超音波診断装置は、被検体内の所定領域に超音波を送信し、当該被検体内から反射・散乱されてくる反射超音波に基づいて組織の構造や血流情報をリアルタイムに映像化する医療用映像装置である。この装置の臨床的な目的の一つに、心臓壁筋の運動（以下、壁運動と称す）の様子を観察したり、壁運動の動きの程度を定量的なパラメータに置き換えて、それを画像として映像化したり、計測したりすることがある。この目的を達成するための手法は、例えば次のように分類することができる。
【０００３】
まず、壁運動形状そのものの形状を映像化し観察するための手法がある。この手法は、例えば組織ドプラ法に従う三次元スキャンにより得られた三次元データ（ボリュームデータ）を、ボリュームレンダリング（VR）法やサーフェイスレンダリング（SR）法等の三次元指表示手法により映像化し、表示するものである。
【０００４】
また、壁運動情報を定量的に観察するための手法がある。これは、いくつかの２次元断面に対して、壁運動情報に関する定量的な画像を得てから、三次元画像を構築するものである。なお、壁運動情報に関する定量的な画像を取得する代表的な手法としては、次の二つの手法がある。一つは、組織ドプラによる手法である。これは、組織ドプラ法により組織の移動速度を得て、上記壁運動の観察と同様にＶＲ等により画像として表示したり、それに対して、処理を行い他のパラメータ（Strain、Strain-Rate、Displacement等）を得てから、同様に画像として表示する。もう一つは、画像のパターンマッチングによる手法である。この手法では、画像のパターンマッチングにより組織の移動量を求め、これを用いて組織の移動速度を得て、上記壁運動の観察と同様にＶＲ法等により画像として表示するものである。
【発明の開示】
【発明が解決しようとする課題】
【０００５】
ところで、心筋は心臓の内部の組織であり、臨床的には特に心臓内膜側の様子が重要である。しかしながら、従来の手法において用いられるＶＲ法等の画像再構成法では、十分に壁運動を観察することができない場合がある。これは、ＶＲ法等は、視点から視線に沿って見た場合の物体（心臓）を、実際に見える様子にできるだけ近づけようとする手法であり、心臓内部の三次元的映像化に必ずしも適さないからである。
【０００６】
また、観察効率性の観点から、心筋の運動情報は、一枚の画像で十分に映像化されていることが好ましい。しかしながら、ＶＲ法等の画像再構成法では、視線から見える部分の影になってしまう部分があり、全体の運動情報を十分に映像化することができない。
【０００７】
また、従来の手法により壁運動を定量化して観察する場合、記述の如く組織ドプラを用いる方法又はパターンマッチングを用いる方法を用いる。しかしながら、リアルタイム三次元像に対して、これらの処理を実行する場合、情報量が多いために多大な時間を必要とし、臨床上実用的でない（この問題は、特にパターンマッチングを用いる方法において顕著である）。
【０００８】
本発明は、上記事情を鑑みてなされたもので、心臓壁筋の全体に関する運動情報を一枚の画像によって観察することができる超音波診断装置、超音波信号処理装置及び超音波信号処理プログラムを提供することを目的としている。
【課題を解決するための手段】
【０００９】
本発明は、上記目的を達成するため、次のような手段を講じている。
【００１０】
請求項１に記載の発明は、被検体内の心臓を含む領域に超音波を送信し、前記送信した超音波に基づく反射波を受信し、前記受信した反射波に用いて前記心臓を含む領域に関するボリュームデータを生成するデータ生成手段と、前記ボリュームデータを用いて、心臓内膜を一覧可能な展開図上に、前記心臓の組織に関する輝度情報分布をマッピングする所定のマッピング処理を実行することで、前記心臓の形態に関する三次元情報を示す第１の展開画像を生成する画像生成手段と、生成された前記第１の展開画像を表示する表示手段と、を具備することを特徴とする超音波診断装置である。
【００１１】
請求項１３に記載の発明は、被検体の心臓を含む領域に超音波を送信して得られる反射波を用いて生成された、前記心臓を含む領域に関するボリュームデータを記憶する記憶手段と、前記ボリュームデータを用いて、心臓内膜を一覧可能な展開図上に、前記心臓の組織に関する輝度情報分布をマッピングする所定のマッピング処理を実行することで、前記心臓の形態に関する三次元情報を示す展開画像を生成する画像生成手段と、生成された前記展開画像を表示する表示手段と、を具備することを特徴とする超音波信号処理装置である。
【００１２】
請求項１４に記載の発明は、コンピュータに、被検体の心臓を含む領域に超音波を送信して得られる反射波に基づいて生成された、前記心臓を含む領域に関するボリュームデータを用いて、心臓内膜を一覧可能な展開図上に、前記心臓の組織に関する輝度情報分布をマッピングする所定のマッピング処理を実行することで、前記心臓の形態に関する三次元情報を示す展開画像を生成させる画像生成機能と、生成された前記展開画像を表示させる表示機能と、を実現させることを超音波信号処理プログラムである。
【発明の効果】
【００１３】
以上本発明によれば、心臓壁筋の全体に関する運動情報を一枚の画像によって観察することができる超音波診断装置、超音波信号処理装置及び超音波信号処理プログラムを実現することができる。
【発明を実施するための最良の形態】
【００１４】
以下、本発明の第１実施形態及び第２実施形態を図面に従って説明する。なお、以下の説明において、略同一の機能及び構成を有する構成要素については、同一符号を付し、重複説明は必要な場合にのみ行う。
【００１５】
（第１実施形態）
図１は、本実施形態に係る超音波診断装置１の構成を示したブロック図である。同図に示すように、本超音波診断装置１は、超音波プローブ１１、送信回路１３、受信回路１５、ビームフォーマ１７、Ｂ／Ｍ処理部１９、ＣＦＭ処理部２１、画像生成部２３、ビデオＩ／Ｆ２７、表示部２９、制御部３０、記憶部３１、入力ユニット３３を具備している。
【００１６】
超音波プローブ１１は、送信回路１３からの駆動信号（駆動パルス）に基づき超音波を発生し且つ被検体からの反射波を電気信号（エコー信号）に変換する複数の圧電振動子、当該圧電振動子に設けられる整合層、当該圧電振動子から後方への超音波の伝播を防止するバッキング材等を有している。当該超音波プローブ１１から被検体に超音波が送信されると、送信超音波は体内組織の音響インピーダンスの境界、微小散乱等により後方散乱され、反射波として超音波プローブ１１に受信される。
【００１７】
送信回路１３は、送信超音波を形成するため、所定のレート周波数ｆｒ  Ｈｚ（周期；１／ｆｒ秒）でレートパルスを繰り返し発生し、チャンネル毎に超音波をビーム状に集束し且つ送信指向性を決定するのに必要な遅延時間を、各レートパルスに与える。また、送信回路１３は、必要な遅延時間が与えられたレートパルスに従って、プローブ１１に対してチャンネル毎に駆動電圧パルスを印加する。
【００１８】
受信回路１５は、超音波プローブ１１が受信したエコー信号をチャンネル毎に増幅し、ビームフォーマ１７に増幅後のエコー信号を出力する。
【００１９】
ビームフォーマ１７は、増幅されたエコー信号をＡ／Ｄ変換し、受信指向性を決定するのに必要な遅延時間を与え、整相加算によるフォーカシングを行う。このフォーカシングにより、エコー信号の受信指向性に応じた方向からの反射成分が強調され、受信指向性と送信指向性とにより超音波送受信の総合的な指向性（走査線）が決定される。
【００２０】
Ｂ／Ｍ処理部１９は、所定の周波数帯域の反射波成分を減衰させ、所望の周波数帯域の反射波成分を抽出するバンドパスフィルタ処理を行った後、包絡線検波を行うことでその包絡線成分を検出する。その後、Ｂ／Ｍ処理部１９は、エコー信号に対して対数圧縮処理及び必要に応じてエッジ強調処理等を実行する。
【００２１】
ＣＦＭ処理部２１は、エコー信号を組織信号と血流信号とに分離するためのハイパスフィルタ処理（「ＭＴＩフィルタ処理」、又は「Ｗａｌｌフィルタ処理」とも呼ばれる。）、組織や血流の移動速度を検出するための自己相関処理を実行する。また、ＣＦＭ処理部２１は、必要に応じてエコー信号中の組織信号を低減・削除するための非線形処理を実行する。
【００２２】
画像生成部２３は、Ｂ／Ｍ処理部１９、ＣＦＭ処理部２１から入力した信号列を、超音波ビームの送受信に対応した位置にマッピングし、画像化する。また、画像生成部２３は、制御部３０の制御のもと、後述する展開画像生成機能を用いた処理（展開画像生成処理）を実行する。
【００２３】
ビデオＩ／Ｆ２７は、Ｂモード／Ｍモード用の画像信号、ＣＦＭからの画像信号、これらの三次元画像、その他の超音波画像に関する種々の情報を各画像と組み合わせ、モニタ表示用にレイアウトされた表示画像を構成し、表示部２９に当該表示画像を表示する。
【００２４】
表示部２９は、画像生成部２３からのビデオ信号に基づいて、生体内の形態学的譲歩や、血流情報を画像として表示する。また、造影剤を用いた場合には、造影剤の空間的分布、つまり血流或いは血液の存在している領域を求めた定量的な情報量に基づいて、輝度画像やカラー画像を表示する。
【００２５】
制御部３０は、ＣＰＵ３０１、メモリ３０３を有しており、システム全体の制御中枢として、本超音波診断装置１を静的又は動的に制御する。当該制御部３０の制御のもと、データ（超音波信号）は、各ユニットで種々の処理を受け次段のユニットに送り出されるとともに、メモリ３０３を介して又は直接記憶部３１等に保存される。なお、すべての段階での超音波信号を記憶部３１に記憶する必要はなく、必要に応じて選択的に記憶してもよい。また、各ユニット内部で処理中の信号が保存可能な構成であってもよい。
【００２６】
また、制御部３０は、後述する展開画像生成処理に関する制御を実行する。より具体的には、記憶部３１に記憶されているプログラムを読み出し、これに基づいてメモリ３０３にシステムを展開することで、本移動ベクトル生成処理を実行する。なお、これに限定されず、ハードウェア構成によって本移動ベクトル生成処理を実行してもよい。
【００２７】
記憶部３１は、過去に撮影したＢモード画像、組織ドプラ画像等をフレーム毎に記憶している。また、記憶部３１は、リアルタイムで撮影されるＢモード画像、組織ドプラ画像、後述する運動情報に関する画像等は、メモリ３０３に一次記憶されると共に、所定の操作に応答して、記憶部３１に記憶される。
【００２８】
入力部３３は、装置１の本体に接続され、オペレータからの各種指示・命令・情報を装置本体にとりこむための入力装置（マウスやトラックボール、モード切替スイッチ、キーボード等）が設けられる。
【００２９】
（展開画像生成機能）
次に、本超音波診断装置１が有する、展開画像生成機能について説明する。この機能は、ボリュームスキャンによって得られたボリュームデータを用いて、所定の投影処理によるマッピングを行うことにより、心筋内膜に関する展開画像を生成するものである。
【００３０】
図２は、本展開画像生成処理の流れを示したフローチャートである。同図に従って、展開画像生成処理において実行される各処理について説明する。
【００３１】
［心臓組織の抽出／軸の設定：ステップＳ１］
まず、制御部３０は、記憶部３１に記憶された所定期間内の各時相のボリュームデータを読み出す。画像生成部２７は、各ボリュームデータにおいて所定の手法（例えば、閾値処理を用いたセグメンテーション等）により、心臓組織領域を抽出する。また、画像生成部２７は、抽出した心臓組織領域の心腔に中心軸を設定する。この中心軸の設定は、例えば心腔を楕円等で近似しその軸を採用することで、設定することができる。
【００３２】
なお、読み出した（或いはスキャンによって得られた）ボリュームデータが直交系ボリュームデータでない場合には、本ステップにおいて所定の座標変換し直交系ボリュームデータを作成すればよい。
【００３３】
［心臓部位の判定：ステップＳ２］
次に、画像生成部２７は、心臓部位を心尖領域及び非心尖領域（すなわち、心腔基部から心腔中部付近にかけての領域）とを判定する。この判定は、例えば設定した中心軸の長さ、或いは心腔容積と中心軸の長さとの比率等に基づいて、判定することができる。
【００３４】
［時相毎の展開画像の生成：ステップＳ３］
次に、画像生成部２７は、各ボリュームデータに対して、設定した中心軸上に複数の視点をおき、各視点から（心臓の）外側に向かう視線方向により、心臓内膜の全体を一覧可能な展開図上の対応する位置に、視点から視線方向に沿った心臓組織の輝度情報分布を一つの値に集約させる所定のマッピング処理（例えば、ボリュームレンダリング処理）を実行し、形態に関する展開画像を時相毎に生成する。
【００３５】
図３は、本ステップにおける視点、視線方向の設定、マッピング処理を説明するための図である。同図に示すように、画像生成部２７は、中心軸Ａの各点に視線方向を設定し、非心尖領域では、視点を含み中心軸に垂直な平面Ｐ上において視線を３６０度回転するようにして、マッピング処理を行う。また、画像生成部２７は、心尖領域では、視点を含み中心軸に垂直な平面Ｐ上において視線を３６０度回転するのに加えて、軸方向への角度（立体角θ）を振っていくようにして、マッピング処理を行う。なお、当然ながら、立体角θは、全ての心臓組織が展開図上にマッピングされるように、各視点毎に所定の値が予め設定される。本ステップにおいて生成される形態に関する展開画像は、表示部２９において所定の形態で表示することができる。
【００３６】
本ステップにおけるマッピング処理は、視点から視線方向に沿った心臓組織の輝度情報分布を一つの値に集約させるものであれば、特に限定はない。典型例としては、ボリュームレンダリング処理（或いは、これと同様のアルゴリズムをもつ処理）、サーフェスレンダリング処理、ＭＩＰ処理、最初に心壁と見なすことができるボクセルの値を閾値により抽出し代表させる処理等を挙げることができる。
【００３７】
また、本ステップにおいて用いられる展開図は、心臓内膜の全体を一覧可能なものであれば、特に限定はない。図４は、心尖部をその中心とし、当該心尖部から心腔基部に向かう深さ方向（すなわち、長軸方向）を半径方向とし、円周方向は水平面上の角度とした極座標系に対応する展開図の例を示している。この例は、いわゆるブルズアイ（Bull’s Eye）表示と言われるものであり、Ｂモードと類似の階調表現をもつ動画像として、基部から心尖部に向かう心筋の運動を一覧可能に表現することができる。また、図４の例の他、図５に示すように、心尖部を頂点とし当該心尖部から心腔基部に向かう方向に沿って心臓を切り開いた展開図等を採用することも可能である。
【００３８】
［展開画像を用いた移動ベクトル／移動速度の計算：ステップＳ４］
次に、画像生成部２７は、隣り合う時相の展開画像を用いて相互相関処理を実行し、各位置の各時相における移動ベクトルを計算する。また、画像生成部２７は、計算された各位置の各時相における移動ベクトルと隣り合う時相の時間間隔とに基づいて、各位置の各時相における移動速度を計算する。
【００３９】
［運動情報の演算：ステップＳ５］
次に、画像生成部２７は、必要に応じて、計算された移動ベクトル、移動速度を用いて運動情報を計算する。計算する運動情報には、特に拘泥されない。典型例としては、例えば壁厚方向への変化に関する運動情報（Radial-strainやRadial-strain rate）、長軸方向への変化に関する運動情報（Longitudinal-strainやLongitudinal-strain rate）、円周方向への変化に関する運動情報（Circumferential-strainやCircumferential-strain rate）、短軸面内での面積重心に関する運動情報（RotationやRotation rate）、異なる短軸面間の回転の差分である運動情報（TwistやTwist rate）、短軸面間の距離でTwist情報を規格化した運動情報（TorsionやTorsion rate）、移動距離に関する運動情報（DisplacementやVelocity ）等を挙げることができる。
【００４０】
なお、本実施形態においては、説明を具体的にするため、例えば運動情報としてLongitudinal-strainが計算されてたものとする。
【００４１】
［運動情報のマッピング：ステップＳ６］
次に、画像生成部２７は、心臓内膜の全体を一覧可能な展開図上の対応する位置に、計算した運動情報をカラーコード化しマッピングするマッピング処理を実行し、運動情報に関する展開画像を時相毎に生成する。
【００４２】
［マッピング画像の表示：ステップＳ７］
次に、表示部２９は、運動情報に関する展開画像を所定の形態にて表示する。
【００４３】
図６は、図４と同様の展開図を用いてマッピングされた運動情報に関する展開画像の一例を示した図である。また、図７は、図５と同様の展開図を用いてマッピングされた運動情報に関する展開画像の一例を示した図である。この運動情報に関する展開画像においては、壁運動異常のある部位は、歪みの小さな部位として表現される。なお、必要に応じて、歪み像のPeakタイミングを画像化しても良い。また、いずれの例においても、運動情報の強弱と色彩との対応を示すカラーバーを同時に表示することが好ましい。
【００４４】
運動情報に関する展開画像の表示形態は、図６、図７の例に拘泥されない。表示形態の変形例としては、次のようなものを挙げることができる。
【００４５】
図８に示すように、形態に関する展開画像（造影剤使用の有無に拘泥されない）及び運動情報に関する展開画像を同時に（又は交互に）表示してもよい。また、係る場合には、形態に関する展開画像と運動情報（例えばStrain）に関する展開画像との間で異なる部分を抽出し、図９に示すように当該部分のみがカラー表示された展開画像を生成し同時に（又は交互に）表示するようにしてもよい。
【００４６】
また、図１０に示すように、時間的に隔たりのある（収集された時刻の異なる）複数の運動情報に関する展開画像を同時に（又は交互に）表示するようにしてもよい。この場合の典型例としては、治療前と治療後の各展開画像、ストレスエコーにおけるストレス前とストレス後の各展開画像を同時に表示する場合等を挙げることができる。さらに、時間的に隔たりのある複数の（例えば二枚）運動情報に関する展開画像を同時に表示する場合には、図１１に示すように、双方の差分画像を生成し同時に（又は交互に）表示するようにしてもよい。
【００４７】
い。なお、図１０、図１１に示す表示形態は、時間的に隔たりのある複数の形態に関する展開画像を同時に（又は交互に）表示する場合にも、用いることができる。
【００４８】
（効果）
以上述べた構成によれば、以下の効果を得ることができる。
【００４９】
本超音波診断装置においては、視点を心腔内の中心付近に設定される軸の各点におき、心腔の基部から中央部付近においては、軸に対して垂直な平面上で視点を中心として視線を３６０°回転させるようにし、心尖付近においては、軸に対して垂直な平面上で視点を中心として視線を３６０°回転させると共に、所定の立体角だけ回転させるようにする。この様に設定された各視点、各視線を用いて、心臓内膜の全体を一覧可能な展開図上の対応する位置に、心臓組織の視点から視線方向に沿った輝度情報分布を一つの値に集約させる所定のマッピング処理（例えば、ボリュームレンダリング処理）を実行し、展開画像を生成し、そのまま表示あるいは展開画像を用いて生成された運動情報のマッピング画像を表示する。従って、操作者は、生成された展開画像等を観察することで、心臓組織全体についての形態や運動情報等を一覧可能な新たな医療情報を把握することができる。また、心臓組織全体についての形態や運動情報等が一覧可能であるため、画像の切り換えや複数画像の同時観察を必要とせず、画像診断時における作業負担を軽減することができる。
【００５０】
また、本超音波診断装置によれば、形態に関する展開画像及び運動情報に関する展開画像、時間的に隔たりのある複数の運動情報に関する展開画像、及びこれらを用いて生成された差分画像等を同時に（又は交互に）表示することができる。従って、操作者は、表示された展開画像を観察することで、心臓組織全体についての形態や運動情報等を一覧可能な新たな医療情報を把握することができる。特に、差分画像によれば、時間経過に伴い変化の生じた部位等を容易に把握することができ、画像診断時における作業負担を軽減することができる。
【００５１】
（第２の実施形態）
次に、本発明の第２の実施形態について説明する。本実施形態に係る超音波診断装置１では、映像化対象以外の領域に対応するボクセルをできる限り取り除くことで、ノイズの少ない好適な診断画像を提供するものである。なお、本実施形態に係る手法の適用は、生成しようとする画像が展開画像であるか否かに拘泥されない。しかしながら、本実施形態では、説明を具体的にするため、第１の実施形態において説明した展開画像を生成する場合への適用を例として説明する。
【００５２】
一般に、ボリュームデータを用いて三次元情報を映像化する場合、ボクセルに含まれるノイズ（スペックル・ノイズやＣＴ画像のフォトン・ノイズ等）は画像の画質を劣化させる。特に、三次元画像生成のための投影面とボクセル内の物体表面との間に存在するボクセルに含まれるノイズ成分が画質に大きく影響する。
【００５３】
図１２は、一般的なボリュームレンダリング処理の形態を説明するための図である。図１３は、図１２内の投影面に直交する１断面上での投影処理を表した模式図である。各図に示すように、投影面と映像化対象である物体と間に存在するボクセルに含まれるノイズを除去するには、Cropping（あるいはClipping）処理により投影面とボクセル内の物体表面との間に存在する不要なボクセル・データを切削したのち、レンダリング処理などを行うのが一般的である。
【００５４】
図１４、図１５は、Cropping処理を説明するための概念図である。Cropping面は平面で設定されるのが一般的であり、図１４に示す例では投影面に対して平行に設定されている。Cropping処理を行う場合に、図１５に示すようにCropping面を三次元空間内に曲面として設定できれば三次元画像の画質には好ましい。しかしながら、一般に曲面を三次元空間内に設定することは操作性、実現性の両面で困難である。また、Cropping面が図１４のように平面であったとしても、平面を三次元空間内で自由に操作することは、操作者にとって煩雑な操作となる。
【００５５】
そこで、本実施形態に係る超音波診断装置１は、物体表面を基準としたノイズ除去機能を有する。この機能は、投影面上の各ピクセルから投影面に垂直な方向に探索を行い、ボクセル内の物体（組織）の表面位置を検出し、不要なボクセル（投影面と物体表面の間のボクセル）を削除した後、マッピング処理を行うものである。なお、物体表面を基準としたノイズ除去機能に従う処理（以下、単にノイズ除去処理）は、画像生成部２７において、例えば図２に示したステップＳ１の心壁抽出処理において実行される。
【００５６】
まず、画像生成部２７は、図１６に示すように、投影面の各ピクセルから投影面に垂直な方向に探索方向を設定する。さらにこの探索方向上に図１７に示すように、重み付き加算区間を設定し、加算区間を探索方向に沿って移動させ、それぞれの位置において、例えば次の様な重み付き加算処理を行う。
【００５７】
Ａ（ｚ）  ＝  Σn｛ｆ(n)×P（ｚ＋ｎ）｝
ここで、Ａは加算値、ｆ(n)は重み係数、P（ｚ＋ｎ）は画素値、ｚは投影面に対する垂線方向の位置をそれぞれ表す。また、ｆ(n)  ≧０，Σn｛ｆ(n)｝＝１を満たす。
【００５８】
また、画像生成部２７は、加算区間を移動させながら重み付き加算処理結果Ａを予め装置にプリセットされている閾値Ｔ（Ｔ＞０）と比較を行う。そして図１８のように、前記加算結果Ａが前記閾値Ｔを越えた加算区間の中間点を物体の表面位置とする。なお、物体の表面位置の決定手法は、当該中間点の例じ拘泥されない。例えば、加算区間の中間点以外に、区間の端点を表面位置とするようにしても良い。
【００５９】
このように一定区間に渡る重み付き加算処理を行うと、図１９のようなノイズに対して加算結果Ａはあまり変動せず、図２０のような組織の境界においてのみ加算結果Ａは大きく変動する。従って、閾値Ｔを適切に設定した加算処理により、所望の組織境界を検出することができる。なお、閾値Ｔの値は、映像化対象に応じた値が初期設定されている。しかしながら、必要に応じて、入力装置３３からの操作により、任意の値を設定するようにしてもよい。
【００６０】
このように投影面上のピクセルから投影面に垂直な方向に探索を行い、三次元画像を生成する物体表面（心臓組織境界）を抽出し、図２１に示すように探索の始点と組織の表面位置の間のボクセルのボクセル値をゼロにした後、図２のステップＳ３以降のレンダリング処理等を行う。あるいは、ボクセル値をゼロにする代わりに、当該のボクセルに対して減弱係数を乗じてもよい。
【００６１】
（変形例）
次に、本実施形態の変形例について説明する。物体表面の検出方法は、上記重み付け加算処理に限定されず、例えば次の様な手法を用いても実現することができる。
【００６２】
図２２は、本変形例に係る物体表面の検出方法を説明するための図である。同図に示すように、まず、画像生成部２７は、投影面上の各ピクセルから投影面に対して垂直な方向に探索方向（探索線）を設定する。次に、画像生成部２７は、探索線上に存在するボクセルに対して投影面に近い側のボクセルからボクセル値の累積加算処理を行う。
【００６３】
図２３は、累積加算処理によって得られるある探索線に関する累積加算値の視点からの距離に対する変化を示した図である。画像生成部２７は、累積加算値を所定の閾値Ｔと比較し、図２３に示すような累積加算値が閾値Ｔを越えたボクセルを、物体表面（組織境界）と判定する。
【００６４】
次に、画像生成部２７は、組織表面位置が検出された後は、図２４に示すように第１の方法と同様に探索の始点と組織の表面位置の間のボクセルのボクセル値をゼロにした後、図２のステップＳ３以降のレンダリング処理等を行う。あるいは、ボクセル値をゼロにする代わりに、当該のボクセルに対して減弱係数を乗じてもよい。
【００６５】
（効果）
以上述べた構成によれば、第１の実施形態で述べた効果に加えて、次の効果を取得することができる。
【００６６】
すなわち、本超音波診断装置では、投影面上の各ピクセルから投影面に垂直な方向に投影を行い、ボクセル内の物体表面（組織境界）の位置を検出し、心壁の映像化に不要なボクセルを削除した後にマッピング処理を行う。従って、特別なCropping面の設定操作を行わずとも、映像化に不要なボクセルに含まれるノイズを簡単且つ高精度で取り除くことができ、画像の高い診断画像を提供することができる。
【００６７】
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。具体的な変形例としては、例えば次のようなものがある。
【００６８】
本実施形態に係る各機能は、当該処理を実行するプログラムを超音波画像処理装置等のワークステーションのコンピュータにインストールし、これらをメモリ上で展開することによっても実現することができる。このとき、コンピュータに当該手法を実行させることのできるプログラムは、磁気ディスク（フロッピー（登録商標）ディスク、ハードディスクなど）、光ディスク（ＣＤ−ＲＯＭ、ＤＶＤなど）、半導体メモリなどの記録媒体に格納して頒布することも可能である。
【００６９】
また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【産業上の利用可能性】
【００７０】
以上本発明によれば、心臓壁筋の全体に関する運動情報を一枚の画像によって観察することができる超音波診断装置、超音波信号処理装置及び超音波信号処理プログラムを実現することができる。
【図面の簡単な説明】
【００７１】
【図１】図１は、第１の実施形態に係る超音波診断装置１の構成を示したブロック図である。
【図２】図２は、本展開画像生成処理の流れを示したフローチャートである。
【図３】図３は、展開画像精製処理における視点、視線方向の設定、マッピング処理を説明するための図である。
【図４】図４は、心臓内膜の全体を一覧可能な展開図の一例を示した図である。
【図５】図５は、心臓内膜の全体を一覧可能な展開図の他の例を示した図である。
【図６】図６は、図４と同様の展開図を用いてマッピングされた運動情報に関する展開画像の一例を示した図である。
【図７】図７は、図５と同様の展開図を用いてマッピングされた運動情報に関する展開画像の一例を示した図である。
【図８】図８は、形態に関する展開画像及び運動情報に関する展開画像を同時に表示する場合の表示形態の一例を示した図である。
【図９】図９は、形態に関する展開画像、運動情報に関する展開画像、さらに双方の間で異なる部分が抽出されカラー表示された展開画像を同時に表示する場合の表示形態の一例を示した図である。
【図１０】図１０は、時間的に隔たりのある運動情報に関する展開画像（二つ）を同時に表示する場合の表示形態の一例を示した図である。
【図１１】図１１は、時間的に隔たりのある運動情報に関する展開画像（二つ）と差分画像とを同時に表示する場合の表示形態の一例を示した図である。
【図１２】図１２は、一般的なボリュームレンダリング処理の形態を説明するための図である。
【図１３】図１３は、図１２内の投影面に直交する１断面上での投影処理を表した模式図である。
【図１４】図１４は、Cropping処理を説明するための概念図である。
【図１５】図１５は、Cropping処理を説明するための概念図である。
【図１６】図１６は、第２の実施形態に係るノイズ除去処理を説明するための図である。
【図１７】図１７は、第２の実施形態に係るノイズ除去処理を説明するための図である。
【図１８】図１８は、第２の実施形態に係るノイズ除去処理を説明するための図である。
【図１９】図１９は、第２の実施形態に係るノイズ除去処理を説明するための図である。
【図２０】図２０は、第２の実施形態に係るノイズ除去処理を説明するための図である。
【図２１】図２１は、第２の実施形態に係るノイズ除去処理を説明するための図である。
【図２２】図２２は、第２の実施形態に係る物体表面の検出処理の他の例を説明するための図である。
【図２３】図２３は、第２の実施形態に係る物体表面の検出処理の他の例を説明するための図である。
【図２４】図２４は、第２の実施形態に係るレンダリング処理を説明するための図である。
【符号の説明】
【００７２】
１…超音波診断装置、１１…超音波プローブ、１３…送信回路、１５…受信回路、１７…ビームフォーマ、１９…Ｂ／Ｍ処理部、２１…ＣＦＭ処理部、２３…画像生成部、２７…ビデオＩ／Ｆ、２９…表示部、３０…制御部、３１…記憶部、３３…入力ユニット
【図１】
2010094287.tif 000003
【図２】
2010094287.tif 000004
【図３】
2010094287.tif 000005
【図４】
2010094287.tif 000006
【図５】
2010094287.tif 000007
【図６】
2010094287.tif 000008
【図７】
2010094287.tif 000009
【図８】
2010094287.tif 000010
【図９】
2010094287.tif 000011
【図１０】
2010094287.tif 000012
【図１１】
2010094287.tif 000013
【図１２】
2010094287.tif 000014
【図１３】
2010094287.tif 000015
【図１４】
2010094287.tif 000016
【図１５】
2010094287.tif 000017
【図１６】
2010094287.tif 000018
【図１７】
2010094287.tif 000019
【図１８】
2010094287.tif 000020
【図１９】
2010094287.tif 000021
【図２０】
2010094287.tif 000022
【図２１】
2010094287.tif 000023
【図２２】
2010094287.tif 000024
【図２３】
2010094287.tif 000025
【図２４】
2010094287.tif 000026