(57)【特許請求の範囲】
【請求項１】
被検者に向けて超音波を送信するとともに、被検者から反射された超音波を受信する超音波送受信手段と、
前記超音波送受信手段が受信した超音波信号の位相成分からスペックル成分を抽出するスペックル成分抽出手段と、
を備えたことを特徴とする超音波画像処理装置。
【請求項２】
請求項１に記載の超音波画像処理装置であって、さらに、前記スペックル成分抽出手段が抽出した抽出結果を表示する表示手段を有することを特徴とする超音波画像処理装置。
【請求項３】
請求項１または２に記載の超音波画像処理装置であって、前記スペックル成分抽出手段は、位相成分の変化量を用いて前記超音波信号の位相成分からスペックル成分を抽出することを特徴とする超音波画像処理装置。
【請求項４】
請求項１または２に記載の超音波画像処理装置であって、前記スペックル成分抽出手段は、位相成分の少なくとも２つ以上の方向の変化量を用いることを特徴とする超音波画像処理装置。
【請求項５】
請求項１〜４のいずれかに記載の超音波画像処理装置であって、前記スペックル成分抽出手段は、スペックルを判別する判別関数を使用することを特徴とする超音波画像処理装置。
【請求項６】
請求項３または４に記載の超音波画像処理装置であって、前記スペックル成分抽出手段は、前記位相成分の変化量と比較してスペックルを判別する閾値を使用することを特徴とする超音波画像処理装置。
【請求項７】
請求項５に記載の超音波画像処理装置であって、前記判別関数は、スペックル及び非スペックルの位置がわかっている位相成分から生成された特徴量によって設定されることを特徴とする超音波画像処理装置。
【請求項８】
請求項７に記載の超音波画像処理装置であって、前記特徴量は、位相成分の変化量であることを特徴とする超音波画像処理装置。
【請求項９】
被検者に向けて超音波を送信するとともに、被検者から反射された超音波を受信する超音波送受信手段と、
前記超音波送受信手段が受信した超音波信号の振幅成分または包絡線成分と位相成分との組み合わせからスペックル成分を抽出するスペックル成分抽出手段と、
を備えたことを特徴とする超音波画像処理装置。
【請求項１０】
請求項１〜９のいずれかに記載の超音波画像処理装置であって、前記位相成分は、前記超音波信号の進行方向に垂直な方向のデータの分解能が、前記超音波送受信手段の前記超音波信号の送受信を行う素子の配列間隔以上であることを特徴とする超音波画像処理装置。
【請求項１１】
被検者に向けて超音波を送信するとともに、被検者から反射された超音波を受信し、
受信した超音波信号の位相成分からスペックル成分を抽出することを特徴とする超音波画像処理方法。
【請求項１２】
被検者に向けて超音波を送信するとともに、被検者から反射された超音波を受信する機能と、
受信して得られた超音波信号の位相成分からスペックル成分を抽出する機能と、
をコンピュータに実現させることを特徴とする超音波画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波画像処理装置及び方法並びにプログラムに係り、特に、超音波診断装置の画像データに対し、そのデータに含まれる位相情報を利用してスペックルノイズを抽出する技術に関する。
【背景技術】
【０００２】
従来より、超音波を用いて被検者の断層画像を取得して医療診断に供することが広く行われているが、超音波診断装置の画像データにおいて、スペックルノイズ（スペックル）と呼ばれる、被検体内の複数の反射体からの反射超音波のランダムな干渉によって生ずる低エコー情報が存在する。このスペックルノイズは、画像診断の妨げになるため、超音波診断画像からスペックルノイズを低減することが求められている。また、その一方で、スペックルノイズを含む信号（スペックルパターン）は、被検者の診断される組織の状態などによって変化し、組織性状を反映していると考えられるため、スペックルパターンを抽出して診断に活用する技術も考えられている。
【０００３】
例えば、超音波画像全体に含まれるスペックル情報を解析し、その解析結果を表示するものが知られている（特許文献１等参照）。これは具体的には、超音波受信信号の振幅情報を画像データに変換し、得られた原画像データに対し極大値、極小値からスペックルノイズの無い構造物画像データを生成し、原画像データからこのスペックルノイズの無い構造物画像データを差し引いた差分からスペックルノイズのみのスペックル画像データを算出するものである。
【特許文献１】特開２００６−２１２０５４号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかしながら、上記従来のスペックル抽出技術では、スペックルノイズの形状やエコーレベルにはバラツキがあり、非常に複雑であるため、振幅情報からスペックルノイズ部を正確に抽出するのは非常に困難であるという問題がある。
【０００５】
本発明は、このような事情に鑑みてなされたもので、スペックルの形状や周囲のエコーレベルに依存せずにスペックルを正確に抽出することのできる超音波画像処理装置及び方法並びにプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００６】
前記目的を達成するために、請求項１に記載の発明は、被検者に向けて超音波を送信するとともに、被検者から反射された超音波を受信する超音波送受信手段と、前記超音波送受信手段が受信した超音波信号の位相成分からスペックル成分を抽出するスペックル成分抽出手段と、を備えたことを特徴とする超音波画像処理装置を提供する。
【０００７】
これにより、振幅情報ではなく位相情報を用いてスペックルを抽出するようにしたため、スペックルの形状や周囲のエコーレベルに依存しないスペックルの抽出が可能となる。
【０００８】
また、請求項２に示すように、請求項１に記載の超音波画像処理装置は、さらに、前記スペックル成分抽出手段が抽出した抽出結果を表示する表示手段を有することを特徴とする。
【０００９】
スペックル抽出結果を表示することにより、診断への利用が促進される。
【００１０】
また、請求項３に示すように、請求項１または２に記載の超音波画像処理装置は、前記スペックル成分抽出手段は、位相成分の変化量を用いて前記超音波信号の位相成分からスペックル成分を抽出することを特徴とする。
【００１１】
位相成分の変化量を用いることによって、スペックルをより正確に抽出することができる。
【００１２】
また、請求項４に示すように、請求項１または２に記載の超音波画像処理装置は、前記スペックル成分抽出手段は、位相成分の少なくとも２つ以上の方向の変化量を用いることを特徴とする。
【００１３】
このように、複数方向の変化量を利用することで、抽出精度を高めることができる。
【００１４】
また、請求項５に示すように、請求項１〜４のいずれかに記載の超音波画像処理装置であって、前記スペックル成分抽出手段は、スペックルを判別する判別関数を使用することを特徴とする。
【００１５】
これにより、スペックル抽出の目的に応じた判別関数を使用することができる。
【００１６】
また、請求項６に示すように、請求項３または４に記載の超音波画像処理装置は、前記スペックル成分抽出手段は、前記位相成分の変化量と比較してスペックルを判別する閾値を使用することを特徴とする。
【００１７】
これにより、閾値と比較してスペックルの判別ができる場合には、スペックル抽出が容易となる。
【００１８】
また、請求項７に示すように、請求項５に記載の超音波画像処理装置であって、前記判別関数は、スペックル及び非スペックルの位置がわかっている位相成分から生成された特徴量によって設定されることを特徴とする。
【００１９】
これにより、スペックル抽出処理に応じた特徴量を用いて判別関数を作成することができる。
【００２０】
また、請求項８に示すように、請求項７に記載の超音波画像処理装置は、前記特徴量は、位相成分の変化量であることを特徴とする。
【００２１】
これによれば、位相成分の変化量をそのまま特徴量としてスペックル抽出に用いるため抽出処理が簡単となる。
【００２２】
また、同様に前記目的を達成するために、請求項９に記載の発明は、被検者に向けて超音波を送信するとともに、被検者から反射された超音波を受信する超音波送受信手段と、前記超音波送受信手段が受信した超音波信号の振幅成分または包絡線成分と位相成分との組み合わせからスペックル成分を抽出するスペックル成分抽出手段と、を備えたことを特徴とする超音波画像処理装置を提供する。
【００２３】
これによれば、位相情報にさらに振幅情報を組み合わせることにより、さらにスペックル抽出精度が良くなる。
【００２４】
また、請求項１０に示すように、請求項１〜９のいずれかに記載の超音波画像処理装置は、前記位相成分は、前記超音波信号の進行方向に垂直な方向のデータの分解能が、前記超音波送受信手段の前記超音波信号の送受信を行う素子の配列間隔以上であることを特徴とする。
【００２５】
これによれば、分解能が良くなり、さらにスペックル抽出の精度が向上する。
【００２６】
また、同様に前記目的を達成するために、請求項１１に記載の発明は、被検者に向けて超音波を送信するとともに、被検者から反射された超音波を受信し、受信した超音波信号の位相成分からスペックル成分を抽出することを特徴とする超音波画像処理方法を提供する。
【００２７】
これによれば、振幅情報ではなく位相情報を用いてスペックルを抽出するようにしたため、スペックルの形状や周囲のエコーレベルに依存しないスペックルの抽出が可能となる。
【００２８】
また、同様に前記目的を達成するために、請求項１２に記載の発明は、被検者に向けて超音波を送信するとともに、被検者から反射された超音波を受信する機能と、受信して得られた超音波信号の位相成分からスペックル成分を抽出する機能と、をコンピュータに実現させることを特徴とする超音波画像処理プログラムを提供する。
【００２９】
これにより、振幅情報ではなく位相情報を用いてスペックルを抽出するようにしたため、スペックルの形状や周囲のエコーレベルに依存しないスペックルの抽出が可能となる。
【発明の効果】
【００３０】
以上説明したように、本発明によれば、振幅情報ではなく位相情報を用いてスペックルを抽出するようにしたため、スペックルの形状や周囲のエコーレベルに依存しないスペックルの抽出が可能となる。
【発明を実施するための最良の形態】
【００３１】
以下、添付図面を参照して、本発明に係る超音波画像処理装置及び方法並びにプログラムについて詳細に説明する。
【００３２】
図１は、本発明に係る超音波画像処理装置を含む超音波診断装置の一実施形態の概略構成を示すブロック図である。この超音波診断装置は、超音波を用いて被検者の診断部位について超音波画像を撮影して表示するものであり、主に、超音波用探触子１０、信号処理部２０、画像処理部３０及び表示部４０を含んで構成されている。
【００３３】
超音波用探触子１０は、被検者の体内の診断部位に向けて超音波を送信するとともに、被検者の体内で反射してきた超音波を受診するものである。すなわち、超音波は構造物の境界のように音響インピーダンスが異なる領域の境界において反射されるため、超音波ビームを人体等の被検体内に送信して被検体内において生じた超音波エコーを受信し、超音波エコーが生じた反射点や反射強度を求めることにより、被検体内に存在する構造物の輪郭を検出することができる。
【００３４】
超音波用探触子１０は、例えば１次元の超音波トランスデューサアレイを構成する複数の超音波トランスデューサを備え、各超音波トランスデューサは、例えばＰＺＴ等の圧電素子の両端に電極を形成した振動子によって構成されるが、特に限定されるものではない。例えば、このように複数の超音波トランスデューサが１次元に配列されたリニアアレイプローブの他、被検体内を扇状に走査するセクタプローブ、複数の超音波トランスデューサが凸面上に配列されたコンベックスアレイプローブや、あるいは複数の超音波トランスデューサが２次元に配列された２次元アレイプローブを用いてもよい。またあるいは超音波内視鏡においてラジアル走査を行うメカニカルラジアルプローブでもよい。
【００３５】
超音波用探触子１０は、図示を省略した制御部の制御により超音波ビームを被検者内に送信し、リニア走査や、セクタ走査や、コンベックス走査や、ラジアル走査等の走査方式で被検者を走査する。超音波用探触子１０が発生した超音波は被検者体内に存在する反射体によって反射され、反射した超音波は超音波用探触子１０で受信される。超音波用探触子１０で受信された超音波受信信号は、電気信号に変換された後、図示を省略した送受信部を介して信号処理部２０に引き渡される。
【００３６】
信号処理部２０は、入力された受信信号に対し所定の信号処理を施すとともに、受け取った超音波受信信号から複数の受信データに遅延をかけて加算する位相整合加算、受信データをＩＱデータ（複素信号）に分離するＩＱ検波などの処理を行いＩＱデータを作成して画像処理部３０に出力する。
【００３７】
なお、ＩＱデータから次の式によって受信データの振幅Ａ及び位相θが算出される。
【００３８】
Ａ＝√（Ｉ２＋Ｑ２）、θ＝arctan（Ｑ／Ｉ）
ここで、記号√（Ｘ）は、Ｘの平方根を表す。また、受信データの波形の包絡線をとる場合には、その振幅は上記値を２倍して２Ａとしなければならない。
【００３９】
なお、信号処理部２０は、１ライン毎にデータを生成する手法でもよいが、例えば特開２００３−１８０６８８号公報に記載されたような、全素子の受信データをメモリに保存しておき、後処理でデータを生成する手法でもよい。この手法であれば、さらに方位分解能（素子方向の分解能）を向上させることができる。例えば、超音波信号の進行方向に垂直な方向に位相情報の分解能が素子間隔以上のデータを利用するように、超音波信号の進行方向に垂直な方向に高分解能な位相情報を利用することで、スペックルのランダムな位相変化をより正確に区別することができる。
【００４０】
また、一般的な超音波Ｂモード画像では、振幅から画像を生成し、位相情報は利用しないため、既存の超音波画像処理技術（Ｂモード）における超音波画像データは振幅データを指している。
【００４１】
次に信号の干渉と位相の関係について図を用いて説明する。干渉には、強め合う干渉と打ち消し合う干渉がある。強め合う干渉は波と波の位相差が小さい場合に生じ、弱め合う干渉は位相差がπに近い場合に生じる。
【００４２】
図２に、強め合う干渉の例を示す。図２（ａ）は干渉前を表し、基準波に対して位相差０．２（ｒａｄ）、位相差０．４（ｒａｄ）及び位相差０．６（ｒａｄ）という位相差の小さな３つの波を干渉させるようにしている。図２（ｂ）は干渉後を表し、破線で表された基準波に対し干渉波が実線で表されている。このように位相差が小さい場合には、強め合う干渉波が得られる。また、図２（ｂ）からわかるように、干渉波の山のピークは基準波の山のピークと近く、位相差が小さい干渉の波は、干渉後も基準波との位相差が小さくなる。
【００４３】
また、図３に、弱め合う干渉の例を示す。図３（ａ）は干渉前を表し、この場合基準波に対して位相差０．２（ｒａｄ）の他に位相差３．０（ｒａｄ）及び位相差３．２（ｒａｄ）という、位相差がπ（ｒａｄ）に近く、大きい波を干渉させている。このとき図３（ｂ）に示すように、破線で表された基準波に対して、実線で表されるように弱め合う干渉波が得られる。このように位相差が大きい波を干渉させると、弱め合う干渉波が得られる。また、図３（ｂ）からわかるように、干渉波の山のピークは基準波の山のピークと離れており、位相差が大きい干渉の波は、干渉後も基準波との位相差が大きい。
【００４４】
前にスペックルノイズは、被検体内の複数の反射体からの反射超音波のランダムな干渉によって生ずると言ったが、図３に示すような位相差が大きい干渉が、スペックルノイズを発生させる干渉となる。
【００４５】
また、ＩＱデータから得られる位相は干渉後の波と検波の波との位相差を表すが、ＩＱ検波の位相は誤差を含むため、干渉後の波から、位相差が大きい干渉の結果か、小さい干渉の結果かを確認するには、ＩＱデータの連続的な位相変化（例えば、隣接する画素間の差分）を見れば良い。これは、位相差が小さい波が連続すると位相変化は小さく、スペックルノイズ部になると位相が急変化し位相変化が大きくなるからである。
【００４６】
このように、位相が急に変化するという特徴を用いることでスペックルノイズを判別することができる。
【００４７】
画像処理部３０は、信号処理部２０で得られたＩＱデータに画像処理を施して、表示部４０に表示する。画像処理部３０における処理の中心は、信号処理部２０のＩＱデータから位相情報を抽出してスペックルの判別を行いスペックルを抽出することである。
【００４８】
図４に、画像処理部３０における、位相情報を利用したスペックル抽出処理の流れを示す。
【００４９】
図４に示すように、画像処理部３０は、位相変化量抽出部３２とスペックル判別部３４を有している。
【００５０】
信号処理部２０で得られたＩＱデータから位相データが位相変化量抽出部３２に入力される。位相は、前述したように、式θ＝arctan（Ｑ／Ｉ）で与えられる。位相変化量抽出部３２は、位相データから位相変化量を抽出する。抽出された位相変化量データはスペックル判別部３４に入力される。スペックル判別部３４は、位相変化量データを用いてスペックルを抽出し、スペックル抽出結果を出力する。
【００５１】
図５に、位相変化量抽出部３２における処理の一例を示す。
【００５２】
位相変化量抽出部３２は、位相データを得ると、まず方位方向（横方向）の位相変化量及び距離方向（深さ方向）の位相変化量、すなわちそれぞれの方向における画素の差分（方位方向差分及び距離方向差分）を算出する。
【００５３】
２次元のＢモード画像の場合、このように方位方向及び距離方向の位相変化量を求めることが望ましい。これは、スペックルノイズが方位方向または距離方向のどちらかに平行に存在する場合、その方向では位相変化が小さいが、それと直交する方向では位相変化が大きくなるためである。
【００５４】
次に、検波のずれによる一定量の位相変化を除去する。なお、スペックルの位相変化は他の部分に比べて急に変化するので、フィルタ（ローパスフィルタ、メディアンフィルタ）との差分データや位相データの２次微分など、変化量が大きくなる成分のみを抽出することによって、よりスペックルの判別をし易い特徴を持つデータを得ることができる。
【００５５】
そして、各方向で位相が急変化する成分が抽出され、各方向での位相変化量データが算出される。変化量は各方向の画素間の差分によって求められる。このとき隣接画素間における差分を用いることが最も好ましいが、適宜間引いた画素間における差分を用いてもよい。
【００５６】
このように、縦方向、横方向それぞれに対し、一定量の位相変化量を取り除いた結果の情報を基にしてスペックル／信号の判定を行うことができる。
【００５７】
なお、位相変化量は、この情報に限らず、複数画素の差分平均や斜め方向の差分を用いてもよい。さらに、別々の変化量ではなく、ベクトル成分のような一つのデータに変換してもよい。
【００５８】
スペックル判別部３４では、得られた位相変化量データに対し、スペックルノイズかどうかという２値的な判別や、どの程度スペックルノイズが含まれているかという多値的な判別を行い、その結果が目的に応じて出力される。もちろん、位相変化量そのものを、スペックルらしさを示すデータとして出力するようにしてもよい。
【００５９】
位相変化量抽出部３２によって抽出された方位方向位相変化量データ及び距離方向位相変化量データはスペックル判別部３４に入力される。スペックル判別部３４は、これらの位相情報データに基づいてスペックルの判別を行う。
【００６０】
図６に、スペックル判別部３４の一構成例を示す。
【００６１】
図６に示すように、スペックル判別部３４は、判別関数作成部３４１を備えている。判別関数作成部３４１は、位相変化量データに対し、スペックルノイズかどうかを判別するための判別関数を作成するものである。
【００６２】
判別関数作成部３４１は、振幅画像においてスペックルまたは非スペックルの位置がわかっている位相データを予め用意しておき、その位相変化量から作成した特徴量を元に判別関数を作成する。すなわち、スペックルであることがわかっている位相変化量データ３４４ａ及び非スペックルであることがわかっている位相変化量データ３４４ｂとから、それぞれ所定の特徴量が特徴量変換部３４２において算出され、これからスペックル判別関数が作成される。
【００６３】
ここで、特徴量は、位相変化量のデータの単一画素でも良い。ただし、単一画素の場合、差分をとって位相変化量を求めた場合には画素の僅かなずれが生じる、縦横方向から位相変化量を求めた場合には十字に交差する部分の中心の位相変化は大きくならないなどの問題があるため、注目画素の近傍画素の値や近傍画素との演算結果等、複数のデータを使用するのが望ましい。このとき、複数のデータは多次元となるため、閾値を設計し易いようにＰＣＡ（主成分分析）などを行って次元を下げるようにしても良い。
【００６４】
図９に、スペックル判別関数の一例を示す。
【００６５】
ここでは、縦方向位相変化量及び横方向位相変化量をそれぞれ特徴量（１）及び特徴量（２）とし、非スペックルノイズを〇で表し、スペックルノイズを×で表している。図９に示す例では、非スペックルノイズ〇とスペックルノイズ×の領域を分離する直線として判別関数が設定される。このように特徴量に変換した結果を基にしてスペックルノイズと非スペックルノイズを判別する関数（あるいは閾値）が設計される。なお、判別関数はこのような線形のものに限定されるものではない。
【００６６】
また、判別関数の設定方法は特に限定されるものではなく、例えば、ＳＶＭ（サポートベクターマシン）などの既知のデータ（学習データ）を利用した統計的手法（例えば、参考文献としてネロ・クリスティアニーニ、ジョン・テーラー著「サポートベクターマシン入門」共立出版などが挙げられる。）などの公知のクラス分類に使用される線形あるいは非線形の判別関数を用いることができる。もちろん、特徴量毎に閾値を与えるだけで判別可能であれば、閾値だけでスペックルを判別してもよい。また、位相変化量に変換することなく、連続的な位相データの画素といった位相変化がわかるデータを特徴量とした判別関数を設定してもよい。
【００６７】
図７に、ＳＶＭ（サポートベクターマシン）を使用した特徴量変換部３４２におけるスペックル抽出の判別関数生成処理の一例を示す。
【００６８】
図７に示すように、まずファントム画像から手作業でスペックル部分及び非スペックル部分をラベリングし、スペックル判別関数を作成するための既知データを作成する。なお、スペックル及び非スペックル部分をラベリングする際、曖昧な箇所についてはラベリングは行わなくてよい。次に、この既知データのスペックル部分及び非スペックル部分からスペックルの判別に用いる特徴量を抽出する。
【００６９】
ここでは特徴量として、図８に示すように、３×３画素の中央の画素ｃを注目画素とし、注目画素ｃとその上下左右の４つの近傍画素ａ、ｂ、ｄ、ｅに関し、それぞれ縦（距離）方向位相変化及び横（方位）方向位相変化の計１０個の特徴量を使用する。
【００７０】
次に、この縦方向及び横方向それぞれのラベリングされた画素（注目画素及びその近傍）の計１０箇所の位相変化量を特徴量としてＳＶＭ（サポートベクターマシン）を適用し、判別関数（スペックル判別器）を生成する。もちろん、ラベリングに使用するデータや特徴量は、判別結果が最適となるように変更しても良い．
【００７１】
このように、判別関数作成部３４１において、予めスペックル及び非スペックルであることがわかっている既知のデータから特徴量を抽出してスペックル判別関数を作成しておく。そして、実際の超音波診断において、超音波用探触子１０の走査によって得られたデータから信号処理部２０によって生成されたＩＱデータから得られた位相情報に基づいてスペックル抽出が行われる。
【００７２】
すなわち、図６において、方位（横）方向位相変化量データ３４６ａ及び距離（縦）方向位相変化量データ３４６ｂが入力されると、特徴量変換部３４３ではこれをスペックル判別に用いる特徴量に変換し、判別関数作成部３４１で予め作成されたスペックル判別関数３４７を用いてスペックルであるかどうか判別し、スペックル抽出が行われる。そして、スペックル抽出結果３４８が出力され表示部４０に表示される。この表示は、スペックルのみを表示してもよいし、原画像である振幅画像に重ねて表示するようにしてもよい。
【００７３】
判別結果は、２値的にスペックルであるかどうかを示すだけでなく、閾値との差を多値的にどの程度スペックルノイズが含まれているかを示すスペックルらしさとして出力してもよい。また、多値的に出力する場合には、ＬＵＴ（ルックアップテーブル）などでさらに値を調整してもよい。
【００７４】
なお、スペックル判別関数は、超音波の送受信の条件などによって変化するため、実際の装置の場合には条件毎に判別関数を設定するのが望ましい。
【００７５】
以上説明したように、本実施形態においては、振幅情報ではなく位相情報を用いてスペックルを抽出するようにしているが、振幅ではスペックルの形状、エコーレベルによっては、スペックルの判別が困難な場合があるため正確にスペックルを抽出できない場合があるが、位相情報を用いることにより、スペックルの形状や周囲のエコーレベルに依存しないスペックルの抽出が可能となった。
【００７６】
このように、位相情報のみでスペックルを抽出することができるが、振幅情報を特徴量に加えるなど、振幅情報（振幅成分あるいはその包絡線成分）と位相情報を組み合わせてスペックルを抽出するようにしてもよい。この場合、例えば図６において、位相情報（方位方向位相変化量データ３４６ａ、距離方向位相変化量データ３４６ｂ）が特徴量変換部３４３に入力されているが、これらのデータとともに振幅情報をも特徴量変換部３４３に入力して、振幅情報（画素値）も加えて特徴量の次元（図９に示す例では２次元）を増加して、スペックル抽出を行うようにしてもよい。
【００７７】
このように、振幅情報と位相情報の両方を用いてスペックル抽出を行うことで、より正確なスペックル抽出が可能となる。
【００７８】
また、本実施形態の超音波画像処理装置は、図示を省略した制御部に付属したメモリに格納された超音波画像処理プログラムによって制御される。すなわち、制御部によってメモリから超音波画像処理プログラムが読み出され、該超音波画像処理プログラムに従って、被検者に向けて超音波を送信するとともに、被検者から反射された超音波を受信する機能と、受信して得られた超音波信号の位相成分からスペックル成分を抽出する機能が実行される。
【００７９】
なお、超音波画像処理プログラムは、このように制御部に付属のメモリに格納されるものに限定されず、該超音波画像処理プログラムを例えばＰＣカードやＣＤ−ＲＯＭなど、本超音波画像処理装置に着脱可能に構成されるメモリ媒体（リムーバブル媒体）に記録しておき、リムーバブル媒体に対応するインターフェイスを介して本装置に読み込むように構成してもよい。
【００８０】
以上、本発明の超音波画像処理装置及び方法並びにプログラムについて詳細に説明したが、本発明は、以上の例には限定されず、本発明の要旨を逸脱しない範囲において、各種の改良や変形を行ってもよいのはもちろんである。
【図面の簡単な説明】
【００８１】
【図１】本発明に係る超音波画像処理装置を含む超音波診断装置の一実施形態の概略構成を示すブロック図である。
【図２】強め合う干渉の例を示すグラフであり、（ａ）は干渉前を表し、（ｂ）は干渉後を表す。
【図３】弱め合う干渉の例を示すグラフであり、（ａ）は干渉前を表し、（ｂ）は干渉後を表す。
【図４】画像処理部における位相情報を利用したスペックル抽出処理の流れを示すブロック図である。
【図５】位相変化量抽出部における処理の一例を示すブロック図である。
【図６】スペックル判別部の一構成例を示すブロック図である。
【図７】特徴量変換部におけるＳＶＭ（サポートベクターマシン）を使用したスペックル抽出の判別関数生成処理の一例を示すブロック図である。
【図８】特徴量として用いる画素の例を示す説明図である。
【図９】スペックル判別関数の例を示す説明図である。
【符号の説明】
【００８２】
１０…超音波用探触子、２０…信号処理部、３０…画像処理部、４０…表示部、３２…位相変化量抽出部、３４…スペックル判別部、３４１…判別関数作成部、３４２、３４３…特徴量変換部
【図１】
0005089442.tif 000002
【図２】
0005089442.tif 000003
【図３】
0005089442.tif 000004
【図４】
0005089442.tif 000005
【図５】
0005089442.tif 000006
【図６】
0005089442.tif 000007
【図７】
0005089442.tif 000008
【図８】
0005089442.tif 000009
【図９】
0005089442.tif 000010