[go: up one dir, main page]

JP2010508881A - 最適なスペックル追跡のための二重経路処理 - Google Patents

最適なスペックル追跡のための二重経路処理 Download PDF

Info

Publication number
JP2010508881A
JP2010508881A JP2009535182A JP2009535182A JP2010508881A JP 2010508881 A JP2010508881 A JP 2010508881A JP 2009535182 A JP2009535182 A JP 2009535182A JP 2009535182 A JP2009535182 A JP 2009535182A JP 2010508881 A JP2010508881 A JP 2010508881A
Authority
JP
Japan
Prior art keywords
processing
speckle
scan line
data
sound waves
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.)
Granted
Application number
JP2009535182A
Other languages
English (en)
Other versions
JP5627890B2 (ja
Inventor
カール イー ディエル
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
Koninklijke Philips Electronics NV
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 Koninklijke Philips NV, Koninklijke Philips Electronics NV filed Critical Koninklijke Philips NV
Publication of JP2010508881A publication Critical patent/JP2010508881A/ja
Application granted granted Critical
Publication of JP5627890B2 publication Critical patent/JP5627890B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Clinical applications
    • 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/52025Details of receivers for pulse systems
    • G01S7/52026Extracting wanted echo signals
    • 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/8995Combining images from different aspect angles, e.g. spatial compounding
    • 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/52053Display arrangements
    • G01S7/52057Cathode ray tube displays
    • G01S7/52071Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric 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
    • 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/52077Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging with means for elimination of unwanted signals, e.g. noise or interference

Landscapes

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

Abstract

本発明は一般に、組織といった標的の時間に対する運動の一連の画像を得るための、スペックル追跡に関する技術を強調すること及び軽減することを結合する改良されたシステム及び方法に関する。上記方法は、人体に音波を送信するステップと、これらの音波のエコーを出力するステップと、スキャンラインデータを生み出すため上記エコーを受信及びビーム形成するステップと、生体構造情報を表示するためスペックルを減らす方法を用いて上記スキャンラインデータを処理するステップと、スペックルを減らさない方法又は手順を用いてスキャンラインデータを処理するステップと、1つのスキャンシーケンスの間、スペックルを減らす処理がされたデータとスペックルを減らすことなく処理されたデータとの2つのスキャンラインデータを同時に取得するステップとを有する。

Description

本発明は、一般に超音波画像診断に関し、より詳細には超音波スペックルパターンの強調及び軽減を両方用いる超音波画像診断に関する。
過去10年にわたり、超音波画像品質に対する有意な改良は、周波数合成及び空間的合成(SonoCT)技術を含む高度な合成技術によりもたらされた。これらの技術は超音波スペックルを軽減することにより機能する。超音波スペックルは、ローリー散乱状態のエコーから得られる建設的/破壊的な干渉パターンに関連付けられる人工的ノイズパターンである。
スペックルは、超音波ビームの分解セルに含まれる多数の小さな生体構造的目標物に関連付けられるランダムな建設的及び破壊的な干渉によりもたらされる。これらの目標物又はローリー散乱は、定義上、検査に使用する音波の波長より非常に短い。送信された音ビームは、広帯域となる傾向がある。これは、このビームがさまざまな波長を備える音波を含む概念に関連する。当業者に知られているように、異なる波長は、異なる建設的及び破壊的な干渉パターンを持ち、従って、異なるスペックルパターンを持つ。プリズムがその成分の波長(カラー)へと白色光を分離する方法に非常に似た態様で、直角位相帯域フィルタは、戻り音エコーを2つのグループに分離する。1つは、より短い波長を持ち、もう1つはより長い波長を持つグループである。従って、2つのグループは、異なる干渉パターンを持ち、従って、異なるスペックルパターンを持つことになる。
近年では、(1D、2D及び3D空間において)血液及び組織の速度及び変位両方を追跡したいという願望があった。超音波画像診断から得られるスペックルパターンは、短い距離の組織間変位を追跡する傾向があるので、組織速度及び変位の正確な測定は、空間にわたり得られるスペックルパターンを時間にわたり得られる類似するスペックルパターンと相互相関させることにより算出されることができる。これらの技術は産業界では、2Dスペックル追跡及び3Dスペックル追跡と呼ばれる。戻りエコーの超音波スペックルを軽減することにより最適な白黒の(BW)画像品質が得られ、超音波スペックルが強調されるとき、最適なスペックル追跡(変位及び速度)が得られる。
本発明は一般に、例えば組織といった目標物のある時間にわたる運動の一連の画像を得るための、スペックル追跡に関する強調及び軽減技術を組み合わせる改良されたシステム及び方法に関連する。
上記方法は、人体に音波を送信し、これらの音波のエコーを出力するステップと、スキャンラインデータを生み出すため上記エコーを受信及びビーム形成するステップと、生体構造情報を表示するためスペックルを減らす方法を用いて上記スキャンラインデータを処理するステップと、スペックルを減らさない方法又は手順を用いてスキャンラインデータを処理するステップと、1つのスキャンシーケンスの間、スペックルを減らす処理がされたデータとスペックルを減らすことなく処理されたデータとの2つのスキャンラインデータを同時に取得するステップとを有する。
低周波直角位相帯域フィルタから生み出される患者の一部(組織)の説明的なスペックル画像を示す図である。 患者の一部(組織)の説明的な合成ファントム又は現実画像を示す図である。 低周波直角位相帯域フィルタから生み出される患者の一部(組織)の別の説明的なスペックル画像を示す図である。 後の時間での図3の組織を示す図である。 患者の一部(組織)の別の説明的な合成ファントム又は現実画像を示す図である。 図5に示される組織の後の時間での現実画像を示す図である。 スペックルパターンを減らすよう構成される従来技術の超音波イメージングシステムの説明的な概略図を示す図である。 最適なスペックル追跡のために構成される従来技術の超音波イメージングシステムの説明的な概略図を示す図である。 本発明のある実施形態による、最適なスペックル追跡のために構成される超音波イメージングシステムの説明的な概略図を示す図である。 本発明の別の実施形態による、最適なスペックル追跡のために構成される超音波イメージングシステムの説明的な概略図を示す図である。
本発明についての前述した及び他の目的、側面、特徴及び利点は、以下の説明及び請求項から一層明らかになるだろう。
図面において、同様な参照文字は一般に、異なる図にわたり同じ部分を参照する。また、図面は必ずしも大きさ通りに描かれておらず、代わりに、一般的に本発明の原理を示す部分が強調される。
本発明は、超音波スペックルパターンの強調及び軽減両方を用いる超音波イメージングに関する。単一のスキャンラインに対応するデジタル化された戻りエコーが複製され、2本の別々の処理経路に送られる。1本の経路は、白黒(BW)画像品質のために最適化される。即ち、スペックルが減じられる。他の経路は、スペックル追跡のために最適化される。即ち、スペックルが強調される。
目標物(例えば、人体組織)が超音波で照射されるとき、その目標物は、建設的又は破壊的な態様で超音波信号と干渉する可能性がある。標的組織の画像は、粒状に見えるか、又はテクスチャを持つように見える。この粒状の外観はスペックルと呼ばれる。スペックルは、画像において基礎をなすデータと何ら関係がない。スペックルは、組織が動くにつれ変化する、データ内の単なる任意のバンプ又はノイズである。従って、スペックルを追跡すること、即ち、時間におけるスペックルデータをキャプチャすることは、時間における組織運動及び/若しくは変位、又は血流の追跡を可能にする。
例えば、スペックルは、以下のようにして、心臓周期を通しての心拍動又は運動を追跡するのに用いられることができる。血液が組織の中を適切に流れるとき、その組織は柔らかい。血液が組織の中を適切に流れないとき、組織は硬くなる。心臓は、スポンジ状であり、収縮性の特性を持つ。心臓が鼓動するとき、それは圧縮及び拡張する。しかしながら、死亡組織又は傷害性組織は、圧縮又は運動をしない。従って、ある時間での心臓の超音波画像からスペックルパターンを追跡すると、心臓の拍動若しくは運動、又はその欠如が追跡されることができる。
図1は、低周波直角位相帯域フィルタから生み出される、患者の一部、即ち組織のスペックル画像を示す。図2は、患者の一部、すなわち組織の合成ファントム又は現実画像を示す。図2において、人工的スペックルは全て取り除かれている。このことは、左側にある点目標、左上にある小さな黒い血管、及び背景のグレーレベル(右下)におけるわずかな変動が容易に確認できる状態にする。従って図2は、2D表示及び/又は生体構造表示にとって最適であると考えられることができる。しかしながら、組織が運動する場合、及び、ある位置におけるこの変位の相対的な変化を検出したい場合、図2は有益でない。なぜなら、それは任意の有意な「テクスチャ」(特に右下において)を欠いているからである。従って、図2を用いて運動を検出することは、非常に困難である。
図3において、中央におけるグレーボックスにより示される組織の任意の領域が、追跡のために特定される。グレーボックスは、関心領域(ROI)として知られている領域を示す。図4は、ROIと共に同じ組織を、しかし、後の時間での組織を示す。ROIが示すように、組織は、図3に示されるその最初の位置から移動した。しかしより重要なことは、スペックル又はグレインとしても知られるテクスチャが、図3及び図4の両方において同じであることである。任意の所与の組織がどれくらい移動したかを決定するためのさまざまな「スペックル追跡」方法を可能にするのが、このテクスチャである。
図5は、図3と同じ時間での同じ組織を示す。しかしながら、この場合、すべてのスペックルが除去されている。この場合も、組織の特定の領域(ROI、グレーボックス)が、追跡のために特定された。図5で使用されたのと同じスペックル低減技術を用いて、図6は、後の時間での同じ組織を示す。しかしながら、図6において、すべてのスペックルが取り除かれたので、図5のROIにおいて所望の組織がどれくらい移動したかを決定するためのいかなるスペックル追跡方法も存在しない。こうして、すべてのスペックルを除去することは、組織の運動の追跡を妨げる。
図7を参照すると、スペックルパターンを減らすよう構成される従来技術の超音波イメージングシステム100の概略図が示される。イメージングシステム100は、超音波トランスデューサ(XD)105、スキャナ110、第1の直角位相帯域フィルタ(QBP1)115、第2の直角位相帯域フィルタ(QBP2)120、ログ検出部125、ログ検出部130、平均化手段135、マルチレート低域フィルタ(LPF)140、SonoCT145及びディスプレイ150を含む。好ましい実施形態において、スキャナ110が戻りエコーをデジタル化することが想定される。その結果、後続の処理ステップが、デジタルハードウエアを用いて又はCPUの一部としてのソフトウェアを用いて処理される。平均化手段135は、ログ検出部125及び130の2つの出力を合計し、及び結果を半分にするものとして単純にすることができる。
動作において、超音波トランスデューサ(XD)105は、電気信号を音波に及びその逆に変換する超音波圧電トランスデューサである。XD105は、被検者(患者)をスキャンし、超音波を生み出し、及びスキャナ110にそれらを出力する。これは、超音波ビームを方向付けし及び焦点合わせするのに用いられる波ビームフォーマに対する位相である。スキャナ110の出力は、QBP1 115及びQBP2 120に入力される。QBP1 115及びQBP2 120は、各々がヒルベルト変換器(1〜3MHz)を含む帯域フィルタである。QBP1 115は2MHzで中心化され、及びQBP2 120は3MHzで中心化される。QBP1 115及びQBP2 120は各々、IQ信号と呼ばれる複素解析信号を出力する。IQ信号は、実部のInphase信号及び虚部のQuadrature信号を持つ。二乗の合計の平方根をとることにより、
Figure 2010508881
としてエコーのエンベロープが計算されることができる。ログ検出部125及びログ検出部130は、それぞれ、QBP1 115及びQBP2 120から複素信号を受信し、受信した複素信号のエンベロープを検出し、及びその後検出された結果の対数をとる。異なる周波数帯域フィルタからの検出信号を結合する方法は、「周波数合成」と呼ばれ、超音波業界において確立した技術であることに留意されたい。
平均化手段135は、ログ検出部125及びログ検出部130から記録されたエンベロープを受信する。ログ検出部125及びログ検出部130からの記録されたエンベロープは、2つの異なる周波数(例えば、それぞれ2及び3MHz)から得られる。スペックルは周波数の関数として変化する。一方、基礎をなす信号は同じままである。記録されたエンベロープが一緒に平均化されるとき、スペックルは平均される。その後、平均化された信号は、マルチレート低域フィルタ140へ入力され、SonoCT 145に出力される。スペックルは基礎をなす平均信号より急速に変化することができるので、このデータを低域フィルタリングすることは、スペックル変動を更に減らすだろう。マルチレート低域フィルタ140は高い空間周波数情報も減らす。それにより、この信号が大幅に縮小されることが可能にされる。これは、スキャンライン当たりのサンプル数を数千からわずか数百にまで減らす。サンプルが少なくなると、下流の処理動作におけるコンピュータの負担が減少する。
SonoCt145は、異なる表示角度から画像を得て、その後それらを単一画像に結合する合成イメージングデバイスである。スペックルパターンは、表示角度とともに変化する。マルチレート低域フィルタ140及びSonoCt 145に平均化手段135の出力を入力する目的は、超音波画像からスペックルを更に取り除くためである。その後、SonoCt 145の出力が、モニタといったディスプレイ150に入力される。
図8を参照すると、最適なスペックル追跡のために構成される従来技術の超音波イメージングシステム200の概略図が示される。イメージングシステム200は、超音波トランスデューサ(XD)105、スキャナ110、直角位相帯域フィルタ(QBP1)115、ログ検出部125、マルチレート低域フィルタ202、スペックル追跡装置205及びディスプレイ150を含む。
動作中、XD 105は、超音波を生み出し、スキャナ110にそれらを出力する。スキャナ110の出力は、QBP1 115に入力される。QBP1 115は、上記の通りIQ信号を出力する。ログ検出部125は、QBP1 115から複素信号を受信し、受信した複素信号のエンベロープを検出する。その後エンベロープは、マルチレート低域フィルタ202へ入力され、及びスペックル追跡装置205に出力される。最適なスペックル低減のため図7において使用されるマルチレート低域フィルタ140とは異なり、このマルチレート低域フィルタ202は、より少ないスムージング及び可能性としてより少ない減少を提供する。最適なスペックル追跡に対して、スペックルが強調されることが要求される。その結果、スペックルをマスキングするのに用いられる従来の技術は、ここでは不利益である。スペックル追跡装置205は、スペックルにおける変動を得るため異なる時点でスペックルを追跡する相互相関デバイスである。即ち、このデバイスは、標的(例えば組織)が運動するとき画像データを記録する。異なる時点でスペックルを相互相関させることにより、スペックル追跡装置は、組織変位、組織運動及び組織圧縮を計算することができる。その後スペックル追跡装置205の出力は、ディスプレイ150に入力される。
「スペックル追跡」に関しては公表された多数の方法が存在する(例えばChenらによる米国特許第5,876,342号)。追跡性能及び速度の両方のために通常使用される方法は、「正規化された相互相関」法である。それは、
Figure 2010508881
として規定される。
ここで、NCCは、正規化された相互相関関数である。
dx、dyは、スペックルがどれくらい移動したかを決定する検索空間である。
Figure 2010508881
は、関心領域(ROI)におけるx及びyに対して、合計が計算されることを示す。
u1は、時間1での画像である。
u2は、時間2での画像である。
この式は、以下のように適用される。
1.まず、第1の画像において追跡のために選択される関心領域(ROI)が特定される。複数のROIが選択されることができる点、及びそのすべてのピクセル(又は、3Dボリュームにおけるすべてのボクセル)が追跡のために選択されることができる点に留意されたい。これは、第1の画像u1においてROI並びにx及びyの範囲を規定する。
2.次に、後の時間において観察される画像u2において同じ大きさのROIを変位させるためにdx及びdyが変化させられる。
3.各dx及びdyに対して、正規化された相互相関(NCC)関数が評価される。
4.NCCのピーク最大値が観察されるまで、ステップ2及び3が繰り返される。NCC値が1.0であるとき、最大の相関を示す。このピーク値でのdx及びdyの値は、所望の組織がROIにおいてどれくらい移動したかを示す。
当業者にとって明らかであるように、(u1における)源ROI内又は(u2における)変位させられたROIにおいて、何のテクスチャ又はスペックル変動もない場合、NCCサーチアルゴリズムは失敗する。dx及びdyのすべての変位させられた値に対して1.0である相関値が観察され、それ故、ピークが特定されることができない。
本発明は、ある時間における組織の運動の一連の画像を得るために、画像強調超音波信号経路から得られるデータ及びスペックル強調超音波信号経路から得られるデータを結合する改良されたシステム及び方法を提供する。
図9を参照すると、最適なスペックル追跡のために構成される超音波イメージングシステム300の好ましい実施形態の概略図が示される。イメージングシステム300は、超音波トランスデューサ(XD)105、スキャナ/ビームフォーマ110、第1の直角位相帯域フィルタ(QBP1)115、第2の直角位相帯域フィルタ(QBP2)120、ログ検出部125、ログ検出部130、平均化手段135、第1のマルチレート低域フィルタ305、第2のマルチレート低域フィルタ310、スペックル追跡装置205、SonoCT145及びディスプレイ150を含む。
動作中、スキャナ110は、超音波トランスデューサXD105に電気信号を送信する。超音波トランスデューサは、この電気信号を音波に変換する。これらの音波は、体に伝播され、さまざまな生体構造において反射する。戻り音波エコーは、同じ超音波トランスデューサXD 105により電気信号に逆変換され、その後スキャナ110に送り返される。スキャナ110は、特定のスキャン方向及び深さからエコーを分離するためにその後これらの信号を処理する。これにより、それらの位置で生体構造が確認される。
スキャナ110の出力は、QBP1 115及びQBP2 120に入力される。ある実施形態において、QBP1 115は2MHzで中心化され、QBP2 120は3MHzで中心化される。QBP1 115及びQBP2 120はそれぞれ、IQ信号を出力する。IQ信号は、信号ノイズが除去される対象となる複素信号である。ログ検出部125及びログ検出部130はそれぞれ、QBP1 115及びQBP2 120から複素信号を受信し、受信した複素信号のエンベロープを検出する。上述したように、平均化手段135は、信号経路320を介してログ検出部125及びログ検出部130から信号エンベロープを受信し、画像のノイズ(スペックル)を平均化する。
その後、平均化された信号は、マルチレート低域フィルタ310に入力される。マルチレート低域フィルタ310の出力はSonoCT 145に入力される。SonoCTは、異なる表示角度からの画像を得て、その後それらを単一画像に結合する。その後SonoCt 145の出力は、ディスプレイ150に入力される。
ログ検出部125からの信号エンベロープも、信号経路315を介してマルチレート低域フィルタ305に入力される。マルチレート低域フィルタ305の出力は、異なる時点でのスペックルを追跡するスペックル追跡装置205に入力される。上記のように、異なる時点でのスペックルを相互相関させることにより、スペックル追跡装置は、組織変位、組織運動及び組織圧縮を計算することができる。その後スペックル追跡装置205の出力は、ディスプレイ150に入力される。
スペックル追跡装置205からのスペックルデータ及びSonoCT 145からの画像データが、ディスプレイ150によって同時に得られる。グラフとして又は第2の画像としてのいずれかとして、このスペックルデータ又は「機能情報」は、生体構造画像データと共に並べて表示されることができる。好ましい実施形態において、この機能情報は、例えば生体構造画像と異なるカラーを用いて、生体構造画像データ上にオーバレイ又は重畳されることができる。斯かる画像は、超音波業界ではしばしば「パラメトリック画像」と呼ばれる。
従って、パラメトリック画像を生み出すためスペックルデータは画像データに重畳されることができる。この重畳は、画像化される組織の運動が観察されることを可能にする。スペックルデータは、その値に基づきさまざまな色で表示されることができる。例えば、ある実施形態において、動く組織を示す「変化する」スペックルデータは緑で表示され、動かない組織を示す「変化しない」スペックルデータは灰色で表示される。有利には、「カラーの」スペックルデータが同時に得られた画像データ(組織)に重畳されるとき、移動する組織及び移動しない組織が観察されることができる。組織運動に加えて、血流量を観察するのに、得られたスペックルデータ及び画像データが使用されることができる。血流が起こるにつれ、時間にわたり組織は膨張及び収縮する。こうして、変化するスペックルデータがもたらされる。血流がない場合、得られるスペックルデータは変化しないだろう。
スペックル追跡装置205の直接出力は、検査される生体構造に関する運動及び変位情報を提供する。この情報は、多数の機能的属性を決定するのに使用されることができる。1つの例において、異なる構造の速度を決定するため、変位フィールドは時間に対して区別されることができる。別の例では、局所ひずみを計算するため、変位における空間差が使用することができる。斯かるひずみの測定は、健康であり及び収縮する心筋の部分と、虚血性で、死んでおり、及び収縮しない部分とを区別するために利用されることができる。更に別の例において、心臓の異なる部分が収縮するときを決定するタイミング解析のため、運動フィールドが使用されることができる。通常の健康な心臓では、左心室のすべての部分が、同時に収縮する傾向がある。しかしながら、逆同期収縮を持つ病気にかかった心臓では、心筋の異なる部分が異なる時間に収縮する。結果としてより効率的でないポンピングがもたらされる。
上記した全ての測定は、専用のハードウェアを用いて又はコンピュータで動作しているソフトウェアを用いて計算されることができる。また、リアルタイム(音波が得られている間)又は非リアルタイム(取得後)のいずれかの態様で斯かる測定を得ることも可能である。
上記した本発明のシステム及び方法は、胸部組織における腫瘍を検出するのに有益である。例えばマンモグラフィといった現在の方法は、例えば40〜50歳の女性といった腫瘍があまり高密度でない組織に囲まれる場合にだけ効果的である。周囲の組織の密度に関係なく、本発明は効果的に腫瘍を検出し、即ち血流又は組織運動のない領域を検出し、こうして20〜40歳の女性における腫瘍を発見することができる。
上記される本発明の方法は、心臓の破られた領域を探すのにも効果的である。斯かる領域は、損傷を受けており、血流が減らされる。従って、運動も減らされる。このことが、追跡及び観察されることができる。
更に、本発明は、電離性放射線又は放射性染料の導入を含む現在の診断方法より速く、より安全であり、及びより非侵襲的である。
図9において示される実施形態の重要な限定は、QBPフィルタのログ検出部処理バンク(例えば、QBPフィルタ115及びログ検出部125)の1つが、減じられたスペックル画像品質経路及び最適なスペックル追跡経路の両方により共有されることである。2つのQBPフィルタログ検出部バンクだけを必要とすることにより、この共有が結果として低コストの実現を生じさせることができる一方で、これは減少したスペックル画像品質経路及び最適なスペックル追跡経路両方の性能を損なう可能性がある。例えば、経路のうちの1つは基本周波数動作のために(QBPフィルタが、送信周波数の近くに中心周波数を持つ)、及び他の経路は組織調和イメージングのために(QBPフィルタが、送信周波数の2倍の中心周波数を持つ)構成されることが望ましい場合がある。
図10を参照すると、図9の性能制限に対処するため別の実施形態において、最適なスペックル追跡のために構成される超音波イメージングシステム400の概略図が示される。イメージングシステム400は、超音波トランスデューサ(XD)105、スキャナ/ビームフォーマ110、第1の直角位相帯域フィルタ(QBP1)115、第2の直角位相帯域フィルタ(QBP2)120、第3の直角位相帯域フィルタ(QBP3)405、ログ検出部125、ログ検出部130、ログ検出部410、平均化手段135、第1のマルチレート低域フィルタ305、第2のマルチレート低域フィルタ310、スペックル追跡装置205、SonoCT145及びディスプレイ150を含む。
動作中、XD 105は、超音波を電気信号に変換し、スキャナ110にそれらを出力する。スキャナ110の出力は、QBP1 115、QBP2 120及びQBP3 405に入力される。ある実施形態において、QBP1 115は2MHzで中心化され、QBP2 120は3MHzで中心化される。これは送信周波数が2.5MHzで中心化されるシナリオに関連する。QBP1 115及びQBP2 120は、送信周波数に近い基本周波数で周波数合成を実行することを試みる。これらの周波数は、最適な画像品質及び最適なスペックル低減のために選択されることができる。この同じシナリオにおいて、最適なスペックル追跡が組織調和イメージングを用いて実行されるべきであると結論されることができる(米国特許第5,879,303号を参照)。この場合、送信された音波の周波数の2倍の5MHzで中心化されるQBP3 405を持つことが適切である。QBP1 115、QBP2 120及びQBP3 405はそれぞれ、IQ信号を出力する。IQ信号は、信号ノイズが除去される対象となる複素信号である。ログ検出部125及びログ検出部130はそれぞれ、QBP1 115及びQBP2 120から複素信号を受信し、受信した複素信号のエンベロープを検出する。ログ検出部410は、QBP3 405から複素信号を受信し、受信した複素信号のエンベロープを検出する。
平均化手段135は、ログ検出部125及びログ検出部130から信号エンベロープを受信し、画像のスペックルを平均化する。その後、平均化された信号は、マルチレート低域フィルタ310に入力される。マルチレート低域フィルタ310の出力は、SonoCT 145に入力される。その後、SonoCt 145の出力は、モニタといったディスプレイ150に入力される。
同時に、ログ検出部410からの信号エンベロープが、マルチレート低域フィルタ305に入力される。マルチレート低域フィルタ305の出力が、スペックル追跡装置205に入力される。その後、スペックル追跡装置205の出力が、ディスプレイ150に入力される。
すべての実施形態及びブロック図は、同じスキャンラインを処理する異なる方法を説明する。その結果、第1の経路は最適な画像品質及びスペックル低減のために最適化され、第2の経路はスペックル追跡及びスペックル強調のために最適化される。スキャンラインは、体における特定の見通し線をなぞる(interrogating)音の単一のビームとして規定される。これは、軸方向の深さ(例えばmm単位)の大きさを持つ。このスキャンラインがどのようにシーケンス化されるかに基づき、異なるイメージングモード及び表示が得られることができる。ある実施形態において、スキャンラインは、同じ見通し線をなぞることができる(Mモードと呼ばれる)。第2の実施形態において、スキャンラインは、体における断層撮影スライスを通りシーケンス化することができる。これは、2D又はBモード動作と呼ばれる。更に別の実施形態において、スキャンラインは方位角方向(横方向)及び上昇方向の大きさの両方において変化することができる。これにより、ボリュームがスキャンされる(3D又は4Dイメージングと呼ばれる)。
また、当業者には明らかであるように、本発明は任意のタイプの超音波トランスデューサに適用できる。例えば、単一要素の機械的なトランスデューサ、フェーズドアレイ、リニア、湾曲リニアアレイ(CLA)、2Dマトリクスアレイ及びフェーズドアレイウォブラを含むが、これらに限定されるものではない。
本発明の更に別の実施形態において、単一の処理経路がライン毎に基づかれ変化するよう、並行処理経路がマルチプレクス化されることが想定される。その結果、1つの受信スキャンイベントの間、経路はスペックル追跡のために最適化され、及び同じ見通し線とすることができる別の受信スキャンイベントに対して、その経路は軽減されたスペックルを持つ最適な画像品質のために最適化される。
本発明の更に別の実施形態において、スペックル追跡のために使用される処理は、帯域フィルタのためにRFフィルタを使用すること、及びログ検出部を持たないことを含む。本発明の別の実施形態は、スペックル追跡経路において使用される遮断周波数以下の遮断周波数で、ポスト検出された低域フィルタを制限することにより、スペックルを減じることを含む。
本願明細書に説明されるもの変形、修正及び他の実現は、本発明の精神及び範囲から逸脱することなく当業者であれば思いつくことができる。従って、本発明は、上記の説明的な記載によってのみ規定されるものではない。

Claims (22)

  1. 超音波システムであって、
    人体に音波を送信し、前記音波のエコーを出力する手段と、
    少なくとも1つのスキャンラインデータを生み出すため、前記エコーを受信し、ビーム形成する手段と、
    生体構造情報を表示するため前記スキャンラインデータの1つを処理する第1の手段であって、スペックルを減らすことを含む処理をする第1の手段と、
    前記スキャンラインデータの1つを処理する第2の手段であって、スペックルを減らすことを含まない処理をする第2の手段と、
    1つのスキャンシーケンスの間、前記第1の手段を用いて処理される前記スキャンラインデータの1つ及び前記第2の手段を用いて処理される前記スキャンラインデータの1つを同時に得る手段とを有する、超音波システム。
  2. 前記第1の処理手段及び前記第2の処理手段が、前記同じスキャンラインを処理する、請求項1に記載のシステム。
  3. 前記第1の処理手段及び前記第2の処理手段が、前記1つのスキャンシーケンスとは異なるスキャンラインデータを処理する、請求項1に記載のシステム。
  4. 前記第1の処理手段が、RF帯域フィルタを含む、請求項1に記載のシステム。
  5. 前記第2の処理手段が、RF帯域フィルタを含む、請求項1に記載のシステム。
  6. 前記第1の処理手段が、前記エコーのエンベロープを検出する手段を含む、請求項1に記載のシステム。
  7. 前記第1の処理手段が、前記検出されたエンベロープの対数をとる手段を含む、請求項6に記載のシステム。
  8. スペックルを減らすことが、周波数合成を用いて実行される、請求項1に記載のシステム。
  9. 2つ又はこれ以上のフィルタバンクを更に有し、各フィルタバンクが、検出器と、各フィルタバンクに対する一意なレスポンスを持つ帯域フィルタとを有する、請求項8に記載のシステム。
  10. スペックルを減らすことが、空間合成により実行される、請求項1に記載のシステム。
  11. 前記スキャンシーケンスが、単一の見通し線、平面及びボリュームの少なくとも1つをなぞる、請求項1に記載のシステム。
  12. 前記音波を送信する手段が、フェーズドアレイ、リニア、湾曲リニアアレイ、機械的なウォブラ及び3Dウォブラからなるグループから選択される、請求項1に記載のシステム。
  13. 前記第1の処理手段及び前記第2の処理手段が、RF帯域フィルタを有する、請求項1に記載のシステム。
  14. 前記第1の処理手段が、専用のハードウェアを用いて、又はCPUにおいて動作するソフトウェアを用いて実現されることができる、請求項1に記載のシステム。
  15. 前記第2の処理手段が、専用のハードウェアを用いて、又はCPUにおいて動作するソフトウェアを用いて実現されることができる、請求項1に記載のシステム。
  16. 前記第2の処理手段を用いて処理される前記スキャンラインデータの1つが、パラメトリックイメージング表示技術を用いて生体構造情報に重畳されることができる、請求項1に記載のシステム。
  17. 前記第2の処理手段を用いて処理される前記スキャンラインデータの1つが、以前のスキャンシーケンスから得られるデータと相互相関される、請求項1に記載のシステム。
  18. 前記相互相関させられたデータが、ひずみ、ひずみ率、エラストグラフィ、壁濃厚化及び収縮タイミングに関して使用されることができる、請求項17に記載のシステム。
  19. ポスト検出された低域フィルタをスペックル追跡経路において使用される遮断周波数以下の遮断周波数に制限することにより、スペックを減らすことが実行される、請求項1に記載のシステム。
  20. 時間に対する組織の空間変位を決定するため、前記1つのスキャンシーケンスが繰り返されることができる、請求項1に記載のシステム。
  21. スペックル追跡を実行する方法において、
    人体に音波を送信し、前記音波のエコーを出力するステップと、
    少なくとも1つのスキャンラインデータを生み出すため前記エコーを受信し、ビーム形成するステップと、
    生体構造情報を表示するため前記スキャンラインデータの1つを処理するステップであって、前記処理がスペックルを減らすことを含む、ステップと、
    前記スキャンラインデータの1つを追加的に処理するステップであって、前記追加的な処理が、スペックルを減らすことを含まない、ステップと、
    1つのスキャンシーケンスの間、生体構造情報を表示するために処理された前記スキャンラインデータの1つと追加的な処理を用いて処理された前記スキャンラインデータの1つとを同時に得るステップとを有する、方法。
  22. スペックル追跡を実行するコンピュータで動作するコンピュータ可読プログラムコードを持つコンピュータ可読媒体であって、
    人体に音波を送信し、前記音波のエコーを出力するステップと、
    少なくとも1つのスキャンラインデータを生み出すため前記エコーを受信し、ビーム形成するステップと、
    生体構造情報を表示するため前記スキャンラインデータの1つを処理するステップであって、前記処理がスペックルを減らすことを含む、ステップと、
    前記スキャンラインデータの1つを追加的に処理するステップであって、前記追加的な処理が、スペックルを減らすことを含まない、ステップと、
    1つのスキャンシーケンスの間、生体構造情報を表示するために処理された前記スキャンラインデータの1つと追加的な処理を用いて処理された前記スキャンラインデータの1つとを同時に得るステップとを有する、コンピュータ可読媒体。
JP2009535182A 2006-11-03 2007-11-02 最適なスペックル追跡のための二重経路処理 Expired - Fee Related JP5627890B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US86425906P 2006-11-03 2006-11-03
US60/864,259 2006-11-03
PCT/IB2007/054466 WO2008053457A2 (en) 2006-11-03 2007-11-02 Dual path processing for optimal speckle tracking

Publications (2)

Publication Number Publication Date
JP2010508881A true JP2010508881A (ja) 2010-03-25
JP5627890B2 JP5627890B2 (ja) 2014-11-19

Family

ID=39301157

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009535182A Expired - Fee Related JP5627890B2 (ja) 2006-11-03 2007-11-02 最適なスペックル追跡のための二重経路処理

Country Status (6)

Country Link
US (1) US20100004540A1 (ja)
EP (1) EP2082261A2 (ja)
JP (1) JP5627890B2 (ja)
KR (1) KR101407425B1 (ja)
CN (1) CN101563626A (ja)
WO (1) WO2008053457A2 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011104361A (ja) * 2009-11-19 2011-06-02 Medison Co Ltd 超音波空間合成映像を提供する超音波システムおよび方法
JP2012176232A (ja) * 2011-02-04 2012-09-13 Toshiba Corp 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム
JP2014528321A (ja) * 2011-10-03 2014-10-27 サーフ テクノロジー アクティーゼルスカブSurf Technology As デュアルバンドパルス合成物による非線形イメージング
US9345446B2 (en) 2011-12-08 2016-05-24 Kabushiki Kaisha Toshiba Ultrasonic diagnosis apparatus, image processing apparatus, and image processing method

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8795181B2 (en) * 2008-11-25 2014-08-05 Mayo Foundation For Medical Education And Research System and method for analyzing carpal tunnel using ultrasound imaging
EP2385391A3 (en) 2010-05-04 2012-08-01 Sony Corporation Active imaging device and method for speckle noise reduction
US10231712B2 (en) 2010-06-09 2019-03-19 Regents Of The University Of Minnesota Dual mode ultrasound transducer (DMUT) system and method for controlling delivery of ultrasound therapy
AU2012242639B2 (en) * 2011-04-14 2016-09-01 Regents Of The University Of Minnesota Vascular characterization using ultrasound imaging
TWI446897B (zh) * 2011-08-19 2014-08-01 Ind Tech Res Inst 超音波影像對齊裝置及其方法
TWI487402B (zh) * 2012-08-10 2015-06-01 Mstar Semiconductor Inc 可用於一無線通訊系統的搜尋方法
EP2944266A4 (en) * 2013-01-11 2016-08-31 Hitachi Aloka Medical Ltd ULTRASOUND IMAGING DEVICE
US10034658B2 (en) 2013-03-05 2018-07-31 Koninklijke Philips N.V. Consistent sequential ultrasound acquisitions for intra-cranial monitoring
WO2015013196A2 (en) 2013-07-23 2015-01-29 Regents Of The University Of Minnesota Ultrasound image formation and/or reconstruction using multiple frequency waveforms
RU2695475C2 (ru) * 2015-01-29 2019-07-23 Конинклейке Филипс Н.В. Оценка инфаркта миокарда с помощью ультразвуковой визуализации деформаций в режиме реального времени
KR102387708B1 (ko) 2015-01-30 2022-04-19 삼성메디슨 주식회사 향상된 hprf 도플러 영상을 위한 가이드를 제공하는 방법 및 초음파 시스템
CN104586433B (zh) * 2015-02-02 2016-08-24 声泰特(成都)科技有限公司 基于变频的基波/谐波融合与空间复合相结合的成像方法
US10675007B2 (en) * 2016-04-19 2020-06-09 Siemens Medical Solutions Usa, Inc. Frequency compounding in elasticity imaging
GB201614950D0 (en) * 2016-09-02 2016-10-19 Ntnu Tech Transfer As Enhanced-resolution ultrasound imaging of fluid paths
WO2018051265A1 (en) 2016-09-15 2018-03-22 Koninklijke Philips N.V. Ultrasonic elastographic pre-load measurement and display
US11458337B2 (en) 2017-11-28 2022-10-04 Regents Of The University Of Minnesota Adaptive refocusing of ultrasound transducer arrays using image data
US11596812B2 (en) 2018-04-06 2023-03-07 Regents Of The University Of Minnesota Wearable transcranial dual-mode ultrasound transducers for neuromodulation
WO2021216723A1 (en) 2020-04-22 2021-10-28 Bfly Operations, Inc. Methods and apparatuses for beamforming in ultrasound systems
WO2024020093A2 (en) * 2022-07-20 2024-01-25 The Johns Hopkins University Methods, systems, and computer readable media for storing and processing ultrasound audio data

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004073620A (ja) * 2002-08-21 2004-03-11 Toshiba Medical System Co Ltd 超音波診断装置
JP2006000618A (ja) * 2004-05-20 2006-01-05 Fuji Photo Film Co Ltd 超音波撮像装置、超音波画像処理方法、及び、超音波画像処理プログラム
US20060241431A1 (en) * 2005-03-31 2006-10-26 Naohisa Kamiyama Ultrasound diagnostic apparatus and ultrasound image processing program

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2005208A (en) * 1934-08-20 1935-06-18 Schuchat Jonas Wrapping
EP0321717B1 (en) * 1987-11-30 1995-07-05 Aloka Co. Ltd. Ultrasonic speckle velocity measurement method and apparatus
US5474070A (en) * 1989-11-17 1995-12-12 The Board Of Regents Of The University Of Texas System Method and apparatus for elastographic measurement and imaging
US5879303A (en) * 1996-09-27 1999-03-09 Atl Ultrasound Ultrasonic diagnostic imaging of response frequency differing from transmit frequency
US5735797A (en) * 1996-12-30 1998-04-07 General Electric Company Method and apparatus for combining topographic flow power imagery with a B-mode anatomical imagery
US5876342A (en) * 1997-06-30 1999-03-02 Siemens Medical Systems, Inc. System and method for 3-D ultrasound imaging and motion estimation
US6213946B1 (en) * 1998-12-24 2001-04-10 Agilent Technologies, Inc. Methods and apparatus for speckle reduction by orthogonal pulse compounding in medical ultrasound imaging
US6142942A (en) * 1999-03-22 2000-11-07 Agilent Technologies, Inc. Ultrasound imaging system and method employing an adaptive filter
US6464637B1 (en) * 2000-06-23 2002-10-15 Koninklijke Philips Electronics N.V. Automatic flow angle correction by ultrasonic vector
CA2445137A1 (en) * 2002-10-15 2004-04-15 Matsushita Electric Industrial Co., Ltd. Image processing apparatus, method and program
US20050053305A1 (en) * 2003-09-10 2005-03-10 Yadong Li Systems and methods for implementing a speckle reduction filter
US20050288589A1 (en) * 2004-06-25 2005-12-29 Siemens Medical Solutions Usa, Inc. Surface model parametric ultrasound imaging
US7678050B2 (en) * 2004-08-24 2010-03-16 General Electric Company Method and apparatus for detecting cardiac events

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004073620A (ja) * 2002-08-21 2004-03-11 Toshiba Medical System Co Ltd 超音波診断装置
JP2006000618A (ja) * 2004-05-20 2006-01-05 Fuji Photo Film Co Ltd 超音波撮像装置、超音波画像処理方法、及び、超音波画像処理プログラム
US20060241431A1 (en) * 2005-03-31 2006-10-26 Naohisa Kamiyama Ultrasound diagnostic apparatus and ultrasound image processing program

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
JPN7013000243; Abd El-Monem El-Sharkawy et al.: 'Improved ultrasound speckle motion tracking using nonlinear diffusion filtering' SPIE Proceedings Vol. 4325, 20010330, pp.453-461 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011104361A (ja) * 2009-11-19 2011-06-02 Medison Co Ltd 超音波空間合成映像を提供する超音波システムおよび方法
US8727990B2 (en) 2009-11-19 2014-05-20 Samsung Medison Co., Ltd. Providing an ultrasound spatial compound image in an ultrasound system
JP2012176232A (ja) * 2011-02-04 2012-09-13 Toshiba Corp 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム
JP2014528321A (ja) * 2011-10-03 2014-10-27 サーフ テクノロジー アクティーゼルスカブSurf Technology As デュアルバンドパルス合成物による非線形イメージング
US9345446B2 (en) 2011-12-08 2016-05-24 Kabushiki Kaisha Toshiba Ultrasonic diagnosis apparatus, image processing apparatus, and image processing method

Also Published As

Publication number Publication date
KR20090084840A (ko) 2009-08-05
WO2008053457A3 (en) 2008-07-03
KR101407425B1 (ko) 2014-06-17
JP5627890B2 (ja) 2014-11-19
WO2008053457A2 (en) 2008-05-08
CN101563626A (zh) 2009-10-21
EP2082261A2 (en) 2009-07-29
US20100004540A1 (en) 2010-01-07

Similar Documents

Publication Publication Date Title
JP5627890B2 (ja) 最適なスペックル追跡のための二重経路処理
Jensen et al. Ultrasound vector flow imaging—Part II: Parallel systems
JP6063553B2 (ja) 超音波イメージング方法及び超音波イメージング装置
CN104272134B (zh) 超声成像系统中的杂波抑制
US11612381B2 (en) Method for tissue characterization by ultrasound wave attenuation measurements and ultrasound system for tissue characterization
JP6063552B2 (ja) 超音波イメージング方法及び超音波イメージング装置
JP2012531933A (ja) 非ecgの生体同期収集による三次元胎児心臓撮像方法を用いた超音波画像診断システム及び方法
CN103202714B (zh) 超声波诊断装置、医用图像处理装置以及医用图像处理方法
US9855025B2 (en) Ultrasonic diagnostic apparatus and ultrasonic image processing apparatus
JP4320392B2 (ja) 高歪みレート除去フィルタ処理のための方法及び装置
JP4870449B2 (ja) 超音波診断装置及び超音波画像処理方法
JP2007020908A (ja) 超音波診断装置及び超音波診断装置の制御プログラム
Hemmsen et al. Tissue harmonic synthetic aperture ultrasound imaging
JP2007518512A (ja) 心筋灌流を表示するための画像分割
WO2019087741A1 (ja) 超音波診断装置、および、生体組織の物性評価方法
CN112842382B (zh) 用于对信道数据进行流处理以应用非线性波束形成的方法和系统
CN112862924B (zh) 多模态成像中图像重建方法、装置和多模态成像技术系统
Byram et al. 3-D phantom and in vivo cardiac speckle tracking using a matrix array and raw echo data
JP4921816B2 (ja) 超音波診断装置及びその制御プログラム
JP2002028160A (ja) 超音波診断装置
Hyun et al. A GPU-based real-time spatial coherence imaging system
JP5299961B2 (ja) 超音波診断装置、画像処理装置及び超音波診断装置の制御プログラム
KR101956460B1 (ko) 초음파 의료영상 장치를 이용한 미세석회화 조직 검출 방법 및 그 초음파 의료영상 장치
KR101610877B1 (ko) 공간 일관성 기초 초음파 신호 처리 모듈 및 그에 의한 초음파 신호 처리 방법
Onyia Displacement Data Processing for ARFI Imaging

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20101029

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20120928

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130124

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20130423

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20130501

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130529

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20131121

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140128

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: 20140902

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20141001

R150 Certificate of patent or registration of utility model

Ref document number: 5627890

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees