(57)【特許請求の範囲】
【請求項１】
被検者に向けて超音波を送信すると共に、被検者から反射された超音波信号を受信することにより受信信号を出力する複数の素子が配列された超音波プローブと、
前記被検者に向けて送信する超音波の実音速に対して予め設定された仮定音速を変化させる手段と、
１回の送信から２種類以上の異なる前記仮定音速を変化させて生成されるＲＦデータ又は振幅画像のスペックルを利用して各フレーム間における各部分の類似度を求め、各フレーム間相互の位置を算出するフレーム位置算出手段と、
を有することを特徴とする超音波診断装置。
【請求項２】
前記フレーム位置算出手段は、前記各フレーム間相互の位置として、前記各部分のフレーム内での移動量及びフレームに垂直な移動量を算出することを特徴とする請求項１に記載の超音波診断装置。
【請求項３】
前記ＲＦデータ又は振幅画像は、前記素子の配列方向の分解能が該素子間隔以上のデータから生成されることを特徴とする請求項１または２に記載の超音波診断装置。
【請求項４】
複数の素子が配列された超音波プローブから被検者に向けて超音波を送信すると共に、被検者から反射された超音波信号を受信し、
前記被検者に向けて送信する超音波の実音速に対して予め設定された仮定音速を変化させ、
１回の送信から２種類以上の異なる前記仮定音速を変化させて生成されるＲＦデータ又は振幅画像のスペックルを利用して各フレーム間における各部分の類似度を求め、各フレーム間相互の位置を算出するようにしたことを特徴とする超音波診断方法。
【請求項５】
前記各フレーム間相互の位置として、前記各部分のフレーム内での移動量及びフレームに垂直な移動量を算出することを特徴とする請求項４に記載の超音波診断方法。
【請求項６】
前記ＲＦデータ又は振幅画像は、前記素子の配列方向の分解能が該素子間隔以上のデータから生成されることを特徴とする請求項４または５記載の超音波診断方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波診断方法及び装置に係り、特に、各フレーム間における各部分の類似度を安定して均一に求める超音波診断方法及び装置に関する。
【背景技術】
【０００２】
従来より、超音波を用いて被検者の断層画像を取得して医療診断に供することが行われているが、例えば、心臓などの臓器の動きを計測するような場合に、各フレーム間において、あるフレーム中のある部分が次のフレームのどの部分に相当するかを求めその位置合わせをするのに「類似度」が用いられている。ここで各フレーム間における各部分の類似度を求める方法として、従来、相互相関法、正規化相互相関法、ＳＡＤ法あるいはＳＳＤ法などの技術が知られており、この技術はパノラマ画像合成、３Ｄ画像合成、動き・回転補正、スペックルや組織のトラッキングなど数多くの分野で使用されている。
【０００３】
例えば、フレーム間でのテクスチャの変化に着目してセンサレスでプローブ位置をトラッキングするものが知られている（例えば、特許文献１等参照）。これは、超音波断層画像に現れるスペックルパターンの、プローブの変位に対する変化が、観察している組織内部の状態に依存しないことに基づき、スペックルパターンの変化を用いて超音波断層画像の３次元的な位置を特定するものである。
【０００４】
すなわち、これにおいては、超音波プローブの移動前後の画像の相関関数のピーク値から断層画像面における移動ベクトルを測定し、さらに移動前後の画像のうちどちらか一方の画像から小さい領域を複数切り出して、その領域と他方画像の相関関数のピーク位置とピーク値とから回転角を測定するようにしている。
【特許文献１】特開２００２−１０２２２３号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
上述したような、従来のフレーム間でのテクスチャの類似性を利用する技術は、一般的に利用できるテクスチャの種類が多い程安定し、テクスチャの分解能が高い程精度向上が望めるが、その反面分解能が高い程フレーム間での変化が大きく、それによる精度低下も大きいという問題がある。
【０００６】
また、画像のある部分のテクスチャの種類を増やすには、異なる周波数やビームの送受信により得られる複数ＲＦデータ又は振幅画像を利用することで可能であるが、余分な送信が必要となってしまう。また、信号を分割することや各送受信間でプローブ位置と反射体の位置関係が変わってしまうことにより、フレームレートや分解能が低下し、精度が犠牲となる。
【０００７】
テクスチャの分解能は、画像生成時の仮定音速と実音速のずれが大きい程低下するが、局所局所で仮定音速と実音速のずれ量が異なることから、部分部分でテクスチャ、スペックルの分解能が異なり、その結果、類似度の精度の劣化や部分部分で精度が異なるという均一性の劣化が生じてしまうという問題があった。
【０００８】
さらに、上記従来技術のようにフリーハンドでの３Ｄ画像構築において、類似度の値からプローブ移動距離を推定する処理が行われるが、スペックルの分解能によってその関係が異なるため推定精度が悪化してしまうという問題があった。
【０００９】
本発明は、このような事情に鑑みてなされたもので、フレームレートを低下させずに安定して、均一、高精度に類似度を求めることのできる超音波診断方法及び装置を提供することを目的とする。
【課題を解決するための手段】
【００１０】
前記目的を達成するために、請求項１に記載の発明は、被検者に向けて超音波を送信すると共に、被検者から反射された超音波信号を受信することにより受信信号を出力する複数の素子が配列された超音波プローブと、前記被検者に向けて送信する超音波の実音速に対して予め設定された仮定音速を変化させる手段と、１回の送信から２種類以上の異なる前記仮定音速を変化させて生成されるＲＦデータ又は振幅画像のスペックルを利用して各フレーム間における各部分の類似度を求め、各フレーム間相互の位置を算出するフレーム位置算出手段と、を有することを特徴とする超音波診断装置を提供する。
【００１１】
これにより、仮定音速によって異なる複数のスペックルパターンを利用するため、類似度を安定して求めることができる。また、分解能の低いパターンから分解能の高いパターンまで均一に含むため、安定して高精度に類似度を求めることができる。さらに、各部分で最適音速（実音速）が異なるにも関わらず、同じように分解能の低いスペックルから高いスペックルまで含むため、どの部分でも均一に類似度を求めることができる。また、同じ送信から得られた受信データから生成することで、フレームレートを低下させずに、また、送信間のずれによる精度低下を防ぐことができる。
【００１２】
また、請求項２に示すように、前記フレーム位置算出手段は、前記各フレーム間相互の位置として、前記各部分のフレーム内での移動量及びフレームに垂直な移動量を算出することを特徴とする。
【００１３】
これにより、算出された各フレームの各部分の断層面に平行な移動量、及び断層面に垂直な移動量に従って各フレームを設置することにより３Ｄ画像を生成することができる。
【００１４】
また、請求項３に示すように、前記ＲＦデータ又は振幅画像は、前記素子の配列方向の分解能が該素子間隔以上のデータから生成されることを特徴とする。
【００１５】
これにより、類似度を高精度に求めることができる。
【００１６】
また、同様に前記目的を達成するために、請求項４に記載の発明は、複数の素子が配列された超音波プローブから被検者に向けて超音波を送信すると共に、被検者から反射された超音波信号を受信し、前記被検者に向けて送信する超音波の実音速に対して予め設定された仮定音速を変化させ、１回の送信から２種類以上の異なる前記仮定音速を変化させて生成されるＲＦデータ又は振幅画像のスペックルを利用して各フレーム間における各部分の類似度を求め、各フレーム間相互の位置を算出するようにしたことを特徴とする超音波診断方法を提供する。
【００１７】
これにより、仮定音速によって異なる複数のスペックルパターンを利用するため、類似度を安定して求めることができる。また、分解能の低いパターンから分解能の高いパターンまで均一に含むため、安定して高精度に類似度を求めることができる。さらに、各部分で最適音速（実音速）が異なるにも関わらず、同じように分解能の低いスペックルから高いスペックルまで含むため、どの部分でも均一に類似度を求めることができる。また、同じ送信から得られた受信データから生成することで、フレームレートを低下させずに、また、送信間のずれによる精度低下を防ぐことができる。
【００１８】
また、請求項５に示すように、前記各フレーム間相互の位置として、前記各部分のフレーム内での移動量及びフレームに垂直な移動量を算出することを特徴とする。
【００１９】
これにより、算出された各フレームの各部分の断層面に平行な移動量、及び断層面に垂直な移動量に従って各フレームを設置することにより３Ｄ画像を生成することができる。
【００２０】
また、請求項６に示すように、前記ＲＦデータ又は振幅画像は、前記素子の配列方向の分解能が該素子間隔以上のデータから生成されることを特徴とする。
【００２１】
これにより、類似度を高精度に求めることができる。
【発明の効果】
【００２２】
以上説明したように、本発明によれば、仮定音速によって異なる複数のスペックルパターンを利用するため、類似度を安定して求めることができる。また、分解能の低いパターンから分解能の高いパターンまで均一に含むため、安定して高精度に類似度を求めることができる。さらに、各部分で最適音速（実音速）が異なるにも関わらず、同じように分解能の低いスペックルから高いスペックルまで含むため、どの部分でも均一に類似度を求めることができる。
【発明を実施するための最良の形態】
【００２３】
以下、添付図面を参照して、本発明に係る超音波診断方法及び装置について詳細に説明する。
【００２４】
本発明は、異なる仮定音速で生成されたＲＦデータ又は振幅画像のスペックルを利用して、フレームレートを低下させずに安定して均一に、高精度に類似度を求め、フレーム位置を算出するものである。
【００２５】
なお、被検者に対して送信する実際の超音波の音速（実音速）に対して、所定量ずつステップ刻みで複数変化させて設定する超音波音速を設定音速あるいは仮定音速と言う。また、実音速は最適音速とも言う。
【００２６】
本発明は、以下の事に着目してなされたものである。すなわち、まず第１に、異なる仮定音速で生成されたスペックルパターンが異なるため、それぞれで類似度を求めることより、安定性が向上すること。第２に、異なる仮定音速で生成されたスペックルパターンは単に異なるというだけではなく、分解能の低いスペックルパターンから分解能が高いスペックルパターンまで均一に含むため、類似度の算出において安定性が向上し、高精度であること。第３に、仮定音速と最適音速とのずれは局所局所で異なるにも関わらず、複数の仮定音速を利用することで、どの部分でもある程度同じように分解能が低いパターンから分解能が高いパターンまで均一に含むため、部分によらず均一に類似度の算出ができ、精度が向上すること。第４に、１回の同じ送信から得られた受信データから生成できるので、余分な送信を必要とせず、送信間のずれが無いこと、である。
【００２７】
以下、フリーハンドでの３Ｄ画像生成を例にとって説明する。
【００２８】
図１は、本発明に係る超音波診断装置の一実施形態の概略構成を示すシステム構成図である。
【００２９】
図１に示すように、本実施形態の超音波診断装置１は、超音波を用いて被検者の診断部位について超音波画像を撮影して表示するものであり、超音波プローブ１０、送受信部１２、走査制御部１４、ＡＤ変換部１６、画像生成部１８、フレーム位置算出部２０、表示画像生成部２２及びモニタ２４を有して構成されている。
【００３０】
超音波プローブ１０は、被検者の体内の診断部位に向けて超音波を送信するとともに体内で反射してきた超音波を受信するものである。本実施形態の超音波プローブ１０は、１次元の超音波トランスデューサアレイを構成する複数の超音波トランスデューサを備えており、各超音波トランスデューサは、例えばＰＺＴ等の圧電素子の両端に電極を形成した振動子によって構成されている。この電極は信号線によって送受信部１２と接続されている。各電極に電圧を印加すると振動子は超音波を発生する。また、振動子は反射してきた超音波を受信すると電気信号を発生し、これが受信信号として出力される。
【００３１】
送受信部１２は、超音波プローブ１０に超音波送信信号を与え振動子から超音波を発生させ、走査制御部１４から与えられた遅延に基づいて送波する。そして、反射した超音波を受信して超音波プローブ１０が出力した各素子の受信信号をそのまま（受波フォーカスをかけず）増幅する。
【００３２】
ＡＤ変換部１６は、送受信部１２から超音波受信信号を受け取りＡＤ変換して画像生成部１８に引き渡す。画像生成部１８はＡＤ変換部１６から受け取った受信データを保存する。画像生成部１８では、保存された各素子の受信データから、詳しくは後述するが、様々に設定される音速（これを上述したように被検者に送波する実際の音速（実音速）に対して仮定音速という。）に基づく遅延で受波フォーカスされ、各仮定音速に基づくＲＦデータが生成される。
【００３３】
フレーム位置算出部２０は、複数の仮定音速を用いてフレームの位置を算出するものである。また、表示画像生成部２２は、画像生成部１８で生成された画像と、フレーム位置算出部２０で算出されたフレーム位置の結果からモニタ２４に表示するための３Ｄ画像を生成するものである。
【００３４】
次に、図１の装置構成における画像生成部１８の作用を図２のフローチャートに沿って説明する。
【００３５】
画像生成部１８は、仮定音速を変化させていろいろな音速で得られたデータから画像を生成するものである。
【００３６】
まず図２のステップＳ１００において、いろいろ変化させる仮定音速の初期値を設定する。この値は特に限定されるものではなく、例えば、１４００［ｍ／ｓ］のように適宜決めればよい。
【００３７】
そして設定された初期値により、走査制御部１４によって制御された送受信部１２から超音波プローブ１０に信号が送られ、その仮定音速初期値によるデータが取得され画像生成部１８に送られる。
【００３８】
次にステップＳ１１０において、仮定音速を所定量１ステップ変更し、変更された仮定音速による超音波データが取得される。この１ステップの所定量は、特に限定されず、例えば２０［ｍ／ｓ］でもよいし、１０［ｍ／ｓ］でも、４０［ｍ／ｓ］でもよく、所定量だけ仮定音速を変化させていく。
【００３９】
次にステップＳ１２０において、得られた各仮定音速によるデータを位相を整合して加算し、ＲＦ（Radio Frequency)データを生成する。このＲＦデータは、振幅情報と位相情報の両方を含むものである。このようにすべての仮定音速での画像でＲＦデータを作成する。
【００４０】
そしてステップＳ１３０において、画像生成が終了したか否か判断し、まだ終了していない場合にはステップＳ１１０に戻り、また仮定音速を１ステップ変更し画像生成を続行する。画像生成の終了は、すべての仮定音速についての処理が終了したか否かで判断する。それは例えば、仮定音速を何ステップ変更したら終了するかを予め決めておき、その回数をカウントして判断するようにすればよい。
【００４１】
次に、フレーム位置算出部２０の作用を説明する。
【００４２】
図３は、フレーム位置算出部２０における処理の流れを示すフローチャートである。
【００４３】
まず、ステップＳ２００において、フレーム１の標準画像についてスペックル領域の探索を行い、検出した各スペックル領域で所定サイズのカーネルを設定する。
【００４４】
なお、ここでスペックル領域の探索に利用される標準画像とは、超音波音速１５４０［ｍ／ｓｅｃ］による画像である。
【００４５】
スペックル領域の探索におけるスペックルか否かを判断する方法は、特には限定されず、例えば、レイリー分布からの逸脱度を調べる周知の方法でもよい。この方法は、超音波画像に現れるスペックルパターンは無数の散乱体が超音波の解像度以下の細かさで分布しているときに散乱体の無数の重畳によって超音波信号の強度に高い部分と低い部分とが生じる現象であり、この散乱体がランダムに分布している場合には、それらの散乱体から反射される超音波信号の強度である振幅置の確率密度分布は、Ｐ（ｘ）＝（ｘ／ｓ２）ｅｘｐ（−ｘ２／２ｓ２）（ただし、ｓ２は分散を表し、平均０と規格化されている。）で表されるレイリー分布に従うのであるが、組織内にある種の構造物が増加していくとそのスペックルパターンはその構造物を反映するようになりランダムとは言えなくなって、その結果、輝度の確率密度関数はレイリー分布から逸脱するようになるので、これによりスペックルを判定するものである。
【００４６】
また、受信した超音波画像からＲＦデータを生成する際に仮定音速を変化させたとき、スペックルの場合には、仮定音速に依らず位相はランダムになるという位相変化特性を利用してスペックルを判定するようにしてもよい。
【００４７】
あるいは、ユーザ自身がスペックル領域を指定するようにしても良い。
【００４８】
次に、ステップＳ２１０において、フレームＮｏ．ｎを初期化して、ｎ＝１とする。また、ステップＳ２２０において、仮定音速を初期化する。これはすでに画像生成部１８における処理で得られているデータを用いればよい。
【００４９】
次に、ステップＳ２３０において、仮定音速を１ステップ変更し、その音速でのデータを取得するが、これも画像生成部１８における処理で得られているデータを使用すればよい。
【００５０】
次に、ステップＳ２４０において、フレームｎの各カーネルにつき、次のフレームｎ＋１上での類似度ピークを探索し、その類似度のピーク位置とそのピーク値とを求める。
【００５１】
類似度を算出する方法は特に限定されず、２つの画像データのフレーム間における相互相関係数として、０〜１間の値として算出する相互相関法や、ＳＡＤ法あるいはＳＳＤ法などの周知の方法で算出することができる。類似度の算出に用いられるデータは、振幅画像とＲＦデータのどちらでも良い。ここで、ＲＦデータとは振幅情報と位相情報の両方を含むデータである。
【００５２】
類似度とは、画像間において変化がない部分であるか否かを示す値であり、類似度が高いということは、画像中のその部分が時間的に変化していないことを示す。すなわち、類似度を求めることにより、あるフレーム中のある部分が次のフレーム中のどの部分に相当するかがわかる。このように類似度を求めることにより、フレームの位置を算出することが可能となる。
【００５３】
次に、ステップＳ２５０において、各カーネルにつき、類似度のピーク位置から断層面に平行な移動ベクトル（断層面内での移動量）を求め、類似度のピーク値から断層面に垂直な移動距離（断層面に垂直な移動量）を求める。
【００５４】
ここで、ピーク値と垂直移動量との関係は予め測定（計算）してテーブル化しておくと良い。
【００５５】
次に、ステップＳ２６０において、以上の処理が全ての仮定音速について終了したか否か判断する。そして、まだ終了していない場合にはステップＳ２３０へもどり、仮定音速を１ステップ変更して次の仮定音速についての処理を行う。
【００５６】
また、全ての仮定音速について処理が終了した場合には、次のステップＳ２７０において、全仮定音速の各カーネルでの断層面に平行な移動ベクトル及び断層面に垂直な移動距離を平均化する。
【００５７】
また、全仮定音速について各カーネルの類似度を平均化してもよい。すなわち、類似度を平均化して、平均化された類似度、類似度ピークからその断層面に平行な移動ベクトル、断層面に垂直な移動距離を求めるようにしてもよい。
【００５８】
次に、ステップＳ２８０において、フレームＮｏ．ｎを１インクリメントしてｎ＋１とする。
【００５９】
そして、ステップＳ２９０において、全てのフレームについての処理が終了したか否か判断する。
【００６０】
まだ、全てのフレームについての処理が終了していない場合には、ステップＳ２２０にもどり、次のフレームｎ＋１についての処理を行う。一方、全てのフレームについての処理が終了した場合には、フレーム位置算出部２０における処理を終了する。
【００６１】
次に、表示画像生成部２２の処理について説明する。
【００６２】
図４は、表示画像生成部２２における処理内容を示すフローチャートである。
【００６３】
まず、図４のステップＳ３００において、各フレームの各カーネル位置で求められた断層面に平行な移動量、及び断層面に垂直な移動量に従って、各フレームを設置することにより３Ｄ画像を生成する。
【００６４】
次に、ステップＳ３１０において、生成された３Ｄ画像に対し一般的な対数圧縮を行い、さらに、ゲイン／ＤＲ（ダイナミックレンジ）／ＳＴＣ（深さ重み付け）／グレーマップ調整等を行う。
【００６５】
そして表示画像生成部２２で生成された表示画像はモニタ２４に表示される。
【００６６】
以上説明したように、本実施形態においては、異なる複数の仮定音速で生成されたＲＦデータ又は振幅画像を使用して各フレームの各部分の類似度を求めるようにしている。
【００６７】
このように、仮定音速によって異なる複数のスペックルパターンを利用するため、安定して類似度を求めることができ、特に分解能の低いパターンから高いパターンまで均一に含むため、安定してしかも高精度に求めることが可能となる。
【００６８】
また、各部分で最適音速が異なるにも関わらず、同じ様に分解能の低いスペックルから高いスペックルまで含むので、どの部分でも均一に求めることができる。さらに、同じ１回の送信から得られた受信データから生成することにより、フレームレートを低下させず、また送信間のずれによる精度低下を防ぐことができる。さらに、フリーハンドでの３Ｄ画像構築で類似度の値とプローブ移動距離との関係を利用する場合には、その関係がスペックルの分解能に依るので、分解能が低いパターンから分解能が高いパターンまでの部分においてもそれを用いる関係がどの部分でも同じであるため、分解能が異なることによる推定精度の低下を防ぐことができる。
【００６９】
なお、フレーム位置算出部２０における処理は上述したものには限らない。以下、フレーム位置算出部２０における他の処理の例を説明する。
【００７０】
まず図５のフローチャートは、フレーム位置算出部２０における他の処理例その１を示すものである。前述した図３のフローチャートに示した処理では、初めのフレームでスペックル位置を設定するようにしていたが、この例は、各フレーム毎においてスペックル位置を設定するようにしたものである。
【００７１】
図５のステップＳ４００において、フレームＮｏ．ｎを初期化し、ｎ＝１とする。
【００７２】
次にステップＳ４１０において、フレームｎの標準画像につきスペックル領域の探索を行い、検出した各スペックル領域で所定サイズのカーネルを設定する。なお、標準画像とは超音波音速１５４０［ｍ／ｓｅｃ］の画像である。
【００７３】
次にステップＳ４２０において、仮定音速を初期化し、次のステップＳ４３０で仮定音速を１ステップ変更する。
【００７４】
そしてステップＳ４４０において、フレームｎの各カーネルについて次のフレームｎ＋１上での類似度ピークを探索する。
【００７５】
次のステップＳ４５０において、各カーネルにつき、断層面に平行な移動ベクトル及び断層面に垂直な移動距離を算出する。
【００７６】
そしてステップＳ４６０で全仮定音速についての処理が終了したか否か判断し、まだ終了していない場合には、ステップＳ４３０にもどり、仮定音速を１ステップ変更して次の仮定音速についての処理を行う。
【００７７】
また、全仮定音速についての処理が終了した場合には、次のステップＳ４７０において、全仮定音速の各カーネルでの断層面に平行な移動ベクトル及び断層面に垂直な移動距離を平均化する。
【００７８】
そして、次のステップＳ４８０において、フレームＮｏ．ｎを１インクリメントしてｎ＋１とし、ステップＳ４９０で全フレームについての処理が終了したか否か判断する。
【００７９】
まだ全フレームについての処理が終了していない場合には、ステップＳ４１０にもどり、次のフレームｎ＋１についての処理を行うこととし、全フレームについての処理が終了した場合には、フレーム位置算出部２０の処理を終了する。
【００８０】
次に図６のフローチャートは、フレーム位置算出部２０における他の処理例その２を示すものである。この例は、初めのフレームでスペックル初期位置を設定し、各フレームでのスペックル位置を前フレームでのスペックル位置から移動ベクトルで移動した位置に設定するようにするものである。
【００８１】
まず図６のステップＳ５００において、フレーム１の標準画像につきスペックル領域の探索を行い、検出した各スペックル領域で所定サイズのカーネルを設定する。なお、標準画像とは超音波音速１５４０［ｍ／ｓｅｃ］の画像である。
【００８２】
次にステップＳ５１０において、フレームＮｏ．ｎを初期化してｎ＝１とする。また、ステップＳ５２０において、仮定音速を初期化し、次のステップＳ５３０において、仮定音速を１ステップ変更する。
【００８３】
ステップＳ５４０において、フレームｎの各カーネルにつき、次のフレームｎ＋１上での類似度ピークを探索する。
【００８４】
ステップＳ５５０において、各カーネルにつき、断層面に平行な移動ベクトル及び断層面に垂直な移動距離を算出する。
【００８５】
そして、次のステップＳ５６０において、全ての仮定音速についての処理が終了したか否か判断し、まだ全仮定音速についての処理が終了していない場合には、ステップＳ５３０に戻り、仮定音速を１ステップ変更して次の仮定音速についての処理を行う。
【００８６】
また、全仮定音速についての処理が終了した場合には、次のステップＳ５７０に進み、全仮定音速の各カーネルでの断層面に平行な移動ベクトル及び断層面に垂直な移動距離を平均化する。
【００８７】
ステップＳ５８０において、フレームｎにおける各カーネル位置に、各カーネルの移動ベクトルを加えたものを、次のフレームｎ＋１における各カーネル位置とする。
【００８８】
すなわち、（次のフレームｎ＋１における各カーネル位置）＝（フレームｎにおける各カーネル位置）＋（各カーネルの移動ベクトル）とする。
【００８９】
そして、次のステップＳ５９０において、フレームＮｏ．ｎを１インクリメントして、ｎ＋１とし、次のステップＳ５９５において、全てのフレームについての処理が終了したか否か判断し、まだ全フレームについての処理が終了していない場合には、ステップＳ５２０に戻り、次のフレームｎ＋１についての処理を行い、全フレームについての処理が終了した場合にはフレーム位置算出部２０における処理を終了する。
【００９０】
このようにフレーム位置算出部２０における処理の流れは一通りに限定されるものではなく、いくつかの方法が有り得る。
【００９１】
最近のソフトウエアベースの超音波装置やアナログベースでも高性能な回路構成によって、同じ１つの送信から得られた受信信号から種々の仮定音速で画像生成することが可能となってきている。本発明における装置構成は、種々の仮定音速でのＲＦデータ又は画像をフレームずれ無しに得るため、またフレームレートを低下させないために必要な構成を有している。
【００９２】
従来の超音波装置は、素子間隔の音線位置で送受信ビームを形成し、ＲＦデータを生成し、振幅画像を生成した後に素子間の音線につき、補間して振幅データを生成し表示画像とする構成が一般的であったが、最近では、素子間の音線についても、送受信ビーム形成し、ＲＦデータ生成する構成が見られる。
【００９３】
本実施形態の装置構成においては、例えばフレーム位置算出部においてスキャン方向に分解能が素子間隔以上のデータを利用するように、スキャン方向に高分解能なデータを利用することで、スペックルのランダムな位相変化をより正確に判別し、スキャン方向に素子間隔以上の分解能で高精度に類似度を算出することが可能である。
【００９４】
以上述べた実施形態には、超音波の送受信周波数が１種類のＲＦデータを利用する場合のみを挙げたが、基本波と高調波など、複数の異なる周波数のＲＦデータを利用する場合も本発明に含まれる。
【００９５】
以上、本発明の超音波診断方法及び装置について詳細に説明したが、本発明は、以上の例には限定されず、本発明の要旨を逸脱しない範囲において、各種の改良や変形を行ってもよいのはもちろんである。
【図面の簡単な説明】
【００９６】
【図１】本発明に係る超音波診断装置の一実施形態の概略構成を示すシステム構成図である。
【図２】画像生成部の作用を示すフローチャートである。
【図３】フレーム位置算出部における処理の流れを示すフローチャートである。
【図４】表示画像生成部における処理内容を示すフローチャートである。
【図５】フレーム位置算出部における他の処理例その１を示すフローチャートである。
【図６】フレーム位置算出部における他の処理例その２を示すフローチャートである。
【符号の説明】
【００９７】
１…超音波診断装置、１０…超音波プローブ、１２…送受信部、１４…走査制御部、１６…ＡＤ変換部、１８…画像生成部、２０…フレーム位置算出部、２２…表示画像生成部、２４…モニタ
【図１】
0005108444.tif 000002
【図２】
0005108444.tif 000003
【図３】
0005108444.tif 000004
【図４】
0005108444.tif 000005
【図５】
0005108444.tif 000006
【図６】
0005108444.tif 000007