JP5134787B2 - 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム - Google Patents

Description

組織歪みイメージング（TSI: Tissue Strain Imaging）法を用いた心臓の映像化に関し、特に、心臓の収縮期における長軸方向の縮み（Shortening）情報の短軸画像へマッピング、又は心臓の心壁厚方向の伸張（Thickening）情報の長軸画像へのマッピングに関する。

一般に、心筋等の生体組織に関して、その運動や機能を客観的かつ定量的に評価することは、その組織の診断にとって非常に重要である。超音波画像処理装置を使用した画像診断においても、主に心臓を例として様々な定量的評価法が試みられている。例えば、正常な心筋は収縮期に壁厚方向（短軸）へは厚みが増し（thickening）、長軸方向へは長さが縮む（shortening）ことがわかっている。一般に、このthickeningとshorteningは、その運動方向が互いに直交し異なる起序を呈すると言われている一方、これらの運動を観察し心筋壁運動を評価することで、例えば心筋梗塞等の心臓疾患に関する診断支援の可能性が示唆されている。

なお、本願に関連する公知文献としては、例えば次のようなものがある。
特開平２００３−１７５０４１号

しかしながら、従来の超音波診断装置においては、例えば次のような問題がある。すなわち、短軸断層像を撮影し観察している場合には、超音波診断装置の特性上当該短軸断層像の走査面に対して垂直な生体情報を取得することはできないため、shortening情報を同時に表現することはできない。また、例えば組織ドプラ法によって例えば心尖長軸断層像を撮影し観察している場合には、ドプラ限界角が存在するため、thickening情報を同時に表現することはできない。

本発明は、上記事情を鑑みてなされたもので、２次元断層画像を表示する場合において、断層画像と直交する方向に関する運動情報を、当該断層画像にマッピングして画像化することができ、より具体的な例としては、短軸断層像を撮影し観察している場合であってもshortening情報を同時に観察することができ、また、長軸断層像を撮影し観察している場合であってもthickening情報を同時に観察することができる超音波診断装置、超音波画像処理装置及び超音波画像処理プログラムを提供することを目的としている。

本発明は、上記目的を達成するため、次のような手段を講じている。

一実施形態に従う超音波診断装置は、３次元空間における、心組織の第１の方向及び第２の方向に関する運動情報を求める運動情報生成ユニットと、超音波送受信により取得された前記３次元空間に関するボリュームデータに基づいて、前記３次元空間における第１の方向、及び前記第１の方向に対し直交する第２の方向に関する２次元画像を生成する２次元画像生成ユニットと、前記第１の方向に関する２次元画像に対して前記第２の方向に関する運動情報をマッピングした第１のマッピング画像と、前記第２の方向に関する２次元画像に対して前記第１の方向に関する運動情報をマッピングした第２のマッピング画像とを生成するマッピング画像生成ユニットと、前記第１あるいは第２のマッピング画像を表示する表示ユニットと、を具備することを特徴とする超音波診断装置である。

以上本発明によれば、短軸断層像を撮影し観察している場合であってもshortening情報を同時に観察することができ、また、長軸断層像を撮影し観察している場合であってもthickening情報を同時に観察することができる超音波診断装置、超音波画像処理装置及び超音波画像処理プログラムを実現することができる。

以下、本発明の第１実施形態乃至第３実施形態を図面に従って説明する。以下の説明において、略同一の機能及び構成を有する構成要素については、同一符号を付し、重複説明は必要な場合にのみ行う。

なお、以下の各実施形態においては、本発明の技術的思想を超音波診断装置に適用する場合を例として説明する。しかしながら、これに拘泥されることなく、本発明の技術的思想は、ワークステーション、パーソナルコンピュータ等の超音波画像処理装置についても適用可能である。

また、本実施形態に係る各構成要素、特に後述するＴＳＩ処理ユニット２０、ボリュームデータ生成ユニット２１、マッピング処理ユニット２４、追跡処理ユニット２７（図１、図１４参照）については、当該各構成要素と同様の処理を実行するソフトウェアプログラムをワークステーション等のコンピュータ、コンピュータ機能を有する超音波診断装置等にインストールし、これらをメモリ上で展開することによっても実現することができる。このとき、コンピュータに当該手法を実行させることのできるプログラムは、磁気ディスク（フロッピー（登録商標）ディスク、ハードディスクなど）、光ディスク（ＣＤ−ＲＯＭ、ＤＶＤなど）、半導体メモリなどの記録媒体に格納して頒布することも可能である。

（第１実施形態）
図１は、第１の実施形態に係る超音波診断装置１０の構成図である。本超音波診断装置１０は、超音波プローブ１１、送信ユニット１２、受信ユニット１３、Ｂモード処理ユニット１４、組織ドプラ処理ユニット１５、表示制御ユニット１７、表示ユニット１８、ＴＳＩ処理ユニット２０、ボリュームデータ生成ユニット２１、記憶部２２、制御ユニット（ＣＰＵ）２３、マッピング処理ユニット２４、入力ユニット２５を具備している。なお、本発明を超音波処理装置に適用する場合には、図１の点線内がその構成要素となる。

超音波プローブ１１は、送信ユニット１２からの駆動信号に基づき超音波を発生し、被検体からの反射波を電気信号に変換する複数の圧電振動子、当該圧電振動子に設けられる整合層、当該圧電振動子から後方への超音波の伝播を防止するバッキング材等を有している。当該超音波プローブ１１から被検体に超音波が送信されると、生体組織の非線形性により、超音波の伝播に伴って種々のハーモニック成分が発生する。送信超音波を構成する基本波とハーモニック成分は、体内組織の音響インピーダンスの境界、微小散乱等により後方散乱され、反射波（エコー）として超音波プローブ１１に受信される。

送信ユニット１２は、図示しない遅延回路およびパルサ回路等を有している。パルサ回路では、所定のレート周波数ｆｒ Ｈｚ（周期；１／ｆｒ秒）で、送信超音波を形成するためのレートパルスが繰り返し発生される。また、遅延回路では、チャンネル毎に超音波をビーム状に集束し且つ送信指向性を決定するのに必要な遅延時間が、各レートパルスに与えられる。送信ユニット１２は、このレートパルスに基づくタイミングで、所定のスキャンラインに向けて超音波ビームが形成されるように振動子毎に駆動パルスを印加する。

受信ユニット１３は、図示していないアンプ回路、Ａ／Ｄ変換器、加算器等を有している。アンプ回路では、プローブ１１を介して取り込まれたエコー信号をチャンネル毎に増幅する。Ａ／Ｄ変換器では、増幅されたエコー信号に対し受信指向性を決定するのに必要な遅延時間を与え、その後加算器において加算処理を行う。この加算により、所定のスキャンラインに対応した超音波エコー信号を生成する。

