JP2006231069A - イメージングシステムの表示プロセッサ - Google Patents

イメージングシステムの表示プロセッサ Download PDF

Info

Publication number
JP2006231069A
JP2006231069A JP2006103636A JP2006103636A JP2006231069A JP 2006231069 A JP2006231069 A JP 2006231069A JP 2006103636 A JP2006103636 A JP 2006103636A JP 2006103636 A JP2006103636 A JP 2006103636A JP 2006231069 A JP2006231069 A JP 2006231069A
Authority
JP
Japan
Prior art keywords
signal
resolution
image signal
image
contrast
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2006103636A
Other languages
English (en)
Inventor
Kutay Ustuner
アースターナー,カーテイ
Matthew I Haller
ハーラー,マセッウ,イサック
Ting-Lan Ji
ジ,ティン−ラン
Pai-Chi Li
リ,パイ−チ
Can Cinbis
チンビス,キャン
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Siemens Medical Solutions USA Inc
Original Assignee
Acuson Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Acuson Corp filed Critical Acuson Corp
Publication of JP2006231069A publication Critical patent/JP2006231069A/ja
Pending legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0202Two or more dimensional filters; Filters for complex signals
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52023Details of receivers
    • G01S7/52025Details of receivers for pulse systems
    • G01S7/52026Extracting wanted echo signals
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52046Techniques for image enhancement involving transmitter or receiver
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8995Combining images from different aspect angles, e.g. spatial compounding
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52053Display arrangements
    • G01S7/52057Cathode ray tube displays
    • G01S7/5206Two-dimensional coordinated display of distance and direction; B-scan display
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52053Display arrangements
    • G01S7/52057Cathode ray tube displays
    • G01S7/52071Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric imaging

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Physics (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Image Processing (AREA)

Abstract

【課題】 超音波イメージングシステムの表示プロセッサの組織種別の差別化能力を改善する。
【解決手段】 このプロセッサには、高空間分解能画像信号IDから平滑化された画像信号IPを発生させる2次元フィルタ(12)が含まれる。信号IPは高いコントラスト分解能に最適化されており、組織に対する良好な差別化を呈し、一方、信号IDは高い細部分解能および精細な構造的細部の表示に最適化されている。信号IPおよびIDはアドレス入力としてルックアップテーブル(14)に与えられ、このテーブルが出力画像I0を発生する。この出力信号は信号IDの細部分解能と信号IPのコントラスト分解能とを結合する。信号I0は細部分解能を輝度として、またコントラスト分解能をカラーとして表示できる。信号I0はまた、信号IPおよびIDの重み付け結合または加算の量として形成することができる。
【選択図】 図1

Description

本発明はイメージングシステムに関し、とくに、このようなシステムにおける多重種類の分解能の同時表示を強調できる表示プロセッサに関する。
超音波画像を診断的に使用することで2つの異なる目的を果たすことができ、その一つは組織を構造的に表示することであり、もう一つは種々の種類の組織を差別化することである。この画像が高い空間分解能(すなわち、細部分解能(detail resolution))を有していると、ユーザの能力を最大限に発揮させて精細な構造的細部を検出することができる。また、この画像のコントラスト分解能(contrast resolution)が高いと、やはりユーザの能力を最大限に発揮させて組織の種類差を検出することができる。
コントラスト分解能に関する基本的な制限の一つはスペックルノイズである。スペックルノイズはコヒーレントなイメージングシステムに固有なものであり、分解できないほどに細かな構造物により散乱された波のコヒーレントな(位相感知の)干渉によって生じる。スペックルノイズは、対象物の構造の詳細に関する情報を有していないが、その低周波成分は対象物の後方散乱の平均強度を示す情報を有している。このため、コントラスト分解能を向上させるには、スペックルノイズの高周波成分を減少させる必要がある。
コントラスト分解能を改善するため様々な信号処理技法が使用されている。ローパス・ビデオフィルタ、ローパス・スキャン変換フィルタ、さらには表示モニタの帯域幅の限定化さえも、スペックルノイズの高周波数の空間変化を減じるために使用することができる。このようにして、ある組織から別の組織への平均後方散乱強度の差は著しく目立たせることができる。しかし、コントラスト分解能を改善するためのこのようなアプローチに従う空間的スムージングによって細部分解能を著しく下げてしまう。
Bamperの米国特許第4783839号には、スペックルノイズを減じるためのシステムが開示されており、このシステムは画像にかけるスムージングの量を、画像に似たスペックルの個々の領域がどの程度接近しているかの関数として変化させる。多くの場合、スペックルノイズは対象物の分解可能な構造的細部に重畳しており、診断的に有益な細部分解能はしたがって失われてしまう。
別の信号処理技法として非線形表示マッピング法があり、この技法もコントラスト分解能を改善するために使用されている。とくに、スムージング用ビデオフィルタにカスケード接続する場合、非線形表示マップを用いると、選択した表示強度レベルではコントラスト分解能を上げることができるが、その代わりに、そのほかの表示強度レベルでの細部分解能及びコントラスト分解能が低下するという犠牲を伴う。
スペックルノイズは画像形成の過程においても、種々の非干渉性アベレージング(合成)法のどれかの方法で減らすことができ、その技法には空間的、周波数的、また時間的合成法が含まれる。
画像形成における空間的合成、横、軸、または深さ方向でトランスデューサをシフトさせることで実現できる。周波数的合成は、入力パルスの通過帯域を周波数ドメイン上の複数のサブ帯域に分割することで行える。時間的合成により画像を生成するには、距離(range)強度サンプルを非干渉的に加算することで行える。
これらの合成法は、スペックル画像の全体平均は元の対象物の非干渉性平均と同じであるという事実を活用している。斑点状の姿で現れるスペックルノイズはイメージング装置についてのみの情報しかもっておらず、イメージングする対象物についての情報は持っていないから、スペックル変動値は、元々の目標コントラストを変化させることなく、独立した測定値を非干渉的に平均化することで減じることができる。スペックル変動を減らした結果として、細部分解能を減少させるという代償を払うが、コントラスト分解能は改善できる。
Lipschultzの米国特許第5224483号には、画像中の血液プール(pool)領域のクラッタを減らすようにして超音波画像を高めるシステムが開示されている。この血液プール領域はローパスフィルタリング及び非線形強度マッピングによって識別され、実質的に組織領域のある値と実質的に血液プール領域の別の値とを持つマスク信号が発生される。画像信号は次いでそのマスク信号で変調されて、血液プールのクラッタが殆ど除去される。このマスク信号は好ましくは、厳密な2値ではなく、マスクする領域とマスクしない領域との間で信号値がスムーズに移動する中間レベルの値を含み、これにより最終画像の不自然な表示が排除される。この技術は血液プール領域のクラッタを抑制する場合のみ適用されるもので、組織領域のコントラスト分解能の改善の問題には言及していない。
したがって、表示プロセッサに対するニーズは、画像上の診断的に重要で精細な構造的細部を失うことなく、コントラスト分解能を高めて組織種別の差別化能(肝臓か腎臓か、健康な組織か病変部か)を改善することである。
本発明によれば、イメージングシステムの表示プロセッサは共通エンティティを描写する第1および第2の画像信号に応答する。この第1の画像信号の細部分解能は第2の画像信号のそれよりも大きく、一方、第2の画像信号のコントラスト分解能は第1の画像信号のそれよりも大きい。表示プロセッサは、第1の画像信号による細部情報および第2の画像信号によるコントラスト情報を合成する出力信号を、この2つの信号に関する関数として形成する。
この表示プロセッサは好適には、2次元マップといった表示発生器を備える。ある実施例では、出力信号は、第1の画像信号による細部分解能を表示する強度と第2の画像信号によるコントラスト分解能を表示するカラーとにより特徴付けられる。ほかの実施例では、出力信号は第1および第2の画像信号の重み付け加算として形成され、その重み付け加算は第1および第2の画像信号の両者に応じて変わる重み付けファクタを有する。
図面を参照すると、図1は本発明の1つの好適な実施例のブロック図であり、任意の適宜な超音波イメージングシステムSが発生する高細部分解能画像信号IDに応ずるものである。このシステムSは従来使用されているものであり、本実施例では、上記信号IDは標準のBモード画像信号で、その強度は組織の後方散乱係数の関数として変わる。勿論、本発明はBモード画像信号の使用に限定されるものはなく、組織のほかのどのような適宜な音響的特性(減衰係数、音速、弾性など。しかし、これらに限定されない。)をモニタする超音波イメージングシステムも同様に本発明に使用できる。この画像信号IDは処理されていない画像信号、すなわち、何がしかの量のビデオフィルタリング(スムージングまたはエッジ強調のいずれか)を要する画像信号であることができる。
図1に示すように、この画像信号IDは入力信号として表示プロセッサ10に加えられ、この表示プロセッサは2次元フィルタ12と、2次元ルックアップテーブルといった表示発生器14を有する。この2次元フィルタ12は、処理された高コントラスト分解能画像信号IPを出力として発生するビデオフィルタで形成できる。この画像信号IPの細部分解能は低いが、上記画像信号IPよりもコントラスト分解能は高い。
フィルタ12で行うスムージングの程度は、コントラスト分解能、すなわち種々の組織の種別を差別化する能力を強調するように選択される。一般に、前記信号IPは組織の種別の変化により変わるマルチレベルの信号で、この信号により固体組織の異なる種類の主要な領域は互いに差別化できる。例えば、血液プールおよび少なくとも第1、第2の種類の固体組織を含む領域の画像信号IPにより、好適に、血液プールを固体組織から差別化することができ、さらに、第1、第2の種類の固体組織を互いに差別化することができる。2次元空間フィルタ12のサイズおよび形状はアプリケーションや患者の体の形状に合うように広範囲に変えることができる。このフィルタとしてマルチ・ユーザの選択型フィルタを提供できるし、シングルフィルタをユーザ・アジャスタブルで使用することもできる。
この2次元空間フィルタ12については多くのインプレメンテーションが当業者には知られている。図2乃至4に示す1つの適した手法は、画像メモリ24、方位フィルタ26、レンジフィルタ28、およびスケーリング回路30を有する。この画像メモリ24は前記細部画像信号IDから得た音響データから成るラインを記憶しており、レンジフィルタ28はそのラインに沿ってフィルタリングし、さらに方位フィルタ26はライン間でフィルタリングする。
図3に示す如く、方位フィルタ26は2個のFIFOレジスタ32、34を含み、この両方のレジスタがメモリ24から同一の画像信号を受ける。このレジスタはクロック信号の制御の元に加算器36に出力信号を供給する。遅延要素38がレジスタ34の出力信号を所望のライン数に等しい量Dだけ時間的にシフトさせる。lnはライン数で、tnはレンジサンプル数とするとき、オリジナルの画像信号はID(ln,tn)であるとすると、加算器36の出力ISは以下のように表される。
Figure 2006231069
この信号ISは次いでアキュムレータ40に供給される。このアキュムレータは加算器42およびラインFIFOレジスタ44を有する。アキュムレータ40はラインデータを積分し、以下のような出力信号IINTを発生する。
Figure 2006231069
このため、信号IINTはD隣接ラインの合計となる。これはサイズDのボックスカーフィルタに相当する。
図4にレンジフィルタ28の詳細を示す。このフィルタはライン数の代わりにレンジサンプル数の関数として差を演算し、累積する。そのほかの点では、レンジフィルタ28は上述した方位フィルタ26と同様の方法で機能する。
2次元フィルタ12のサイズは上限値および下限値との間で選択することができる。この上限値は最小の病変部または被検出種類のほかの組織のサイズに等しい。下限値は平均スペックルサイズおよびコントラスト・ノイズ比(CNR)の所望の改善度に左右される。平均スペックルサイズは点広がり関数の自己相関の範囲(例えば、-6.8dB)に等しい。CNRは以下のように定義される。
Figure 2006231069
ΔIは病変部と背景との間の平均強度(後方散乱係数)の差で、σβ およびσ はそれぞれ背景と病変部の変数である。CNRを3dBだけ増大させるとすると、フィルタサイズは平均スペックルサイズに等しく設定する必要がある。
図1の如くの実施例の場合、IPはIDから発生させるが、IDの高細部分解能はIPのコントラスト分解能に寄与する。これは、細部分解能が高くなると、それにより、処理(平均)に使用できる、より独立したサンプルが得られるからである。これによる改善度は、独立サンプル数の二乗平方根に比例する。上述したように、高い細部分解能が高いコントラスト分解能にとっての十分な条件ではない。というのは、画像信号IDの高い空間周波数成分が、画像のコントラスト差を知覚するユーザの能力に不利な影響を与えるからである。
高コントラスト分解能画像信号IPは画像信号IDからビデオフィルタを使って得なければならないということは、本質的事項ではない。ほかの種類の空間フィルタリング、周波数合成、時間合成、空間合成、これらの技法の組み合わせ、さらには、ほかのスペックル除去法などの各種の技法も使用できる。また、ビーム形成法を使用して直接に画像信号IPおよびIDを発生させることが好適な場合もある。この場合、画像プロセッサにいかなる種類の2次元フィルタを含める必要がないことも多々あると想定できる。一般に、同一のビーム形成した画像信号を二度処理することは容易である(一度は細部分解能のためのIDを生成する処理で、もう一度はコントラスト分解能のためのIPを生成する処理)。
IDおよびIPは一度上述したように形成されると、出力画像信号IOを発生する表示発生器14の第1および第2の入力16および18に供給される。この信号IOはIDの細部情報とIPのコントラスト情報とを合成または重畳するもので、これにより重畳された細部およびコントラスト情報は共に保持され、信号IOがモニタ20に表示されたとき人間の目はその両方の情報を区別できる。高い空間分解能のとき、IOは IDにしたがって変化する。反対に、低い空間分解能のとき、IOは IPにしたがって変化する。
このため、出力画像信号IOは、収集され処理されて細部分解能を最大化させた信号IDと収集され処理されてコントラスト分解能を最大化させた信号IPとの合成信号になる。
この実施例では、表示発生器14は2次元のルックアップテーブルであり、その中でIDおよびIPは、テーブルから対応する格納値IOを選択するアドレスとして使用されている。IDおよびIPは、時間軸上のどの点においても、画像化する対象物の同一の物理領域に関係するように同期されている。
この実施例の表示発生器のルックアップテーブルはカラーテーブルで、以下の式(1)に示すようにIDおよびIPの両方に基づいて画素値をIOに割り当てる。
Figure 2006231069
式(1)の関数Fはベクトル関数として教示でき、このベクトル関数の各コンポーネントを表示カラーに対応させることができる。従来のRGB(赤、緑、青)モニタ20を用いたときには、Fは好適には3次元ベクトル関数である。この赤、緑、および青のコンポーネントが互いに等しく設定されたとき、IOの対応する値はグレースケール値になる。関数Fの形状を正確に設定すると、画質の多様な面を強調するように設定できる。例えば、IDが高いときでさえも、IPが低いならば、Fの値を低くすることでノイズを減らすようにFを設定することができる。IOの特定のある画素の輝度がIDとともに変わり、その一方で、その画素のカラーがIPとともに変わるようにFを設定することで、コントラスト分解能と細部分解能の両方をIPに保持させることができる。これが行われた場合、コントラストレベルが異なる2つの違う組織が異なるカラーで表示される。図5は、IDの所定値およびIPの変化値に対する3次元のFの3つのコンポーネントに対する可能な1つの形を示す。
以下の式(2)乃至(5)を使ってIOの赤、緑、および青のコンポーネントをIDおよびIPの関数として決めることができる。
Figure 2006231069
Figure 2006231069
Figure 2006231069
Figure 2006231069
以下の定数を式(2)乃至(5)に挿入することができ、ここで使用している添字R,G,Bは赤、緑、および青コンポーネントをそれぞれ示す。
Figure 2006231069
このアプローチによれば、最適化された細部分解能を有する画像信号IDを、カラーシェーディングによって、最適化されたコントラスト分解能を有する画像信号IPに重ねることができる。このように、強調されたコントラスト情報および強調された細部情報は信号IOを介してモニタ20上に同時に提示され、その両方の種類の情報がユーザによって知覚できるようになる。
図5および6は表示発生器14の別の好適な実施例に関する。この実施異例では、表示発生器は出力画像信号IOをIDおよびIPの単色の2次元マッピングとして発生させる。すなわち、
Figure 2006231069
このF’(IP,ID)を明示する形は式(7)により与えられる。すなわち、
Figure 2006231069
であり、ここでα(IP,ID)はIPおよびIDの重み付けファクタである。このため、IOはIPおよびIDの重み付けされた組み合わせ信号になる。例えば、α=0のときIO=IPであり、α=1のときIO=IDである。αを0と1の間で変化させることで、出力画像信号IOは、このような2つの限度値間でIPおよびIDのいかなる重み付けの組み合わせにも等しくすることができる。
出力画像IOの平均輝度は式(8)により与えられ、<>は空間アベレージング処理を示す。
Figure 2006231069
実際上、式(8)の第2項は零または略零となる。したがって、画像信号IOの平均輝度はIPのそれと略等しくなる。
関数αは超音波画像の周知の特性に基づいて設定することができる。例えば、IPが低い場合、αを零に近い値に設定することにより、暗い背景中のノイズを除去または実質的に減少させることができる。また暗い背景中にある小さい目標体の抑制を望まず、この抑制を回避するには、IPおよびID間の差が小さい場合、より大きい値をαに与えることができる。目標体が小さく、そのサイズがノイズスペックルのサイズよりも大きい限り、すなわち、その目標体がノイズよりも明るい限り、その目標体の(ID-IP)の値はノイズよりも小さくなる傾向にある。この理由によって、IPが小さく、(ID-IP)が大きい場合、αは好適には零に設定される。同様に、背景の強度が高く(IPが大きい)かつ(IP-ID)が大きい場合、αを零に近い値に設定して、明るい背景中の暗い穴を抑制し、スペックルノイズを減少させ、そして比較的明るい組織領域中に良好なコントラスト分解能を達成することができる。
IDおよびIPが殆ど等しい場合、IDは主に画像の細部を表し、とくにIPが中間レベルから高いレベルに及んでいる場合にそうである。この場合、αを高い値に設定してそれらの画像の細部を保持させる。
図6はα(IP,ID)の1つの好適なフォームを図示しており、このフォームは先に説明した考察に大略したがうものである。図6中の対角線状に延びているリッジの幅を調整し、またその対角線に沿ったリッジの高さを変えることにより、ノイズ抑制と細部保持とを良好にバランスさせることができる。
表示発生器14は、入力アドレスとしてのIPおよびIDを有する2次元ルックアップテーブルとしてインプリメントすることができる。これはα(IP,ID)およびその結果としてのF’(IP,ID)を事前に演算してテーブルに格納しておくことができることに因る。図7はIO=F’(IP,ID)を示し、これは図6に示すα(IP,ID)を使って計算されたものである。勿論、種々のルックアップテーブルを、特定のアプリケーションや種類の組織または患者で使用するために格納するようにしてもよい。
図6および図7に関連して先に議論した手法により、重要な利点が得られる。つまり、この手法により暗い領域中のノイズを抑制し(例えば血液プール中のクラッタ)、明るい組織領域中のスペックルノイズを減少させることができるということである。その結果、小さい構造やエッジ部分の細部を殆ど失う又は失うことがない。さらに、画像の平均強度も保持される。
表示プロセッサ10の2次元フィルタ12は全ての実施例で必要なわけではない。また、処理された画像信号IPは任意の適宜な手法で発生させて、異なる種別の組織間の有効なコントラスト情報を提供するようにもできる。
アプリケーションの中には、コントラストを高め、境界を強調するため、表示画像中の血液プールを抑制することが強く望まれる場合がある。運動情報は、以下に述べるように、時間コヒーレンス、すなわち連続するフレーム間の相関を測定することにより容易に得ることができる。この運動情報を用いて、組織の細部情報に影響を及ぼすことなく、血液プールを検出し、その血液プール画像を抑制することができる。
図8はそのような実施例の1つを示し、図1の2次元フィルタ12を運動推定器22で置換したものである。この運動推定器22は処理された画像信号IPを2次元ルックアップテーブル14’に供給するようになっている。
運動推定器22は、画像信号IDの連続するフレーム間で相関係数(すなわち、零次遅れにおける規格化された相互相関)を計算する。好適には、各フレームがサブフレームに分割されていることで、各サブフレームが20λ×20λ(λはイメージングプローブの中心周波数に対する波長)のサイズを有することである。このようなサブフレームを用いると、計算上の必要条件を減らすことができ、その一方で運動推定の有効性を確保することができる。細部情報は高分解能画像信号IDによっても提供されるので、この細部情報の喪失は重大なものにはならない。
相関係数は次いで補間され、しきい値処理されて、異なる運動速度間で区別される。ローパスフィルタリングにより異なる速度の領域間の移行が平滑化される。この結果生じる信号IDは画素毎の信号で、この信号がルックアップテーブル14’に対するアドレス用信号として使用される。このテーブル14’を適宜に選択して運動速度が速い領域(血液プールなど)を抑制することができ、これにより、組織のコントラストを増強することができる。このテーブル14’により発生される出力IOは、上述したようにマルチカラー画像または単色画像のいずれか一方をエンコードすることができる。
超音波イメージングシステムの標準動作モードの1つは、スペクトルドップラーモードである。このモードの場合、サンプルボリュームから得られるドップラーシフト周波数fdのパワースペクトラムの変動が時間の関数として表示される。時間は水平方向の軸に沿ってプロットされ、周波数が垂直方向の軸に沿ってプロットされる。(このスペクトラムの振幅は強度に変調される。)ドップラースペクトラムはサンプルボリューム内の速度分布に比例する。平坦な速度プロファイルからのスペクトラムは狭い帯域幅となり、一方、放物線状のプロファイルからのそれは広い帯域幅となる。乱流の場合、スペクトラムは時間の経過に伴って素早く変動する。
Bモード画像の増強について説明した上述の方法は同様に、スペクトラムドップラーモードの画像を改善するために用いることができる。その例の場合、第1の画像信号IDは従来のスペクトラムドップラー信号ID(fd,t)であり、第2の画像信号IPは処理されたスペクトラムドップラー信号IP(fd,t)であり、ここで、
Figure 2006231069
である。式(8)において、*はコンボリューションを意味し、h(fd)およびh(t)はそれぞれ周波数と時間のフィルタである。このフィルタh(fd)およびh(t)は増強する特徴に基づいて選ばれる。
例えば、パワースペクトラム振幅のゆっくりした変化を増強し、かつ、ノイズを減らすためには(すなわち、画像のコントラストを増加させるには)、周波数または時間のドメインにおいて、フィルタh(fd)またはh(t)はローパスフィルタとして設定される。このため、IPは式(2)または(7)を使って振幅をカラーエンコードするのに用いられる。反対に、速い変化を増強させるには、そのフィルタh(fd)およびh(t)はハイパスフィルタになるように設定できる。例えば、IPを生成するハイパスフィルタh(t)を使用し、かつ、式(2)で説明したように出力をカラーコーディングすることによって、乱流またはジェットを強調することができる。速く、時間変動し、かつ高速の流れ、または推定の変動といった特別な特徴を抽出するため、IPを生成するプロセッサは非線形のプロセッサ(または、線形フィルタおよび非線形デシジョンメーキングの組み合わせ)にも形成することができる。
この実施例では、画像信号ID(fd,t)は画像信号IP(fd,t)よりも広い帯域幅を有し、この両方の信号により描写されるエンティティ(entity)は生体組織よりむしろスペクトラムドップラーのプロットである。この結果生じる出力信号は、ID(fd,t)にしたがって変化する広い帯域幅の部分と、IP(fd,t)にしたがって変化する狭い帯域幅の部分とを有する。
上述のことから、先に説明した実施例に広範囲にわたる変形または修正を加えることができるということは明白である。F(IP,ID),F’(IP,ID)およびα(IP,ID)は全て特定のアプリケーションに適宜になるように適応させることができる。表示発生器はIPおよびIDから直接にIOを決めるようにしてもよく(例えば、計算によって)、これにより、ルックアップテーブルを用いる必要性をなくするようにしてもよい。最も広範囲に及ぶ手法を用いる場合、分離したビーム形成法またはIDからIPを発生させる手法のいずれかによって、IPおよびIDを発生させることができる。
本発明は超音波イメージングシステムに使用する場合に限定されない。そのシステムの代わりに、例えば走査型電子顕微鏡などの多様なイメージングシステムに使用することができる。加えて、本発明は物理領域または目標体以外のエンティティの画像を強調するのに使用することができる。例えば、波形やスペクトラムドップラーデータ表示などのエンティティの画像も上述した手法を使って強調することができる。したがって、ここで意図したいことは、先の説明は限定をなすものではなく、説明として見做されるということである。本発明の範囲を定義することを意図しているのは、全ての均等物を含む、以下のクレームである。
本発明の好適な実施例を具体化する超音波イメージングシステムのブロック図である。 図1に示すフィルタのブロック図である。 図2に示す方位フィルタのブロック図である。 図2に示すレンジフィルタのブロック図である。 図1に示す表示発生器の1つの好適な実施例に関するグラフである。 図1に示す表示発生器の別の好適な実施例に関するα(IP,ID)を示すグラフである。 図1に示す表示発生器の前記別の好適な実施例に関するF’(IP,ID)を示すグラフである。 本発明の別の実施例を具体化する超音波イメージングシステムのブロック図である。

Claims (1)

  1. 共通のエンティティを描写する第1および第2の画像信号に応答し、前記第1の画像信号は前記第2の画像信号よりも大きな細部分解能を有し、かつ前記第2の画像信号は前記第1の画像信号よりも大きなコントラスト分解能を有するイメージングシステムの表示プロセッサであって、この表示プロセッサは、
    前記第1の画像信号を受信する第1の入力と、
    前記第2の画像信号を受信する第2の入力と、
    前記第1および第2の入力に結合され、少なくとも前記第1および第2の画像信号の重み付けされた組み合わせとして出力信号を形成し、この重み付けされた組み合わせは前記第1および第2の画像信号の両方にしたがって変化する重み付けファクタを有するように形成した手段と、を備えたイメージングシステムの表示プロセッサ。
JP2006103636A 1995-03-10 2006-04-04 イメージングシステムの表示プロセッサ Pending JP2006231069A (ja)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US08/401,715 US5479926A (en) 1995-03-10 1995-03-10 Imaging system display processor

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP8527714A Division JPH11501841A (ja) 1995-03-10 1996-03-07 イメージングシステム用表示プロセッサ

Publications (1)

Publication Number Publication Date
JP2006231069A true JP2006231069A (ja) 2006-09-07

Family

ID=23588912

Family Applications (2)

Application Number Title Priority Date Filing Date
JP8527714A Pending JPH11501841A (ja) 1995-03-10 1996-03-07 イメージングシステム用表示プロセッサ
JP2006103636A Pending JP2006231069A (ja) 1995-03-10 2006-04-04 イメージングシステムの表示プロセッサ

Family Applications Before (1)

Application Number Title Priority Date Filing Date
JP8527714A Pending JPH11501841A (ja) 1995-03-10 1996-03-07 イメージングシステム用表示プロセッサ

Country Status (5)

Country Link
US (1) US5479926A (ja)
JP (2) JPH11501841A (ja)
AU (1) AU5093896A (ja)
DE (1) DE19681290T1 (ja)
WO (1) WO1996028095A1 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013533050A (ja) * 2010-07-29 2013-08-22 ビー−ケー メディカル エーピーエス 動き補償処理のためのシステムおよび方法
WO2013129185A1 (ja) * 2012-03-01 2013-09-06 日立アロカメディカル株式会社 超音波診断装置
JP2013180036A (ja) * 2012-03-01 2013-09-12 Hitachi Aloka Medical Ltd 超音波診断装置
JP2013180035A (ja) * 2012-03-01 2013-09-12 Hitachi Aloka Medical Ltd 超音波診断装置

Families Citing this family (79)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6005827A (en) * 1995-03-02 1999-12-21 Acuson Corporation Ultrasonic harmonic imaging system and method
US5608690A (en) * 1995-03-02 1997-03-04 Acuson Corporation Transmit beamformer with frequency dependent focus
US6009046A (en) * 1995-03-02 1999-12-28 Acuson Corporation Ultrasonic harmonic imaging system and method
US6027448A (en) * 1995-03-02 2000-02-22 Acuson Corporation Ultrasonic transducer and method for harmonic imaging
US6104670A (en) * 1995-03-02 2000-08-15 Acuson Corporation Ultrasonic harmonic imaging system and method
US5595179A (en) * 1995-05-02 1997-01-21 Acuson Corporation Adaptive persistence processing
US5544658A (en) * 1995-09-18 1996-08-13 Siemens Medical Systems, Inc. Doppler ultrasound velocity estimation
US5582173A (en) * 1995-09-18 1996-12-10 Siemens Medical Systems, Inc. System and method for 3-D medical imaging using 2-D scan data
US5653234A (en) * 1995-09-29 1997-08-05 Siemens Medical Systems, Inc. Method and apparatus for adaptive spatial image filtering
JP3542858B2 (ja) * 1995-09-29 2004-07-14 株式会社日立メディコ 超音波診断装置
US5653235A (en) * 1995-12-21 1997-08-05 Siemens Medical Systems, Inc. Speckle reduction in ultrasound imaging
US6553327B2 (en) * 1998-09-16 2003-04-22 Yeda Research & Development Co., Ltd. Apparatus for monitoring a system with time in space and method therefor
US6353803B1 (en) * 1996-01-18 2002-03-05 Yeda Research And Development Co., Ltd. At The Welzmann Institute Of Science Apparatus for monitoring a system in which a fluid flows
EP0883860B1 (en) * 1996-02-29 2006-08-23 Acuson Corporation Multiple ultrasound image registration system, method and transducer
US5724974A (en) * 1996-03-22 1998-03-10 Acuson Corporation Two-dimensional ultrasound display system
US5846202A (en) * 1996-07-30 1998-12-08 Acuson Corporation Ultrasound method and system for imaging
US5797397A (en) * 1996-11-25 1998-08-25 Hewlett-Packard Company Ultrasound imaging system and method using intensity highlighting to facilitate tissue differentiation
US6030344A (en) * 1996-12-04 2000-02-29 Acuson Corporation Methods and apparatus for ultrasound image quantification
US5882306A (en) * 1997-04-11 1999-03-16 Acuson Corporation Ultrasound imaging methods and systems
US5961460A (en) * 1997-04-11 1999-10-05 Acuson Corporation Ultrasound imaging enhancement methods and systems
US6110120A (en) * 1997-04-11 2000-08-29 Acuson Corporation Gated ultrasound imaging apparatus and method
US6050944A (en) 1997-06-17 2000-04-18 Acuson Corporation Method and apparatus for frequency control of an ultrasound system
US5913823A (en) * 1997-07-15 1999-06-22 Acuson Corporation Ultrasound imaging method and system for transmit signal generation for an ultrasonic imaging system capable of harmonic imaging
US6193659B1 (en) 1997-07-15 2001-02-27 Acuson Corporation Medical ultrasonic diagnostic imaging method and apparatus
US5833614A (en) * 1997-07-15 1998-11-10 Acuson Corporation Ultrasonic imaging method and apparatus for generating pulse width modulated waveforms with reduced harmonic response
US6023977A (en) * 1997-08-01 2000-02-15 Acuson Corporation Ultrasonic imaging aberration correction system and method
US6132374A (en) * 1997-08-01 2000-10-17 Acuson Corporation Ultrasonic imaging method and system
US5860928A (en) * 1997-08-07 1999-01-19 Siemens Medical Systems, Inc. Method and system for selectively smoothing color flow images in an ultrasound system
US6312379B1 (en) 1997-08-15 2001-11-06 Acuson Corporation Ultrasonic harmonic imaging system and method using waveform pre-distortion
US5944666A (en) * 1997-08-21 1999-08-31 Acuson Corporation Ultrasonic method for imaging blood flow including disruption or activation of contrast agent
US5928151A (en) * 1997-08-22 1999-07-27 Acuson Corporation Ultrasonic system and method for harmonic imaging in three dimensions
US6106465A (en) * 1997-08-22 2000-08-22 Acuson Corporation Ultrasonic method and system for boundary detection of an object of interest in an ultrasound image
US5873830A (en) * 1997-08-22 1999-02-23 Acuson Corporation Ultrasound imaging system and method for improving resolution and operation
US5876343A (en) * 1997-09-23 1999-03-02 Scimed Life Systems, Inc. Methods and apparatus for blood speckle detection in an intravascular ultrasound imaging system
US5860931A (en) * 1997-10-10 1999-01-19 Acuson Corporation Ultrasound method and system for measuring perfusion
US5935069A (en) * 1997-10-10 1999-08-10 Acuson Corporation Ultrasound system and method for variable transmission of ultrasonic signals
US5897500A (en) * 1997-12-18 1999-04-27 Acuson Corporation Ultrasonic imaging system and method for displaying composite fundamental and harmonic images
US5902242A (en) * 1998-01-22 1999-05-11 Acuson Corporation System and method for forming a combined ultrasonic image
US6511426B1 (en) 1998-06-02 2003-01-28 Acuson Corporation Medical diagnostic ultrasound system and method for versatile processing
US6116244A (en) * 1998-06-02 2000-09-12 Acuson Corporation Ultrasonic system and method for three-dimensional imaging with opacity control
US5957852A (en) * 1998-06-02 1999-09-28 Acuson Corporation Ultrasonic harmonic imaging system and method
US6126599A (en) * 1998-10-01 2000-10-03 Atl Ultrasound, Inc. Ultrasonic diagnostic imaging system with real time spatial compounding processor
US6117081A (en) * 1998-10-01 2000-09-12 Atl Ultrasound, Inc. Method for correcting blurring of spatially compounded ultrasonic diagnostic images
US6048316A (en) * 1998-10-16 2000-04-11 Acuson Corporation Medical diagnostic ultrasonic imaging system and method for displaying composite fundamental and harmonic images
US7837624B1 (en) 1998-11-20 2010-11-23 Siemens Medical Solutions Usa, Inc. Medical diagnostic ultrasound imaging methods for extended field of view
US6554770B1 (en) * 1998-11-20 2003-04-29 Acuson Corporation Medical diagnostic ultrasound imaging methods for extended field of view
US6042545A (en) * 1998-11-25 2000-03-28 Acuson Corporation Medical diagnostic ultrasound system and method for transform ultrasound processing
US6120446A (en) * 1998-12-17 2000-09-19 Acuson Corporation Diagnostic medical ultrasonic imaging system and method with adaptive gain
US6213946B1 (en) 1998-12-24 2001-04-10 Agilent Technologies, Inc. Methods and apparatus for speckle reduction by orthogonal pulse compounding in medical ultrasound imaging
US6322505B1 (en) 1999-06-08 2001-11-27 Acuson Corporation Medical diagnostic ultrasound system and method for post processing
US6484191B1 (en) 1999-07-02 2002-11-19 Aloka Co., Ltd. Apparatus and method for the real-time calculation of local variance in images
US6454715B2 (en) 2000-04-11 2002-09-24 Scimed Life Systems, Inc. Methods and apparatus for blood speckle detection in an intravascular ultrasound imaging system
US6726631B2 (en) 2000-08-08 2004-04-27 Ge Parallel Designs, Inc. Frequency and amplitude apodization of transducers
US6416477B1 (en) 2000-08-22 2002-07-09 Koninklijke Philips Electronics N.V. Ultrasonic diagnostic systems with spatial compounded panoramic imaging
US6464638B1 (en) * 2000-10-05 2002-10-15 Koninklijke Philips Electronics N.V. Ultrasound imaging system and method for spatial compounding
US6663567B2 (en) * 2002-03-19 2003-12-16 Zonare Medical Systems, Inc. System and method for post-processing ultrasound color doppler imaging
US6858010B2 (en) * 2003-05-06 2005-02-22 Siemens Medical Solutions Usa, Inc. Identifying clinical markers in spatial compounding ultrasound imaging
KR101089426B1 (ko) * 2003-09-11 2011-12-07 파나소닉 주식회사 시각 처리 장치, 시각 처리 방법, 시각 처리 프로그램 및 반도체 장치
KR20060120643A (ko) * 2003-09-11 2006-11-27 마쯔시다덴기산교 가부시키가이샤 화상처리 장치, 화상처리 방법 및 화상처리 프로그램
EP1667066B1 (en) * 2003-09-11 2020-01-08 Panasonic Intellectual Property Corporation of America Visual processing device, visual processing method, visual processing program, integrated circuit, display device, imaging device, and mobile information terminal
CN100417175C (zh) * 2003-09-11 2008-09-03 松下电器产业株式会社 视觉处理装置、视觉处理方法和半导体器件
KR101089394B1 (ko) * 2003-09-11 2011-12-07 파나소닉 주식회사 시각 처리 장치, 시각 처리 방법, 시각 처리 프로그램 및 반도체 장치
US20090216124A1 (en) * 2005-05-19 2009-08-27 Hitachi Medical Corporation Ultrasonic diagnostic apparatus and image processing method thereof
US7983456B2 (en) * 2005-09-23 2011-07-19 Siemens Medical Solutions Usa, Inc. Speckle adaptive medical image processing
WO2008010375A1 (fr) * 2006-07-20 2008-01-24 Hitachi Medical Corporation Dispositif ultrasonographique
US7780601B2 (en) 2007-06-05 2010-08-24 Siemens Medical Solutions Usa, Inc. Adaptive clinical marker preservation in spatial compound ultrasound imaging
JP2009011711A (ja) * 2007-07-09 2009-01-22 Toshiba Corp 超音波診断装置
JP5441795B2 (ja) * 2010-03-31 2014-03-12 キヤノン株式会社 イメージング装置及びイメージング方法
GB2493396B (en) * 2011-08-05 2018-08-08 Snell Advanced Media Ltd Multidimensional sampled-signal filter
GB2502047B (en) * 2012-04-04 2019-06-05 Snell Advanced Media Ltd Video sequence processing
US8696577B2 (en) * 2012-04-20 2014-04-15 Siemens Medical Solutions Usa, Inc. Tongue imaging in medical diagnostic ultrasound
US9081097B2 (en) 2012-05-01 2015-07-14 Siemens Medical Solutions Usa, Inc. Component frame enhancement for spatial compounding in ultrasound imaging
KR102025756B1 (ko) * 2013-01-04 2019-09-27 삼성전자주식회사 영상에서 스펙클을 제거하는 방법, 장치 및 시스템.
CN110338751A (zh) 2013-03-15 2019-10-18 瓦索普蒂克医疗公司 使用多种照明形态的眼科检查及疾病管理
CN104323794B (zh) * 2013-07-22 2016-12-28 华中科技大学 一种用于三维超声成像的拆分式行列寻址方法
EP2989986B1 (en) * 2014-09-01 2019-12-18 Samsung Medison Co., Ltd. Ultrasound diagnosis apparatus and method of operating the same
US10430688B2 (en) 2015-05-27 2019-10-01 Siemens Medical Solutions Usa, Inc. Knowledge-based ultrasound image enhancement
CN109152564A (zh) * 2015-12-31 2019-01-04 皇家飞利浦有限公司 自动血池识别***及其操作方法
US11185305B2 (en) * 2016-06-30 2021-11-30 Koninklijke Philips N.V. Intertial device tracking system and method of operation thereof

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB8431374D0 (en) * 1984-12-12 1985-01-23 Bamber J C Adaptive filtering
US4787393A (en) * 1985-11-20 1988-11-29 Matsushita Electric Industrial Co., Ltd. Ultrasonic tomographic with alternate image scaling
JPH03224552A (ja) * 1990-01-31 1991-10-03 Toshiba Corp 超音波診断装置
JP2938123B2 (ja) * 1990-03-30 1999-08-23 株式会社東芝 多機能デジタルカメラ
US5224483A (en) * 1992-05-19 1993-07-06 Hewlett-Packard Company Adaptive contrast enhancement for scanned ultrasonic image

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013533050A (ja) * 2010-07-29 2013-08-22 ビー−ケー メディカル エーピーエス 動き補償処理のためのシステムおよび方法
WO2013129185A1 (ja) * 2012-03-01 2013-09-06 日立アロカメディカル株式会社 超音波診断装置
JP2013180036A (ja) * 2012-03-01 2013-09-12 Hitachi Aloka Medical Ltd 超音波診断装置
JP2013180035A (ja) * 2012-03-01 2013-09-12 Hitachi Aloka Medical Ltd 超音波診断装置
US9514521B2 (en) 2012-03-01 2016-12-06 Hitachi, Ltd. Diagnostic ultrasound apparatus

Also Published As

Publication number Publication date
US5479926A (en) 1996-01-02
AU5093896A (en) 1996-10-02
WO1996028095A1 (en) 1996-09-19
DE19681290T1 (de) 1998-03-19
JPH11501841A (ja) 1999-02-16

Similar Documents

Publication Publication Date Title
JP2006231069A (ja) イメージングシステムの表示プロセッサ
US5961460A (en) Ultrasound imaging enhancement methods and systems
US6679843B2 (en) Adaptive ultrasound image fusion
US8021300B2 (en) Three-dimensional fly-through systems and methods using ultrasound data
Michailovich et al. Despeckling of medical ultrasound images
JP2558108B2 (ja) 影像ろ過装置及び影像ろ過方法
JP4757307B2 (ja) 超音波画像処理装置
US5931784A (en) Ultrasonic diagnostic apparatus
KR101205107B1 (ko) 스페클 감소 필터의 구현 방법, 스페클 감소 필터링 장치 및 초음파 촬상 시스템
KR100742466B1 (ko) 이미징 방법 및 장치와 그를 이용한 혈액 움직임 이미징 시스템
US5793883A (en) Method for enhancing ultrasound image
US6508766B2 (en) Ultrasound diagnostic apparatus
US6036643A (en) Ultrasonic harmonic doppler imaging
JP2002177275A5 (ja)
US9123139B2 (en) Ultrasonic image processing with directional interpolation in order to increase the resolution of an image
CN110832343B (zh) 具有改善的杂波抑制的能量多普勒成像***和方法
US20120101384A1 (en) Ultrasonic diagnostic device
JP2011087965A (ja) 空間合成による超音波診断造影画像
Bamber et al. Real-time implementation of coherent speckle suppression in B-scan images
JP4732345B2 (ja) 超音波撮像方法および装置
Lovstakken et al. Blood flow imaging-a new real-time, flow imaging technique
JP2004154567A (ja) 画像処理装置、方法およびプログラム
US5865752A (en) Method and apparatus for ultrasound imaging using normalized difference between successive frames
JP4808373B2 (ja) Bモード画像のバンディング抑制に関連する応用のための方法及び装置
JP2950632B2 (ja) 超音波診断装置

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20061114

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20070306

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20081222