JP2002045359A - フィルター係数決定方法、超音波撮像システム及びその方法 - Google Patents

フィルター係数決定方法、超音波撮像システム及びその方法

Info

Publication number
JP2002045359A
JP2002045359A JP2001180172A JP2001180172A JP2002045359A JP 2002045359 A JP2002045359 A JP 2002045359A JP 2001180172 A JP2001180172 A JP 2001180172A JP 2001180172 A JP2001180172 A JP 2001180172A JP 2002045359 A JP2002045359 A JP 2002045359A
Authority
JP
Japan
Prior art keywords
signal
ultrasonic imaging
spread spectrum
ultrasonic
pulse compression
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2001180172A
Other languages
English (en)
Inventor
Te Kyon Son
テ キョン ソン
Yan Mo Yuu
ヤン モ ユウ
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Samsung Medison Co Ltd
Original Assignee
Medison Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Medison Co Ltd filed Critical Medison Co Ltd
Publication of JP2002045359A publication Critical patent/JP2002045359A/ja
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/895Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques characterised by the transmitted frequency spectrum
    • G01S15/8954Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques characterised by the transmitted frequency spectrum using a broad-band spectrum
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Acoustics & Sound (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radiology & Medical Imaging (AREA)
  • Veterinary Medicine (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Medical Informatics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biophysics (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

(57)【要約】 【課題】 許容サイドローブを有するスペクトル拡散信
号を用いるパルス圧縮技法に基づいた超音波撮像システ
ム及び方法、このシステムで使用されるFIRフィルタ
ー係数決定方法を提供する。 【解決手段】 予め決められた第1スペクトル拡散信号
を超音波信号に変換してオブジェクトに発射し、該オブ
ジェクトから反射された反射信号を受け取って反射信号
に対しパルス圧縮を行ってパルス圧縮信号を生成し、こ
のパルス圧縮信号を処理して受信集束信号を生成し、こ
の受信集束信号からオブジェクトの映像を生成する。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は、超音波撮像システ
ムに関し、特に、スペクトル拡散信号と効率的なハード
ウェア構造を有する有限インパルス応答(FIR)フィル
ターとを利用するパルス圧縮技法に基づく超音波撮像シ
ステムに関する。
【0002】
【従来の技術】通常、医療用超音波撮像システムは、短
パルスの超音波を人体へ発射し、それから反射される信
号を受け取って人体関連情報を得る。図1は、従来の短
パルス超音波撮像システム100のブロック図であっ
て、複数のトランスジューサよりなるトランスジューサ
列1、パルサー11、送信(Tx)集束遅延メモリ14、
送受信(Tx/Rx)スイッチ21、レシーバ31、ビー
ム形成器37、受信(Rx)集束遅延調整器36、信号処
理器41及びスキャンコンバータ42を含む。
【0003】Tx集束遅延メモリ14は、トランスジュ
ーサ列1を通じてオブジェクト(例えば、人体)内に伝送
されるべき超音波パルスの遅延パターンを格納し、該遅
延パターンに対応する2進系列をパルサー11に供給す
る。
【0004】各トランスジューサに対する遅延パターン
を決める方法として、通常、人体内の所定地点上に超音
波パルスのエネルギーを集束する固定集束技法が用いら
れている。最近、送信時には固定集束技法を、受信時に
は動的集束技法を使用して生じる分解能低下の問題を解
決するための技法の一つとして、合成開口技法(synthet
ic aperture technique)に対する研究が行われてい
る。この合成開口技法によれば、超音波の送信時に一つ
またはそれ以上のトランスジューサを用いることがで
き、送受信両側で両方向動的集束が可能である。このよ
うな合成開口技法を用いることによって、分解能は高
め、信号対雑音比(SNR)は低下する。
【0005】パルサー11はTx集束遅延メモリ14か
らの2進系列に応じて、増幅信号(例えば、+80V又は
-80V)をトランスジューサ列1へ供給するバイポーラ
形パルサーである。パルサー11からの出力は予め定め
られた大きさの電圧であって、該遅延パターンによって
決定される時刻にてトランスジューサ列1の各トランス
ジューサに印加される。
【0006】トランスジューサ列1はパルサー11から
の出力電圧に応じて、超音波パルスをオブジェクト内に
発射する。任意の時刻にて超音波パルスの送信の際、ト
ランスジューサ列1のうちの一部が選択的に用いられ得
る。例えば、トランスジューサ列1がN個(128個)の
トランスジューサよりなる場合、一回送信の際に、64
個のトランスジューサのみが開口を通じて超音波を送信
するようにする。その後、トランスジューサ列1は超音
波パルスの送信の際該オブジェクトから反射されるパル
ス信号を受け取る。
【0007】Tx/Rxスイッチ21は、パルサー11
から出る高電圧に対してレシーバ31を分離させる送受
切換え器として機能を果たし、送信の際にはトランスジ
ューサ列1をパルサー11に接続させ、受信の際にはト
ランスジューサ列1をレシーバ31に接続させる。
【0008】レシーバ31は受信信号を増幅する前置増
幅器と、超音波の伝搬の際生じる減衰を補償するための
時間利得補償器(TGC)と、増幅された受信信号をディ
ジタル信号に変換するアナログディジタル変換器とから
構成される。
【0009】ビーム形成器37はRx集束遅延調整器3
6からの遅延パターンに基づいて、レシーバ31から入
力された該当ディジタル信号に対して受信集束を行う。
【0010】信号処理器41は包絡線検波、対数補償な
どの信号処理を行い、Bモードの映像信号を生成する。
【0011】スキャンコンバータ42は、Bモードの映
像信号を表示装置(図示せず)上に表示可能な信号に変換
する。
【0012】超音波をゴム、軟組織のように減衰の激し
い媒体内へ伝搬する時、超音波の電力が減少することと
なり、短パルス方式の撮像システムは人体内の目標物に
対する正確な情報を得ることが困難である。
【0013】医療用超音波撮像システム100では、送
信される短パルスのピーク電圧が増加すれば、人体に損
傷をもたらし得るため、送信パルスの電力を増加させ受
信信号の電力を増加させることができない。
【0014】一方、レーダー装置で用いられるパルス圧
縮技法は、送信パルスのピーク電圧を増加させる代わり
に、送信パルスの平均電力を増加させて超音波撮像シス
テムの信号対雑音比を向上させることができる。通常、
そのようなパルス圧縮技法を用いる撮像システムは短パ
ルスの代わりに、長期波形信号(即ち、長いパルス)を人
体に送信してSNRを向上させる。
【0015】短パルスを用いる従来の医療用撮像システ
ム100では、高電圧の短パルスを用いるため、超音波
伝搬方向への分解能は選択された超音波トランスジュー
サのインパルス応答によって左右される。しかし、パル
ス圧縮技法を用いる撮像システムの場合に、分解能は超
音波トランスジューサと長パルスとの間のたたみ込みに
より決定される。
【0016】パルス圧縮方式の撮像システムにおいて、
超音波レシーバはFIRフィルターを有するパルス圧縮
器を用いて、短パルス方式のピーク電圧に比べ低い電圧
を有する長パルス信号を送信することによって信号対雑
音比を効果的に増加させることができる。
【0017】長パルスを用いる超音波撮像システムの性
能は、用いられた長パルス信号の特性によって左右され
る。詳記すると、最終の画質は長パルス信号と超音波ト
ランスジューサとの間の周波数特性の関係によって決定
される。また、システムの性能はFIRフィルターまた
はパルス圧縮器をどのように具現するかによって左右さ
れる。
【0018】また、動的受信集束のためには、パルス圧
縮器がすべてのチャネルに使用されなければならないの
で、受信部のハードウェア複雑さはパルス圧縮器のハー
ドウェア構造によって左右される。
【0019】そのような長パルス信号を用いる超音波撮
像システムでは、例えば、チャープ信号(線形周波数変
調信号)のようなスペクトル拡散信号を長パルス信号と
して使用することができる。そのようなチャープ信号
は、帯域幅の制限された超音波撮像システムのトランス
ジューサのスペクトルと整合する周波数特性を有する。
また、スペクトル拡散信号が従来のFIRフィルターを
通過すると、メインローブより略-13dB位低い尖頭
サイドローブが生成される。しかしながら、医療用超音
波撮像システムに用いるためには、パルス圧縮器を通し
た出力信号のスペクトルにおけるサイドローブがメイン
ローブに比べて-50dB以下でなければならないの
で、従来のスペクトル拡散信号は医療用超音波撮像シス
テムに不適合であるという不都合がある。
【0020】
【発明が解決しようとする課題】従って、本発明の目的
は、許容サイドローブを有するスペクトル拡散信号を用
いるパルス圧縮技法に基づいた超音波撮像システム及び
その方法を提供することにある。
【0021】本発明の他の目的は、効率的な構造を有す
る有限インパルス応答(FIR)フィルターとスペクトル
拡散信号を用いる超音波撮像システムで使用されるFI
Rフィルター係数を決定する方法を提供することにあ
る。
【0022】
【課題を解決するための手段】上記の目的を達成するた
めに、本発明の一好適実施例によれば、オブジェクトの
映像を撮像する超音波撮像方法であって、予め決められ
た第1スペクトル拡散信号を超音波信号に変換して前記
オブジェクトに発射する第1ステップと、前記オブジェ
クトから反射された前記超音波信号の反射信号を受け取
って、前記反射信号に対してパルス圧縮を行ってパルス
圧縮信号を生成する第2ステップと、前記パルス圧縮信
号を処理して受信集束信号を生成し、前記受信集束信号
から前記オブジェクトの映像を生成する第3ステップと
を含む超音波撮像方法が提供される。
【0023】本発明の他の好適実施例によれば、オブジ
ェクトの映像を撮像する超音波撮像システムであって、
予め決められたスペクトル拡散信号に応じて、前記オブ
ジェクトに超音波信号を発射する一つまたはそれ以上の
トランスジューサと、前記オブジェクトから反射される
反射信号を受け取って、 前記反射信号に対してパルス
圧縮を行ってパルス圧縮信号を生成するパルス圧縮手段
と、前記パルス圧縮信号を処理して受信集束信号を生成
し、前記受信集束信号に基づいて前記オブジェクトの映
像を生成する手段とを含む超音波撮像システムが提供さ
れる。
【0024】
【発明の実施の形態】以下、本発明の好適実施例につい
て、図面を参照しながらより詳しく説明する。
【0025】図2は、本発明による超音波撮像システム
200の概略的なブロック図である。本発明の好適実施
例によれば、超音波の送信の際、重みつきチャープ信号
を用いているが、それと類似した特性を有するスペクト
ル拡散信号などの信号を用いてもよい。また、本発明の
超音波撮像システム200は受信ビーム形成以前にRF
データに基づいたパルス圧縮を施して受信集束を行うこ
とを特徴とする。従って、本発明によれば、受信集束以
後にパルス圧縮を行ったとき、近距離フィールドで超音
波伝搬方向に生じ得るサイドローブを除去することがで
き、メインローブの幅が増加することを防止することが
できる。
【0026】図2において、超音波撮像システム200
はパルサー12、Txパターンメモリ13、Tx集束遅
延メモリ14、フィルター係数制御器15、フィルター
係数メモリ16、Tx/Rxスイッチ21、レシーバ3
1、パルス圧縮器35、Rx集束遅延調整器36、ビー
ム形成器37、信号処理器41及びスキャンコンバータ
42から成る。パルス圧縮器35はRFメモリ及びFI
Rフィルター34を有する。
【0027】Txパターンメモリ13はTx信号パター
ンを格納し、Tx集束遅延メモリはTx信号に対する遅
延パターンを格納する。Txパターンメモリ13に格納
されたTx信号パターンはパルサー12に伝送され、こ
のパルサー12は遅延パターンによって遅延されたパル
スを有する。
【0028】本発明の一実施例によれば、Txパターン
メモリ13及びTx集束遅延メモリ14は各々予め定め
られたTx信号パターン及び予め定められた遅延パター
ンを格納する。パルサ12への入力を供給するものとし
て、遅延されたTx信号パターンを生成し得る回路また
はソフトウェアを用いてもよい。
【0029】使用ウィンドウによって、パルサー12は
Txパターンメモリ13から受け取ったTx信号パター
ンを増幅する。四角形ウィンドウが用いられる場合、パ
ルサー12の代わりに、Tx/Rxスイッチ21を通じ
てバイポーラパルスをトランスジューサ列1に供給する
バイポーラパルサーを用いることができる。四角形ウィ
ンドウ以外のものを用いる場合は、パルサー12は増幅
されたバイポーラパルスをトランスジューサ列1に供給
するためにバイポーラパルサー及び線形増幅器を有し得
る。図1に示したように、固定集束技法または合成開口
技法がチャープ信号の伝送の際に用いられる。
【0030】Tx/Rxスイッチ21は図1中のスイッ
チと同様に、パルサー12から出る高電圧に対してレシ
ーバ31を分離させる送受切換え器として機能を果た
す。
【0031】レシーバ31が受信した信号は、サイドロ
ーブが大きすぎて映像形成に適合でない。従って、結果
映像の分解能が従来の短パルス方式の撮像システムの分
解能と互換されるように、パルス圧縮器35を用いてレ
シーバ31からの信号を処理する。
【0032】詳記すると、パルス圧縮器35は、レシー
バ31から受信したRF信号をRFメモリ33に格納す
る。各チャネルに対するRF信号はRFメモリ33に格
納されているため、FIRフィルター34はハードウェ
アまたはソフトウェアによって具現されることができ
る。また、スペクトル拡散信号の代わりに通常の短パル
スを使用する場合、図2中の撮像システム200はFI
Rフィルター34をバイパスすることによって従来の超
音波撮像システムの機能を行うこともできる。
【0033】FIRフィルター34はフィルター係数制
御器15から供給されるフィルター係数に応じて、レシ
ーバ31から受信した信号、即ち、RFメモリ33に格
納された信号をフィルタリングしてパルス圧縮信号を出
力する。即ち、長パルスが送信されれば、FIRフィル
ター34は短パルスが送信されたことと同一の出力信号
を発生する。
【0034】フィルター係数メモリ16は複数組のフィ
ルター係数を格納する。各係数組は原チャープ信号及び
ウィンドウ関数の特性と超音波映像形成の際に用いられ
る周波数成分(例えば、基本周波数成分または高調波成
分)とによって予め定められる。このような周波数成分
は外部から使用者の入力によって選択され得る。
【0035】フィルター係数制御器15は、フィルター
係数メモリ16に格納された係数組のうちのいずれかを
選択してFIRフィルター34に供給する。
【0036】ビーム形成器37はRx集束遅延調整器3
6から供給されたRx集束遅延を用いてRx集束を行
う。Rx集束の後、信号処理器41は包落線検波及びロ
グ補償を行ってBモードの映像信号を生成する。スキャ
ンコンバータ42はBモードの映像信号を変換して表示
装置(図示せず)上にディスプレイされるようにする。
【0037】前述したように、本発明の好適実施例によ
れば、超音波送信のために重みつきチャープ信号が用い
られる。重みつきチャープ信号は実際上一定な振幅及び
周波数を有することで、制限された周波数帯域で線形的
に時変する。チャープ信号のスペクトルは医療用撮像シ
ステムに用いられる狭帯域トランスジューサの特性とよ
く整合する。チャープ信号の周波数はドップラーシフト
によって変化され、このような周波数変化はFIRフィ
ルターにて瞬間測定され目標オブジェクトの速度を表
す。時変する角周波数がω=ω0+μtであるとき、チャー
プ信号は下記式のように表現され得る。
【数1】 ここで、Aは正の整数であり、ω0はチャープ信号の中
心角周波数を示し、定数μは時間間隙Tの間角周波数の
変化量Δω/Tを示し、w1 (t)は送信されるチャープ信
号の包落線を示す。
【0038】通常のチャープ信号の場合、w1 (t)は矩形
波信号、即ち、w1 (t)=rect(t/T)であって、通常のF
IRフィルターを通過するとメインローブ周りに略−1
3dBの尖頭サイドローブを生成するので、医療用超音
波撮像システムには適合でない。即ち、そのような通常
のチャープ信号を用いる場合、医療用映像を得るのに必
要な最小信号対雑音比を得ることができなくなる。
【0039】従って、本発明はウィンドウ関数を用いて
チャープ信号の包落線に重みを付けて、フィルター係数
を変更させてデータ送受信を行うシステムを提供する。
チャープ信号は上記式(1)の四角形ウィンドウ関数以外
にも、低いサイドローブ及び狭いメインローブを生成し
得るウィンドウ関数を用いて重み付きが加えられる。
【0040】図3は従来のチャープ信号のグラフを示
し、図4は本発明で使用されるハミング・ウィンドウ(H
amming-window)、ハンニング・ウィンドウ(Hanning-win
dow)、ブラックマン・ウィンドウ(Blackman-window)、
矩形ウィンドウ(Rectangular-window)を示す模式図であ
り、図5はハンニング・ウィンドウによって重み付きさ
れたチャープ信号を示すグラフである。図4及び図5に
おいて、Lはフィルター係数の個数を示す。重みつきチ
ャープ信号をトランスジューサに供給するためには、二
つの離散値だけを供給するパルサ11よりは、前述のよ
うに連続的に信号を増幅し得るパルサ12が必要であ
る。
【0041】重みつきチャープ信号のフィルター出力信
号は下記式(2)の通りチャープ信号を変更するウィンド
ウ関数w1 (t)とフィルター係数を変更するウィンドウ関
数w2(t)とによって決定される。
【数2】
【0042】その後、フィルター係数は、FIRフィル
ターに使用されたウィンドウ関数及び原チャープ信号の
特性と、受信信号に用いられた周波数成分が基本周波数
成分か高調波成分かとによって計算される。
【0043】重みつきチャープ信号は通常のチャープ信
号の包落線を変更して得られるため、送信電力効率が低
下され、パルス圧縮技法を通じて得るようになる信号対
雑音比の増加を減少させる。従って、トランスジューサ
1における入力電力に対する出力電力の比で定義される
伝送電力効率(TPE)を考慮してウィンドウ関数を選択
しなければならない。また、超音波映像の分解能はパル
ス圧縮信号のサイドローブのレベル及びメインローブの
幅によって決定されるので、最低レベルのサイドローブ
及び最低幅のメインローブを齎すウィンドウ関数を選択
することが望ましい。即ち、本発明の医療用超音波撮像
システムに採用される送受信用ウィンドウ関数は、フィ
ルター出力信号のメインローブ幅、最大サイドローブ及
び送信電力効率を考慮して選択されなければならない。
【0044】一方、前述したFIRフィルター34は各
チャネル別に備えなければならない。重みつきチャープ
信号を利用する場合や変形フィルター係数を使用する場
合に、FIRフィルター34の係数は連続的な値を有す
ることとなる。このため、乗算器及び加算器を有するF
IRフィルター34の構造は複雑になって、多数のフィ
ルターを含む撮像システムの全体的な構造が複雑にな
る。
【0045】図6は、従来の整合フィルターの形態で具
現されたFIRフィルターの構造を示した模式図であ
る。送信信号の持続時間が長いほど整合フィルタの長さ
も増加し、FIRフィルターのハードウェア構造が相当
に複雑になる。図7は、本発明によるFIRフィルター
の構造を示した模式図である。フィルター係数の連続的
な値によるFIRフィルター34のハードウェア複雑さ
は主に乗算器によって齎される。従って、本発明では発
生的アルゴリズムを利用して乗算器が除去されたFIR
フィルターを設計することによって受信部のハードウェ
ア複雑さを減らすことができる。相関値計算に使用され
るフィルター係数を2の自乗またはその和で表現できる
値に量子化される。そのような量子化によって、FIR
フィルター34は図6及び図7に示すように加算器及び
シフターだけで具現され得る。
【0046】図8は、本発明の実施例によってフィルタ
ー係数の決定に用いられる発生的アルゴリズムを示すフ
ロー図である。発生的アルゴリズムは生物学的進化をモ
デリングするアルゴリズムであって、目的関数によって
個体群のうちで染色体を選択し発生的演算を適用して次
の個体群を生成する技法である(D.E.Goldbe
rg,「サーチ、最適化、マシン学習における発生的ア
ルゴリズム」、Addison−Wesley Pub
lishingCompany、Inc.(198
9))。この発生的アルゴリズムにおいて、個体群は染色
体の集合を示し、一般的に染色体は固定長さの列であ
る。初期にランダムに発生された各染色体集合は目的関
数により評価されて、その中の一部が選択演算によって
次の染色体の親として選択されることによって、次の世
代の良い染色体が生成される。選択された染色体は互い
に有している染色体の一部を交換(交配)するかまたは交
配により得ることが難しい染色体を得るための突然変異
を通じて次の個体群の新たな染色体を生成する。
【0047】例えば、フィルター係数の個数が512個
で、一個の係数は8ビットで表現される時、512個の
係数全体を表す4096ビットで構成されるビット列は
一つの染色体に対応する。染色体はエラー自乗平均(M
SE)、最小局所最大(minimum local m
axmum)などの目的関数により評価される。本発明
の一実施例によれば、目的関数としてのMSEは次の通
り与えられる。
【数3】 ここで、Rmm(i)は通常のFIRフィルターを使用し
た場合のフィルター出力であり、Rmd(i)は本発明によ
って設計されたFIRフィルターを使用した場合のフィ
ルター出力である。Nはフィルター係数の個数を示す。
【0048】図8を参照すれば、発生的アルゴリズムを
用いてFIRフィルター34の係数を求める過程が述べ
られている。ステップS101にて、各ビットが乗算器
を加算器及びシフターで代替させるために量子化された
一セットのビット列の初期値がランダムに発生される。
ステップS102にて、発生されたビット列の各々に対
するMSEが上記式(3)を使用して求められる。
【0049】ステップS103にて、各ビット列に対し
て、求められたMSEが医療分野に適合であるかを判断
する。ステップS104にて、MSEが医療分野に適合
であれば該当ビット列を相関係数として選択し、プロセ
スを終了する。一方、MSEが医療分野に適合でなけれ
ば、ステップS105にて発生されたビット列の一部を
選択し、選択されたビット列はステップS106及びS
107にて交配または変異過程を経ることとなる。続け
て、制御プロセスはステップS102に戻して次のビッ
ト列に対するMSEを計算する。
【0050】前述のように、各フィルター係数を2の自
乗またはその和で表現されることができる値に量子化す
れば、FIRフィルター34を加算器及びシフターだけ
で具現でき、システムの構造の単純化が可能である。一
方、送信信号の中心周波数が変更される度に、その変更
された周波数に対して新しく計算された係数をフィルタ
ー係数として使用しなければならない。本発明におい
て、同一の送信信号に対して同一の相関係数を使用する
ためには、送信信号の中心周波数が変更されるとき、サ
ンプリング周波数を変更して送信パターンメモリ13か
らパルサー12へ印加される送信パターンを変更する。
即ち、中心周波数とサンプリング周波数との比を一定に
して、このサンプリング周波数に基づいてRFメモリ3
3からFIRフィルター34への受信信号の供給速度を
変更する。
【0051】本発明の他の実施例では、連続的な重みつ
きチャープ信号の送信の代わりに、前述の発生的アルゴ
リズムを使用して量子化された信号を送信することもで
きる。この場合、パルサ12は高価の線形増幅器の代わ
りに安価のデジタル−アナログ変換器で具現されること
ができる。
【0052】図9〜図13において、原チャープ信号の
中心周波数、サンプリング周波数及び帯域幅は各々5M
Hz、40MHz及び3MHzである。また、この原チ
ャープ信号は8ビットで量子化され、フィルター係数は
ハンニングウィンドウによって変更される。図9〜図1
3中で、Lはフィルター係数の個数を表す。
【0053】図9は、ハンニングウィンドウ関数によっ
て重みつきされたチャープ信号の基本周波数成分に対す
るFIRフィルターの出力を示したグラフである。図9
から分かるように、メインローブの幅は変化しなく、サ
イドローブは医療用超音波撮像システムに適合な水準で
ある。従って、本発明では、乗算器が除去されたFIR
フィルター34はシステムのハードウェア複雑さをより
一層減らし、望ましいシステム性能が提供できることが
分かる。
【0054】図10は、ハンニングウィンドウ関数によ
って重みつきされたチャープ信号の第2高調波成分に対
するフィルターの出力を示したグラフである。図10の
フィルター出力は図9のフィルター出力に比べて、メイ
ンローブの幅が狭くサイドローブのレベルも低い。従っ
て、場合によって高調波成分を考慮してフィルター係数
を選択すればより良い超音波映像を形成できる。
【0055】図11は、FIRフィルター34に供給さ
れる基本周波数及び第2高調波成分を有する受信信号の
例を示すグラフである。高調波成分は例えば、超音波の
非線形的伝搬等によって生じる。このような高調波成分
を用いて、目標オブジェクトによってより良い画質の超
音波映像を得ることができる。
【0056】図12は、図11の受信信号の基本周波数
成分に対するフィルターの出力を示すグラフである。図
12から、受信信号の基本周波数成分によって得られた
メインローブの幅及びサイドローブのレベルが医療用撮
像システムに適合であることが分かる。図12中でグラ
フの右側に示されたサイドローブは第2高調波成分によ
り発生したものであるが、このサイドローブは−60d
B以下であってシステムの特性を大きく低下させない。
【0057】図13は、図11の受信信号の第2高調波
成分に対するフィルターの出力を示すグラフである。図
13から、受信信号の基本周波数成分によって得られた
メインローブの幅及びサイドローブのレベルが医療用撮
像システムに適合であることが分かる。図13中でグラ
フの左側に示されたサイドローブは基本周波数成分によ
って生じたものであるが、このサイドローブは−60d
B以下であってシステムの特性を大きく低下させない。
【0058】上記において、本発明の好適な実施の形態
について説明したが、本発明の請求範囲を逸脱すること
なく、当業者は種々の改変をなし得るであろう。
【0059】
【発明の効果】従って、本発明によれば、重みつきスペ
クトル拡散信号を利用したパルス圧縮システムはメイン
ローブの幅及び最大サイドローブの側面で医療用超音波
撮像システムに適合である。特に、スペクトル拡散信号
をウィンドウ関数で重みつき処理し、受信集束以前にパ
ルス圧縮を行うことによって、受信信号のサイドローブ
を低下させて超音波映像の画質を向上させることができ
る。さらに、本発明によれば、発生的アルゴリズムを利
用して乗算器が除去されたフィルターを設計することに
よって、システムのハードウェア構造を単純化すること
ができる。
【図面の簡単な説明】
【図1】 従来の短パルス超音波撮像システムのブロッ
ク図である。
【図2】 本発明による超音波撮像システムの概略的な
ブロック図である。
【図3】 従来のチャープ信号を示したグラフである。
【図4】 本発明で使用されるハミングウィンドウ、ハ
ンニングウィンドウ、ブラックマンウィンドウ及び矩形
ウィンドウを示す模式図である。
【図5】 ハンニングウィンドウによって重み付きされ
たチャープ信号を示すグラフである。
【図6】 従来の整合フィルターの形態で具現されたF
IRフィルターの構造を示した模式図である。
【図7】 本発明によるFIRフィルターの構造を示し
た模式図である。
【図8】 本発明の実施例によってフィルター係数の決
定に用いられる発生的アルゴリズムを示すフロー図であ
る。
【図9】 ハンニングウィンドウ関数によって重みつき
されたチャープ信号の基本周波数成分に対するFIRフ
ィルターの出力を示したグラフである。
【図10】 ハンニングウィンドウ関数によって重みつ
きされたチャープ信号の第2高調波成分に対するフィル
ターの出力を示したグラフである。
【図11】 FIRフィルターに供給される基本周波数
及び第2高調波成分を有する受信信号の例を示すグラフ
である。
【図12】 図11の受信信号の基本周波数成分に対す
るフィルターの出力を示すグラフである。
【図13】 図11の受信信号の第2高調波成分に対す
るフィルターの出力を示すグラフである。
【符号の説明】
12 パルサー 13 Txパターンメモリ 14 Tx集束遅延メモリ 15 フィルター係数制御器 16 フィルター係数メモリ 21 Tx/Rxスイッチ 31 レシーバ 34 FIRフィルター 35 パルス圧縮器 36 Rx集束遅延調整器 37 ビーム形成器 41 信号処理器 42 スキャンコンバータ 200 超音波撮像システム
───────────────────────────────────────────────────── フロントページの続き (72)発明者 ユウ ヤン モ 大韓民国 ソウル トクビョルシ ヨンド ンポクダンサンドン5ガ カンナムアパー ト25ドン205ホ Fターム(参考) 2G047 BA03 BC13 CA01 DB02 EA07 EA14 GB02 GF07 GF12 GF18 GG17 4C301 AA01 BB22 EE07 EE11 EE15 GB02 HH01 HH24 HH43 JB02 JB29 JB33 JB35 5J083 AA02 AB17 AC15 AC28 AD13 AE10 BA03 BC01 BE53 BE57 BE58

Claims (16)

    【特許請求の範囲】
  1. 【請求項1】 オブジェクトの映像を撮像する超音波撮
    像方法であって、 予め決められた第1スペクトル拡散信号を超音波信号に
    変換して前記オブジェクトに発射する第1ステップと、 前記オブジェクトから反射された反射信号を受け取っ
    て、前記反射信号に対してパルス圧縮を行ってパルス圧
    縮信号を生成する第2ステップと、 前記パルス圧縮信号を処理して受信集束信号を生成し、
    前記受信集束信号から前記オブジェクトの映像を生成す
    る第3ステップとを含むことを特徴とする超音波撮像方
    法。
  2. 【請求項2】 前記予め決められた第1スペクトル拡散
    信号が、少なくとも一つ以上のトランスジューサで前記
    超音波信号に変換されることを特徴とする請求項1記載
    の超音波撮像方法。
  3. 【請求項3】 前記予め決められた第1スペクトル拡散
    信号が、チャープ信号に対して第1ウィンドウ関数を適
    用することによって生成されることを特徴とする請求項
    1記載の超音波撮像方法。
  4. 【請求項4】 前記第1ウィンドウ関数が、ハンニング
    ウィンドウ、ハミングウィンドウ、ブラックマンウィン
    ドウのうちのいずれか一つであることを特徴とする請求
    項3記載の超音波撮像方法。
  5. 【請求項5】 前記第2ステップが、前記反射信号を予
    め決められたフィルター係数でフィルタリングするステ
    ップを含むことを特徴とする請求項1記載の超音波撮像
    方法。
  6. 【請求項6】 前記予め決められた第1スペクトル拡散
    信号が、予め決められた第2スペクトル拡散信号に対し
    て第2ウィンドウ関数を適用して生成される重みつきス
    ペクトル拡散信号であり、前記予め決められたフィルタ
    ー係数が、前記第2スペクトル拡散信号の有限インパル
    ス応答(FIR)フィルター係数を変更し、前記変更され
    たFIRフィルター係数に対して発生的アルゴリズムを
    適用して、該係数が2の自乗またはその和で表現される
    よう量子化することによって生成されることを特徴とす
    る請求項5記載の超音波撮像方法。
  7. 【請求項7】 前記FIRフィルター係数がウィンドウ
    関数の適用によって変更されることを特徴とする請求項
    6に記載の超音波撮像方法。
  8. 【請求項8】 前記第1スペクトル拡散信号が、発生的
    アルゴリズムによって求められる量子化信号であること
    を特徴とする請求項1記載の超音波撮像方法。
  9. 【請求項9】 前記パルス圧縮信号が、前記反射信号の
    基本周波数成分信号及び高調波成分信号のうちのいずれ
    か一つに対して提供されることを特徴とする請求項1記
    載の超音波撮像方法。
  10. 【請求項10】 超音波撮像システム及びオブジェクト
    撮像方法で用いられるフィルター係数を決定する方法で
    あって、 予め決められたスペクトル拡散信号のFIRフィルター
    係数を求める第1ステップと、 前記FIRフィルター係数に対して発生的アルゴリズム
    を適用して、該係数が2の自乗またはその和で表現され
    るよう量子化する第2ステップとを含むことを特徴とす
    るフィルター係数決定方法。
  11. 【請求項11】 オブジェクトの映像を撮像する超音波
    撮像システムであって、 予め決められたスペクトル拡散信号に応じて、前記オブ
    ジェクトに超音波信号を発射する一つまたはそれ以上の
    トランスジューサと、 前記オブジェクトから反射される反射信号を受け取っ
    て、前記反射信号に対してパルス圧縮を行ってパルス圧
    縮信号を生成するパルス圧縮手段と、 前記パルス圧縮信号を処理して受信集束信号を生成し、
    前記受信集束信号に基づいて前記オブジェクトの映像を
    生成する手段とを含むことを特徴とする超音波撮像シス
    テム。
  12. 【請求項12】 前記スペクトル拡散信号が、チャープ
    信号に対してウィンドウ関数を適用することで生成され
    る重みつきチャープ信号であることを特徴とする請求項
    11記載の超音波撮像システム。
  13. 【請求項13】 前記パルス圧縮手段が、予め決められ
    たフィルター係数で前記反射信号をフィルタリングする
    手段を含むことを特徴とする請求項11記載の超音波撮
    像システム。
  14. 【請求項14】 前記スペクトル拡散信号が、発生的ア
    ルゴリズムによって求められる量子化信号であることを
    特徴とする請求項11記載の超音波撮像システム。
  15. 【請求項15】 前記パルス圧縮信号が、前記反射信号
    の基本周波数成分信号及び高調波成分信号のうちのいず
    れか一つに対して提供されることを特徴とする請求項1
    1記載の超音波撮像システム。
  16. 【請求項16】 前記パルス圧縮手段が、フィルター係
    数を用いて前記反射信号をフィルタリングするフィルタ
    リング手段を含み、 前記フィルター係数が、前記スペクトル拡散信号のFI
    Rフィルター係数を計算し、前記FIRフィルター係数
    に対して発生的アルゴリズムを適用して該係数が2の自
    乗またはその和で表現されるよう量子化することによっ
    て決められ、 前記フィルタリング手段が加算器及びシフターを用いて
    構成されることを特徴とする請求項11記載の超音波撮
    像システム。
JP2001180172A 2000-06-17 2001-06-14 フィルター係数決定方法、超音波撮像システム及びその方法 Pending JP2002045359A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
KR1020000033417A KR100350026B1 (ko) 2000-06-17 2000-06-17 확산 대역 신호를 이용한 펄스 압축 방식에 기초한 초음파영상 형성 방법 및 장치
KR2000-33417 2000-06-17

Publications (1)

Publication Number Publication Date
JP2002045359A true JP2002045359A (ja) 2002-02-12

Family

ID=36314059

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2001180172A Pending JP2002045359A (ja) 2000-06-17 2001-06-14 フィルター係数決定方法、超音波撮像システム及びその方法

Country Status (5)

Country Link
US (1) US6622560B2 (ja)
EP (1) EP1164385B1 (ja)
JP (1) JP2002045359A (ja)
KR (1) KR100350026B1 (ja)
DE (1) DE60118837T2 (ja)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007508544A (ja) * 2003-10-01 2007-04-05 フローライン インコーポレイテッド 深さ測定システム
JP2008511408A (ja) * 2004-08-27 2008-04-17 ベラソン インコーポレイテッド 超音波画像内の体腔液を定量化および分類化するためのシステムおよび方法
JP2009156694A (ja) * 2007-12-26 2009-07-16 Kanazawa Inst Of Technology 計測方法および装置
JP2010142442A (ja) * 2008-12-19 2010-07-01 Konica Minolta Medical & Graphic Inc 超音波画像診断装置
JP2012502692A (ja) * 2008-09-15 2012-02-02 テラテク・コーポレーシヨン 3次元超音波画像形成システム
KR101643304B1 (ko) * 2015-08-20 2016-07-28 대진대학교 산학협력단 비선형 필터를 이용한 초음파 영상의 부엽 저감장치

Families Citing this family (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100406097B1 (ko) * 2001-12-26 2003-11-14 주식회사 메디슨 가중 쳐프 신호를 이용한 초음파 영상 형성 장치 및 방법
US7819806B2 (en) * 2002-06-07 2010-10-26 Verathon Inc. System and method to identify and measure organ wall boundaries
GB2391625A (en) 2002-08-09 2004-02-11 Diagnostic Ultrasound Europ B Instantaneous ultrasonic echo measurement of bladder urine volume with a limited number of ultrasound beams
US20090112089A1 (en) * 2007-10-27 2009-04-30 Bill Barnard System and method for measuring bladder wall thickness and presenting a bladder virtual image
WO2006026605A2 (en) * 2002-06-07 2006-03-09 Diagnostic Ultrasound Corporation Systems and methods for quantification and classification of fluids in human cavities in ultrasound images
US8221322B2 (en) * 2002-06-07 2012-07-17 Verathon Inc. Systems and methods to improve clarity in ultrasound images
US7520857B2 (en) * 2002-06-07 2009-04-21 Verathon Inc. 3D ultrasound-based instrument for non-invasive measurement of amniotic fluid volume
US20100036252A1 (en) * 2002-06-07 2010-02-11 Vikram Chalana Ultrasound system and method for measuring bladder wall thickness and mass
US20060025689A1 (en) * 2002-06-07 2006-02-02 Vikram Chalana System and method to measure cardiac ejection fraction
US20090062644A1 (en) * 2002-06-07 2009-03-05 Mcmorrow Gerald System and method for ultrasound harmonic imaging
US20080262356A1 (en) * 2002-06-07 2008-10-23 Vikram Chalana Systems and methods for ultrasound imaging using an inertial reference unit
US20040127797A1 (en) * 2002-06-07 2004-07-01 Bill Barnard System and method for measuring bladder wall thickness and presenting a bladder virtual image
US8221321B2 (en) * 2002-06-07 2012-07-17 Verathon Inc. Systems and methods for quantification and classification of fluids in human cavities in ultrasound images
US7094204B2 (en) * 2002-08-23 2006-08-22 Siemens Medical Solutions Usa, Inc. Coded excitation imaging for use with bipolar, unipolar and other waveforms
EP1532433A2 (en) * 2002-08-28 2005-05-25 Wayne State University System for infrared imaging by inducing acoustic chaos
KR100579715B1 (ko) * 2002-09-16 2006-05-15 학교법인 인하학원 유전자 알고리즘을 이용한 영상의 적응적인 전처리 기법
US20050050483A1 (en) * 2003-08-25 2005-03-03 Keller S. Brandon System and method analyzing design elements in computer aided design tools
US7151484B2 (en) 2003-09-30 2006-12-19 Kabushiki Kaisha Toshiba Pulse compression processor
US20060160516A1 (en) * 2005-01-14 2006-07-20 Dongsong Zeng Pulse shaping optimizer in UWB receiver
JP4163733B2 (ja) * 2006-07-18 2008-10-08 アロカ株式会社 超音波診断装置
JP5283888B2 (ja) * 2006-11-02 2013-09-04 株式会社東芝 超音波診断装置
US8167803B2 (en) * 2007-05-16 2012-05-01 Verathon Inc. System and method for bladder detection using harmonic imaging
US8225998B2 (en) * 2008-07-11 2012-07-24 Es&S Innovations Llc Secure ballot box
CA2957778C (en) * 2008-08-07 2022-12-13 Jongtae Yuk Device, system, and method to measure abdominal aortic aneurysm diameter
KR101126182B1 (ko) * 2010-04-01 2012-03-22 서강대학교산학협력단 비선형 쳐프 신호를 이용한 고조파 영상 장치 및 방법
WO2012160541A2 (en) 2011-05-25 2012-11-29 Orcasonix Ltd. Ultrasound imaging system and method
KR101336757B1 (ko) * 2012-05-21 2013-12-04 동국대학교 산학협력단 주파수 변조된 초음파를 이용한 세포 또는 입자 제어 장치 및 방법
GB201222882D0 (en) 2012-12-19 2013-01-30 Univ Leeds Ultrasound generation
SG11201509789WA (en) * 2013-05-31 2015-12-30 Nestec Sa Systems and methods for detecting water/product interfaces during food processing
KR102146374B1 (ko) * 2013-11-18 2020-08-20 삼성전자주식회사 초음파 영상장치 및 그 제어방법
CN112162266B (zh) * 2020-09-28 2022-07-22 中国电子科技集团公司第五十四研究所 一种基于凸优化理论的共形阵二维波束优化方法

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NO831718L (no) * 1983-05-13 1984-11-14 Vingmed As Fremgangsmaate og apparat ved blodstroem-hastighetsmaaling med ultralyd for dannelse av todimensjonal avbildning av blodets hastighet
US5345426A (en) 1993-05-12 1994-09-06 Hewlett-Packard Company Delay interpolator for digital phased array ultrasound beamformers
JPH095310A (ja) 1995-06-16 1997-01-10 Nkk Corp 超音波検査方法
JP3482039B2 (ja) 1995-07-31 2003-12-22 アルパイン株式会社 ディジタルフィルタの係数決定方法
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
US5938611A (en) * 1998-03-26 1999-08-17 General Electric Company Method and apparatus for color flow imaging using coded excitation with single codes
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
US6074348A (en) 1998-03-31 2000-06-13 General Electric Company Method and apparatus for enhanced flow imaging in B-mode ultrasound
US6010456A (en) * 1998-12-30 2000-01-04 General Electric Company Method and apparatus for acoustic subtraction imaging using linear and nonlinear ultrasonic images
US6155980A (en) * 1999-03-16 2000-12-05 General Electric Company Ultrasonic imaging system with beamforming using unipolar or bipolar coded excitation
US6343510B1 (en) * 1999-04-22 2002-02-05 Vn Instruments Limited Ultrasonic testing using synthetic impulses
KR100362000B1 (ko) * 2000-02-01 2002-11-22 주식회사 메디슨 변형된 골레이 코드를 이용한 펄스 압축 방식의 초음파영상 형성 방법 및 장치
US6309356B1 (en) * 2000-03-06 2001-10-30 Acuson Corporation Method and apparatus for forming medical ultrasound images

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007508544A (ja) * 2003-10-01 2007-04-05 フローライン インコーポレイテッド 深さ測定システム
JP2008511408A (ja) * 2004-08-27 2008-04-17 ベラソン インコーポレイテッド 超音波画像内の体腔液を定量化および分類化するためのシステムおよび方法
JP2012101103A (ja) * 2004-08-27 2012-05-31 Verathon Inc 超音波画像内の体腔液を定量化および分類化するためのシステムおよび方法
JP2009156694A (ja) * 2007-12-26 2009-07-16 Kanazawa Inst Of Technology 計測方法および装置
JP2012502692A (ja) * 2008-09-15 2012-02-02 テラテク・コーポレーシヨン 3次元超音波画像形成システム
JP2010142442A (ja) * 2008-12-19 2010-07-01 Konica Minolta Medical & Graphic Inc 超音波画像診断装置
KR101643304B1 (ko) * 2015-08-20 2016-07-28 대진대학교 산학협력단 비선형 필터를 이용한 초음파 영상의 부엽 저감장치

Also Published As

Publication number Publication date
KR20010113215A (ko) 2001-12-28
DE60118837T2 (de) 2007-05-16
EP1164385A3 (en) 2004-06-02
KR100350026B1 (ko) 2002-08-24
DE60118837D1 (de) 2006-05-24
US6622560B2 (en) 2003-09-23
EP1164385A2 (en) 2001-12-19
EP1164385B1 (en) 2006-04-19
US20020005071A1 (en) 2002-01-17

Similar Documents

Publication Publication Date Title
JP2002045359A (ja) フィルター係数決定方法、超音波撮像システム及びその方法
JP3423935B2 (ja) 超音波撮像方法及び改良ゴールレイコードを用いたパルス圧縮方式に基づく装置
EP1122556B1 (en) Enhanced tissue-generated harmonic imaging using coded excitation
US6547733B2 (en) Ultrasound imaging apparatus and method using Golay codes with orthogonal property
EP0985936B1 (en) Nonlinear imaging using orthogonal transmit and receive codes
US7393324B2 (en) Ultrasonic transmitting and receiving apparatus
EP1121901A2 (en) Angle independent ultrasound volume flow measurement
US6050947A (en) Method and apparatus for harmonic tissue imaging and contrast imaging using coded transmission
US20040034305A1 (en) Ultrasound imaging system and method based on simultaneous multiple transmit-focusing using weighted orthogonal chirp signals
JP4227783B2 (ja) 超音波撮像システム
JP2002534184A (ja) コヒーレンス・イメージングのための方法および装置
JPH11197151A (ja) 超音波イメージング・システム用のbモード処理装置および検出後像処理方法
JP2000157548A (ja) 超音波散乱体をイメ―ジングするための方法及びシステム
KR20010014579A (ko) 유니폴라 또는 바이폴라 코딩 여기를 이용하여 빔을형성하는 초음파 이미징 시스템
US6786097B2 (en) Ultrasound imaging system and method using weighted chirp signals
JP4469031B2 (ja) 超音波散乱体をイメージングするシステム及び方法
JP4490128B2 (ja) 超音波診断装置
JP3222648B2 (ja) 医用超音波分散圧縮送受信方法及び医用超音波分散圧縮送受信装置
JP4499477B2 (ja) 超音波診断装置
CN113768542A (zh) 超声血流成像装置及超声设备
JPH08131435A (ja) 超音波診断装置
JPH0723959A (ja) 超音波診断装置
WO2004055941A1 (fr) Dispositif de production de retards pour appareil de synthese de faisceaux
JPH01303133A (ja) 分散圧縮方式パルスエコーシステム