Ｂモード処理ユニット１４は、受信ユニット１３から受け取った超音波エコー信号に対して包絡線検波処理を施すことにより、超音波エコーの振幅強度に対応したＢモード信号を生成する。

組織ドプラ処理ユニット１５は、受信ユニット１３から受け取ったエコー信号に対して直交検波処理、自己相関処理等を行い、遅延加算処理された超音波エコー信号のドプラ偏移成分に基づいて、被検体内で移動している組織の速度、分散、パワーに対応した組織ドプラ信号を求める。

表示制御ユニット１７は、Ｂモード信号の所定断層に係る次元分布を表したＢモード超音波像を生成する。また、表示制御ユニット１７は、組織ドプラ信号に基づいて、速度、分散、パワー値の所定断層にかかる２次元分布を表した組織ドプラ超音波像を生成する。さらに、表示制御ユニット１７は、必要に応じてＢモード超音波像と組織ドプラ超音波像の重畳画像、Ｂモード超音波像と組織の運動情報に関する画像との重畳画像等を生成する。ここで、組織の運動情報とは、組織の歪み、歪み率、移動距離、速度その他の組織の運動に関して取得可能な物理情報である。以下、このような組織の運動情報を含む画像の総称を「ＴＳＩ画像」と呼ぶ。

表示部１８は、表示制御ユニット１７からのビデオ信号に基づいて、生体内の形態学的情報や、血流情報を画像として表示する。また、造影剤を用いた場合には、造影剤の空間的分布、すなわち血流或いは血液の存在している領域を求めた定量的な情報量に基づいて、輝度画像やカラー画像として表示する。

ＴＳＩ処理ユニット２０は、Ｂモード処理ユニット１４が出力するＢモード信号、組織ドプラ処理ユニット１６が出力するドプラ信号、記憶部２２に記憶されたＢモード画像データ、ドプラ画像データ、速度分布画像等を用いてＴＳＩ処理を実行し、組織の歪みに関する画像であるＴＳＩ画像を生成する。ここで、速度分布画像とは、診断対象組織上の複数位置における速度を時相毎に表す画像である。

ボリュームデータ生成ユニット２１は、必要に応じて、ＴＳＩ処理ユニット２０において生成された所定の対象に関する複数の断層に関するＴＳＩ画像を用いて、空間的な補間処理を実行する。これにより、当該所定の対象に関するＴＳＩ画像のボリュームデータが生成される。

記憶部２２は、各時相に対応する超音波画像データ（例えば、組織ドプラモード、Ｂモード等によって撮影された組織画像データ）、ＴＳＩ処理ユニット２０によって生成された各時相に対応する速度分布画像及びＴＳＩ画像等を記憶する。また、記憶部２２は、ボリュームデータ生成ユニット２１において生成されたＴＳＩ画像のボリュームデータを、必要に応じて記憶する。なお、記憶部２２が記憶する超音波画像データは、スキャンコンバート前の所謂Raw画像データであってもよい。

制御ユニット（ＣＰＵ）２３は、情報処理装置（計算機）としての機能を持ち、本超音波診断装置本体の動作を静的又は動的に制御する。

マッピング処理ユニット２４は、ボリュームデータ生成ユニット２１において生成されたＴＳＩ画像のボリュームデータに基づいて、選択された短軸断層上にShortening情報をマッピングするマッピング処理を行う。また、マッピング処理ユニット２４は、同じくＴＳＩ画像のボリュームデータに基づいて、選択された長軸断層上にThickening情報をマッピングするマッピング処理を行う。

入力部２５は、装置本体に接続され、オペレータからの各種指示、関心領域（ＲＯＩ）の設定指示、種々の画質条件設定指示等を装置本体にとりこむためのマウスやトラックボール、モード切替スイッチ、キーボード等を有している。

（組織歪みイメージング：ＴＳＩ）
次に、本実施形態の前提となる技術である組織歪みイメージング法について、簡単に説明する。この組織歪みイメージング法は、組織の運動情報として、運動に伴う組織位置を追跡しながら、速度情報に由来する信号を積分することで得られる、局所の変位、歪み等のパラメータを画像化するものである。当該手法によれば、心臓の局所心筋の歪みや変位に関する画像を作成・表示することができ、画像出力値の局所領域に対する時間変化の解析が支援される。

例えば短軸断層像を用いてThickeningを観察する場合、組織歪みイメージング法では、このThickeningに関わる成分を角度補正によって検出して画像化するために、収縮中心に向かう運動の場（コントラクション場：Contraction Motion Field）の概念や設定を用いている。また、組織歪みイメージング法では、心臓全体の並進運動（「translation」とも呼ばれる）の影響も考慮して収縮中心位置を時間的に移動させ、時間的に可変な運動場に対しても適用可能となっている。従って、並進運動による収縮中心位置の変動には追随することができる。この組織歪みイメージング法のさらなる詳細については、例えば特開平２００３−１７５０４１号に説明されている。

なお、本組織歪みイメージング法では、複数の時相に関する速度分布画像が必要とされる。この速度分布画像は、組織ドプラ法によって収集された複数の時相に関する２次元又は３次元超音波画像データから生成するか、Ｂモード等によって収集された複数の時相に関する複数の２次元又は３次元組織画像に対してパターンマッチング処理を施すことで得られる。本実施形態では、説明を具体的にするため、時相毎のＢモード画像をパターンマッチング処理することによって生成された２次元の速度分布画像を用いるものであるとする。しかしながら、これに限定する趣旨ではなく、例えば組織ドプラ法によって生成された２次元又は３次元の速度分布画像を利用するものであってもよい。

（ボリュームデータの取得）
後述するマッピング処理においては、例えば診断対象（今の場合、心臓）に関するＴＳＩ画像のボリュームデータ等が必要である。そのため、本超音波診断装置１０は、ＴＳＩ画像のボリュームデータを取得するためのいくつかの機能を有している。以下、この機能について説明する。

図２は、ＴＳＩ画像のボリュームデータを取得するための第１の手法を示しており、実際に診断部位を含む３次元領域を走査するボリュームスキャンによってボリュームデータを取得する手法を示した図である。すなわち、実際に診断部位を含む３次元領域を、２次元アレイプローブ（すなわち、超音波振動子がマトリックス状に配列されたプローブ）を用いたボリュームスキャン、又は一次元アレイプローブを機械的又は手技等により煽り動作させることによるボリュームスキャンすることで、ボリュームデータを取得することができる。また、ボリュームスキャンを連続的に行うことで、ボリュームデータを時相毎に取得でき、これらを用いたＴＳＩ処理を行うことで、診断対象に関する時相毎のＴＳＩ画像のボリュームデータを取得することができる。

