JP4037689B2

JP4037689B2 - 超音波画像処理装置

Info

Publication number: JP4037689B2
Application number: JP2002154519A
Authority: JP
Inventors: 剛望月; 昌紀広瀬
Original assignee: Aloka Co Ltd
Current assignee: Hitachi Ltd
Priority date: 2002-05-28
Filing date: 2002-05-28
Publication date: 2008-01-23
Anticipated expiration: 2022-05-28
Also published as: JP2003339705A

Description

【０００１】
【発明の属する技術分野】
本発明は、超音波ビームにより走査された空間にて得られたエコーデータから超音波画像を生成する超音波画像処理装置に関し、関心部位の識別を容易とすることに関する。
【０００２】
【従来の技術】
超音波診断装置は、超音波ビームの走査により、走査面上や三次元空間に配置されたサンプリング点でのエコーデータを取得する。この二次元的、又は三次元的に配置されたエコーデータ（エコーデータ空間）から超音波画像が生成される。例えば、エコーデータ空間が二次元、すなわち面である場合には、その面に沿ったＢモード画像等の断層画像を生成することができる。また、エコーデータ空間が三次元である場合には、ボリュームレンダリング（Volume Rendering）法などを用いて三次元画像を生成することができる。
【０００３】
例えば、Ｂモード画像では、画素近傍の数点でのエコーデータを加重平均した補間値に応じて画素値が定められる。また、三次元画像の生成においては、その画像の各画素に対応して互いに平行な複数の視線が設定され、各視線に沿ったエコーデータを用いて所定の累積演算を行い、画素値が求められる。このように、従来の超音波画像では、エコー信号の振幅情報であるエコーデータに応じてその画素値が定められる。すなわち、基本的にエコーデータ、又は平均演算に基づく補間値や積算演算に基づく値が画素値とされる。
【０００４】
【発明が解決しようとする課題】
超音波画像は、生体内部の組織構造を可視化する目的で幅広く用いられている。その目的の１つに肝臓などの実質臓器内に生じ得る腫瘍の診断がある。この診断では、Ｂモード画像や三次元画像により腫瘍の広がり具合を観察しながら、その良悪性が判断される。
【０００５】
ここで、エコー信号は生体内の音響インピーダンスが相異する境界で発生するが、腫瘍と実質臓器とでは音響インピーダンスの相違が比較的小さく、その境界では大きなエコーデータが生じない場合がある。そのため、エコーデータの大きさに応じて画素値を定めて超音波画像を形成する従来の画像処理では、腫瘍と実質臓器とでの画素値の差が小さく、それらを画像上で明確に識別することが難しいという問題点があった。
【０００６】
本発明は上記問題点を解決するためになされたもので、異なる組織領域間で音響インピーダンスの差異が小さい場合においても、それら組織領域を容易に識別可能な超音波画像が得られる超音波画像処理装置を提供することを目的とする。
【０００７】
【課題を解決するための手段】
本発明に係る超音波画像処理装置は、エコーデータ空間内の各エコーデータを対象エコーデータとして、その対象エコーデータごとに、当該対象エコーデータを含む標本空間を設定する標本空間設定手段と、前記対象エコーデータごとに、前記標本空間内の複数のエコーデータを用いて、互いに異なる性質を持った複数の統計量からなる統計量セットを求める統計演算手段と、前記統計量セットに基づいて変換関数を選択し、当該変換関数により前記対象エコーデータを変換して変換エコーデータを生成するエコーデータ変換手段と、前記変換エコーデータを用いて超音波画像を生成する画像生成手段とを有し、前記変換関数は、前記統計量セットに基づいて選択された重み付け係数を乗じる演算であり、前記統計量セットを構成する複数の統計量は、それぞれ、所定の互いに異なる次数ｎを有するｎ次積率であり、当該ｎ次積率には、前記標本空間内の前記エコーデータの標準偏差値が含まれ、前記標準偏差値が小さい場合よりも前記標準偏差値が大きい場合における重み付け係数を小さくしてエコーデータのばらつきの幅を圧縮し、前記ｎ次積率には、標準偏差値に加え、前記標本空間内の前記エコーデータの平均値が含まれ、平均値と標準偏差値が共に比較的小さい血管に対して他の部位よりも小さな重み付け係数を選択し、血管に対応するエコーデータに当該小さな重み付け係数を乗じた後にエコーデータの大小関係を反転して前記変換エコーデータを生成することを特徴とする。
【０００８】
本発明によれば、対象エコーデータに対応した標本空間内の複数のエコーデータから得られる統計量セットに応じて変換関数が選択される。そして、各対象エコーデータに対して選択された変換関数で当該対象エコーデータを変換して変換エコーデータが生成される。例えば、ある２つの対象エコーデータが同じ値であっても、統計量セットを構成する統計量の値がそれら対象エコーデータで相違する場合には、互いに異なる変換関数が割り当てられ、互いに異なる変換エコーデータが得られる。変換関数は、識別したい組織の統計的性質に対応してあらかじめ定められる。このように、識別したい組織間でのエコーデータの相違が小さい場合であっても、比較的大きな差異を有する変換エコーデータを用いて画素値が定められ、超音波画像が生成される。
【０００９】
本発明の好適な態様は、前記変換関数が、前記統計量セットに基づいて選択された重み付け係数を乗じる演算である超音波画像処理装置である。
【００１０】
本発明の他の好適な態様は、前記統計量セットを構成する複数の統計量がそれぞれ、所定の互いに異なる次数ｎを有するｎ次積率である超音波画像処理装置である。
【００１１】
本発明の別の好適な態様は、前記ｎ次積率が、前記標本空間内の前記エコーデータの平均及び分散である超音波画像処理装置である。
【００１２】
本発明のさらに別の好適な態様は、前記標本空間が、中央部に前記対象エコーデータを有する超音波画像処理装置である。
【００１３】
他の本発明に係る超音波画像処理装置は、三次元空間に配列された各ボクセルでのエコーデータに対し、視線に沿った演算を行い、前記視線に対応する画素の画素値を求めるものであって、前記各ボクセルを対象ボクセルとして、その対象ボクセルごとに、当該対象ボクセルを含み前記視線方向に一列に並んだ複数のボクセルからなる標本空間を設定する標本空間設定手段と、前記対象ボクセルごとに、前記標本空間内の複数のボクセルのエコーデータを用いて、互いに異なる性質を持った複数の統計量からなる統計量セットを求める統計演算手段と、前記統計量セットに基づいて変換関数を選択し、当該変換関数により前記対象ボクセルのエコーデータを変換して変換エコーデータを生成するエコーデータ変換手段と、前記変換エコーデータを用いて前記画素値を求め超音波画像を生成する画像生成手段とを有し、前記変換関数は、前記統計量セットに基づいて選択された重み付け係数を乗じる演算であり、前記統計量セットを構成する複数の統計量は、それぞれ、所定の互いに異なる次数ｎを有するｎ次積率であり、当該ｎ次積率には、前記標本空間内の前記エコーデータの標準偏差値が含まれ、前記標準偏差値が小さい場合よりも前記標準偏差値が大きい場合における重み付け係数を小さくしてエコーデータのばらつきの幅を圧縮し、前記ｎ次積率には、標準偏差値に加え、前記標本空間内の前記エコーデータの平均値が含まれ、平均値と標準偏差値が共に比較的小さい血管に対して他の部位よりも小さな重み付け係数を選択し、血管に対応するエコーデータに当該小さな重み付け係数を乗じた後にエコーデータの大小関係を反転して前記変換エコーデータを生成することを特徴とする。
【００１４】
この場合にはエコーデータ空間は三次元空間であり、その中に設定される複数の視線に沿って画素値を定める演算が行われる。本発明によれば、その視線に沿って並んだエコーデータ列の一部分が標本空間とされる。
【００１５】
【発明の実施の形態】
次に、本発明の実施形態である超音波診断装置について図面を参照して説明する。本装置は、三次元空間にてエコーデータを取得し、視線に沿ってレンダリング処理を行い三次元画像を生成し表示することができる。
【００１６】
図１は、エコーデータが得られる三次元空間と、視線及び投影面と、複数の統計量が求められる標本空間であるカーネルとを示す模式図である。
【００１７】
エコーデータが得られる三次元空間２はボクセルに分割され、各ボクセルにエコーデータが割り付けられる。また、三次元空間２には、投影面４に形成される画像の各画素６に対応して平行な視線８が複数本、設定される。各ボクセルを対象ボクセルとして、その対象ボクセルごとにカーネル１０が設定される。カーネル１０は、例えば、高さ、幅、奥行き方向にそれぞれボクセルがＮｘ，Ｎｙ，Ｎｚ個並んだ直方体形状であり、奥行き方向ｚが視線方向に一致する。カーネル１０は、統計量を算出するために複数のボクセルを内包し、また一般にはカーネル１０は、そのうちの所定位置、例えば中央のボクセルが対象ボクセルとなるように設定される。以下の説明では、処理を簡単とするためにＮｘ＝Ｎｙ＝１、すなわち、カーネル１０は視線８に沿って一列に並んだ複数個のボクセルで構成されるものとする。
【００１８】
カーネル１０内の複数のエコーデータから、互いに異なる性質を持った複数の統計量からなる統計量セットが求められる。そして、あらかじめ用意された複数の変換関数のいずれかが、統計量セットに基づいて選択され、対象ボクセルのエコーデータはその選択された変換関数により変換され、変換エコーデータが生成される。
【００１９】
三次元画像のレンダリングでは、視線に沿ったレイキャスティング等の演算処理により画素６の画素値が決定される。この画素値を求める演算処理では通常、エコーデータが視線方向の位置（深さ）に応じて重み付けされ、視線方向に積算される。本装置では、この画素値を求める演算は変換エコーデータを用いて行われる。すなわち、画素値には、積算処理における深さに応じた重み付けとは別に、上述の統計量セットに応じた変換がさらに加味され反映される。
【００２０】
図２は、本発明の実施形態である超音波診断装置の概略のブロック構成図である。超音波探触子２０は、例えばリニアアレイ型の超音波振動子を有し、その超音波振動子を電子走査（リニア走査、セクタ走査）することによって、走査面が形成される。この走査面を当該面の法線方向へ例えば機械的に走査することによって、三次元空間でのエコーデータが取得される。ここで、超音波探触子２０の機械的な走査は図示しない駆動部によって行われている。
【００２１】
送受信部２２は超音波探触子２０に対して送信信号を供給すると共に、超音波探触子２０から出力された受信信号を処理するものである。その送受信部２２は受信信号をデジタル信号（エコーデータ）に変換した後に、データメモリ部２４へ出力する。
【００２２】
データメモリ部２４は、次に説明する前処理部２６での処理で用いられるエコーデータを一時的に保持するためのメモリであり、例えば三次元のエコーデータ空間を格納可能に構成される。
【００２３】
前処理部２６は、エコーデータから、上述した変換エコーデータを求める処理を行う。その詳細は後述する。
【００２４】
カーネルサイズ発生部２８は、ユーザがトラックボール等の入力機器３０を操作して指定した値Ｎを受けて、これをカーネルの透視線方向のサイズＮｚとして前処理部２６に設定する。すなわちＮｚ＝Ｎに設定される。
【００２５】
濃度変換部３２は、基本的には輝度反転処理を行う。例えば、血管はエコー強度がかなり小さく、そのまま表示すると、画像上暗くなり見にくい。また、単純に、エコーデータを増倍する変換ではノイズまで増幅される。そこで、そのようなエコーデータが小さい部分に関心がある場合に、濃度変換部３２は、所定の基準レベルＥsを用いて、例えば次式で表される処理を行い、エコーデータの大小関係を反転させる。ここで、Ｅoは前処理部２６から出力された変換エコーデータ、Ｅrは反転処理後のエコーデータを表す。
【００２６】
【数１】
Ｅr ＝Ｅs−Ｅo ………（１）
三次元データメモリ部３４は、濃度変換部３２から出力される反転処理後のエコーデータを保持するためのメモリであり、三次元のエコーデータ空間を格納可能に構成される。
【００２７】
三次元データレンダリング演算部３６は、三次元データメモリ部３４に格納されたエコーデータを用いて、各視線に対するレイキャスティング処理等のレンダリング処理を行って、当該視線に対応する画素値を算出する。
【００２８】
フレームバッファメモリ３８は、三次元データレンダリング演算部３６で算出された各画素値を格納する。その結果、フレームバッファメモリ３８には、１枚分の三次元画像の画像データが格納される。
【００２９】
フレームバッファメモリ３８から読み出された画像データは表示部４０に表示される。
【００３０】
以上、各部の動作は制御部４２により制御される。例えば、制御部４２は送受信部２２に対して送信同期信号、またデータメモリ部２４にはエコーデータの読み出し／書き込みを制御するＲ／Ｗ信号を供給する。また、制御部４２は、データメモリ部２４からのエコーデータの読み出し動作に同期して、前処理部２６の内部動作のタイミング制御を行う。さらに、制御部４２は、三次元データメモリ部３４に三次元のエコーデータ空間に対応するデータが格納されると、三次元データレンダリング演算部３６に対して演算開始信号を与える。
【００３１】
次に、本発明の主要な特徴を有する前処理部２６をより詳細に説明する。図３が、前処理部２６のブロック構成図である。制御部４２は、データメモリ部２４から視線に沿って順次、エコーデータを読み出す。このエコーデータが前処理部２６の統計演算部５０及び遅延時間部５２に入力される。
【００３２】
統計演算部５０は、処理対象とされるエコーデータ（対象エコーデータ）に対して、カーネルサイズ発生部２８により指定されるサイズのカーネルを設定する。ここではカーネルは視線に沿って一列に並ぶＮ個のエコーデータであり、対象エコーデータがその中央に位置するように設定される。統計演算部５０は、互いに異なる複数の統計量として、カーネル内に含まれるＮ個のエコーデータについて次数の異なるｎ次積率を計算し、それらを統計量セットとして出力する。ここでは統計量セットを構成するｎ次積率として、平均値Ｍ（次数ｎ＝１に相当）及び標準偏差値Ｓ（次数ｎ＝２に相当）が計算される。
【００３３】
乗算係数選択部５４は、統計演算部５０から与えられる統計量セットに基づいて、乗算係数を選択する。例えば、乗算係数選択部５４は、平均値Ｍと標準偏差値Ｓとをパラメータとし、乗算係数を割り当てる二次元のテーブルであり、例えばメモリを用いて構成される。
【００３４】
乗算部５６は、遅延時間部５２から順次出力されるエコーデータに、乗算係数選択部５４にて選択された乗算係数を乗じて、その乗算結果を出力する。この出力が前処理部２６の出力である変換エコーデータとなる。なお、本装置において、前処理部２６に入力されたエコーデータに対する変換関数は、入力データに係数を乗じた値を求めるものであり、乗算係数選択部５４がその係数を選択することにより、当該変換関数が具体的に選択されることとなる。
【００３５】
遅延時間部５２は、統計演算部５０及び乗算係数選択部５４での処理時間に応じた遅延を発生させる。すなわち、乗算部５６では、遅延時間部５２から入力されるエコーデータに、当該エコーデータを対象エコーデータとして統計演算部５０及び乗算係数選択部５４にて決定された乗算係数が乗じられる。
【００３６】
図４は、統計演算部５０の概略のブロック構成図である。統計演算部５０は、データメモリ部２４から入力されたエコーデータを順次、シフトレジスタ６０に入力する。入力されたエコーデータはシフトレジスタ６０内を順送りされる。シフトレジスタ６０に保持される所定個数のエコーデータのうち、カーネルサイズで指定されたＮ個の連続するデータ（ｅ₁，ｅ₂，…，ｅ_Nと表す）が、平均演算部６２及び標準偏差演算部６４に読み出されて、それらＮ個のエコーデータに基づいて平均値Ｍ及び標準偏差値Ｓが計算される。平均演算部６２は、シフトレジスタ６０から読み出されたエコーデータを加算回路７０で合計した後、この合計値を除算回路７２にて加算されたデータ数Ｎで除算して、次式で表される平均値Ｍを出力する。
【００３７】
【数２】

