JP5905080B2 - オーバーラップする送信ビームにおける適格と評価された領域を使用する増強された超音波画像形成 - Google Patents

オーバーラップする送信ビームにおける適格と評価された領域を使用する増強された超音波画像形成 Download PDF

Info

Publication number
JP5905080B2
JP5905080B2 JP2014509393A JP2014509393A JP5905080B2 JP 5905080 B2 JP5905080 B2 JP 5905080B2 JP 2014509393 A JP2014509393 A JP 2014509393A JP 2014509393 A JP2014509393 A JP 2014509393A JP 5905080 B2 JP5905080 B2 JP 5905080B2
Authority
JP
Japan
Prior art keywords
transmit
field
beams
field point
transmit beam
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.)
Active
Application number
JP2014509393A
Other languages
English (en)
Other versions
JP2014512923A (ja
Inventor
デイグル,ロナルド,エルヴィン
Original Assignee
ヴェラゾニックス,インコーポレーテッド
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 ヴェラゾニックス,インコーポレーテッド filed Critical ヴェラゾニックス,インコーポレーテッド
Publication of JP2014512923A publication Critical patent/JP2014512923A/ja
Application granted granted Critical
Publication of JP5905080B2 publication Critical patent/JP5905080B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52019Details of transmitters
    • G01S7/5202Details of transmitters for pulse systems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4483Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
    • A61B8/4488Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer the transducer being a phased array
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • 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
    • 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/8909Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
    • G01S15/8915Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
    • G01S15/8927Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array using simultaneously or sequentially two or more subarrays or subapertures
    • 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/8977Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using special techniques for image reconstruction, e.g. FFT, geometrical transformations, spatial deconvolution, time deconvolution
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52023Details of receivers
    • G01S7/52033Gain control of receivers
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52046Techniques for image enhancement involving transmitter or receiver
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52085Details related to the ultrasound signal acquisition, e.g. scan sequences
    • 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/8979Combined Doppler and pulse-echo imaging systems
    • 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/8979Combined Doppler and pulse-echo imaging systems
    • G01S15/8984Measuring the velocity vector

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Veterinary Medicine (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Acoustics & Sound (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Gynecology & Obstetrics (AREA)
  • Theoretical Computer Science (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Description

関連出願の相互参照
本出願は、2011年5月2日に出願した米国特許仮出願第61/481,476号の利益を、米国特許法第119条(e)に従って主張するものであり、この仮出願は、参照により全体として本明細書に組み込まれる。
本開示は、超音波撮像に関するものであり、より詳細には、画像形成処理のために画像フィールド点特性(image field point characteristics)を決定し利用するための方法及びシステムに関するものである。
超音波撮像では、高周波音波を利用して、物体の内側を可視化する。例えば、超音波探傷検査機器などの音響検知デバイスは、人体、溶接部の周りの領域、及び木材ベースのパネルなどの製造品を含む、さまざまな物体の内側を検査する際に使用される。本開示は、医療用超音波検査などの組織撮像の文脈において説明されているが、応用はこの分野以外にもある。
医療用超音波撮像は、さまざまな疾病状態及び病状を診断するための効果的なツールに発展している。音波機器の市場は、画質及びさまざまな種類の組織を区別する機能の改善に支えられて、年月を重ね着実に成長してきている。超音波撮像は、常に、広範な信号及び画像処理方法を必要とし、特に、それぞれが固有の信号処理要件を有する128個以上ものトランスデューサ素子を採用するアレイシステムにはそのような方法を必要としてきた。この10年間に、市場の最下層のものを除く、ほとんどすべてのシステムにおけるデジタル信号処理の精度及び柔軟性の改善への移行が見られた。この移行は、送信された音波及び戻った超音波エコー信号における情報の多くを利用することができる画像形成の方法が改善される可能性をもたらした。
商業用の超音波システムは、典型的には、画像形成のために集束送信ビームを利用する。二次元画像フィールドは、典型的には、最良の画像分解能が望まれているフィールド内のある深さにそれぞれが集束される、フィールドの幅にわたって均一な間隔で並ぶ送信ビームの集合による高周波音波を当てられる。それぞれの順次送信ビームから戻る超音波エコーが受信され、処理されて、さらに画像データのより多くのラインのうちの一つが得られ、それらのラインはこの集合内のそれぞれの送信ビームの軸に対応する。次いで、複数の画像データラインが一つのピクセルアレイに補間され、一つの画像を生成する。
前述の方法は、図1に例示されている。典型的な順次ライン走査30は、128本の送信ビーム32を含むこともあり、その結果、画像ライン34は128となり、これらは次いで表示できるようにピクセル36に補間される。一般に、送信ビーム32は、特に送信ビームの焦点深度とは異なる深さで、再構成された画像ライン34に比べてかなり広くなる。この結果、横方向分解能は深さとともに変化し、図1に示されているように送信集束ゾーンで最もシャープな分解能となる。ビームの軸にそったフィールド点のみが、画像パラメータの再構成に使用され、その結果、画像ライン34の集合は走査で使用される送信ビーム32の数に等しくなる。この画像ライン34の集合は、一般的に、ディスプレイ38内のピクセル36の間隔に関してまばらな間隔で並ぶので、画像ライン34内の点は、ディスプレイ内のそれぞれのピクセル36について補間されなければならない。
現代的な商用システムでは、フィールドの幅にわたってそれぞれの順次走査位置での複数の送信を利用することによってより大きな被写界深度(depth of field)上で横方向分解能を改善することを試みている。それぞれの位置において、複数の送信ビームは、注目する深さの上で異なる集束ゾーンを利用する。それぞれのゾーンからの画像ラインデータが組み合わされて、より大きな実行被写界深度が得られる。この技術は、収集時間の増加又はフレームレートの低下と引き換えに横方向分解能を改善する。完全な画像フレームを生成する時間は、個別の送信ビームのそれぞれからエコー信号を収集するのに要する時間の総和であり、これは一般的に音速と撮像される媒体の注目する最大の深さとによって制限される。それぞれの画像フレームに対して利用される送信ビームの数が多ければ多いほど、画像フレームを収集する時間が長くなり、フレームレートは低速になる。
フレームレートは、多くの超音波用途において重要なファクタであるため、画像収集時間を短縮するために別の技術が利用されることが多い。送信ビームは、それぞれのビームにおいて複数の受信ラインが再構成されうるように弱集束され、したがって、並ぶ送信ビームの間隔をフィールドの幅にわたって広げることができ、また画像フィールドを被覆するのに必要な送信ビームの総数を減らせる。この技術は、画像収集をスピードアップし、複数の集束ゾーンの使用に対して失われるフレームレートの一部を回復することができるが、横方向分解能は、典型的には、それぞれの送信ビームが広がるため低下する。
本開示の一態様によれば、多素子トランスデューサから得られた超音波画像の分解能及びフレームレートの両方を改善するための方法及びシステムが提供される。
この方法は、
a.画像フレームを収集するために使用される送信ビームの集合内のそれぞれの送信ビームに対する送信開口及び個別の素子送信特性を指定するステップであって、前記送信ビームの集合は、所望の視野の少なくとも一部に高周波音波を当てる少なくとも一つの送信ビームを含む、ステップと、
b.測定又はシミュレーションを通じて、視野を覆う複数のフィールド点で前記送信ビームの集合内のそれぞれの送信ビームのさまざまな属性を決定するステップと、
c.それぞれの新しい画像フレームを処理する際の繰り返し使用に対するそれぞれのフィールド点の前記属性を格納するステップと、
d.前記送信ビームの集合を使用してエコー信号を送受信し、前記エコー信号をメモリに格納するステップと、
e.前記格納されている送信ビーム属性のうちの1つ以上を使用して、前記格納されているエコー信号を処理し、画像形成で使用するため与えられた送信ビームが寄与するフィールド点から受信された超音波エコー信号の適格性評価(qualify)を行い、信号がどのように処理されるべきかを指定するステップと、
f.フィールド点に対してオーバーラップする、適格と評価された送信ビーム(qualified transmit beam)の集合から複数の画像パラメータを生成して組み合わせて、前記画像フレームに対するフィールド点の値を構成する最終画像パラメータを生成するステップとを含む。
本開示のさらなる態様により、送信ビームの測定された、又はシミュレートされた属性は、
a.送信イベントの持続時間にわたる送信ビームのピーク強度、
b.送信ビームのピーク強度がフィールド点で生じた時間、
c.パルス強度がパルスピークの典型的には−20dB未満である閾値を超える時間によって決定される送信パルス持続時間、
d.前記ピーク強度と、フィールド点で送信イベントの発生中に生じるビーム強度の他のピークとの比、
e.フィールド点の位置を指定するために使用される選択された軸を持つ送信ビームパルスの波面の入射角、
のうちの少なくとも一つからなる。
本開示の一実施形態において、送信ビームの集合は、集束しないビームを使用して送受信収集の少ない視野に高周波音波を当てることによってフレームレートを最大化するように設計される。別の実施形態では、送信ビームの集合は、空間的及びコントラスト分解能を最大化するように設計され、その一方で、従来の走査システムの場合と比べて大きなレートでフレームをそのまま収集する。
本開示の別の態様によれば、組織又は血流速度に対する受信エコー信号の処理では、入射送信ビーム属性の格納されている角度を使用して、ドップラー周波数偏移に基づき速度計算を補正し、それぞれのフィールド点における運動の大きさと運動の方向を求める。
本開示の前述の特徴及び利点並びに他の特徴及び利点は、添付の図面に関して、以下の詳細な説明を読むことでよりよく理解できるため明白であろう。
集束ビームを使用する従来の走査及び画像形成の図である。 カットオフ閾値によって決定される境界を示すフィールド点のグリッド上の送信ビームの最大強度の図である。 異なるフィールド点における送信波形の図である。 三つ又は四つの受信されたエコー信号が画像形成のため組み合わされうるオーバーラップした領域内のフィールド点を示す四つのオーバーラップする送信ビームの図である。 ピクセル指向処理の一実施形態のソフトウェアベースのアーキテクチャの概略図である。 ピクセル指向処理に従って形成されたプラグインモジュールの線図である。 ピクセル指向処理に従って形成された128素子線形配列に対する収集データの概略図である。 ピクセル指向処理で使用されるピクセルマッピングプロセスの図である。 本開示のプロセスに対するシステムアーキテクチャの高水準表現の図である。
以下の説明では、さまざまな開示されている実施形態を完全に理解できるように、いくつかの具体的詳細が述べられている。しかし、当業者であれば、実施形態は、これらの具体的詳細の1つ以上がなくても、又は他の方法、コンポーネント、材料などを使用しても、実施できることを理解するであろう。
文脈上他の意味に解すべき場合を除き、本明細書及び添付の特許請求の範囲全体を通して、「含む」、「備える」、及びこれらの活用形は、開いた包含的意味で解釈されるべき、つまり、「限定はしないが・・・を含む」という意味で解釈すべきである。
本明細書全体を通して「一実施形態」又は「実施形態」と記述されている場合、これは、その実施形態に関して説明されている特定の特徴、構造、又は特性が少なくとも一つの実施形態に含まれることを意味する。そのため、「一実施形態では」又は「実施形態では」という語句が本明細書全体のさまざまな箇所に記載されていても、必ずしもすべて同じ実施形態を指しているとは限らない。さらに、特定の特徴、構造、又は特性は、1以上の実施形態において適当な任意の方法で組み合わせられうる。
本明細書及び添付の特許請求の範囲で使用されているように、「一つの(又は使わない場合もある)」及び「その(使わない場合もある)」(英語原文の冠詞「a」、「an」、及び「the」に対応する)で示される単数形は、文脈上明らかにそうでないことを示していない限り、複数形の指示対象を含む。「又は」という言い回しは、一般的に、内容が明らかに別のことを示していない限り最も広い意味で、つまり「及び/又は」を意味するものとして使用されることにも留意されたい。
本明細書の見出し及び「要約」は、便宜上のものにすぎず、実施形態の範囲又は意味を解釈することはしない。
媒体に高周波音波を当てるために、深さで集束又は弱集束(weakly focused)されるか、又は完全に集束されない送信ビームを含む、さまざまなビーム特性を持つ部分的に又は完全にオーバーラップする送信ビームの集合を使用することができる超音波画像形成の方法及び対応するシステムが説明される。本開示の一態様によれば、方法は、送信ビームの関連する特性を予めシミュレーション又は測定を通じてあらかじめ決定することを伴い、ここで、送信ビームは、この集合内の与えられたビームに対する送信パルスがフィールド点を通過するときにそれぞれの画像フィールド点において生成される。
このような特性として、フィールド点において結果として生じるピーク音響強度、フィールド点においてピーク音響強度が生じる時間、特定のレベルより大きい強度の持続時間(パルス持続時間)、ピーク強度とフィールド点において送信イベントの発生時に生じるビーム強度の他のピークとの間の比、及びフィールド点とのピーク送信強度波の入射角が挙げられる。次いで、画像形成に適している送信フィールドの領域の適格性評価(qualify)を行うためにこれらの特性の1以上が使用されうる。これらの送信フィールド特性は、コンピューティングシステムのメモリに格納され、次いで、すでに特徴付けられている送信ビーム集合を使用する超音波走査のために画像形成処理時に取り出されうる。どのフィールド点が送信ビームの適格と評価される領域(qualified region)内に入るか、またそれぞれのフィールド点において与えられた送信ビームについて戻されるエコー信号がどのように処理されるべきかを決定するためにビーム特性が使用されうる。フィールド点における最終音響画像パラメータは、集合内の送信ビームの1以上を基に処理済み信号の組み合わせから導き出される。
本開示では、それぞれの送信ビームから生成される複数のエコー信号を利用し、画像形成プロセスで複数のオーバーラップする送信ビームからの信号を組み合わせて画像分解能を改善し、及び/又はフレーム収集時間を短縮する。図1に示されているように、個別の送信ビーム32は、トランスデューサの視野のかなり大きな領域に高周波音波を当てることができる。さらに、典型的な順次走査で高周波音波を当てられた領域はかなりオーバーラップし、同じフィールド点がこの集合内の複数の送信ビームによる高周波音波で当てられる。それぞれのビームの軸にそってのみ画像パラメータを再構成する代わり、送信ビームに含まれる他のフィールド点を処理し、フィールド点におけるそれぞれのオーバーラップするビームからの結果を組み合わせることが可能である。このような再構成法を実行するために、われわれは、送信ビームの高周波音波を当てられる領域を決定し、領域内のフィールド点に対する最良の種類の処理を決定するのを補助する領域のパラメータを測定する何らかの方法を必要としている。
送信ビームによる高周波音波を当てられる領域の範囲を決定する一方法では、図2に示されているように、グリッド42内のそれぞれのフィールド点40における送信ビームの最大強度を、ビームシミュレーション技術を使用して、又は直接測定によって、計算する。媒体中の音速に依存する、ある知られている時刻において、パルス送信ビームは、ビームの高周波音波を当てられた領域内のフィールド点で強度ピークを発生する。この強度ピークは、ビームの軸に関してフィールド点の位置とともに変化し、典型的には、ビーム軸からの距離が長くなるにつれ減少する。フィールド点における強度ピークの大きさは、図2のフィールド点の灰色の陰によって示される。ビーム軸からある距離のところで、ピーク強度は、フィールド点から戻るエコー信号が、バックグラウンドノイズ及び音波散乱と混じるので、検知不可能又は使用不可能になるくらい弱くなる。したがって、カットオフ閾値は、最大ピーク強度のある比で設定することができ、これに対して、その閾値より低い強度を持つフィールド点は送信ビームの高周波音波を当てられる領域の外部にあると指定される。異なるフィールド点で異なるカットオフ閾値を有するように選択することができるが、それは、フィールド点における送信パルスの全強度が深さ、組織減衰、及び他のファクタとともに変化しうるからである。次いで、さまざまな閾値が、ビームの主要な高周波音波を当てられる領域の境界43を画成するが、これは、適格と評価される領域(qualified region)とここでは称している(図2を参照)。送信ビームの高周波音波を当てられる領域内のフィールド点に対し、画像形成操作でピーク強度値を使用してそれぞれのフィールド点から戻るエコーを正規化することもできる。これにより、送信ビームの主要な高周波音波を当てられる領域の上に画像パラメータの均一な強度再構成を行うことができる。
送信パルスのピーク強度と異なる送信ビーム特性は、画像形成にどのフィールド点を使用できるのかという適格性評価を行うために重要である場合がある。ビームの焦点から移動されたいくつかの距離のところにあるフィールド点に対して、送信パルスの持続時間は延長されうるが、それは、トランスデューサ内の個別の素子からの送信パルスについての到来時間が同時でないからである。焦点のビームの軸上のフィールド点53におけるパルス強度の持続時間を見た場合、焦点から遠い、軸外れのフィールド点54におけるパルス強度の持続時間と比較すると、図3に示されているように、波形50及び波形51であることがわかるであろう。
焦点におけるフィールド点に対する強度波形51は、適切な振る舞いをし、特異ピーク55を形成する。軸外れ点における強度波形50は、より広がっており、二次ピーク52を有する。中央のピークの広がり、及びフィールド点における強度波形の二次ピークの存在により、フィールド点からの画像形成に対するエコー信号の処理において分解能が低下し、散乱が増大する。したがって、強度波形の中央のピークが指定された量を超えて広がっているか、又は中央のピークのある指定された比より高い二次ピークがあるフィールド点を除外すると都合がよい。図3において、二次ピークは、主ピークの比αであり、高い二次ピークを有するフィールド点を除去するために使用されうるαの最大値を設定することができる。二次ピークの低い寄与を保証するのに望ましいαの値は、0.1(−20dB)未満となるであろう。
次いで、送信ビームの上述の特性の一部又は全部を使用してどのフィールド点がその特定のビームによる画像形成に使用されうるかの適格性評価を行うことで、フィールド点について計算された画像パラメータが正確であり、不要信号によって損なわれていないことを確認することができる。したがって、この集合内のそれぞれの送信ビームについて、ビーム軸のみに沿うのではなく、送信ビームの実質的領域を覆うフィールド点にわたって画像形成操作を実行することが可能である。これにより、ビーム軸のみが再構成されるときに必要な多数のビームではなく、少数の部分的に、又は完全にオーバーラップする送信ビームのみで注目する所望の画像フィールドを覆う送信ビームの集合を設計することが可能になる。64本より少ないビームを使用することで、収集時間を短縮し、フレームレートを高めることができる。それに加えて、同じフィールド点は、集合内の複数の送信ビームに対する適格と評価されたフィールド点であってよいので、画像形成処理では、1本だけなくそれ以上の数の送信ビームからのエコー信号を利用することができる。
特定の一例について、図4に示されているように、複数の送信ビームT1−Nの集合によって高周波音波を当てられる与えられたフィールド点P(63)を考える。図4には、四つのオーバーラップするビーム、三つのオーバーラップするビーム、及び四つのオーバーラップするビームを持つフィールド点を含む三つの領域60、61、62が示されている。さらに完全に説明されているように、複数のオーバーラップする送信ビーム領域を使用する画像形成では、空間分解能が高められる。それぞれの送信ビームについて、以下の特性が、シミュレーション又は直接測定によってフィールド点Pにおいて決定される。
I(T)−Pにおけるピーク音響強度。
tp(T)−Pにおける送信開始からピーク音響強度の発生までの時間。
tb(T)−tp(T)の指定された数分の一よりも大きい音響強度の持続時間。
r(T)−二次強度ピークの強度とI(T)との比。
θ(T)−トランスデューサの軸に関する音波の波面の入射角。
エコー画像形成のために、個別のトランスデューサ素子受信信号を遅延させ、総和してフィールド点で反射体によって生成されるエコー信号を再構成することによって、集合内の送信ビーム毎にフィールド点Pに対する複素信号値Snを求めることができる。これは、ビーム形成のよく知られている方法であるけれども、この場合、信号再構成は、それぞれの送信ビームの軸に沿って配置されている点ではなく、特定のフィールド点で実行されている。信号再構成で使用するさまざまなトランスデューサ素子信号に対する時間遅延を計算することを目的として、送信の開始からフィールド点におけるピーク音響強度tp(Tn)の発生までの時間が使用されうる。
この事例では、われわれのフィールド点Pにおいてある程度の音響エネルギーを発生する8本ほどの送信ビームを有し、八つの複数信号値S1〜8を計算することが可能である。送信ビームの集合に対してPにおける全部の複素信号値を得るために、以下のように個別の信号値を組み合わせることができる。
=(S*N+S*N+S*N+S*N+S*N+S*N+S*N+S*N)/8
ただしN1〜8は、ピーク音響強度値I(T1〜8)の逆数を計算することによって得られる正規化定数である。例えば、送信ビームTのPにおける強度が0.5(任意の何らかの単位)であり、送信ビームTのPにおける強度が1.0である場合、N=2及びN=1と計算することが可能である。この場合、信号値Sは信号値Sの強度の半分であると予想し、信号値Sに2を乗算して、これに信号値Sと同等の総和への寄与を与える。
上の組み合わせる式では、複数のファクタを考慮していない。特定の送信ビームに対するピーク音響強度I(T)は、組み合わされた信号にノイズを加えるだけであればそのように弱い場合がある。これを防ぐために、Sについて正規化定数Nをゼロに設定するが、ただし、I(T)は、Pにおける最大音響強度の0.05倍などの、指定された閾値より小さい。それに加えて、点Pにおける送信パルスの持続時間tb(T)は、いくつかの送信ビームTについて延長されるものとしてよく、ただし、フィールド点Pはビーム軸に沿って配置されない。これらの送信ビームからの信号を使用することで、画像分解能は低下し、そこでtb(T)の値が最短の送信パルス持続時間の1.2倍などの、指定された量より大きい場合にこれらのビームに対する正規化定数をゼロに設定する。最後に、いくつかの送信ビームについて、Pにおける送信パルスは、一次ピークからしばらくして発生する二次ピークなどの、望ましくない特性を発現させている可能性がある。この場合、二次ピークの強度の比r(T)を使用してこれらのビームの適格性評価を行うことができる。r(T)の値が0.02などの、特定の閾値より高い(二次ピークは、一次ピークの強度の2%超である)場合、正規化定数はゼロに設定されうる。これらの制約条件が追加されると、われわれの組み合わせ公式は以下のように修正されうる。
I(T)に適用される強度閾値によって適格性評価がなされる場合:
Figure 0005905080
持続時間閾値tb(T)及び/又は二次ピークの比r(T)によってさらに適格性評価がなされる場合:
Figure 0005905080
そこで、Pに対するこれらの組み合わされた複素信号は、寄与する8本のビームのうちの4本に基づくが、選択に対するわれわれの基準では、十分な強度及び適切に形成された送信パルスを有するビームのみが寄与することを許容している。適格と評価された送信ビーム領域が、それらの適格と評価された領域がオーバーラップしている領域に沿った形で図4に示されている。Pでその結果得られる組み合わされた信号は、典型的には、単一のビームからの信号と比較したフィールド点から戻るエコーのより正確な推定をもたらす。次いで、複素信号値Sがさらに処理され、これによりエコー強度などの音響画像パラメータ、又はドップラー速度測定に対する位相情報を得ることができる。
次いで、トランスデューサの視野内のさまざまなフィールド点における音響画像パラメータを集めて画像フレームを組み立てる。ディスプレイデバイスへの表示、プリンタでのプリントアウト、他のコンピューティングシステムへの送信、及び同様の操作のため、コンピューティングシステム内でフレームが処理されうる。この方法は、それぞれのフィールド点に高周波音波を当てる複数の送信ビームからさまざまなフィールド点における追加の音響情報を抽出するが、これは従来のビーム形成方法に勝る利点をもたらす。つまり、従来の処理で使用されるビームの数以下のビームの数で改善された画質が得られるか、又は画質を落とすことなくより少ない送信ビームで超音波画像フレームが形成され、これにより、より高い収集フレームレートを達成することができるということである。
より少ない送信ビームを使用してより高いフレームレートで注目するフィールドを覆う代わりにより小さな強度フィールドを持つ送信ビームのより大きな集合を使用することもできる。強度フィールドが小さければ、送信ビームは画像形成における散乱増大に寄与しうる大きな軸外れ反射面からのより少ないエコーを含む、トランスデューサの視野全体からのより少ないエコーを発生する。強度フィールドの小さい送信ビームは、それでも、オーバーラップする適格と評価された画像形成領域を有し、空間分解能が高められるという利点を維持しうる。散乱のレベルが低ければ、コントラスト分解能の増強も得られる。
複数の送信ビームからのフィールド点に対する信号を組み合わせるために、われわれは、それぞれの送信ビームに対する送信パルス波面が到来する正確な時間を知る必要がある。この時間とフィールド点で生成されるエコーから個別のアレイ素子に戻るまでの時間を知ることで、位相同期して個別の素子信号を組み合わせてフィールド点からの信号を再構成することができる。与えられた送信ビームに対するフィールド点における送信パルスの到来時間は、ピーク強度及びピーク強度持続時間と同様に、シミュレーションによって、又は実際の測定によって、決定されうる。送信ビームのそれぞれの適格と評価されたフィールド点に対する送信パルスの到来時間は、送信ビームによって生成されるエコー信号の画像形成処理時に使用するため予め決定され、メモリデバイス内に格納されうる。画像形成時に、パルスの到来時間が、フィールド点から特定のトランスデューサ素子までの移動時間に加えられ、これにより素子がエコー信号を受信した時刻を決定し、その信号と他の送信ビームから得られた信号と組み合わせるために使用する。
複数のオーバーラップする送信ビームによって高周波音波を当てられたフィールド点に対するエコー信号を組み合わせることで、フィールド点における画像空間及びコントラスト分解能を改善することができる。これは、それぞれのビームがフィールド点上でのわずかに異なる入射角を有することがあるという事実によるものであり、この結果、焦点での送信ビームの幾何学的集束によく似た合成送信集束の一形態が得られる。したがって、オーバーラップする送信ビームからの信号を組み合わせることは、画像の焦点が合わされる被写界深度を広げる効果を有し、複数の送信焦点深度を組み合わせたのと同じ結果をもたらすが、走査に沿ったそれぞれのビーム位置で複数回送信しなくてよい。これにより、フレーム収集時間を延ばさなくても画像分解能を改善できる。
それに加えて、画像形成において複数の送信ビームを組み合わせることで、コントラスト分解能を改善することができる。コントラスト分解能は、散乱及びスペックルのアーチファクトの低減により改善される。散乱は、フィールド点からのエコーと同時に受信素子に到来するフィールド点におけるエコーとは別のターゲットからのエコーによる。これらのエコーは、他の発生源に由来し、送信ビームのサイドローブ又はグレーティングローブによって生じることが多い。
異なる発生源(origin)及び入射角の送信ビームからのいくつかのエコー信号をフィールド点で組み合わせることによって、フィールド点における散乱信号は変動し、非コヒーレントに組み合わさる傾向を有する。フィールド点からのエコーは、コヒーレントに組み合わされ、したがって、複数のビームからの信号を加えることで信号対散乱比を増大させうる。超音波画像中のスペックルパターンは、送信ビーム特性にも部分的に依存し、複数のビームからの信号の組み合わせは、これらの変動も同様に平均する傾向を有する。
オーバーラップする送信ビームの集合による画像形成は、実装するにはかなり複雑な場合があり、好ましくは、以下でさらに詳しく説明されている、ソフトウェアによるピクセル指向の処理を使用して実行される。フィールド点は、画像表示のピクセル位置にあるか、又は仮想画像のピクセル位置にあるものとしてよく、最終的には、補間されて高分解能表示を形成する。それぞれのフィールド点において、受信されたエコー信号が画像形成に使用されうるオーバーラップする適格と評価された送信ビーム領域の数を決定しなければならない。集合内のそれぞれの送信ビームは、与えられたフィールド点で異なる強度を発生することができるため、戻される信号は、組み合わせる前に強度の差を補正するため正規化されるべきである。組み合わせるプロセスは、これが信号強度にも影響を及ぼすので、それぞれの送信ビームについて受信プロセスに加わるトランスデューサ素子の数を追跡し、正規化することもすべきである。これらの正規化ファクタは、注目する画像フィールド内のそれぞれのフィールド点について無関係であり、画像形成処理が進行するとともにリアルタイムで計算されるか、又はそれぞれの画像フレームの処理中にパラメータを取り出せる格納テーブルに保持されるかのいずれかとしなければならない。
画像形成処理が注目する画像フィールド内の運動を測定することを目指している場合、オーバーラップする適格と評価された送信ビーム領域の集合を使用することには、いくつかの独自の利点がある。ビームのシミュレーション又は直接測定を通じて、それぞれのフィールド点とともに送信パルス波面の入射角を計算することができるため、フィールド点における運動ベクトルから戻されるドップラー周波数偏移の大きさを知ることができる。オーバーラップする集合内の異なる送信ビームは、それぞれのフィールド点で異なる入射角を生成するように設計することができ、それぞれの入射角による測定された運動ベクトル成分の変化から、運動ベクトルの大きさの絶対値及び方向を計算することができる。次いで、これは、単一の送信ビームの入射角を使用する従来の運動撮像技術に比べてかなり正確な血流速度及び組織運動検知を行える。
従来のドップラー撮像法では、走査におけるそれぞれの送信ビーム位置に対する複数の送信及び受信信号収集を利用して、ビーム軸に沿った運動を検知するために使用されうるそれぞれの位置における受信信号の集合体を生成する。従来のアプローチでは、絶対速度ではなく、ビームの軸に沿った速度ベクトルの成分を測定する。従来のドップラー撮像方法で使用される送信ビームが、注目するフィールド点で何らかのオーバーラップをもたらすように十分に広いものである場合、同じ走査順序で絶対速度情報が得られる。送信ビーム位置における収集のそれぞれの集合体(ensemble)により、それぞれのフィールド点における運動ベクトルの単一成分を、送信ビームのすでに決定されている入射角とともに測定することができる。次の送信ビーム位置における収集の次の集合体は、前のビームとオーバーラップするフィールド点における運動ベクトルの追加の成分をもたらす。次いで、知られている送信ビームの入射角とともにそれぞれのフィールド点における複数のオーバーラップする送信ビームから得られる複数の成分を組み合わせることによって、運動ベクトルの大きさの絶対値及び方向を計算することができる。
上記の説明において、ここでは与えられた用途に対するオーバーラップする送信ビームの集合の幾何学的形状を明示的に指定してはいない。これは、オーバーラップするビームの集合の設計が、超音波用途と測定される所望の画像形成パラメータの両方に依存するからである。典型的には、送信ビームは、送信開口(transmit aperture)のサイズ、その開口に収容される送信素子の重み付け又はアポダイゼーション、及びそれぞれの送信素子の波形の開始までの時間遅延を制御することによって整形される。個別の送信機の波形も、ほとんどのシステムにおいて制御されうる。
増強された画像形成に対する送信ビームの集合を設計するための一般的ルールは、以下のとおりである。a)注目する画像フィールド内に運動がほとんど、又は全くなく、明るい鏡面反射体が最小数である場合、広く、実質的なオーバーラップのあるビームの集合を設計し、それぞれのフィールド点が5本を超えるビームによって触れられるようにする。アポダイゼーションを通じて送信ビームの焦点をソフトにするか、又は非常に深い焦点を使用して、ビームの最も狭い部分に十分なオーバーラップが生じるようにする。それに加えて、送信ビームの集合を、フィールド点の過半数において送信パルス波面の異なる入射角を与えるように設計する。b)注目するフィールド内に実質的な運動があり、及び/又は鏡面反射体が大きい場合、運動が予想される領域に少量のオーバーラップを有する比較的狭い送信ビームの集合を設計し、領域内のフィールド点におけるオーバーラップしたビーム数を2本又は3本のビームに減らす。これにより、媒体の運動による複数のビームを使用するコヒーレントな画像形成における誤差を最小にする。
上記の画像形成方法の主な利点は、これにより、超音波画像エンジニアが異なる撮像用途及び測定に合わせてシステムの撮像性能を最適化することができる点である。フレーム毎にわずかな送信ビームのみを使用して、非常に高いフレームレートを有することができ(毎秒100フレームを超える)、それでも、妥当な画質を得られるか、又はより多くのビームで典型的なフレームレート(毎秒約30フレーム)を得られ、最良の画像空間及びコントラスト分解能を得られる。ドップラー撮像では、血流又は組織運動の大きさ及び方向は、フレームレートを犠牲にすることなく得ることができる。
ピクセル指向の処理
前記の方法は、本開示の一実施形態によるソフトウェアベースの方法及びシステムアーキテクチャを利用するコンピューティングシステムで実装されうる。このシステムでは、すべてのリアルタイム処理機能をソフトウェアで実行する。提案されているアーキテクチャは、図5に概略として示されている。
ソフトウェアベースのシステムにおける唯一のカスタムハードウェアコンポーネントは、パルス発生及び信号収集回路、並びに信号データを格納するために使用される拡張メモリの大きなブロックを収容するコンピュータの拡張バスに接続する収集モジュールである。信号収集プロセスは、送信パルスに続いてトランスデューサ素子のそれぞれから戻された信号を増幅し、2値化することからなる。典型的には、トランスデューサそれ自体が備える自然なバンドパスフィルタリング以外の、2値化の前の信号のフィルタリングだけが、A/D変換を行うためのローパスアンチエイリアス処理フィルタリングである。信号は、関わっている超音波周波数と一致する一定のレートでサンプリングされ、2値化されたデータは、最小の信号処理によりメモリに格納される。信号収集の直截的(straight-forward)な設計により、回路は、比較的少ないボード面積で市販のコンポーネント使って実装することができる。
収集モジュールのより詳細な外観は、図6に概略が示されている。複数の収集チャネルが図示されており、それぞれ送信器、受信器プリアンプ、A/Dコンバータ、及びメモリブロックから構成されている。受信時に、トランスデューサ信号が2値化され、個別のメモリブロックに直接書き込まれる。メモリブロックはデュアルポートであり、これは、収集データがA/Dコンバータ側から書き込まれるのと同時にコンピュータ側から読み出されうることを意味する。メモリブロックは、システムのCPU(複数可)には通常の拡張メモリとして見える。システムは好ましくはカスタムエンクロージャ内に収納されるので、収集モジュールのサイズは、標準のコンピュータ拡張カードの通常のサイズに限定されないことに留意されたい。また、多数のトランスデューサ素子に対応できるように複数の収集モジュールが使用され、それぞれのモジュールはトランスデューサ開口の部分集合を処理することができる。
増幅器、A/Dコンバータ及び関連するインターフェース回路、並びに送信パルス発生及び信号収集に必要なコンポーネントを含む、収集モジュール用のコンポーネントは、容易に入手できる市販のコンポーネントであり、本明細書では詳しく説明しない。受信されたエコーから得られるエコー信号のRFデータの格納に必要なメモリブロックは、市販の拡張メモリカードに見られるのと本質的に同じ回路であるが、ただし、2値化された信号データを書き込むために第2のダイレクトメモリアクセスポートが加えられている。(受信されたエコー信号データは、トランスデューサによって生成される高周波電気振動からなるため、RFデータと一般的に称される。)
メモリは、中央演算処理装置のアドレス空間内にマッピングされ、コンピュータのマザーボード上に配置されている他のCPUメモリと同様にしてアクセスされうる。或いは、RFデータは、ダイレクトメモリアクセスを使って収集モジュールからホストコンピュータに転送することもできる。収集モジュール上のメモリのサイズは、最大256以上の別々の送信/受信サイクルに対して個別のチャネル受信データを受け入れることができるようなサイズである。人体内の超音波パルスの往復に対する最大の実用的な浸透深さは、波長約500個分なので、中心周波数の4倍の標準サンプリングレートは、個別のトランスデューサ素子からの4000サンプル程度のストレージを必要とする。16ビット及び128個のトランスデューサチャネルのサンプリング精度では、最大深さの受信データ収集は、それぞれの送信/受信イベントについて約1メガバイトのストレージを必要とする。したがって、256個のイベントを格納するには、256MBのストレージが必要であり、総計で、128個のチャネルシステムを一つ又は二つと少ない収集モジュール上に構築することが可能である。
ソフトウェアベースの超音波システムの別の態様は、コンピュータのマザーボード及びその関連するコンポーネントである。提案されている設計に対するマザーボードは、好ましくは、必要な処理能力を得るために、マルチプロセッサCPU構成をサポートすべきである。電源、メモリ、ハードディスクストレージ、DVD/CD−RWドライブ、及びモニターを完備した、完全なマルチプロセッサコンピュータシステムは、当業者に知られており、容易に市販のものを購入できるので、これ以上詳しく説明することはしない。
ソフトウェアベースの超音波システムは、「高性能」、つまり、健康産業に著しいメリットをもたらすために、既存のハイエンドシステムに匹敵する画質を真に達成しなければならない。このレベルの性能は、現在のシステムのフロースルー処理方法をソフトウェア実装に単純に変換することでは達成できないが、それは、フロースルーアーキテクチャでは1秒分のリアルタイム撮像に必要なすべての処理演算を単純に追加すると、数は複数の汎用プロセッサで現在達成可能な1秒あたりの演算数を超えるからである。その結果、フロースルー方法に比べてはるかに効率のよい新しい処理方法が求められている。
本発明のソフトウェアベースの超音波システムアーキテクチャの一実施形態では、信号及び画像処理の入力データは、1以上の送信イベントに続く個別のトランスデューサチャネルから得られたRFサンプルの集合からなる。そこで、例えば、図8に示されているように、128素子線形トランスデューサアレイによる典型的な2D撮像走査モードを考える。
この場合、「送信イベント」は、媒体内で組み合わされる複数の音波を生成し特定の素子位置においてトランスデューサ上の原点から外へ放射される集束超音波ビームを形成する複数のトランスデューサ素子からの、時間が指定されているパルス(timed pulse)からなる。複数の送信イベント(全部で128個)は、トランスデューサ面の幅にわたって徐々に順次放射される超音波ビームを生成し、これにより、画像フレーム全体のインテロゲーションを行う。これらの送信ビームのそれぞれについて、受信エコーデータは、トランスデューサ内の128個の受信器素子のそれぞれから集められて、それぞれの列が対応するトランスデューサ素子によって受信されるサンプリングされたエコー信号を表すデータ配列にまとめられる。こうして、それぞれの配列は、128個のトランスデューサ素子に対応する128個の列と、取り出された深さのサンプルの数に対応する多数の行(この場合、4096個の行を仮定するとサンプルは4096となる)を有する。次いで、これらの128個のデータ配列は、一つの完全な画像フレームを生成するのに十分なRFデータセットを構成する。
フロースルーアーキテクチャでは、上述のRFデータセットは、ビーム及び画像形成がトランスデューサからのデータストリームとして行われるので、存在することすらしない(少なくとも一度には全くない)ことに注目すべきである。言い換えると、データが送信イベントの後にそれぞれの素子に戻ると、これらは処理され、組み合わされて(ビーム形成と称される)、単一のビーム(走査ライン)にそって集束戻り信号を表す単一のRF信号を生成する。このRF信号は、エコー振幅サンプルに処理され(ここでもリアルタイムで)、メモリアレイ内に格納される。すべてのビーム方向が処理されたときに、エコー振幅データは表示できるように補間されてピクセル画像にフォーマットされる。すべての処理がリアルタイムで行われるので、処理回路は、トランスデューサ素子から入るデータストリーミングに「遅れずについてくる」ことができなければならない。
本発明のソフトウェアベースのアーキテクチャでは、すべての入力データは、処理前に格納される。これは、収集レートを処理レートから切り離し、したがって、必要ならば、処理時間を収集時間よりも長くとることができる。これは、収集の深さが短く、サンプルレートが高い、高周波走査でははっきりと有利である。例えば、10MHzの走査ヘッドは、約4センチメートルの撮像の使用可能な深さを有するものとしてよい。この場合、組織中の音速により、128個の送信/受信イベントのそれぞれは、52マイクロ秒でデータを取得して格納することを命令することになるが、これは非常に高速な収集データ転送速度である。フロースルーアーキテクチャでは、これらの収集データは、高処理速度でリアルタイムで走査ラインに形成される。本発明のソフトウェアベースのアーキテクチャでは、RFデータを格納することにより、処理の時間を、組織移動のリアルタイムの可視化のために典型的には33ミリ秒(30フレーム/秒)である、表示のフレーム期間と同程度の時間とすることを可能にする。128個のピクセル列(走査ラインにおおよそ類似している)について、これは、フロースルーアーキテクチャの52マイクロ秒ではなく、1列あたり258マイクロ秒の処理を可能にする。この格納戦略は、典型的な走査深さに対してフロースルーアーキテクチャと比較して処理の最大速度を実質的に下げる効果を有する。
入力データの格納により、最大処理速度は低下するが、必ずしも処理ステップの数が減るわけではない。これを達成するために、超音波データ処理に対する新しいアプローチがとられる。第1のステップは、撮像モードに入っているときのシステムの最終的な目標は、出力ディスプレイ上に画像を生成することであることを認識することである。超音波画像は、周波数及び配列次元などの、収集システムの物理的パラメータに依存する基本的分解能を有し、エコー振幅又は他の何らかの組織(音響)特性を符号化したピクセル値の矩形配列として表されうる。この矩形ピクセルの密度は、画像分解能の適切な空間サンプリングを行えるものでなければならない。表示画像は、ピクセルの矩形配列のみからなる必要はないが、異なる幾何学的形状を表す、ピクセルの任意の集合からなるものとしてよいことは理解されるであろう。
次のステップは、この画像配列内のピクセルの一つから始めて、RFデータセット内のどのサンプル点がこのピクセルの強度の計算に寄与するかを考慮し、これらにアクセスし処理する最も効率のよい方法を決定することである。このアプローチは、現在のフロースルーアーキテクチャによって利用されているものと完全に異なるアプローチであるが、それは、ディスプレイ上のピクセルに寄与する情報のみが処理されればよいからである。本発明のアプローチでは、表示画像上の小さな領域では大きな画像領域に比べて全体的処理時間が短縮されるが、それは、小さな領域に含まれるピクセルの数が少ないからである。対照的に、フロースルー処理方法は、画像領域のサイズと無関係に、最大データストリーム帯域幅を取り扱えるように設計されなければならない。
超音波画像を適切に表すために必要なピクセル配列を処理した後、この配列をコンピュータのディスプレイに対して、表示に適したサイズでレンダリングすることができる。追加CPU処理を必要としない、コンピュータのグラフィックスプロセッサは、典型的には、この演算を実行することができ、これは単純なスケーリングと補間からなる。
次に、われわれの超音波画像の単一ピクセルに対する処理戦略を考える。この説明では、トランスデューサアレイに関してピクセルの対応する空間位置でエコー強度を得ることが目的であると仮定する。他の音響パラメータも同様に得ることができる。われわれの第1のステップは、エコー強度計算に寄与するサンプルを含む収集RFデータの領域を見つけることである。図7の操作方法に対してこれを遂行するために、最初に、ピクセル位置との交差点に最も近くなる収集走査ラインを見つけて、対応する個別の素子データ配列を使用する。
図8は、超音波画像内の例示的なピクセルに対するこのマッピングプロセスを示している。図8では、指し示されているピクセルは、走査の最も近い収集ラインにマッピングされ、この場合これは走査ライン4であり、RFデータは第4の個別の素子RFデータ配列(第4の送信/受信イベントから集められたデータを表す)内に置かれている。1以上のRFデータ配列が、ピクセル信号に寄与するものとして選択されうるが、この例では、単一のデータ配列のみを考察する。
われわれの次のステップは、ピクセルの強度計算に寄与するサンプルを含む個別の素子配列内の領域をマッピングすることである。このマッピングプロセスは、かなり複雑であり、複数のファクタに依存する。トランスデューサ素子はそれぞれ、画像フィールド内の特定の点から戻る信号にどのように応答するかを決定する感度(sensitivity)の領域を有する。与えられた画像点に対して、所定の閾値より高い感度を有する素子のみを考察すればよいが、それは、感度が低すぎると、素子はピクセルの量に有用な情報を付与しない。次いで、この感度閾値は、マッピングされた領域に含める素子データ列の数を決定する。
マッピングされたデータ領域の始まりの深さは、それぞれの個別のトランスデューサ素子における戻るエコーの到着時間によって決定される。図8に示されているように、画像点からさらに離れる素子に対する画像点信号は、時間的に後で捕捉され、したがって、データセットの開始点は、メモリ内のより深くにある。最後に、マッピングされるデータ領域に必要な深さ範囲は、生成される送信パルスの持続時間に依存する。より長い送信パルスは、より長い期間に画像点を励起し、RFメモリのより大きな深さ範囲にわたるエコー信号を発生する。
幸いなことに、マッピングされたデータの領域を決定することになるファクタの多くは、与えられたピクセルグリッドついて事前に計算されうるが、それは、このグリッドがリアルタイムの画像列の複数のフレームにわたって変化しないからである。事前に計算されたファクタを使用することで、与えられたピクセルに対するマッピングされたデータ領域は、素早く、効率的に決定されうるため、リアルタイムの撮像の際にかなりの計算を削減できる。
ピクセルマッピングRFデータを選び出した後、これを以下に示されているように行列RFPnmに構成することができる。
「Pnm」という記法は、第n行、第m列の画像ピクセルを指す。行列の列は、図11の垂直バーであり、
Figure 0005905080
であると仮定し、それぞれの垂直バー内のサンプルの数jは、同じである。サンプルの数jは、送信パルスによって生成される信号を捕捉するのに必要な時間内でRFデータの範囲に依存する。インデックスkは、強度計算に関わる画像点からの適切な信号強度を有するRFデータ配列内のチャネルの数である。ピクセルPnmの信号強度値を計算するプロセスは、一連の行列演算からなり、最終的に単一の値が得られる。
図9は、本開示のプロセスを実装するように適合可能なシステムアーキテクチャの高水準表現の図であり、図10は、一実施形態のソフトウェアベースのアーキテクチャの概略図である。それに加えて、図11は、上で説明されているピクセル指向処理に従って形成されたプラグインモジュールの線図である。
より具体的には、図9は、本開示のプロセスを実装する高水準システムアーキテクチャ70を表すシステムレベルのブロック図を示している。これは、単に代表的な一実施形態にすぎず、例示されているアーキテクチャ70は、本開示のすべての実施形態に対する要件ではないことは理解されるであろう。
アーキテクチャ70は、ホストコンピュータ72を備え、ホストコンピュータ72は、PCI−express74を介してマルチチャネルトランシーバ及びデータ収集システム76に結合されている。ホストコンピュータ72は、ユーザーインターフェース及び制御装置78、並びにディスプレイ80を有し、両方ともピクセルベースのアプリケーション処理ソフトウェア84を利用するプロセッサ82に結合されている。マルチチャネルトランシーバ及びデータ収集システム76のハードウェアは、超音波トランスデューサ86に結合されており、この超音波トランスデューサ86は、モニター、プロジェクターなどのディスプレイ80上に表示するか、又はデバイスの表示若しくは操作若しくはその両方のため別のデバイスに送信することを目的として音響媒体内の領域88を撮像するために使用されるものである。これらのコンポーネントは、容易に市販のものが入手可能であるため、本明細書では詳しく説明しない。
ピクセル指向の処理を使用することで、オーバーラップする送信ビームを利用するすでに述べられている種類の複素エコー信号再構成を行うことができる。この方法では、トランスデューサの視野に関して位置決めされた点のグリッドのそれぞれにおいて送信ビームの集合内のそれぞれの送信ビームの計算された又は測定されたビーム特性を格納するためにルックアップテーブルメモリが使用される。ピクセル指向の信号再構成は、点のグリッド内のそれぞれの点で実行され、それぞれの送信ビームのエコー信号の寄与は、再構成点に対応するルックアップテーブルパラメータを使用して計算され、組み合わされる。
上で説明されているさまざまな実施形態は、さらなる実施形態を実現するために組み合わされうる。実施形態の態様は、さらなる実施形態を実現するためにさまざまな特許、出願、及び公開の概念を採用する必要がある場合に、修正されうる。
上述の説明に照らしてこれら及び他の変更を実施形態に加えることができる。一般に、添付の特許請求の範囲において、使用されている用語は、特許請求の範囲を明細書及び特許請求の範囲で開示されている特定の実施形態に限定するものとして解釈すべきでなく、そのような特許請求の範囲が関わる均等物の全範囲とともにすべての可能な実施形態を含むものとして解釈すべきである。したがって、特許請求の範囲は、開示によって限定されない。

Claims (14)

  1. 多素子トランスデューサから得られた超音波画像フレームの分解能及びフレームレートを改善する方法であって、
    a.画像フレームを収集するために使用される送信ビームの集合内のそれぞれの送信ビームに対する送信開口及び個別の素子送信特性を指定するステップであって、前記送信ビームの集合が、所望の視野の少なくとも一部に高周波音波を当てる少なくとも一つの送信ビームを含む、ステップと、
    b.測定又はシミュレーションを通じて、前記視野を覆う複数のフィールド点で前記送信ビームの集合内のそれぞれの送信ビームのさまざまな属性を決定するステップと、
    c.それぞれの新しい画像フレームを処理する際の繰り返し使用に対するそれぞれのフィールド点の前記属性を格納するステップと、
    d.前記送信ビームの集合を使用してエコー信号を送受信し、前記エコー信号をメモリに格納するステップと、
    e.前記格納されている送信ビームの属性のうちの1つ以上を使用して前記格納されているエコー信号を処理し、画像形成で使用するため与えられた送信ビームが寄与するフィールド点から受信された前記超音波エコー信号の適格性評価を行い、前記信号がどのように処理されるべきかを指定するステップと、
    f.フィールド点に対してオーバーラップする、適格と評価された前記送信ビームの集合から複数の画像パラメータを生成して組み合わせ、前記画像フレームに対するフィールド点の値を構成する最終画像パラメータを生成するステップと
    を含む方法。
  2. フィールド点における前記送信ビームの前記測定された、又はシミュレートされた属性は、
    a.送信イベントの持続時間にわたる前記送信ビームのピーク強度、
    b.前記送信ビームの前記ピーク強度が前記フィールド点で生じた時刻、
    c.パルス強度が前記パルスピークの−20dB未満である閾値を超える時間によって決定される送信パルス持続時間、
    d.前記ピーク強度と、前記フィールド点で送信イベントの発生中に生じるビーム強度の他のピークとの比、
    e.前記フィールド点の位置を指定するために使用される選択された軸を持つ前記送信ビームパルスの波面の入射角
    のうちの少なくとも一つを含む、請求項1に記載の方法。
  3. 前記格納されているエコー信号に対する前記処理は、前記フィールド点におけるそれぞれの送信ビームのピーク送信強度を使用することによってそれぞれの格納されているエコー信号の利得を補償することを含む、請求項1に記載の方法。
  4. 前記フィールド点は、表示される画像フレーム内のピクセルの対応する位置における矩形配列で配置される、請求項1に記載の方法。
  5. 前記送信ビームの集合は、集束しない送信ビームを使用して送受信収集が64未満である視野に高周波音波を当てることによってフレームレートを最大化するように構成される、請求項1に記載の方法。
  6. 前記送信ビームの集合は、前記集合内の前記送信ビームのオーバーラップの量を制御することによって最適な空間及びコントラスト分解能をもたらすように構成される、請求項1に記載の方法。
  7. 組織又は血流速度に対する受信したエコー信号の処理では、入射送信ビーム属性の格納されている角度を使用して、ドップラー周波数偏移に基づき速度計算を補正し、それぞれのフィールド点における運動の大きさと運動の方向を求めるように構成される、請求項1に記載の方法。
  8. 多素子トランスデューサから得られた超音波画像の分解能及びフレームレートを改善するためのコンピュータで実装されるシステムであって、
    プロセッサ及びメモリを備え、前記プロセッサ及び前記メモリは、
    a.画像フレームを収集するために使用される送信ビームの集合内のそれぞれの送信ビームに対する送信開口及び個別の素子送信特性を指定する処理であって、前記送信ビームの集合が、所望の視野の少なくとも一部に高周波音波を当てる少なくとも一つの送信ビームを含む、処理と、
    b.測定又はシミュレーションを通じて、前記視野を覆う複数のフィールド点で前記送信ビームの集合内のそれぞれの送信ビームのさまざまな属性を決定する処理と、
    c.それぞれの新しい画像フレームを処理する際の繰り返し使用に対するそれぞれのフィールド点の前記属性を格納する処理と、
    d.前記送信ビームの集合を使用してエコー信号を送受信し、前記エコー信号をメモリに格納する処理と、
    e.前記格納されている送信ビーム属性のうちの1つ以上を使用して前記格納されているエコー信号を処理し、画像形成で使用するため与えられた送信ビームが寄与するフィールド点から受信された前記超音波エコー信号の適格性評価を行い、前記信号がどのように処理されるべきかを指定する処理と、
    f.フィールド点に対してオーバーラップする、適格と評価された前記送信ビームの集合から複数の画像パラメータを生成して組み合わせ、前記画像フレームに対するフィールド点の値を構成する最終画像パラメータを生成する処理と
    を実行するように構成されている、システム。
  9. 前記送信ビームの前記測定された、又はシミュレートされた属性は、
    a.送信イベントの持続時間にわたる前記送信ビームのピーク強度、
    b.前記送信ビームの前記ピーク強度が前記フィールド点で生じた時刻、
    c.パルス強度が前記パルスピークの−20dB未満である閾値を超える時間によって決定される送信パルス持続時間、
    d.前記ピーク強度と、前記フィールド点で送信イベントの発生中に生じるビーム強度の他のピークとの比、
    e.前記フィールド点の位置を指定するために使用される選択された軸を持つ前記送信ビームパルスの波面の入射角
    のうちの少なくとも一つを含む、請求項8に記載のシステム。
  10. 前記格納されているエコー信号に対する前記処理は、前記フィールド点におけるそれぞれの送信ビームのピーク送信強度を使用することによってそれぞれの格納されているエコー信号の利得を補償することを含む、請求項8に記載のシステム。
  11. 前記フィールド点は、表示される画像フレーム内のピクセルの対応する位置における矩形配列で配置される、請求項8に記載のシステム。
  12. 前記送信ビームの集合は、集束しない送信ビームを使用して送受信収集が64未満である視野に高周波音波を当てることによってフレームレートを最大化するように構成される、請求項8に記載のシステム。
  13. 前記送信ビームの集合は、前記集合内の前記送信ビームのオーバーラップの量を制御することによって最適な空間及びコントラスト分解能をもたらすように構成される、請求項8に記載のシステム。
  14. 組織又は血流速度に対する受信したエコー信号の処理では、入射送信ビーム属性の格納されている角度を使用して、ドップラー周波数偏移に基づき速度計算を補正し、それぞれのフィールド点における運動の大きさと運動の方向を求めるように構成される、請求項8に記載のシステム。
JP2014509393A 2011-05-02 2012-05-02 オーバーラップする送信ビームにおける適格と評価された領域を使用する増強された超音波画像形成 Active JP5905080B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201161481476P 2011-05-02 2011-05-02
US61/481,476 2011-05-02
PCT/US2012/036155 WO2012151300A2 (en) 2011-05-02 2012-05-02 Enhanced ultrasound image formation using qualified regions of overlapping transmit beams

Publications (2)

Publication Number Publication Date
JP2014512923A JP2014512923A (ja) 2014-05-29
JP5905080B2 true JP5905080B2 (ja) 2016-04-20

Family

ID=47108208

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014509393A Active JP5905080B2 (ja) 2011-05-02 2012-05-02 オーバーラップする送信ビームにおける適格と評価された領域を使用する増強された超音波画像形成

Country Status (7)

Country Link
US (1) US9384530B2 (ja)
EP (2) EP3570068A1 (ja)
JP (1) JP5905080B2 (ja)
KR (1) KR102014946B1 (ja)
CN (1) CN103765240B (ja)
CA (1) CA2834993C (ja)
WO (1) WO2012151300A2 (ja)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5905856B2 (ja) * 2012-07-13 2016-04-20 富士フイルム株式会社 超音波検査装置
JP6006249B2 (ja) * 2014-03-24 2016-10-12 富士フイルム株式会社 音響波処理装置、音響波処理装置の信号処理方法およびプログラム
JP6129108B2 (ja) * 2014-03-31 2017-05-17 富士フイルム株式会社 音響波処理装置、音響波処理装置の信号処理方法およびプログラム
WO2015180069A1 (zh) 2014-05-28 2015-12-03 深圳迈瑞生物医疗电子股份有限公司 一种超声成像方法和***
US9945946B2 (en) 2014-09-11 2018-04-17 Microsoft Technology Licensing, Llc Ultrasonic depth imaging
US9479694B2 (en) * 2014-10-28 2016-10-25 Google Inc. Systems and methods for autonomously generating photo summaries
US11583254B2 (en) 2017-06-28 2023-02-21 Duke University Systems and methods for intra-beam compounding for ultrasound imaging
CN117064447B (zh) * 2023-10-13 2024-01-05 深圳华声医疗技术股份有限公司 超声穿刺针显影方法、装置、超声设备及存储介质

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4159462A (en) * 1977-08-18 1979-06-26 General Electric Company Ultrasonic multi-sector scanner
US5908391A (en) * 1997-05-07 1999-06-01 General Electric Company Method and apparatus for enhancing resolution and sensitivity in color flow ultrasound imaging using multiple transmit focal zones
KR100330855B1 (ko) * 1999-02-09 2002-04-03 이민화 초단수신신호의 저장과 저장된 데이터의 재사용이 가능한 디지털 초음파영상장치
US6056693A (en) * 1999-08-16 2000-05-02 General Electric Company Ultrasound imaging with synthetic transmit focusing
US6537217B1 (en) * 2001-08-24 2003-03-25 Ge Medical Systems Global Technology Company, Llc Method and apparatus for improved spatial and temporal resolution in ultrasound imaging
KR100419806B1 (ko) * 2001-12-31 2004-02-21 주식회사 메디슨 평면파를 이용하는 초음파 영상의 합성 구경 집속 방법
JP4232975B2 (ja) 2004-07-01 2009-03-04 株式会社リコー 画像形成方法、画像形成装置および画像形成装置用プロセスカートリッジ
US8287456B2 (en) 2005-04-14 2012-10-16 Verasonics, Inc. Ultrasound imaging system with pixel oriented processing
KR100868483B1 (ko) * 2005-10-07 2008-11-12 주식회사 메디슨 초음파영상 디스플레이 방법
US20090048440A1 (en) 2005-11-03 2009-02-19 Neose Technologies, Inc. Nucleotide Sugar Purification Using Membranes
WO2007058895A2 (en) * 2005-11-11 2007-05-24 Visualsonics Inc. Overlay image contrast enhancement
BRPI0711595A2 (pt) * 2006-05-12 2011-11-16 Koninkl Philips Electronics Nv sistema de formação de imagem de diagnóstico ultra-sÈnico, e, método para produzir uma imagem de ultra-som com uma gama focal estendida.
US8317712B2 (en) * 2006-05-12 2012-11-27 Koninklijke Philips Electronics N.V. Eindhoven Retrospective dynamic transmit focusing for spatial compounding
WO2009064891A2 (en) 2007-11-13 2009-05-22 Wisconsin Alumni Research Foundation A method for producing highly constrained ultrasound images
US8241216B2 (en) * 2008-06-06 2012-08-14 Siemens Medical Solutions Usa, Inc. Coherent image formation for dynamic transmit beamformation
DK2303131T3 (en) 2008-06-26 2015-05-04 Verasonics Inc Quantitative Doppler flow imaging with high frame rate using unfocused transmission rays

Also Published As

Publication number Publication date
WO2012151300A2 (en) 2012-11-08
CN103765240A (zh) 2014-04-30
CN103765240B (zh) 2017-02-15
CA2834993C (en) 2019-04-30
JP2014512923A (ja) 2014-05-29
US20140140600A1 (en) 2014-05-22
WO2012151300A3 (en) 2013-01-03
KR20140031287A (ko) 2014-03-12
EP3570068A1 (en) 2019-11-20
EP2705383A2 (en) 2014-03-12
KR102014946B1 (ko) 2019-08-27
US9384530B2 (en) 2016-07-05
CA2834993A1 (en) 2012-11-08
EP2705383A4 (en) 2014-11-05

Similar Documents

Publication Publication Date Title
JP5905080B2 (ja) オーバーラップする送信ビームにおける適格と評価された領域を使用する増強された超音波画像形成
US9649094B2 (en) Ultrasound imaging system with pixel oriented processing
KR101820422B1 (ko) 포커싱되지 않은 송신 빔들을 이용한 높은 프레임 레이트의 정량적 도플러 흐름 촬영
KR101581227B1 (ko) 로컬 미디어 움직임의 감지에 기초한 적응적 초음파 이미지 복원
Montaldo et al. Coherent plane-wave compounding for very high frame rate ultrasonography and transient elastography
CN111819467B (zh) 波传播和散射参数的估计方法和仪器
CN103505243A (zh) 测量超声波的声吸收或衰减
CA3127900A1 (en) Method and system for ultrasonic characterization of a medium
CA3127822A1 (en) Method and system for ultrasonic characterization of a medium
CA3127924A1 (en) Method and system for ultrasonic characterization of a medium
US10792518B2 (en) System and device for improved ultrasound cavitation mapping
US9759690B2 (en) Method and system for nondestructive ultrasound testing
CN108784744A (zh) 用于剪切波成像的可变焦点
KR20220036337A (ko) 매질의 초음파 특성화 방법 및 시스템
Guenther et al. Synthetic aperture methods for angular scatter imaging

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20140428

RD01 Notification of change of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7426

Effective date: 20140711

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20140711

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20150226

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150303

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20150526

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150903

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20160216

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160315

R150 Certificate of patent or registration of utility model

Ref document number: 5905080

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250