図３、図４は、ボリュームデータを取得するための第２の手法を示している。各図に示すように、例えば心腔２断層（2-CH）、心尖長軸像（3-CH）、心腔４断層（4-CH）の三つの長軸断層像を時相毎に取得し、これらを用いたＴＳＩ処理を行う。得られた三つの長軸断層像に関するＴＳＩ画像を補間処理することで、診断対象に関するＴＳＩ画像のボリュームデータを時相毎に取得することができる。

図５、図６は、ボリュームデータを取得するための第３の手法を示している。各図に示すように、例えば心基部レベル（Basal）、中央レベル（Mid）、心尖部レベル（Apex）の三つの断軸像を時相毎に取得し、これを用いてこれらを用いたＴＳＩ処理を行う。得られた三つの短軸断層像に関するＴＳＩ画像を補間処理することで、診断対象に関するＴＳＩ画像のボリュームデータを時相毎に取得することができる。

なお、第１の手法と第２及び第３の手法とを比較した場合、装置の仕様によっては、第２及び第３の手法は第１の手法に比して空間分解能が劣る一方、時間分解能が勝る場合がある。係る場合には、診断の目的に応じていずれの手法を採用すべきかを決定する必要がある。また、第２及び第３の手法においては、三断層に基づく補間処理を例に説明した。しかしながら、これに拘泥されず、補間可能であれば、撮影にて取得する（補間処理の基礎とする）断層像の数は、いくつであってもよい。

（直交方向のマッピング処理）
次に、本超音波診断装置１が有する直交方向のマッピング機能について説明する。この機能は、２次元断層画像を表示する場合において、ＴＳＩ画像のボリュームデータを用いて、断層画像と直交する方向に関する運動情報を当該断層画像に色づけしてマッピング（カラーマッピング）し、画像化するものである。これにより、断層画像上の走査線に沿った方向の運動と当該断層画像に直交する方向に関する運動とを、同一の２次元断層画像上において観察することが可能となる。

図７、図８は、直交方向のマッピング機能を説明するための概念図の一例である。図７に示すように、例えば異なる三つの長軸断層画像（図３、図４参照）に関するＴＳＩ画像を用いた補間処理により、診断対象に関するボリュームデータＶ１が生成されているとする。このボリュームデータＶ１に対して任意の短軸断層面Ｅ１を指定すると共に、ボリュームデータＶ１を用いて当該短軸断層面Ｅ１上の組織上の各点のShortening情報を例えばその歪みの大きさに応じた輝度によってカラー変換し、マッピング画像Ｆ１を生成する。このマッピング画像Ｆ１を、例えば短軸断層面Ｅ１の通常の超音波画像と重畳して表示することで、短軸断層画像においてShortening情報を観察することができる。

また、図９、図１０は、直交方向のマッピング機能を説明するための概念図の他の例である。図９に示すように、例えば異なる三つの短軸断層画像（図５、図６参照）に関するＴＳＩ画像を用いた補間処理により、診断対象に関するボリュームデータＶ２が生成されているとする。このボリュームデータＶ２に対して任意の長軸断層面Ｅ２を指定すると共に、ボリュームデータＶ２を用いて当該長軸断層面Ｅ２上の組織上の各点のThickening情報を例えばその歪みの大きさに応じた輝度によってカラー変換し、マッピング画像Ｆ２を生成する。このマッピング画像Ｆ２を、例えば長軸断層面Ｅ２の通常の超音波画像と重畳して表示することで、長軸断層画像においてThickening情報を観察することができる。

上記各例に従う直交方向のマッピング処理は、診断対象のＴＳＩ画像に関するボリュームデータを用いたマッピングである。このＴＳＩ画像に関するボリュームデータの取得は、補間処理を含む既述の手法によって可能である。従って、本手法によれば、例えば図１１上段に示すように複数の長軸（複数の短軸）断層画像を撮影するのみで、図１１下段に示すように短軸（長軸）断層画像上にShortening情報をマッピングして映像化することができる。

（実施例１）
次に、本超音波診断装置１の直交方向のマッピング処理における動作について説明する。本実施例１は、短軸断層画像上に歪みのShortening情報をマッピングする場合を例として説明する。

図１２は、本実施例に係る直交方向のマッピング処理において実行される各処理の流れを示したフローチャートである。同図に示すように、まず、ＴＳＩ画像のボリュームデータに対して、任意の（所望の）位置における短軸断層が選択される（ステップＳ１）。

次に、マッピング処理ユニット２４は、ボリュームデータを用いて、指定された短軸断層上に存在する組織の各位置に関する長軸方向の縮み（Shortening情報）を取得し、これを歪みの強度とその濃度とが対応するようなカラーマッピングを行う。このマッピング処理により、マッピング画像が生成される（ステップＳ２）。

次に、表示制御ユニット１７は、指定された短軸断層に関する通常の超音波画像（短軸断層像）と生成されたマッピング画像とを重畳させて、表示ユニット１８の画面に表示する（ステップＳ３）。

（実施例２）
次に、本超音波診断装置１の直交方向のマッピング処理における動作について説明する。本実施例１は、長軸断層画像上に歪みのThickening情報をマッピングする場合を例として説明する。

図１３は、本実施例に係る直交方向のマッピング処理において実行される各処理の流れを示したフローチャートである。同図に示すように、まず、ＴＳＩ画像のボリュームデータに対して、任意の（所望の）位置における長軸断層が選択される（ステップＳ１１）。

次に、マッピング処理ユニット２４は、ボリュームデータを用いて、指定された長軸断層上に存在する組織の各位置に関する短軸方向の伸張（Thickening情報）を取得し、これを歪みの強度とその濃度とが対応するようなカラーマッピングを行う。このマッピング処理により、マッピング画像が生成される（ステップＳ１２）。

次に、表示制御ユニット１７は、指定された長軸断層に関する通常の超音波画像（長軸断層像）と生成されたマッピング画像とを重畳させて、表示ユニット１８の画面に表示する（ステップＳ１３）。

以上述べた構成によれば、以下の効果を得ることができる。

本超音波診断装置では、２次元断層画像を表示する場合において、断層画像と直交する方向に関する運動情報を当該断層画像にマッピングし、画像化する。これにより、従来取得できなかった表示情報が得られ、互いに直交し異なる起序を呈すると言われるShortening情報とThickening情報とを複合的に判断することが可能となる。その結果、心臓の画像診断において、局所的な心筋壁運動の把握が容易となり、心筋梗塞等の診断を支援することができる。

また、本超音波診断装置によれば、例えば長軸（短軸）断層画像のＴＳＩ画像を用いた補間処理により、診断対象のＴＳＩ画像に関するボリュームデータを生成する。従って、例えば複数の長軸（複数の短軸）断層画像を撮影するのみで、２次元画像上にShortening（Thickening）情報をマッピングし、Shortening情報とThickening情報とを複合的に判断することが可能となる。その結果、長軸断層画像及び短軸断層画像の双方を撮影する必要がなく、画像診断における作業時間を短縮することができ、術者及び被験者の負担を軽減させることができる。

