(57)【要約】【課題】精度の高い、組織の移動若しくは変形に関するタイミングを映像化すること。【解決手段】超音波を被検体Ｓ内にスキャンしたときの受信信号に基づいて制御部２１により被検体Ｓの動きの速度の時間的変化を求め、この速度の時間的変化から被検体Ｓにおける動きのタイミングの違いをモニター２０にカラー表示する。【選択図】図１
2006289067.tif 000002
【特許請求の範囲】
【請求項１】
被検体内に超音波を送信し、前記被検体からの反射波を受信して受信信号を出力するスキャン部と、
前記スキャン部から出力される前記受信信号に基づいて前記被検体の動きの速度の時間的変化を求め、この速度の時間的変化から前記被検体における動きのタイミングの違いを表示する演算表示部と、
を具備することを特徴とする超音波診断装置。
【請求項２】
前記スキャン部から出力される前記受信信号に基づいて前記被検体内の画像データを逐次作成する画像データ作成部を有し、
前記演算表示部は、前記画像データ作成部により逐次作成された前記画像データ上に制御点を付し、逐次作成される前記画像データ上の前記制御点の移動を追跡して前記被検体の動きの速度の時間的変化を求め、この速度の時間的変化から前記被検体における動きのタイミングの違いを表示する、
ことを特徴とする請求項１記載の超音波診断装置。
【請求項３】
前記画像データ作成部は、前記画像データを時間の経過と共に複数フレーム作成し、
前記演算表示部は、前記画像データにおける画素を前記制御点とし、後続して作成される前記各フレームの前記画像データ上における前記制御点の移動を追跡して前記制御点の速度の時間的変化を求め、この速度の時間的変化に基づいて前記被検体における動きのタイミングの違いを表示することを特徴とする請求項１記載の超音波診断装置。
【請求項４】
前記演算部は、前記制御点を前記画像データにおける全ての前記画素又は一部の前記画素に設定することを特徴とする請求項１記載の超音波診断装置。
【請求項５】
前記被検体は、組織から成り、
前記演算表示部は、前記被検体の前記組織の動きの速度の時間的変化の各ピーク位置を求める第１の処理部と、
前記第１の処理部により求められた前記各ピーク位置に基づいて前記組織の動きのタイミングの違いを求め、この動きのタイミングの違いを表示する第２の処理部と、
を有することを特徴とする請求項１記載の超音波診断装置。
【請求項６】
前記被検体の前記組織の動きのタイミングの違いをカラー表示するカラー表示部を有することを特徴とする請求項５記載の超音波診断装置。
【請求項７】
前記被検体は、収縮運動し、
前記カラー表示部は、前記被検体の収縮運動による前記被検体の各部位における前記組織の動きのタイミングの違いをカラー表示することを特徴とする請求項６記載の超音波診断装置。
【請求項８】
前記第１の処理部は、前記被検体の前記組織の動きの速度の時間的変化の重心又は標準偏差を求めることを特徴とする請求項５記載の超音波診断装置。
【請求項９】
ディスプレイ装置を有し、
前記カラー表示部は、前記画像データ作成部により逐次作成される前記画像データから超音波画像を前記ディスプレイ装置に表示し、かつ前記組織の動きのタイミングの違いの前記カラー表示を前記超音波画像に重ねて表示する、
ことを特徴とする請求項６記載の超音波診断装置。
【請求項１０】
被検体内に超音波を送信し、前記被検体からの反射波を受信して受信信号を出力するスキャン部と、
前記スキャン部から出力される前記受信信号に基づいて前記被検体の動きの速度を検出し、前記被検体の動きの速度によりstrain値、velocity gradient値又はdisplacement値の少なくとも１つの値を求め、この値の時間的変化から前記被検体における動きのタイミングの違いを表示する演算表示部と、
を具備することを特徴とする超音波診断装置。
【請求項１１】
前記スキャン部から出力される前記受信信号に基づいて前記被検体内の画像データを複数フレーム逐次作成する画像データ作成部を有し、
前記演算表示部は、前記画像データにおける画素を前記制御点とし、後続して作成される前記各フレームの前記画像データ上における前記制御点の前記strain値、前記velocity gradient値又は前記displacement値の少なくとも１つの値を求める、
ことを特徴とする請求項１０記載の超音波診断装置。
【請求項１２】
前記演算表示部は、前記制御点を前記画像データにおける全ての前記画素又は一部の前記画素に設定することを特徴とする請求項１０記載の超音波診断装置。
【請求項１３】
前記演算表示部は、前記strain値、前記velocity gradient値又は前記displacement値の少なくとも１つの値の変化のピーク位置を求める第１の処理部と、
前記第１の処理部により求められた前記ピーク位置に基づいて前記被検体における動きのタイミングの違いを求める第２の処理部と、
を有することを特徴とする請求項１０記載の超音波診断装置。
【請求項１４】
前記第１の処理部は、前記strain値、前記velocity gradient値又は前記displacement値の少なくとも１つの値の変化の重心又は標準偏差を求めることを特徴とする請求項１３記載の超音波診断装置。
【請求項１５】
前記被検体の前記組織の動きのタイミングの違いをカラー表示するカラー表示部を有することを特徴とする請求項１３記載の超音波診断装置。
【請求項１６】
前記被検体は、収縮運動し、
前記カラー表示部は、前記被検体の収縮運動による前記被検体の各部位における前記組織の動きのタイミングの違いをカラー表示することを特徴とする請求項１５記載の超音波診断装置。
【請求項１７】
ディスプレイ装置を有し、
前記カラー表示部は、前記画像データ作成部により逐次作成される前記画像データから超音波画像を前記ディスプレイ装置に表示し、かつ前記組織の動きのタイミングの違いの前記カラー表示を前記超音波画像に重ねて表示する、
ことを特徴とする請求項１５記載の超音波診断装置。
【請求項１８】
スキャン部から超音波を被検体内に送信し、前記被検体からの反射波を受信して受信信号を出力し、
前記スキャン部から出力される前記受信信号に基づいて前記被検体の動きの速度の時間的変化を求め、この速度の時間的変化から前記被検体における動きのタイミングの違いを表示する、
ことを特徴とする超音波診断装置の制御プログラム。
【請求項１９】
スキャン部から超音波を被検体内に送信し、前記被検体からの反射波を受信して受信信号を出力し、
前記受信信号に基づいて前記被検体の動きの速度を検出し、
前記検出された前記被検体の動きの速度によりstrain値、velocity gradient値又はdisplacement値の少なくとも１つの値を求め、この値の時間的変化から前記被検体における動きのタイミングの違いを表示する、
ことを特徴とする超音波診断装置の制御プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体内に超音波パルスを送信してスキャンし、被検体からの反射波を受信し、被検体内の画像に基づいて組織の運動情報を求める超音波診断装置及びその制御プログラムに関する。
【背景技術】
【０００２】
超音波診断装置は、超音波パルス反射法により体表から生体内の軟組織の断層像を無侵襲に得る。超音波診断装置は、超音波探触子から被検体内に超音波を送信し、被検体内部で音響インピーダンスの不整合により生じる反射波を超音波探触子で受信して受信信号を生成し、この受信信号に基づいて被検体の内部を画像化する。
【０００３】
この超音波診断装置は、他の医療用画像機器に比べて次の特長を有する。例えば小型で安価である。Ｘ線等の被爆がなく安全性が高い。血流イメージングが可能である。超音波診断装置は、例えば心臓、腹部の診断に用いられ、かつ泌尿器科及び産婦人科等で広く利用されている。超音波診断装置は、心臓の診断に有用であることが知られている。心筋等の生体組織の機能を客観的かつ定量的に評価することは、非常に有用である。
【０００４】
最近注目を集めている治療法として、重度心不全患者に対する心臓再同期療法（ＣＲＴ）が挙げられる。この心臓再同期療法の適用可否の判定や治療の効果判定のための定量的評価に対しても超音波診断装置を用いることが試みられている。
【０００５】
心臓再同期療法について簡単に説明する。重度心不全患者は、心壁運動の収縮同期不全を併発することが多い。心臓は、電気信号の伝導によって動いている。重度心疾患の患者には、心室内伝導障害が発生することがある。心室内伝導障害は、心筋を運動させる電気信号の伝わる順序にずれが起きることである。心室内伝導障害は、このずれにより本来ほぼ同時に全体に電気信号が伝わる心室において、早く伝わる部分と遅く伝わる部分ができてしまうことがある。この結果、心壁の収縮が同期せず、血液を充分に送り出せなくなり、心不全の状態となる。
【０００６】
心臓再同期療法は、このような障害に対して人工的に電気信号を出して心臓に伝わる電気信号の順序を整え、心臓のポンプ機能を助ける治療である。具体的には、ペースメーカーを胸の皮膚の下に埋め込むことによりなされる。心臓再同期療法は、既に多数の患者に施されており、劇的な症状の改善が確認されている。
【０００７】
一方、心臓再同期療法を適応しても症状の改善が見られない症例がある。心不全症例である。心不全症例の患者は、全体の約３割確認されている。これは、心不全症例の原因が収縮同期不全であるか否かを正確に判断できないためである。
【０００８】
現状、心臓再同期療法の適応基準は、心電波形のＱＲＳ幅が１３０ｍｓｅｃ未満、左室駆出率（ＥＦ）が３５％であることとされている。この基準では、心不全であるが収縮同期不全ではない患者も含まれてしまう。
【０００９】
そこで、超音波診断装置を用いた定量的評価法により収縮同期不全のみを抽出しようという技術が開発されている。この技術は、例えば特許文献１に開示されている。特許文献１は、心筋（心臓壁）の運動速度をドプラー法により検出し、この運動速度を演算・解析することを開示する。この技術は、心筋の複数部位における運動速度や変位などの経時変化のピークを自動的に検出できる。そして、この技術は、所定心時相から各ピークに到達するまでの時間を算出し、この時間に応じて心筋の超音波画像に色付ける。これにより、心筋全体の運動状態の違いがカラー画像として出力される。心筋の各部位における運動タイミングの違いが画像化できる。
【００１０】
超音波診断装置は、生体内の構造だけでなく組織の移動速度を組織ドプラー法により映像化することが実現されている。組織ドプラー法は、心臓の壁などの反射が強く比較的動きの遅い部位の速度を計測して二次元カラー表示する。近年では、組織の移動速度を組織ドプラー法により映像化する機能を応用して、さらに臨床上重要な情報を得る手法が提案されている。例えば、組織の移動速度分布像を用いた移動タイミングの影像法と、Velocity Gradient影像法と、Displacement影像法と、Strain影像法と、角度補正を用いたティッシュトラッキング法とである。
【００１１】
組織の移動速度分布像を用いた移動タイミングの影像法は、組織ドプラー法によって得られた組織の移動速度分布像を用い、各画素の組織移動速度がある閾値に達する時間やピークとなる時間をカラーコーディングする。これにより、組織の移動タイミングを映像化する、この影像法は、例えば特許文献２に開示されている。
【００１２】
Velocity Gradient影像法は、組織ドプラー法で得られた移動速度分布に対して一定距離離れた２点の速度差をこれら２点間の距離で除算し、局所のVelocity Gradientを取得する。そして、局所のVelocity Gradientの取得を画像の多数の点で行い、画像として表示する。
【００１３】
Displacement影像法は、組織ドブラ法で得られた移動速度分布に対して一定時間に含まれる各フレームの速度値にフレーム間時間差を乗算した値を積分し、組織の移動量を計算する。
【００１４】
Strain影像法は、Displacement影像法で得られた移動量分布像若しくはVelocity Gradient像を用いて組織の歪(Strain)分布像を得る。
【００１５】
Displacement像若しくはStrain像を得る際には、複数フレームの値を用いた積分処理が必要となる。通常考えられるように、複数フレームの同一画素の値を積分したのでは、移動している組織に対して同一組織の値を積算したことにならない。角度補正を用いたティッシュトラッキング法は、組織の移動を追いかけながら積算する必要がある。
【００１６】
しかしながら、上記移動タイミングの影像法では、以下の問題を生じる。  
第１に、組織は動いている。このため、各フレームの同一画素の速度値を用いたのでは、同一組織の速度の最大時刻を表示しない。例えば、心臓の短軸像を考える。図１１は例えばモニターに表示された被検体Ｓとしての例えば心臓の左室短軸像の模式図を示す。ＵＳ１は心臓の左室短軸像の外膜を表す。ＵＳ２は心臓の拡張末期における左室短軸像の内膜を表す。このＵＳ３は、心臓の収縮末期における左室短軸像の内膜を表す。
【００１７】
この心臓の短軸像に点Ａを設定する。左室の心筋が短軸像で動いた場合、点Ａは拡張末期において心腔内にある。点Ａは収縮末期において心筋内である。図１２は左室の心筋が短軸像で動いたときの収縮運動の速度変化を示す。点Ａの速度の時間的変化は、後半において心筋の速度を反映している。これに対して点Ａの速度の時間的変化は、前半において心腔内の速度となる。このため、心筋の移動速度のピークは正しく検出されない。
【００１８】
第２に、心筋の複数部位における運動速度や変位などの経時変化のピークを自動的に検出する。図１３は時間経過に対する心筋のある部位の速度の変化を示す。この速度の変化は、頂上付近で平坦を示す。このような速度の変化が頂上付近で平坦な場合、ちょっとした変動でピーク位置の時刻が大きく動いてしまうことがある。
【００１９】
第３に、移動速度の変化の仕方が臨床的に有用である可能性もある。心筋が短時間で急激に収縮したのか、比較的長い時間でゆっくり収縮したのかを、ピーク時刻からは分からない。
【００２０】
第４に、速度は、周囲の心筋の影響を受けてしまう。速度は、心筋の移動をみる。速度は、心筋のその部分が収縮したかを見ているわけではない。例えば、心尖像の場合、心筋局所の移動速度は、心尖部からその部分に至る部分全体の収縮によって決まる。心筋局所の移動速度は、その部分の収縮の強さは表さない。よって、局所の収縮タイミングを表す指標としては、速度を用いることは必ずしも適切でない場合がある。
【００２１】
第５に、速度のピークは、最もActiveな時刻を表し完結時刻は表さない。例えば心筋が最も活発に収縮している時刻はわかる。ところが、いつ収縮期が終了したかはわからない。
【００２２】
従来、Velocity Gradient影像法、Displacement影像法、Strain影像法、又は角度補正を用いたティッシュトラッキング法のいずれかを用いる場合がある。これら場合についても、精度の高い、組織の移動若しくは変形に関するタイミングを映像化できるようにしたものが強く望まれている。
【特許文献１】特開平１０−２６２９７０号公報
【特許文献２】特許第３１８７００８号公報
【発明の開示】
【発明が解決しようとする課題】
【００２３】
本発明の目的は、精度の高い、組織の移動若しくは変形に関するタイミングを映像化し得る超音波診断装置及びその制御プログラムを提供することにある。
【課題を解決するための手段】
【００２４】
本発明の第１の局面に係る超音波診断装置は、被検体内に超音波を送信し、被検体からの反射波を受信して受信信号を出力するスキャン部と、スキャン部から出力される受信信号に基づいて被検体の動きの速度の時間的変化を求め、この速度の時間的変化から被検体における動きのタイミングの違いを表示する演算表示部とを具備する。
【００２５】
本発明の第２の局面に係る超音波診断装置は、被検体内に超音波を送信し、被検体からの反射波を受信して受信信号を出力するスキャン部と、スキャン部から出力される受信信号に基づいて被検体の動きの速度を検出する検出部と、検出部により検出された被検体の動きの速度によりstrain値、velocity gradient値又はdisplacement値の少なくとも１つの値を求め、この値の変化から被検体における動きのタイミングの違いを表示する演算表示部とを具備する。
【発明の効果】
【００２６】
本発明によれば、精度の高い、組織の移動若しくは変形に関するタイミングを映像化し得る超音波診断装置及びその制御プログラムを提供できる。
【発明を実施するための最良の形態】
【００２７】
以下、本発明の第１の実施の形態について図面を参照して説明する。
【００２８】
図１は超音波診断装置のブロック構成図を示す。超音波プローブ１１には、送信部１２と受信部１３とが接続されている。超音波プローブ１１は、複数の圧電振動子と、各圧電振動子に設けられる整合層と、各圧縮振動子から後方への超音波の伝播を防止するバッキング材等とを有する。各圧電振動子は、送信部１２からの駆動信号に基づいて超音波パルスを発生し、被検体Ｓからの反射波を受信して電気信号に変換する。従って、超音波プローブ１１は、例えば被検体Ｓ内に超音波パルスを送信する。被検体Ｓ内に送信された超音波パルスは、被検体Ｓ内の体内組織や血流等の音響インピーダンスの不連続面で次々と反射し、反射波として超音波プローブ１１に入射する。超音波プローブ１１は、反射波を受信して電気信号に変換出力する。
【００２９】
送信部１２は、超音波プローブ１１の各圧縮振動子に駆動信号を送出し、超音波パルスを発生させる。送信部１２は、各圧縮振動子に送出する駆動信号の位相を予め設定された一定範囲で制御する。送信部１２は、各圧縮振動子の駆動を位相制御し、送信する超音波パルスをスキャンさせる。
【００３０】
受信部１３は、超音波プローブ１１からの電気信号を増幅して出力する。
【００３１】
ビーム形成部１４は、受信部１３から出力される電気信号をデジタル変換し、このデジタル電気信号を所定の時間だけ遅延し、この後、整相加算し、フォーカシングを行う。ビーム形成部１４は、フォーカシングした信号をＢ／Ｍ（Ｂモード／Ｍモード）処理部１５とＣＦＭ（カラードプラーモード）処理部１６とに供給する。
【００３２】
Ｂ／Ｍ処理部１５は、ビーム形成部１４の出力信号に対してバンドパスフィルタ処理を施し、この後、その包絡線成分を検出し、ＬＯＧ圧縮処理を行う。Ｂ／Ｍ処理部１５は、エッジ強調等の処理を行ってもよい。
【００３３】
ＣＦＭ処理部１６は、ビーム形成部１４の出力信号に対してハイパスフィルタ処理を行い、この後、自己相関処理を行う。ハイパスフィルタ処理は、組織信号と血流信号とを分離する。ハイパスフィルタ処理は、例えばＭＴＩフィルタ若しくはＷａｌｌフィルタを有する。本発明は、組織ドプラー法を前提とする。この場合、ＭＴＩフィルタは、全周波数通過若しくはローパスフィルタとなる。これにより、ＭＴＩフィルタは、組織信号を通過させる。その他、ハイパスフィルタ処理は、組織信号を低減・削除するための非線形処理を行う場合もある。自己相関処理は、血流や組織の移動速度を検出する。
【００３４】
第１のスキャンコンバータ１７は、Ｂ／Ｍ処理部１５の出力信号を超音波パルスの送受信に対応した位置にマッピングし、１フレームずつ複数の画像データを逐次作成する。
【００３５】
第２のスキャンコンバータ１８は、ＣＦＭ処理部１６の出力信号を超音波パルスの送受信に対応した位置にマッピングし、１フレームずつ複数の画像データを逐次作成する。
【００３６】
ビデオインタフェース部１９は、第１のスキャンコンバータ１７から出力されるＢ／Ｍ用の画像データ、又は第２のスキャンコンバータ１８から出力されるＣＦＭ用の画像データを入力し、これら画像データと画像に関する様々な情報とを組み合わせ、レイアウトした超音波画像をモニター２０に表示する。
【００３７】
制御部２１は、送信部１２と、受信部１３と、ビーム形成部１４と、Ｂ／Ｍ処理部１５と、ＣＦＭ処理部１６と、第１のスキャンコンバータ１７と、第２のスキャンコンバータ１８と、ビデオインタフェース部１９とを制御する。制御部２１は、Ｂ／Ｍ処理部１５の出力信号、ＣＦＭ処理部１６の出力信号、第１と第２のスキャンコンバータ１７、１８からそれぞれ出力される各画像データを取り込む。
【００３８】
制御部２１は、第１のスキャンコンバータ１７と第２のスキャンコンバータ１８とにより逐次作成される画像データ上における画素を制御点として設定し、後続して作成される各フレームの画像データ上における制御点の移動を追跡し、各フレームの画像データ間の相関処理を行って制御点の速度を求め、この速度の変化から被検体Ｓにおける動きのタイミングの違いをモニター２０に表示する。
【００３９】
制御部２１は、制御点を画像データにおける全ての画素又は一部の画素に設定する。
【００４０】
制御部２１は、演算表示部としての第１の処理部２１ａと、第２の処理部２１ｂと、カラー表示部２１ｃとを有する。第１の処理部２１ａは、時間経過に対する被検体Ｓの組織の動きの速度の変化の各ピーク位置を求める。
【００４１】
第２の処理部２１ｂは、第１の処理部２１ａにより求められた各ピーク位置に基づいて被検体Ｓの収縮運動による組織の動きのタイミングの違いを求める。
【００４２】
第２の処理部２１ｂは、第１のスキャンコンバータ１７と第２のスキャンコンバータ１８とにより逐次作成される各画像データから超音波画像をモニター２０に表示する。第２の処理部２１ｂは、モニター２０に表示する超音波画像に対して被検体Ｓの収縮運動による組織の動きのタイミングの違いを例えば並べて表示する。これにより、第２の処理部２１ｂは、被検体Ｓの組織の動きの急峻、緩慢を表示する。
【００４３】
カラー表示部２１ｃは、第２の処理部２１ｂにより求められた被検体Ｓの組織の動きのタイミングの違いをモニター２０にカラー表示する。
【００４４】
制御部２１は、Ｂ／Ｍ処理部１５及びＣＦＭ処理部１６の各出力信号を入力し、スキャンコンバーションを行って被検体Ｓの組織の動きのタイミングの違いの解析を行ってもよい。
【００４５】
なお、制御部２１には、記憶部２２が接続されている。第１のスキャンコンバータ１７と第２のスキャンコンバータ１８とにより逐次作成される画像データは、パーソナルコンピュータ２３に送られている。
【００４６】
次に、上記の如く構成された装置の動作について説明する。
【００４７】
先ず、被検体Ｓの診断処理が開始されると、送信部１２は、超音波プローブ１１の各圧縮振動子に駆動信号を送出する。このとき、送信部１２は、各圧縮振動子に送出する駆動信号の位相を予め設定された一定範囲で制御する。これにより、超音波プローブ１１は、例えば被検体Ｓ内に超音波パルスをスキャンして送信する。被検体Ｓ内に送信された超音波パルスは、被検体Ｓ内の体内組織や血流等の音響インピーダンスの不連続面で次々と反射し、反射波として超音波プローブ１１に入射する。超音波プローブ１１は、反射波を受信して電気信号に変換出力する。受信部１３は、超音波プローブ１１からの電気信号を増幅して出力する。
【００４８】
ビーム形成部１４は、受信部１３から出力される電気信号をデジタル変換し、このデジタル電気信号を所定の時間だけ遅延し、この後、整相加算し、フォーカシングを行う。
【００４９】
Ｂ／Ｍ処理部１５は、ビーム形成部１４の出力信号に対してバンドパスフィルタ処理を施し、この後、その包絡線成分を検出し、ＬＯＧ圧縮処理を行う。
【００５０】
これと共にＣＦＭ処理部１６は、ビーム形成部１４の出力信号に対して例えばＭＴＩフィルタによりハイパスフィルタ処理を行って組織信号と血流信号とを分離する。次に、ＣＦＭ処理部１６は、組織信号と血流信号とに対して自己相関処理を行って血流や組織の移動速度を検出する。
【００５１】
第１のスキャンコンバータ１７は、Ｂ／Ｍ処理部１５の出力信号を超音波パルスの送受信に対応した位置にマッピングし、１フレームずつ複数の画像データを逐次作成する。
【００５２】
これと共に第２のスキャンコンバータ１８は、ＣＦＭ処理部１６の出力信号を超音波パルスの送受信に対応した位置にマッピングし、１フレームずつ複数の画像データを逐次作成する。
【００５３】
ビデオインタフェース部１９は、第１のスキャンコンバータ１７から出力されるＢ／Ｍ用の画像データ、又は第２のスキャンコンバータ１８から出力されるＣＦＭ用の画像データを入力し、これら画像データと画像に関する様々な情報とを組み合わせ、レイアウトした被検体Ｓ内の超音波画像をモニター２０に表示する。
【００５４】
図１１はモニター２０に表示された被検体Ｓとしての例えば心臓の左室短軸像の模式図を示す。ＵＳ１は心臓の左室短軸像の外膜を表す。ＵＳ２は心臓の拡張末期における左室短軸像の内膜を表す。このＵＳ３は、心臓の収縮末期における左室短軸像の内膜を表す。この心臓の短軸像に点Ａを設定する。左室の心筋が短軸像で動いた場合、点Ａは拡張末期において心腔内にある。点Ａは収縮末期において心筋内である。図１２は左室の心筋が短軸像で動いたときの収縮運動の速度変化を示す。点Ａの速度の時間的変化は、後半において心筋の速度を反映している。これに対して点Ａの速度の時間的変化は、前半において心腔内の速度となる。このため、心筋の移動速度のピークは正しく検出されない。
【００５５】
これに対して本実施の形態の制御部２１は、図２に示す制御処理フローチャートを実行する。
【００５６】
制御部２１は、第１と第２のスキャンコンバータ１７、１８により作成される各フレームの画像データを取り込んで例えば記憶部２２に一時記憶する。制御部２１は、ステップ＃１において、第１と第２のスキャンコンバータ１７、１８から取り込んだ１フレームの画像データの全ての画素に各制御点を自動的に設定する。なお、制御部２１は、１フレームの画像データの少なくとも一部を構成する各画素に制御点を設定してもよい。但し、これ以下、説明の都合上により２つの制御点についてのみ説明する。
【００５７】
図３はフレームＦ１の画像データ上に設定した２つの制御点Ｃ１、Ｃ２を示す。これら制御点Ｃ１、Ｃ２は、心臓の収縮運動による組織の移動する速度変化を求めるために設定される。
【００５８】
次に、制御部２１は、ステップ＃２において、フレームＦ１よりも後続して作成されるフレームＦ２の画像データ上における各制御点Ｃ１、Ｃ２の移動を追跡する。各制御点Ｃ１、Ｃ２は、各フレームＦ１、Ｆ２の各画像データ上における心臓の同一組織を追跡する。追跡は、組織ドプラー法により得られる速度情報を用いて行われる。これにより、各フレームＦ１とＦ２の各画像データ間における各制御点Ｃ１、Ｃ２の移動が追跡される。
【００５９】
制御部２１は、各フレームＦ１とＦ２の各画像データ上における各制御点Ｃ１、Ｃ２の各移動距離と各フレームＦ１とＦ２の各画像データの各間隔とに基づいて各制御点Ｃ１、Ｃ２の各速度値を求める。なお、各フレーム間の移動距離は、次式により求められる。
【００６０】
次のフレームまでの移動距離＝速度値×フレーム間時間（一定）    …（１）
制御部２１は、各制御点Ｃ１、Ｃ２の速度値を例えば記憶部２２に記憶する。
【００６１】
次に、制御部２１は、ステップ＃３において、全フレームＦ１〜Ｆｎの画像データ上における各制御点Ｃ１、Ｃ２の速度値を求めることを終了したか否かを判断する。この判断の結果、終了していなければ、制御部２１は、ステップ＃４からステップ＃２に戻り、フレームＦ２と次のフレームＦ３の画像データ上における各制御点Ｃ１、Ｃ２を追跡する。これにより、制御部２１は、各フレームＦ２とＦ３の各画像データ上における各制御点Ｃ１、Ｃ２の各移動距離と各フレームＦ１〜Ｆｎの各画像データの各間隔とに基づいて各制御点Ｃ１、Ｃ２の各速度値を求める。制御部２１は、各制御点Ｃ１、Ｃ２の速度値を例えば記憶部２２に記憶する。
【００６２】
以下、同様に、制御部２１は、ステップ＃２からステップ＃４を繰り返し、全フレームＦ１〜Ｆｎの各画像データ間における各制御点Ｃ１、Ｃ２を追跡し、これら制御点Ｃ１、Ｃ２の速度値を求める。制御部２１は、各制御点Ｃ１、Ｃ２の各速度値を例えば記憶部２２に記憶する。なお、記憶部２２には、全フレームＦ１〜Ｆｎの各画像データにおける全ての画素における各制御点の各速度値が記憶される。
【００６３】
全フレームＦ１〜Ｆｎの各画像データ間における各制御点Ｃ１、Ｃ２の速度値が求められると、制御部２１は、ステップ＃５において、記憶部２２に記憶した全フレームＦ１〜Ｆｎの各画像データ間の各制御点Ｃ１、Ｃ２の各速度値を読み出し、各制御点Ｃ１、Ｃ２の各速度の時間的変化を求める。
【００６４】
図４は時間経過に対する各制御点Ｃ１、Ｃ２の各速度の各変化曲線を示す。これら変化曲線は、心臓に各制御点Ｃ１、Ｃ２おける動きの速度の変化を示す。同図から心臓の収縮運動において制御点Ｃ１にて観測された組織よりも制御点Ｃ２にて観測された組織の方が遅れて動いていることが分かる。
【００６５】
制御点の変化曲線は、全フレームＦ１〜Ｆｎの各画像データにおける全ての画素において求められている。これにより、被検体Ｓである心臓全体における各制御点の各変化曲線が求められる。
【００６６】
第１の処理部２１ａは、２つの制御点Ｃ１、Ｃ２の各変化曲線における各ピーク位置を求める。図４は各制御点Ｃ１、Ｃ２の各変化曲線における各ピーク位置Ｐ１、Ｐ２を示す。同図１３各ピーク位置Ｐ１、Ｐ２が異なる。
【００６７】
第２の処理部２１ｂは、第１の処理部２１ａにより求められた各ピーク位置に基づいて被検体Ｓの収縮運動による組織の動きのタイミングの違いを求める。第２の処理部２１ｂは、図４に示すような各制御点Ｃ１、Ｃ２の各変化曲線における各ピーク位置Ｐ１、Ｐ２を
モニター２０に表示してもよい。各ピーク位置Ｐ１、Ｐ２が異なることから各制御点Ｃ１、Ｃ２における被検体Ｓの収縮運動による組織の動きのタイミングが異なることが分かる。心臓の収縮運動による組織の動きのタイミングが異なるのは、心壁運動の収縮同期不全であると診断できる。
【００６８】
なお、図５は正常な心臓における２つの制御点Ｃ１、Ｃ２の各変化曲線における各ピーク位置Ｐ１、Ｐ２を示す。同図１３各ピーク位置Ｐ１、Ｐ２が略同一である。これにより、各制御点Ｃ１、Ｃ２における被検体Ｓの収縮運動による組織の動きのタイミングが略同一であることが分かる。
【００６９】
次に、制御部２１は、ステップ＃６において、画像データの全ての画素において求められた各制御点Ｃ１、Ｃ２を含む全変化曲線に基づいて例えば各制御点Ｃ１、Ｃ２の移動の画像データを作成する。
【００７０】
次に、カラー表示部２１ｃは、ステップ＃７において、画像データの全ての画素において求められた各制御点Ｃ１、Ｃ２を含む全変化曲線に基づいて被検体Ｓである心臓の収縮運動による組織の動きのタイミングの違いをカラーコーディングしてモニター２０に表示する。
【００７１】
具体的にカラー表示部２１ｃは、例えば図６に示すように各制御点Ｃ１、Ｃ２を含む全変化曲線に対してカラーバーＣＢを設定する。カラーバーＣＢは、例えば時間経過が遅くなると共に例えば緑色から赤色に変化する配色になっている。カラーバーＣＢの配色は、緑色から赤色に限らず、他の色を用いてもよい。カラーバーＣＢは、各制御点Ｃ１、Ｃ２を含む全変化曲線の各ピーク位置Ｐ１、Ｐ２を含む長さに設定されている。すなわち、カラーバーＣＢは、心臓の収縮運動による組織の動きのタイミングの相違により生じる各ピーク位置Ｐ１、Ｐ２の違いを区分するために設定される。
【００７２】
カラー表示部２１ｃは、各ピーク位置Ｐ１、Ｐ２に対応する色をカラーバーＣＢから求め、これら色を各制御点Ｃ１、Ｃ２のカラーとしてコーディングする。例えば、制御点Ｃ１の変化曲線のピーク位置Ｐ１は、緑色にコーディングする。制御点Ｃ２の変化曲線のピーク位置Ｐ２は、赤色にコーディングする。このようにして画像データの全ての画素において求められた各制御点の各変化曲線のピーク位置の色をコーディングする。
【００７３】
次に、制御部２１は、ステップ＃８において、各制御点Ｃ１、Ｃ２の移動の画像データ、カラーコーディングの画像データをビデオインタフェース部１９に送る。このビデオインタフェース部１９は、第１のスキャンコンバータ１７から出力されるＢ／Ｍ用の画像データ、又は第２のスキャンコンバータ１８から出力されるＣＦＭ用の画像データを入力し、これら画像データと各制御点Ｃ１、Ｃ２の移動の画像データ、カラーコーディングの画像データ、さらには画像に関する様々な情報とを組み合わせ、レイアウトした被検体Ｓ内の超音波画像をモニター２０に表示する。
【００７４】
図７はモニター２０に表示された超音波画像の一例を示す。この超音波画像は、心臓の収縮運動を動画で表示する。この超音波画像は、心臓の収縮末期の画像を示す。全制御点のうち例えば制御点Ｃ１の変化曲線のピーク位置Ｐ１は、緑色にコーディングされて表示される。制御点Ｃ２の変化曲線のピーク位置Ｐ２は、赤色にコーディングされて表示される。従って、心臓の超音波画像において緑色と赤色とがコーディングされて表示されている。コーディングされた色が緑色と赤色とで異なることは、心臓の収縮運動による組織の動きのタイミングが異なることを示す。この結果、モニター２０に表示された超音波画像のコーディングの色から心壁運動の収縮同期不全であることが診断できる。
【００７５】
一方、図８は正常な心臓における各制御点Ｃ１、Ｃ２の各変化曲線に設定したカラーバーＣＢを示す。各ピーク位置Ｐ１、Ｐ２は、略同一位置となっている。これにより、全制御点のうち例えば制御点Ｃ１の変化曲線のピーク位置Ｐ１は、緑色にコーディングされて表示される。制御点Ｃ２の変化曲線のピーク位置Ｐ２も緑色にコーディングされて表示される。従って、モニター２０に表示された超音波画像は、例えば図９に示すように緑色一色でコーディングされて表示される。この結果、モニター２０に表示された超音波画像のコーディングの色から正常な心臓であることが診断できる。
【００７６】
なお、上記処理において、組織ドプラー法で得られた速度像をそのまま用いてもよいが、角度補正を最初に行い、角度補正された速度像を用いるようにすれば正確な結果が得られる。
【００７７】
このように上記第１の実施の形態によれば、画像データ上に例えば制御点Ｃ１、Ｃ２を設定し、全フレームＦ１〜Ｆｎの各画像データ間における各制御点Ｃ１、Ｃ２を追跡し、これら制御点Ｃ１、Ｃ２の速度値を求め、これら各制御点Ｃ１、Ｃ２の各変化曲線を求める。そして、これら制御点Ｃ１、Ｃ２の各変化曲線の時間経過に対する各ピーク位置Ｐ１、Ｐ２に対応する色をコーディングする。これにより、例えば心臓の収縮運動による組織の動きのタイミングの違いが色の違いにより表示される。
【００７８】
この結果、正常な心臓であれば、心臓の全体が略同一タイミングで移動するので、超音波画像は、例えば緑色一色でコーディングされて表示される。これに対して心臓が心壁運動の収縮同期不全であれば、心臓の収縮運動による組織の動きのタイミングに違いが生じる。これにより、超音波画像は、例えば緑色と赤色とでコーディングされて表示される。従って、心臓の収縮運動による組織の動きのタイミングの違いが視覚的に認識できる。例えば収縮同期不全が診断できる。
【００７９】
なお、第１の処理部２１ａは、図４に示すように例えば２つの制御点Ｃ１、Ｃ２の各変化曲線における各ピーク位置Ｐ１、Ｐ２を求めているが、これに限らず、各変化曲線の重心位置又は標準偏差の位置を求めてもよい。
【００８０】
次に、本発明の第２の実施の形態について図面を参照して説明する。なお、上記第１の実施の形態と同一部分の説明は省略し、相違するところを説明する。
【００８１】
制御部２１は、第１のスキャンコンバータ１７と第２のスキャンコンバータ１８とにより逐次作成される画像データ上における全画素を制御点として設定する。制御部２１は、後続して作成される各フレームの画像データ上における各制御点から被検体Sの動きの速度を検出する。
【００８２】
制御部２１は、被検体Ｓの動きの速度によりstrain値、velocity gradient値又はdisplacement値の少なくとも１つの値を求める。制御部２１は、strain値、velocity gradient値又はdisplacement値の少なくとも１つの値の時間的変化から被検体Ｓにおける動きのタイミングの違いをモニター２０にカラー表示する。
【００８３】
次に、上記の如く構成された装置の動作について図１０に示す制御処理フローチャートに従って説明する。
【００８４】
制御部２１は、第１と第２のスキャンコンバータ１７、１８により作成される各フレームの画像データを取り込んで例えば記憶部２２に一時記憶する。制御部２１は、ステップ＃１０において、第１と第２のスキャンコンバータ１７、１８から取り込んだ１フレームの画像データの全画素に各制御点を設定する。制御点は、画像データの少なくとも一部を構成する各画素に設定してもよい。但し、ここでも説明の都合上、２つの制御点についてのみ説明する。
【００８５】
制御部２１は、ステップ＃１１において、フレームＦ１よりも後続して作成されるフレームＦ２の画像データ上における２つの制御点Ｃ１、Ｃ２の移動を追跡する。制御部２１は、各フレームＦ１とＦ２の各画像データ上における各制御点Ｃ１、Ｃ２を追跡しながら各制御点Ｃ１、Ｃ２の移動速度を検出する。制御部２１は、各制御点Ｃ１、Ｃ２の移動速度から各strain値を求める。strain値は、各制御点Ｃ１、Ｃ２間の距離が例えば各フレームＦ１、Ｆ２間でどれだけ変化したかを示す。
【００８６】
従って、フレームＦ１の画像データ上の制御点Ｃ１とフレームＦ２の画像データ上の制御点Ｃ１と間のstrain値が求められる。これと共にフレームＦ１の画像データ上の制御点Ｃ２とフレームＦ２の画像データ上の制御点Ｃ２と間のstrain値が求められる。制御部２１は、各制御点Ｃ１、Ｃ２の各strain値を例えば記憶部２２に記憶する。
【００８７】
次に、制御部２１は、ステップ＃１２において、全フレームＦ１〜Ｆｎの画像データ上における各制御点Ｃ１、Ｃ２の各strain値を求めることを終了したか否かを判断する。この判断の結果、終了していなければ、制御部２１は、ステップ＃１３からステップ＃１１に戻り、フレームＦ２と次のフレームＦ３の画像データ上における各制御点Ｃ１、Ｃ２を追跡する。これにより、制御部２１は、各フレームＦ２とＦ３の各画像データ上における各制御点Ｃ１、Ｃ２の各strain値を求める。
【００８８】
以下、同様に、制御部２１は、ステップ＃１１からステップ＃１３を繰り返し、全フレームＦ１〜Ｆｎの各画像データ間における各制御点Ｃ１、Ｃ２を追跡し、これら制御点Ｃ１、Ｃ２間の各strain値を求める。
【００８９】
全フレームＦ１〜Ｆｎの各画像データ間における各制御点Ｃ１、Ｃ２間の各strain値が求められると、制御部２１は、ステップ＃１４において、記憶部２２に記憶した全フレームＦ１〜Ｆｎの各画像データ間の各制御点Ｃ１、Ｃ２間の各strain値を読み出し、各制御点Ｃ１、Ｃ２における各strain値の時間的変化を示す各変化曲線を求める。制御点の変化曲線は、全フレームＦ１〜Ｆｎの各画像データにおける全ての画素において求められているので、心臓全体における制御点の変化曲線が求められる。
【００９０】
次に、制御部２１は、ステップ＃１５において、画像データの全ての画素において求められた各制御点Ｃ１、Ｃ２を含む全変化曲線に基づいて例えば各制御点Ｃ１、Ｃ２の移動の画像データを作成する。
【００９１】
次に、第１の処理部２１ａは、各制御点Ｃ１、Ｃ２を含む全制御点の各変化曲線における各ピーク位置を求める。
【００９２】
次に、第２の処理部２１ｂは、ステップ＃１６において、画像データの全ての画素において求められた各制御点Ｃ１、Ｃ２を含む全変化曲線に基づいて心臓等の被検体Ｓの収縮運動による組織の動きのタイミングの違いをモニター２０に表示する。
【００９３】
具体的にカラー表示部２１ｃは、例えば上記図６に示すのと同様に、各制御点Ｃ１、Ｃ２を含む全変化曲線に対してカラーバーＣＢを設定する。
【００９４】
カラー表示部２１ｃは、各ピーク位置Ｐ１、Ｐ２に対応する色をカラーバーＣＢから求め、これら色を各制御点Ｃ１、Ｃ２のカラーとしてコーディングする。例えば、制御点Ｃ１の変化曲線のピーク位置Ｐ１は、緑色にコーディングする。制御点Ｃ２の変化曲線のピーク位置Ｐ２は、赤色にコーディングする。このようにして画像データの全ての画素において求められた各制御点の各変化曲線のピーク位置の色をコーディングする。
【００９５】
次に、制御部２１は、ステップ＃１７において、各制御点Ｃ１、Ｃ２の移動の画像データ、カラーコーディングの画像データをビデオインタフェース部１９に送る。このビデオインタフェース部１９は、第１のスキャンコンバータ１７から出力されるＢ／Ｍ用の画像データ、又は第２のスキャンコンバータ１８から出力されるＣＦＭ用の画像データを入力し、これら画像データと各制御点Ｃ１、Ｃ２の移動の画像データ、カラーコーディングの画像データ、さらには画像に関する様々な情報とを組み合わせ、レイアウトした被検体Ｓ内の超音波画像をモニター２０に表示する。
【００９６】
この結果、被検体Ｓが心壁運動の収縮同期不全であれば、上記図７に示す超音波画像と同様に、モニター２０に表示される超音波画像は、緑色と赤色とがコーディングされる。従って、コーディングされた色が緑色と赤色とで異なることから心壁運動の収縮同期不全であることが認識できる。
【００９７】
一方、被検体Ｓが正常な心臓であれば、モニター２０に表示される超音波画像は、例えば図９に示すように緑色一色でコーディングされて表示される。これにより、正常な心臓であることが認識できる。
【００９８】
このように上記第２の実施の形態によれば、画像データ上に例えば制御点Ｃ１、Ｃ２を設定し、全フレームＦ１〜Ｆｎの各画像データ間における各制御点Ｃ１、Ｃ２を追跡し、各フレームの各制御点Ｃ１、Ｃ２における各strain値を求め、各制御点Ｃ１、Ｃ２における各strain値の各変化曲線を求める。これにより、上記第１の実施の形態と同様の効果を奏することができる。
【００９９】
本発明は、上記各実施の形態に限定されるものでなく次の通り変形してもよい。
【０１００】
例えば、パーソナルコンピュータ２３は、第１及び第２のスキャンコンバータ１７、１８から出力される各画像データを取り込む。パーソナルコンピュータ２３は、各画像データを処理して例えば被検体Ｓの組織の運動を解析してもよい。
【０１０１】
制御部２１は、２つの制御点Ｃ１、Ｃ２の速度値又はstrain値に限らず、velocity gradient値又はdisplacement値の少なくとも１つの値を求めてもよい。Velocity Gradient値は、組織ドプラー法で得られた移動速度分布に対して一定距離離れた２点の速度差をこれら２点間の距離で除算して取得する。Velocity Gradient値の取得は、各画像データの全画素で行う。Displacement値は、組織ドブラ法で得られた移動速度分布に対して一定時間に含まれる各フレームの速度値にフレーム間時間差を乗算した値を積分し、組織の移動量として求める。Displacement値は、被検体Ｓ内における組織が任意の時間内にどれだけ移動したかを示す。制御部２１は、全ての制御点における各velocity gradient値又は各displacement値の時間的変化を示す各変化曲線を求める。
【０１０２】
上記各実施の形態は、制御点の速度値、strain値、velocity gradient値又はdisplacement値を求めて記憶部２２に記憶し、この後に自動的に速度値、strain値、velocity gradient値又はdisplacement値の変化曲線を求めている。これに限らず、制御部２１は、制御点の速度値、strain値、velocity gradient値又はdisplacement値を一旦記憶部２２に記憶する。この後、制御部２１は、例えばユーザの指示の入力を待つ。ユーザの指示が制御部２１に入力されると、制御部２１は、速度値、strain値、velocity gradient値又はdisplacement値の変化曲線を求めるようにしてもよい。
【０１０３】
その他、超音波診断装置の構成、制御部の機能、タイミング画像生成手順及びその内容、制御点の変化量の種類等についても、本発明の要旨を逸脱しない範囲で種々変形して実施できる。
【図面の簡単な説明】
【０１０４】
【図１】本発明に係わる超音波診断装置の第１の実施の形態を示すブロック構成図。
【図２】同装置における制御処理フローチャート。
【図３】同装置において画像データ上に設定した２つの制御点の移動例を示す図。
【図４】同装置により求められる収縮同期不全の心臓における時間経過に対する各制御点の各変化曲線における各ピーク位置を示す図。
【図５】同装置により求められる正常な心臓における時間経過に対する各制御点の各変化曲線における各ピーク位置を示す図。
【図６】同装置による収縮同期不全の心臓に対するカラーバーを用いたカラーコーディングを示す図。
【図７】同装置のモニターに表示された収縮同期不全の心臓のカラーコーディングされた超音波画像の一例を示す図。
【図８】同装置による正常な心臓に対するカラーバーを用いたカラーコーディングを示す図。
【図９】同装置のモニターに表示された正常な心臓のカラーコーディングされた超音波画像の一例を示す図。
【図１０】本発明に係わる超音波診断装置の第２の実施の形態における制御処理フローチャート。
【図１１】心臓の左室短軸像を示すモニター表示例の模式図。
【図１２】左室の心筋が短軸像で動いたときの収縮運動の速度変化を示す図。
【図１３】時間経過に対する心筋のある部位の速度の変化を示す図。
【符号の説明】
【０１０５】
１１：超音波プローブ、１２：送信部、１３：受信部、１４：ビーム形成部、１５：Ｂ／Ｍ（Ｂモード／Ｍモード）処理部、１６：ＣＦＭ（カラードプラーモード）処理部、１７：第１のスキャンコンバータ、１８：第２のスキャンコンバータ、１９：ビデオインタフェース部、２０：モニター、２１：制御部、２１ａ：第１の処理部、２１ｂ：第２の処理部、２１ｃ：カラー表示部、２２：記憶部、２３：パーソナルコンピュータ。
【図１】
2006289067.tif 000003
【図２】
2006289067.tif 000004
【図３】
2006289067.tif 000005
【図４】
2006289067.tif 000006
【図５】
2006289067.tif 000007
【図６】
2006289067.tif 000008
【図７】
2006289067.tif 000009
【図８】
2006289067.tif 000010
【図９】
2006289067.tif 000011
【図１０】
2006289067.tif 000012
【図１１】
2006289067.tif 000013
【図１２】
2006289067.tif 000014
【図１３】
2006289067.tif 000015