JP2002539881A - パルス圧縮を使用する超音波イメージング - Google Patents
パルス圧縮を使用する超音波イメージングInfo
- Publication number
- JP2002539881A JP2002539881A JP2000607546A JP2000607546A JP2002539881A JP 2002539881 A JP2002539881 A JP 2002539881A JP 2000607546 A JP2000607546 A JP 2000607546A JP 2000607546 A JP2000607546 A JP 2000607546A JP 2002539881 A JP2002539881 A JP 2002539881A
- Authority
- JP
- Japan
- Prior art keywords
- pulse
- function
- phase
- time
- envelope
- 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
- 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/8959—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using coded signals for correlation purposes
-
- 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/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
- G01S7/52038—Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
Abstract
Description
特許出願第09/283,346号(代理人ドケット5050/540)に関連する。
超音波イメージが、医療イメージング(像形成)における貴重な新しい診断モー
ドとして医療用超音波分野に出現してきた。好ましい組織調波イメージは、送信
された基本波信号の第2高調波を取り込んで形成される。非線形コントラスト剤
を注入して使用すると、基本周波数の第2高調波、または基本周波数の半分、即
ち低調波のような、送信された基本波周波数以外の周波数帯からの信号を更に増
加させる。最近の数年は、診断検査を向上させるために、非線形コントラスト剤
の注入使用に関心が集まっている。より近年になって、コントラスト剤の分裂を
最小にし、コントラスト剤と組織との間のコントラストを改善し、そして小血管
を視覚化する目的で、低振幅励振の使用に関心が寄せられてきた。
、そしてコントラストを改善するので、従来は患者をイメージングするのが困難
であった組織調波イメージが好まれてきた。組織によって生成される調波の量は
、実質的に基本波よりも少なく、調波エネルギは変換器の面から徐々に累積され
、そして組織減衰は周波数が高い程大きいから、調波エネルギから生成されるイ
メージは、基本波エネルギだけから形成されるイメージに比して信号対雑音比(
SNR)が劣る。従って、調波イメージは、変換器面付近の浅い領域において、
基本波イメージの浸透限界付近の深い領域において、走査フォーマットの外側エ
ッジにおいて、及び個々の超音波線及び個々の変換器素子のための大きい舵取り
角において、臨界的な診断情報に欠ける可能性がある。送信電圧を増加させると
、戻り信号レベル及び屡々SNRが増加するが、実際の電気的限界によってある
最大値に到達する。
ントラスト剤破壊のようなコントラスト剤分裂が最小であり、より多くのコント
ラスト剤がより長い検査時間のために利用可能になり、そして小血管の流れの検
出を増加させることができる。しかしながら、低振幅励振は、コントラスト剤が
存在しないか、または低濃度で存在する場合には、受入れ難い貧弱なイメージ品
質を発生し得る。更に、コントラスト剤信号は第2高調波周波数帯域で戻り、低
調波及び超高調波周波数のような他の帯域は普通の送信パルスによって制限され
得る。コントラスト剤を分裂させずに送信電圧を増加させることはできないから
、低振幅励振を用いるイメージは貧弱なSNRを呈する。
の平均電力を増加させるための公知の技術である。この技術は、1940年代及び19
50年代にレーダー応用のために開発され、近年になって医療用超音波の分野に示
唆されてきた。基本波イメージングに関しては、M. O'Donnell, Coded Excitati
on System for Improving the Penetration of Real-Time Phased Array Imagin
g Systems, IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequen
cy Control, Vol. 39, No. 3, pp. 341-51, May 1992を参照されたい。またコン
トラスト剤イメージングに関しては、Y. Takeuchi, Coded Excitation for Harm
onic Imaging, Ultrasonics, PH-3, 1996を参照されたい。
向細部分解能に重大な(または、何等かの)損失を生ずることなく、送信パルス
を符号化し、組織調波イメージ形成プロセスにパルス圧縮受信フィルタを適用す
ることによってSNRを向上させる方法を実現している。
ス圧縮受信フィルタリング技術と共に使用する装置は、基本波信号を使用する装
置に比して、SNRの改善及びクラッタ雑音の減少についてより多くの妥協を示
すことができる。詳述すれば、組織調波イメージは本質的に、浸透が小さく、視
野のエッジにおけるSNRが低いが、普通の基本波イメージと比較した場合には
クラッタ雑音アーチファクトが少なく、より多くの診断情報を供給することが多
い。即ち、あるイメージング環境においては、基本波イメージにおけるSNR及
び浸透を増加させても、強いクラッタ雑音アーチファクトのために診断情報を増
加させることはないが、組織調波イメージにおいてはSNRを同じように、また
は等価的に増加させると診断情報を増加させることができる。
たは分数調波エネルギを検出するための独特な非線形位相変調符号化計画を実現
する。SNRの改善及びコントラスト剤の特異性の増加によって、診断情報を増
加させることができる。
ーク時間平均限界、またはそれ以外によって制限されてはいない場合に、基本波
イメージと比較した時の調波イメージのSNRが劣るのを解決するために、送信
される信号は、合計エネルギを増加させるように長くする(圧縮していない、ま
たは伸張させる)ことができる。信号処理鎖の後半における軸方向細部分解能を
最大にするために、どの次数調波信号も再圧縮することができる。この方法は、
総合帯域幅(軸方向細部分解能に比例する)を維持しながら、送信される信号の
時間的な幅を増加させることによって、時間・帯域幅積を効果的に増加させる。
公知の符号化送信技術と、以下に開示する独特な受信フィルタリングとを組合せ
ることによって、どの次数調波イメージをも生成することが可能になり、超音波
イメージの診断値を改善することができる。
に、組織調波イメージングにおけるSNRを改善するために使用することができ
る。1つの好ましい実施の形態では、コントラスト剤を使用せずに、組織の第2
高調波イメージングがSNRを改善し、それが浸透の増加及び診断情報の改善に
変換される。
システムのブロック図である。図1のイメージングシステム10は変換器(トラ
ンスデューサ)12を含んでおり、この実施の形態では、変換器12は一次元フ
ェーズドアレイ変換器である。変換器12は、送受信スイッチ14を介して、送
信ビーム形成器16及び受信ビーム形成器18の両方に結合されている個々の変
換器素子のアレイを含んでいる。送信ビーム形成器16は、入力として符号化送
信機20a、20b、…20Ntから符号化された送信パルスを受け、送信機2
0a、20b、…20Ntからの符号化された送信パルスに適当な遅延及び/ま
たは位相変化を与え、イメージ中の組織の所定の領域内にコヒーレントに追加す
るように時間調整され、位相調整された超音波信号(関連する変換器素子からの
)に変換する。
え、個々の変換器素子から受信したエコー信号を所望のビーム方向に沿ってコヒ
ーレントに合計する。受信ビーム形成器18によって生成され、ビーム形成され
たエコー信号は受信復調器20へ印加され、受信復調器20はエコー信号を所望
周波数範囲に復調する。復調されたエコー信号はデコーディング受信機22へ印
加され、受信機22はビーム形成され、復調されたエコー信号にパルス圧縮関数
を適用する。デコードされたエコー信号は、検出及び表示のためにイメージプロ
セッサ24へ印加される。
、図1のシステムに適用される。超音波線の形成は以下のようにモデル化するこ
とができる。
幅a(t)及び変調周波数fmを有する普通の送信波形Xtx(t)は次のように表すこ
とができる。
の実数部分を表す。一般的に言えば、新しい送信波形は、
、約1の時間・帯域幅積を有するコンパクトな固定帯域幅時間パルスをより大き
い時間・帯域幅積に時間的に広げ、より大きい送信信号エネルギを可能にするこ
とができる。組織からの散乱した信号を受信した後に、適切に設計された受信機
が時間的に広げられたパルスを圧縮し、コンパクトなパルスを回復することがで
きる。このプロセスは、広げられた送信パルスから重大な劣化を生ずることなく
、指定された帯域幅に関連する好ましい軸方向細部分解能を保持する。優れた軸
方向細部分解能は、適切なデコーディング受信機を用い、帯域外周波数成分を十
分に排除することによってのみ維持される。所望の第2、低次または高次調波エ
ネルギを選択的に取り込むように設計された受信機を用いて適切にデコーディン
グすると、時間・帯域幅を約1の普通の値を越えて増加させ、SNRを改善する
ことができる。当分野においては公知のように、SNRを最大にすることができ
る基本波イメージングのための好ましいデコーダは整合(マッチド)フィルタで
ある。従って、優勢な基本周波数において復調した後に、受信機デコーダは、
ニズムに依存する。この場合も、関心調波に整合した整合フィルタが好ましいデ
コーダである。第2高調波組織イメージングの場合、公知のリーマンの波動方程
式を解くことによって合理的なモデルを求めることができる(Nonlinear Acoust ics , Mark F. Hamilton, et al. 1998, pg 75)。
しい境界条件、即ち強制関数
ckstock, copyright 1998の75ページを参照されたい。
参照)を使用すると、一般モデルは
高調波信号が得られる。上式のKtは、組織材料特性に等しい。第2高調波エネル
ギのコントラスト剤イメージングに適用することができるモデルの例は時間導関
数を排除するが、それでも基本波圧力の自乗に依存している。
信伝達関数が関心周波数帯域にわたって重大な影響を及ぼさないものとすれば、
組織において生成される第2高調波信号は次のように示すことができる。
的に周波数変調されており、最後の項
い。この最後の項は、非線形音響伝播信号ピークが、有限帯域幅信号のための信
号トラフよりも速く走行する、即ち正の相対圧力(ピーク)が負の相対圧力(ト
ラフ)よりも大きい組織音速度を発生するから現れるのである。コントラスト剤
の第2高調波イメージングの場合、
延、位相調整、及びアポジゼーションが適用された後に、さらなる処理のために
関心調波信号を所望の周波数帯域に復調することができる。好ましい帯域はベー
スバンドである。もし復調関数d(t)を次式のように第2高調波周波数の中心に
選択し、
方向細部分解能を維持する好ましい実施の形態を実現することができる。なお、
式(13)において、 組織第2高調波エネルギの場合には、
ス圧縮受信機Xrx (t)のために、SNRを最大にすることができる整合フィルタ
を使用している。普通の受信機は、デコーディング、及び所望の軸方向分解能を
回復するために必要な付加的な非線形位相変調が欠如している。復調の後の好ま
しいパルス圧縮整合受信機は、
帯域を好ましい周波数帯域へ効果的にシフトさせる。整合フィルタ、またはより
一般的にはパルス圧縮フィルタは、直流に近い低周波数項、及び基本波エネルギ
または他の望ましくないエネルギを効果的に抑圧することができる。もし望むな
らば、整合フィルタ包絡線以外の振幅変調を受信機内で使用することもできる。
詳述すれば、付加的なパルス成形を使用して、不要範囲のローブを抑圧するのを
援助することができる。
線形周波数変調(FM)チャープパルス圧縮コードの例である。送信信号は、以
下の特性を有している。 1)振幅は時間幅
場合には、瞬時周波数(ラジアン/秒)は、
復調後の受信パルス圧縮整合フィルタは、
ことができる典型的な時間・帯域幅積は、
は、Ispta限界、熱限界、及びチャンネル毎の電力によって決定されることが考
えられる。図2乃至9は、あるシステムを実現する上での異なるステージにおけ
る信号の例である。図2乃至5は、時間・帯域幅積が4であるシステムを表し、
一方図6乃至9は、時間・帯域幅積が8、即ち後者の4つの図面を作成するため
に使用した送信信号時間長が、前者の4つの図面を作成するために使用した送信
信号長の2倍であるシステムを表している。4つの図面の2組のために送信され
た帯域幅は同一である。一般的には、時間・帯域幅積が50より小さい、より好ま
しくは20より小さい、最も好ましくは10より小さいシステムが好ましい。
1MHzである。図2a及び6aは、普通の送信される信号、及び対応する送信
包絡線を示しており、各図においては、時間依存・非線形位相変調は使用せず、
これら2つの普通の送信信号は同一であって、各ケースのために参考として示し
たものである。図2b及び6bはそれぞれ、二次PMまたは線形FM符号化され
た送信信号、及び対応する送信包絡線を示している。4つの図面においては、電
気的限界または所望の空間ピーク音響限界に基づいて、典型的なように、ピーク
信号レベルは同一である。図2c及び6cは、先行図に示すように、普通の場合
の線形送信位相、及び符号化信号の場合の非線形位相を示している。図2d及び
6dは、これらの位相の時間導関数によって決定される送信位相に対応する瞬時
周波数を示している。符号化送信信号位相が、普通の送信信号位相とは明白に異
なっており、非線形であることに注目されたい。
絡線を示しており、組織減衰及び音響回折は排除されている。また、これらの図
には、参考として組織基本波信号包絡線が示されている。第2高調波信号が、改
善された軸方向細部分解能を与えることに注目されたい。図3b及び7bはそれ
ぞれ、符号化送信信号に関連する組織第2高調波信号及びその包絡線を示してい
る。これらの図面にも、参考として組織基本波信号包絡線が示されている。図3
c及び7cは、普通の送信信号を用いた時の第2高調波信号、符号化送信信号を
用いた時の第2高調波信号、及び符号化送信信号を用いた時の基本波信号を示し
ている。図3d及び7dは、これらの組織信号位相に対応する瞬時周波数を示し
ている。各型の組織信号毎の、及び各型の送信信号毎の組織信号位相及び瞬時周
波数が独特であることに注目されたい。最適化されたパルス圧縮(または、デコ
ーディング)受信フィルタを設計する場合、これらの差を斟酌する必要がある。
即ち、位相は、符号化送信信号が変換器から発射された時の組織基本波信号に比
して、組織第2高調波信号の場合には、時間の関数として2倍程度速く変化して
いる。図(4a、4b、4c)及び(8a、8b、8c)は、第2高調波及び基
本波で受信した信号をフィルタするための受信ベースバンドフィルタ成分のパル
ス圧縮、またはデコーディングを示している。2組の各々において、これらの図
はそれぞれ、普通の送信信号を用いた時の第2高調波デコーディングフィルタ、
符号化送信信号を用いた時の第2高調波デコーディングフィルタ、及び符号化送
信信号を用いた時の基本波デコーディングフィルタを示している。3つの図の最
後の組は比較のために示したものであり、基本波信号を用いて組織をイメージン
グする場合の典型である。各図は、各フィルタのレスポンスの実数部分、虚数部
分、及び振幅を示している。第2高調波信号のための受信デコーディングフィル
タは、符号化された送信信号が変換器から発射される場合は、基本波信号のため
の受信デコーディングフィルタとは異なることに注目されたい。これらの差は、
受信フィルタの位相を時間の関数としてプロットしても明らかになる。図4d及
び8dは、3つの異なる型のフィルタのためのこれらの位相を示している。これ
らの位相が、符号化送信信号を使用する時の基本波に比して、第2高調波の場合
には時間の関数として2倍程度速く変化することに注目されたい。より高次の、
または低次の調波のイメージングの場合には、デコーディングフィルタの瞬時位
相は、所望の調波位相変化と整合させることが好ましい。例えば、1/2次の低調
波は、その位相が送信された信号の位相の半分のレートで変化するフィルタを用
いて優先的にデコードされる。図4e及び8eは、優勢第2高調波周波数(4M
Hz)における適切な復調、第2高調波パルス圧縮フィルタを用いたベースバン
ド受信フィルタリング(2つの場合に、それぞれ図4aまたは4b、及び図8a
または8bから)、及びディスプレイ検出の後の最終出力信号を示している。以
下のことに注目されたい。
軸方向分解能は同一である。符号化された送信信号及び適切な受信フィルタリン
グを用いると、軸方向細部分解能が維持される。
は、普通の符号化されてない技術に関連する戻りエネルギよりも大きい。時間的
に広げられた送信波形によって、付加的なエネルギが利用可能にされている。2
つの別個のケース(ケース1は時間・帯域幅積=4、ケース2は時間・帯域幅積
=8)の間の戻されて検出される信号エネルギの差は、予測した通り、時間・帯
域幅積がより大きい8の場合の方が大きい。
ために残留基本波信号が存在する。符号化信号処理技術は、優勢残留エネルギを
普通の非符号化信号処理技術に対して効果的に、時間的にシフトさせる。この残
留エネルギは、信号帯域幅が適切に選択されている臨床設定ではそれ程重大では
ない。周波数スペクトルの重畳のために、帯域幅を大きくする程基本波抑圧が小
さくなる。貧弱な基本波エネルギ抑圧、または受入れ難い残留基本波エネルギは
、異なって送信された包絡線位相に関連して予備検出された受信パルスの加算、
またはフィルタリングによって更に減少させることができる。複数の受信したパ
ルスを単純に加算またはフィルタリングによって組合せると、所望の関心調波を
効果的に強調しながら、望ましくない基本波または他の調波エネルギを打ち消す
ことができる。例えば、同一方向に逆相(即ち、0°及び180°)で送信された
2つのパルスは、使用可能な軸方向細部分解能及び第2高調波信号のSNRを同
時に増加させ、基本波信号を抑圧するためのパルス圧縮受信フィルタの後に加算
することができる。別の例では、異なる初期包絡線を有する2つのパルスを異な
る方向に沿って、または異なる原点から送信し、次いで受信してパルス圧縮フィ
ルタでデコーディングの後に加算することができる。非共線対送信パルスは、共
線送信パルスに比してフレームレートを改善することができ、位相及び振幅変調
符号化及びデコーディングのためにSNRが改善される。図4e及び8eにおい
ては、これらの複数のパルス組合せが残留基本波信号内のエネルギを減少させて
いる。
織信号、及び検出された出力信号を示している。これらの図は、送信される信号
の指定された時間幅には無関係に、組織基本波及び組織第2高調波信号の帯域幅
が−6.28dBに一定に保たれていることを示している。また、普通の送信信号か
らの戻りと比較した時の、符号化された送信信号からの戻りの検出された信号エ
ネルギの差が2つのケースについて示されている。2つの時間・帯域幅積の大き
い方の出力信号は、提示されている2つの符号化のケースの間の同一ピーク送信
信号レベルの大きい方であるが、図5及び9の全てのスペクトルは、各ケース毎
に、符号化パルス圧縮技術の調波組織信号ピークによって独立的に正規化されて
いる。従って、信号出力の最大スペクトル振幅は、常に0デシベルである。この
正規化によって、普通の送信信号の組織調波信号レベルと、調波受信フィルタリ
ングの後の残留基本波信号レベルとの差が強調されている。
り、それらは以下のものを含む。
調包絡線の例を説明したが、多くの異なる型の振幅及び位相変調関数を送信機及
びパルス圧縮受信フィルタ内に使用することができる。図10は、線形、非線形
、離散、及び連続瞬時FM関数対時間の幾つかの例を示している。これらの関数
は、位相変調(PM)関数の時間導関数である。これらの関数は、パルス時間幅
TFM及び所望の周波数スパンWFMの範囲内に定義することが好ましい。線30が
包絡線持続時間の半分に位置していることに注目されたい。これらは例であり、
全てを網羅しているのではない。線形FM(二次PM)関数の例が、a、b、及
びcとラベル付けされている。2つの区分的線形FM関数は、関数f及びcであ
る。2つの非線形FM関数は、d及びeとラベル付けされている。変調関数c及
びeは線30に関して対称であり、他の関数は線30に関して非対称であること
に注目されたい。
振幅スペクトルと共に示す。破線で示すスペクトルA及びBは、送信機における
、またはパルス圧縮フィルタの入力における信号スペクトルを表している。実線
曲線で示す位相関数a及びbは、対称及び非対称関数の例である。関数aは振幅
スペクトルAに関して対称であるか、または振幅Bに関して非対称である。関数
bは、振幅スペクトルA及びBに関して非対称である。
振幅変調関数対時間の幾つかの例を示している。関数a及びbは、対称振幅プロ
ファイルであり、関数c及びdは非対称プロファイルである。これらの型の関数
は、所望のシステム総合スペクトルレスポンスに依存して、送信機及びパルス圧
縮受信機に独立的に適用することができる。非対称関数は、組織減衰、音響回折
、及びシステム及び変換器設計からのインパルスレスポンスによって発生するス
ペクトルの成形を補償することによって、SNRを最大にするのに特に有用であ
る。
計のために好ましいものであることができる。例を以下に示す。
時周波数変調関数のような、線形に、または二次的に変化するもの以外の位相変
調関数対時間が好ましいかも知れない。二次位相変調(または、線形FM)を定
振幅変調と共に使用すると、臨床的には受入れ難い範囲のローブを発生し得るが
、定振幅以外の振幅変調関数は送信機において使用することができない。定振幅
または均一振幅変調は、これらの型の送信機のコストが、Acuson Corporationか
らSequoiaの商品名で販売されているプログラマブル波形発生器(PWG)のよ
うな非均一振幅変調を行うより精緻な送信機と比較して低いために、商用超音波
製造業者の間で人気がある。定振幅・非線形FMは、高レンジローブを抑圧する
のを援助するように、パルス圧縮出力スペクトルを優先的に成形するために使用
することができる。代替として、または更に、パルス圧縮受信機における適切な
非均一対称または非対称振幅変調によって、不要なレンジローブを減少させるこ
とができる。受信機における非均一振幅変調の潜在的な欠陥は、SNRが低下す
ることである。非線形FMは、付加的な受信振幅変調の必要性を排除し、従って
定振幅送信されるパルスに伴うSNRを最大にする。
の型には無関係に非線形FMが使用される別の例は、近フィールド調波イメージ
ングである。少なくとも同一の変換器を使用して超音波エネルギを送受信するシ
ステムの場合には、送信機がターンオフされるまでは受信回路はターンオンでき
ないから、変換器の面付近のイメージングはパルス長の拡張に制限を課すことに
なる。また、もし組織非線形伝播によって単独に調波エネルギが生成されるので
あれば、SNRは貧弱になり得る。もし最短の時間窓で最大エネルギを組織内へ
送信するのであれば、均一に振幅変調された波形が好ましい。受入可能なレンジ
ローブレベルを維持するために、非線形FMを使用することができる。
益であるようにランダム化することができる。送信波形の瞬時周波数を制御する
ことによって、信号が0交差する間の時間を精密に制御して疑似ランダムシーケ
ンスを生成することができる。整合フィルタレスポンスを用いて受信フィルタし
た後には時間的にコンパクト出力が発生し、これが優れた軸方向細部分解能、及
び改善されたSNRを与える。更に、もしレンジローブが受入れ難く、且つ非均
一振幅変調が使用できなければ、フレームレートの損失が予測される第2相補送
信発射を使用して最初の送信発射に関連する極性とは逆の極性の副レンジローブ
を生成させ、出力信号の主時間ローブを強調することができる。逆極性のレンジ
ローブは打ち消し合うから、不要レンジローブを抑圧するのを援助するための表
示処理の前にこれら2つの分離出力を加算することができる。これらの型の相補
コードの例は、1961年10月のIRE Transactions on Information Theory, Volume
IT-7, Number 4に所載のM.J.E. GolayのComplementary Seriesに記述されてい
るGolay codesである。
うな、これらの型のコードのいろいろな形状を調波イメージングへ適用すること
に関して2つの重要な面に注目すべきである。第1は、これらのコードの殆どが
要求する第2高調波パルスの反転バージョンを発生させるために、基本波送信パ
ルスを±90°だけシフトさせなければならないことである。例えば、
状の第2高調波信号を実現するためには、送信される信号は次の近似形状を有し
ているべきである。
るためには、送信される信号の2つの成分を時間Tだけ十分に分離すべきことで
ある。そのようにしなければ、調波信号内に第3の望ましくないクロス乗積項が
現れるようになる。説明の都合で示した調波波形は、相補Golayコードを使用し
た符号化の例である。更に、2つの波形の第1の波形は、サイズ2の最も簡単な
Bakerコードを使用して符号化された送信波形の例である。
Rを最大にするためには、非対称位相関数対周波数、即ち図11の関数bのよう
な関数が好ましい。周波数依存減衰は、低めの周波数成分よりも高めの周波数成
分の方を優先的に大きく減衰させるから、戻されたスペクトルエネルギをスペク
トル幅及び中心周波数の両方で追跡する受信機は、SNRを最大にすることがで
きる。深さ依存スペクトルエネルギを所与のベースバンドフィルタ内に選択的に
位置決めする1つの方法は、復調周波数を深さの関数として変化させることであ
る。もし、復調周波数を変化させながら、不要なイメージアーチファクトを発生
させずにSNRを最大にさせれば、受信機フィルタに非対称位相関数を組み込み
、フィルタインパルスレスポンスを深さと共に変化させることができる。この動
的に変化する受信フィルタが好ましいであろう。
機に伴う不正確なデコーディングまたはイメージアーチファクトを最小にするこ
とができる。
を送信変調周波数のような特定周波数に関して対称に保つことである。
ィルタする前に、深さの関数として受信信号を時間遅延させることである。深さ
が増加するにつれて遅延が追加されて増加し、正確なデコーディングが維持され
る。
維持し、可能なアーチファクトを受入れることである。幾つかの位相関数または
対応するFM関数の場合には、イメージアーチファクトは許容できるものである
。
散したセグメントを有する送信パルスが好ましい。この型の特性を有する例は、
図10に関数fとラベル付けされて示されている。関数fは、2つの独立非線形
二次PM関数(または、図示のように、区分的線形FM関数)を組み入れている
。実効的にサブコードからなるこの型のコードは、2つの密に離間する時間間隔
で組織から情報を抽出するために、またはコントラスト剤調波を優先的に励振す
るために使用することができる。独特な振幅及び位相関数を有する1つの時間セ
グメントは、他の符号化されたセグメントの存在に依存する振動モードを優先的
に励振する、または核を作ることができる。この型のコードは、低調波のような
分数調波エネルギを最適に検出するために使用することができる。勿論、複数の
コードを時間的に重畳するように(これは、実効的に2つのコードの加算である
)設計することができる。2つの異なるセグメントを用いているこの例では、パ
ルス圧縮受信機は一方の、または両方のセグメントに関連するエネルギをデコー
ドすることができる。もし一方のセグメントが特定の振動モードを励振するため
に使用されていれば、別のセグメントをイメージング提示のために使用すること
ができる。これらのコードが有用であることができる特定例は、送信されたパル
ス内に含まれる高めの振幅基本波エネルギから、低振幅送信された低調波エネル
ギの成長の核を作るために使用する場合である。この例は、本明細書と同一の日
付で出願され、本明細書が参照しているコペンディング米国特許出願第09/282,6
03号(代理人ドケット5050/515)に詳細に開示されている。
変換器に限定されない。例えば、本発明は静電変換器を用いて実現することがで
きる。
積を深さ依存とすることができる。良好な近フィールドイメージングを維持する
ためには、最も浅い焦点のために小さめの時間・帯域幅積が好ましい。
ることができる。送信機及び受信機は、ディジタル及び/またはアナログ回路を
使用することができ、また送信波形は必ずしもミクサーを必要とはしない。適当
な送信機の例が、米国特許第5,675,554号“Method and Apparatus for Transmit
Beamformer”に開示されている。例えば、簡単なディジタル的に格納された波
形はミクサーを使用せずに増幅することができ、またはアナログ発振器は振幅を
変化させながら異なる周波数にわたって掃引させることができる。フィルタを実
現し、所望の送信コードを生成させるために適切な複素フィルタと共に、インパ
ルス発生器を使用することができる。受信フィルタリングは、上述したアーキテ
クチャの例に限定されるものではない。例えば、復調ステージはビーム形成が完
了した後に発生させる必要はなく、代わりとして、ビーム形成が完了する前に各
受信チャンネル上に、または受信チャンネルのサブグループ上に発生させること
ができる。更に、受信フィルタリングは、ベースバンドに発生させる必要はない
。複素フィルタリングは、中間周波数において実現することができる。
整数調波または分数調波もデコードすることができる。
圧縮調波Bモードイメージを、別の型のイメージ、カラー、またはBモードと共に
、同時または順次取得及び処理と共に使用することができる。これらの例を以下
に示す。 a.符号化パルス圧縮調波イメージと共に基本波イメージ。 b.符号化パルス圧縮調波イメージと共に基本波符号化パルス圧縮イメージ。 c.符号化パルス圧縮調波イメージと共に基本波または調波符号化パルス圧縮
または非符号化カラーイメージ。
分離したビームのどのような予備検出組合せとも使用することができる。組合せ
内の個々のビームは、同一の振幅、周波数、及び位相を有する遅延プロファイル
を有することも、または異なる振幅、周波数、または位相を有することもできる
。この実施の形態の例は、改善された基本波信号排除/打ち消しを得るために、
従って改善された信号帯域幅及び軸方向細部分解能を得るために、逆極性送信位
相調整を用いてパルス圧縮受信フィルタリングした後に、2つのビームを加算す
ることを含む。関連付けられた送信ビームは、普通のパルス反転イメージング(
例えば、Chapmanの米国特許第5,632,277号)のように空間的に整列させることも
、または本明細書が全文を参照しているコペンディング米国特許出願第09/282,3
96号(代理人ドケット5050/542)に開示されているように空間的に別個であるこ
ともできる。
、コントラスト剤を用いるイメージングとに適用されることを記述している。以
下に代替実施の形態を、コントラスト剤イメージングに関連して説明する。1つ
の代替実施の形態は、コントラスト剤を破壊せずにSNRを改善するために、低
圧力振幅で送信されるパルスの時間幅を増加させる。コントラスト剤の整数また
は分数調波を励振するために、時間的に伸張された励振を使用することができる
。第2の代替実施の形態は、幾らかのコントラスト剤は破壊され得るが、それで
もSNRを改善するために、第1の代替実施の形態の圧力振幅よりも大きいが、
それでもある最大値よりは小さい圧力振幅を有するように、送信されるパルス幅
を増加させる。検出された非線形コントラスト剤信号と組織信号との間の最適イ
メージコントラストは、前記最大値よりも小さい送信圧力において得ることがで
きるから、改善されたSNRはイメージ品質及びコントラストを改善する。
積及びパルス幅が上述したように増加された場合には、分解能を損なわずにSN
Rを増加させる。しかしながら、低送信電圧、従って低音響圧力におけるイメー
ジングのようなコントラスト剤イメージングに特定の付加的長所が存在する。第
1に、低送信電圧を使用すると、コントラスト剤の破壊が減少する。長い持続時
間の(大きい時間・帯域幅積の)低音響圧力送信パルスは、圧縮した高音響圧力
送信パルス程コントラスト剤を破壊しないものと考えられる。第2に、音響圧力
を減少させると、非線形に伝播する調波信号の組織散乱に対して、コントラスト
剤の非線形散乱を増加させることができる。考え得る理由の1つの説明は、より
高い音響圧力においては、入射エネルギが散逸するか、(第2高調波散乱に変換
されるのではなく)低調波に、またはより高次の調波に変換されるので、コント
ラスト剤からの散乱に飽和効果が存在するというものである。正味の結果として
、より高い音響圧力よりも低い音響圧力における方が、コントラスト剤調波信号
と組織調波信号との間のコントラストが大きくなり得る。
ト剤が有利であろう。しかしながら、これらの低めの音響圧力においては、SN
Rがかなり劣化する。本発明が概要開示しているパルス伸張及び圧縮は、軸方向
分解能を犠牲にすることなく、そのSNRの幾らかを回復するために使用するこ
とができる。実際に、ピーク振幅を減少させると、合計パルスエネルギを時間平
均パルス電力に関するFDA限界のような限界以下に留めたまま送信パルスをよ
り長くすることができるので、より大きい時間・帯域幅積は、組織高調波イメー
ジングよりもコントラスト剤イメージングの方を支持できることになる。100よ
りも大きい時間・帯域幅積は、多くの理由から実現不能である。当業界で広く行
われているように普通のダイナミック受信集束を使用した場合、これらの長いパ
ルスを正確にデコードしようとすると、必ず受入れ難い程大きいレンジローブが
発生する。ダイナミック受信集束を用いる長いパルスのデコードに伴う不正確さ
を回避するために、パルス圧縮フィルタを多数回反復し、動的な時間遅延を適用
する前に各反復を各受信ビーム形成チャンネルの背後に配置することができる。
しかしながら、多くの素子及びシステムチャンネルを有する殆どの実際のフェー
ズドアレイ変換器システムの場合、コスト及び複雑さが莫大になる。普通の受信
機は、送信機が活動を停止するまで受信を行うことはできないから、スタンドオ
フパッドを使用しない限り、長いパルスは近フィールドイメージングも禁止する
可能性がある。
あり、この散乱のためのモデルは多くのx2(t)次の項を含み得る。ここにx(t)
は入射送信パルスである。例えば、正確なモデルは、
理的な帯域幅である場合には、コントラスト剤からの第2高調波散乱をx2(t)に
よって表す近似を使用すると、有益な結果を得ることができる。これは散乱した
第2高調波信号を、式8のように
フィルタを調波散乱パルスと整合するように設計することによって、幾分良好な
結果を得ることができる。
するために、非線形位相変調送信パルスを適切に設計された受信パルス圧縮フィ
ルタと共に使用することによって、さらなる有益な結果が得られる。交互送信極
性技術に関しては、特に図2乃至9の説明のアイテム3、及び好ましい実施の形
態の例のアイテム7において既述した。一般的には、送信パルス電圧を低下させ
ると、組織非線形伝播による基本波信号レベルに対して第2高調波信号レベルが
低下し、基本波信号の抑圧を効果的に減少させる。更に、送信及び受信に大きい
時間・帯域幅積を使用すると、送信及び受信における有限フィルタ長のために、
基本波信号抑圧が低下するという結果がもたらされる。何れの場合も、交互極性
送信パルスからの空間的に共線の、または空間的に別個のビームの予備検出組合
せを使用して、組織からの基本波散乱の抑圧を増加させることができる。これは
、コントラスト剤の特異性を改善し、受入れ難い程高い基本波信号レベルからの
重大なアーチファクトを生ずることなく、軸方向細部分解能を増加させる。上述
した非線形位相変調された送信パルスは、
スの数まで交互位相項を用いて変更することができる。
の説明の殆どは、受信機が、受信した調波パルスの整合フィルタリングを遂行す
るような好ましい実施の形態に集中していた。ある場合には、送信パルスの帯域
幅よりも狭い、または広い帯域幅の何れかを使用することが望ましいかも知れな
い。その場合、受信機は、受信機の位相ひずみ(位相対周波数)が受信した調波
パルスの位相ひずみ(位相対周波数)を補償するように選択することが好ましい
。これは、受信の時間的周波数掃引(調波パルスの周波数掃引とは幾分異なる)
を含むことができる。例えば、もし前記ガウス送信パルスを使用すれば、調波信
号は、幅(持続時間)T=(2α)1/2、帯域幅W=[2(α2+γ2)/α]1/2、及び式
19によって与えられる時間・帯域幅を用いて、
解能を得るために伸張された調波信号を最適に圧縮するためには、到来する調波
信号の分散(位相対周波数)を最良に取り消すように受信機位相対周波数を決定
することが好ましい。特に受信機帯域幅を調波信号帯域幅よりも狭くするような
パルス圧縮受信機を設計する別の方法は、得られる出力信号が帯域幅を減少させ
るように、適切な整合受信機及びフィルタを設計することである。これは、最適
な圧縮を、従って調波信号の軸方向分解能を保証する。
きである。これは、特に、受信した信号への基本波または望ましくない高調波及
び低調波の何れかからの寄与を最小にしながら、最大の軸方向分解能及び信号エ
ネルギが得られるように、送信及び受信スペクトルが好ましく設計されている整
数調波または低調波イメージングの場合に真である。
時間・帯域幅積及びモノトニック瞬時周波数変調に関する次の近似を使用する。
応用の場合に有用である2または4程度に小さい時間・帯域幅積にとっても有用
な結果を得ることができる。上述したように成形されたパルスに関連する調波ス
ペクトルは、
い。この理由から、この詳細な説明は単なる例示に過ぎず、限定を意図してはい
ない。本発明の範囲を限定するのは特許請求の範囲の独立及び従属請求項だけで
ある。
テムの概要図である。
織基本波信号の包絡線も示す。
波信号の包絡線も示す。
ルタのグラフである。
ルタのグラフである。
タのグラフである。
なる。
なる。
なる。
ていることが異なる。
ていることが異なる。
ていることが異なる。
ィングフィルタ及びエコー信号に関していることが異なる。
ィングフィルタ及びエコー信号に関していることが異なる。
ィングフィルタ及びエコー信号に関していることが異なる。
ィングフィルタ及びエコー信号に関していることが異なる。
を形成するために使用できる代替周波数変調関数対時間の図である。
トルの図である。
Claims (72)
- 【請求項1】 医療診断用超音波イメージング方法において、 (a)時間・帯域幅積が100より小さく1より大きいことを特色としている符
号化された基本超音波パルスを組織内へ送信するステップと、 (b)前記組織からN次調波エコー信号を受信するステップと、 (c)前記符号化された基本超音波パルスの速さの約N倍で変化する位相を有
する圧縮関数を用いて前記調波エコー信号を圧縮するステップと、 を含むことを特徴とする方法。 - 【請求項2】 前記Nは、1以外の整数に等しいことを特徴とする請求項1
に記載の方法。 - 【請求項3】 前記Nは、分数に等しいことを特徴とする請求項1に記載の
方法。 - 【請求項4】 (d)前記ステップ(a)及び(b)中、前記組織にコント
ラスト剤を添加しないように維持するステップ、 を更に含むことを特徴とする請求項1に記載の方法。 - 【請求項5】 前記ステップ(a)及び(b)中、前記組織内にコントラス
ト剤を供給するステップを更に含むことを特徴とする請求項1に記載の方法。 - 【請求項6】 前記ステップ(d)は、前記ステップ(a)、(b)、及び
(c)が遂行される全患者検査セッションを通して、前記組織にコントラスト剤
を添加しないように維持するステップからなることを特徴とする請求項1に記載
の方法。 - 【請求項7】 前記ステップ(a)は、フェーズドアレイ内に含まれる複数
の変換器素子からパルスを発射するステップからなることを特徴とする請求項1
に記載の方法。 - 【請求項8】 前記パルスは、時間・帯域幅積が50より小さいことを特徴と
する請求項1に記載の方法。 - 【請求項9】 前記パルスは、時間・帯域幅積が20より小さいことを特徴と
する請求項1に記載の方法。 - 【請求項10】 前記パルスは、時間・帯域幅積が10より小さいことを特徴
とする請求項1に記載の方法。 - 【請求項11】 前記パルスは、包絡線持続時間と、前記包絡線持続時間の
半分を中心として時間が対称である振幅変調関数とを更に特徴とする請求項1に
記載の方法。 - 【請求項12】 前記パルスは、包絡線持続時間と、前記包絡線持続時間の
半分を中心として時間が非対称である振幅変調関数とを更に特徴とする請求項1
に記載の方法。 - 【請求項13】 前記パルスは、実質的に均一である振幅変調関数を更に特
徴とする請求項1に記載の方法。 - 【請求項14】 前記パルスは、時間的に連続である振幅変調関数を更に特
徴とする請求項1に記載の方法。 - 【請求項15】 前記パルスは、時間的に不連続である振幅変調関数を更に
特徴とする請求項1に記載の方法。 - 【請求項16】 前記パルスは、実質的に非均一である振幅変調関数を更に
特徴とする請求項1に記載の方法。 - 【請求項17】 前記パルスは、スペクトルピークと、前記スペクトルピー
クを中心として周波数が対称である位相関数対周波数スペクトルとを更に特徴と
する請求項1に記載の方法。 - 【請求項18】 前記パルスは、スペクトルピークと、前記スペクトルピー
クを中心として周波数が非対称である位相関数対周波数スペクトルとを更に特徴
とする請求項1に記載の方法。 - 【請求項19】 前記パルスは、包絡線持続時間と、前記包絡線持続時間の
半分を中心として時間が対称である位相関数対周波数スペクトルとを更に特徴と
する請求項1に記載の方法。 - 【請求項20】 前記パルスは、包絡線持続時間と、前記包絡線持続時間の
半分を中心として時間が非対称である位相関数対周波数スペクトルとを更に特徴
とする請求項1に記載の方法。 - 【請求項21】 前記圧縮関数は、包絡線持続時間と、前記包絡線持続時間
の半分を中心として時間が対称である振幅変調関数とを更に特徴とする請求項1
に記載の方法。 - 【請求項22】 前記圧縮関数は、包絡線持続時間と、前記包絡線持続時間
の半分を中心として時間が非対称である振幅変調関数とを更に特徴とする請求項
1に記載の方法。 - 【請求項23】 前記圧縮関数は、実質的に均一である振幅変調関数を更に
特徴とする請求項1に記載の方法。 - 【請求項24】 前記圧縮関数は、時間的に連続である振幅変調関数を更に
特徴とする請求項1に記載の方法。 - 【請求項25】 前記圧縮関数は、時間的に不連続である振幅変調関数を更
に特徴とする請求項1に記載の方法。 - 【請求項26】 前記圧縮関数は、実質的に非均一である振幅変調関数を更
に特徴とする請求項1に記載の方法。 - 【請求項27】 前記圧縮関数は、スペクトルピークと、前記スペクトルピ
ークを中心として周波数が対称である位相関数対周波数スペクトルとを更に特徴
とする請求項1に記載の方法。 - 【請求項28】 前記圧縮関数は、スペクトルピークと、前記スペクトルピ
ークを中心として周波数が非対称である位相関数対周波数スペクトルとを更に特
徴とする請求項1に記載の方法。 - 【請求項29】 前記圧縮関数は、包絡線持続時間と、前記包絡線持続時間
の半分を中心として時間が対称である位相関数対周波数スペクトルとを更に特徴
とする請求項1に記載の方法。 - 【請求項30】 前記圧縮関数は、包絡線持続時間と、前記包絡線持続時間
の半分を中心として時間が非対称である位相関数対周波数スペクトルとを更に特
徴とする請求項1に記載の方法。 - 【請求項31】 前記圧縮関数は、時間的に連続である振幅変調関数を更に
特徴とする請求項1に記載の方法。 - 【請求項32】 前記圧縮関数は、時間的に不連続である振幅変調関数を更
に特徴とする請求項1に記載の方法。 - 【請求項33】 前記圧縮関数は、前記圧縮関数の入力における前記N次調
波エコー信号と振幅及び位相が整合していることを更に特徴とする請求項1に記
載の方法。 - 【請求項34】 前記圧縮関数は、前記圧縮関数の入力における前記N次調
波エコー信号と位相が整合し、前記送信された包絡線とは異なる包絡線を有する
ことを更に特徴とする請求項1に記載の方法。 - 【請求項35】 前記圧縮関数は、前記圧縮関数の入力における前記N次調
波エコー信号と位相が整合し、前記圧縮関数の入力における前記信号包絡線とは
異なる包絡線を有することを更に特徴とする請求項1に記載の方法。 - 【請求項36】 前記圧縮関数は、組織内の異なる範囲に対応するパルス受
信時間の関数として変化することを特徴とする請求項1に記載の方法。 - 【請求項37】 複数の関心周波数帯をデコードするために、複数のパルス
圧縮関数を使用することを特徴とする請求項1に記載の方法。 - 【請求項38】 (d)前記ステップ(a)のパルスと空間的に整列し、前
記ステップ(a)のパルスとは送信包絡線の位相及び極性の少なくとも一方が異
なる第2の符号化された基本超音波パルスを前記組織内へ送信するステップと、 (e)前記第2のパルスに応答して前記組織から第2のN次調波エコー信号を
受信するステップと、 (f)前記第2のパルスの速さの約N倍で変化する位相を有する圧縮関数を用
いて、前記第2の調波エコー信号を圧縮するステップと、 (g)前記ステップ(c)及び(f)からの圧縮された調波信号を組合せるス
テップと、 を更に含むことを特徴とする請求項1に記載の方法。 - 【請求項39】 (d)前記ステップ(a)のパルスと空間的に整列し、前
記ステップ(a)のパルスとは送信包絡線の位相及び極性の少なくとも一方が異
なる少なくとも2つの付加的な符号化された基本超音波パルスを前記組織内へ送
信するステップと、 (e)前記ステップ(d)の付加的なパルスに応答して前記組織から付加的な
第2のN次調波エコー信号を受信するステップと、 (f)前記付加的なパルスの速さの約N倍で変化する位相を有する圧縮関数を
用いて、前記付加的なN次調波エコー信号を圧縮するステップと、 (g)前記ステップ(c)及び(f)からの圧縮された調波信号を組合せるス
テップと、 を更に含むことを特徴とする請求項1に記載の方法。 - 【請求項40】 (d)前記ステップ(a)のパルスとは空間的に別個で、
前記ステップ(a)のパルスとは送信包絡線の位相及び極性の少なくとも一方が
異なる第2の符号化された基本超音波パルスを前記組織内へ送信するステップと
、 (e)前記第2のパルスに応答して前記組織から第2のN次調波エコー信号を
受信するステップと、 (f)前記第2のパルスの速さの約N倍で変化する位相を有する圧縮関数を用
いて、前記第2の調波エコー信号を圧縮するステップと、 (g)前記ステップ(c)及び(f)からの圧縮された調波信号を組合せるス
テップと、 を更に含むことを特徴とする請求項1に記載の方法。 - 【請求項41】 (d)前記ステップ(a)のパルスとは空間的に別個で、
前記ステップ(a)のパルスとは送信包絡線の位相及び極性の少なくとも一方が
異なる少なくとも2つの付加的な符号化された基本超音波パルスを前記組織内へ
送信するステップと、 (e)前記ステップ(d)の付加的なパルスに応答して前記組織から付加的な
第2のN次調波エコー信号を受信するステップと、 (f)前記付加的なパルスの速さの約N倍で変化する位相を有する圧縮関数を
用いて、前記付加的なN次調波エコー信号を圧縮するステップと、 (g)前記ステップ(c)及び(f)からの圧縮された調波信号を組合せるス
テップと、 を更に含むことを特徴とする請求項1に記載の方法。 - 【請求項42】 医療診断用超音波イメージングシステムにおいて、 フェーズドアレイ変換器プローブと、 モードと結合され、100より小さく1より大きい時間・帯域幅積を特色として
いる符号化された基本超音波パルスを前記プローブから組織内へ送信させるよう
な送信波形を前記モードに供給するように動作する送信機と、 前記プローブに結合され、前記組織からN次調波エコー信号を受信するように
動作する受信機と、 前記符号化された基本超音波パルスの約N倍で変化する位相を有する圧縮関数
を用いて前記調波エコー信号を圧縮するように動作する圧縮フィルタと、 を備えていることを特徴とするシステム。 - 【請求項43】 前記Nは、1以外の整数に等しいことを特徴とする請求項
42に記載の発明。 - 【請求項44】 前記Nは、分数に等しいことを特徴とする請求項42に記
載の発明。 - 【請求項45】 前記パルスは、時間・帯域幅積が50より小さいことを特徴
とする請求項42に記載の発明。 - 【請求項46】 前記パルスは、時間・帯域幅積が20より小さいことを特徴
とする請求項42に記載の発明。 - 【請求項47】 前記パルスは、時間・帯域幅積が10より小さいことを特徴
とする請求項42に記載の発明。 - 【請求項48】 前記パルスは、包絡線持続時間と、前記包絡線持続時間の
半分を中心として時間が対称である振幅変調関数とを更に特徴とする請求項42
に記載の発明。 - 【請求項49】 前記パルスは、包絡線持続時間と、前記包絡線持続時間の
半分を中心として時間が非対称である振幅変調関数とを更に特徴とする請求項4
2に記載の発明。 - 【請求項50】 前記パルスは、実質的に均一である振幅変調関数を更に特
徴とする請求項42に記載の発明。 - 【請求項51】 前記パルスは、時間的に連続である振幅変調関数を更に特
徴とする請求項42に記載の発明。 - 【請求項52】 前記パルスは、時間的に不連続である振幅変調関数を更に
特徴とする請求項42に記載の発明。 - 【請求項53】 前記パルスは、実質的に非均一である振幅変調関数を更に
特徴とする請求項42に記載の発明。 - 【請求項54】 前記パルスは、スペクトルピークと、前記スペクトルピー
クを中心として周波数が対称である位相関数対周波数スペクトルとを更に特徴と
する請求項42に記載の発明。 - 【請求項55】 前記パルスは、スペクトルピークと、前記スペクトルピー
クを中心として周波数が非対称である位相関数対周波数スペクトルとを更に特徴
とする請求項42に記載の発明。 - 【請求項56】 前記パルスは、包絡線持続時間と、前記包絡線持続時間の
半分を中心として時間が対称である位相関数対周波数スペクトルとを更に特徴と
する請求項42に記載の発明。 - 【請求項57】 前記パルスは、包絡線持続時間と、前記包絡線持続時間の
半分を中心として時間が非対称である位相関数対周波数スペクトルとを更に特徴
とする請求項42に記載の発明。 - 【請求項58】 前記圧縮関数は、包絡線持続時間と、前記包絡線持続時間
の半分を中心として時間が対称である振幅変調関数とを更に特徴とする請求項4
2に記載の発明。 - 【請求項59】 前記圧縮関数は、包絡線持続時間と、前記包絡線持続時間
の半分を中心として時間が非対称である振幅変調関数とを更に特徴とする請求項
42に記載の発明。 - 【請求項60】 前記圧縮関数は、実質的に均一である振幅変調関数を更に
特徴とする請求項42に記載の発明。 - 【請求項61】 前記圧縮関数は、実質的に非均一である振幅変調関数を更
に特徴とする請求項42に記載の発明。 - 【請求項62】 前記圧縮関数は、スペクトルピークと、前記スペクトルピ
ークを中心として周波数が対称である位相関数対周波数スペクトルとを更に特徴
とする請求項42に記載の発明。 - 【請求項63】 前記圧縮関数は、スペクトルピークと、前記スペクトルピ
ークを中心として周波数が非対称である位相関数対周波数スペクトルとを更に特
徴とする請求項42に記載の発明。 - 【請求項64】 前記圧縮関数は、包絡線持続時間と、前記包絡線持続時間
の半分を中心として時間が対称である位相関数対周波数スペクトルとを更に特徴
とする請求項42に記載の発明。 - 【請求項65】 前記圧縮関数は、包絡線持続時間と、前記包絡線持続時間
の半分を中心として時間が非対称である位相関数対周波数スペクトルとを更に特
徴とする請求項42に記載の発明。 - 【請求項66】 前記圧縮関数は、時間的に連続である振幅変調関数を更に
特徴とする請求項42に記載の発明。 - 【請求項67】 前記圧縮関数は、時間的に不連続である振幅変調関数を更
に特徴とする請求項42に記載の発明。 - 【請求項68】 前記圧縮関数は、前記圧縮関数の入力における前記N次調
波エコー信号と振幅及び位相が整合していることを更に特徴とする請求項42に
記載の発明。 - 【請求項69】 前記圧縮関数は、前記圧縮関数の入力における前記N次調
波エコー信号と位相が整合し、前記送信された包絡線とは異なる包絡線を有する
ことを更に特徴とする請求項42に記載の発明。 - 【請求項70】 前記圧縮関数は、前記圧縮関数の入力における前記N次調
波エコー信号と位相が整合し、前記圧縮関数の入力における前記信号包絡線とは
異なる包絡線を有することを更に特徴とする請求項42に記載の発明。 - 【請求項71】 前記圧縮関数は、組織内の異なる範囲に対応するパルス受
信時間の関数として変化することを特徴とする請求項42に記載の発明。 - 【請求項72】 複数の関心周波数帯をデコードするために、複数のパルス
圧縮関数を使用することを特徴とする請求項42に記載の発明。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/282,510 | 1999-03-31 | ||
US09/282,510 US6241674B1 (en) | 1999-03-31 | 1999-03-31 | Medical ultrasound diagnostic imaging method and system with nonlinear phase modulation pulse compression |
PCT/US2000/008133 WO2000057791A1 (en) | 1999-03-31 | 2000-03-28 | Ultrasonic imaging using pulse compression |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2002539881A true JP2002539881A (ja) | 2002-11-26 |
JP2002539881A5 JP2002539881A5 (ja) | 2007-05-10 |
Family
ID=23081827
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2000607546A Pending JP2002539881A (ja) | 1999-03-31 | 2000-03-28 | パルス圧縮を使用する超音波イメージング |
Country Status (5)
Country | Link |
---|---|
US (1) | US6241674B1 (ja) |
JP (1) | JP2002539881A (ja) |
AU (1) | AU3926400A (ja) |
DE (1) | DE10084427T1 (ja) |
WO (1) | WO2000057791A1 (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005334197A (ja) * | 2004-05-26 | 2005-12-08 | Aloka Co Ltd | 超音波診断装置 |
JP2010246640A (ja) * | 2009-04-13 | 2010-11-04 | Aloka Co Ltd | 超音波診断装置 |
Families Citing this family (39)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6375618B1 (en) * | 2000-01-31 | 2002-04-23 | General Electric Company | Enhanced tissue-generated harmonic imaging using coded excitation |
KR100362000B1 (ko) * | 2000-02-01 | 2002-11-22 | 주식회사 메디슨 | 변형된 골레이 코드를 이용한 펄스 압축 방식의 초음파영상 형성 방법 및 장치 |
US6599245B1 (en) * | 2000-06-27 | 2003-07-29 | Siemens Medical Solutions Usa, Inc. | Ultrasound transmission method and system for simulating a transmit apodization |
US6454714B1 (en) * | 2000-10-20 | 2002-09-24 | Koninklijke Philips Electronics N.V. | Ultrasonic harmonic flash suppression |
JP2002306477A (ja) * | 2001-04-11 | 2002-10-22 | Ge Medical Systems Global Technology Co Llc | 超音波送受信方法、超音波送受信装置、超音波撮影方法および超音波撮影装置 |
JP4666815B2 (ja) * | 2001-06-06 | 2011-04-06 | アロカ株式会社 | 超音波診断装置 |
KR100406099B1 (ko) * | 2001-09-05 | 2003-11-14 | 주식회사 메디슨 | 다단계 구조의 펄스 압축기를 이용한 초음파 영상 형성 장치 및 방법 |
KR100406097B1 (ko) | 2001-12-26 | 2003-11-14 | 주식회사 메디슨 | 가중 쳐프 신호를 이용한 초음파 영상 형성 장치 및 방법 |
JP3910860B2 (ja) * | 2002-02-05 | 2007-04-25 | 株式会社日立メディコ | 超音波撮像装置 |
US6632177B1 (en) | 2002-05-01 | 2003-10-14 | Acuson Corporation | Dual process ultrasound contrast agent imaging |
US20090062644A1 (en) * | 2002-06-07 | 2009-03-05 | Mcmorrow Gerald | System and method for ultrasound harmonic imaging |
US7094204B2 (en) * | 2002-08-23 | 2006-08-22 | Siemens Medical Solutions Usa, Inc. | Coded excitation imaging for use with bipolar, unipolar and other waveforms |
US6589177B1 (en) * | 2002-11-15 | 2003-07-08 | Koninklijke Philips Electronics N.V. | Method and apparatus for obtaining B-flow and B-mode data from multiline beams in an ultrasound imaging system |
US7035040B2 (en) * | 2003-05-16 | 2006-04-25 | Imation Corp. | Sequenced time-based servo techniques |
US6960169B2 (en) * | 2003-05-19 | 2005-11-01 | Siemens Medical Solutions Usa, Inc. | Spread spectrum coding for ultrasound contrast agent imaging |
US6952317B2 (en) * | 2003-06-17 | 2005-10-04 | Imation Corp. | Amplitude-based servo patterns for magnetic media |
EP1660908A1 (en) * | 2003-08-26 | 2006-05-31 | Koninklijke Philips Electronics N.V. | System and method for ultrasound pulse shaping and output power adjustment using multiple drive pulses |
US7037265B2 (en) * | 2003-09-09 | 2006-05-02 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for tissue harmonic imaging with natural (tissue) decoded coded excitation |
US7038872B2 (en) * | 2003-11-10 | 2006-05-02 | Imation Corp. | Servo patterns with inherent track ID |
US7038871B2 (en) * | 2003-11-10 | 2006-05-02 | Imation Corp. | Multi-band servo patterns with inherent track ID |
US7142381B2 (en) * | 2003-11-10 | 2006-11-28 | Imation Corp. | Servo writing devices for creating servo patterns with inherent track ID |
JP4602993B2 (ja) * | 2004-01-16 | 2010-12-22 | ボストン サイエンティフィック リミテッド | 医用撮像のための方法及び装置 |
US7095583B2 (en) * | 2004-06-02 | 2006-08-22 | Imation Corp. | Dual mode servo pattern |
US7578790B2 (en) * | 2004-07-20 | 2009-08-25 | Boston Scientific Scimed, Inc. | Systems and methods for detecting and presenting textural information from medical images |
US20060036147A1 (en) * | 2004-07-20 | 2006-02-16 | Scimed Life Systems, Inc. | Systems and methods for detecting and presenting textural information from medical images |
US20060173318A1 (en) * | 2004-07-20 | 2006-08-03 | Scimed Life Systems Inc. | Systems and methods for detecting and presenting textural information from medical images |
US7199958B2 (en) * | 2004-08-25 | 2007-04-03 | Imation Corp. | Servo head with varying write gap width |
US7466510B2 (en) * | 2005-06-03 | 2008-12-16 | Imation Corp. | Distributed servo patterns for data storage media |
US7379254B2 (en) * | 2005-06-29 | 2008-05-27 | Imation Corp. | Mixed frequency amplitude-based servo pattern |
US7946988B2 (en) * | 2005-12-13 | 2011-05-24 | Siemens Medical Solutions Usa, Inc. | Medical diagnostic ultrasound contrast agent destruction with reduced biological effect |
US20070196282A1 (en) * | 2006-02-21 | 2007-08-23 | Siemens Medical Solutions Usa, Inc. | Medical diagnostic ultrasound with temperature-dependent contrast agents |
US7535797B2 (en) * | 2006-06-20 | 2009-05-19 | Rehabtek | High-resolution ultrasound displacement measurement apparatus and method |
US7436622B2 (en) * | 2006-07-31 | 2008-10-14 | Imation Corp. | Concurrent servo and data track writing |
US7967753B2 (en) * | 2006-08-01 | 2011-06-28 | Stichting Voor de Technische Wetenschappen of Van Vollenhovenlaan | Pulse inversion sequences for nonlinear imaging |
JP4928989B2 (ja) * | 2007-03-07 | 2012-05-09 | 株式会社東芝 | 超音波診断装置および超音波診断装置の制御プログラム |
GB2457240B (en) * | 2008-02-05 | 2013-04-10 | Fujitsu Ltd | Ultrasound probe device and method of operation |
TWI521223B (zh) | 2012-02-13 | 2016-02-11 | 國立臺灣科技大學 | 超音波成像裝置及其方法 |
US20130281858A1 (en) * | 2012-04-19 | 2013-10-24 | Industrial Technology Research Institute | Method for increasing depth of field and ultrasound imaging system using the same |
GB201222882D0 (en) | 2012-12-19 | 2013-01-30 | Univ Leeds | Ultrasound generation |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0779973A (ja) * | 1993-09-20 | 1995-03-28 | Ge Yokogawa Medical Syst Ltd | 超音波分散圧縮送受信方法及び超音波分散圧縮送受信装置 |
JPH09187457A (ja) * | 1996-01-12 | 1997-07-22 | Ge Yokogawa Medical Syst Ltd | 超音波撮像方法および装置 |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2478822A1 (fr) | 1980-03-18 | 1981-09-25 | Thomson Csf | Systeme de detection active au moyen d'emissions multiples simultanees |
FR2478824A1 (fr) | 1980-03-21 | 1981-09-25 | Thomson Csf | Systeme d'imagerie acoustique |
FR2492109B1 (fr) | 1980-10-10 | 1985-07-05 | Thomson Csf | Systeme d'imagerie a emissions multiples et simultanees |
FR2493528A1 (fr) | 1980-11-04 | 1982-05-07 | Thomson Csf | Systeme de detection multivoies a emission diversifiee |
US5014712A (en) | 1989-12-26 | 1991-05-14 | General Electric Company | Coded excitation for transmission dynamic focusing of vibratory energy beam |
US5224482A (en) | 1991-04-08 | 1993-07-06 | Kabushiki Kaisha Toshiba | Ultrasound high velocity flow correlation measurement using coded pulses |
US5142649A (en) | 1991-08-07 | 1992-08-25 | General Electric Company | Ultrasonic imaging system with multiple, dynamically focused transmit beams |
US5522393A (en) | 1994-05-24 | 1996-06-04 | Duke University | Multi-dimensional real-time ultrasonic blood flow imaging apparatus and method |
US5454372A (en) | 1994-06-17 | 1995-10-03 | Siemens Medical Systems, Inc. | Angle independent doppler in ultrasound imaging |
US5675554A (en) | 1994-08-05 | 1997-10-07 | Acuson Corporation | Method and apparatus for transmit beamformer |
US5608690A (en) | 1995-03-02 | 1997-03-04 | Acuson Corporation | Transmit beamformer with frequency dependent focus |
US5851187A (en) | 1997-10-01 | 1998-12-22 | General Electric Company | Method and apparatus for ultrasonic beamforming with spatially encoded transmits |
US5964706A (en) * | 1998-03-18 | 1999-10-12 | General Electric Company | Method and apparatus for pulsed doppler imaging using coded excitation on transmit and pulse compression on receive |
US6095977A (en) * | 1998-03-26 | 2000-08-01 | Hall; Anne Lindsay | Method and apparatus for color flow imaging using Golay-coded excitation on transmit and pulse compression on receive |
US5980459A (en) * | 1998-03-31 | 1999-11-09 | General Electric Company | Ultrasound imaging using coded excitation on transmit and selective filtering of fundamental and (sub)harmonic signals on receive |
US6050947A (en) | 1998-04-20 | 2000-04-18 | General Electric Company | Method and apparatus for harmonic tissue imaging and contrast imaging using coded transmission |
US5961463A (en) * | 1998-08-24 | 1999-10-05 | General Electric Company | Nonlinear imaging using orthogonal transmit and receive codes |
US6010456A (en) | 1998-12-30 | 2000-01-04 | General Electric Company | Method and apparatus for acoustic subtraction imaging using linear and nonlinear ultrasonic images |
-
1999
- 1999-03-31 US US09/282,510 patent/US6241674B1/en not_active Expired - Lifetime
-
2000
- 2000-03-28 JP JP2000607546A patent/JP2002539881A/ja active Pending
- 2000-03-28 WO PCT/US2000/008133 patent/WO2000057791A1/en active Application Filing
- 2000-03-28 DE DE10084427T patent/DE10084427T1/de not_active Withdrawn
- 2000-03-28 AU AU39264/00A patent/AU3926400A/en not_active Abandoned
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0779973A (ja) * | 1993-09-20 | 1995-03-28 | Ge Yokogawa Medical Syst Ltd | 超音波分散圧縮送受信方法及び超音波分散圧縮送受信装置 |
JPH09187457A (ja) * | 1996-01-12 | 1997-07-22 | Ge Yokogawa Medical Syst Ltd | 超音波撮像方法および装置 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005334197A (ja) * | 2004-05-26 | 2005-12-08 | Aloka Co Ltd | 超音波診断装置 |
JP4499477B2 (ja) * | 2004-05-26 | 2010-07-07 | アロカ株式会社 | 超音波診断装置 |
JP2010246640A (ja) * | 2009-04-13 | 2010-11-04 | Aloka Co Ltd | 超音波診断装置 |
Also Published As
Publication number | Publication date |
---|---|
AU3926400A (en) | 2000-10-16 |
US6241674B1 (en) | 2001-06-05 |
DE10084427T1 (de) | 2002-03-07 |
WO2000057791A1 (en) | 2000-10-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2002539881A (ja) | パルス圧縮を使用する超音波イメージング | |
US6491631B2 (en) | Harmonic golay-coded excitation with differential pulsing for diagnostic ultrasound imaging | |
US6213947B1 (en) | Medical diagnostic ultrasonic imaging system using coded transmit pulses | |
JP3423935B2 (ja) | 超音波撮像方法及び改良ゴールレイコードを用いたパルス圧縮方式に基づく装置 | |
Chiao et al. | Coded excitation for diagnostic ultrasound: A system developer's perspective | |
EP1122556B1 (en) | Enhanced tissue-generated harmonic imaging using coded excitation | |
AU688116B2 (en) | Ultrasonic spectral contrast imaging | |
US6050944A (en) | Method and apparatus for frequency control of an ultrasound system | |
US6440075B1 (en) | Ultrasonic diagnostic imaging of nonlinearly intermodulated and harmonic frequency components | |
EP1121901A2 (en) | Angle independent ultrasound volume flow measurement | |
MXPA96001754A (en) | Images by spectro ultrason contrast | |
JPH10295694A (ja) | 超音波イメージング・システムの作動方法 | |
JP4642977B2 (ja) | 超音波診断装置及び超音波イメージング方法 | |
US7037265B2 (en) | Method and apparatus for tissue harmonic imaging with natural (tissue) decoded coded excitation | |
KR100406097B1 (ko) | 가중 쳐프 신호를 이용한 초음파 영상 형성 장치 및 방법 | |
WO2003090625A1 (en) | Ultrasonograph | |
US5957852A (en) | Ultrasonic harmonic imaging system and method | |
US6213946B1 (en) | Methods and apparatus for speckle reduction by orthogonal pulse compounding in medical ultrasound imaging | |
JP2003265466A (ja) | 超音波診断装置 | |
JP2005131407A (ja) | 単一送信ゴレイ符号化励起のための方法及び装置 | |
JP4728466B2 (ja) | 物質をイメージングする方法及びイメージング・システム | |
US6440074B1 (en) | Ultrasonic diagnostic imaging with nonlinearly intermodulated frequency components | |
Ramalli et al. | Real-time base-band pulse compression imaging | |
Trucco et al. | Harmonic beamforming: a new approach to removing the linear contribution from harmonic imaging | |
Trucco et al. | Improving harmonic imaging by a novel beamforming process for ultrasound medical scanners |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070307 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20070307 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20100316 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100401 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20100909 |