（第２実施形態）
次に、本発明の第２実施形態について説明する。本実施形態に係る超音波診断装置１は、任意の断層画像上に、Shortening情報とThickening情報とを同時にマッピングし映像化しうるものである。

図１４は、第２の実施形態に係る超音波診断装置１０の構成図である。図１の超音波診断装置１０と比較した場合、追跡処理ユニット２７をさらに具備する点が異なる。この追跡処理ユニット２７は、例えばＢモード撮影によって取得された時相毎の組織像において、２次元又は３次元のスペックル追跡（トラッキング処理）を行い、心筋の壁厚方向のThickening情報等を演算する。

図１５は、本実施形態のマッピング処理の概念を説明するための図である。図１５に示すように、既述の方法にて診断対象のＴＳＩ画像に関するボリュームデータＶを生成し、これを用いてShortening情報を演算する。

次に、任意の短軸断層像Ｅ１を指定し、当該短軸断層内での白黒像のスペックル追跡により、壁方向のThickening情報を演算する。また、これと並列して、任意の長軸断層像Ｅ２を指定し、当該長軸断層内での白黒像のスペックル追跡により、壁方向のThickening情報を演算する。

次に、演算されたShortening情報の伸張成分を赤色、収縮成分を青色、Thickening情報の厚さ増加成分を緑色といった具合に、運動に応じた色を割り当てた後カラー変換することで、マッピング画像Ｆ１、Ｆ２を生成する。生成された各マッピング画像は、背景の白黒画像（短軸断層像又は長軸断層像）に重畳させて表示される。

なお、以上の手順を図１６のフローチャートに示した。

以上述べた構成によれば、必要に応じてShortening情報がマッピングされた短軸断層像、及びThickening情報がマッピングされた長軸断層像を同時に表示することができる。従って、観察者は、短軸断層像及び長軸断層像の双方でThickening情報とShortening情報とを観察することができ、心臓の画像診断において、局所的な心筋壁運動の把握が容易となり、心筋梗塞等の診断を好適に支援することができる。

（第３の実施形態）
第１の実施形態、第２の実施形態では、Thickening情報、Shortening情報を運動情報とする例を示した。しかしながら、本超音波診断装置のマッピング機能に用いることができる運動情報は、Thickening情報、Shortening情報に限定されない。そこで、本実施形態では、Thickening情報及びShortening情報以外の運動情報を用いて、第１の実施形態で述べたマッピング処理、及び第２の実施形態で述べたマッピング処理を実行する例について説明する。

本実施形態に係る超音波診断装置１０は、図１４と略同様の構成を有する。以下、各運動情報を用いたマッピング機能について説明する。

（運動情報：捻れ運動）
まず、運動情報として心臓壁の捻れ運動の状態を採用する場合の例を図２１、図２２を参照しつつ説明する。表示ユニット１８に表示された断層画像上に計測画像領域が指定されると、追跡処理ユニット２７は、時系列に沿った複数のボリュームデータを構成する断層画像を用いて内膜ｍ１の２次元トラッキングを実行することにより、心尖部短軸像Ｇ１の断面に直交する方向（心臓の長軸方向）を軸とする内膜ｍ１の回転角度（局所運動情報）を演算する（ステップＳ３１）。同様に、乳頭筋短軸像Ｇ２及び心基部短軸像Ｇ３についても、心臓の長軸方向を軸とする内膜ｍ２、ｍ３の回転角度（局所運動情報）を演算する（ステップＳ３２））。なお、内膜ｍ１、ｍ２、ｍ３に代えて、外膜Ｍ１、Ｍ２、Ｍ３の回転角度を演算してもよい。

このとき、追跡処理ユニット２７は、たとえば各時相について、ステップＳ３１で内膜ｍ１等が入力された時相（基準時相）に対する回転角度として内膜ｍ１、ｍ２、ｍ３の回転角度を演算する。また、時系列に沿って隣接するフレーム（つまり連続するフレーム）における内膜ｍ１、ｍ２、ｍ３の回転角度を順次演算していくようにしてもよい。

ＴＳＩ処理ユニット２０は、内膜ｍ１の回転角度と、内膜ｍ２の回転角度との差（相対回転角度）を演算する（ステップＳ３３）。同様に、内膜ｍ２の回転角度と、内膜ｍ３の回転角度との差（相対回転角度）を演算する（ステップＳ３４）。

ステップＳ３３、ステップＳ３４の処理について、図２２を参照しつつ具体的に説明する。入力ユニット２５からの所定の操作により、例えば反時計回り方向を正の回転方向（＋θ方向）と定義する。また、内膜ｍ１の回転角度をθ１、内膜ｍ２の回転角度をθ２、内膜ｍ３の回転角度をθ３とする。

このとき、ステップＳ３３で演算される相対回転角度Δθ１２は、Δθ１２＝θ１−θ２（又は、θ２−θ１）によって算出される。また、ステップＳ３４で演算される相対回転角度Δθ２３は、Δθ２３＝θ２−θ３（又は、θ３−θ２）によって算出される。

ステップＳ３３で得られる相対回転角度Δθ１２は、心尖部短軸像Ｇ１の断面位置と乳頭筋短軸像Ｇ２の断面位置との間における心臓壁の捻れ運動の状態（大きさ）を反映する情報である。つまり、相対回転角度Δθ１２＝０（θ１＝θ２）である場合、これらの断面位置の間の任意の位置において、心臓壁は同じ方向に同じ角度だけ回転しており、回転方向への捻れは無いものと考えることができる。

一方、｜Δθ１２｜≠０である場合には、これらの断面位置の間において回転角度に差異があり、心臓壁が回転角度方向に捻れていることになる。この心臓壁の捻れは、相対回転角度Δθ１２の絶対値が大きいほど大きくなる。たとえば、θ１の符号とθ２の符号とが異なる場合、すなわち内膜ｍ１の回転方向と内膜ｍ２の回転方向とが逆である場合には、相対回転角度Δθ１２の絶対値は比較的大きくなる。

ステップＳ３４で得られる相対回転角度Δθ２３は、同様に、乳頭筋短軸像Ｇ２の断面位置と心基部短軸像Ｇ３の断面位置との間における心臓壁の捻れ運動の大きさを反映する情報である。

