(57)【特許請求の範囲】
【請求項１】
超音波プローブの直下に位置する生体組織に変形を加えながら、異なる時間に同一方向に送信した超音波の２つ以上の反射波に基づく受信信号を用いて、超音波送信方向にて所定間隔で設定された各ウィンドウ幅にある受信信号間の自己相関関数を計算し、当該計算した自己相関関数を用いて位相差を計算することで、前記各ウィンドウ幅における変位を計算する変位計算部と、
前記自己相関関数を計算した時と略同一のウィンドウ幅にて各受信信号の振幅値の積の重心位置を計算し、当該計算した重心位置を前記変位計算部により計算された前記変位を与える位置として補正する補正部と、
前記各ウィンドウ幅の変位及び重心位置から傾きを計算することで、前記生体組織の歪みを計算する歪み計算部と、
前記歪み計算部が計算した前記歪みを用いてエラストグラフィーを生成する画像生成部と、
を備えたことを特徴とする超音波診断装置。
【請求項２】
前記補正部は、更に、各ウィンドウ幅の中心位置の変位を、他のウィンドウ幅にて計算された重心位置の変位を用いて補間することで計算し、
前記歪み計算部は、前記各ウィンドウ幅の補間された変位及び中心位置から傾きを計算することで、歪みを計算することを特徴とする請求項１に記載の超音波診断装置。
【請求項３】
超音波プローブの直下に位置する生体組織に変形を加えながら、異なる時間に同一方向に送信した超音波の２つ以上の反射波に基づく受信信号を用いて、超音波送信方向にて所定間隔で設定された各ウィンドウ幅にある受信信号間の相互相関係数がピークを取るずれ量を計算するずれ量計算部と、
前記ずれ量計算部が前記ずれ量の取得時に用いた相互相関係数を計算した時と略同一のウィンドウ幅にて各受信信号の振幅値の積の重心位置である第１の重心位置を前記ずれ量を補正した補正ずれ量として計算する第１の補正部と、
前記補正ずれ量によりずらした受信信号間の相互相関関数の位相差から微小変位を計算し、当該算出した微小変位に前記補正ずれ量を加算することで、前記各ウィンドウ幅における変位を計算する変位計算部と、
前記微小変位を計算した時と略同一のウィンドウ幅にて各受信信号の振幅値の積の重心位置である第２の重心位置を計算し、当該計算した第２の重心位置を前記変位計算部により計算された前記変位を与える位置として補正する第２の補正部と、
前記所定間隔で設定された各ウィンドウ幅の変位及び重心位置から傾きを計算することで、歪みを計算する歪み計算部と、
前記歪み計算部が計算した前記歪みを用いてエラストグラフィーを生成する画像生成部と、
を備えたことを特徴とする超音波診断装置。
【請求項４】
前記補正部は、更に、各ウィンドウ幅の中心位置の変位を、他のウィンドウ幅にて計算された第２の重心位置の変位を用いて補間することで計算し、
前記歪み計算部は、前記各ウィンドウ幅の補間された変位及び中心位置から傾きを計算することで、歪みを計算することを特徴とする請求項３に記載の超音波診断装置。
【請求項５】
超音波プローブの直下に位置する生体組織に変形を加えながら、異なる時間に同一方向に送信した超音波の２つ以上の反射波に基づく受信信号を用いて、超音波送信方向にて所定間隔で設定された各ウィンドウ幅にある受信信号間の自己相関関数を計算し、当該計算した自己相関関数を用いて位相差を計算することで、前記各ウィンドウ幅における変位を計算する変位計算部と、
前記自己相関関数を計算した時と略同一のウィンドウ幅にて各受信信号の振幅値の積の重心位置を計算し、当該計算した重心位置を前記変位計算部により計算された前記変位を与える位置として補正する補正部と、
前記各ウィンドウ幅の変位及び重心位置から傾きを計算することで、前記生体組織の歪みを計算する歪み計算部と、
を備えたことを特徴とする画像処理装置。
【請求項６】
超音波プローブの直下に位置する生体組織に変形を加えながら、異なる時間に同一方向に送信した超音波の２つ以上の反射波に基づく受信信号を用いて、超音波送信方向にて所定間隔で設定された各ウィンドウ幅にある受信信号間の相互相関係数がピークを取るずれ量を計算するずれ量計算部と、
前記ずれ量計算部が前記ずれ量の取得時に用いた相互相関係数を計算した時と略同一のウィンドウ幅にて各受信信号の振幅値の積の重心位置である第１の重心位置を前記ずれ量を補正した補正ずれ量として計算する第１の補正部と、
前記補正ずれ量によりずらした受信信号間の相互相関関数の位相差から微小変位を計算し、当該算出した微小変位に前記補正ずれ量を加算することで、前記各ウィンドウ幅における変位を計算する変位計算部と、
前記微小変位を計算した時と略同一のウィンドウ幅にて各受信信号の振幅値の積の重心位置である第２の重心位置を計算し、当該計算した第２の重心位置を前記変位計算部により計算された前記変位を与える位置として補正する第２の補正部と、
前記所定間隔で設定された各ウィンドウ幅の変位及び重心位置から傾きを計算することで、歪みを計算する歪み計算部と、
を備えたことを特徴とする画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、超音波診断装置及び画像処理装置に関する。
【背景技術】
【０００２】
近年、超音波診断装置において実時間で生体組織の硬さを表示するエラストグラフィーが実用化されている。被検体の体表に当てた状態で超音波プローブを手で軽く押し付けたり離したりすると、超音波プローブの直下にある生体組織は、変形する。エラストグラフィーでは、変形の前後、あるいは変形の程度の異なる２つ以上の状態から同一部位の位置の変化（変位と呼ばれる）を計測し、変位を空間的に微分することで歪み（ストレインとも呼ばれる）を計算して表示する。硬い生体組織ほど変形しにくいので、硬い生体組織のストレイン値は小さくなり、軟らかい生体組織のストレイン値は大きくなる。
【０００３】
エラストグラフィーにおいて変位を計測する方法には、様々な方法がある。すなわち、変位を計測する方法としては、「自己相関法」や、「相互相関法」、「複合自己相関法」、「ゼロ位相法」などがある。「自己相関法」は、組織ドプラ法を用いて短い時間差での変位を計測する方法である。また、「相互相関法」は、連続する２フレーム間で相互相関係数のピーク位置を計測する方法である。また、「複合自己相関法」は、連続する２フレーム間で相互相関法と自己相関法とを組み合わせた方法である。また、「ゼロ位相法」は、連続する２フレーム間で相互相関関数の位相が「０」となる点を検索する方法である。なお、変位を計測する方法には、「複合自己相関法」を応用した「拡張複合自己相関法」も知られている。
【０００４】
上記の変位を計測する方法は、全て、所定の範囲（ウィンドウ幅と呼ばれる）内の積算値を使用するため、ウィンドウ内でエコー強度の大きい場所があると、エコー強度の大きい場所の重みが大きくなる。しかし、上記の変位を計測する方法では、ウィンドウ内の中心を変位が計測された位置とするために、ウィンドウ内でエコー強度の大きい場所があると、位置には、誤差が生じる。そのため、エラストグラフィーで表示される画像には、均一なストレインを持つ生体組織でも、エコー強度（Ｂモード像の輝度）の変化に対応するかのような模様が見えてしまう場合がある。
【０００５】
このため、従来のエラストグラフィーでは、エコー強度の影響を受けてしまっていた。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００４−５７６５２号公報
【特許文献２】国際公開第２００７／１１０６６９号
【非特許文献】
【０００７】
【非特許文献１】Tsuyoshi  Shiina  ,  Naotaka  Nitta,  Ei  Ueno,  Jeffrey  C.  Bamber,  "Real  time  tissue  elasticity  imaging  using  the  combined  autocorrelation  method",  Journal  of  Medical  Ultrasonics,  Vol.26,  No.2,  pp.57-66,  1999
【発明の概要】
【発明が解決しようとする課題】
【０００８】
本発明が解決しようとする課題は、エコー強度の影響を低減したエラストグラフィーを表示することが可能となる超音波診断装置及び画像処理装置を提供することである。
【課題を解決するための手段】
【０００９】
実施形態の超音波診断装置は、変位計算部と、補正部と、歪み計算部と、画像生成部とを備える。変位計算部は、超音波プローブの直下に位置する生体組織に変形を加えながら、異なる時間に同一方向に送信した超音波の２つ以上の反射波に基づく受信信号を用いて、超音波送信方向にて所定間隔で設定された各ウィンドウ幅にある受信信号間の自己相関関数を計算し、当該計算した自己相関関数を用いて位相差を計算することで、前記各ウィンドウ幅における変位を計算する。補正部は、前記自己相関関数を計算した時と略同一のウィンドウ幅にて各受信信号の振幅値の積の重心位置を計算し、当該計算した重心位置を前記変位計算部により計算された前記変位を与える位置として補正する。歪み計算部は、前記各ウィンドウ幅の変位及び重心位置から傾きを計算することで、前記生体組織の歪みを計算する。画像生成部は、前記歪み計算部が計算した前記歪みを用いてエラストグラフィーを生成する。
【図面の簡単な説明】
【００１０】
【図１】図１は、第１の実施形態に係る超音波診断装置の全体構成を説明するための図である。
【図２】図２は、ウィンドウ幅の一例を説明するための図である。
【図３】図３は、従来のエラストグラフィーにおける課題を説明するための図である。
【図４】図４は、第１の実施形態に係る歪みデータ処理部の構成を説明するための図である。
【図５】図５は、第１の実施形態に係る歪みデータ処理部の概要を説明するための図である。
【図６】図６は、第１の実施形態における補正処理を説明するための図である。
【図７】図７は、補正処理の結果を説明するための図である。
【図８】図８は、第１の実施形態に係る超音波診断装置の歪み計算処理を説明するためのフローチャートである。
【図９】図９は、第２の実施形態に係る歪みデータ処理部の概要を説明するための図である。
【図１０】図１０は、第２の実施形態に係る歪みデータ処理部の概要を説明するための図である。
【図１１】図１１は、第２の実施形態に係る超音波診断装置の歪み計算処理を説明するためのフローチャートである。
【発明を実施するための形態】
【００１１】
以下、添付図面を参照して、超音波診断装置の実施形態を詳細に説明する。
【００１２】
（第１の実施形態）
まず、第１の実施形態に係る超音波診断装置の構成について説明する。図１は、第１の実施形態に係る超音波診断装置の全体構成を説明するための図である。図１に示すように、本実施形態に係る超音波診断装置は、超音波プローブ１と、モニタ２と、入力装置３と、装置本体１０とを有する。
【００１３】
超音波プローブ１は、複数の圧電振動子を有し、これら複数の圧電振動子は、後述する装置本体１０が有する送受信部１１から供給される駆動信号に基づき超音波を発生する。また、超音波プローブ１は、被検体Ｐからの反射波を受信して電気信号に変換する。また、超音波プローブ１は、圧電振動子に設けられる整合層及び音響レンズと、圧電振動子から後方への超音波の伝播を防止するバッキング材などを有する。超音波プローブ１は、装置本体１０と着脱自在に接続される。
【００１４】
超音波プローブ１から被検体Ｐに超音波が送信されると、送信された超音波は、被検体Ｐの体内組織における音響インピーダンスの不連続面で次々と反射され、反射波信号として超音波プローブ１が有する複数の圧電振動子にて受信される。受信される反射波信号の振幅は、超音波が反射される不連続面における音響インピーダンスの差に依存する。なお、送信された超音波パルスが、移動している血流や心臓壁などの表面で反射された場合の反射波信号は、ドプラ効果により、移動体の超音波送信方向に対する速度成分に依存して、周波数偏移（ドプラ偏移）を受ける。
【００１５】
なお、本実施形態は、複数の圧電振動子が一列で配置された１次元超音波プローブである超音波プローブ１により、被検体Ｐを２次元でスキャンする場合であっても、１次元超音波プローブの複数の圧電振動子を機械的に揺動する超音波プローブ１や複数の圧電振動子が格子状に２次元で配置された２次元超音波プローブである超音波プローブ１により、被検体Ｐを３次元でスキャンする場合であっても、適用可能である。
【００１６】
入力装置３は、マウス、キーボード、ボタン、パネルスイッチ、タッチコマンドスクリーン、フットスイッチ、トラックボールなどを有し、超音波診断装置の操作者からの各種設定要求を受け付け、装置本体１０に対して受け付けた各種設定要求を転送する。
【００１７】
モニタ２は、超音波診断装置の操作者が入力装置３を用いて各種設定要求を入力するためのＧＵＩ（Graphical  User  Interface）を表示したり、装置本体１０において生成された超音波画像などを表示したりする。
【００１８】
装置本体１０は、超音波プローブ１が受信した反射波に基づいて超音波画像を生成する装置である。装置本体１０は、図１に示すように、送受信部１１と、位相検波回路１２と、フレームバッファ１３と、Ｂモード処理部１４と、ドプラ処理部１５と、歪みデータ処理部１６と、画像処理部１７と、画像メモリ１８と、制御部１９と、内部記憶部２０とを有する。
【００１９】
送受信部１１は、トリガ発生回路、送信遅延回路およびパルサ回路などを有し、超音波プローブ１に駆動信号を供給する。パルサ回路は、所定の繰り返し周波数（ＰＲＦ：Pulse  Repetition  Frequency）送信超音波を形成するためのレートパルスを繰り返し発生する。なお、ＰＲＦは、レート周波数とも呼ばれる。また、送信遅延回路は、超音波プローブ１から発生される超音波をビーム状に集束して送信指向性を決定するために必要な圧電振動子ごとの送信遅延時間を、パルサ回路が発生する各レートパルスに対し与える。また、トリガ発生回路は、レートパルスに基づくタイミングで、超音波プローブ１に駆動信号（駆動パルス）を印加する。すなわち、送信遅延回路は、各レートパルスに対し与える送信遅延時間を変化させることで、圧電振動子面からの送信方向を任意に調整する。
【００２０】
なお、送受信部１１は、後述する制御部１９の指示に基づいて、所定のスキャンシーケンスを実行するために、送信周波数、送信駆動電圧などを瞬時に変更可能な機能を有している。特に、送信駆動電圧の変更は、瞬間にその値を切り替え可能なリニアアンプ型の発信回路、または、複数の電源ユニットを電気的に切り替える機構によって実現される。
【００２１】
また、送受信部１１は、アンプ回路、Ａ／Ｄ変換器、受信遅延回路、加算器などを有し、超音波プローブ１が受信した反射波信号に対して各種処理を行なって反射波データを生成する。アンプ回路は、反射波信号をチャンネルごとに増幅してゲイン補正処理を行なう。Ａ／Ｄ変換器は、ゲイン補正された反射波信号をＡ／Ｄ変換する。受信遅延回路は、デジタルデータに受信指向性を決定するのに必要な受信遅延時間を与える。加算器は、受信遅延回路により受信遅延時間が与えられたデジタルデータの加算処理を行なう。加算器の加算処理により、反射波信号の受信指向性に応じた方向からの反射成分が強調される。
【００２２】
位相検波回路１２は、送受信部１１が有する加算器の出力信号をベースバンド帯域の同相信号（Ｉ信号、Ｉ：In-pahse）と直交信号（Ｑ信号、Ｑ：Quadrature-phase）とに変換する。そして、位相検波回路１２は、Ｉ信号及びＱ信号（以下、ＩＱ信号と記載する）を反射波データとして後段のフレームバッファ１３に格納する。なお、位相検波回路１２は、加算器の出力信号を、ＲＦ信号や解析信号に変換した上で、フレームバッファ１３に格納することも可能である。以下では、超音波の反射波に基づいて生成された「ＩＱ信号、ＲＦ信号、解析信号」をまとめて、「反射波データ」、又は、「受信信号」と記載する。
【００２３】
フレームバッファ１３は、位相検波回路１２から出力された反射波データを一時的に記憶するバッファである。具体的には、フレームバッファ１３は、数フレーム分の反射波データを記憶する。例えば、フレームバッファ１３は、ＦＩＦＯ（First-In/First-Out）メモリであり、所定フレーム分の反射波データを記憶し、新たに１フレーム分の反射波データが出力された場合、生成時間が最も古い１フレーム分の反射波データを破棄して、新たに出力された１フレーム分の反射波データを記憶する。
【００２４】
なお、１フレーム分の反射波データとは、１枚の超音波画像を生成するための反射波データのことであり、送受信部１１は、複数の走査線（スキャンライン）で形成される走査範囲にて超音波送受信を超音波プローブ１に行なわせることで、位相検波回路１２から１フレーム分の反射波データを出力する。
【００２５】
Ｂモード処理部１４は、フレームバッファ１３から送受信部１１が生成した反射波データを読み出し、読み出した反射波データに対して、対数増幅、包絡線検波処理、対数圧縮などを行なって、信号強度が輝度の明るさで表現されるデータ（Ｂモードデータ）を生成する。
【００２６】
ドプラ処理部１５は、フレームバッファ１３から送受信部１１が生成した反射波データ（ＩＱ信号）を読み出し、読み出した反射波データを周波数解析することで、走査範囲内にある移動体のドプラ効果に基づく運動情報を抽出したデータ（ドプラデータ）を生成する。具体的には、ドプラ処理部１５は、移動体の運動情報として、平均速度、分散値、パワー値などを多点に渡り抽出したドプラデータを生成する。より具体的には、ドプラ処理部１５は、血流の動態を示すカラードプラ画像を生成するためのカラードプラデータや、組織の動態を示す組織ドプラ画像を生成するための組織ドプラデータを生成する。
【００２７】
歪みデータ処理部１６は、反射波データから、生体組織の硬さを表す歪み（ストレイン値）がカラー表示されるエラストグラフィーを生成するためのデータを生成する。具体的には、歪みデータ処理部１６は、被検体Ｐの体表に当てた状態で超音波プローブ１を操作者が手で軽く押し付けたり離したりすることで変形した生体組織の変位を計算する。より具体的には、歪みデータ処理部１６は、変形の前後、あるいは変形の程度の異なる２つ以上の反射波データから同一部位（ウィンドウ幅内）の位置の変位を計算する。そして、歪みデータ処理部１６は、変位を深さ方向に微分することで歪みを計算する。具体的には、歪みデータ処理部１６は、複数のウィンドウ幅それぞれのストレイン値を計算する。図２は、ウィンドウ幅の一例を説明するための図である。
【００２８】
図２に示すように、１本のスキャンライン上には、深さ方向に沿って、複数のサンプル点が設定される。エラストグラフィーでは、複数のサンプル点が一つのウィンドウ幅として設定され、複数のウィンドウ幅が重複して等間隔で設定される。図２に示す一例では、サンプル点の間隔を「λ／２」とすると、「８×λ」がウィンドウ幅（Ｌ）として設定されている。また、図２に示す一例では、深さ方向に沿って、各ウィンドウ幅が「４×λ」分重複した状態で、等間隔で設定されている。このように、ウィンドウ幅は、サンプル点の間隔から、深さ方向の距離に換算して表記することが出来る。なお、図２に示すウィンドウ幅の設定例は、あくまでも、一例であり、ウィンドウ幅を形成するためのサンプル点の数や、ウィンドウ幅の間隔、ウィンドウ幅の重複距離は、任意に設定可能である。
【００２９】
なお、以下では、歪みデータ処理部１６が生成するデータを歪みデータと記載する場合がある。第１の実施形態に係る歪みデータ処理部１６が行なう処理については、後に詳述する。
【００３０】
画像処理部１７は、Ｂモード処理部１４、ドプラ処理部１５及び歪みデータ処理部１６が生成したデータを用いて表示用の超音波画像を生成したり、生成した超音波画像に対して画像処理を行なったりする処理部である。図１に示す画像処理部１７は、画像生成部１７ａ及び画像合成部１７ｂを有する。画像生成部１７ａは、Ｂモード処理部１４、ドプラ処理部１５及び歪みデータ処理部１６が生成したデータから超音波画像を生成する。すなわち、画像生成部１７ａは、Ｂモード処理部１４が生成したＢモードデータから反射波の強度を輝度にて表したＢモード画像を生成する。また、画像生成部１７ａは、ドプラ処理部１５が生成したドプラデータから移動体情報（血流運動情報や組織運動情報）を表す平均速度画像、分散画像、パワー画像、又は、これらの組み合わせ画像としてのドプラ画像（カラードプラ画や組織ドプラ画像）を生成する。
【００３１】
また、画像生成部１７ａは、歪みデータ処理部１６が生成した歪みデータを用いてエラストグラフィーを生成する。具体的には、画像生成部１７ａは、Ｂモード画像上に、各ウィンドウ幅にて計算されたストレイン値の大きさに応じて階調を変化させたストレイン画像を生成し、生成したストレイン画像をＢモード画像に重畳させることで、エラストグラフィーを生成する。
【００３２】
ここで、画像生成部１７ａは、一般的には、超音波走査の走査線信号列を、テレビなどに代表されるビデオフォーマットの走査線信号列に変換（スキャンコンバート）し、表示用画像としての超音波画像を生成する。具体的には、画像生成部１７ａは、超音波プローブ１による超音波の走査形態に応じて座標変換を行なうことで、表示用画像としての超音波画像を生成する。また、画像生成部１７ａは、スキャンコンバート以外に種々の画像処理として、例えば、スキャンコンバート後の複数の画像フレームを用いて、輝度の平均値画像を再生成する画像処理（平滑化処理）や、画像内で微分フィルタを用いる画像処理（エッジ強調処理）などを行なう。
【００３３】
画像合成部１７ｂは、画像生成部１７ａが生成した各種超音波画像に、種々のパラメータの文字情報、目盛り、ボディーマークなどを合成した合成画像を生成する。
【００３４】
画像メモリ１８は、画像生成部１７ａが生成した超音波画像や、画像合成部１７ｂが生成した合成画像を記憶するメモリである。また、画像メモリ１８は、Ｂモード処理部１４やドプラ処理部１５、歪みデータ処理部１６が生成したデータ（生データ）を記憶することも可能である。
【００３５】
内部記憶部２０は、超音波送受信、画像処理及び表示処理を行なうための制御プログラムや、診断情報（例えば、患者ＩＤ、医師の所見など）や、診断プロトコルや各種ボディーマークなどの各種データを記憶する。また、内部記憶部２０は、必要に応じて、画像メモリ１８が記憶する画像の保管などにも使用される。また、内部記憶部２０が記憶するデータは、図示しないインターフェースを経由して、外部の周辺装置へ転送することができる。
【００３６】
制御部１９は、超音波診断装置の処理全体を制御する。具体的には、制御部１９は、入力装置３を介して操作者から入力された各種設定要求や、内部記憶部２０から読込んだ各種制御プログラム及び各種データに基づき、送受信部１１、位相検波回路１２、Ｂモード処理部１４、ドプラ処理部１５、歪みデータ処理部１６、画像処理部１７の処理を制御する。また、制御部１９は、画像メモリ１８が記憶する超音波画像及び合成画像や、操作者が各種処理を指定するためのＧＵＩなどをモニタ２にて表示するように制御する。
【００３７】
以上、第１の実施形態に係る超音波診断装置の全体構成について説明した。かかる構成のもと、第１の実施形態に係る超音波診断装置は、エラストグラフィーの生成及び表示を行なう。より具体的には、第１の実施形態に係る超音波診断装置は、組織ドプラ法を用いて短い時間差での変位を計測する自己相関法によりエラストグラフィーを生成する。
【００３８】
ここで、変位を計算する際には、ウィンドウ幅内の積算値が用いられる。そのために、ウィンドウ内でエコー強度の大きい場所があると、その大きい場所の重みが大きくなる。しかしウィンドウ内の中心を変位の位置とするために、位置に誤差が生じる。そのために、均一なストレインを持つ組織でも、エコー強度つまりＢモード像の輝度の変化に対応するかのような模様が見えてしまう場合がある。
【００３９】
上記の現象を、以下、自己相関法を例に、図３などを用いて詳細に説明する。図３は、従来のエラストグラフィーにおける課題を説明するための図である。なお、自己相関法では、一般的な組織ドプラ法と同様に、１フレームのエラストグラフィーを生成するために、１本のスキャンラインにて２回以上に渡って超音波の送受信が行なわれる。すなわち、自己相関法では、１フレームのスキャンを行なう際に行なわれた同一スキャンラインにおける２つ以上の受信信号を用いて変位を計算する。以下では、２つのＩＱ信号を用いて自己相関法により変位が計算される場合について説明する。
【００４０】
図３の（Ａ）は、深さ方向（x）に「x0＝５ｍｍ」ごとに振幅値が変化する２つのＩＱ信号（IQ1(x)及びIQ2(x)）を示している。図３の（Ａ）に示すグラフでは、縦軸が振幅値（ＩＱ信号の絶対値）であり、横軸が深さ方向（x）となっている。また、図３の（Ａ）に示すグラフでは、IQ1(x)の振幅値を「absIQ1」とし、IQ2(x)の振幅値を「absIQ2」として示している。IQ1(x)を生成するために行なわれた超音波送信と、IQ2(x)を生成するために行なわれた超音波送信との時間差は、小さいため、図３の（Ａ）に示すグラフでは、「absIQ1」と「absIQ2」とが略同じとなっている。なお、図３に示すＩＱ信号では、超音波送信の基本周波数「f0」が「６．６７ＭＨｚ」である。また、図３では、ウィンドウ幅「Ｌ」は、波長を「λ」とすると、「１６×λ＝３．７ｍｍ」となっている。
【００４１】
図３の（Ａ）に示すIQ1(x)及びIQ2(x)は、以下の式（１）で表される。式（１）において、「x」は、深さ方向の距離、「s」は、ストレイン値、「λ」は、超音波の波長、「θ1」及び「θ2」それぞれは、IQ1(x)及びIQ2(x)の初期位相を示す。IQ1(x)及びIQ2(x)の振幅値は、上述したように、同じものとする。また、式（１）では、深さ方向に沿って同一のストレイン値を有する生体組織が存在すると仮定しているために、「s」が「x」に依存しない値となっている。
【００４２】
【数１】
0005509058.tif 000002
【００４３】
図３の（Ｂ）は、式（１）に示す２つのＩＱ信号に対して、自己相関法によって変位を計算した結果を示す。図３の（Ｂ）では、「s」として、４つの値「０．５％、１．０％、１．５％、２．０％」を与えた場合の結果を示している。
【００４４】
本来は、「x」に比例した直線で変位が深さ方向に沿って変化する結果が得られることが期待される。しかし、実際に自己相関法で計算すると、図３の（Ｂ）に示すように、直線にはならずに、振幅値が「０」となる「x」で段差を有する線となってしまう。
【００４５】
図３の（Ｂ）に示す結果をストレイン画像として示したものを図３の（Ｃ）に示す。図３の（Ｃ）では、横方向に４種類のストレイン値「０．５％、１．０％、１．５％、２．０％」それぞれのストレイン画像を示している。なお、図３の（Ｃ）では、下方向が深さ方向であり、下に向かって、「x」が大きくなっている。理想的には、同一のストレイン値を有する生体組織ならば、変位から計算されたストレイン値をマッピングしたストレイン画像は、同じ階調になる。しかし、実際には、図３の（Ｃ）に示すように、同一のストレイン値を有する生体組織にも関わらず、ストレイン画像は、振幅値が「０」となる「x」付近のストレイン値が大きくなっている。つまり、変位から計算されたストレイン値は、エコーの振幅の影響を受けている。
【００４６】
なお、３つ以上のＩＱ信号を用いて自己相関法により変位を計算する場合であっても、上記の現象は発生する。また、自己相関法以外の方法を用いて変位を計算する場合であっても、上記の現象は発生する。
【００４７】
このように、従来のエラストグラフィーでは、エコー強度の影響を受けていた。
【００４８】
そこで、自己相関法によりエラストグラフィーを生成する第１の実施形態に係る超音波診断装置は、エコー強度の影響を低減したエラストグラフィーを表示するために、以下、詳細に説明する歪みデータ処理部１６の処理を実行する。第１の実施形態に係る歪みデータ処理部１６は、自己相関法によって組織ドプラ法と同様な手法を用いて変位を計算し、計算した変位に対して、エコーの振幅による誤差を除去する補正を行なう処理部である。
【００４９】
図４は、第１の実施形態に係る歪みデータ処理部の構成を説明するための図である。図４に示すように、第１の実施形態に係る歪みデータ処理部１６は、変位計算部１６ａと、補正部１６ｂと、歪み計算部１６ｃとを有する。以下、図５を参照して、各部の処理の概要を説明し、その後、数式を用いて、各部の処理を詳細に説明する。図５は、第１の実施形態に係る歪みデータ処理部の概要を説明するための図である。
【００５０】
変位計算部１６ａは、自己相関法により変位を計算する処理部である。変位計算部１６ａは、超音波プローブ１の直下に位置する生体組織に変形を加えながら、異なる時間に同一方向に送信した超音波の２つ以上の反射波に基づく受信信号（ＩＱ信号）を用いて変位を計算する。すなわち、変位計算部１６ａは、１本のスキャンラインにて２回以上に渡って送信された超音波の反射波から生成されたＩＱ信号を用いて変位を計算する。
【００５１】
具体的には、変位計算部１６ａは、図５に示すように、超音波送信方向（深さ方向）にて所定間隔で設定された各ウィンドウ幅にある受信信号間（ＩＱ信号間）の自己相関関数を計算する。そして、変位計算部１６ａは、図５に示すように、計算した自己相関関数を用いて位相差を計算することで、各ウィンドウ幅における変位を計算する。
【００５２】
そして、補正部１６ｂは、図５に示すように、変位補正として、重心位置計算処理と、補間処理とを行なう。すなわち、補正部１６ｂは、自己相関関数を計算した時と略同一のウィンドウ幅にて各受信信号（各ＩＱ信号）の振幅値の積の重心位置を計算し、当該計算した重心位置を変位計算部１６ａにより計算された変位を与える位置として補正する。更に、補正部１６ｂは、各ウィンドウ幅の中心位置の変位を、他のウィンドウ幅にて計算された重心位置の変位を用いて補間することで計算する。
【００５３】
そして、歪み計算部１６ｃは、図５に示すように、各ウィンドウ幅の補間された変位及び中心位置から傾きを計算することで、歪み（ストレイン値）を計算する。
【００５４】
以下では、第１の実施形態に係る歪みデータ処理部１６が２つのＩＱ信号を用いて上記の処理を行なう場合について、数式を用いて詳細に説明する。なお、第１の実施形態に係る歪みデータ処理部１６は、３つ以上のＩＱ信号を用いて上記の処理を行なうことも可能である。
【００５５】
まず、式（１）に示す２つのＩＱ信号（IQ1(x)及びIQ2(x)）を一般化して、以下の式（２）とする。なお、式（２）に示す「A1(x)」及び「A2(x)」それぞれは、IQ1(x)及びIQ2(x)それぞれの振幅値である。また、式（２）に示す「φ1(x)」及び「φ2(x)」それぞれは、IQ1(x)及びIQ2(x)それぞれの位相である。
【００５６】
【数２】
0005509058.tif 000003
【００５７】
ここで、「x」は、上述したように、深さ方向の距離を表わす。また、IQ1(x)とIQ2(x)との時間差を「T」とする。ここで、「T」は、同一地点におけるIQ1(x)とIQ2(x)との変位が半波長（-λ/2〜λ/2）の範囲にあることが保証されるような十分に短い時間（例えば、０．２ｍｓｅｃ〜１ｍｓｅｃの範囲）とする。
【００５８】
理想的な条件では、以下の式（３）に示すような等式が成立する。
【００５９】
【数３】
0005509058.tif 000004
【００６０】
ここで、「s(x)」は、位置「x」におけるストレイン値である。また、上述したように、「λ」は、超音波の波長、「θ1」及び「θ2」それぞれは、IQ1(x)及びIQ2(x)の初期位相を示す。以下では、式（３）に示す理想的な条件で処理が行なわれるが、実際の計算方法を明確化するために式（２）のようにIQ1(x)及びIQ2(x)を記述する。式（１）を用いると、自己相関関数「ac(x)」は、以下の式（４）で計算することが出来る。
【００６１】
【数４】
0005509058.tif 000005
【００６２】
ここで、式（４）に示す「アスタリスク：＊」は、共役複素数を表し、「Ｌ」は、上述したようにウィンドウ幅を表す。
【００６３】
式（４）における複素数「ac(x)」のアークタンジェント（偏角）、すなわち、「angle｛ac(x)｝」がIQ1(x)とIQ2(x)との位相差となり、変位「D(x)」は、以下の式（５）により計算される。
【００６４】
【数５】
0005509058.tif 000006
【００６５】
すなわち、変位計算部１６ａは、式（４）及び式（５）を用いて、各ウィンドウ幅の変位を計算する。
【００６６】
なお、式（４）は、式（３）に示す理想的な条件では、以下の式（６）となる。
【００６７】
【数６】
0005509058.tif 000007
【００６８】
ここで、式（５）に示すように、変位を計算するために必要となるのは、自己相関関数の位相差「angle｛ac(x)｝」であるが、式（６）で明らかなように、振幅値及び変位の双方が「i」によって変化するので、変位は、振幅値の重みの影響を受ける。つまり、変位は、振幅値の大きい位置の位相に引っ張られる。従来では、式（４）及び式（５）により算出された変位をウィンドウ幅の中央位置に与えていた。変位をウィンドウ幅の中央位置に与えることが、従来のエラストグラフィーで、均一なストレインを持つ生体組織でも、エコー強度の変化に対応するかのような模様が見えてしまう原因となっていた。
【００６９】
そこで、上述したように、補正部１６ｂの補正処理が行なわれる。まず、補正部１６ｂは、式（４）及び式（５）により算出された変位を与える位置を、以下の式（７）で算出される位置とする。すなわち、補正部１６ｂは、変位を与える位置を、振幅値「A1(x)」及び「A2(x)」の積の重心位置「xc(x)」とする。
【００７０】
【数７】
0005509058.tif 000008
【００７１】
図６は、第１の実施形態における補正処理を説明するための図である。すなわち、補正部１６ｂは、図６の（Ａ）に示すように、ウィンドウ幅の中心位置「x」に対して与えられていた変位「D(x)」を、式（７）で算出された重心位置「xc」に与えると補正する。
【００７２】
ここで、歪み計算部１６ｃは、各ウィンドウ幅で計算された重心位置を用いて変位の傾きを計算することで、ストレイン値を算出することが出来る。しかし、第１の実施形態では、ストレイン値の算出処理を容易にするために、変位の補間処理を行なう。
【００７３】
中心位置「x」は、等間隔にサンプリングされる位置となるが、「xc(x)」は、重心位置であるため、等間隔でなくなる。等間隔である方が、以降の処理が容易となるので、補正部１６ｂは、位置と変位を表わす座標「｛xc(x),D(x)｝」の点を補間により「｛x,Dc(x)｝」という等間隔の中心位置「x」で表わされた座標に変換する。
【００７４】
かかる処理は、数値解析ソフトウェアであるＭＡＴＬＡＢ（登録商標）を用いると、以下の式（８）で表される。式（８）に示す「interp1」は、座標列「｛xc,D｝」に対して、Ｘ座標の値が「xc」の代わりに「x」であった場合に、Ｙ座標の値「Dc」を「D」を補間することにより求める関数である。すなわち、補正部１６ｂは、他のウィンドウ幅にて計算された重心位置の変位を用いて、式（８）から、各ウィンドウ幅の中心位置の変位を補間処理により算出する。
【００７５】
【数８】
0005509058.tif 000009
【００７６】
式（８）に示すように、補正部１６ｂは、重心位置「xc」に与えた変位「D(x)」を、他のウィンドウ幅にて計算した重心位置「xc」に与えた変位「D(x)」を用いて、補間変位「Dc(x)」を計算する。そして、補正部１６ｂは、図６の（Ｂ）に示すように、計算した補間変位「Dc(x)」をウィンドウ幅の中心位置「x」に対して与える。
【００７７】
このように、補正部１６ｂは、式（７）及び式（８）を用いて補正処理を行なう。
【００７８】
そして、歪み計算部１６ｃは、補正部１６ｂにより計算された「Dc(x)」を用いて、以下の式（９）により歪み、すなわち、ストレイン値「S(x)」を計算する。
【００７９】
【数９】
0005509058.tif 000010
【００８０】
なお、歪み計算部１６ｃは、各中心位置にて、「Dc(x)」の傾きを最小２乗法によって求めることで、ストレイン値「S(x)」を計算しても良い。
【００８１】
以上の処理の後、画像生成部１７ａは、歪み計算部１６ｃが計算した歪み（ストレイン値）を用いてエラストグラフィーを生成する。具体的には、画像生成部１７ａは、Ｂモード処理部１４が生成したＢモードデータと、歪み計算部１６ｃが計算したストレイン値を用いてエラストグラフィーを生成する。そして、制御部１９の制御により、エラストグラフィーがモニタ２にて表示される。
【００８２】
図７は、補正処理の結果を説明するための図である。従来では、図３の（Ａ）に示すように、変位が深さ方向に沿って、振幅値が「０」となる「x」で段差を有する線となっていた。しかし、補正部１６ｂの補正処理が行なわれることで、図７の（Ａ）に示すように、変位は、深さ方向に沿って「x」に比例した直線として得られる。また、従来では、図３の（Ｃ）に示すストレイン画像のように、同一のストレイン値を有する生体組織であっても、均一の階調とはならなかった。しかし、補正部１６ｂの補正処理が行なわれることで、図７の（Ｂ）に示すように、同一のストレイン値を有する生体組織のストレイン画像は、同一階調で表現されるようになる。
【００８３】
以上、まとめると、変位計算部１６ａは、フレームバッファ１３から２つのＩＱ信号を読み出して、式（４）により自己相関関数を計算し、式（５）により変位計算を行なう。そして、補正部１６ｂは、式（７）により重心位置を計算し、式（８）により補間変位の計算を行なう。そして、歪み計算部１６ｃは、式（９）により、歪み（ストレイン値）を計算する。
【００８４】
なお、第１の実施形態は、補正部１６ｂが式（７）を用いて重心位置を計算し、歪み計算部１６ｃが重心位置と変位とを用いて歪みを計算する場合であっても良い。かかる場合でも、図７の（Ｂ）に示すように、同一のストレイン値を有する生体組織のストレイン画像は、同一階調で表現されるようになる。
【００８５】
次に、図８を用いて、第１の実施形態に係る超音波診断装置の歪み計算処理について説明する。図８は、第１の実施形態に係る超音波診断装置の歪み計算処理を説明するためのフローチャートである。なお、以下に示す歪み計算処理は、１フレーム分のスキャンが行なわれるごとに、同一フレーム内の受信信号を用いて繰り返し実行される。
【００８６】
図８に示すように、第１の実施形態に係る超音波診断装置は、フレームバッファ１３にＩＱ信号が格納されたか否かを判定する（ステップＳ１０１）。具体的には、第１の実施形態に係る超音波診断装置は、１フレーム分のスキャンにより収集されたＩＱ信号が格納されたか否かを判定する。ここで、ＩＱ信号が格納されていない場合（ステップＳ１０１否定）、第１の実施形態に係る超音波診断装置は、待機状態となる。
【００８７】
一方、ＩＱ信号が格納された場合（ステップＳ１０１肯定）、変位計算部１６ａは、自己相関関数を計算し（ステップＳ１０２）、変位を計算する（ステップＳ１０３）。すなわち、変位計算部１６ａは、各スキャンラインの各ウィンドウ幅の変位を計算する。例えば、変位計算部１６ａは、同一のスキャンラインにて２回超音波送信が行なわれることで生成された２つのＩＱ信号を用いて、式（４）により自己相関関数を計算し、式（５）により変位を計算する。
【００８８】
そして、補正部１６ｂは、自己相関関数を計算した時と略同一のウィンドウ幅にて各ＩＱ信号の振幅値の積の重心位置を計算し（ステップＳ１０４）、中心位置の変位を補間により計算する（ステップＳ１０５）。例えば、補正部１６ｂは、式（７）により各ウィンドウ幅の重心位置を計算する。そして、補正部１６ｂは、式（８）により各ウィンドウ幅の補間変位を計算する。
【００８９】
そして、歪み計算部１６ｃは、各ウィンドウ幅の中心位置と補間変異とから傾きを計算することで、各ウィンドウ幅の歪み（ストレイン値）を計算し（ステップＳ１０６）、処理を終了する。すなわち、歪み計算部１６ｃは、式（９）を用いて歪み（ストレイン値）を計算する。なお、歪み計算処理の後に、画像生成部１７ａは、エラストグラフィーを生成し、制御部１９の制御により、モニタ２は、エラストグラフィーを表示する。
【００９０】
上述してきたように、第１の実施形態では、変位計算部１６ａは、超音波プローブ１の直下に位置する生体組織に変形を加えながら、異なる時間に同一方向に送信した超音波の２つ以上の反射波に基づく受信信号を用いて、超音波送信方向にて所定間隔で設定された各ウィンドウ幅にある受信信号間の自己相関関数を計算し、当該計算した自己相関関数を用いて位相差を計算することで、前記各ウィンドウ幅における変位を計算する。そして、補正部１６ｂは、自己相関関数を計算した時と略同一のウィンドウ幅にて各受信信号の振幅値の積の重心位置を計算し、当該計算した重心位置を変位計算部１６ａにより計算された変位を与える位置として補正する。更に、補正部１６ｂは、各ウィンドウ幅の中心位置の変位を、他のウィンドウ幅にて計算された重心位置の変位を用いて補間することで計算する。
【００９１】
そして、歪み計算部１６ｃは、各ウィンドウ幅の補間された変位及び中心位置から傾きを計算することで、歪み（ストレイン値）を計算する。そして、画像生成部１７ａは、歪み計算部１６ｃが計算した歪みを用いてエラストグラフィーを生成する。
【００９２】
従って、第１の実施形態では、計算した変位を与える位置を補正し、更に、中心位置の補間変位を計算することで、エコーの振幅による誤差を除去することが出来る。その結果、第１の実施形態では、自己相関法を適用する場合に、エコー強度の影響を低減したエラストグラフィーを表示することが可能となる。また、第１の実施形態では、等間隔の中心位置における補間変位を計算するので、ストレイン値を計算する際の処理に要する負荷を軽減することが出来る。
【００９３】
（第２の実施形態）
第２の実施形態では、複合自己相関法を用いて変位を計算し、計算した変位に対するエコーの振幅による誤差を除去する方法について説明する。
【００９４】
複合自己相関法でも、超音波プローブ１の直下に位置する生体組織に変形を加えながら、異なる時間に同一方向に送信した超音波の２つ以上の反射波に基づく受信信号が用いられる。ただし、複合自己相関法では、一般的には、２フレーム間の同一スキャンライン上の同一位置の受信信号が用いられる。例えば、複合自己相関法では、解析信号が用いられる。なお、複合自己相関法では、自己相関法と同様に、１フレームでの同一スキャンライン上の複数の受信信号を用いても良い。
【００９５】
ここで、解析信号について、以下、数式を用いて説明する。時間「t」のＩＱ信号「IQ(t)」を、振幅値「A(t)」及び位相「φ(t)」を用いて式（１０）により定義すると、ＲＦ信号「rf(t)」は、以下の式（１１）により定義される。また、解析信号「RF(t)」は、以下の式（１２）により定義される。ここで、式（１１）の「f0」は、上述したように、基本周波数である。
【００９６】
【数１０】
0005509058.tif 000011
【００９７】
【数１１】
0005509058.tif 000012
【００９８】
【数１２】
0005509058.tif 000013
【００９９】
複合自己相関法では、２つの解析信号に対して相互相関法を用いて相互相関係数がピークとなる位置が得られる「ずれ量」を計算する。しかし、計算された「ずれ量」は、サンプル点単位である。そこで、複合自己相関法では、変位の計算精度を向上させるために、相互相関関数で位相差を調べ、ずれ量に位相差を足し合わせることで、１サンプル以下の変位まで計算する。ただし、複合自己相関法でも、第１の実施形態で説明したように、エコー強度の影響を受けることが知られている。
【０１００】
そこで、第２の実施形態に係る歪みデータ処理部１６は、エコー強度の影響を低減するために、図９に示すように構成される。すなわち、第２の実施形態に係る歪みデータ処理部１６は、ずれ量計算部１６ｄと、第１の補正部１６ｅと、変位計算部１６ｆと、第２の補正部１６ｇと、歪み計算部１６ｈとを有する。以下、図１０を参照して、各部の処理の概要を説明し、その後、数式を用いて、各部の処理を詳細に説明する。図１０は、第２の実施形態に係る歪みデータ処理部の概要を説明するための図である。
【０１０１】
ずれ量計算部１６ｄは、超音波プローブ１の直下に位置する生体組織に変形を加えながら、異なる時間に同一方向に送信した超音波の２つ以上の反射波に基づく受信信号（解析信号）を用いてずれ量を計算する。具体的には、ずれ量計算部１６ｄは、図１０に示すように、超音波送信方向にて所定間隔で設定された各ウィンドウ幅にある受信信号間（解析信号間）の相互相関係数がピークを取るずれ量を計算する。
【０１０２】
そして、第１の補正部１６ｅは、図１０に示すように、ずれ量計算部１６ｄがずれ量の取得時に用いた相互相関係数を計算した時と略同一のウィンドウ幅にて各受信信号（各解析信号）の振幅値の積の重心位置である第１の重心位置をずれ量を補正した補正ずれ量として計算する。
【０１０３】
そして、変位計算部１６ｆは、補正ずれ量によりずらした受信信号間（解析信号間）の相互相関関数の位相差から微小変位を計算する。そして、変位計算部１６ｆは、算出した微小変位に補正ずれ量を加算することで、図１０に示すように、各ウィンドウ幅における変位を計算する。
【０１０４】
そして、第２の補正部１６ｇは、図１０に示すように、第２の重心位置計算処理及び補間処理を行なうことで、変位補正を行なう。具体的には、第２の補正部１６ｇは、微小変位を計算した時と略同一のウィンドウ幅にて各受信信号（各解析信号）の振幅値の積の重心位置である第２の重心位置を計算する。そして、第２の補正部１６ｇは、計算した第２の重心位置を変位計算部１６ｆにより計算された変位を与える位置として補正する。そして、補正部１６ｇは、更に、各ウィンドウ幅の中心位置の変位を、他のウィンドウ幅にて計算された第２の重心位置の変位を用いて補間することで計算する。
【０１０５】
そして、歪み計算部１６ｈは、図１０に示すように、各ウィンドウ幅の補間された変位及び中心位置から傾きを計算することで、歪み（ストレイン値）を計算する。
【０１０６】
以下では、第２の実施形態に係る歪みデータ処理部１６が２つの解析信号を用いて上記の処理を行なう場合について、数式を用いて詳細に説明する。なお、第２の実施形態に係る歪みデータ処理部１６は、３つ以上の解析信号を用いて上記の処理を行なうことも可能である。また、第２の実施形態に係る歪みデータ処理部１６は、ＩＱ信号やＲＦ信号を用いる場合であっても適用可能である。
【０１０７】
ここで、以下では、２つの解析信号を、「RF1(x)」及び「RF2(x)」と定義する。ここでは、「RF」という記号を用いているがＲＦ信号ではなく解析信号である。通常、「RF1(x)」と「RF2(x)」との間には、１６ｍｓｅｃ〜５０ｍｓｅｃといった１フレームの時間差があるので、自己相関法のように同一地点での変位が半波長以内に収まる保証がない。
【０１０８】
相互相関法では、１サンプル点単位で相互相関係数を計算し、相互相関係数が最大となる変位を求める。ただし、相互相関法で計算される変位は、サンプル点間の距離精度しかないので、複合自己相関法では、相互相関法で計算された変位（上記の「ずれ量」に対応する）での相互相関関数の位相差を利用してサンプル点間の距離以下の精度を有する変位を求める。
【０１０９】
まず、ずれ量計算部１６ｄは、以下の式（１３）を用いて、解析信号間の相互相関係数「C(x,ξ)」を計算する。
【０１１０】
【数１３】
0005509058.tif 000014
【０１１１】
そして、ずれ量計算部１６ｄは、式（１３）から、位置「x」において、相互相関係数がピークとなる１サンプル単位の変位（単位：ｍ）を「ξ(x)」として計算する。
【０１１２】
ここで、従来の複合自己相関法では、ずれ量計算部１６ｄが計算した「ξ(x)」を、以下の式（１４）に代入することで相互相関関数「R(x)」を計算していた。
【０１１３】
【数１４】
0005509058.tif 000015
【０１１４】
そして、従来の複合自己相関法では、相互相関関数「R(x)」の位相差「angle｛R(x)｝」に「λ／２π」を乗算することで微小変位を計算していた。そして、従来の複合自己相関法では、微小変位にずれ量計算部１６ｄが計算した「ξ(x)」を加算することで変位「D(x)」を計算していた。すなわち、従来の複合自己相関法では、以下の式（１５）により変位「D(x)」を計算していた。
【０１１５】
【数１５】
0005509058.tif 000016
【０１１６】
しかし、式（１３）及び式（１４）には、式（１２）で定義される解析信号を参照しても明らかなように、振幅値の重みが掛かっている。
【０１１７】
そこで、第２の実施形態では、第１の補正部１６ｅは、振幅値の重みを除去するために「ξ(x)」を以下の式（１６）を代入することで、第１の重心位置である補正ずれ量「ξc(x)」を計算する。式（１６）により計算される補正ずれ量は、「RF1(x)」の振幅値「A1(x)」と、「RF2(x)」の振幅値「A2(x)」との重心となる。
【０１１８】
【数１６】
0005509058.tif 000017
【０１１９】
そして、変位計算部１６ｆは、第１の補正部１６ｅが式（１６）を用いて計算した「ξc(x)」を、式（１４）及び式（１５）に代入することで、変位を計算する。すなわち、変位計算部１６ｆは、以下に示す式（１７）及び式（１８）により変位を計算する。なお、式（１７）では、式（１４）の「R(x)」と区別するために、相互相関関数を「R2(x)」として示している。また、式（１８）では、式（１５）の「D(x)」と区別するために、変位を「D2(x)」として示している。
【０１２０】
【数１７】
0005509058.tif 000018
【０１２１】
【数１８】
0005509058.tif 000019
【０１２２】
ただし、補正ずれ量「ξc(x)」に相当するサンプル点がない時は、相互相関関数を計算することが出来ないので、第１の補正部１６ｅは、補正ずれ量「ξc(x)」を一番近いサンプル点の値に丸め込みを行う。例えば、第１の補正部１６ｅは、以下の式（１９）に示すように、四捨五入の値を計算する関数「round」を用いて、補正ずれ量「ξc(x)」を最寄りのサンプル点に補正された「ξ’c(x)」を計算する。例えば、「ξc(x)」が「４．３」であった場合、第１の補正部１６ｅは、式（１９）を用いて、「ξ’c(x)」を「４」として計算する。
【０１２３】
【数１９】
0005509058.tif 000020
【０１２４】
そして、変位計算部１６ｆは、「ξc(x)」の代わりに「ξ’c(x)」を式（１７）及び式（１８）を代入することで、変位を計算する。なお、「ξ’c(x)」により計算された変位に含まれる誤差は、以下の補正処理により除去される。
【０１２５】
ここで、式（１８）に示す位相差「angle｛R2(x)｝」には、第１の実施形態（自己相関法）で説明したように、振幅値による誤差が含まれている。そこで、第２の補正部１６ｇは、以下に示す式（２０）により、第２の重心位置「xc2(x)」を計算する。すなわち、第２の補正部１６ｇは、変位計算部１６ｆが計算した変位を与える位置として、第２の重心位置を計算する。これにより、変位計算部１６ｆが計算した変位「D2(x)」を与える位置が、ウィンドウ幅の中心位置「x  (x)」から第２の重心位置「xc2(x)」に補正される。
【０１２６】
【数２０】
0005509058.tif 000021
【０１２７】
ここで、歪み計算部１６ｈは、各ウィンドウ幅で計算された第２の重心位置を用いて変位の傾きを計算することで、ストレイン値を算出することが出来る。しかし、第２の実施形態では、第１の実施形態と同様に、ストレイン値の算出処理を容易にするために、変位の補間処理を行なう。
【０１２８】
すなわち、第２の補正部１６ｇは、ＭＡＴＬＡＢ（登録商標）の関数「interp1」により、位置と変位を表わす座標「｛xc2(x),D2(x)｝」の点を補間により「｛x, Dc(x)｝」という等間隔の中心位置「x」で表わされた座標に変換する。すなわち、第２の補正部１６ｇは、他のウィンドウ幅にて計算された第２の重心位置の変位を用いて、式（２１）から、各ウィンドウ幅の中心位置の変位を補間処理により算出する。
【０１２９】
【数２１】
0005509058.tif 000022
【０１３０】
式（２１）に示す「interp1」は、式（８）に示す「interp1」と同様に、座標列「｛xc2,D2｝」に対して、Ｘ座標の値が「xc2」の代わりに「x」であった場合に、Ｙ座標の値「Dc」を「D2」を補間することにより求める関数である。
【０１３１】
そして、歪み計算部１６ｈは、第２の補正部１６ｇにより計算された「Dc(x)」を用いて、以下の式（２２）により歪み、すなわち、ストレイン値「S(x)」を計算する。
【０１３２】
【数２２】
0005509058.tif 000023
【０１３３】
なお、歪み計算部１６ｃは、各中心位置にて、「Dc(x)」の傾きを最小２乗法によって求めることで、ストレイン値「S(x)」を計算しても良い。
【０１３４】
以上の処理の後、画像生成部１７ａは、歪み計算部１６ｃが計算した歪み（ストレイン値）を用いてエラストグラフィーを生成する。具体的には、画像生成部１７ａは、Ｂモード処理部１４が生成したＢモードデータと、歪み計算部１６ｃが計算したストレイン値を用いてエラストグラフィーを生成する。そして、制御部１９の制御により、エラストグラフィーがモニタ２にて表示される。
【０１３５】
以上、まとめると、ずれ量計算部１６ｄは、フレームバッファ１３から２つの解析信号を読み出して、式（１３）により「ずれ量」を計算する。第１の補正部１６ｅは、式（１６）を用いて、補正ずれ量（第１の重心位置）を計算する。また、第１の補正部１６ｅは、補正ずれ量がサンプル点にないときは、式（１９）を用いて補正ずれ量の近似値を計算する。そして、第１の補正部１６ｅは、補正ずれ量、又は、補正ずれ量の近似値と、式（１７）及び式（１８）を用いて、変位を計算する。そして、第２の補正部１６ｇは、式（２０）により第２の重心位置を計算し、式（２１）により補間変位の計算を行なう。そして、歪み計算部１６ｃは、式（２２）により、歪み（ストレイン値）を計算する。
【０１３６】
次に、図１１を用いて、第２の実施形態に係る超音波診断装置の歪み計算処理について説明する。図１１は、第２の実施形態に係る超音波診断装置の歪み計算処理を説明するためのフローチャートである。なお、以下に示す歪み計算処理は、１フレーム分のスキャンが行なわれるごとに、フレーム間の受信信号を用いて繰り返し実行される。
【０１３７】
図１１に示すように、第２の実施形態に係る超音波診断装置は、フレームバッファ１３に解析信号が格納されたか否かを判定する（ステップＳ２０１）。具体的には、第２の実施形態に係る超音波診断装置は、１フレーム分のスキャンにより収集された解析信号が格納されたか否かを判定する。ここで、解析信号が格納されていない場合（ステップＳ２０１否定）、第２の実施形態に係る超音波診断装置は、待機状態となる。
【０１３８】
一方、解析信号が格納された場合（ステップＳ２０１肯定）、ずれ量計算部１６ｄは、相互相関係数がピークをとるずれ量を計算する（ステップＳ２０２）。すなわち、ずれ量計算部１６ｄは、各スキャンラインの各ウィンドウ幅のずれ量を計算する。例えば、変位計算部１６ａは、隣接するフレーム間の同一のスキャンラインにて生成された２つの解析信号を用いて、式（１３）により自己相関係数を計算することで、各ウィンドウ幅のずれ量を計算する。
【０１３９】
そして、第１の補正部１６ｅは、ずれ量の取得時に用いた相互相関係数を計算した時と略同一のウィンドウ幅にて各解析信号の振幅値の積の重心位置である第１の重心位置（補正ずれ量）を計算する。例えば、第１の補正部１６ｅは、式（１６）により各ウィンドウ幅の第１の重心位置（補正ずれ量）を計算する（ステップＳ２０３）。なお、第１の補正部１６ｅは、補正ずれ量に相当するサンプル点がない時は、式（１９）により補正ずれ量の近似値を計算する。
【０１４０】
その後、変位計算部１６ｆは、各ウィンドウ幅の変位を計算する（ステップＳ２０４）。すなわち、変位計算部１６ｆは、補正ずれ量（又は、補正ずれ量の近似値）によりずらした受信信号間（解析信号間）の相互相関関数を計算する（式（１７））。そして、変位計算部１６ｆは、相互相関関数の位相差から微小変位を計算し、算出した微小変位に補正ずれ量（又は、補正ずれ量の近似値）を加算することで、各ウィンドウ幅における変位を計算する（式（１８））。これにより、変位計算部１６ｆは、各スキャンラインの各ウィンドウ幅の変位を計算する。
【０１４１】
続いて、第２の補正部１６ｇは、微小変位を計算した時と略同一のウィンドウ幅にて各解析信号の振幅値の積の重心位置である第２の重心位置を計算する（ステップＳ２０５）。例えば、第２の補正部１６ｇは、式（２０）により第２の重心位置を計算する。そして、第２の補正部１６ｇは、各ウィンドウ幅の中心位置の変位を補間により計算する（ステップＳ２０６）。すなわち、第２の補正部１６ｇは、各ウィンドウ幅の中心位置の変位を、他のウィンドウ幅にて計算された第２の重心位置の変位を用いて補間することで計算する。例えば、第２の補正部１６ｇは、式（２１）を用いて第２の重心位置を計算する。
【０１４２】
そして、歪み計算部１６ｈは、各ウィンドウ幅の中心位置と補間変異とから傾きを計算することで、各ウィンドウ幅の歪み（ストレイン値）を計算し（ステップＳ２０７）、処理を終了する。すなわち、歪み計算部１６ｈは、式（２２）を用いて歪み（ストレイン値）を計算する。なお、歪み計算処理の後に、画像生成部１７ａは、エラストグラフィーを生成し、制御部１９の制御により、モニタ２は、エラストグラフィーを表示する。
【０１４３】
上述してきたように、第２の実施形態では、ずれ量計算部１６ｄは、超音波プローブ１の直下に位置する生体組織に変形を加えながら、異なる時間に同一方向に送信した超音波の２つ以上の反射波に基づく受信信号を用いて、超音波送信方向にて所定間隔で設定された各ウィンドウ幅にある受信信号間の相互相関係数がピークを取るずれ量を計算する。そして、第１の補正部１６ｅは、ずれ量計算部１６ｄがずれ量の取得時に用いた相互相関係数を計算した時と略同一のウィンドウ幅にて各受信信号の振幅値の積の重心位置である第１の重心位置をずれ量を補正した補正ずれ量として計算する。変位計算部１６ｆは、補正ずれ量によりずらした受信信号間の相互相関関数の位相差から微小変位を計算し、当該算出した微小変位に補正ずれ量を加算することで、各ウィンドウ幅における変位を計算する。
【０１４４】
そして、第２の補正部１６ｇは、微小変位を計算した時と略同一のウィンドウ幅にて各受信信号の振幅値の積の重心位置である第２の重心位置を計算し、当該計算した第２の重心位置を変位計算部１６ｆにより計算された変位を与える位置として補正する。更に、第２の補正部１６ｇは、更に、各ウィンドウ幅の中心位置の変位を、他のウィンドウ幅にて計算された第２の重心位置の変位を用いて補間することで計算する。そして、歪み計算部１６ｈは、各ウィンドウ幅の補間された変位及び中心位置から傾きを計算することで、歪みを計算する。そして、画像生成部１７ａは、歪み計算部１６ｃが計算した歪みを用いてエラストグラフィーを生成する。
【０１４５】
従って、第２の実施形態では、補正ずれ量を計算することで、計算したずれ量のエコーの振幅による誤差を除去することが出来る。また、第２の実施形態では、計算した変位を与える位置を第２の重心位置に補正し、更に、中心位置の補間変位を計算することで、エコーの振幅による誤差を除去することが出来る。その結果、第２の実施形態では、複合自己相関法を適用する場合に、エコー強度の影響を低減したエラストグラフィーを表示することが可能となる。また、第２の実施形態では、等間隔の中心位置における補間変位を計算するので、ストレイン値を計算する際の処理に要する負荷を軽減することが出来る。
【０１４６】
なお、第２の実施形態は、第２の補正部１６ｇが第２の重心位置を計算し、歪み計算部１６ｈが第２の重心位置と変位とを用いて歪みを計算する場合でも適用可能である。かかる場合でも、エコーの振幅による誤差を除去することが出来、複合自己相関法を適用する場合に、エコー強度の影響を低減したエラストグラフィーを表示することが可能となる。
【０１４７】
なお、上記では、超音波診断装置内で自己相関法や複合自己相関法を用いたエラストグラフィーの生成処理が行なわれる場合について説明した。しかし、上記の第１の実施形態、又は、第２の実施形態は、超音波診断装置とは独立に設置された画像処理装置内で、エコーの振幅による誤差を除去したデータを生成する場合であっても適用可能である。具体的には、第１の実施形態に係る歪みデータ処理部１６や第２の実施形態に係る歪みデータ処理部１６の機能を有する画像処理装置が、受信信号を受信した上で、ストレイン値を計算する場合であっても良い。かかる場合、超音波診断装置は、画像処理装置により計算されたストレイン値を用いて、エラストグラフィーを生成及び表示する。
【０１４８】
また、上記の画像処理装置は、歪みデータ処理部１６とともに、Ｂモード処理部１４及び画像生成部１７ａを有する場合であっても良い。かかる場合は、画像処理装置内で、ストレイン値の計算とともに、エラストグラフィーの生成及び表示を実行する。
【０１４９】
また、図示した各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部または一部を、各種の負荷や使用状況などに応じて、任意の単位で機能的または物理的に分散・統合して構成することができる。さらに、各装置にて行なわれる各処理機能は、その全部または任意の一部が、ＣＰＵおよび当該ＣＰＵにて解析実行されるプログラムにて実現され、あるいは、ワイヤードロジックによるハードウェアとして実現され得る。
【０１５０】
以上、説明したとおり、第１の実施形態及び第２の実施形態によれば、エコー強度の影響を低減したエラストグラフィーを表示することが可能となる。
【０１５１】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【０１５２】
１  超音波プローブ
２  モニタ
３  入力装置
１０  装置本体
１１  送受信部
１２  位相検波回路
１３  フレームバッファ
１４  Ｂモード処理部
１５  ドプラ処理部
１６  歪みデータ処理部
１６ａ  変位計算部
１６ｂ  補正部
１６ｃ  歪み計算部
１７  画像処理部
１７ａ  画像生成部
１７ｂ  画像合成部
１８  画像メモリ
１９  制御部
２０  内部記憶部
【図１】
0005509058.tif 000024
【図２】
0005509058.tif 000025
【図３】
0005509058.tif 000026
【図４】
0005509058.tif 000027
【図５】
0005509058.tif 000028
【図６】
0005509058.tif 000029
【図７】
0005509058.tif 000030
【図８】
0005509058.tif 000031
【図９】
0005509058.tif 000032
【図１０】
0005509058.tif 000033
【図１１】
0005509058.tif 000034