標準偏差演算部６４では、シフトレジスタ６０から読み出されたデータそれぞれが合成回路８０にて平均値Ｍを減算された後、その差分値が乗算回路８２にて二乗される。さらに、それら各二乗値を加算回路８４で合計し、この合計値を除算回路８６にてデータ数Ｎで除算した後、平方根演算回路８８にて平方根を求め、次式で表される標準偏差値Ｓを出力する。
【００３８】
【数３】

なお、標準偏差値Ｓとして、次式で表される簡易計算式による値を求めるように標準偏差演算部６４を構成してもよい。
【００３９】
【数４】

図５は、乗算係数選択部５４での処理の一例を説明する模式図である。この図は、乗算係数選択部５４を、平均値Ｍと標準偏差値Ｓとをパラメータとした二次元のルックアップテーブル（ＬＵＴ）９０として構成した場合の当該テーブルの概念図でもある。ここで、平均値Ｍ及び標準偏差値Ｓそれぞれは複数の範囲に区分され、それらＭ及びＳの範囲に基づいてＭ−Ｓ平面が格子状に区分される。図５に示す例では、Ｍ及びＳそれぞれが２つの範囲に区分され、Ｍ−Ｓ平面が４つの格子領域９２-1〜９２-4に区分される。これら各格子領域９２-1〜９２-4にそれぞれ乗算係数が設定される。ちなみに、この例では、Ｍは値“１１０”を境界として２つの範囲に区分され、Ｓは値“５０”を境界として２つの範囲に区分されている。
【００４０】
図５は、例えば、肝臓などの実質臓器内に生じた腫瘍及び血管を他の部位と識別容易に表示する場合に対応している。腫瘍は、実質臓器と同様、比較的大きな平均値Ｍを生じる点で類似するが、実質臓器より大きな標準偏差値Ｓとなるといった差異を有し得る。この場合には、実質臓器に設定されたカーネルから得られる統計量セットは、Ｍが大きくＳが小さい領域９２-1に属し、一方、腫瘍は当該領域９２-1よりＳが大きな領域９２-2に属すると考えられる。よって、領域９２-1に設定される乗算係数“１”に対し、腫瘍に対応する領域９２-2には相対的に小さい乗算係数“０．５”が設定される。
【００４１】
もう１つの関心部位である血管の内部は血液であり、これは一様な液体であるので、実質臓器に比べて平均値Ｍが小さくなり、また標準偏差値Ｓは実質臓器と程度の小さな値を取る。つまり、血管が占める領域に設定されたカーネルから得られる統計量セットは、Ｍ，Ｓ共に小さい領域９２-3に属すると考えられる。よって、領域９２-1に設定される乗算係数“１”に対し、血管に対応する領域９２-3には相対的に小さい乗算係数“０”が設定される。
【００４２】
また領域９２-4も実質臓器以外の部位が属すると考えられ、領域９２-1に設定される乗算係数“１”よりも相対的に小さい乗算係数、例えば“０．５”が設定される。
【００４３】
以上のように腫瘍及び血管が属する領域９２-2，９２-3に、実質臓器が属する領域９２-1より小さい乗算係数を割り当てることにより、前処理部２６から出力される変換エコーデータは、実質臓器より腫瘍、血管において比較的小さな値を取り得る。しかし、濃度変換部３２での反転処理により、三次元データメモリ部３４に格納される値は、実質臓器より腫瘍、血管において比較的大きな値となり、画像上で関心部位である腫瘍及び臓器が実質臓器より明るく表示される。
【００４４】
なお、標準偏差値Ｓが大きい領域９２-2に属する部位ではエコーデータの標準偏差、すなわちばらつきが大きいが、その領域９２-2の乗算係数を領域９２-1の乗算係数より小さくすることで、そのエコーデータのばらつきの幅が圧縮がされる。これにより領域９２-2に属する関心部位である腫瘍が、ざらつきの低減された見やすい画質で表示される。また、領域９２-3に関しては、大きな乗算係数をエコーデータに乗じて直接、明るい輝度値に対応した変換エコーデータを生成すると、ノイズまでも増幅されてしまい画質が劣化する。これに対して、小さい乗算係数を乗じた後、反転して明るい輝度値に対応するデータを生成する本装置の構成は、関心部位である血管が明るいだけでなく、ノイズの少ない良好な画質で表示される。
【００４５】
図６は、上述した前処理部２６の処理フロー図である。ユーザはトラックボール等の入力機器３０を操作してカーネルサイズを指定する（Ｓ１００）。統計演算部５０は、各視線に対する初期状態として、その視線に沿ったエコーデータ列の先頭にカーネルを位置させる（Ｓ１０５）。統計演算部５０は、設定したカーネル内のエコーデータの平均値Ｍ及び標準偏差値Ｓを計算する（Ｓ１１０）。得られたＭ及びＳに基づいて、乗算係数選択部５４は１つの乗算係数を選択する。遅延時間部５２は、乗算係数選択部５４からの乗算係数の出力タイミングに合わせて、その乗算係数に対応するカーネル内の対象エコーデータを乗算部５６へ出力する。乗算部５６では、この対象エコーデータと乗算係数とが乗じられる（Ｓ１１５）。その乗算結果が変換エコーデータとして前処理部２６から出力される（Ｓ１２０）。
【００４６】
統計演算部５０は視線に沿ってカーネルを順に移動させ（Ｓ１２５）、上記ステップＳ１１０〜Ｓ１２０が繰り返される。カーネルがある視線の終わりまで移動すると（Ｓ１３０）、視線を移動させ（Ｓ１３５）、以上のステップＳ１０５〜Ｓ１３０が全ての視線について繰り返される（Ｓ１４０）。
【００４７】
図７は、前処理部２６及び濃度変換部３２での処理例を示すエコーデータの模式図である。図において、縦軸がデータ値Ｅ、横軸が視線上での位置を表す。図７（ａ）は、データメモリ部２４から読み出されるエコーデータを表し、実質臓器に対応し、エコーデータのばらつきが比較的小さい（例えばＳ＝４０の）実質臓器区間２００、腫瘍に対応し、エコーデータのばらつきが比較的大きい（例えばＳ＝８０の）腫瘍区間２０２及び、血管に対応し、エコーデータが小さい血管区間２０４を含んでいる。ここで、実質臓器区間２００と腫瘍区間２０２とはエコーデータの平均値が同程度（例えばそれぞれＭ＝１２０，１００）である。
【００４８】
図７（ｂ）は、図５に示す乗算係数を乗じて得られた前処理部２６の出力でのエコーデータ（（１）式におけるＥo）を表す。前処理部２６での変換により、腫瘍区間２０２の平均値は実質臓器区間２００の半分程度（例えばＭ＝５０）となる。また腫瘍区間２０２でのエコーデータのばらつきが半分に低減される（例えばＳ＝４０となる）。血管区間２０４のエコーデータは０とされる。
【００４９】
前処理部２６の出力に対し、濃度変換部３２は、実質臓器区間２００でのエコーデータの平均値付近に設定する基準レベルＥs（図７（ｂ）の水平の点線２１０）に基づいて（１）式による反転演算を行い、（１）式におけるＥrを出力する。図７（ｃ）は、濃度変換部３２の出力値Ｅrを表す。反転することにより、実質臓器区間２００より腫瘍区間２０２及び血管区間２０４でのデータが大きくなり、これを用いて関心部位の腫瘍及び血管が明るく見やすく表示される。
【００５０】
【発明の効果】
本発明の超音波画像処理装置によれば、エコーデータそのものでは区別しにくい被検体内の異なる領域が、それら領域でのエコーデータの統計量セットの相違に基づいて、超音波画像上にて識別容易に表示される。
【図面の簡単な説明】
【図１】エコーデータが得られる三次元空間と、視線及び投影面と、複数の統計量が求められる標本空間であるカーネルとを示す模式図である。
【図２】本発明の実施形態である超音波診断装置の概略のブロック構成図である。
【図３】前処理部の概略のブロック構成図である。
【図４】統計演算部の概略のブロック構成図である。
【図５】乗算係数選択部での処理の一例を説明する模式図である。
【図６】前処理部での処理の概略のフロー図である。
【図７】前処理部及び濃度変換部での処理例を示すエコーデータの模式図である。
【符号の説明】
２０超音波探触子、２２送受信部、２４データメモリ部、２６前処理部、２８カーネルサイズ発生部、３０入力機器、３２濃度変換部、３４三次元データメモリ部、３６三次元データレンダリング演算部、３８フレームバッファメモリ、４０表示部、５０統計演算部、５２遅延時間部、５４乗算係数選択部、５６乗算部。