マッピング処理ユニット２４は、ステップＳ３３、ステップＳ３４で演算された相対回転角度Δθ１２、Δθ２３を、心臓壁の捻れ運動の大きさを示す運動情報とし、これを例えば長軸断層像等にカラーマッピングすることでマッピング画像を生成する。表示制御ユニット１７は、このマッピング画像を、組織の形態画像又は所定の運動情報がマッピングされたマッピング画像と重畳させて表示する（ステップＳ３５）。この表示された相対回転角度Δθ１２、Δθ２３を参照することにより、ユーザは、心臓壁の捻れ運動の大きさを把握することができる。ここで、心臓壁の内膜及び外膜についてそれぞれ相対回転角度を演算し、この２つの相対回転角度に基づいて捻れ運動の大きさを評価することもできる（例えば２つの相対回転角度の平均値を取るなど。）。

なお、相対回転角度Δθ１２を時間で微分することにより、内膜ｍ１、ｍ２の間における心臓壁の捻れ運動の速度を求めることができる。同様に、相対回転角度Δθ２３を時間で微分することにより、内膜ｍ２、ｍ３の間における心臓壁の捻れ運動の速度を求めることができる。そして、これらの速度を表示部８１に表示させるように構成することが可能である。ここで、「微分」とは、通常の微分演算とともに、相対回転角度を求めたフレーム間の時間間隔で当該相対回転角度を除算する処理も含むものとする。

（運動情報：相対回転勾配）
心臓壁の相対回転勾配を運動情報として取得するときの処理について、図２３、図２４を参照しつつ説明する。この相対回転勾配は、心臓壁の捻れ運動の度合いを示す運動情報である。

まず、追跡処理ユニット２７は、心尖部短軸像Ｇ１の内膜ｍ１の回転角度θ１と、乳頭筋短軸像Ｇ２の内膜ｍ２の回転角度θ２と、心基部短軸像Ｇ３の内膜ｍ３の回転角度θ３とをそれぞれ演算する（ステップＳ４１）。

次に、追跡処理ユニット２７は、内膜ｍ１の回転角度θ１と、内膜ｍ２の回転角度θ２との相対回転角度Δθ１２を演算し（ステップＳ４２）、内膜ｍ２の回転角度θ２と、内膜ｍ３の回転角度θ３と相対回転角度Δθ２３を演算する（ステップＳ４３）。

追跡処理ユニット２７は、心尖部短軸像Ｇ１と乳頭筋短軸像Ｇ２との間の距離ｄ１２を演算し（ステップＳ４４）、乳頭筋短軸像Ｇ２と心基部短軸像Ｇ３との間の距離ｄ２３を演算する（ステップＳ４５）。この距離ｄ１２、ｄ２３は、たとえば心尖部短軸像Ｇ１、乳頭筋短軸像Ｇ２、心基部短軸像Ｇ３の断面位置の座標に基づいて演算することができる。

次に、ＴＳＩ処理ユニット２０は、図２４に示すように、ステップＳ４２で求めた相対回転角度Δθ１２を、ステップＳ４４で求めた距離ｄ１２で除算して、内膜ｍ１と内膜ｍ２との間における相対回転勾配δθ１２＝Δθ１２／ｄ１２を演算する（ステップＳ４６）。同様に、ＴＳＩ処理ユニット２０は、ステップＳ４３で求めた相対回転角度Δθ２３を、ステップＳ４５で求めた距離ｄ２３で除算して、内膜ｍ２と内膜ｍ３との間における相対回転勾配δθ２３＝Δθ２３／ｄ２３を演算する（ステップＳ４７）。

マッピング処理ユニット２０は、ステップＳ４６、ステップＳ４７で演算された相対回転勾配δθ１２、δθ２３を、心臓壁の捻れ運動の度合いを示す運動情報とし、これを例えば長軸断層像等にカラーマッピングすることでマッピング画像を生成する。表示制御ユニット１７は、このマッピング画像を、組織の形態画像又は所定の運動情報がマッピングされたマッピング画像と重畳させて表示する（ステップＳ４８）。

相対回転勾配δθ１２は、心尖部レベルの内膜と乳頭筋レベルの内膜の間における単位距離当たりの捻れの大きさを示している。また、相対回転勾配δθ２３は、乳頭筋レベルの内膜の心基部レベルの内膜と間における単位距離当たりの捻れの大きさを示している。つまり、相対回転勾配δθ１２、δθ２３は、心臓壁（内膜）の捻れの度合いを反映した運動情報である。ユーザは、表示された相対回転勾配δθ１２、δθ２３を参照することにより、心臓壁の捻れ運動の度合いを把握することができる。なお、心臓壁の内膜及び外膜についてそれぞれ相対回転勾配を演算し、この２つの相対回転勾配に基づいて捻れ運動の度合いを評価することもできる（例えば、２つの相対回転勾配の平均値を取るなど。）。

（運動情報：長軸方向のストレイン）
心臓壁の長軸方向のストレインを運動情報として取得するときの処理について、図２５、図２６を参照しつつ説明する。このストレインは、心臓壁の歪みの大きさの度合いを示す情報であり、心臓壁の歪み状態を示すものである。

まず、追跡処理ユニット２７は、内膜ｍ１、内膜ｍ２、内膜ｍ３のそれぞれについて、計測画像領域が指定された断層画像３次元的な変位（Δｘ１、Δｙ１、Δｚ１）、（Δｘ２、Δｙ２、Δｚ２）、（Δｘ３、Δｙ３、Δｚ３）を演算し（ステップＳ５１）、これらの３次元的な変位から、Ｚ方向（長軸方向）への変位Δｚ１、Δｚ２、Δｚ３をそれぞれ抽出する（ステップＳ５２）。

次に、追跡処理ユニット２７は、心尖部レベルと乳頭筋レベルとの間における心臓壁の伸縮Δｚ１２＝Δｚ１−Δｚ２を演算し（ステップＳ５３）、乳頭筋レベルと心基部レベルとの間における心臓壁の伸縮Δｚ２３＝Δｚ２−Δｚ３を演算する（ステップＳ５４）。

また、追跡処理ユニット２７は、計測画像領域が指定された心尖部短軸像Ｇ１、乳頭筋短軸像Ｇ２、心基部短軸像Ｇ３について、心尖部短軸像Ｇ１と乳頭筋短軸像Ｇ２との間の距離ｄ１２を演算し（ステップＳ５５）、乳頭筋短軸像Ｇ２と心基部短軸像Ｇ３との間の距離ｄ２３を演算する（ステップＳ５６）。

ＴＳＩ処理ユニット２０は、ステップＳ５３で演算した伸縮Δｚ１２をステップＳ４５で演算した距離ｄ１２で除算することにより、心尖部レベルと乳頭筋レベルとの間における長軸方向のストレインδｚ１２＝Δｚ１２／ｄ１２を演算する（ステップＳ５７）。また、ＴＳＩ処理ユニット２０は、ステップＳ４４で演算した伸縮Δｚ２３をステップＳ４６で演算した距離ｄ２３で除算することにより、乳頭筋レベルと心基部レベルとの間における長軸方向のストレインδｚ２３＝Δｚ２３／ｄ２３を演算する（ステップＳ５８）。

