(57)【要約】      （修正有）【課題】周波数領域でデジタルスキャン変換を実施して超音波映像を形成する方法を提供する。【解決手段】プローブから入力された受信信号で円筒座標系の空間映像データを形成し、前記空間映像データに周波数領域でスキャン変換されたフーリエ変換式を適用して直交座標系にスキャン変換された周波数領域の映像データを形成し、前記直交座標系にスキャン変換された周波数領域の映像データを逆フーリエ変換して直交座標系の空間映像データを形成し、前記直交座標系の空間映像データから超音波映像を形成する。【選択図】図３
2008212678.tif 000002
【特許請求の範囲】
【請求項１】
プローブから入力された受信信号から円筒座標系の空間映像データを形成し、
前記空間映像データに周波数領域でスキャン変換されたフーリエ変換式を適用して直交座標系にスキャン変換された周波数領域の映像データを形成し、
前記直交座標系にスキャン変換された周波数領域の映像データを逆フーリエ変換して直交座標系の空間映像データを形成し、
前記直交座標系の空間映像データから超音波映像を形成することを特徴とする超音波映像形成方法。
【請求項２】
前記空間映像データがｆ（ｒ、θ）で表現される時、前記空間映像データｆ（ｒ、θ）に周波数領域でスキャン変換されたフーリエ変換式を次の式のように適用して前記スキャン変換された周波数領域の映像データを形成し、
2008212678.tif 000013
前記式内のｕ、ｖ、ρ及びφは周波数領域の座標系で
2008212678.tif 000014
及び
2008212678.tif 000015
を満たすことを特徴とする請求項１に記載の超音波映像形成方法。
【請求項３】
前記プローブはコンベックスプローブ(ｃｏｎｖｅｘ  ｐｒｏｂｅ)、位相配列プローブ(ｐｈａｓｅｄ  ａｒｒａｙ  ｐｒｏｂｅ)またはセクタプローブ(ｓｅｃｔｏｒ  ｐｒｏｂｅ)であることを特徴とする請求項１または請求項２に記載の超音波映像形成方法。
【請求項４】
前記スキャン変換された周波数領域の映像データをゼロ挿入(ｚｅｒｏ  ｐａｄｄｉｎｇ)した後、前記逆フーリエ変換を実施することを特徴とする請求項３に記載の超音波映像形成方法。
【請求項５】
前記逆フーリエ変換は、高速逆フーリエ変換であることを特徴とする請求項４に記載の超音波映像形成方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波映像形成分野に関するものであって、特に周波数領域でデジタルスキャン変換(digital scan conversion)を実施して超音波映像を形成する方法に関する。
【背景技術】
【０００２】
超音波診断システムは超音波を用いて非破壊、非侵襲方式で対象体の内部構造を示す。超音波診断システムは超音波信号を送受信するためのプローブ、ビームフォーマ、アナログ-デジタル変換器、映像処理部を含む。スキャン方法または類型に応じて、プローブはコンベックスプローブ(convex probe)、位相配列プローブ(phased array probe)またはセクタプローブ(sector probe)になり得る。プローブは超音波信号と電気信号を相互変換するための多数の変換子(transducer)を含む。1回の超音波送信時にプローブの各変換子は別途に超音波信号を発生したり、複数の変換子が同時に超音波信号を発生させることもある。各変換子から送信された超音波信号は対象体内部の音響インピーダンス(acoustic impedance)の不連続面(反射体表面)で反射される。各変換子は個別に反射された超音波信号を電気的受信信号に変換する。アナログ-デジタル変換器は、受信信号をデジタル信号に変換し、ビームフォーマは超音波信号の集束点と各変換子の位置を考慮して送信集束及び受信集束する。映像処理部は集束された受信信号に基づいて超音波映像データを形成する。
【０００３】
コンベックスプローブ(convex probe)、位相配列プローブ(phased array probe)、セクタプローブ(sector probe)から得られる受信信号は円筒座標系で表現される信号である。映像処理部は受信信号を処理して空間領域(spatial domain)の円筒座標系映像データを形成する。直交座標系の画面に超音波映像をディスプレイするためには円筒座標系の空間映像データを直交座標系の空間映像データに変換しなければならない。このようなデータの座標変換をデジタルスキャン変換(digital scan conversion)という。
【０００４】
図1を参照すると、従来、デジタルスキャン変換器(digital scan converter)は方眼紙のような格子状の格納領域を有する半導体メモリで具現される。普通、デジタルスキャン変換器は水平に512、垂直に512の格納領域を有するメモリで具現される。各格納領域は画面の一つの画素に該当する。各格納領域には任意の大きさのビットで表現される映像データが格納される。例えば、8ビット映像の場合、一つの格納領域に格納される映像データの表現範囲は0から255に至る。
【０００５】
図2を参照すると、従来のスキャン変換ではメモリに格納されている円筒座標系の空間映像データf(r、θ)を直交座標系の空間映像データf(x、y)に変換する。この時、直交座標系の一つの画素Rの空間映像データは、画素R周辺の円筒座標系上の四つの画素A、B、C、Dのデータを双線形補間(bilinear interpolation)して得る。
【０００６】
このように、従来のデジタルスキャン変換では、円筒座標系の空間映像データを補間して直交座標系の空間映像データを得るので非常に多量のデータ計算が必要である。また、円筒座標系で表現される扇状映像の場合、プローブから遠くなるほど(即ち、映像の深さが深くなるほど)スキャンライン(SL)間の間隔が広くなり(d1<d2)、広い間隔に任意の値を作って埋めていかなければならないので、計算誤差が増加する問題がある。
【発明の開示】
【発明が解決しようとする課題】
【０００７】
本発明は、周波数領域でデジタルスキャン変換を実施して超音波映像を形成する方法を提供する。
【課題を解決するための手段】
【０００８】
本発明は、プローブから入力された受信信号から円筒座標系の空間映像データを形成し、前記空間映像データに周波数領域でスキャン変換されたフーリエ変換式を適用して直交座標系にスキャン変換された周波数領域の映像データを形成し、前記直交座標系にスキャン変換された周波数領域の映像データを逆フーリエ変換して直交座標系の空間映像データを形成し、前記直交座標系の空間映像データから超音波映像を形成する超音波映像形成方法を提供する。
【発明の効果】
【０００９】
本発明は、周波数領域でスキャン変換することによってスキャン変換時に映像の位置情報に対する依存性を減らすことができる。これによって、映像深さによって計算誤差が増加する問題を解決することができる。
【発明を実施するための最良の形態】
【００１０】
以下、添付した図面を参照して本発明の実施例を説明する。本発明の実施例で、空間映像データf(r、θ)はコンベックスプローブ(convex probe)、位相配列プローブ(phased array probe)またはセクタプローブ(sector probe)から得られた受信信号をデジタル変換して集束して得られたデータである。本発明は、円筒座標系の空間映像データf(r、θ)に周波数領域(frequency domain)でスキャン変換された(直交座標系に変換された)フーリエ変換式Fsc(u、v)を適用してスキャン変換された周波数領域の映像データを形成し、スキャン変換された周波数領域の映像データを逆フーリエ変換(Inverse Fourier Transformation, IFT)して直交座標系の空間映像データf(x、y)を得ることにその特徴がある。
【００１１】
図3は、本発明によるスキャン変換過程を概略的に示す。本発明は図3の(A)の通り扇状で表現される円筒座標系の空間領域(spatial domain)を長方形状の周波数領域に変換してスキャン変換を実施することによって、プローブから遠くなるほど(即ち、映像の深さが深くなるほど)スキャンライン間の間隔が広くなるのに伴って必要なデータ計算量を映像深さに関係なく一定にすることができ、計算誤差を減少させることができる。
【００１２】
以下、周波数領域でスキャン変換されたフーリエ変換式Fsc(u、v)を得るための方法を説明する。
直交座標系のフーリエ変換(FT)式は、次の式1のように表現される。
【００１３】
【数１】
2008212678.tif 000003
図4a及び図4bを参照すると、一点Pの円筒座標系と直交座標系の座標変換式は、次の式2の通りである。
【００１４】
【数２】
2008212678.tif 000004
式2から式1の指数部分は、数式3のように表現される。
【００１５】
【数３】
2008212678.tif 000005
式2及び3を用いて円筒座標系の空間映像データf(r、θ)をフーリエ変換(FT)すると、周波数領域(frequency domain)のフーリエ変換式F(ρ、φ)は、式4のように表現される。
【００１６】
【数４】
2008212678.tif 000006
F(ρ、φ)をスキャン変換し、Fsc(u、v)(スキャン変換されたフーリエ変換式)を得るために、式4で位置(u、v)に対応する(ρ、φ)でフーリエ係数を直接計算する。次の式5のような座標変換に従えば、周波数領域のスキャン変換された(直交座標系の)フーリエ変換式Fsc(u、v)は式6のように表現される。
【００１７】
【数５】
2008212678.tif 000007
【００１８】
【数６】
2008212678.tif 000008
前述した過程に従って円筒座標系の空間映像データf(r、θ)に周波数領域(frequency domain)のスキャン変換された(直交座標系に変換された)フーリエ変換式Fsc(u、v)を適用する。式6のように表現されるFsc(u、v)は、周波数領域の変換式であるので、空間領域(spatial domain)の映像データf(x、y)を得るためにフーリエ変換式Fsc(u、v)を逆変換する。本発明の実施例では式7のようにフーリエ逆変換の例として高速計算のためのIFFT(inverse fast Fourier transformation)を用いる。
【００１９】
【数７】
2008212678.tif 000009
一方、数式7でFsc(u、v)をゼロ挿入(zero padding)した後、IFFTを実施すると、周期が拡張されて任意の倍率で拡大された空間映像データf(x、y)を得ることができる。
【００２０】
空間映像データf(r、θ)が実数である場合、周波数領域でFsc(u、v)は、次の式8を満たすので、対称関係を用いて周波数領域映像データの1/2のみ計算すればいいので、計算量を減少させることができる。
【００２１】
【数８】
2008212678.tif 000010
また、ノイズ(noise)を除去するために周波数領域で次の式9のような低域通過フィルタを適用することができる。
【００２２】
【数９】
2008212678.tif 000011
式9で
2008212678.tif 000012
は、遮断周波数を示す。
【００２３】
以下、図5〜図9を参照して従来技術と本発明を具体的に対比して説明する。
図5は、各スキャンライン別にメモリに格納された円筒座標系の空間映像データｆ(ｒ、θ)で形成された映像(Ａ乃至Ｅ)を示す。図6に示したデジタルスキャン変換映像(Ａ１乃至Ｅ１)は従来技術によって映像データｆ(ｒ、θ)を空間領域(spatial domain)で双線形補間(bilinear transform)して得た映像である。図7は、本発明によって空間領域の映像データをフーリエ変換して周波数領域のデータを得、周波数領域でデジタルスキャン変換し、逆フーリエ変換して得た映像(Ａ２乃至Ｅ２)を示す。図8は、本発明によって前述した式7のフーリエ変換Fsc(u、v)で得たスペクトルである。図8に示したスペクトルは双線形補間でスキャン変換した後、フーリエ変換して得たスペクトル(図9)と同一または非常に類似することが分かる。図8と図9との比較から、本発明による周波数領域でのスキャン変換によって情報の損失なく超音波映像を得ることができることを予測することができる。
上述した実施例は、本発明の原理を応用した多様な実施例の一部を示したものに過ぎないことを理解しなければならない。本技術分野で通常の知識を有する者は、本発明の本質から逸脱することなく様々な変形が可能であることが明確に分かる。
【図面の簡単な説明】
【００２４】
【図１】スキャン変換器(digital scan converter)の格納領域を示す概略図である。
【図２】従来技術による円筒座標系と直交座標系との座標変換を説明するための概略図である。
【図３】本発明による周波数領域でのスキャン変換を用いた超音波映像形成過程を示す概略図である。
【図４ａ】図４ａは座標変換のための空間座標系及び周波数座標系を示す概略図である。
【図４ｂ】図４ｂは座標変換のための空間座標系及び周波数座標系を示す概略図である。
【図５】円筒座標系の空間映像データがスキャンライン別にメモリに格納される例を示す例示図である。
【図６】従来技術によって映像データを空間領域で双線形補間して得たデジタルスキャン変換映像の例示図である。
【図７】本発明によって空間領域の映像データをフーリエ変換して周波数領域のデータを得、周波数領域でデジタルスキャン変換し、逆フーリエ変換して得た映像の例示図である。
【図８】本発明によってスキャン変換されたフーリエ変換で得たスペクトルを示す例示図である。
【図９】双線形補間でスキャン変換した後、フーリエ変換して得たスペクトルを示す例示図である。
【符号の説明】
【００２５】
A〜E:円筒座標系の映像データf(r、θ)で形成された空間映像
A1〜E1:従来技術によってスキャン変換された映像
A2〜E2:本発明によってスキャン変換された映像
【図１】
2008212678.tif 000016
【図２】
2008212678.tif 000017
【図３】
2008212678.tif 000018
【図４ａ】
2008212678.tif 000019
【図４ｂ】
2008212678.tif 000020
【図５】
2008212678.tif 000021
【図６】
2008212678.tif 000022
【図７】
2008212678.tif 000023
【図８】
2008212678.tif 000024
【図９】
2008212678.tif 000025