Claims

エコーデータ空間内の各エコーデータを対象エコーデータとして、その対象エコーデータごとに、当該対象エコーデータを含む標本空間を設定する標本空間設定手段と、
前記対象エコーデータごとに、前記標本空間内の複数のエコーデータを用いて、互いに異なる性質を持った複数の統計量からなる統計量セットを求める統計演算手段と、
前記統計量セットに基づいて変換関数を選択し、当該変換関数により前記対象エコーデータを変換して変換エコーデータを生成するエコーデータ変換手段と、
前記変換エコーデータを用いて超音波画像を生成する画像生成手段と、
を有し、
前記変換関数は、前記統計量セットに基づいて選択された重み付け係数を乗じる演算であり、
前記統計量セットを構成する複数の統計量は、それぞれ、所定の互いに異なる次数ｎを有するｎ次積率であり、当該ｎ次積率には、前記標本空間内の前記エコーデータの標準偏差値が含まれ、
前記標準偏差値が小さい場合よりも前記標準偏差値が大きい場合における重み付け係数を小さくしてエコーデータのばらつきの幅を圧縮し、
前記ｎ次積率には、標準偏差値に加え、前記標本空間内の前記エコーデータの平均値が含まれ、
平均値と標準偏差値が共に比較的小さい血管に対して他の部位よりも小さな重み付け係数を選択し、血管に対応するエコーデータに当該小さな重み付け係数を乗じた後にエコーデータの大小関係を反転して前記変換エコーデータを生成する、
ことを特徴とする超音波画像処理装置。
請求項１に記載の超音波画像処理装置において、
前記標本空間は、中央部に前記対象エコーデータを有することを特徴とする超音波画像処理装置。
三次元空間に配列された各ボクセルでのエコーデータに対し、視線に沿った演算を行い、前記視線に対応する画素の画素値を求める超音波画像処理装置において、
前記各ボクセルを対象ボクセルとして、その対象ボクセルごとに、当該対象ボクセルを含み前記視線方向に一列に並んだ複数のボクセルからなる標本空間を設定する標本空間設定手段と、
前記対象ボクセルごとに、前記標本空間内の複数のボクセルのエコーデータを用いて、互いに異なる性質を持った複数の統計量からなる統計量セットを求める統計演算手段と、
前記統計量セットに基づいて変換関数を選択し、当該変換関数により前記対象ボクセルのエコーデータを変換して変換エコーデータを生成するエコーデータ変換手段と、
前記変換エコーデータを用いて前記画素値を求め超音波画像を生成する画像生成手段と、
を有し、
前記変換関数は、前記統計量セットに基づいて選択された重み付け係数を乗じる演算であり、
前記統計量セットを構成する複数の統計量は、それぞれ、所定の互いに異なる次数ｎを有するｎ次積率であり、当該ｎ次積率には、前記標本空間内の前記エコーデータの標準偏差値が含まれ、
前記標準偏差値が小さい場合よりも前記標準偏差値が大きい場合における重み付け係数を小さくしてエコーデータのばらつきの幅を圧縮し、
前記ｎ次積率には、標準偏差値に加え、前記標本空間内の前記エコーデータの平均値が含まれ、
平均値と標準偏差値が共に比較的小さい血管に対して他の部位よりも小さな重み付け係数を選択し、血管に対応するエコーデータに当該小さな重み付け係数を乗じた後にエコーデータの大小関係を反転して前記変換エコーデータを生成する、
ことを特徴とする超音波画像処理装置。