マッピング処理ユニット２０は、ステップＳ５７、ステップＳ５８で演算された心臓壁のストレインδｚ１２、δｚ２３を、心臓壁のストレインの大きさを示す運動情報とし、これを例えば短軸断層像等にカラーマッピングすることでマッピング画像を生成する。表示制御ユニット１７は、このマッピング画像を、組織の形態画像又は所定の運動情報がマッピングされたマッピング画像と重畳させて表示する（ステップＳ５９）。ユーザは、表示された心臓壁のストレインδｚ１２、δｚ２３を参照することにより、心臓壁の歪みの大きさを把握することができる。

なお、心臓壁の内膜及び外膜についてそれぞれストレインを演算し、この２つのストレインの値に基づいて歪みの大きさを評価することもできる（例えば２つのストレインの値の平均値を取るなど。）。

（運動情報：長軸方向のストレインレート）
心臓壁の長軸方向のストレインレートを運動情報として取得するときの処理について説明する。このストレインレートは、心臓壁の歪み（ストレイン）の時間変化率を示す情報であり、心臓壁の歪み状態を示すものである。

ストレインレートを求める場合、図２５のフローチャートのステップＳ５１〜Ｓ５８と同様の処理を行って、心尖部レベルと乳頭筋レベルとの間における長軸方向のストレインδｚ１２と、乳頭筋レベルと心基部レベルとの間における長軸方向のストレインδｚ２３とを演算する。

ここで、ストレインδｚ１２、ストレインδｚ２３は、２つの時相ｔ１、ｔ２における心尖部短軸像Ｇ１、乳頭筋短軸像Ｇ２、心基部短軸像Ｇ３について演算されるものである。ＴＳＩ処理ユニット２０は、ストレインδｚ１２を時間間隔Δｔ＝｜ｔ１−ｔ２｜で除算することにより、心尖部レベルと乳頭筋レベルとの間における長軸方向のストレインレートを演算する。また、ストレインδｚ２３を時間間隔Δｔで除算することにより、乳頭筋レベルと心基部レベルとの間における長軸方向のストレインレートを演算する。なお、通常の微分演算を実行することにより、ストレインからストレインレートを算出するようにしてもよい。

マッピング処理ユニット２０は、演算された心臓壁のストレインレートδｚ１２／Δｔ、δｚ２３／Δｔを、心臓壁のストレインの時間変化率を示す運動情報とし、これを例えば短軸断層像等にカラーマッピングすることでマッピング画像を生成する。表示制御ユニット１７は、このマッピング画像を、組織の形態画像又は所定の運動情報がマッピングされたマッピング画像と重畳させて表示する。ユーザは、表示された心臓壁のストレインレートを参照することにより、心臓壁の歪みの時間変化率を把握することができる。

なお、心臓壁の内膜及び外膜についてそれぞれストレインレートを演算し、この２つのストレインレートの値に基づいて歪みの時間変化率を評価することもできる（例えば、２つのストレインレートの値の平均値を取るなど。）。

以上述べた構成によれば、捻れ情報等の運動情報を用いた場合であっても、第１の実施形態で述べたマッピング処理、及び第２の実施形態で述べたマッピング処理を実現することができる。

なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。具体的な変形例としては、例えば次のようなものがある。

（１）上記各実施形態においては、マッピング画像を当該断面の白黒像と重畳させて表示する構成とした。これに対し、一枚の短軸断層像又は長軸断層像にThickening情報とShortening情報との双方をマッピングしたマッピング画像を生成し、これを表示するようにしてもよい。この場合には、例えば図１７に示すように、Thickening情報とShortening情報とで異なる色相を割り当て、その歪みの強度に応じた濃度によってこれを表示するようにすればよい。

（２）上記各実施形態においては、Thickening情報又はShortening情報を色相表示するためのマッピング画像を生成し表示する構成であった。これに対し、例えば図１８、図１９に示す様に、マッピング画像上の心臓組織を小領域に分割し、当該小領域毎にその運動情報を数値化して表示するようにしてもよい。

また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。

（３）上記各実施形態においては、説明を具体的にするため、３次元スキャン等によって取得されてボリュームデータを用いて診断対象に関するＴＳＩ画像のボリュームデータを生成し、これを用いた直交方向のマッピング処理を説明した。

た。しかしながら、直交方向のマッピング処理は、必ずしも診断対象に関するＴＳＩ画像のボリュームデータを必要としない。すなわち、本直交方向のマッピング処理は、マッピング画像を生成するために選択される短軸断層像又は長軸断層像に表示された組織の各点において、当該断層面と直交する方向についての運動情報を取得しうる画像データがあればよい。

従って、例えば図２０に示すように、マッピング画像を生成するために選択される断層像Ｅ１を含み、当該断層像と直交する方向に関して一定の厚さを有する画像データｖによっても、本直交方向のマッピング処理を実現することができる。なお、この画像データｖを取得する場合、断層像Ｅ１に略直交する複数の断層像を取得し既述の補間処理を行うか、少なくとも当該画像データｖに対応する領域を超音波走査すれば十分である。

（４）上記各実施形態においては、ＴＳＩ画像を利用する場合を例を説明した。しかしながら、これに拘泥されず、組織の運動情報に関する画像であれば本直交方向のマッピング処理を適用することができる。従って、例えば組織ドプラ画像を用いても本直交方向のマッピング処理は実施可能であり、同様の効果を得ることができる。

（５）上記各実施形態においては、マッピング処理によって得られたマッピング画像を他の画像と重畳させて表示する例を説明した。しかしながら、これに拘泥されず、マッピング処理によって得られる複数のマッピング画像を用いて投影画像（例えば、ボリュームレンダリング画像、サーフェスレンダリング画像、極座標等による展開画像）を生成し、これを表示するようにしてもよい。

また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。

以上本発明によれば、短軸断層像を撮影し観察している場合であってもshortening情報を同時に観察することができ、また、長軸断層像を撮影し観察している場合であってもthickening情報を同時に観察することができる超音波診断装置、超音波画像処理装置及び超音波画像処理プログラムを実現することができる。

