JP4306864B2 - 3次元イメージング・システム並びに走査平面の動きを追跡する方法及び装置 - Google Patents
3次元イメージング・システム並びに走査平面の動きを追跡する方法及び装置 Download PDFInfo
- Publication number
- JP4306864B2 JP4306864B2 JP07470999A JP7470999A JP4306864B2 JP 4306864 B2 JP4306864 B2 JP 4306864B2 JP 07470999 A JP07470999 A JP 07470999A JP 7470999 A JP7470999 A JP 7470999A JP 4306864 B2 JP4306864 B2 JP 4306864B2
- Authority
- JP
- Japan
- Prior art keywords
- correlation
- image frames
- interval
- degree
- data
- 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.)
- Expired - Lifetime
Links
- 230000033001 locomotion Effects 0.000 title claims description 14
- 238000003384 imaging method Methods 0.000 title claims description 11
- 238000000034 method Methods 0.000 title description 29
- 238000012935 Averaging Methods 0.000 claims description 13
- 238000002604 ultrasonography Methods 0.000 claims description 12
- 238000012216 screening Methods 0.000 claims description 7
- 238000002592 echocardiography Methods 0.000 claims description 5
- 230000000149 penetrating effect Effects 0.000 claims 1
- 230000006835 compression Effects 0.000 description 9
- 238000007906 compression Methods 0.000 description 9
- 239000000523 sample Substances 0.000 description 9
- 230000006870 function Effects 0.000 description 8
- 238000012805 post-processing Methods 0.000 description 8
- 230000000875 corresponding effect Effects 0.000 description 7
- 210000003484 anatomy Anatomy 0.000 description 5
- 238000005314 correlation function Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000013459 approach Methods 0.000 description 4
- 238000003491 array Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 230000002596 correlated effect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000017531 blood circulation Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012285 ultrasound imaging Methods 0.000 description 3
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 238000005266 casting Methods 0.000 description 2
- 230000006837 decompression Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 239000013598 vector Substances 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000013329 compounding Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/02—Affine transformations
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03B—APPARATUS OR ARRANGEMENTS FOR TAKING PHOTOGRAPHS OR FOR PROJECTING OR VIEWING THEM; APPARATUS OR ARRANGEMENTS EMPLOYING ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ACCESSORIES THEREFOR
- G03B42/00—Obtaining records using waves other than optical waves; Visualisation of such records by using optical means
- G03B42/06—Obtaining records using waves other than optical waves; Visualisation of such records by using optical means using ultrasonic, sonic or infrasonic waves
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/32—Determination of transform parameters for the alignment of images, i.e. image registration using correlation-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/97—Determining parameters from multiple pictures
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/13—Tomography
- A61B8/14—Echo-tomography
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8934—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a dynamic transducer configuration
- G01S15/8936—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a dynamic transducer configuration using transducers mounted for mechanical movement in three dimensions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S128/00—Surgery
- Y10S128/916—Ultrasound 3-D imaging
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Description
【産業上の利用分野】
本発明は一般的には、医療診断を目的とした人体の解剖学的構造の超音波イメージングに関する。具体的には、本発明は、人体内の走査される領域によって反射された超音波エコーの強度を検出することにより人体の解剖学的構造をイメージングする方法に関する。
【0002】
【従来の技術】
従来の超音波スキャナは、ピクセル(画素)の輝度がエコー反射の強度に基づくようにして組織の2次元Bモード画像を形成する。従来のBモードでの基本的な信号処理チェーンを図1に示す。超音波トランスデューサ・アレイ2が駆動されて、走査線に沿って音波バーストを送信する。反射RF信号は、トランスデューサ素子によって検出され、次いで、ビームフォーマ4によって受信ビームとして形成される。各々の走査線ごとのビームフォーマの出力データ(I/Q又はRF)は、等化フィルタリング、包絡線検出及び対数圧縮を含むBモード処理チェーン6に供給される。走査形状に応じて、数百までのベクトルを用いて、単一の音波画像フレームを形成することができる。1つの超音波フレームから次の音波フレームへの時間的な移行を平滑化するために、走査変換の前に、何らかの音波フレーム平均化(ブロック8で示す)を行ってもよい。フレーム平均は、FIRフィルタ又はIIRフィルタによって実現することができる。一般に、圧縮された画像は、R−θフォーマット(セクタ走査の場合)を有しており、R−θフォーマットは、走査変換器10によって変換されて、ビデオ表示用のX−Yフォーマットとなる。システムによっては、走査変換の前に音波フレームに対してフレーム平均化を行うのではなく、ビデオX−Yデータに対して行ってもよく(破線のブロック12で示す)、所与のビデオ表示フレーム・レート(典型的には、30Hz)を達成するために、重複したビデオ・フレームが音波フレームの間に挿入されることもある。ビデオ・フレームはビデオ・プロセッサ14に伝送され、ビデオ・プロセッサ14は基本的にはビデオ・データをビデオ表示のための表示用グレイ・マップにマッピングする。
【0003】
システムの制御は、ホスト・コンピュータ20に集中化されており、コンピュータ20は、オペレータ・インタフェイス22(例えば、キーボード)を介してオペレータの入力を受け付け、次いで、様々なサブシステムを制御する。(図1には、画像データの転送経路のみが示されている。)Bモード・イメージングの際には、シネ・メモリ16において、最新の画像の長い系列(シーケンス)が記憶されていると共に自動的に連続して更新されている。R−θ音波画像を保存するように設計されているシステムもあり(このデータ経路を図1の破線で示す)、X−Yビデオ画像を記憶するシステムもある。シネ・メモリ16に記憶されている画像ループはトラック・ボール制御を介して通覧することができ、画像ループの一部分は選択してハード・ディスクに記憶させることもできる。フリー・ハンド式3次元イメージング性能を持つ超音波スキャナでは、シネ・メモリ16に記憶されている選択された画像系列が、3次元再構成のためにホスト・コンピュータ20に転送される。この結果は、シネ・メモリの他の部分に書き戻され、ここから、ビデオ・プロセッサ14を介して表示システム18に伝送される。
【0004】
図2を見ると、走査変換器10は、音響線メモリ24とX−Yメモリ26とを含んでいる。音響線メモリ24内に極座標(R−θ)セクタ・フォーマットで記憶されているBモード・データが、適当にスケーリングされたデカルト座標の強度データに変換されて、X−Yメモリ26に記憶される。Bモード・データの相次ぐ多数のフレームが、先入れ先出し方式でシネ・メモリ16に記憶されている。シネ・メモリは、バックグラウンドで稼働する循環的画像バッファのようなものであり、実時間でユーザに対して表示されている画像データを連続的に取り込む。ユーザがシステムをフリーズさせると、ユーザは、シネ・メモリに以前に取り込まれている画像データを見ることができるようになる。
【0005】
ホスト・コンピュータ20は、中央処理装置(CPU)28とランダム・アクセス・メモリ(RAM)30とを含んでいる。CPU28は読み出し専用メモリ(ROM)を内蔵しており、このメモリは、収集された強度データを異なる角度から見た多数の3次元投影画像に変換するのに用いられるルーチンを記憶している。CPU28は、システム制御バス32を介してX−Yメモリ26及びシネ・メモリ16を制御する。具体的には、CPU28は、走査変換器10の音響線メモリ24又はX−Yメモリ26からビデオ・プロセッサ14及びシネ・メモリ16へのデータの流れ、並びにシネ・メモリからビデオ・プロセッサ14及びCPU28自体へのデータの流れを制御する。各々のフレームのイメージング・データは、被検体を貫通する多数の走査又はスライスのうちの1つを表すものであり、音響線メモリ24、X−Yメモリ26及びビデオ・プロセッサ14に順次記憶される。それと並行して、音響線メモリ又はX−Yメモリのいずれかからの画像フレームが、シネ・メモリ16に記憶される。走査された物体領域を表す一群のフレームが、シネ・メモリ16の部分16Aに記憶される。
【0006】
2次元超音波画像は、しばしば解釈が困難なものであり、これは、走査中の解剖学的構造の2次元表現を観測者が思い描くことができないことに起因している。しかしながら、超音波プローブが関心領域にわたって掃引されて、2次元画像が3次元領域を形成するように蓄積されると、解剖学的構造は、熟練した観測者及び熟練していない観測者の両者にとって、遥かに思い描き易いものになる。
【0007】
3次元画像を形成するために、CPU28は、米国特許第5,226,113号又は同第5,485,842号に開示されているようなレイ・キャスティング(ray casting )アルゴリズムを用いた一連の変換を実行することができる。レイ・キャスティング法は、シネ・メモリ16の部分16Aから検索された関心のあるソース・データ・ボリュームのデータに適用される。続いての変換は、ある角度範囲、例えば、+90°から−90°までの角度範囲内のある角度増分、例えば、10°の間隔で行われる最大投影、最小投影、合成投影、表面投影又は平均投影等の様々な投影手法を含み得る。投影後の画像の各々のピクセルは、所与の画像平面への投影によって得られる変換後のデータを含んでいる。加えて、シネ・メモリがオペレータによってフリーズされた時点で、CPU28は選択により、シネ・メモリ16の部分16B内の多数の相次ぐアドレスにX−Yメモリ28からの最後のフレームを記憶する。最初に投影されたビュー角度についての投影画像データが、シネ・メモリの部分16Bの最初のアドレスに書き込まれているので、関心領域の投影画像データは、背景フレームにスーパインポーズ(重ね表示)される。この過程は、すべての投影画像がシネ・メモリの部分16Bに記憶されるまで各々の角度増分ごとに繰り返される。このとき、各々の投影画像フレームは、変換後の強度データを含んでいるデータ関心領域と、随意選択により関心領域を包囲する背景の周縁領域であって、変換後の強度データによって上書きされることのない背景強度データから成っている背景の周縁領域とで構成されている。背景画像は、各々の表示されている投影が、どこから見たものであるかをより明確にする。次いで、オペレータは、投影画像から任意の画像を選択して表示することができる。加えて、これら一連の投影画像を表示モニタ上で再生して、物体空間が観察者の前で回転しているかのように表現することもできる。
【0008】
多数行型(multi-row )トランスデューサ・アレイには、いわゆる「1.25D」アレイ及び「1.5D」アレイを含めて様々な種類があり、従来の単一行型(「1D」)アレイの限定された仰角(elevation )性能を向上させるために開発されている。ここで用いる場合には、上の各用語は次の意味を有するものとする。即ち、1Dでは、仰角開口が固定されており、焦点が固定されたレンジ(距離)にある。1.25Dでは、仰角開口が可変であるが、焦点合わせが依然静的である。1.5Dでは、仰角開口、シェーディング及び焦点合わせが動的に可変であるが、但し、アレイの中心線に関して対称である。
【0009】
フリー・ハンド式3次元超音波走査の場合には、トランスデューサ・アレイ(1D〜1.5D)を仰角方向に平行移動させて、関心のある解剖学的構造の全体にわたる実質的に平行な1組の画像平面を収集する。これらの画像は、シネ・メモリに記憶されると共に、後には3次元再構成のためにシステム・コンピュータによって検索することができる。画像フレーム相互の間の間隔が既知であるならば、3次元領域は、平面外寸法と走査平面寸法との間の正しいアスペクト比によって再構成され得る。しかしながら、スライス相互の間の間隔の推定値が貧弱なものであると、3次元物体の形状が大幅に歪む可能性がある。
【0010】
従来技術では、プローブの仰角方向の動き(モーション)をそれぞれ制御する又は追跡するために、さまざまなモーション制御方法及び位置検知方法が提案されている。しかしながら、これらのシステムは、臨床環境で用いるにはしばしば経費が高く、面倒なものである。従って、3次元画像を仰角方向に良好な分解能で再構成するためには、相次ぐ画像フレーム相互の間のスペックル相関外れ(decorrelation )の度合から走査平面の変位を直接的に推定し得るようにすることが極めて望ましい。
【0011】
国際特許WO97/00482号において、Fowlkes 等は、画像フレーム相互の間の相関の計算に基づく走査平面の動きを追跡する方法を提案している。この相関方法は、血流を監視するのに用いられる相関外れ手法を改造したものであると述べられている。このような従来技術を再検討すると、以下のような大きく2つのアプローチがあることがわかる。
【0012】
(1) IEEE Trans. Ultrasonics, Ferroelec. and Freq. Control 誌、第UFFC−33巻(1986年)、第257頁〜第264頁所載のTrahey等の論文「横方向の開口の平行移動によるスペックル・パターン相関:実験結果と空間合成における意味(Speckle pattern correlation with lateral aperture translation: experimental results and implications for spatial compounding)」には、超音波画像内の強度の完全相関関数を用いた最初の研究が報告されている。この手法は、圧縮前のRF画像データ又は検出された画像データを用いており、このことは、相関関数が全エコー強度によって正規化されている(即ち、システム利得依存性である)という事実から明らかである。Int. J. Imaging Syst. Technol.誌、第8巻(1997年)、第38頁〜第44頁所載のChen等の論文「スペックル相関外れを用いた走査平面の動きの決定:理論的考察と初期試験(Determination of scan-plane motion using speckle decorrelation: theoretical considerations and intial test)」では、3次元掃引距離の推定のためにこの相関方法を用いたファントムの研究が報告されている。
【0013】
(2) IEEE Trans. Biomed. Eng.誌、第38巻(1991年)、第280頁〜第286頁所載のBohs等の論文「血流及び組織の動きの角度非依存的超音波イメージングの新しい方法(A novel method for angle independent ultrasonic imaging of blood flow and tissue motion )」には、圧縮後の超音波画像に適用されるべき一層単純な相関方法が提案されている。この方法は、相関させられる2つのカーネルの対応するピクセル相互の間での絶対差(差の絶対値)の和を算出することに基づいているので、SAD(sum of absolute differences )法と呼ばれている。計算上効率的なこの方法は、Trahey等の完全相関関数とほとんど同じ程度に良好に機能することが判明している。IEEE Trans. Ultrasonics, Ferroelec. and Freq. Control 誌、第UFFC−41巻(1994年)、第96頁〜第104頁所載のShehada 等の論文「全血の非ニュートン特性を調べるための超音波方法(Ultrasound methods for investigating the non-Newtonian characteristics of whole blood )」にもまた、超音波画像のSAD相関に基づく流れの測定研究が報告されている。流れの測定のためには、SAD法は基本的に、SADを最小化する走査平面内での変位ベクトルを求めることから成っているので、SADの相対的な変化のみが関係してくる。しかしながら、一般的には、所与のカーネル寸法について、SADはダイナミック・レンジの設定及び後処理フィルタリングによって大幅に変化する可能性がある。
【0014】
【発明が解決しようとする課題】
3次元フリー・ハンド式走査性能を持つ超音波スキャナでは、シネ・メモリに記憶されている画像は典型的には、表示(典型的には、8ビットの振幅表示)のために対数圧縮又は他の何らかの高次非線形圧縮を既に経ている。これらの画像はまた、平滑化又はエッジ強調等のための何らかの後処理フィルタを経てきている。圧縮及びフィルタリングの操作はしばしば非可逆的なものであるので、近似的な「逆圧縮(decompression )」であっても、行おうとすれば画像に量子化雑音を導入する可能性がある。この理由から、上述した1番目の方法は、圧縮前の画像向けに設計されており、掃引速度の推定に最も適しているとは言えない。
【0015】
SAD法は、圧縮後の画像に適用されるものであり、計算速度についての利点がある。しかしながら、3次元再構成については、掃引距離を推定するために、フレームからフレームへの実際の相関外れを定量化する必要がある。SADのみを用いたとすると、表示ダイナミック・レンジ(これは、深さ依存性であり得る)、フィルタ、カーネル寸法及びカーネルの深さ位置のすべての可能な組み合わせについて較正を行う必要が出てくる。
【0016】
従って、異なるダイナミック・レンジ設定及び後処理フィルタに適応する新たな相関インデクスが必要とされる。
【0017】
【課題を解決するための手段】
本発明は、適応スペックル相関を用いてフリー・ハンド式3次元超音波走査における走査平面の動き(モーション)を追跡する方法及び装置である。この方法は、異なる表示ダイナミック・レンジ及び後処理フィルタに適応する相関インデクスを採用する。1つの好ましい態様による装置は、相関計算のために各々の画像フレーム内でカーネルを選択する手段と、重複した画像フレームを除去する手段と、相次ぐ画像フレーム相互の間の相関度を測定する手段と、手ぶれ及び他のアーティファクトに関連している可能性のある相関推定値を除去する手段と、平均相関推定値に基づいて平均フレーム間(即ち、スライス間)間隔を算出する手段とを含んでいる。これらの各手段は、シネ・メモリに結合されているホスト・コンピュータに組み込まれている。画像に基づくこのモーション追跡方法の主な長所は、外部の位置検知装置を一切用いずに、良好な形状忠実度で3次元再構成を行うことが可能になることである。
【0018】
【発明の実施の形態】
図3は、ユーザによって選択された1つの画像フレーム36(この例ではセクタ走査)内の典型的な関心領域(ROI)ボックス34を示している。フリー・ハンド式走査によって形成されたN個のフレームがシネ・メモリに記憶されていると仮定する。図4に、3次元再構成(これは、シネ・メモリから検索されたデータに基づいてホスト・コンピュータによって実行される)のために平均フレーム間間隔dを推定する処理工程を示す。先ず、本発明における相関法が拡散散乱媒体によって生ずる純粋なスペックルの統計に基づいているので、比較的純粋なスペックル・パターン(巨視的構造を有さないもの)を示すROI34内のカーネル38(この例ではM×Mピクセルとしているが、一般的には、カーネルは正方形でなくてもよい)を同定(identify)しなければならない(工程40)。カーネル38は、1つ又はそれ以上の画像フレームの目視判定に基づいて、手動で(例えば、トラックボールを用いて)選択することができる。代替的には、何らかの自動化された方法を用いて、そのピクセルの振幅ヒストグラムが理論的なスペックル分布に一致しているようなカーネルを探してもよい。例えば、このような試験は、通常のダイナミック・レンジ設定に対するヒストグラム幅の測度に基づいたものであり得る。その平均ピクセル値が過度に低いカーネル(信号が存在しないもの)又はその分散が過度に大きいカーネル(均一でないもの)は排除すべきである。図3に示すように、ユーザが関心のある構造をROIの中央に配置する傾向にあるとすると、試験に好適な初期カーネル38は、ROI34の四隅の任意位置にある。
【0019】
相次ぐ画像フレーム内のカーネル相互の間の相関の計算に移る前に、N個のソース・フレームに重複した画像フレームが存在しているならば、これらの画像フレームを同定して除去する(図4の工程42)。ビデオ・モニタの表示速度(典型的には、30Hz)に合わせるために、重複したフレームが音波画像フレーム相互の間に挿入されていることがある。重複したフレームの対が厳密に同一であれば、SADが実質的にゼロに等しくなるという規準に基づいて、これらの対を容易に検出することができる。重複したフレームが、例えばフレーム平均化効果によって、厳密には同一でない場合は、何らかの閾値法を用いて近似的に同一なフレームを検出すればよい。例えば、新たな画像フレームのカーネル内のピクセルの或る特定のパーセント(例えば、25%)よりも多いピクセルが或る特定の値よりも大きく前のフレームのカーネル内のピクセルと異なっていた場合、この新たな画像フレームを新たな音波フレームとする。そうでない場合は、この新たな画像フレームを、除去すべき重複した画像フレームと見做す。代替的には、重複したフレームをスクリーニング(篩い分け)して除去するときの信頼性を向上させるために、(N−1)個のソース・フレームの対すべてについてのSAD値の集団的統計を考察して、通常の統計的な偏差より外側に位置するSAD値を有するフレームを除去する必要があろう。
【0020】
重複している可能性のあるフレームをスクリーニングして除去した後、次の工程(図4の工程44)は、音波フレームの残りのセット内のすべての隣接するフレームの対について相関インデクス(CI)を算出することである。本発明では、異なるカーネル寸法、表示ダイナミック・レンジ及び後処理フィルタに適応し得る正規化されたSADと見做すことのできる相関インデクスを用いる。このインデクスは、画像の逆圧縮を必要とし且つ各々の相関関数推定値ごとに最大10個までのフレームを用いる可能性のある上記のChen等によって開示された完全相関関数よりも計算効率が遥かによいので、有利である。本発明の相関インデクスの範囲は、同等なカーネルについての100%から、完全に独立なスペックル・パターンについてのゼロまでにわたる。原理的には、この相関インデクスはまた、2つのカーネルが異なる構造をピックアップする場合には負となる可能性もある。
【0021】
一般に、1つのフレームに基づいて選択されたカーネルがまた、他のフレームにおいて均一なスペックル・パターンを含んでいるという保証はない。従って、相関インデクス推定値のスクリーニング試験(図4の工程46)を行うことが望ましい。例えば、軸方向及び/若しくは横方向の走査平面のぶれ、仰角方向の動きの滑り(skid)又はフレーム平均(これは、相関していないフレームを弱い相関があるように見せる可能性がある)を示すものである一定の信頼性閾値(例えば、20%)を下回るすべての相関インデクスのサンプルを除去するように選択することができる。また、残りの良好な相関インデクス値を数えて(図4の工程47)、これらの相関インデクス値がN個のフレームのうちの有意な部分を構成しているかどうかを知ることも有用である(例えば、「少なくとも10%の相関インデクス値が良好でなければならない」等)。信頼性のあるフレーム(CI>20%)があまりにも少ないようであれば、より遅い又はより一定した速度で再走査を行うように、或いは全掃引距離の推定値を手入力するように、ユーザに促さなければならない(工程48)。
【0022】
良好な相関インデクスのサンプルが十分に存在しているならば、これらの平均を求めて(図4の工程50)、統計的なばらつきを減少させるべきである。この結果を用いて、各々のプローブの種類及びカーネルの深さごとに予め較正されているCI対dモデル(図2に示すCPU28のメモリに記憶されている)に基づいて、対応する平均スライス間間隔dを算出することができる(工程52)。良好な相関インデクスのサンプルが十分に存在しているならば、対応する平均dは、3次元再構成のために信頼性が非常に高いものとなるはずである。
【0023】
本発明の核心にあるのは、結果として得られる相関インデクスが妥当な限度内で表示ダイナミック・レンジ、後処理フィルタ及びカーネル寸法に独立であるように、2つの画像カーネルのSADを正規化する適応方法である。主眼の着想は、理論的なスペックル統計から、画像のカーネルが非常にかけ離れているため統計的に独立となった場合(フレーム平均化がない場合)にピクセル当たりの平均SADが接近していく値を決定することである。
【0024】
均一な拡散散乱媒体について検出されるスペックル振幅は、レイリー(Rayleigh)分布によって記述されることは周知である。表示の前の画像の圧縮が、以下の単純な対数関数によってモデル化され得るものとする。
y=10 log[x+1] (1)
標準的な統計的演算が示すところによると、xがレイリー分布しているならば、yの確率密度関数(pdf)は、
fy(y)=(a/m)exp(ay)exp(−[1/m][exp(ay)−1]) (2)
となり、ここで、a=(0.1)ln(10)は定数であり、m=E[x]はシステム利得に依存するxの期待値である。8ビットの線形グレイ・マップについては、yが表示のために[0,255]にマッピングされ、m=100についてのサンプルのヒストグラムが図5にプロットされている。
【0025】
y1 及びy2 が、相関する2つのカーネル内の対応するピクセルの振幅を表す同等に分布した独立確率変数であるとする。ここで必要なのは、SAD=|y1 −y2 |のpdfを決定することである。先ず、y2 のpdfがy1 のpdfと同じであると仮定すると、(−y2 )のpdfは単純にy2 のpdfの鏡像である。y1 +(−y2 )のpdfは、これらのそれぞれのpdfのコンボリューション(畳み込み)によって与えられ、これを図6に示す。この和のpdfは、ゼロに関して対称な関数であるので、その絶対値(SAD)のpdfは単純に、分布の正の半分(の2倍)である。
【0026】
(y1 −y2 )のpdfの半値幅(図6)が、個別のyのpdfの半値幅(図5)の約1.5倍であることに留意されたい。実用上は、圧縮関数は式(1)と異なっていてもよいが、常に同じアプローチに従ってSADのpdfを導出することができる。圧縮関数が式(1)から大幅に逸脱しておらず、また、後処理の空間フィルタリング効果が穏やかなものであれば、SADヒストグラムの幅を各々の画像カーネルのSADヒストグラム幅の定数γ倍によって近似することができ、ここではγ≒1.5である。空間フィルタが極めて強いものである場合は、γの値をこれに応じて調節しなければならないかも知れない。当業者には明らかなように、k番目の画像フレームにおける(M×M)ピクセルのカーネルにわたるピクセル振幅分布の幅を、そのrmsピクセル値によって定義することもできるし、又は本発明の好ましい実施例に従って、単純に以下のような平均絶対偏差によって定義することもできる。
【0027】
【数1】
【0028】
ここで、yk (i)はi番目のピクセルの振幅であり、また、meank はすべてのM2 個のピクセルについての平均振幅である。一般的には、sk は、ダイナミック・レンジ設定及び後処理用空間フィルタの関数である。
SADk 画像が、以下のようにしてk番目及び(k+1)番目のフレームの対応する2つのカーネルについて算出されたとする。
【0029】
【数2】
【0030】
フレーム離隔距離が増大するにつれて、2つのカーネルは統計的に独立になってゆくので、ピクセル当たりの平均絶対差、即ちSADk /M2 は、ゼロから限度値γsk まで増大する。従って、適切な相関インデクスを以下のように定義することができる。
【0031】
【数3】
【0032】
ここで、対数圧縮されたスペックル画像についてはγ≒1.5である。
異なる種類のイメージング用アレイを用いて行われた実験によると、この相関インデクスは、スライス間間隔dについての以下のような指数減衰関数によって極めて良好に記述され得ることがわかった。
CI= exp(−d/Dz ) (6)
ここで、Dz はプローブの仰角方向のビーム・プロファイルの特性値である相関外れ長さ(decorrelation length)である。仰角方向のビーム・プロファイルが回折及び組織減衰の影響から深さに応じて変わるので、Dz は一般的には所与の送信周波数について深さzの関数となる。ビーム・プロファイルは一般的には近距離音場では干渉性(コヒーレンス)がより低く且つ複雑さがより少ないので、Dz は中距離音場及び遠距離音場におけるよりも近距離音場で小さくなると予測される。
【0033】
ある相関インデクス推定値が与えられたとき、対応するスライス間間隔は、以下のようにして算出することができる。
d=−Dz ln(CI) (7)
CIがdによって指数的に減衰するという事実は、3次元走査にとっての利点となることがわかるであろう。即ち、CIは、小さな変位又は遅い動きに極めて敏感なものとなるからである。dが仰角方向のスライス厚みを超えて増大するにつれて、CIはゼロに向かって緩やかに減少する。式(7)から、CI<20%の場合には、CIの小さな変化がdの大きな変化に変換する可能性がある。従って、不良な相関インデクスのサンプルを除去するための妥当な信頼性閾値(図4の工程46)は、例えばCI=20%とし、即ち、この閾値を下回る値はすべて、手ぶれ又はプローブの滑りによって生じている可能性がある。式(7)を用いて、信頼性閾値よりも大きいすべての相関インデクス値の平均値に基づいて、N個のフレームについての平均スライス間間隔を算出することができる。
【0034】
走査中にプロセッサを用いて平均スライス間間隔dを算出するのではなく、ルックアップ・テーブル等のその他の手段によってdとCIとの間の関係を指定しておくこともできる。但し、式(6)の指数モデルは、各々の深さにおいて、この関係が相関外れ長さによって完全に特定されるという簡便さを提供する。所与のプローブについて、モータ駆動式プローブ・ホルダを用いてプローブを均一な散乱ファントムの全体にわたって一定の速度で平行移動させるような制御された実験を行うことにより、異なる深さ又はz間隔について相関外れ長さを較正することができる。
【0035】
CIの導出に当たっては、フレーム平均化は想定されていない。実用上は、フレーム平均化は、隣接する音波フレームがどれだけ遠くに離隔していようがこれらのフレーム相互の間に付加的な相関を生成する傾向を有する。このことは、フレーム平均化を行う結果、フレーム平均化の程度に応じて、相関インデクスがdの増大に伴ってより緩やかに減衰する(より大きな有効Dz )と共に、或る非ゼロの基準線のレベルに近付くということを意味している。このことは、実験研究によって確認されており、これによると、より大きな有効Dz が用いられると共に、相関インデクスについての信頼性閾値がフレーム平均化効果に関連する非ゼロの基準線の相関値よりも大きな値に選択される限り、式(6)の指数減衰モデルが依然として成立することがわかった。
【0036】
以上の好ましい実施例は、説明の目的のために開示されたものである。本発明の基本的な概念の変形及び改変は、当業者には容易に明らかであろう。具体的には、本発明を用いて、平均値ではない、隣接する走査平面相互の間の間隔を算出し得ることが認められよう。これらのような変形及び改変はすべて、特許請求の範囲によって包含されているものとする。
【図面の簡単な説明】
【図1】実時間超音波イメージング・システム内の主要な機能的サブシステムを示すブロック図である。
【図2】ピクセル強度データの相次ぐ3次元投影を含んでいるフレームを再構成するための手段を示すブロック図である。
【図3】1つの画像フレーム内での典型的な関心領域ボックスを示すと共に、関心領域内の選択されたカーネルを示している概略図である。
【図4】3次元再構成に用いるために平均フレーム間間隔dを推定する好ましい実施例による方法の工程を示す流れ図である。
【図5】m=E[x]がxの期待値であるとして、m=100について、対数圧縮されたノイズ・スペクトル・パワーの確率密度分布fy (y)を示すグラフである。
【図6】y1 及びy2 が、相関する2つのカーネルにおける対応するピクセルの振幅を表す同等に分布した独立確率変数であるとして、確率密度分布fy (SAD=|y1 −y2 |)を示すグラフである。
【符号の説明】
2 超音波トランスデューサ・アレイ
4 ビームフォーマ
6 Bモード処理チェーン
8 音波フレーム平均化ブロック
10 走査変換器
12 ビデオ・フレーム平均化ブロック
14 ビデオ・プロセッサ
16 シネ・メモリ
18 表示システム
20 ホスト・コンピュータ
22 オペレータ・インタフェイス
24 音響線メモリ
26 X−Yメモリ
28 中央処理装置(CPU)
30 ランダム・アクセス・メモリ
32 システム制御バス
Claims (16)
- 物体領域の3次元イメージング・システムにおいて、超音波ビームを送信すると共に前記物体領域から反射した超音波エコーを走査平面内の多数のサンプル領域において検出する超音波トランスデューサ・アレイと、該超音波トランスデューサ・アレイに結合されていて、前記物体領域を貫通する多数の平面の各々から反射した超音波エコーから導出されるイメージング・データを収集する手段と、前記多数の走査平面に対応する多数の画像フレームの各々についてのデータを記憶する手段と、相次ぐ画像フレームのそれぞれの対のそれぞれのカーネルのデータ相互の間の相関度を決定する手段であって、それぞれの相関度が前記それぞれのカーネルの対応するデータ相互の間の絶対差の正規化された和の関数である手段と、前記多数の走査平面のうちの相次ぐ走査平面相互の間の間隔を前記相関度の関数として算出する手段と、表示モニタと、前記表示モニタに前記間隔の関数である投影画像を表示させる手段とを有している3次元イメージング・システム。
- 前記間隔は平均値である請求項1に記載のシステム。
- 前記の相関度を決定する手段は、実質的に重複した連続した画像フレームをスクリーニングして除去する手段を含んでいる請求項1に記載のシステム。
- 前記の実質的に重複した相次ぐ画像フレームをスクリーニングして除去する手段は、前記絶対差の和を閾値と比較する手段を含んでいる請求項3に記載のシステム。
- 前記の相関度を決定する手段は、相次ぐ画像フレームの各々の対ごとにそれぞれの相関インデクスを推定する手段と、該相関インデクスの推定値を平均する手段とを含んでいる請求項1に記載のシステム。
- 前記間隔は、以下の方程式d=−Dz ln(CI)
に従って決定される値dを有しており、ここで、Dz は前記超音波トランスデューサ・アレイの仰角方向のビーム・プロファイルの特性値である相関外れ長さであり、またCIは前記相関インデクスの推定値の平均である請求項5に記載のシステム。 - 前記の間隔を算出する手段は、前記値dを記憶しているルックアップ・テーブルを含んでいる請求項6に記載のシステム。
- 前記の間隔を算出する手段は、プロセッサと、該プロセッサに前記相関外れ長さDz を入力する手段とを含んでいる請求項6に記載のシステム。
- 走査平面の動きを追跡するシステムであって、超音波ビームを送信すると共に多数の走査平面から反射した超音波エコーを検出する超音波スキャナと、前記超音波スキャナに結合されていて、前記走査平面のそれぞれから反射した超音波エコーから導出される多数の画像フレームの各々についてのデータを収集する手段と、相次ぐ画像フレームのそれぞれの対のそれぞれのカーネルのデータ相互の間の相関度を決定する手段であって、それぞれの相関度が前記それぞれのカーネルの対応するデータ相互の間の絶対差の正規化された和の関数である手段と、前記多数の走査平面のうちの相次ぐ走査平面相互の間の間隔を前記相関度の関数として算出する手段とを有しているシステム。
- 前記間隔は平均値である請求項9に記載のシステム。
- 前記の相関度を決定する手段は、実質的に重複した相次ぐ画像フレームをスクリーニングして除去する手段を含んでいる請求項9に記載のシステム。
- 前記の実質的に重複した相次ぐ画像フレームをスクリーニングして除去する手段は、前記絶対差の和を閾値と比較する手段を含んでいる請求項11に記載のシステム。
- 前記の相関度を決定する手段は、相次ぐ画像フレームの各々の対ごとにそれぞれの相関インデクスを推定する手段と、該相関インデクスの推定値を平均する手段とを含んでいる請求項9に記載のシステム。
- 前記間隔は、以下の方程式d=−Dz ln(CI)
に従って決定される値dを有しており、ここで、Dz は前記超音波スキャナの仰角方向のビーム・プロファイルの特性値である相関外れ長さであり、またCIは前記相関インデクスの推定値の平均である請求項13に記載のシステム。 - 前記の間隔を算出する手段は、前記値dを記憶しているルック・アップ・テーブルを含んでいる請求項14に記載のシステム。
- 前記の間隔を算出する手段は、プロセッサと、該プロセッサに前記相関外れ長さDz を入力する手段とを含んでいる請求項14に記載のシステム。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/045780 | 1998-03-20 | ||
US09/045,780 US6012458A (en) | 1998-03-20 | 1998-03-20 | Method and apparatus for tracking scan plane motion in free-hand three-dimensional ultrasound scanning using adaptive speckle correlation |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH11313822A JPH11313822A (ja) | 1999-11-16 |
JP4306864B2 true JP4306864B2 (ja) | 2009-08-05 |
Family
ID=21939843
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP07470999A Expired - Lifetime JP4306864B2 (ja) | 1998-03-20 | 1999-03-19 | 3次元イメージング・システム並びに走査平面の動きを追跡する方法及び装置 |
Country Status (4)
Country | Link |
---|---|
US (2) | US6012458A (ja) |
JP (1) | JP4306864B2 (ja) |
DE (1) | DE19910771A1 (ja) |
IL (1) | IL128866A (ja) |
Families Citing this family (72)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0839497A1 (en) * | 1996-11-01 | 1998-05-06 | EndoSonics Corporation | A method for measuring volumetric fluid flow and its velocity profile in a lumen or other body cavity |
US6269176B1 (en) * | 1998-12-21 | 2001-07-31 | Eastman Kodak Company | Method for x-ray antiscatter grid detection and suppression in digital radiography |
US6618493B1 (en) * | 1999-11-26 | 2003-09-09 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for visualization of motion in ultrasound flow imaging using packet data acquisition |
US6402693B1 (en) * | 2000-01-13 | 2002-06-11 | Siemens Medical Solutions Usa, Inc. | Ultrasonic transducer aligning system to replicate a previously obtained image |
US6728394B1 (en) * | 2000-02-14 | 2004-04-27 | Siemens Medical Solutions Usa, Inc. | Dynamic measurement of object parameters |
US6530885B1 (en) * | 2000-03-17 | 2003-03-11 | Atl Ultrasound, Inc. | Spatially compounded three dimensional ultrasonic images |
US6656119B2 (en) * | 2000-03-17 | 2003-12-02 | Kabushiki Kaisha Toshiba | Imaging diagnostic apparatus and maintenance method of the same |
KR100355718B1 (ko) * | 2000-06-10 | 2002-10-11 | 주식회사 메디슨 | 스티어링이 가능한 프로브를 사용한 3차원 초음파 영상시스템 및 영상 형성 방법 |
CA2314794A1 (en) * | 2000-08-01 | 2002-02-01 | Dimitre Hristov | Apparatus for lesion or organ localization |
JP2002102223A (ja) * | 2000-10-03 | 2002-04-09 | Mitani Sangyo Co Ltd | 超音波断層画像における面座標検出方法ならびにシステムおよび同方法がプログラムされ記録された記録媒体 |
JP4659974B2 (ja) * | 2000-12-12 | 2011-03-30 | 株式会社東芝 | 超音波診断装置 |
US6599248B1 (en) | 2001-03-20 | 2003-07-29 | Aloka | Method and apparatus for ultrasound diagnostic imaging |
US7846096B2 (en) | 2001-05-29 | 2010-12-07 | Ethicon Endo-Surgery, Inc. | Method for monitoring of medical treatment using pulse-echo ultrasound |
US20030069502A1 (en) | 2001-05-29 | 2003-04-10 | Makin Inder Raj. S. | Ultrasound feedback in medically-treated patients |
US6605042B2 (en) | 2001-08-10 | 2003-08-12 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for rotation registration of extended field of view ultrasound images |
US6517486B1 (en) * | 2001-08-16 | 2003-02-11 | Computed Ultrasound Global, Inc. | Compounding method for reducing speckle noise |
GB0206413D0 (en) * | 2002-03-19 | 2002-05-01 | Refrigerant Products Ltd | Refrigerant for centrifugal compressors |
KR100437974B1 (ko) * | 2002-05-11 | 2004-07-02 | 주식회사 메디슨 | 측면거리 상관함수를 이용한 3차원 초음파 영상 형성 방법및 장치 |
CN100393283C (zh) * | 2002-09-12 | 2008-06-11 | 株式会社日立医药 | 生物体组织动状态跟踪方法、使用该方法的图像诊断装置 |
US7356172B2 (en) * | 2002-09-26 | 2008-04-08 | Siemens Medical Solutions Usa, Inc. | Methods and systems for motion tracking |
JP3748848B2 (ja) * | 2002-11-11 | 2006-02-22 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 超音波診断装置 |
JP2005058551A (ja) * | 2003-08-18 | 2005-03-10 | Ge Medical Systems Global Technology Co Llc | 超音波画像表示方法および超音波診断装置 |
JP4729491B2 (ja) * | 2003-09-09 | 2011-07-20 | ゼネラル・エレクトリック・カンパニイ | 超音波ドップラー・カラーフロー・イメージングのための運動適応型フレーム平均化 |
US7658714B2 (en) * | 2003-10-31 | 2010-02-09 | Siemens Medical Solutions Usa, Inc. | Intelligent ultrasound examination storage system |
US7596255B2 (en) * | 2003-11-26 | 2009-09-29 | Ge Medical Systems Information Technologies, Inc. | Image navigation system and method |
US7806839B2 (en) | 2004-06-14 | 2010-10-05 | Ethicon Endo-Surgery, Inc. | System and method for ultrasound therapy using grating lobes |
US7128712B2 (en) * | 2004-06-21 | 2006-10-31 | General Electric Company | Adaptive ultrasound imaging system |
US7672705B2 (en) * | 2004-07-19 | 2010-03-02 | Resonant Medical, Inc. | Weighted surface-to-surface mapping |
US7729744B2 (en) * | 2004-07-20 | 2010-06-01 | Resonant Medical, Inc. | Verifying lesion characteristics using beam shapes |
US8000442B2 (en) * | 2004-07-20 | 2011-08-16 | Resonant Medical, Inc. | Calibrating imaging devices |
ES2365063T3 (es) * | 2004-09-20 | 2011-09-21 | Resonant Medical Inc. | Vigilancia de un tratamiento por radioterapia por medios de ultrasonidos. |
WO2006103603A1 (en) * | 2005-03-28 | 2006-10-05 | Koninklijke Philips Electronics, N.V. | Adaptive parallel artifact mitigation |
US20090306509A1 (en) * | 2005-03-30 | 2009-12-10 | Worcester Polytechnic Institute | Free-hand three-dimensional ultrasound diagnostic imaging with position and angle determination sensors |
US20080219405A1 (en) * | 2005-09-06 | 2008-09-11 | Tony Falco | System and method for patient setup for radiotherapy treatment |
US8929621B2 (en) * | 2005-12-20 | 2015-01-06 | Elekta, Ltd. | Methods and systems for segmentation and surface matching |
US8035637B2 (en) * | 2006-01-20 | 2011-10-11 | 3M Innovative Properties Company | Three-dimensional scan recovery |
AU2007272373B2 (en) * | 2006-07-13 | 2011-10-27 | The Regents Of The University Of Colorado | Echo particle image velocity (EPIV) and echo particle tracking velocimetry (EPTV) system and method |
US9451928B2 (en) | 2006-09-13 | 2016-09-27 | Elekta Ltd. | Incorporating internal anatomy in clinical radiotherapy setups |
US8699765B2 (en) * | 2006-11-13 | 2014-04-15 | Siemens Medical Solutions Usa, Inc. | Reducing jittering in medical diagnostic ultrasound imaging |
US7925068B2 (en) * | 2007-02-01 | 2011-04-12 | General Electric Company | Method and apparatus for forming a guide image for an ultrasound image scanner |
JP5624258B2 (ja) * | 2007-04-26 | 2014-11-12 | 株式会社東芝 | 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム |
WO2009012577A1 (en) | 2007-07-20 | 2009-01-29 | Resonant Medical Inc. | Methods and systems for compensating for changes in anatomy of radiotherapy patients |
US10531858B2 (en) | 2007-07-20 | 2020-01-14 | Elekta, LTD | Methods and systems for guiding the acquisition of ultrasound images |
JP5394620B2 (ja) * | 2007-07-23 | 2014-01-22 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 超音波撮像装置および画像処理装置 |
US8135198B2 (en) | 2007-08-08 | 2012-03-13 | Resonant Medical, Inc. | Systems and methods for constructing images |
JP5108444B2 (ja) * | 2007-10-11 | 2012-12-26 | 富士フイルム株式会社 | 超音波診断方法及び装置 |
US8189738B2 (en) | 2008-06-02 | 2012-05-29 | Elekta Ltd. | Methods and systems for guiding clinical radiotherapy setups |
US9826959B2 (en) * | 2008-11-04 | 2017-11-28 | Fujifilm Corporation | Ultrasonic diagnostic device |
JP5355074B2 (ja) * | 2008-12-26 | 2013-11-27 | キヤノン株式会社 | 3次元形状データ処理装置、3次元形状データ処理方法及びプログラム |
US10542962B2 (en) | 2009-07-10 | 2020-01-28 | Elekta, LTD | Adaptive radiotherapy treatment using ultrasound |
US9248316B2 (en) | 2010-01-12 | 2016-02-02 | Elekta Ltd. | Feature tracking using ultrasound |
US20110172526A1 (en) * | 2010-01-12 | 2011-07-14 | Martin Lachaine | Feature Tracking Using Ultrasound |
JP4709937B2 (ja) * | 2010-10-01 | 2011-06-29 | 株式会社東芝 | 超音波診断装置及び画像処理装置 |
US9208570B2 (en) * | 2012-03-28 | 2015-12-08 | Sony Corporation | System and method for performing depth estimation by utilizing an adaptive kernel |
US20130131510A1 (en) * | 2011-05-30 | 2013-05-23 | Tadamasa Toma | Ultrasound image generation apparatus and ultrasound image generation method |
JP5389118B2 (ja) * | 2011-08-19 | 2014-01-15 | 株式会社東芝 | 超音波診断装置 |
CA2851590A1 (en) * | 2011-10-10 | 2013-04-18 | Tractus Corporation | Method, apparatus and system for complete examination of tissue with hand-held imaging devices |
US9020294B2 (en) * | 2012-01-18 | 2015-04-28 | Dolby Laboratories Licensing Corporation | Spatiotemporal metrics for rate distortion optimization |
US8891881B2 (en) | 2012-01-25 | 2014-11-18 | General Electric Company | System and method for identifying an optimal image frame for ultrasound imaging |
EP2807978A1 (en) | 2013-05-28 | 2014-12-03 | Universität Bern | Method and system for 3D acquisition of ultrasound images |
WO2014210430A1 (en) | 2013-06-27 | 2014-12-31 | Tractus Corporation | Systems and methods for tissue mapping |
JP6835587B2 (ja) * | 2014-02-19 | 2021-02-24 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 医用4dイメージングにおける動き適応型可視化 |
EP3247281B1 (en) | 2015-01-23 | 2020-12-02 | The University of North Carolina at Chapel Hill | Apparatuses, systems, and methods for preclinical ultrasound imaging of subjects |
HU231249B1 (hu) | 2015-06-26 | 2022-05-28 | Dermus Kft. | Eljárás ultrahangkép előállítására és számítógépes adathordozó |
JP6997508B2 (ja) * | 2016-08-31 | 2022-01-17 | ゼネラル・エレクトリック・カンパニイ | 超音波診断装置及びその制御プログラム |
WO2019048482A1 (en) | 2017-09-07 | 2019-03-14 | Piur Imaging Gmbh | APPARATUS AND METHOD FOR DETERMINING THE MOTION OF AN ULTRASONIC PROBE |
US11504095B2 (en) | 2018-01-08 | 2022-11-22 | Rivanna Medical, Inc. | Three-dimensional imaging and modeling of ultrasound image data |
CN109064491A (zh) * | 2018-04-12 | 2018-12-21 | 江苏省基础地理信息中心 | 一种自适应分块的核相关滤波跟踪方法 |
SG10201810322YA (en) * | 2018-11-19 | 2020-06-29 | Medo Dx Pte Ltd | Method and system for generating a three-dimensional ultrasound image of a tissue volume |
EP3705049A1 (en) | 2019-03-06 | 2020-09-09 | Piur Imaging GmbH | Apparatus and method for determining motion of an ultrasound probe including a forward-backward directedness |
US11911213B2 (en) | 2019-06-03 | 2024-02-27 | General Electric Company | Techniques for determining ultrasound probe motion |
CN113812978B (zh) * | 2021-10-25 | 2023-08-18 | 深圳市德力凯医疗设备股份有限公司 | 数据采样方法、病灶部位检查方法及智能终端 |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5226113A (en) * | 1989-10-30 | 1993-07-06 | General Electric Company | Method and apparatus for volumetric projection rendering using reverse ray casting |
US5331964A (en) * | 1993-05-14 | 1994-07-26 | Duke University | Ultrasonic phased array imaging system with high speed adaptive processing using selected elements |
US5485842A (en) * | 1994-11-30 | 1996-01-23 | Advanced Technology Laboratories, Inc. | Ultrasonic diagnostic scan conversion for three dimensional display processing |
US5538004A (en) * | 1995-02-28 | 1996-07-23 | Hewlett-Packard Company | Method and apparatus for tissue-centered scan conversion in an ultrasound imaging system |
US5817022A (en) * | 1995-03-28 | 1998-10-06 | Sonometrics Corporation | System for displaying a 2-D ultrasound image within a 3-D viewing environment |
US5655535A (en) * | 1996-03-29 | 1997-08-12 | Siemens Medical Systems, Inc. | 3-Dimensional compound ultrasound field of view |
US5782766A (en) * | 1995-03-31 | 1998-07-21 | Siemens Medical Systems, Inc. | Method and apparatus for generating and displaying panoramic ultrasound images |
US5899861A (en) * | 1995-03-31 | 1999-05-04 | Siemens Medical Systems, Inc. | 3-dimensional volume by aggregating ultrasound fields of view |
AU6276696A (en) * | 1995-06-15 | 1997-01-15 | Regents Of The University Of Michigan, The | Method and apparatus for composition and display of three-dimensional image from two-dimensional ultrasound |
US5582173A (en) * | 1995-09-18 | 1996-12-10 | Siemens Medical Systems, Inc. | System and method for 3-D medical imaging using 2-D scan data |
US5787889A (en) * | 1996-12-18 | 1998-08-04 | University Of Washington | Ultrasound imaging with real time 3D image reconstruction and visualization |
US5876342A (en) * | 1997-06-30 | 1999-03-02 | Siemens Medical Systems, Inc. | System and method for 3-D ultrasound imaging and motion estimation |
-
1998
- 1998-03-20 US US09/045,780 patent/US6012458A/en not_active Expired - Lifetime
-
1999
- 1999-03-08 IL IL12886699A patent/IL128866A/en not_active IP Right Cessation
- 1999-03-11 DE DE19910771A patent/DE19910771A1/de not_active Withdrawn
- 1999-03-19 JP JP07470999A patent/JP4306864B2/ja not_active Expired - Lifetime
- 1999-04-23 US US09/299,032 patent/US6106470A/en not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
US6106470A (en) | 2000-08-22 |
IL128866A (en) | 2002-12-01 |
US6012458A (en) | 2000-01-11 |
JPH11313822A (ja) | 1999-11-16 |
DE19910771A1 (de) | 1999-09-23 |
IL128866A0 (en) | 2000-01-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4306864B2 (ja) | 3次元イメージング・システム並びに走査平面の動きを追跡する方法及び装置 | |
JP4722283B2 (ja) | 連続データ獲得を用いた超音波フロー・イメージングにおける運動の可視化のための方法および装置 | |
US7223241B2 (en) | Method and apparatus for elasticity imaging | |
EP1046929B1 (en) | Method and apparatus for three-dimensional ultrasound imaging using surface-enhanced volume rendering | |
US6733455B2 (en) | System and method for adaptive clutter filtering in ultrasound color flow imaging | |
EP0882426B1 (en) | Ultrasonic imaging apparatus and method for breast examination | |
JP3935433B2 (ja) | 空間的複合のための適応画像処理 | |
JP4204095B2 (ja) | 被検体ボリュームの三次元イメージング・システムおよび方法 | |
US20070093716A1 (en) | Method and apparatus for elasticity imaging | |
US4817015A (en) | High speed texture discriminator for ultrasonic imaging | |
JP4627366B2 (ja) | パケット・データ獲得を用いた超音波フロー撮像における運動の可視化のための方法および装置 | |
US20130343627A1 (en) | Suppression of reverberations and/or clutter in ultrasonic imaging systems | |
JP2000152939A (ja) | 物体ボリュ―ムの3次元イメ―ジングのためのシステム及び方法 | |
JP7405950B2 (ja) | 微小脈管の高時空分解能超音波イメージングを行うための方法 | |
Skorton et al. | Range-and azimuth-dependent variability of image texture in two-dimensional echocardiograms. | |
JPH1128214A (ja) | 三次元イメージング・システムおよび方法 | |
US20070230756A1 (en) | Method and apparatus for processing ultrasound image signals | |
WO2007110669A1 (en) | Image data processing systems | |
US6500125B1 (en) | Ultrasound b/color priority threshold calculation | |
JP2000149015A (ja) | 画像のエッジ強調方法及びイメ―ジング装置 | |
EP0493909A2 (en) | Image processing system for producing an attenuation map of a scanned image | |
WO2022051609A1 (en) | Time-aligned plane wave compounding of ultrasound data | |
Insana et al. | Analysis of ultrasound image texture via generalized Rician statistics | |
Watkin et al. | Three-dimensional reconstruction and enhancement of arbitrarily oriented and positioned 2D medical ultrasonic images | |
Soler et al. | Superresolution spatial compounding techniques with application to 3D breast ultrasound imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060315 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20090204 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090224 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090227 |
|
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: 20090331 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20090402 |
|
RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20090402 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20090428 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120515 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120515 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130515 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140515 Year of fee payment: 5 |
|
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 |
|
EXPY | Cancellation because of completion of term |