JPH07503552A - 周波数領域障壁フィルタを用いたカラー・フロー作像システム - Google Patents
周波数領域障壁フィルタを用いたカラー・フロー作像システムInfo
- Publication number
- JPH07503552A JPH07503552A JP6516228A JP51622894A JPH07503552A JP H07503552 A JPH07503552 A JP H07503552A JP 6516228 A JP6516228 A JP 6516228A JP 51622894 A JP51622894 A JP 51622894A JP H07503552 A JPH07503552 A JP H07503552A
- Authority
- JP
- Japan
- Prior art keywords
- signal
- echo signal
- frequency
- output
- frequency domain
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details 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/52053—Display arrangements
- G01S7/52057—Cathode ray tube displays
- G01S7/52071—Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8979—Combined Doppler and pulse-echo imaging systems
- G01S15/8981—Discriminating between fixed and moving objects or between objects moving at different speeds, e.g. wall clutter filter
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるため要約のデータは記録されません。
Description
【発明の詳細な説明】
周波数領域障壁フィルタを用いたカラー・フロー作像システム
関連出願の記載
本出願は、1992年に出願され、本出願人に譲渡されたダフト等による出願(
社内控え番号RD−21580)に関連する。
発明の背景
本発明は、音響作像に関し、特に、超音波ドツプラ信号からのカラー・フロー画
像の発生に関する。
超音波のような振動エネルギを用いて物体の像を作成することのできる多数のモ
ードがある。超音波送信器を物体の一方の側に配置し、物体を透過した音波を反
対側に配置されている超音波受信器で受け取る(「送信モード」)。
この送信モードの方法では、画素の各々の輝度が受信器に達する超音波の振幅の
関数である(「減衰」モード)か、又は画素の各々の輝度が音波が受信器に到達
するのに要した時間の関数である(「飛行時間」又は「音速」モード)ような像
を発生することができる。この代わりに、送信器と同じ物体の側に受信器を配置
し、画素の各々の輝度が物体から受信器に反射された超音波の振幅又は飛行時間
の関数である(「屈折」、「後方散乱」又は「エコー」モード)ような像を発生
することができる。本発明は、超音波像を作成する後方散乱方法に関連する。
超音波データを収集する多数の後方散乱方法がよく知られている。所謂「八−走
査」方法では、送信器によって超音波パルスを物体に送り込み、ある期間にわた
って反射された音波の振幅を記録する。エコー信号の振幅は、物体内にある散乱
体の散乱強度に比例し、遅延時間は、送信器からの散乱体の距離に比例する。所
謂「B−走査」方法では、送信器が1つの運動軸線に沿って物体を走査するとき
に一連の超音波パルスを送信する。この結果得られるエコー信号が八−走査方法
の場合のように記録され、エコー信号の振幅を用いて、遅延時間における表示装
置の画素の輝度を変調する。B−走査方法では、走査された視野の画像を再生す
ることができるようにするのに十分なデータが収集される。
所謂C−走査方法では、変換器を物体の上方の平面にわたって走査し、変換器の
焦点深度から反射されるエコーのみを記録する。CRT表示装置の電子ビームの
掃引を変換器の走査と同期させて、変換器のX及びy座標が像のX及びy座標と
対応するようにする。
医療用の超音波変換器は、一対の電極の間に挟まれている1つ又は更に多くの圧
電素子から構成されている。このような圧電素子は典型的にはチタン酸ジルコン
酸鉛(PZT)、二弗化ポリビニリデン(PVDF)又はPZTセラミック/ポ
リマ複合体で構成されている。電極は電圧源に接続されており、電圧波形が印加
されたときに、圧電素子は印加電圧の周波数に対応する周波数で、寸法が変化す
る。
電圧波形が印加されたとき、圧電素子は、圧電素子が結合されている媒質内に超
音波を放出する。逆に、超音波が圧電素子に入射すると、素子はその電極の間に
対応する電圧を発生する。典型的には、素子の前側は音響整合層で覆われており
、音響整合層は、超音波がその中を伝搬する媒質との結合を改善する。更に、圧
電素子の背側には支持材料が接続されており、支持材料は素子の裏側から出て来
る超音波を吸収し、超音波が干渉しないようにする。このような超音波変換器の
多数の構造が、いずれも本出願人に譲渡された米国特許番号第4217684号
、同第4425525号、同第4441503号、同第4470305号及び同
第4569231号に記載されている。
超音波作像に用いるとき、典型的には、変換器は列に配置されている多数の圧電
素子を有しており、これら多数の圧電素子は別々の電圧で駆動される(アポダイ
ゼーション)。印加電圧の遅延時間(又は位相)及び振幅を制御することにより
、圧電素子によって発生される超音波が組み合わさって、選択されたビーム方向
に沿って進むと共にそのビームに沿った選択された点に焦点を合わされた正味の
超音波を発生する。印加電圧の遅延時間及び振幅を制御することにより、ビーム
をその焦点と共に平面内で移動させて、被検体を走査することができる。
変換器を用いて、反射された音波を受信するときにも(受信モード)、同じ原理
が働く。即ち、配列内の変換器素子に発生される電圧は、正味の信号が被検体の
1つの焦点から反射された音波を表すように合計される。送信モードの場合と同
じく、超音波エネルギをこのように焦点を合わせて受信することは、配列の各々
の変換器素子からの信号に別々の遅延時間(及び/又は位相シフト)並びに利得
を与えることにより達成される。
このような形式の超音波作像は、[位相調整アレイ扇形走査]又はPASS”
と呼ばれる。この走査は一連の測定で構成されており、この測定で、進路を決め
られた超音波を送信し、短い期間の後にシステムが受信モードに切り換わり、反
射された超音波を受信して記憶する。典型的には、音響ビーム又は走査線に沿っ
た一連の点からのデータを収集するために、毎回の測定の間、送信及び受信は同
じ方向(θ)に方向法めされる。反射された超音波を受信するとき、受信器は走
査線に沿った一連の距離(R)の所に動的に焦点合わせされる。走査全体を行う
のに必要な時間は、毎回の測定を行うのに要する時間と、所望の分解能及び信号
対雑音比で関心のある領域全体をカバーするのに必要な測定の回数との関数であ
る。例えば、各々の走査線を0.70°の増分で進め、合計128本の走査線を
90゜の扇形にわたって取得することができる。このような多数の超音波作像シ
ステムが、本出願人に譲渡された米国特許番号第4155258号、同第415
5260号、同第4154113号、同第4155259号、同第418079
0号、同第4470303号、同第4662223号、同第4669314号及
び同第4809184号に記載されている。
ドツプラ・シフトを利用して心臓及び血管の血液流を測定することは、周知であ
る。反射波の振幅は組織のグレー画像を発生するために用いられているが、反射
組織の速度を測定するために、反射波の周波数シフトを用いてもよい。
カラー・フロー画像は、流れる組織のカラー画像を黒白の大きさ画像に重ねるこ
とにより発生される。各画素で測定された流れの速度によって、そのカラー(色
)が決定される。
反射された超音波のドツプラ測定をする際の困難さの主な点は、エコー信号には
典型的には、静止した又はゆっくりと移動する組織によって発生された多くの成
分が含まれていることにある。静止した組織は反射波におけるいかなる周波数シ
フトも発生せず、これらの成分は、流れの測定に影響を及ぼすことなく、容易に
フィルタで取り除くことができる。しかしながら、心臓及び血管の運動する壁に
よって発生された反射波は、周波数シフトされ、ゆっくりと流れる血液と区別し
にくい。従来のシステムには、「障壁(つオール・)フィルタ」が設けられてお
り、このフィルタは、搬送周波数の中心に位置しているエコー信号の狭帯周波数
を取り除くように、オペレータによって手動で調節されている。オペレータのす
るべきことは、反射された障壁信号を、血液流の測定をゆがめることなく除去す
るような方法で、このフィルタの幅を調節することである。このような静的なフ
ィルタは、画像の視野内の一部の位置では十分に作用し得るが、他の位置に対し
ては適切な設定ではないことがある。
発明の要約
本発明は、超音波作像システム用のカラー・フロー・プロセッサに関し、特に、
その中心周波数及び帯域幅を、受は取ったエコー信号の関数として自動的に調節
する適応障壁フィルタに関する。更に具体的に言うと、カラー・フロー・プロセ
ッサは、第1の自動相関評価装置を含んでいる。
第1の自動相関評価装置は、復調されたベースバンド・エコー信号を受け取ると
共に、このベースバンド・エコー信号の平均周波数及びこの平均周波数からの変
動を示す出力信号を発生する。フーリエ変換回路がベースバンド・エコー信号を
受け取ると共に、そのベースバンド・エコー信号を周波数領域エコー信号に変換
する。削除制御手段が、平均周波数信号及び変動信号を受け取ると共にフィルタ
にかけられるべき周波数のバンドを示す削除信号を発生するために用いられてい
る。フィルタ手段が、周波数領域エコー信号及び削除信号を受け取ると共に周波
数領域エコー信号から削除信号によって示される周波数を除去するために設けら
れている。重心手段か、フィルタにかけられたエコー信号を受け取ると共にその
平均周波数を示す出力信号を発生するように接続されている。重心手段の出力信
号に応答する手段が、流れの速度を示すカラー信号を発生するために用いられて
いる。
本発明の全般的な目的は、超音波ドツプラ信号に基づくカラー・フロー測定の精
度を改良することにある。フィルタにかけられていないベースバンド・エコー信
号の測定された平均周波数及び変動によって、血管の壁のような静止した又はゆ
っくりと運動する組織により発生されるドツプラ・シフトしたエコー信号成分の
中心周波数及びバンド幅の正確な指示が得られる。これらの測定された信号は、
戻って来た周波数シフトされたベースバンド・エコー信号の周波数領域表示から
削除された周波数を動的に制御するために用いられる。このフィルタにかけられ
たエコー信号における残りの周波数成分は主として、流れる血液によるものであ
るので、その平均され、ドツプラ・シフトされた周波数のより正確な測定が可能
になる。
本発明の他の目的は、動作条件の広範な変化に合致するように動的に変化する障
壁フィルタを提供することである。
エコー信号の障壁成分が典型的には、平均周波数及び変動量の測定を左右するの
で、これらの測定は、実時間で行われ得ると共に、障壁信号成分の周波数又はバ
ンド幅に拘らず、障壁フィルタの特性を動的に制御するために用いられ得る。
本発明の更に他の目的は、焦点合わせされ復調されたエコー信号に対する動的障
壁フィルタを提供することにある。
エコー信号が復調されて、単一の受信ビームに動的に焦点合わせされた後に受信
器の中間プロセッサ部において本発明を実施することにより、マルチチャンネル
PASSシステムの各チャンネルにおいて障壁フィルタを繰り返す必要がない。
本発明の上述及びその他の目的及び利点は、以下説明するところから明らかにな
ろう。以下の説明では、本発明の実施例を示す図面を参照する。しかしながら、
この実施例は必ずしも本発明の範囲全体を表すものではなく、本発明の要旨を解
釈するに当たっては、請求の範囲によるべきである。
図面の簡単な説明
第1図は本発明を用いた超音波作像システムのブロック図である。
第2図は第1図のシステムの一部を形成している送信器のブロック図である。
第2A図及び第2B図は第2図の送信器5oの任意のチャンネルの信号を表示す
るグラフである。
第3図は第1図のシステムの一部を形成している受信器のブロック図である。
第4図は第1図のシステムの一部を形成している表示システムのブロック図であ
る。
第5図は第3図の受信器の一部を形成している受信チャンネルのブロック図であ
る。
第5A図〜第5E図は第5図の受信チャンネルの種々の点に現れる信号を示すグ
ラフである。
第6図は第3図の受信器の一部を形成している本発明によるカラー・フロー・プ
ロセッサの概略ブロック図である。
第7図は第6図のカラー・フロー・プロセッサの動作説明を補助するためのグラ
フである。
好ましい実施例の説明
第1図について具体的に説明すると、振動エネルギ作像システムは変換器配列1
1を含んでおり、変換器配列11は別々に駆動される複数の素子12で構成され
ている。各々の素子は、送信器13によって発生されたパルス波形によりエネル
ギを与えられたときに、超音波エネルギのバーストを発生する。被検体から変換
器配列11に反射された超音波エネルギは、変換器素子12の各々によって電気
信号に変換され、−組の送信/受信(T/R)スイッチ15を介して受信器14
に別々に印加される。送信器13、受信器14及びスイッチ15は、オペレータ
の指令に応答して、ディジタル制御装置16の制御の下に作動される。一連のエ
コーを収集することにより完全な走査が行われる。
このとき、スイッチ15を送信位置に設定し、送信器13を一時的にオンにゲー
ト駆動して、変換器素子12の各々にエネルギを与え、その後スイッチ15を受
信位置に設定して、変換器素子12の各々によって発生された後続のエコー信号
を受信器14に印加する。変換器素子12の各々からのこれらの別々のエコー信
号を受信器14で組み合わせて、1つのエコー信号を発生し、そのエコー信号を
用いて、表示システム17の機内の線を発生する。
送信器13は、発生される振動エネルギ、例えば超音波エネルギがビームとして
差し向けられ又は方向決めされるように、変換器配列11を駆動する。従って、
変換器配列11を物理的に動かす代わりに、このビームを一組の角度にわたって
点から点へ移動させることによりB−走査を行うことかできる。このためには、
送信器13が相次ぐ変換器素子12に印加されるそれぞれのパルス波形20に遅
延時間(T、)を加える。遅延時間がO(T、=0)であれば、すべての変換器
素子12は同時にエネルギを与えられ、その結果生ずる超音波ビームは、変換器
の面に対して法線方向で、変換器配列11の中心から発する軸21に沿った向き
になる。遅延時間(T、)が増加するにつれて、第1図に示すように、超音波ビ
ームは中心軸21から角度θ分だけ下向きに向けられる。変換器配列11の一端
(i=1)から他端(i−n)までのi番目の信号の各々に加えられる遅延時間
T は、次の式で表される。
ここで、Xは変換器配列の中心から変換器素子12の中心までの距離、θは送信
ビーム角度、Cは被検体内での音速、及びR1は送信ビームが焦点合わせされる
距離である。
式(1)における遅延時間T は、ビームを所望の角度θに方向決めすると共に
ビームを固定された距離R1の所で焦点合わせさせる効果を有している。相次ぐ
励振で遅延時間T を漸進的に変化させることにより、扇形走査が行われる。こ
のため、角度θが増分に分けて変えられて、送信ビームを一連の方向に方向決め
する。ビームの方向が中心軸21より上方であるときに、パルス20のタイミン
グは反対になるが、それでも式(1)は依然として成立する。
更に第1図の説明を続けると、超音波エネルギのパーストの各々によって発生さ
れるエコー信号は、超音波ビームに沿って相次ぐ位置(R)の所にある反射物体
から来る。
これらのエコー信号は変換器配列11のセグメント12の各々によって別々に感
知され、特定の時点におけるエコー信号の大きさのサンプルが、特定の距離(R
)で起こった反射の大きさを表す。しかしながら、反射点Pと変換器素子12の
各々との間の伝搬経路の違いのため、これらのエコー信号は同時には発生せず、
又、それらの振幅は等しくない。受信器14の動作は、これらの別々のエコー信
号を増幅して復調し、各々に正しい遅延時間及び位相変化(シフト)を加え、こ
れらを加算して、角度θに向いている超音波ビームに沿って距離Rの所にある点
Pから反射された全超音波エネルギを正確に表す1つのエコー信号を発生するこ
とである。
各々の変換器素子からのエコーによって発生される電気信号を同時に加算するた
め、受信器14の各々別個の変換器素子のチャンネルに遅延時間及び位相シフト
が導入される。受信時のビーム遅延時間は、前に述べた送信時の遅延と同じ遅延
(T、)である。しかしながら、動的に焦点合わせするために、各受信チャンネ
ルの遅延時間及び位相シフトは、エコーを受信する間、連続的に変化して、エコ
ー信号が出てきた距離Rに受信ビームを動的に焦点合わせする。各変換器素子に
よって受け取られた信号に課される遅延時間の正確な式は、次の通りである。
ここで、tは変換器配列の中心からの音の伝播後の経過時間、Cは被検体内の音
速、θはビーム角度、及びXは変換器配列の中心から素子の中心までの距離であ
る。
適当な測定による同様な計算によって、正確な位相シフトがめられる。
ディジタル制御装置16の指示の下に、受信器14は、受信器14の方向決めが
送信器13によって方向決めされたビームの方向を追跡するように、走査の間遅
延を発生すると共に、一連の距離(R)の所でエコー信号をサンプリング(標本
化)して、ビームに沿った点Pに動的に焦点合わせするのに適当な遅延及び位相
シフトを発生する。このため、超音波パルス波形が発せられる度に、その超音波
ビームに沿った所にある対応する一連の点Pからの反射された音波の大きさを表
す一連のデータ点が収集される。
表示システム17は、受信器14によって発生された一連のデータ点を受け取り
、このデータを所望の像を発生する形式に変換する。例えば、A−走査を所望す
る場合には、一連のデータ点の大きさは単に時間の関数として図示される。B−
走査を所望する場合には、この一連の中の各々のデータ点は像の画素の輝度を制
御するために用いられ、相次ぐ方向決め角(θ)における一連の走査で構成され
た走査が実施されて、表示に必要なデータを発生する。
第1図と共に第2図を参照して説明すると、送信器13は、包括的にメモリ50
として示す一組のチャンネル・パルス・コード・メモリを含んでいる。好ましい
実施例では、128個の別個の変換器素子12があり、従って、128個の別個
のチャンネル−パルス・コード−メモリ50がある。パルス・コード・メモリ5
0の各々は典型的には、1ビツト×512ビツト・メモリであり、発生すべき超
音波パルス52の周波数を決定するビット・パターン51を記憶している。好ま
しい実施例では、このビット・パターンを40MHzのマスククロックによって
、各々のパルス・コード・メモリ50から読み出し、駆動器53に印加する。
この駆動器はこの信号を、変換器11を駆動するのに適した電力レベルまで増幅
する。第2A図に示す例では、ビット・パターンは4つの「1」ビットと4つの
「0」ビットとが交互に並んだ順序であって、5MHzの超音波パルス52を発
生している。しかしながら、他の搬送波周波数(Fo)として、2.5MHz、
3.75MHz、6.25MHz、7.5MHz、8.75MHz及び10MH
zのような周波数を好ましい実施例において用いてもよい。
これらの超音波パルス52が印加される変換器素子12は、それに応答して超音
波エネルギを発生する。512個のビットすべてを用いた場合、(5M Hzの
搬送波周波数では)その搬送波周波数を中心とする40kHzという幅の狭い帯
域幅のパルスが放出される。
前に述べたように、超音波エネルギの送信ビームを所望の方向(θ)に方向決め
するためには、第2B図に示すようなnチャンネルの各々に対するパルス52を
適正量だけ遅延させなければならない。これらの遅延は送信制御装置54によっ
て行われる。この制御装置は、ディジタル制御装置16(第1図)から4つの制
御信号(開始、マスタクロック、R1及びθ)を受け取る。入力制御信号θ、一
定の送信焦点R1及び前掲の式(1)を用いて、送信制御装置54は、相次ぐ送
信チャンネルの間に必要な遅延増分子を計算する。開始制御信号を受け取ったと
き、送信制御装置54は40MHzのマスタクロック信号が取り得る4つの位相
のうちの1つを第1の送信チャンネル50にゲートする。その後、相次ぐ遅延期
間(T、)をおいて、40MHzのマスタクロック信号が、次のチャンネルのパ
ルス・コード・メモリ50を介してゲートされ、これがn=128個のチャンネ
ル全部がその超音波パルス52を発生するまで続けられる。各々の送信チャンネ
ル50は、そのビットパターン51全体が送信された後にリセットされ、その後
、送信器13はディジタル制御装置16からの次の入力制御信号θ及び次の開始
制御信号を待つ。前に述べたように、本発明の好ましい実施例では、完全なり一
走査は変換器11の中心軸21 (第1図)を中心とする90°の扇形にわたっ
て、0.70°のΔθ増分に分けて方向決めされる128個の超音波パルスで構
成されている。
送信器13についての詳しいことは、1991年5月14日に付与された発明の
名称「振動エネルギビームの送信の動的な焦点合わせを行うための符号化された
励振(Coded Excitation Fo+ Transmission
Dynamic Focusing of Vib+ata+y Energ
y Beam ) Jという本出願人に譲渡された米国特許番号第501471
2号を参照されたい。
特に第3図を参照すると、受信器14は3つの部分、即ち、時間−利得制御部1
00と、受信ビーム形成部101と、中心プロセッサ102とで構成されている
。時間−利得制御(TGC)部100は、n=128の受信チャンネルの各々に
対する増幅器105と、時間−利得制御回路106とを含んでいる。各々の増幅
器105の入力はそれぞれ1つの変換器素子12に接続されており、増幅器10
5は変換器素子12が受信したエコー信号を受け取って増幅する。増幅器105
による増幅の程度は、時間−利得制御回路106によって駆動される制御ライン
107を介して制御される。周知のように、エコー信号の距離が増加するにつれ
て、エコー信号の振幅は減少する。その結果、より遠い反射体から出るエコー信
号が近くの反射体からのエコー信号よりも一層大きく増幅されなければ、像の輝
度は距離(R)の関数として急速に低下する。この増幅はオペレータによって制
御される。オペレータは手作業で(典型的には)8つのTGC線形ポテンショメ
ータ108を、扇形走査の範囲全体にわたって比較的一様な輝度が得られるよう
な値に設定する。エコー信号が収集される期間が、エコー信号が出てきた距離を
決定し、この期間をTGC回路106によって8つのセグメントに分割する。8
つのポテンショメータの設定により、8つのそれぞれの期間の各々の間、増幅器
105の利得を設定して、エコー信号がエコー信号収集期間にわたって次第に増
加する分だけ増幅されるようにする。
受信器14の受信ビーム形成部101は、n=128個の別個の受信チャンネル
110を含んでいる。後で更に詳しく説明するが、各々の受信チャンネル110
は、その人力111に1つのTGC増幅器105からのアナログ・エコー信号を
受け取り、■母線112及びQ母線113にディジタル化された出力値のストリ
ームを発生する。これらの■及びQの値の各々が、特定の距離(R)におけるエ
コー信号の包絡線の復調されたサンプルを表す。これらのサンプルは遅延及び位
相シフトさせられて、サンプルを加算点114及び115で他の各々の受信チャ
ンネル110からのI及びQサンプルと加算したとき、加算されたサンプルが方
向決めされたビーム(θ)上の距離Rの所にある点Pから反射されたエコー信号
の大きさ及び位相を示すようにする。好ましい実施例では、各々のエコー信号は
、走査線の範囲全体(典型的には40ミリから200ミリ)にわたって約150
ミクロンの増分の所でサンプリングされる。
受信器14についての更に詳しいことは、1991年1月8日に付与された発明
の名称「ディジタル位相配列作像の方法及び装置(Method And Ap
puajus lot Digital t’hase A++ay Imag
ing) Jという米国特許番号第4983970号を参照されたい。
第3図の説明を続けると、中心プロセッサ部102が加算点114及び115か
らの受信ビーム・サンプルを受け取る。ビーム・サンプルの各々の1及びQの値
は、点(R。
θ)から反射された音波の大きさの同相及び直角位相成分を表す32ビツトのデ
ィジタル数である。中心プロセッサ102は、これらのビーム・サンプルに対し
て種々の計算を行うことができ、どれを選択するかは、再生する像の種類によっ
て決定される。例えば、普通の大きさの像を発生する場合には、検出プロセッサ
120で検波過程が行われ、各々の受信ビーム信号サンプルからディジタルの大
きさCM)を次の式に従って計算し、反射点の座標R1θと共に出力121に供
給する。
検出プロセッサ120は又、1989年5月30日に付与された発明の名称「位
相結合を用いた適応可能なコヒーレント・エネルギ・ビームの形成(Adapt
ive Coherent Enetgy Beam Formation U
sing Phase Conjuguion ) Jという本出願人に譲渡さ
れた米国特許番号第4835689号に開示されているような補正方法をも実施
し得る。このような補正方法は、受は取ったビーム・サンプルを検査し、ヒーム
ノ焦点合わせ及び方向決めを向上させるように、送信器13及び受信器14によ
るその後の測定において使用され得る補正値を算出する。このような補正は、例
えば各変換器素子からの音波が走査の間に進む媒体の非均質性を説明するために
必要である。
本発明は、中心プロセッサ102に配置されているカラー・フロー・プロセッサ
123に特に関する。カラー・フロー・プロセッサ123は、第6図に関して後
で詳しく説明するように、加算点114及び115からのそれぞれ復調及び焦点
合わせされた信号サンプルの値I及びQを受け取っており、中心プロセッサの出
力121に12ビツトの流れ値を発生する。この流れ信号は、同じ点に対してデ
ィジタルの大きさMで表示システム17に印加され、表示システム17において
赤、緑及び青の表示カラーを制御する。
後で詳しく説明するように、この流れ値によって示されるカラー(色)は、カラ
ー・フロー・プロセッサ123によって測定される流れの速度及び方向の関数で
ある。
特に第1図及び第4図について説明すると、受信器14はその出力121にディ
ジタル数のストリームを発生し、このストリームは表示システム17の入力に印
加される。
これらの「走査データ」が配列としてメモリ150に記憶され、走査データ配列
150の行が、収集されたそれぞれのビーム角度(θ)に対応し、走査データ配
列150の列が、各々のビームに沿ってサンプルが収集されたそれぞれの距離(
R)に対応する。受信器14からのR及びθ制御信号151及び152は、配列
150内のどこに各々の入力の値を記憶すべきかを示し、メモリ制御回路153
がその値を配列150内の正しいメモリ位置に書き込む。この走査を連続的に繰
り返すことができ、受信器14からの値の流れが走査データ配列150を連続的
に更新する。
更に第4図について説明すると、配列150内にある走査データは、ディジタル
走査変換器154によって読み取られ、所望の像を発生する形式に変換される。
例えば、通常のB−走査像を発生するとき、走査データ配列150に記憶されて
いる大きさの値M (R,θ)は、機内の画素位置(x、y)における大きさを
示す大きさの値M (x、y)に変換される。このような超音波像データの極座
標からデカルト座標への変換は、例えばヒユーレットパラカード・ジャーナル誌
、1983年10月号、第30頁〜第33頁所載のステイーブンC,リービット
等の論文「超音波像を表示するための走査変換アルゴリズム(A 5can C
onve+si。
n A1go+i+bm lot Displaying Ult+a+oun
d Images) J lこ記載されている。
ディジタル走査変換器154によって行われる特定の変換が何であっても、こう
して得られた像データはメモリ155に書き込まれる。このメモリは、変換され
た走査データの2次元配列を記憶する。メモリ制御回路156がメモリ155に
アクセスする二重ポートを有しており、ディジタル走査変換器154がメモリ1
55内の値を新しいデータで連続的に更新することができるようにすると共に、
表示プロセッサ157が更新されたデータを読み取るようにする。表示プロセッ
サ157は制御パネル158から受け取ったオペレータの指令に応答して、メモ
リ155内の変換された走査データに対する通常の像処理作用を行う。例えば、
メモリ155内の変換された走査データが示す輝度レベルの範囲は、表示装置1
60の輝度範囲をずっと超えることがある。実際、メモリ155内の変換された
走査データの輝度の分解能は、人間の目の輝度の分解能を遥かに超えることがあ
り、典型的には、オペレータが、最大の像コントラストを達成しようとする輝度
の値のウィンドウを選択することができるようにする手動で操作し得る制御装置
が設けられている。表示プロセッサはメモリ155から、変換された走査データ
を読み込んで、所望の像の強化を行い、強化された値を表示メモリ161に書き
込む。
表示メモリ161はメモリ制御回路163を介して、表示制御回路162と共有
されており、表示メモリ161内の値が表示装置160の対応する画素の輝度及
びカラーを制御するために写像される。表示制御回路162は、用いられる特定
の型の表示装置160を作動するように設計された市場で入手し得る集積回路で
ある。例えば、表示装置160はCRT (陰画線管)であってもよく、この場
合、表示制御回路162は、水平及び垂直掃引回路に対して必要な同期パルスを
供給すると共に掃引中の適当な時刻に表示データをCRTに写像するCRT制御
チップである。
流れ信号は、適当な座標変換の後にメモリ155に記憶される。これらのデータ
は、表示装置160に適当なカラー出力信号を発生するために、表示プロセッサ
157によって用いられる。種々の表示装置のうちのいずれの装置を用いてもよ
い。例えば、それぞれが8ビツトである赤、緑及び青のカラーに対する信号から
成る出力信号が通常、このような装置によって用いられる。従って、表示プロセ
ッサ157は、第6図に関して後で詳しく説明するように、メモリ155に記憶
されているカラー・フロー信号と、赤(R)、緑(G)及び青(B)のカラー信
号との間のマツピング(写像)を作成する。表示プロセッサ157は又、エコー
信号から導出される表示カラー・フロー情報又は表示グレー・スケール情報の選
択を判定する。
特定の超音波システムの能力及び融通性に応じて、表示システム17がいろいろ
な形式のどの形式にしてもよいことは当業者に明らかであろう。上に述べた好ま
しい実施例では、プログラムされたマイクロプロセッサを用いて、ディジタル走
査変換器及び表示プロセッサの機能を実行するが、このため、この結果として得
られる表示システムは非常に融通性があって強力である。
前に第3図について述べたように、受信器14のビーム形成部Lotは、変換器
11(第1図)の各々の素子12に1つずつ、−組の受信チャンネル110で構
成されている。特に第5図について説明すると、各々の受信チャンネルはディジ
タル制御装置16(第1図)からの開始指令、40MHzのマスタクロック、及
びビーム角度信号(θ)に応答して、ディジタルビーム形成機能を実施する。そ
の機能の中には、アナログ/ディジタル変換器200でアナログ入力信号111
をサンプリングし、サンプリングされた信号を復調器201で復調し、復調器2
01によって発生された高周波数の和信号を低域フィルタ202のフィルタ作用
にかけ、減数器203でデータ速度を下げ、その結果得られたディジタル・デー
タ・ストリームを遅延FIFO(即ち、先入れ/先出しメモIJ)204及び位
相回転装置205で時間的に遅延させると共に位相調節をすることが含まれてい
る。これらのすべての素子は、ディジタル制御装置16からの指令に応答して、
必要なりロック及び制御信号を発生する受信チャンネル制御装置206によって
制御される。好ましい実施例では、これらのすべての素子は1つの集積回路に収
容されている。
更に第5図について説明すると、アナログ/ディジタル変換器200は、第5A
図の波形210のグラフで示すアナログ入力信号を、受信チャンネル制御装置2
06からの線211上の遅延サンプル・クロック信号の前縁によって決定された
規則的な間隔でサンプリングする。好ましい実施例では、サンプル・クロック信
号は4QMHzのクロック信号であり、20MHzまでの超音波周波数を用いて
も、ナイキストのサンプリング判断基準に違反しない。例えば、5 M Hzの
超音波搬送波周波数を用いるとき、搬送波の1サイクルごとに8回サンプリング
され、アナログ/ディジタル変換器の出力には40MHzの速度で10ビツトの
ディジタル・サンプルが発生される。これらのサンプルは復調器201に供給さ
れ、この復調器は、各々のサンプルを送信された超音波搬送波と同相の基準、及
び送信された超音波搬送波と直角位相の基準の両方と混合する。復調器の基準信
号は、記憶されている正弦及び余弦テーブルから発生される。これらのテーブル
は受信チャンネル制御信号206からの40MHzの基準クロック信号によって
、それぞれのROM (リード・オンリ・メモリ)メモリから読み出される。正
弦値にサンプリングされた入力信号をディジタル乗算して、低域フィルタ202
に供給される復調同相値■信号を発生し、余弦値に同じサンプリングされた入力
信号をディジタル乗算して、別個の低域フィルタ202に対して復調直角位相値
Q出力信号を発生する。低域フィルタ202は、復調器201から供給される差
周波数を通過させるが、それよりも高い和周波数を阻止するように同調された有
限インパルス応答フィルタである。第5B図の波形215で示すように、このた
め、各々の低域フィルタの出力信号は40MHzのディジタル値のストリームで
あり、これがベースバンド・エコー信号の包絡線の■又はQ成分の大きさを示す
。
アナログ/ディジタル変換器、復調器及び低域フィルタ回路についての詳しいこ
とは、1989年6月13日に付与された発明の名称「高速ディジタル位相調整
アレイコヒーレント作像システム用の方法及び装置(Method and A
pparalus Fo「High 5peed Digital Phase
d A+ra7 Cohe+enlImaging SHtem) Jという本
出願人に譲渡された米国特許番号第4839652号を参照されたい。
更に第5図の説明を続けると、エコー信号の復調されたI及びQ成分をサンプリ
ングする速度を減数器203によって下げる。12ビツトのディジタル・サンプ
ルが40MHzの速度で減数器203に供給されるが、この速度はバンド幅の観
点からすると不必要に高く、これはシステム全体にわたって維持するのに困難な
データ速度である。このため、減数器203は8番目ごとのディジタル・サンプ
ルを選択して、データ速度を5MHzの速度に下げる。このことは、受信チャン
ネル制御装置206によって発生され、受信チャンネルの他の素子を動作するの
に用いられるベースバンド・クロックの周波数に対応する。こうして、減数器2
03の!及びQ出力が、第5C図のグラフの破線220で示すエコー信号の包絡
線のディジタル化されたサンプル219である。減数比及びベースバンド・クロ
ック周波数は、8:1及び5MHz以外の値に変えてもよく、例えば、4:1及
び10MHzのベースバンドを用いてもよい。
この後、復調されて減数されたディジタル・サンプルにより表されるエコー信号
の包絡線は、遅延FIFO204によって遅延され、位相回転装置205によっ
て位相がシフトされて、所望のビームの方向決め及びビームの焦点合わせを行う
。遅延FIFO204は、相次ぐディジタル・サンプルの値が5MHzの速度で
減数器203によって発生されたときに書き込まれるメモリ装置である。このよ
うに記憶された値は、メモリの相次ぐアドレスに書き込まれ、その後メモリ装置
から読み出されて、位相回転装置205に供給される。初期の遅延の大きさは、
そこからディジタル・サンプルを現在供給しているメモリ位置と、現在受は取っ
たディジタル・サンプルが記憶されているメモリ位置との間の差によって決定さ
れる。5MHzのベースバンド・クロック信号が、記憶されているディジタル・
サンプルの間の200ナノ秒の間隔を定め、従って、FIFO204は、200
ナノ秒の増分で測って、最大25.6マイクロ秒までの遅延時間を発生すること
ができる。
遅延FIFO204によって発生される遅延時間は、第5D図にグラフで示され
るように、減数器203によってサンプリングされたデータ点を前進させること
により、エコー信号を受け取っている開動的に変化し得る。サンプリングされた
各データの前進により、遅延FIFO204に供給されているデータ・ストリー
ムは、更に25ナノ秒(1/ 40 M Hz )だけ遅延させられる。受信チ
ャンネル制御装置206からの信号による制御は、サンプリングされたデータの
前進が起こるべきときに行われる。
位相回転装置205はエコー信号のディジタル表示を、位相回転させることがで
きるようにする。位相回転装置205に供給されたI及びQディジタル・サンプ
ルは、第5E図に示すようにフェーサ221によって表すことができ、位相回転
装置205によって発生された回転I及びQディジタル−サンプルは、フェーサ
222によって表すことができる。フェーサの大きさく即ち、それぞれのI及び
Q成分のベクトル和)は変化していないが、■及びQの値は互いに変化しており
、出カフェーザ222は入カフェーザ221から量Δφだけ回転している。受信
チャンネル制御装置206から母線223を介して受け取る位相制御信号に応答
して、位相は(+Δφ)進めることも(−Δφ)遅らせることもできる。位相回
転装置205についての詳しいことは、1990年1月23日に付与された発明
の名称「コルディック複素数乗算器(Co+dic Complex Mull
iplie+)Jという本出願人に譲渡された米国特許番号第4896287号
を参照されたい。
受信チャンネル110の一般的な説明、並びに各受信チャンネル110の■及び
Q出力信号がビーム信号を形成するためにどのように加算されているかについて
の詳しい説明については、1991年1月8日に付与され本出願人に譲渡された
前述の米国特許番号第4983970号を参照されたい。受信チャンネル制御装
置206の詳しい説明については、1992年4月13日に出願された発明の名
称「改良された動的焦点合わせを用いた超音波作像システム(L:ll+aso
und Imaging system With Improved Dyn
amic F。
(uting) Jという本出願人に譲渡された米国特許出願番号第86759
7号を参照されたい。
第6図について説明すると、第3図のカラー・フロー・プロセッサ123は、市
販され入手可能な集積回路から構成されていてもよい。プロセッサ123はバッ
ファ・メモリ300を含んでいる。バッファ・メモリ300は、変換器配列11
が繰り返し点弧される間に受信器のビーム形成部によって発生された復調エコー
信号の■及びQサンプルを記憶していると共に分類しており、これら■及びQサ
ンプルを自動相関評価装置301の1及びQ入力に供給している。メモリ300
は、「コーナ・ターニング・メモリ」として当業者に知られているものである。
自動相関評価装置301は、アイ・イー・イー・イー・トランザクション・オン
・ソエックス・アンド・ウルトラソニックス(IEEET「ansaction
s on 5onics and Ull+a+onic+ ) 、ボリューム
5U−32、第3巻1985年5月に掲載されたカサイ等による論文「自動相関
法を用いた血液流の実時間2次元作像(Real−Time Two−Dime
nsional Blood Flow Imaging Using Alo
co++elalion Technique) Jに記載されているような、
例えば複合結合倍率器(マルチプライア)、一対の遅延器及び一対の積分器から
構成されていてもよい。この回路は、その入力に供給されるべきn個のサンプル
の周波数の平均値を表す出力信号φ(T)を発生する。反射体の移動がなければ
、エコー・サンプルの周波数におけるドツプラ変化は存在せず、この出力信号は
ゼロである。ある方向への移動があれば、出力信号φ(T)は正の値を有し、移
動が反対方向へのものであれば、出力信号φ(T)は負の値を有する。
自動相関評価装置回路301は、更に2つの出力R(0)及びR(T)を有して
おり、これらの出力から、フィルタのかけられていないベースバンド−エコー信
号の変動が計算され得る。この計算は、乱流計算器回路302によって次の式に
従って実行される。
σ2=(1−IR(T)l/R(0))/T2 (4)ここで、Tは超音波パル
スの放出間隔である。この変動σ2は、障壁信号成分の周波数幅を示すものであ
り、自動相関評価装置回路301から平均速度測定値φ(T)と共に、削除制御
LUT (ルック・アップ・テーブル)303のアドレス端子に供給される。
特に第7A図を参照すると、典型的なベースバンド・エコー信号の周波数スペク
トルが示されており、同図において、相対的に大きな障壁信号成分305が、(
ベースバンドに復調された後にゼロである)搬送周波数f。から一方の方向に周
波数シフトされている。又、非常に小さな流れ成分306が、他方の方向に周波
数シフトされている。全信号の平均周波数測定値φ(T)は本質的には、望まし
くない障壁成分305の中心周波数を示すものであり、変動測定値σ2は本質的
には、その幅を示すものである。
再び第6図を参照すると、本発明の適応可能な障壁フィルタは、復調されたベー
スバンド・エコー信号(即ち、搬送周波数f。=0における信号)が時間領域か
ら周波数領域に変換されることを要する。このことは、高速フーリエ変換(FF
T)回路310によって達成される。この回路は、バッファ・メモリ300から
復調されたエコー信号サンプルの■成分及びQ成分を受け取る。FFT回路31
0は16個から64個のサンプルを受け取ったときに、これらのサンプルの複素
フーリエ変換を実行して、ベースバンド・エコー信号の周波数成分を表す16個
から64個の1組のI値及びQ値を発生する。このことが、第7A図に例示した
ベースバンド・エコー信号に対して、第7B図にグラフで示されている。同図に
おいて、曲線311の下方の垂直な各線は、FFT回路310によって発生され
た特定の周波数成分の大きさくマグニチュード)を表している。
破線で表された曲線311は、これらの大きさの包絡線を表しており、第7A図
に示すベースバンド・エコー信号305に対応している。
再び第6図を参照すると、FFT回路310の各周波数成分出力信号、即ちrb
inJは、マグニチュード(大きさ)計算器313に印加される。この計算器は
、I値及びQ値の実成分及び虚成分の二乗の和の平方根を決定する。
(繰り返し点弧の数に関係する)これら16個から64個の周波数成分の大きさ
は、1組のフィルタ・ゲート314を介して順次クロックされて、重心計算器3
15及び出力計算器316に供給される。フィルタ・ゲート314は削除制御L
UT303の出力信号によって、測定された平均周波数φ(T)におけるその周
波数成分、及び測定された変動値σ2によって決定された成分の周囲のバンドが
阻止されるか又はフィルタにかけられるように制御される。このことが第7C図
にグラフで示されている。同図において、障壁信号305の大きさを含んでいる
周波数成分は、フィルタ・ゲート3140作用によって除去されている。残って
いるこれらの成分は、流れによりドツプラ・シフトされたエコー信号成分の正確
な指示を与える。
フィルタ・ゲート314によって除去された周波数成分は、障壁信号成分の平均
周波数及び変動値の関数として変化するであろう。その結果、本発明による障壁
フィルタは、異なる動作条件に自動的に適合する。例えば、障壁信号成分が比較
的狭い周波数のバンドを占める場合には、変動測定値σ2は比較的小さくなるで
あろうが、比較的小さな周波数のバンドがベースバンド・エコー信号から削除さ
れることになる。
重心計算器315は、フィルタにかけられ変換されたベースバンド・エコー信号
の平均周波数を決定する。この計算は、次の式に従って出カスベクトルの第1の
瞬間(モーメント)を計算することにより実行され得る。
f mean=Σf:Power+/ΣPowe ri (5)ここで、iはF
FT信号の各binに対応する値、f+はi番目の各binの周波数、及びPo
we riはi番目の各binにおける信号の二乗である。
重心計算器315によって発生された平均周波数は、流れる反射体の平均速度V
に比例しており、カラー・フロー・プロセッサ123の出力信号の一部を形成し
ている。重心計算器315、乱流計算器318及びマグニチュード計算器313
の出力信号は、カラー・70−申プロセツサ123の出力121に発生されたカ
ラー・フロー信号を含んでいる。これらの信号は、表示システム17に供給され
た入力信号を含んでいる。表示システム17のディジタル走査変換器154(第
4図)の内部で、これらの信号は座標変換されて、上に述べたように極座標から
デカルト座標へデータが変換される。
座標変換の後、重心計算器315からの出力信号Vは、表示システム17の一部
を形成している色度差(クロミナンス)制御ルック・アップ・テーブル・メモリ
317の最少有効数のアドレス入力に印加される。このメモリに設けられている
アドレス可能な線の各々は、上に述べたように、各画素に対して表示制御器16
2(第4図)に読み出される24ビツトの赤、緑及び青のデータを記憶している
。これらのビットのうちの8ビツトが赤の強度を制御し、他の8ビツトが緑の強
度を制御し、残りの8ビツトが青の強度を制御する。これらのビット・パターン
は、信号Vによって表されているような流れ速度が方向又は大きさにおいて変化
するにつれて、位置R1θにおける画素のカラーが変化するように予め選択され
ている。例えば、変換器に向かう流れは赤として示されていてもよ(、変換器か
ら離れる流れは青として示されていてもよい。流れが速ければ速いほど、カラー
は明るくなる。
流れる反射体の性質を示すために、乱流計算器318が出力計算器316、フィ
ルタ314及び重心計算器315の出力に接続されており、ベースバンド・エコ
ー信号の流れ信号成分の周波数の拡がりを示す出力信号を発生する。
乱流が多くの速度を合成したものであるのに対して層流は非常に狭い範囲の速度
を有しているので、この値は流れの乱れを示すものである。乱流値は、カラー・
フロー・プロセッサの出力121に供給されて、色度差制御LUT317の最大
有効数のアドレス端子に印加され、このLUT内の選択された線に対するアドレ
スを完了する。例えば、速度入力信号Vによって色度差制御LUT317にアド
レス可能な赤及び青の各々の値に対して、異なる量の緑を有している1組の16
個の入力(エントリ)がLUT317に設けられている。緑の量は乱流測定値に
よって選択されている。その結果、画像内の一点における平均流れ速度は赤又は
青の量によって示され、乱流の程度は緑の量によって示される。
ベースバンド・エコー信号における流れ信号成分が所定のレベル以下に低下する
と、流れをバックグラウンド・ノイズから正確に識別することが不可能になる。
こうなると、色度差制御LUT317からのカラー信号がディセーブルされて、
カラーは発生されなくなり、画素はエコー信号の大きさのみによって決定される
グレーのレベル又は程度になる。このことは、出力計算器316の出力に接続さ
れていると共に色度差制御LUT317に論理信号を供給する出力しきい値検出
器320によって達成される。流れ信号出力が予め調節された(プリセット)レ
ベルを超えると、色度差制御LUT317は出力しきい値検出器320によって
イネーブルされ、流れ信号出力がプリセット・レベル以下に低下すると、LUT
317はディセーブルされて、表示制御器162(第4図)はエコー信号の大き
さのみによって制御される。
出力計算器316及び乱流計算器318は、フィルタ・ゲート314及び重心計
算器315によって発生されたデータのストリームに次の計算を施す個別のディ
ジタル信号処理部から成っている。出力計算器316に関して、出力される全出
力は次のように表され得る。
ΣPower。
乱流計算器318に関して、
乱れσ2=Σ(fl −f、、、++) 2Powe rl/ΣPower+
=Σf12POwer1/ΣP o w e r lf men。2 (6)
ここで、flはフィルタ・ゲート314の出力binの周波数、及びf41.。
は重心計算器315の出力信号である。
各チャンネル内にベースバンド信号を発生する、ベースバンド・ビーム形成受信
器を用いた超音波作像装置に関連して本発明を説明したが、他の形式の振動エネ
ルギ受信器を用いてもよい。例えば、各信号が適当な量だけ遅延しているチャン
ネル信号の総和によって、所与の方向にビームが形成される無線周波数(r、f
、)システムを用いてもよい。このようにして形成されたr、f、ビームは、中
心プロセッサ部102内でベースバンドに変換される。この代わりに第2の例と
して、受は取ったチャンネルデータを、受信器部101でビーム形成するために
中間周波数に復調してもよい。この例では、ビームも中心プロセッサ部102(
第3図)内でベースバンドに変換される。どちらの場合にも、ビーム・データは
、それらがメモリ300(第6図)に記憶される前にベースバンドに変換される
。
本発明に関連して用いることのできるカラー・フロー作像システムの更に他の形
式は、当業者には明らかであろう。
電気的に走査される配列に関して本発明を説明したが、焦点合わせされた超音波
変換器プローブを機械的に位置決めすることにより、ビーム方向を設定してもよ
い。この場合には、ビームを形成するために、更にビーム形成部を必要としない
。復調は中心プロセッサ102(第3図)で行われ、ベースバンド信号は中心プ
ロセッサ102で、それらがメモリ300(第6図)に記憶される前に形成され
る。
ビーム方向信号151及び152が、プローブの機械的位置を示すインジケータ
から導出される。
更に、信号のディジタル処理に関して本発明を説明したが、アナログ処理、又は
アナログ処理及びディジタル処理の組み合わせも又、可能である。
本発明のある好ましい特徴のみを図面に示して説明したが、当業者には多くの改
変及び変更が考えられよう。従って、請求の範囲は、本発明の要旨の範囲内に属
するこのようなすへての改変及び変更を包括するものであることを承知されたい
。
フロントページの続き
(51) Int、 C1,6識別記号 庁内整理番号GOIS 15158
8113−5J(72)発明者 ベルナルディ、リチャード・ブルースアメリカ
合衆国、19087、ペンシルバニア州、ストラフオード、ウッドフレスト・ロ
ード、440番
(72)発明者 オドネル、マチュー
アメリカ合衆国、48104、ミシガン州、アン・アーバー、プルツクリン・ア
ベニュー、1607番
I
(72)発明者 トーマス、ルイス・ジョーンズ、ザ・サード
アメリカ合衆国、12308、ニューヨーク州、スケネクタデイ、ラグビー・ロ
ード、1142番
(72)発明者 ウェルス、ケニース・ブレイクリイ、ザ・セカンド
アメリカ合衆国、12302、ニューヨーク州、スコティア、ヘラトチエルタウ
ン・ロード、203#
(72)発明者 エンシェラ−、ウィリアム・アーネストアメリカ合衆国、12
302、ニューヨーク州、スコティア、エステイ−・ステイーブンズ・レーン、
88#
Claims (13)
- 1.振動エネルギ変換器から受け取ったエコー信号を復調してベースバンド・エ コー信号を発生する受信器と、カラー表示を発生するカラー・フロー・プロセッ サとを含んでいる振動エネルギ作像システムであって、前記カラー・フロー・プ ロセッサは、 前記ベースバンド・エコー信号を受け取るように接続されており、該ベースバン ド・エコー信号の平均周波数を示す出力信号を発生するように動作可能な自動相 関評価器と、前記ベースバンド・エコー信号を受け取ると共に該ベースバンド・ エコー信号を実成分と虚成分とを有している周波数領域エコー信号に変換するフ ーリエ変換手段と、前記自動相関評価器に応答して、前記周波数領域エコー信号 からフィルタにかけるべき平均周波数を示す削除信号を発生する削除制御手段と 、 フィルタ・ゲート手段と、 前記フーリエ変換手段を前記フィルタ・ゲート手段に接続する手段であって、前 記フィルタ・ゲート手段は、前記削除信号により示される周波数を前記周波数領 域エコー信号から除去するように前記削除制御手段に応答している、接続する手 段と、 前記フィルタ・ゲート手段に接続されており、フィルタにかけられた前記周波数 領域エコー信号を受け取ると共に該フィルタにかけられた周波数領域エコー信号 の平均周波数に対応する出力信号を発生する重心計算器手段と、色度差制御手段 と、 前記重心計算器手段を前記色度差制御手段に接続する手段であって、前記色度差 制御手段は、前記表示のカラーを制御する前記重心計算器手段の平均周波数出力 信号に応答して動作可能である、接続する手段とを備えた振動エネルギ作像シス テム。
- 2.前記ベースバンド・エコー信号の障壁信号成分の周波数幅を示す変動出力信 号を発生する乱流計算器を含んでおり、 前記削除制御手段は、前記乱流計算器から前記変動出力信号を受け取ると共に該 信号に応答して前記周波数領域エコー信号からフィルタにかけるべき前記平均周 波数の周りの周波数のバンドを前記削除信号に示す手段を含んでいる請求項1に 記載のシステム。
- 3.前記フィルタ手段の出力及び前記重心計算器手段の出力を前記色度差制御手 段に結合すると共に、フィルタにかけられた前記周波数傾城エコー信号の周波数 の拡がりに対応する出力信号を前記色度差制御手段に供給するように動作可能な 乱流計算器を含んでいる請求項1に記載のシステム。
- 4.前記フィルタ・ゲートの出力に接続されており、前記フィルタ・ゲートによ り供給された全出力に対応する出力信号を発生する出力計算器手段と、該出力計 算器手段を前記色度差制御手段に接続しており、前記フィルタにかけられた周波 数領域エコー信号の出力レベルがプリセット・レベルを超えたときに前記色度差 制御手段をイネーブルすると共に、前記出力レベルが前記プリセット・レベル以 下に低下したときに前記色度差制御手段をディセーブルする出力しきい値検出手 段とを含んでいる請求項3に記載のシステム。
- 5.前記ベースバンド・エコー信号の障壁信号成分の周波数幅を示す変動出力信 号を発生する第2の乱流計算器を含んでおり、 前記削除制御手段は、前記乱流計算器から前記変動出力信号を受け取ると共に該 信号に応答して前記周波数領域エコー信号からフィルタにかけるべき前記平均周 波数の周りの周波数のバンドを前記削除信号に示す手段を含んでいる請求項3に 記載のシステム。
- 6.前記ベースバンド・エコー信号の障壁信号成分の周波数幅を示す変動出力信 号を発生する第2の乱流計算器を含んでおり、 前記削除制御手段は、前記乱流計算器から前記変動出力信号を受け取ると共に該 信号に応答して前記周波数領域エコー信号からフィルタにかけるべき前記平均周 波数の周りの周波数のバンドを前記削除信号に示す手段を含んでいる請求項4に 記載のシステム。
- 7.前記フーリエ変換手段を前記フィルタ・ゲート手段に接続する前記手段は、 前記フーリエ変換手段により発生された前記周波数領域エコー信号の実成分及び 虚成分の二乗の和の平方根を計算するマグニチュード計算器手段を含んでいる請 求項2に記載のシステム。
- 8.前記振動エネルギは、超音波エネルギを含んでおり、前記振動エネルギ変換 器は、超音波変換器を含んでいる請求項1に記載のシステム。
- 9.振動エネルギ変換器から受け取ったエコー信号を復調してベースバンド・エ コー信号を発生する受信器を含んでいる振動エネルギ作像システムにおいて、表 示システムを制御するカラー・フロー・プロセッサを動作させる方法であって、 前記ベースバンド・エコー信号の平均周波数を示す出力信号を発生する工程と、 前記ベースバンド・エコー信号を周波数領域エコー信号にフーリエ変換する工程 と、 前記周波数領域信号から、フィルタにかけるべき平均周波数を示す削除信号を発 生する工程と、前記削除信号により示される周波数を前記周波数領域エコー信号 からフィルタ除去する工程と、フィルタにかけられた前記周波数領域エコー信号 の平均周波数に対応する出力信号を発生する工程と、前記フィルタにかけられた 周波数領域エコー信号の平均周波数に従って表示のカラーを制御する工程とを備 えたカラー・フロー・プロセッサを動作させる方法。
- 10.前記ベースバンド・エコー信号の障壁信号成分の周波数幅を示す変動出力 信号を発生する工程と、前記周波数領域エコー信号からフィルタにかけるべき前 記平均周波数の周りの周波数のバンドを前記削除信号に示す工程とを含んでいる 請求項9に記載の方法。
- 11.前記フィルタにかけられた周波数領域エコー信号の周波数の拡がりに従っ て前記表示のカラーを制御する工程を含んでいる請求項9に記載の方法。
- 12.前記フィルタにかけられた周波数領域エコー信号の全出力に対応する信号 を発生する工程を含んでおり、前記フィルタにかけられた周波数領域エコー信号 の平均周波数に従って表示のカラーを制御する工程は、前記フィルタにかけられ た周波数領域エコー信号の出力レベルがプリセット・レベルを超えたときに前記 表示におけるカラーの発生をイネーブルすると共に前記出力レベルが前記プリセ ット・レベル以下に低下したときに前記表示におけるカラーの発生をディセーブ ルする工程を含んでいる請求項11に記載の方法。
- 13.前記振動エネルギは、超音波エネルギを含んでおり、前記振動エネルギ変 換器は、超音波変換器を含んでいる請求項9に記載の方法。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US002,000 | 1993-01-08 | ||
| US08/002,000 US5349525A (en) | 1993-01-08 | 1993-01-08 | Color flow imaging system utilizing a frequency domain wall filter |
| PCT/US1994/000236 WO1994016342A1 (en) | 1993-01-08 | 1994-01-06 | Color flow imaging system utilizing a frequency domain wall filter |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH07503552A true JPH07503552A (ja) | 1995-04-13 |
Family
ID=21698785
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP6516228A Pending JPH07503552A (ja) | 1993-01-08 | 1994-01-06 | 周波数領域障壁フィルタを用いたカラー・フロー作像システム |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US5349525A (ja) |
| EP (1) | EP0630484A1 (ja) |
| JP (1) | JPH07503552A (ja) |
| WO (1) | WO1994016342A1 (ja) |
Families Citing this family (27)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5446800A (en) * | 1994-06-13 | 1995-08-29 | Diasonics Ultrasound, Inc. | Method and apparatus for displaying angiographic data in a topographic format |
| AU4523796A (en) * | 1994-12-30 | 1996-07-24 | Acuson Corporation | Imaging modality showing energy and velocity |
| US6177923B1 (en) | 1994-12-30 | 2001-01-23 | Acuson Corporation | Imaging modality showing energy and velocity |
| US5568446A (en) * | 1995-08-17 | 1996-10-22 | General Electric Company | Dual mode ultrasonic imager system |
| JP3682104B2 (ja) * | 1995-12-27 | 2005-08-10 | 株式会社東芝 | 超音波診断装置 |
| US5823964A (en) * | 1996-12-30 | 1998-10-20 | Siemens Medical Systems, Inc. | Ultrasound doppler wall filter |
| US6001063A (en) * | 1998-06-23 | 1999-12-14 | Acuson Corporation | Ultrasonic imaging method and apparatus for providing doppler energy correction |
| US6146331A (en) * | 1998-09-30 | 2000-11-14 | Siemens Medical Systems, Inc. | Method for improved clutter suppression for ultrasonic color doppler imaging |
| US6210334B1 (en) | 1999-03-31 | 2001-04-03 | Acuson Corporation | Medical diagnostic ultrasound method and apparatus for harmonic detection using doppler processing |
| US6685645B1 (en) | 2001-10-20 | 2004-02-03 | Zonare Medical Systems, Inc. | Broad-beam imaging |
| US6733455B2 (en) * | 1999-08-20 | 2004-05-11 | Zonare Medical Systems, Inc. | System and method for adaptive clutter filtering in ultrasound color flow imaging |
| US6267725B1 (en) * | 1999-09-13 | 2001-07-31 | General Electric Company | Individual channel analog wall filtering to improve flow sensitivity in ultrasound imaging |
| US6287258B1 (en) | 1999-10-06 | 2001-09-11 | Acuson Corporation | Method and apparatus for medical ultrasound flash suppression |
| US7269750B1 (en) * | 2001-06-15 | 2007-09-11 | Silicon Motion, Inc. | Method and apparatus for reducing power consumption in a graphics controller |
| US7297118B2 (en) * | 2003-06-12 | 2007-11-20 | Ge Medical Systems Global Technology Company | Ultrasound method and apparatus for multi-line acquisition |
| US20050054931A1 (en) * | 2003-09-09 | 2005-03-10 | Clark David W. | Tracking clutter filter for spectral & audio doppler |
| USD546837S1 (en) * | 2005-10-13 | 2007-07-17 | Samsung Electronics Co., Ltd. | MP3 player |
| FR2913171A1 (fr) | 2007-02-28 | 2008-08-29 | France Telecom | Procede de determination de la presence d'un signal de telecommunications sur une bande de frequences. |
| US20100004542A1 (en) * | 2008-07-03 | 2010-01-07 | Texas Instruments Incorporated | System and method for ultrasound color doppler imaging |
| JP5653057B2 (ja) * | 2009-05-27 | 2015-01-14 | キヤノン株式会社 | 測定装置 |
| US8795179B2 (en) * | 2011-04-12 | 2014-08-05 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Methods, modules, and systems for gain control in B-mode ultrasonic imaging |
| US9420997B2 (en) | 2012-06-14 | 2016-08-23 | Siemens Medical Solutions Usa, Inc. | Motion artifact suppression in ultrasound diagnostic imaging |
| US9943289B2 (en) | 2013-05-29 | 2018-04-17 | B-K Medical Aps | Color flow ultrasound imaging |
| US10908130B2 (en) * | 2014-10-01 | 2021-02-02 | Sensor Networks, Inc. | Asset integrity monitoring using cellular networks |
| US10247705B2 (en) | 2014-10-01 | 2019-04-02 | Sensor Networks, Inc. | Asset-condition monitoring system |
| EP3187896A1 (en) * | 2015-12-29 | 2017-07-05 | Esaote S.p.A. | Method and system for performing baseband digital receiver beamforming on ultrasound signals |
| US12493067B1 (en) * | 2024-06-10 | 2025-12-09 | Bae Systems Information And Electronic Systems Integration Inc. | Non-linear flitering for pulse detection and pulse width discrimination |
Family Cites Families (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4839652A (en) * | 1987-06-01 | 1989-06-13 | General Electric Company | Method and apparatus for high speed digital phased array coherent imaging system |
| US4835689A (en) * | 1987-09-21 | 1989-05-30 | General Electric Company | Adaptive coherent energy beam formation using phase conjugation |
| US4896287A (en) * | 1988-05-31 | 1990-01-23 | General Electric Company | Cordic complex multiplier |
| US4930513A (en) * | 1988-07-26 | 1990-06-05 | U.S. Philips Corporation | Two dimensional processing of pulsed Doppler signals |
| JPH0614932B2 (ja) * | 1989-05-29 | 1994-03-02 | 株式会社東芝 | 超音波診断装置 |
| US5014712A (en) * | 1989-12-26 | 1991-05-14 | General Electric Company | Coded excitation for transmission dynamic focusing of vibratory energy beam |
| US4983970A (en) * | 1990-03-28 | 1991-01-08 | General Electric Company | Method and apparatus for digital phased array imaging |
| DE4134724C2 (de) * | 1990-10-24 | 1995-11-16 | Hitachi Medical Corp | Einrichtung zur farbigen Strömungsaufzeichnung mit Ultraschall |
| US5086775A (en) * | 1990-11-02 | 1992-02-11 | University Of Rochester | Method and apparatus for using Doppler modulation parameters for estimation of vibration amplitude |
| US5048528A (en) * | 1990-11-15 | 1991-09-17 | North American Philips Corp. | Alias suppression in pulsed doppler systems |
| US5226420A (en) * | 1991-06-07 | 1993-07-13 | Advanced Technology Laboratories, Inc. | Ultrasonic color flow imaging using autoregressive processing |
| US5228009A (en) * | 1992-04-10 | 1993-07-13 | Diasonics, Inc. | Parametric clutter elimination |
-
1993
- 1993-01-08 US US08/002,000 patent/US5349525A/en not_active Expired - Fee Related
-
1994
- 1994-01-06 WO PCT/US1994/000236 patent/WO1994016342A1/en not_active Ceased
- 1994-01-06 EP EP94906548A patent/EP0630484A1/en not_active Withdrawn
- 1994-01-06 JP JP6516228A patent/JPH07503552A/ja active Pending
Also Published As
| Publication number | Publication date |
|---|---|
| WO1994016342A1 (en) | 1994-07-21 |
| EP0630484A1 (en) | 1994-12-28 |
| US5349525A (en) | 1994-09-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5349524A (en) | Color flow imaging system utilizing a time domain adaptive wall filter | |
| US5349525A (en) | Color flow imaging system utilizing a frequency domain wall filter | |
| US5142649A (en) | Ultrasonic imaging system with multiple, dynamically focused transmit beams | |
| US5230340A (en) | Ultrasound imaging system with improved dynamic focusing | |
| US6193665B1 (en) | Doppler angle unfolding in ultrasound color flow and Doppler | |
| US5235982A (en) | Dynamic transmit focusing of a steered ultrasonic beam | |
| EP1041395B1 (en) | Method and apparatus for positioning region of interest in image | |
| US5291892A (en) | Ultrasonic flow imaging | |
| US5121364A (en) | Time frequency control filter for an ultrasonic imaging system | |
| JP2000083958A (ja) | 鏡面状超音波反射体をイメ―ジングするシステム及び方法 | |
| US6267725B1 (en) | Individual channel analog wall filtering to improve flow sensitivity in ultrasound imaging | |
| US5706818A (en) | Ultrasonic diagnosing apparatus | |
| JPH05237117A (ja) | コヒーレント振動エネルギビーム作像装置及び伝達媒質内の収差を補正する方法 | |
| JPH1156837A (ja) | 超音波イメージング・システム及び超音波イメージング・システムの開口をアポダイズする方法 | |
| US5327894A (en) | Wall filter using circular convolution for a color flow imaging system | |
| JPH0347A (ja) | 超音波診断装置 | |
| KR100483783B1 (ko) | 도플러 신호 처리 장치 및 초음파 진단 장치 | |
| EP1016879A2 (en) | Ultrasound color flow imaging utilizing a plurality of algorithms | |
| US6123672A (en) | Color flow imaging for enhancing segmentation and flow dynamics | |
| US6174284B1 (en) | Ultrasound harmonic imaging techniques | |
| KR910009411B1 (ko) | 초음파 촬상장치 | |
| JPH0595948A (ja) | 超音波診断装置 | |
| EP1774360B1 (en) | High resolution images from reflected wave energy | |
| JPH069556B2 (ja) | 血流速度表示装置 | |
| JPH06114058A (ja) | 超音波診断装置及び超音波プローブの走査方法 |