図１は、第１の実施形態に係る超音波診断装置１０の構成図である。 図２は、ＴＳＩ画像のボリュームデータを取得するための第１の手法を示している。 図３は、ボリュームデータを取得するための第２の手法を説明するための図である。 図４は、ボリュームデータを取得するための第２の手法を説明するための図である。 図５は、ボリュームデータを取得するための第３の手法を説明するための図である。 図６は、ボリュームデータを取得するための第３の手法を説明するための図である。 図７は、直交方向のマッピング機能を説明するための概念図の一例である。 図８は、直交方向のマッピング機能を説明するための概念図の一例である。 図９は、直交方向のマッピング機能を説明するための概念図の他の例である。 図１０は、直交方向のマッピング機能を説明するための概念図の他の例である。 図１１は、直交方向のマッピング機能を説明するための概念図である。 図１２は、実施例１に係る直交方向のマッピング処理において実行される各処理の流れを示したフローチャートである。 図１３は、実施例２に係る直交方向のマッピング処理において実行される各処理の流れを示したフローチャートである。 図１４は、第２の実施形態に係る超音波診断装置１０の構成図である。 図１５は、第２実施形態のマッピング処理の概念を説明するための図である。 図１６は、第２の実施形態に係るマッピング処理の流れを示したフローチャートである。 図１７は、マッピング画像におけるThickening情報又はShortening情報の表示形態の変形例を示した図である。 図１８は、マッピング画像におけるThickening情報又はShortening情報の表示形態の変形例を示した図である。 図１９は、マッピング画像におけるThickening情報又はShortening情報の表示形態の変形例を示した図である。 図２０は、本発明の実施形態の変形例を説明するための図である。 図２１は、捻れ運動を運動情報とした場合のマッピング処理の流れを示したフローチャートである。 図２２は、捻れ運動を運動情報とした場合のマッピング処理を説明するための図である。 図２３は、相対回転勾配を運動情報とした場合のマッピング処理の流れを示したフローチャートである。 図２４は、相対回転勾配を運動情報とした場合のマッピング処理を説明するための図である。 図２５は、長軸方向のストレインを運動情報とした場合のマッピング処理の流れを示したフローチャートである。 図２６は、長軸方向のストレインを運動情報とした場合のマッピング処理を説明するための図である。

符号の説明

１０…超音波診断装置、１１…超音波プローブ、１２…送信ユニット、１３…受信ユニット、１４…Ｂモード処理ユニット、１５…組織ドプラ処理ユニット、１７…表示制御ユニット、１８…表示ユニット、２０…ＴＳＩ処理ユニット、２１…ボリュームデータ生成ユニット、２２…記憶部、２３…制御ユニット（ＣＰＵ）、２４…マッピング処理ユニット、２５…入力ユニット２５

Claims (18)

1. ３次元空間における、心組織の第１の方向及び第２の方向に関する運動情報を求める運動情報生成ユニットと、
   超音波送受信により取得された前記３次元空間に関するボリュームデータに基づいて、前記３次元空間における第１の方向、及び前記第１の方向に対し直交する第２の方向に関する２次元画像を生成する２次元画像生成ユニットと、
   前記第１の方向に関する２次元画像に対して前記第２の方向に関する運動情報をマッピングした第１のマッピング画像と、前記第２の方向に関する２次元画像に対して前記第１の方向に関する運動情報をマッピングした第２のマッピング画像とを生成するマッピング画像生成ユニットと、
   前記第１あるいは第２のマッピング画像を表示する表示ユニットと、
   を具備することを特徴とする超音波診断装置。
2. 前記第１の方向を設定する設定ユニットを更に備え、
   前記第１の方向は、心臓の長軸方向に対応し、前記第２の方向は、心臓の短軸方向に対応することを特徴とする請求項１記載の超音波診断装置。
3. ３次元空間における心組織の長軸方向に関する運動情報と、前記長軸方向に対して直交する短軸方向に関する運動情報とを求める運動情報生成ユニットと、
   超音波送受信により取得された前記３次元空間に関するボリュームデータに基づいて、前記短軸断層に関する２次元画像である短軸画像と、前記長軸断層に関する２次元画像である長軸画像と、を生成する２次元画像生成ユニットと、
   前記短軸画像に前記心組織の長軸方向に関する運動情報をマッピングすることで第１のマッピング画像を生成し、前記長軸画像に前記心組織の短軸方向に関する運動情報をマッピングすることで第２のマッピング画像を生成するマッピング画像生成ユニットと、
   前記第１及び第２のマッピング画像を表示する表示ユニットと、
   を具備することを特徴とする超音波診断装置。
4. 前記運動情報生成ユニットは、前記ボリュームデータに対してスペックル追跡処理を施すことにより、前記３次元空間に関する運動情報を求めることを特徴とする請求項１乃至３のいずれか１項に記載の超音波診断装置。
5. 前記運動情報生成ユニットは、前記運動情報から前記心組織の短軸方向に関する運動情報をさらに求め、
   前記マッピング画像生成ユニットは、前記短軸断層上の心組織の長軸方向に関する運動情報及び前記短軸断層上の心組織の短軸方向に関する運動情報のそれぞれについて、運動の度合いに応じて異なる色彩が割り当てられた前記マッピング画像を生成すること、
   を特徴とする請求項２記載の超音波診断装置。
6. ３次元空間における心組織の長軸方向に関する運動情報と、前記運動情報から前記心組織の短軸方向に関する運動情報と、を求める運動情報生成ユニットと、
   超音波送受信により取得された前記３次元空間に関するボリュームデータに基づいて前記長軸方向に対して直交する短軸断層に関する２次元画像を生成する２次元画像生成ユニットと、
   前記２次元画像に対して、前記短軸断層上の心組織の長軸方向に関する運動情報及び前記短軸断層上の心組織の短軸方向に関する運動情報のそれぞれについて、運動の度合いに応じて異なる色彩が割り当てられたマッピング画像を生成するマッピング画像生成ユニットと、
   前記マッピング画像を表示する表示ユニットと、
   を具備することを特徴とする超音波診断装置。
7. 前記運動情報生成ユニットは、前記運動情報から前記心組織の短軸方向に関する運動情報をさらに求め、
   前記マッピング画像生成ユニットは、前記短軸断層上の心組織の長軸方向に関する運動情報については運動の度合いに応じて異なる色彩が割り当てられ、前記短軸断層上の心組織の短軸方向に関する運動情報については運動の度合いに応じた数値が割り当てられた前記マッピング画像を生成すること、
   を特徴とする請求項２記載の超音波診断装置。
8. ３次元空間における心組織の長軸方向に関する運動情報と、前記運動情報から前記心組織の短軸方向に関する運動情報と、を求める運動情報生成ユニットと、
   超音波送受信により取得された前記３次元空間に関するボリュームデータに基づいて前記長軸方向に対して直交する短軸断層に関する２次元画像を生成する２次元画像生成ユニットと、
   前記２次元画像に対して、前記短軸断層上の心組織の長軸方向に関する運動情報については運動の度合いに応じて異なる色彩が割り当てられ、前記短軸断層上の心組織の短軸方向に関する運動情報については運動の度合いに応じた数値が割り当てられたマッピング画像を生成するマッピング画像生成ユニットと、
   前記マッピング画像を表示する表示ユニットと、
   を具備することを特徴とする超音波診断装置。
