(57)【特許請求の範囲】
【請求項１】
線形および非線形超音波エコー信号成分を微分的に位相変調する超音波診断画像処理装置の動作方法であって、
ビーム形成器制御器が発信ビーム形成器を制御して、発生する超音波エコーに位相差を生じさせるように、体内に変調された超音波シーケンスを発信する段階、
該ビーム形成器制御器が受信ビーム形成器を制御して、該発信シーケンスに応答した超音波エコー信号の集合を受信する段階、および、
線形および非線形信号成分の位相シフト情報を分離するために該集合を分析する段階、
を含む方法。
【請求項２】
該発信する段階が、パルス毎に変調された超音波の発信を含む請求項１の方法。
【請求項３】
該パルス毎に変調された超音波の発信が、線形信号成分と非線形信号成分の間に位相の分離を生じさせる請求項２の方法。
【請求項４】
該変調が、該発信パルス速度のナイキスト限界（最大測定可能速度）で生じる請求項２の方法。
【請求項５】
該変調が、パルス毎の１８０゜位相変化、パルス毎の振幅変調、またはパルス毎の位相および振幅変調の組合せである請求項２の方法。
【請求項６】
該分析する段階が、線形および非線形の信号成分を微分的に検出する段階を含む請求項１の方法。
【請求項７】
微分的に検出する段階が、該集合をパワードップラー処理することを含む請求項６の方法。
【請求項８】
微分的に検出する段階が、該集合を色流れ処理することを含む請求項６の方法。
【請求項９】
線形および非線形信号情報を分離するように、ドップラー領域に無線周波数超音波エコー信号情報をマッピングする超音波診断画像処理装置の動作方法であって、
ビーム形成器制御器が発信ビーム形成器を制御して、位相変化超音波のシーケンスを発信する段階、
該ビーム形成器制御器が受信ビーム形成器を制御して、該位相変化超音波に応答した超音波エコー信号の集合を受信する段階、および、
エコー信号の該集合をドップラー処理し、エコー信号の線形性に従って分離されたドップラー位相シフト情報を作成する段階、
を含む方法。
【請求項１０】
線形および非線形の超音波信号成分を分離する超音波診断画像処理装置の動作方法であって、
ビーム形成器制御器が発信ビーム形成器を制御して、一つの波から次へと異なって振幅変調された超音波シーケンスを体内に発信する段階、
該ビーム形成器制御器が受信ビーム形成器を制御して、該振幅変調超音波に応答した超音波エコー信号の集合を受信する段階、および、
エコーの該集合の位相シフト特性に基づいて、受信超音波エコー信号を処理し、分離された線形および非線形信号成分を作成する段階、
を含む方法。
【請求項１１】
線形および非線形超音波エコー信号成分を識別する超音波診断画像処理装置であって、
被検者に超音波の変調シーケンスを発信する超音波発信器、
超音波の該変調発信シーケンスに応答した超音波エコー信号の集合を受信する受信器、および、
該エコーの位相シフト情報を線形および非線形信号成分に分離するプロセッサ、
を含む超音波診断画像処理装置。
【請求項１２】
該発信器が、パルス毎に変調された超音波を発信する請求項１１の超音波診断画像処理装置。
【請求項１３】
該発信パルス毎の変調が、線形信号成分と非線形信号成分の間に位相分離を導入することを含む請求項１２の超音波診断画像処理装置。
【請求項１４】
該変調が、該発信パルス速度(rate)のナイキスト限界で生じる請求項１２の超音波診断画像処理装置。
【請求項１５】
該変調が、パルス毎の１８０゜位相変化、パルス毎の振幅変調、またはパルス毎の位相および振幅変調の組合せを含む請求項１２の超音波診断画像処理装置。
【請求項１６】
該プロセッサが、線形および非線形信号成分を微分的に検出する微分プロセッサを含む請求項１１の超音波診断画像処理装置。
【請求項１７】
該プロセッサが、パワードップラープロセッサを含む請求項１６の超音波診断画像処理装置。
【請求項１８】
該プロセッサが、色流れプロセッサを含む請求項１６の超音波診断画像処理装置。
【請求項１９】
線形および非線形の超音波信号情報を分離する超音波診断画像処理装置であって、
被検者に位相変化超音波のシーケンスを発信する超音波発信器、
該位相変化超音波に応答した超音波エコー信号の集合を受信する受信器、および、
エコー信号の該集合を処理し、該超音波エコー信号の線形性に従って描出されたドップラー位相シフト情報を作成するドップラープロセッサ、
を含む超音波診断画像処理装置。
【請求項２０】
動いている組織の高調波信号効果を識別する超音波診断画像処理装置であって、
被検者に位相変化超音波のシーケンスを発信する発信器、
該超音波に応答したエコーの集合を受信する受信器、および、
該集合を処理し、動いている組織の奇数次および偶数次高調波信号の効果のスペクトルを識別する組織ドップラープロセッサ、
を含む超音波診断画像処理装置。
【請求項２１】
動いている組織または高調波造影剤の高調波信号効果を識別する超音波診断画像処理装置であって、
被検者に位相変化超音波のシーケンスを発信する発信器、
該超音波に応答したエコーの集合を受信する受信器、および、
該集合を処理し、動いている組織または高調波造影剤の奇数次および偶数次高調波信号効果のスペクトルを識別するパワードップラープロセッサ、
を含む超音波診断画像処理装置。
【請求項２２】
奇数次および偶数次高調波の超音波信号情報を分離する超音波診断画像処理装置であって、
被検者に位相変化超音波のシーケンスを発信する超音波発信器、
該位相変化超音波に応答した超音波エコー信号の集合を受信する受信器、
エコー信号の該集合を処理し、ドップラー位相シフト情報を作成するドップラープロセッサ、および、
奇数次高調波ドップラー位相シフト情報および偶数次高調波ドップラー位相シフト情報を、それぞれ通過させる第１および第２フィルタ、
を含む超音波診断画像処理装置。
【請求項２３】
さらに、奇数次および偶数次高調波位相シフト情報の比を計算するための該第１および第２のフィルタに結合する比回路を含む請求項２２の超音波診断画像処理装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、超音波診断画像処理の方法および装置に関し、特にパルス反転ドップラーの新技術による線形(linear)と非線形(nonlinear)超音波散乱物体(scatterer)からのエコーの検出と処理に関する。
【０００２】
【従来の技術】
近年、従来の線形信号に代わって非線形信号を利用する２つの新しい診断超音波画像処理技術が出現した。従来の超音波画像処理においては、発信波の反射したエコーが受信され、検出され、画像を形成する。これらのエコーは、音響エネルギーの吸収により生じる振幅減衰のみが異なるが、周波数と波形を含む実質的にあらゆる点で発信波形を写しているので線形エコーと呼ばれる。エコーの振幅減少の対応は、エコーが発信波形と一次関数的に関係していることを意味し、このため線形画像処理なる用語で呼ばれる。
【０００３】
受信エコーは、いくつかの出所からの非線形信号成分も含む場合がある。これらの構成部分は、発信波形に対して非線形関係にある。非線形効果の１つの出所は、発信波形が組織を通過するときに受ける歪みである。音波が組織を通って伝わる中で、波の圧縮(compressional)位相の速度は、 希薄(rarefactional)位相の速度を越える。透過波の二相の速度の相違は、組織中を進行するにつれて段階的に波形の歪みを引き起こし、元の波の周波数の高調波成分を生じさせる。これらの歪められた波形が反射するとき、戻るエコーにはこれらの歪みが含まれ、受信エコーには基本周波数と、基本周波数のより高次の高調波エネルギーが含まれる。これら自然に発生する発信波形の高調波の検出と画像処理は、組織高調波画像処理と呼ばれる。
【０００４】
非線形エコー構成部分の第２の出所は、微小気泡(microbubble)造影剤である。今日の新世代の造影剤は、十分なエネルギー準位で照射(insonified)されるとき、不規則で非線形に共鳴し、従って非線形音響エネルギー応答の原因となる。基本周波数エネルギーと同様にこの非線形エネルギーは、それが受信される変換器に戻り、加工され、表示される。さらに高いエネルギー準位が微小気泡に照射されると、微小気泡は照射された音波により破壊され、パルス毎に明確に変化した強いエコーを作成し、容易に検出し、表示することができる。いずれの場合にも、これらの強い非線形な微小気泡の応答は、微小気泡を含む体液の非常に明瞭な画像の形成に使用することができる。
【０００５】
【発明が解決しようとする課題】
線形と非線形の信号を分離する従来の方法は、フィルタを使用する周波数を基礎とするものである。基本および高調波信号を含む受信エコーは、低域通過または帯域通過フィルタ処理され、基本周波数信号を通過させ、第二高調波信号を排除することができる。信号は、高域通過または、より高い周波数帯域通過フィルタ処理し、基本周波数を除外して第二高調波信号を通過させることもできる。しかし、基本および高調波エネルギーのフィルタ分離は、使用される帯域通過フィルタの通過帯域により帯域制限され、基本と高調波エネルギーに周波数の重複部分があると、完全に有効とならない。
【０００６】
フィルタの必要性とそれらの制限条件を回避する、大きく改良された基本と高調波信号の分離技術がアメリカ合衆国特許第５,７０６, ８１９号に記載されている。本発明者の一人も共同発明者であるこの特許は、信号処理によって基本と、高調波の信号成分を分離する「パルス反転(pulse inversion)」と呼ばれる技術を開示する。パルス反転法においては、２つの反対位相の波形が発信され、各発信後にエコーが受信される。２つのエコーのシーケンスは、次いで空間基準で加算されるか減算される。ある位置で受信されたエコーが加算されると、発信波形と線形関係にあり、位相のみ相違する線形成分は、それらの位相が反対であるため消去される。連続パルスからの非線形エコー成分の大部分は、反対位相ではなく、加算的に組合わされる。対応して、ある位置からの２つのエコーが減算されると、反対の位相の基本（または線形）成分は互いに増強され、一方非線形成分の大部分は減算によって消去される。
【０００７】
パルス反転分離がフィルタに依存しないので、基本と高調波成分の周波数が重なっている時にも、広帯域信号の分離に有効である。しかしながら、パルス反転が２−パルス方式であるので、運動虚像の対象となりがちである。特に応答がいずれの場合にも作成されるので、非線形エコーと運動とをパルス反転法で識別することができない。従って、エコーの線形と非線形成分のドップラーシフトを検出し、分離できることが望ましい。さらに線形と、非線形の信号成分を十分に分離できることが望ましく、さらに非線形エコー成分を運動の効果から識別できることが望ましい。そのような信号分離は、広帯域操作で、ｒｆ（無線周波数）フィルタ回路の制約によって妨げられないことが望ましい。
【０００８】
【課題を解決するための手段】
本発明により、運動の効果を識別しつつ、同時に広帯域の線形および非線形エコー信号成分を分離する超音波診断画像処理方法および装置が提供される。本発明の方法では、無線周波数（時間）領域の超音波エコー信号を、ドップラー領域のドップラーシフト周波数に、エコーの直線性に従って、マッピングする。本発明方法では、変調(modulated)発信シーケンスに応答して受信された連続エコーの間の位相シフトを分析することによりこれを実施する。本発明方法によると、周期的に位相変化する波形のシーケンスの１つの発信により１つのドップラーエコー集合が取得される。ドップラー処理中に、エコーの集合が組合わされると、結果として生ずるドップラースペクトルは、１つの領域に存在する奇数次高調波（例えば線形信号成分）と、その他の領域に存在する偶数次高調波（例えば第二高調波信号成分）の、別個の領域に分けられる。各領域は、運動の存在なしで体内の静止している物質からの線形および非線形エコーに含まれる画像化される物質の運動に対応する周波数範囲に広がる。好適例においては、組織の運動、高調波造影剤、または組織高調波画像処理情報などの、特定種類の処理用信号情報を同時に描出するのにドップラーフィルタが使用される。本発明方法を実施するための装置についても説明する。
【０００９】
図１は、本発明による超音波診断画像処理装置のブロックダイヤグラムである。図２は、本発明によるエコー集合を取得するための好適な発信波形シーケンスを示す。図３は、図１の集合記憶装置の説明図である。図４ａと４ｂは、図１のウォ−ルフィルタの２つの例を示す。図５ａは、従来のドップラースペクトルを説明し、図５ｂは本発明により作成されたドップラースペクトルを示す。図６は、本発明による異なる高調波スペクトルの分離を示す。図７−１０は、本発明により分析され、画像処理される超音波情報の種々のスペクトルを示す。図１１−１５は、図７−１０のスペクトルの種々の超音波情報を描出するためのウォ−ルフィルタ特性を示す。
【００１０】
【発明の実施の態様】
最初に図１に、本発明の超音波診断画像処理装置のブロックダイヤグラムを示す。アレー変換器１２を有するプローブ１０は、体内に超音波エネルギーを発信し、体内組織、細胞、体液、超音波造影剤を使用する場合には造影剤から戻るエコーを受信する。アレー変換器は線形またはカーブドアレーであり、フェーズドアレイまたは線形アレイとして動作する。ドップラー研究には通常、フェーズドアレイ操作が望ましい。アレー変換器による発信と受信のタイミングは、発信ビーム形成器１６と受信ビーム形成器１８に接続するビーム形成器制御器１４により同調される。各ビーム形成器のチャネルはアレー変換器の個々の素子に接続され、個々の素子からの信号の発信と受信を別々に制御する。ビーム形成器制御器の制御下に、ビーム形成器１６はアレーの各素子が活性化され波もしくはパルスを発信する時期を決定する。この発信の時間制御は、全アレーにより発信される波が、所定の方向、即ち前もって決められた走査線に沿って操作され、所望の焦点深度に焦点を合わせることを可能にする。２つのビーム形成器のチャネルは、高発信電圧から受信ビーム形成器チャネルの入力を保護する発信／受信スイッチ１７により、アレーの素子に接続する。
【００１１】
個々の変換器素子によって受信されたエコーは、発信／受信スイッチ１７により受信ビーム形成器１８の個々のチャネルに転送される。これらの入力経路には、受信エコー信号を増幅するための前置増幅器(preamplifier)や、深さ依存減衰効果を相殺するための時間利得代償回路を含めてもよい。好適例にあるように、受信ビーム形成器１８がデジタルビーム形成器であるとき、ビーム形成器の各チャネルにはアナログ−デジタル変換器が前に置かれるかまたは含まれる。ビーム形成器のチャネルは、走査線からの各変換器素子により受信されたエコーを連続的に適切に遅延させ、走査線に沿って共通の点（試料容積(sample volume)）から受信された信号を時間一致させる。連続的な遅延変位により、走査線に沿った受信エコー信号の動的焦点合わせが可能となる。次いでチャネルの出力での信号が組合わされ、コヒレントエコー信号のシーケンスが形成され
る。
【００１２】
受信ビーム形成器は、動的開口変化の効果を相殺する信号振幅の正規化など、他の処理操作も従来通り実行する。受信ビーム形成器は、複数線受信を実施するために、それぞれがそれ独自のビーム形成器制御器からの遅延プログラミングを有するチャネルの２あるいはそれ以上の群に分割することもできる。複数線受信において、チャネルの各群は、それ自身の走査線に沿って受信ビームを操作し、焦点を合わせ、これによって２あるいはそれ以上の受信走査線を同時に形成する。複数線受信により、より短時間で一画像の全走査線が取得されるので、フレーム速度を上げることができるが、一方、受信走査線の全てが発信開口と一致しているとは限らないので、収差(aberration)効果を生じる傾向が出てくる。
【００１３】
走査線に沿って受信されたコヒ−レントエコーのシーケンスを検出することができ、グレースケール値の範囲に調節され、所望の画像記録方式に走査変換され、表示され、このようにしてＢモード画像を形成する。図１の装置において、コヒレントエコーは直角位相帯域通過（ＱＢＰ）フィルタ２０により、同相（Ｉ）と直角位相（Ｑ）試料に復調される。このＩ，Ｑ試料はドップラー処理され、ドップラーパワー、速度、加速式、分散、および流れまたは運動の方向を決定し、計算式、(Ｉ２＋Ｑ２)１/２によりエコー信号の振幅の検出に使用することもできる。図１の具体例におけるＢモード画像処理では、Ｉ，Ｑ試料は処理され、周波数混合回路２２によりスペックルを除去される。エコー振幅は、検出器２４により検出され、検出されたエコー信号は、エコーデータバス３２を通じてグレースケールプロセッサ６０に転送され、そこでエコーは対数圧縮され、グレースケールマッピングされる。対数圧縮とグレースケール処理の好適技術の詳細は、アメリカ合衆国特許出願第０８/９５５,８１９号に見出される。グレースケール信号は、画像データバス８２を通じて走査変換器８０に転送され、そこでＲ−θ走査線データは所望の表示形式に変換される。走査変換画像は、表示装置９０に表示される。
【００１４】
ＱＢＰフィルタ２０、周波数混合回路２２および検出器２４の詳細な説明は、アメリカ合衆国特許出願第０８/８９３,４２６に見出される。
ドップラー画像処理については、走査線はある時間間隔で連続的に走査され、走査線に沿った各試料容積で、時間的にエコーシーケンスが集められる。この時間的エコーシーケンスは、集合(ensemble)と呼ばれ、発信波のシーケンスにより取得され、その繰返し周波数はパルス繰返周波数、ＰＲＦと呼ばれる。それぞれの個々に発信された波もしくはパルスは、診断用超音波の通常の無線周波数(r.f.)範囲にあるドップラー周波数と呼ばれる名目上の周波数を示す。ＰＲＦは通常、キロヘルツかそれ以下の範囲である。エコー集合は集合記憶装置２６に蓄積され、完成された集合がドップラー処理用に作成される。
【００１５】
ドップラー処理の従来からの第一段階は、ウォ−ルフィルタ処理である。心臓と血管中の血流の画像処理もしくは測定を実施するとき、血液細胞からの相対的に低いレベルのエコーは、血管や心臓壁などの近接組織により反射された強いエコーに圧倒される場合がある。処理の目的は血流の画像処理もしくは測定であるから、この場合には組織エコーは、除去すべきクラッタとなる。このフィルタの目的が心臓や血管壁からのエコーを除去することにあるので、これらの不要な信号を除去する回路はウォ−ルフィルタと呼ばれている。組織信号は通常血流信号よりも大きい振幅と低い周波数を有するので、これらの信号は振幅、周波数、または両者の組合せによって識別することができる。「フラッシュ(flash)」として知られるドップラー虚像と共に組織信号を除去する好適技術がアメリカ合衆国特許第５,１９７,４７７号に開示されている。
【００１６】
ウォ−ルフィルタは、血流ドップラー信号を排除し、組織のドップラー信号を通すよう、逆特性で機能させることもできる。組織のこれらの信号がドップラー処理されるとき、心筋や弁膜などの動いている組織の画像が作成される。この画像処理技術は組織ドップラー画像処理として知られている。
【００１７】
血流または組織の、フィルタ処理されたドップラー信号は、ドップラープロセッサ４０に転送され、そこでそれらはドップラー位相シフトまたは信号強度（パワードップラー）のドップラー計算の実施に使用される。従来よりこれはドップラー信号データのフーリエ変換または自己相関(autocorrelation)により実施される。好適技術は、ドップラー位相シフトとドップラー信号の参照もしくは中心周波数を同時に計算する二次元自己相関の実施である。後者は、位相シフト計算における深さ依存周波数減衰効果の修正に有用である。そのような二次元ドップラープロセッサが、アメリカ合衆国特許５，３８６，８３０に記載されている。ドップラー周波数もしくは位相シフトは、エコーを返す血流または組織の速度に比例するので、速度、加速または分散数の延長(production)は、直線的である。色流れドップラーにおいて、血流速度は色スケールにマッピングされ、画像データバス８２を通じて走査変換器８０に転送され、血流を含む組織構造のグレースケール画像上に重ねられる。パワードップラー画像処理において、ドップラー信号強度は同様にマッピングされ、グレースケール画像上に表示される。アメリカ合衆国特許５,４７４,０７３および５,７２０,２９１に記載されているように、ドップラーとグレースケール画像データを３Ｄプロセッサ７０で処理し、血流および／または組織の三次元画像表示を形成することもできる。
【００１８】
本発明によると、エコー集合は波毎に異なる位相を示す発信波のシーケンスにより取得される。好適例において位相は、図２のパルスシーケンスにより示されるように波毎に１８０゜位相変化する。説明の容易さのために、画像処理の目的に応じ、より長いパルス周期もしくは単極性パルスを採用することもできるが、これらのパルスはサイン波形の一周期として示される。図中、パルス１０２が時間Ｔ１で発信され、その後走査線からのエコーのシーケンスが取得され、各集合中の第１の試料となる。時間Ｔ１からパルス速度(rate)間隔（ＰＲI）と呼ばれる間隔で隔てられている時間Ｔ２で、パルス１０２の反転として示されるパルス１０４が発信され、パルス １０４は、パルス１０２の有する位相から１８０゜変化する。同様に、時間Ｔ３、Ｔ４およびＴ５におけるパルス１０６,１０８および１１０の位相は、パルス毎に１８０゜変化する。
【００１９】
各パルスの後にエコーが受信され、図３に示すように他の試料が各集合に加えられ、それは集合記憶装置２６の目的を説明している。発信パルス１０２に応答して、走査線に沿った連続試料容積（ＳＶ#）から受信されたエコーは、集合記憶装置の第１の列に、Ｔ１と記された列に＋記号で示されるように記憶される。発信パルス１０４に応答して、走査線に沿った試料容積から受信されたエコーは、集合記憶装置の第２の列にＴ２と記された列に−記号で示されたように記憶される。列Ｔ３−Ｔ５の残りは、パルス１０６,１０８および１１０に応答して受信されるエコーにより満たされる。
【００２０】
各試料容積ＳＶの集合は次いで、ウォ−ルフィルタ３０中で、ウォールフィルタ処理される。好適例において、ウォ−ルフィルタはマトリックス(matrix)フィルタである。マトリックスウォールフィルタ処理は、適切な積算係数を選択して実施され、該マトリックスフィルタ処理は下式で表わされる。
【００２１】
【数１】
0004365909.tif 000002
【００２２】
ここに、ｘ１・・・・ｘnは試料の集合であり、 ｙ１・・・・・ｙnはフィルタ処理された出力値である。エコー値と係数の積ｙnは、係数により決定された特性に従ってフィルタ処理されたエコー値である。係数は、所望のフィルタ特性、特に所望のフィルタインパルス応答に関するものを考慮して選択される。
【００２３】
本発明のさらなる面によると、マトリックスウォ−ルフィルタの組合せ処理は、エコー信号の線形と非線形成分を分離する。以下に論ずるように、この分離は発信信号の位相交替によりドップラー値に導入される位相シフトによりもたらされる。
【００２４】
図４ａと４ｂは、ウォ−ルフィルタ３０の異なる２つの例を示す。図４ａの例は、ウォ−ルフィルタとしての使用に適する有限インパルス応答（ＦＩＲ）フィルタを示す。集合のエコーは、一連の遅延段階τにシフトされる。集合が遅延段階を通過するとき、試料は、マルチプライヤ１２０で重み付け係数を積算され、信号に所望のフィルタ特性、およびフィルタ出力信号のシーケンスを作成する加算段階１２２で加算された結果を与える。図４ｂは、マルチプライヤ１２０とアキュムレータ１２４により形成されたフィルタを示す。この例において、集合のエコーには順次マルチプライヤ１２０に適用される重み付け係数が積算され、その積はアキュムレータ１２４に蓄積される。完全な出力値がアキュムレータに蓄積される前に、一連の積算命令が連続して行なわれなければならないので、マルチプライヤ−アキュムレータの例は縦列し、より高速の出力信号を作成することができる。これらのフィルタの例のいずれかの出力は、定期的に跳び越す(skip)ことができ、得られた出力に間引き(decimation)フィルタ特性を与える。集合が有限長さのシーケンスなので、処理を開始するときにフィルタが「振動的に歪む(ring)」ことのないよう、初期化に注意しなければならない。これは直流値を使用するか、または、シーケンスのいくつかの初期値を複製して、エコー値の前に初期化値を加える(pad)ことにより、実施することができる。有限フィルタ長はまた、フィルタ応答がフィルタ処理される集合長によって影響されることを意味し、試料のより長いシーケンスは、より多くの係数、重み付け係数、および積を有し、より複雑で、詳細に調整されたフィルタ特性を与える。
【００２５】
従来のドップラー処理により作成されたドップラースペクトルを図５ａに示す。ドップラー応答が発信パルスによって効果的にサンプリングされるので、パルス繰返周波数は、エイリアシングなしで検出することができるドップラー周波数の範囲としてナイキスト限界(Nyquist limit; 最大測定可能速度))を設定する。運動または流れは変換器に向かう（＋）か、離れる（−）かのいずれかであるから、スペクトルの限界は、−(ＰＲＦ／２）と＋(ＰＲＦ／２)であり、ゼロの周波数基準に中心がある。縦座標は、ドップラーパワーの大きさである。ドップラーシフト周波数は、運動の速度と直接対応するので、横座標はドップラー周波数と速度の両方に置くことができる。静止物体から戻るドップラー信号は、中央ゼロの速度軸に位置するであろう。動いている血液または組織からの信号のドップラーシフトは、動きの方向により、より遅い速度はより中央軸に近く、より大きい速度は横座標の末端よりに位置付けるように、中央軸の左または右に位置する。
【００２６】
本発明の実施により生ずるドップラースペクトルを図５ｂに示す。従来のスペクトルと同様、図５ｂのドップラースペクトルは−(ＰＲＦ／２)と＋(ＰＲＦ／２)のナイキスト限界により境界づけられ、ゼロの速度座標付近に中心がある。しかしながら、パルス毎の発信信号の位相交替は、ドップラーシフトスペクトルを、奇数（線形＋より高次の非線形）と偶数（第二高調波＋より高次の非線形）成分領域に区分し、静止組織から戻るものなどの静的線形成分の、スペクトルの両端への変位を生じさせる。２次高調波成分が存する非線形領域は、−(ＰＲＦ／４)と＋(ＰＲＦ／４)により境界づけられ、中心軸付近に中心がある。静止組織から戻る組織高調波成分は、中心軸に見出され、流れている高調波造影剤などの動いている物質から戻る高調波成分は中心軸から外側に、−Ｖhと＋Ｖhと表示された矢印により示されるように、−(ＰＲＦ／４)と＋(ＰＲＦ／４)の境界内に見出される。動いていないおよび動いている組織や血液から戻る基本周波数成分など、線形成分が存在する領域は、−(ＰＲＦ／４)および＋(ＰＲＦ／４)境界から外側に位置する。静止物質からの基本成分は、＋Ｖfおよび−Ｖfと表示された矢印に示されたように、両末端の内側に位置する速度が増加する物質からの信号成分を有し、スペクトルの−(ＰＲＦ／２)および＋(ＰＲＦ／２)末端に位置付けられるであろう。従って、適切な領域からの信号成分を選択することによって、使用者は静止目標からの線形信号、静止目標からの非線形信号、動いている目標からの線形信号、または動いている目標からの非線形信号など、所望の分割された信号成分の画像処理または分析を実施することができる。
【００２７】
本発明の実施と結果の数学的表現を以下に示す。一連のＮ個の超音波パルスがＴ秒のパルス繰返し間隔で発信される。各発信パルスは、前のパルスの反転されたコピーである。
【００２８】
【数２】
0004365909.tif 000003
【００２９】
ここにＰk(t)は発信パルスｋの圧力波形である。受信エコー、ek(t)は、典型的には変換器通過帯域の中心近く（広帯域発信パルス用）、または発信周波数の２次高調波もしくはその近傍（狭帯域発信パルス用）のいずれかに選択される周波数ｆcの搬送波(carrier)の直角位相にそれらを混合することにより復調される。復調された信号は次いで処理されドップラースペクトルを作成する。線形分布構造には、連続エコーは下式で表される。
【００３０】
【数３】
0004365909.tif 000004
【００３１】
ここに、Δτｄ＝２Δｄ／Ｃｓ  であり、
【００３２】
それは、パルス間で変換器に関係する分散構造の軸運動Δｄによる時間遅れである（Ｃsは、媒体中の音速）。ｆｃと関連して測定された、連続エコーの間の対応する位相シフトは下式で表される。
【００３３】
【数４】
0004365909.tif 000005
【００３４】
付加されたπ項は、数式２のekの前の負記号の結果である。計算されたドップラー周波数は、従って下式となる。
【００３５】
【数５】
0004365909.tif 000006
【００３６】
ここにvは変換器に関する散乱物体の相対軸速度である。そしてｆｐｒ＝１／Ｔは 、 パルス繰返周波数（ＰＲＦ）である 。従って、線形散乱物体からの信号は、ナイキスト周波数と等しい量だけ移動して、ドップラースペクトル中に現れるであろう。非線形散乱からのエコーは、下式のように、奇数と偶数成分に分解することができる。
【００３７】
【数６】
0004365909.tif 000007
【００３８】
連続ＰIＤエコーにおいて、eodd(t)成分はＰ(t)と正負符号を交換するが、
eeven(t)成分は交換しない。この分解は、一般的であり、分散過程の正確な性質に依存しない。下式の簡単なモデルにとって、
【００３９】
【数７】
0004365909.tif 000008
【００４０】
eodd(t)は、奇数次の分散項によるエコーを含み、eeven(t)は偶数次の分散項によるエコーを含む。これら２成分によるドップラーシフトは下式となる。
【００４１】
【数８】
0004365909.tif 000009
【００４２】
fd,evenが、従来のドップラーにより作成されるものと同一のドップラーシフトを有することに注意しなければならない。全ての散乱物体速度が、修正(revised)ナイキスト限界に従うとすると、
【００４３】
【数９】
0004365909.tif 000010
【００４４】
それは従来の限界の半分であり、−fpr/４とfpr/４の間のドップラースペクトルの部分は、非線形散乱プロセスから生じる、唯一のドップラー信号を含む（これをここに、スペクトルの低周波数または非線形部分という。）。ドップラースペクトルの残り半分（高周波数または線形部分と呼ぶ。）は、線形散乱と、おそらくより高次の非線形散乱から生じるドップラー信号を含む。
【００４５】
図６に示すように本発明による発信波の位相交替により作成された周波数分離は、相補的である。図６ａは、２ＭＨｚに中心がある名目上(nominal)の発信周波数に応答して受信される２ＭＨｚの名目上の周波数付近の基本周波数の帯域、それに続いてそれぞれ約４、６、８および１０ＭＨｚに中心がある第２、第３、第４および第５次高調波帯域を示す。パルス反転発信と、取得信号の減算の組合せにより実施されるスペクトル分離が、図６ｂに示される。基本周波数を含む奇数次の高調波は、この処理によって分離される。同様に、加算的組合せによるパルス反転発信は、図６ｃに示すように、２次高調波を含む偶数次の高調波を分離する。
【００４６】
従って、図６ｂのスペクトルは、図５ｂで「線形」と標識された領域に存在し、図６ｃのスペクトルは図５ｂに「非線形」と標識された領域に存在する。しかしながら、超音波画像処理の変換器は有限の通過帯域を有することを思い起こされなければならない。ワシントン州、ボゼルのエイティーエル・ウルトラサウンド・インコーポレイテッド(ATL Ultrasound Incorporated)によって製造されたＰ４−２フェーズドアレー変換器により示された通過帯域は、ほぼ図６ｄの通過帯域１３０によって示されたものであり、２ＭＨｚから４ＭＨｚまで実質的に平坦である。従ってＰ４−２フェーズドアレー変換器は、上記通過帯域１３０に示された基本と２次高調波帯域を再現し、より高次の帯域は変換器の上部遮断(upper cutoff)により除去される。Ｐ４−２変換器によって画像化された非線形および線形スペクトルは従って、高次高調波の大きなエネルギーなしで、基本と２次高調波成分から構成される。これは、Ｐ４−２変換器が２ＭＨｚの発信波の基本または２次高調波エコーのいずれか単独の画像処理に使用しうることを意味している。
【００４７】
図７−１０は、本発明による様々な診断用画像処理方法で取得しうる種々のドップラー信号のスペクトルの例を示す。図７ａは、低強度波の発信に伴う実質的に静止している組織からのエコーのスペクトルを示す。実質的に全てのスペクトル成分は線形であり、ドップラースペクトルのナイキスト末端の帯域２０２,２０４に位置している。発信されたエネルギーの強度が増加すると、図７ｂのスペクトルが結果として生ずる。同じ静止線形帯域２０２,２０４が存在する。しかし、発信エネルギーの増加は、発信波中の高調波成分のより速やかな増強を生じさせ、２次高調波エネルギーは、スペクトルの非線形領域の中央（静止）軸付近に位置する帯域２０６により示されるように、静止組織から戻るエコー中に含まれる。
【００４８】
もし図７ｂのスペクトル中の２次高調波組織成分だけの画像形成を望むならば、図１２のフィルタ通過帯域３０６に示すものなどの、単純応答特性をフィルタに付与するように、マトリックスフィルタの係数またはＦＩＲフィルタの重みを調整する。このフィルタ通過帯域は、通過帯域の中心の組織高調波成分を通過させるが、通過帯域３０６の外側の基本成分は通過させないことが分かる。従って、ウォ−ルフィルタはさらなる処理のために使用される信号成分の種類の選択に使用される。
【００４９】
図８ａは、灌流中の心筋層など、灌流されている組織内の高調波造影剤の低発信強度による検出から得られたスペクトルの内容を示す。造影剤は、帯域２１２,２１４に示されたように、発信波の線形成分を、そして非線形帯域２１６によって示されたように剤の非線形特性による非線形信号成分も反射する。これらの帯域の全ては、ゼロ速度中心およびドップラー帯域の両末端付近に位置し、これは造影剤の組織内の遅い流速から予想される。発信エネルギーが増加し、造影剤の微小気泡を破壊させると、これらの結果のスペクトルは相当に広がり、通常図８ｂの帯域２１８に示すように、ドップラースペクトル全体を満たす。帯域２１８の信号成分をパワードップラー記録方式に効果的に表示することができ、またはアメリカ合衆国特許５,４５６, ２５７に記載されたように微分し、振幅検出することができる。
【００５０】
図９ａは、血管内を流れる造影剤のスペクトルを示す。予想されるように、造影剤のエコーの基本と高調波帯域の両方は、帯域２２２, ２２４の位置により示されるように、周波数（速度）Δfだけ静止軸から離れるように動く。ここに、
【００５１】
Δｆ＝２ｖｆｃ／ｃである（式９参照）。
【００５２】
流れている造影剤が、より高い発信エネルギーに曝されると、微小気泡が破裂し、図９ｂに示すスペクトルを返す。このスペクトルは、微小気泡の破裂により放出された広帯域エネルギーのために広がり、流れる剤の進行速度の基本および２次高調波の戻りに対応する周波数２２６,２２８の頂点を有する。
【００５３】
図１０ａは、いかなる造影剤も不存在時に、中程度の発信エネルギーに応答して静止組織から戻るエネルギーのスペクトルを示す。エコーは、帯域２３２,２３４に示される基本周波数エネルギーと中程度量の帯域２３６の組織の高調波エネルギーを戻す。発信エネルギーが実質的に増加すると、図１０ｂの帯域２３６’に示したように振幅が大きくなるにつれて、発信波の高調波エネルギーは大きい速度で蓄積する。
【００５４】
図７−１０のスペクトルと帯域は、発信波のサンプリング速度（ＰＲＦ）が、遭遇する流れまたは運動速度のナイキスト基準を満足していると仮定している。もし高速度域がサンプリング不足になると、生じるスペクトルの高調波は、サンプリング速度基準が満足される領域で部分的に分離されるのみで、サンプリング不足のより高速度域ではスペクトルの他の高調波領域と重なるであろう。
【００５５】
図１１−１５は、ウォ−ルフィルタ用の係数および重みの選択により生じうる多数の通過帯域を示す。図７−１０と同様に、図１１−１５のスペクトル図の横座標は、標識されていないがＰＲＦ／４分画で示されている。より複雑な機能は通常、より高次のフィルタとより大きい集合長さを必要とする。図１１は、図８のスペクトルである、心筋層内の高調波造影剤の検出に有用なウォ−ルフィルタ通過帯域を示す。低発信エネルギーで、図１１ａの通過帯域３１６は、図８ａの帯域２１６からの戻りを極大化するように、非線形領域での中心の低速度利得を極大化するのに使用される。微小気泡の破裂が始まったときに、より穏やかなエネルギーで、低速非線形領域により大きい利得を与えるために、図１１ｂの通過帯域３１４を使用することができる。より高エネルギーで、図１１ｃの通過帯域３１８は（スペクトルの末端で）線形成分と、（スペクトルの中心で）静止組織からの高調波成分を拒絶するために使用することができ、一方造影剤の微小気泡の破裂から戻る広帯域の戻りを極大化する。簡単にいうと、スペクトルのある領域に対して他の領域中の信号利得を減少、増加させるために、ドップラースペクトルの異なる領域中のフィルタ利得を利用することができる。
【００５６】
図１３は、本発明による、高調波造影剤の色流れ画像処理用の詳細設定された通過帯域３２０を示す。組織高調波からの成分は、スペクトルの中央部では拒絶され、全通過帯域利得がＰＲＦ／４境界までの全高調波周波数（速度）範囲について適用される。
【００５７】
図１４と１５は、 ２次高調波（および、他の偶数次の非線形）および基本（および奇数次の非線形）信号用の理想化された通過帯域を示す。通過帯域３３０は、ＰＲＦ／４でのその鋭い遮断により、全ての出所、高調波造影剤、組織高調波からの２次高調波非線形信号を通過させる。図１５の相補的通過帯域３４０は、全ての原因からの線形（基本）信号成分を通過させる。
【００５８】
１つの画像または分析処理用のスペクトルの２つの領域を分離するのに、複数のフィルタを使用するいくつかの処理を実施することも予想される。例えばそれぞれが、図９の奇数次高調波（線形）と偶数次高調波帯域を別々に通過させ、そして別々に各組の信号をドップラー処理するために、２つのフィルタを使用することができる。これは、基本信号情報の組織ドップラー処理により、偶数次高調波情報から、拍動している心臓の、動いている心壁内を流れている造影剤の表示を可能にする。第二の例として、（図１４に示されたような応答特性を使用して）非線形信号強度を得るために１つのフィルタを使用し、（図１５に示された応答特性を使用して）線形信号強度を得るのに他のフィルタを使用し、その両方をパワードップラー処理法により処理することができる。次いで画像平面内の各画素位置での線形信号パワーに対する非線形信号パワーの割合が計算され、表示される。好適例において、この線形と非線形の比は、結果として生じる係数のパワー依存をより小さくするよう、分母の線形値を２乗することにより計算することができる。拡張された係数の組が、別々に蓄積され、次いでドップラー処理される２つの通過帯域について、結果を算出するので、マトリックスウォ−ルフィルタはそのような「２経路」実行によく役立つ。ドップラースペクトルの線形領域を処理するのに使用される係数は、乗法プロセスの間、集合中の１つおきのエコーを反転させる極性を有することができ、奇数次チャネル中のほぼ±（ＰＲＦ）/２からゼロまでに分布する線形成分を再変調する。
【００５９】
本発明者らは、パワードップラーがほとんどの用途に好適ドップラー技術であり、流れの方向でなく、流れの存在が必要な多くの場合に優れた画像を作成することを見出した。潅流中の心臓組織の、潅流している心筋のパワードップラー画像が優れた解像度で作成されている。本発明による色ドップラー画像処理は、例えば、腎臓移植の画像化に非常に有用である。グレースケール表示も可能である。ドップラー集合は通常、２から１６の試料の長さに渡りうる一方、より短い集合長さで優れた結果を達成することができることが見出され、２試料集合を有するパワードップラー画像は、高解像度画像を作成し、高い実時間フレーム速度をも与えることが見出された。より長い発信パルスが、ある状況では望ましい。本発明者らは、図１のパワーモ−ションプロセッサ５０を使用するパワーモーション画像処理技術により、パルス反転技術により動いている組織から得られた２つの試料を効果的に画像化できることも見出した。パワーモーション画像処理技術とプロセッサ５０の詳細が、アメリカ合衆国特許５,７１８, ２２９に見出される。
【００６０】
受信された集合のエコーからエコーへの分散を導入する他の変調技術を、許容可能な効果をもって使用することができる。例えば、 ２つの異なったレベルの間で、連続発信パルスの振幅を交互にし、パルス毎の１８０゜位相シフトをもって、またはシフトなしで、そして受信エコーが発信パルス振幅に従って正規化され(normalize)、次いで上記したように処理される。振幅変調は、本発明技術の実行に十分な程度の線形と非線形成分の分離を生じさせるのに十分な変位(variation)を生じさせる。もしも１８０゜位相シフトが使用されないならば、スペクトルの奇数次部分と偶数次部分は逆になる。これが本発明の次善の実行であると判断される一方、それはある環境においては有効であるか、望ましい。
【００６１】
本発明は、ドップラースペクトルが図形的にまたは数値的に表示されるＰＷドップラーの応用に容易に適用することができる。従って、所望のウォ−ルフィルタの選択により、線形または非線形信号成分のドップラースペクトルを、選択的に処理し、表示することができる。
【００６２】
線形散乱から非線形散乱を明瞭に分離するので、パルス反転ドップラーは線形組織中の微小気泡造影剤の検出によく適合している。さらに本願技術は、広帯域、高解像度パルスをドップラー用途に使用することを可能にする。
【００６３】
本発明は、パルス反転ドップラー超音波診断画像処理に関し、広帯域線形と非線形エコー信号成分を分離しつつ、同時に運動の効果を同時に識別する超音波診断画像処理方法と装置が提供される。本発明方法は、エコーの直線性に依存するようにドップラー領域中のドップラーシフト周波数に、無線周波数（時間）領域中の超音波エコー信号をマッピングする。本発明方法は、これを連続するエコー間の位相シフトを分析することにより実施する。
【００６４】
周期的に位相変化する波形の発信シーケンスによりドップラーエコー集合が取得される装置が開示される。エコー集合がドップラー処理中に組合わされると、結果として生ずるドップラースペクトルは、奇数次高調波（例えば線形信号成分）が存在する１つの領域と、偶数次高調波（例えば、２次高調波信号成分）が存在する他の領域というように、別個の領域に分離される。好適例において、ドップラーウォ−ルフィルタは、組織運動、高調波造影剤、組織高調波画像処理情報など、特定種類の処理用信号情報の描出に使用される。
【００６５】
本発明の主たる態様を以下に示す。
１、発生する超音波エコーに位相差を生じさせるよう、体内に変調された超音波シーケンスを発信する段階、該発信シーケンスに応答した超音波エコー信号の集合を受信する段階、そして、線形と非線形信号成分の位相シフト情報を分離するために該集合を分析する段階、以上の段階からなる線形と非線形超音波エコー信号成分の微分的位相変調方法。２、該発信する段階が、パルス毎に変調された超音波の発信からなる前記１の方法。３、該パルス毎に変調された超音波の発信が、線形と非線形信号成分の間に位相の分離を生じさせる前記２の方法。４、該変調が、該発信パルス速度のナイキスト限界（最大測定可能速度）で生じる前記２の方法。５、該変調が、パルス毎の１８０゜位相変化、パルス毎の振幅変調、またはパルス毎の位相と振幅変調の組合せである前記２の方法。６、該分析する段階が、線形と非線形の信号成分を微分的に検出する段階からなる前記１の方法。７、微分的に検出する段階が、該集合をパワードップラー処理することからなる前記６の方法。８、微分的に検出する段階が、該集合を色流れ処理することからなる前記６の方法。
【００６６】
９、位相変化超音波のシーケンスを発信する段階、該位相変化超音波に応答した超音波エコー信号の集合を受信する段階、そして、エコー信号の該集合をドップラー処理し、エコー信号の直線性に従って分離されたドップラー位相シフト情報を作成する段階、以上の段階からなる、線形と非線形信号情報を分離するようにドップラー領域に、無線周波数超音波エコー信号情報をマッピングする方法。１０、一つの波から次へと異なって振幅変調された超音波シーケンスを体内に発信する段階、該振幅変調超音波に応答した超音波エコー信号の集合を受信する段階、そして、エコーの該集合の位相シフト特性を基準に、受信超音波エコー信号を処理し、分離された線形と非線形信号成分を作成する段階、以上の段階からなる、線形と非線形の超音波信号成分を分離する方法。１１、被検者の体内に超音波の変調シーケンスを発信する超音波発信器、超音波の該変調発信シーケンスに応答した超音波エコー信号の集合を受信する受信器、および、該エコーの位相シフト情報を線形と非線形信号成分に分離するプロセッサ、からなる、線形と非線形超音波エコー信号成分を識別する超音波診断画像処理装置。１２、該発信器が、パルス毎に変調された超音波を発信する１１の超音波診断画像処理装置。１３、該発信パルス毎の変調が、線形と非線形信号成分の間に位相分離を導入することからなる前記１２の超音波診断画像処理装置。１４、該変調が、該発信パルス速度(rate)のナイキスト限界で生じる前記１２の超音波診断画像処理装置。１５、該変調が、パルス毎の１８０゜位相変化、パルス毎の振幅変調、またはパルス毎の位相と振幅変調の組合せからなる前記１２の超音波診断画像処理装置。１６、該プロセッサが、線形と非線形信号成分を微分的に検出する微分プロセッサからなる前記１１の超音波診断画像処理装置。１７、該プロセッサが、パワードップラープロセッサからなる前記１６の超音波診断画像処理装置。１８、該プロセッサが、色流れプロセッサからなる前記１６の超音波診断画像処理装置。
【００６７】
１９、被検者の体内に位相変化超音波のシーケンスを発信する超音波発信器、該位相変化超音波に応答した超音波エコー信号の集合を受信する受信器、および、エコー信号の該集合を処理し、該超音波エコー信号の直線性に従って描出されたドップラー位相シフト情報を作成するドップラープロセッサ、からなる線形と非線形の超音波信号情報を分離する超音波診断画像処理装置。２０、被検者の体内に位相変化超音波のシーケンスを発信する発信器、該超音波に応答するしたエコーの集合を受信する受信器、および該集合を処理し、動いている組織の奇数次または偶数次高調波信号の効果のスペクトルを識別する組織ドップラープロセッサ、からなる動いている組織の高調波信号効果を識別する超音波診断画像処理装置。２１、被検者の体内に位相変化超音波のシーケンスを発信する発信器、該超音波に応答したエコーの集合を受信する受信器、および該集合を処理し、動いている組織または高調波造影剤の奇数次または偶数次高調波信号効果のスペクトルを識別するパワードップラープロセッサ、からなる動いている組織または高調波造影剤の高調波信号効果を識別する超音波診断画像処理装置。２２、被検者の体内に位相変化超音波のシーケンスを発信する超音波発信器、該位相変化超音波に応答した超音波エコー信号の集合を受信する受信器、およびエコー信号の該集合を処理し、ドップラー位相シフト情報を作成するドップラープロセッサ、および、奇数次高調波ドップラー位相シフト情報と偶数次高調波ドップラー位相シフト情報を、それぞれ通過させる第１と第２フィルタ、からなる奇数次と偶数次高調波の超音波信号情報を分離する超音波診断画像処理装置。２３、さらに、奇数次と偶数次高調波位相シフト情報の比を計算するための該第１と第２のフィルタに接続する比回路からなる前記２２の超音波診断画像処理装置。
【００６８】
【発明の効果】
運動の効果を識別しつつ、同時に広帯域の線形および非線形エコー信号成分を分離する超音波診断画像処理方法および装置が提供される。該方法は、超音波エコーの直線性により、無線周波数領域の超音波エコー信号を、ドップラー領域のドップラーシフト周波数にマッピングする。変調発信シーケンスに応答して受信された連続エコー間の位相シフトを分析することによりこれを実施する。周期的に位相変化する波形のシーケンスの１つの発信により１つのドップラーエコー集合が取得される。エコーの集合が組合わされたドップラースペクトルは、奇数次高調波（例えば線形信号成分）が存在する１つの領域と、偶数次高調波（例えば第二高調波信号成分）が存在するその他の領域という、別々の領域に分けられる。各領域は、運動していない体内の静止物質からの線形と非線形エコー中に含まれる画像化される物質の運動に対応する周波数範囲を含んでいる。好適例においては、組織の運動、高調波造影剤、組織高調波画像処理情報などの、特定種類の処理用信号情報を同時に描出するのにドップラーフィルタが使用される。
【図面の簡単な説明】
【図１】  本発明超音波診断画像処理装置のブロックダイヤグラムである。
【図２】  エコー集合を取得するための好適発信波形シーケンスを示す。
【図３】  図１の集合記憶装置の説明図である。
【図４】  図１のウォ−ルフィルタの２つの相互の例を示す。
【図５】  従来の及び本発明のドップラースペクトルを示す。
【図６】  本発明による異なる高調波スペクトルの分離を示す。
【図７】  分析、画像処理される超音波情報スペクトルの第１の例である。
【図８】  分析、画像処理される超音波情報スペクトルの第２の例である。
【図９】  分析、画像処理される超音波情報スペクトルの第３の例である。
【図１０】  分析、画像処理される超音波情報スペクトルの第４の例である。
【図１１】  超音波情報の描出におけるウォ−ルフィルタ特性の第１の例である。
【図１２】  超音波情報の描出におけるウォ−ルフィルタ特性の第２の例である。
【図１３】  超音波情報の描出におけるウォ−ルフィルタ特性の第３の例である。
【図１４】  超音波情報の描出におけるウォ−ルフィルタ特性の第４の例である。
【図１５】  超音波情報の描出におけるウォ−ルフィルタ特性の第５の例である。
【符号の説明】
１０・・・プローブ、１２・・・アレー変換器、１４・・・ビーム形成器制御器、１６・・・発信ビーム形成器、１７・・・発信／受信スイッチ、１８・・・受信ビーム形成器、２０・・・ＱＢＰフィルタ、２２・・・周波数混合回路、２４・・・検出器、３０・・・ウォールフィルタ、３２・・・エコーデータバス、５０・・・パワーモーションプロセッサ、６０・・・グレースケールプロセッサ、７０・・・３Ｄプロセッサ、８０・・・走査変換器、８２・・・画像データバス、９０・・・表示装置、１２０・・・マルチプライヤ、１２４・・・アキュムレータ。
【図１】
0004365909.tif 000011
【図２】
0004365909.tif 000012
【図３】
0004365909.tif 000013
【図４】
0004365909.tif 000014
【図５】
0004365909.tif 000015
【図６】
0004365909.tif 000016
【図７】
0004365909.tif 000017
【図８】
0004365909.tif 000018
【図９】
0004365909.tif 000019
【図１０】
0004365909.tif 000020
【図１１】
0004365909.tif 000021
【図１２】
0004365909.tif 000022
【図１３】
0004365909.tif 000023
【図１４】
0004365909.tif 000024
【図１５】
0004365909.tif 000025