9. ３次元空間における心組織の短軸方向に関する運動情報と、前記運動情報から前記心組織の長軸方向に関する運動情報と、を求める運動情報生成ユニットと、
   超音波送受信により取得された前記３次元空間に関するボリュームデータに基づいて前記短軸方向に対して直交する長軸断層に関する２次元画像を生成する２次元画像生成ユニットと、
   前記２次元画像に対して、前記長軸断層上の心組織の短軸方向に関する運動情報及び前記長軸断層上の心組織の長軸方向に関する運動情報のそれぞれについて、運動の度合いに応じて異なる色彩が割り当てられたマッピング画像を生成するマッピング画像生成ユニットと、
   前記マッピング画像を表示する表示ユニットと、
   を具備することを特徴とする超音波診断装置。
10. ３次元空間における心組織の短軸方向に関する運動情報と、前記運動情報から前記心組織の長軸方向に関する運動情報と、を求める運動情報生成ユニットと、
    超音波送受信により取得された前記３次元空間に関するボリュームデータに基づいて前記短軸方向に対して直交する長軸断層に関する２次元画像を生成する２次元画像生成ユニットと、
    前記２次元画像に対して、前記長軸断層上の心組織の長軸方向に関する運動情報については運動の度合いに応じて異なる色彩が割り当てられ、前記長軸断層上の心組織の長軸方向に関する運動情報については運動の度合いに応じた数値が割り当てられたマッピング画像を生成するマッピング画像生成ユニットと、
    前記マッピング画像を表示する表示ユニットと、
    を具備することを特徴とする超音波診断装置。
11. 前記運動情報生成ユニットは、前記心臓の複数の短軸断層上のそれぞれにおいて、前記心組織の長軸方向に関する運動情報を生成し、
    前記２次元画像生成ユニットは、前記複数の短軸断層のそれぞれに対応する前記２次元画像を生成し、
    前記マッピング画像生成ユニットは、前記各２次元画像に対して、対応する前記心組織の長軸方向に関する運動情報が少なくともマッピングされた複数の前記マッピング画像を生成し、
    前記表示ユニットは、前記複数のマッピング画像を表示すること、
    を特徴とする請求項６又は８記載の超音波診断装置。
12. 前記運動情報生成ユニットは、前記心臓の複数の長軸断層上のそれぞれにおいて、前記心組織の短軸方向に関する運動情報を生成し、
    前記２次元画像生成ユニットは、前記複数の長軸断層のそれぞれに対応する前記２次元画像を生成し、
    前記マッピング画像生成ユニットは、前記各２次元画像に対して、対応する前記心組織の短軸方向に関する運動情報が少なくともマッピングされた複数の前記マッピング画像を生成し、
    前記表示ユニットは、前記複数のマッピング画像を表示すること、
    を特徴とする請求項９又は１０記載の超音波診断装置。
13. 前記運動情報は、心臓組織の捻れ運動に関する情報、心臓組織の相対回転勾配に関する情報、心臓組織の歪み又は歪み率に関する情報、心臓組織の変位に関する情報のいずれかであることを特徴とする請求項１乃至１２のうちいずれか一項記載の超音波診断装置。
14. 前記運情報生成ユニットは、被検体を３次元スキャンすることで生成された３次元の運動情報、又は複数の２次元スキャンに基づく補間処理によって生成された運動情報のいずれかを生成することを特徴とする請求項１乃至１３のうちいずれか一項記載の超音波診断装置。
15. ３次元空間における、心組織の第１の方向及び第２の方向に関する運動情報を求める運動情報生成ユニットと、
    超音波送受信により取得された前記３次元空間に関するボリュームデータに基づいて、前記３次元空間における第１の方向、及び前記第１の方向に対し直交する第２の方向に関する２次元画像を生成する２次元画像生成ユニットと、
    前記第１の方向に関する２次元画像に対して前記第２の方向に関する運動情報をマッピングした第１のマッピング画像と、前記第２の方向に関する２次元画像に対して前記第１の方向に関する運動情報をマッピングした第２のマッピング画像とを生成するマッピング画像生成ユニットと、
    前記第１あるいは第２のマッピング画像を表示する表示ユニットと、
    を具備することを特徴とする超音波画像処理装置。
16. ３次元空間における心組織の長軸方向に関する運動情報と、前記長軸方向に対して直交する短軸方向に関する運動情報とを求める運動情報生成ユニットと、
    超音波送受信により取得された前記３次元空間に関するボリュームデータに基づいて、前記短軸断層に関する２次元画像である短軸画像と、前記長軸断層に関する２次元画像である長軸画像と、を生成する２次元画像生成ユニットと、
    前記短軸画像に前記心組織の長軸方向に関する運動情報をマッピングすることで第１のマッピング画像を生成し、前記長軸画像に前記心組織の短軸方向に関する運動情報をマッピングすることで第２のマッピング画像を生成するマッピング画像生成ユニットと、
    前記第１及び第２のマッピング画像を表示する表示ユニットと、
    を具備することを特徴とする超音波画像処理装置。
17. コンピュータに、
    ３次元空間における、心組織の第１の方向及び第２の方向に関する運動情報を求めさせる運動情報生成機能と、
    超音波送受信により取得された前記３次元空間に関するボリュームデータに基づいて、前記３次元空間における第１の方向、及び前記第１の方向に対し直交する第２の方向に関する２次元画像を生成させる２次元画像生成機能と、
    前記第１の方向に関する２次元画像に対して前記第２の方向に関する運動情報をマッピングした第１のマッピング画像と、前記第２の方向に関する２次元画像に対して前記第１の方向に関する運動情報をマッピングした第２のマッピング画像とを生成させるマッピング画像生成機能と、
    前記第１あるいは第２のマッピング画像を表示させる表示機能と、
    を実現させることを特徴とする超音波画像処理プログラム。
18. コンピュータに、
    ３次元空間における心組織の長軸方向に関する運動情報と、前記長軸方向に対して直交する短軸方向に関する運動情報とを求めさせる運動情報生成機能と、
    超音波送受信により取得された前記３次元空間に関するボリュームデータに基づいて、前記短軸断層に関する２次元画像である短軸画像と、前記長軸断層に関する２次元画像である長軸画像と、を生成させる２次元画像生成機能と、
    前記短軸画像に前記心組織の長軸方向に関する運動情報をマッピングすることで第１のマッピング画像を生成させ、前記長軸画像に前記心組織の短軸方向に関する運動情報をマッピングすることで第２のマッピング画像を生成させるマッピング画像生成機能と、
    前記第１及び第２のマッピング画像を表示させる表示機能と、
    を実現させることを特徴とする超音波画像処理プログラム。

Images