JP2009172381A - 造影エンハンス超音波診断イメージング方法、コンピュータ読み取り可能な記憶媒体および造影エンハンス超音波診断イメージング装置 - Google Patents
造影エンハンス超音波診断イメージング方法、コンピュータ読み取り可能な記憶媒体および造影エンハンス超音波診断イメージング装置 Download PDFInfo
- Publication number
- JP2009172381A JP2009172381A JP2009012590A JP2009012590A JP2009172381A JP 2009172381 A JP2009172381 A JP 2009172381A JP 2009012590 A JP2009012590 A JP 2009012590A JP 2009012590 A JP2009012590 A JP 2009012590A JP 2009172381 A JP2009172381 A JP 2009172381A
- Authority
- JP
- Japan
- Prior art keywords
- data
- contrast agent
- contrast
- information
- frames
- 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
Links
- 239000002872 contrast media Substances 0.000 title claims abstract description 103
- 238000000034 method Methods 0.000 title claims description 39
- 238000002604 ultrasonography Methods 0.000 title claims description 34
- 238000002059 diagnostic imaging Methods 0.000 title claims description 27
- 230000000747 cardiac effect Effects 0.000 claims abstract description 57
- 230000033001 locomotion Effects 0.000 claims description 73
- 238000002607 contrast-enhanced ultrasound Methods 0.000 claims description 19
- 230000005540 biological transmission Effects 0.000 claims description 16
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims 1
- 230000017531 blood circulation Effects 0.000 abstract description 6
- 230000004087 circulation Effects 0.000 abstract description 6
- 210000000056 organ Anatomy 0.000 abstract description 4
- 238000003384 imaging method Methods 0.000 description 21
- 230000000875 corresponding effect Effects 0.000 description 17
- 230000015654 memory Effects 0.000 description 16
- 230000004044 response Effects 0.000 description 15
- 230000006870 function Effects 0.000 description 7
- 239000000872 buffer Substances 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 6
- 238000002592 echocardiography Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 6
- 230000002107 myocardial effect Effects 0.000 description 6
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 4
- 238000012935 Averaging Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000002861 ventricular Effects 0.000 description 3
- 239000002131 composite material Substances 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 210000001174 endocardium Anatomy 0.000 description 2
- 238000005562 fading Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 210000005240 left ventricle Anatomy 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000001960 triggered effect Effects 0.000 description 2
- 241001116459 Sequoia Species 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 230000003139 buffering effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 229940039231 contrast media Drugs 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 239000002961 echo contrast media Substances 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000008155 medical solution Substances 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/02—Measuring pulse or heart rate
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
- A61B8/5276—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts due to motion
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/481—Diagnostic techniques involving the use of contrast agent, e.g. microbubbles introduced into the bloodstream
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61K—PREPARATIONS FOR MEDICAL, DENTAL OR TOILETRY PURPOSES
- A61K49/00—Preparations for testing in vivo
- A61K49/22—Echographic preparations; Ultrasound imaging preparations ; Optoacoustic imaging preparations
-
- 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/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/06—Systems determining the position data of a target
- G01S15/08—Systems for measuring distance only
- G01S15/10—Systems for measuring distance only using transmission of interrupted, pulse-modulated waves
- G01S15/102—Systems for measuring distance only using transmission of interrupted, pulse-modulated waves using transmission of pulses having some particular characteristics
- G01S15/108—Systems for measuring distance only using transmission of interrupted, pulse-modulated waves using transmission of pulses having some particular characteristics using more than one pulse per sonar period
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8959—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using coded signals for correlation purposes
- G01S15/8963—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using coded signals for correlation purposes using pulse inversion
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52085—Details related to the ultrasound signal acquisition, e.g. scan sequences
- G01S7/52087—Details related to the ultrasound signal acquisition, e.g. scan sequences using synchronization techniques
- G01S7/52088—Details related to the ultrasound signal acquisition, e.g. scan sequences using synchronization techniques involving retrospective scan line rearrangements
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0883—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the heart
-
- 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/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Biophysics (AREA)
- General Physics & Mathematics (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Acoustics & Sound (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Computer Networks & Wireless Communication (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physiology (AREA)
- Hematology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Cardiology (AREA)
- Epidemiology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
【課題】画像の暈けなしに器官または組織の循環流ないし血流をイメージングできるようにする。
【解決手段】少なくとも部分的に造影剤情報を含む複数の心周期にわたる超音波データフレームのシーケンスが形成され、該シーケンスのうち心周期のほぼ同じ相に相応する複数の超音波データフレームから成る第１のサブセットが選択され、前記第１のサブセットの複数の超音波データフレームについて動き補償が行われ、動き補償された前記第１のサブセットの前記超音波データフレームに基づいて、少なくとも部分的に当該の心周期の相での造影剤情報を表す画像が形成される。
【選択図】図２
【解決手段】少なくとも部分的に造影剤情報を含む複数の心周期にわたる超音波データフレームのシーケンスが形成され、該シーケンスのうち心周期のほぼ同じ相に相応する複数の超音波データフレームから成る第１のサブセットが選択され、前記第１のサブセットの複数の超音波データフレームについて動き補償が行われ、動き補償された前記第１のサブセットの前記超音波データフレームに基づいて、少なくとも部分的に当該の心周期の相での造影剤情報を表す画像が形成される。
【選択図】図２
Description
本発明は医療分野の造影エンハンス超音波診断イメージング装置に関する。特に造影剤イメージング情報が所定の時間にわたって結合される。
造影心エコー法は心臓のイメージングに広く利用されている技術である。この技術の基本的な適用分野の１つに左心室不透明化ＬＶＯが挙げられる。ここでは造影剤のミクロ泡によって左心室を不透明化し、心内膜の境界を良好に視覚化することができる。造影エコーを発展させた心筋造影エコー法ＭＣＥによれば、造影剤のミクロ泡によって循環流および心筋血流のトレースが行われる。器官または組織の血流をイメージングできると有意義である。取得されたデータフレームを所定の時間にわたって積算すると有利な適用分野も存在する。取得された画像から毛細血管や血流路などの診断に有用な情報が提供される。
種々の手法を組み合わせた例として、最大強度処理法ＭＩＰ、最小強度保持法、および、時間強度曲線構成法ＴＩＣなどが挙げられる。最大強度処理法によれば、所定の時間にわたって、高い輝度コントラスト部分が結合される。時間強度曲線構成法ＴＩＣによれば、ピクセルごとまたは関心領域ＲＯＩごとの強度（例えば造影剤フレーム強度）が時間の関数として特性マップに表される。当該の特性マップにより、造影剤のインフロー、アウトフローまたはその双方が所定の時間にわたってデータフレームの成分に関連して表される。ただし、オペレータの運動または内的運動のために、種々のフレームの情報を結合すると画像の暈けや不正確な情報が生じる。フレーム間での動き補償のみではこうした暈けを完全に補正することはできない。なぜならイメージングされた組織が他の組織に対して運動しているからである。
本発明の基礎とする課題は、画像の暈けなしに器官または組織の循環流ないし血流をイメージングできるようにすることである。
本発明の第１の特徴は造影エンハンス超音波診断イメージング方法に関しており、ここでは、少なくとも部分的に造影剤情報を含む複数の心周期にわたる超音波データフレームのシーケンスが形成され、該シーケンスのうち心周期のほぼ同じ相に相応する複数の超音波データフレームから成る第１のサブセットが選択され、前記第１のサブセットの複数の超音波データフレームについて動き補償が行われ、動き補償された前記第１のサブセットの前記超音波データフレームに基づいて、少なくとも部分的に当該の心周期の相での造影剤情報を表す画像が形成される。
本発明の第２の特徴は造影エンハンス超音波診断イメージングのためにプログラミングされたプロセッサによって実行される命令を表すデータを記憶したコンピュータ読み取り可能な記憶媒体に関しており、ここでは、複数の心周期にわたる心周期の第１の相に相応する造影剤データフレームおよびＢモードデータフレームを取得する命令、前記造影剤データフレーム間の動きをＢモードデータフレームに基づいて動き補償する命令、ならびに、第１の相に相応する造影剤データを結合して造影剤データに基づく画像を形成する命令を表すデータが記憶されている。
本発明の第３の特徴は造影エンハンス超音波診断イメージング装置に関しており、ここでは、心周期情報を受信するＥＣＧ入力端子、造影剤検出器、Ｂモード検出器、および、前記造影剤検出器から出力されたデータフレームを前記Ｂモード検出器から出力されたデータフレームに基づいて空間的にアライメントするプロセッサが設けられており、前記造影剤検出器から出力されたデータフレームは１心周期内のほぼ同じ時点を表す心周期情報により空間的にアライメントされ、前記プロセッサにより、１心周期内のほぼ同じ時点を表す心周期情報を有する空間的にアライメントされたデータフレームが結合される。
最初に説明する以下の有利な実施例は、造影エンハンス超音波診断イメージング方法、造影エンハンス超音波診断イメージング装置、および、コンピュータで読み取りな可能な記憶媒体を含んでいる。共通の相の期間で得られた複数のデータフレームがグループ化される。動き補償は共通の相のグループの内部で行われる。造影剤を表す画像は共通の相のグループ内部の動き補償された各フレームを結合することによって形成される。
本発明は特許請求の範囲によって規定される。ここでの説明は特許請求の範囲を限定するものでないことに注意されたい。以下では、有利な実施例に則して、本発明のさらなる特徴や利点について説明する。
図の各構成要素は必ずしも縮尺通りに描かれてはおらず、本発明の原理を示すために拡大されていることがあることに注意されたい。なお、全図を通して、対応する部分には同じ参照番号を付してある。
プロセッサは種々の心周期にわたって同じ相での造影剤データフレームを抽出する。所定の相での複数の造影剤フレームがアライメントされる。造影剤データはＢモードデータを用いてアライメントすることができる。種々の心周期にわたる同じ相でアライメントされた造影剤データフレームが結合される。動き平均、最小強度保持法ＭＩＰまたはピーク時点などのパラメータ計算が動き補償された造影剤データフレームに適用される。
左心室不透明化ＬＶＯおよび心筋造影エコー法ＭＣＥは造影エコー法の２つの基本的な適用分野である。これらの適用分野は心周期の共通相を有する複数のフレームを動き補償してイメージングすることにより実現される。この方法によればスペックルノイズが低減され、信号雑音比が上昇する。コントラスト分解能および毛細血管の視認性も向上する。
図１には造影エンハンス超音波診断イメージング装置１０が示されている。イメージング装置１０は、送信ビーム形成装置１２，トランスデューサ１４，受信ビーム形成装置１６，イメージプロセッサ１８，相選択プロセッサ２０，ディスプレイ２２，ＥＣＧ検出器２４，ＥＣＧ入力端子２６およびメモリ２８を含む。コンポーネントを加えたり、別のものと置換したり、減らしたりすることができる。例えば、データフレームのバッファリングのための別のメモリを設けることができる。別の例として、相選択プロセッサ２０がイメージプロセッサ１８に結合されるかまたはその一部として構成されていてもよい。さらに別の例として、ＥＣＧ検出器２４が設けられず、ＥＣＧ入力端子２６がイメージプロセッサ１８または他のプロセッサの出力端子であってもよい。
イメージング装置１０は医療用の超音波診断イメージング装置であるが、他のモダリティを用いたイメージング装置であってもよい。イメージング装置１０はリアルタイム動作またはオフライン動作を行う。他の実施例として、装置１０の一部または全てがコンピュータまたはワークステーションとして実現されてもよい。例えば、以前に取得されたデータフレームはビーム形成装置１２，１６またはトランスデューサ１４なしで処理される。コンピュータまたはワークステーションのオフラインソフトウェアによって造影イメージングに対する共通相の動き補償を行うこともできる。
ＥＣＧ検出器２４は１つのプロセッサおよび１つまたは複数の電極リードを含む。患者の心周期が検出される。心周期、トリガイベントその他の心周期の特性がＥＣＧ入力端子２６へ出力される。受信された心周期情報は当該の心周期内の特定の時点に関連するデータフレームを識別するために用いられる。他の実施例では、ＥＣＧ入力端子２６は超音波データを分析するプロセッサ１８，２０の一部として構成される。
送信ビーム形成装置１２は超音波送信器、メモリ、パルス発生器、アナログ回路、ディジタル回路またはこれらの組み合わせである。送信ビーム形成装置１２は種々の振幅、遅延量および／またはフェージングを有する複数のチャネルに対する波形を形成する。形成された波に応じてトランスデューサ１４から音波が送信されると、１つまたは複数のビームが形成される。送信ビーム形成装置１２は特定の相および／または振幅を有するビームを形成する。例えば、送信ビーム形成装置１２は所定の走査線または隣接する走査線に関するパルスシーケンスを送信する。各パルスの振幅および／または相は各ビームのそれに相応する。他の実施例として、単一のビームを所定の走査線に対しておよび／または同じ振幅ないし相対的な相を有するビームに対して用いてもよい。
トランスデューサ１４は圧電素子または容量性メンブレイン素子の１次元、１．２５次元、１．５次元、１．７５次元または２次元のアレイである。トランスデューサ１４は音響エネルギと電気エネルギとのあいだの変換のための複数の変換素子を有している。これらの変換素子は送信ビーム形成装置１２および受信ビーム形成装置１６のチャネルに接続される。
受信ビーム形成装置１６は増幅器、遅延素子、位相回転器、１つまたは複数の累算器などを有する複数のチャネルを含む。各チャネルはトランスデューサの１つまたは複数の変換素子に接続されている。受信ビーム形成装置１６は、相対遅延量、位相および／またはアポダイゼーション量に基づき、各送信に応じて１つまたは複数の受信ビームを形成する。他の実施例として、受信ビーム形成装置１６はフーリエ変換または他の変換を用いてサンプルを形成するプロセッサであってもよい。
受信ビーム形成装置１６は２次高調波または送信周波数帯域とは異なる周波数帯域での情報を分離するためのフィルタを含むことができる。こうした情報は所望の組織、造影剤および／またはフロー情報を含む。他の実施例として、受信ビーム形成装置１６がメモリ、バッファおよびフィルタあるいは加算器を含んでいてもよい。２つまたはそれ以上の受信ビームが結合され、これにより、２次高調波、立方基本波または他の帯域などの所望の周波数帯域での情報が分離される。
送信動作および受信動作の所望のシーケンスは超音波情報の取得に用いられる。例えば、Ｂモードのデータは１つの領域を一度走査することにより取得される。Ｂモードデータは組織のイメージングに用いられる。相関または運動の追跡はＢモードデータから流体情報を導出するために用いられる。Ｂモード動作から造影剤情報が形成される。ドップラー情報は各走査線に沿ってビームシーケンスを送信することにより取得される。コーナーターニングメモリは組織、造影剤および／またはフロー情報をドップラー信号から分離するために用いられる。公知の他のモードまたは将来開発される他のモードを利用することもできる。
１つの実施例では、造影剤イメージングモードが行われる。造影剤イメージングは典型的にはＢモード技術またはドップラー技術によって行われる。２次調波、偶数調波、奇数調波、下調波その他での情報を分離することにより、造影剤情報が識別される。例えば２つのパルス技術が用いられる。各パルスは同じ振幅で異なる相を有する。応答を加算することにより、偶数調波に関連する情報が識別される。フィルタリングは交互に利用される。これに代えてまたはこれに加えて、相対的なフェージングが受信処理において形成される。
１つの実施例では、送信シーケンスは立方基本波に応じたエコー信号が形成されるように制御される。送信ビーム形成装置１２は少なくとも２つのそれぞれ異なる振幅レベルおよび少なくとも２つの逆相または異なる相を有する複数のパルスを送信する。送信出力は適切な方式にしたがって変更可能である。この変更は、例えば、トランスデューサの個々の変換素子に印加される電圧を調整することにより、または、特定のパルスを形成するトランスデューサの変換素子数または送信アパーチャを調整することにより行われる。
立方基本波での超音波データを取得するために、受信ビーム形成装置１６はラインメモリおよび加算器またはフィルタを含んでおり、当該の送信に相応する複数の信号が結合される。ラインメモリまたはバッファは物理的に分離されたメモリとして構成される。また、これらを共通の物理的デバイスの選択位置として構成してもよい。ビーム形成信号はラインメモリまたはバッファに記憶され、重みづけ加算器において重みづけ加算される。振幅および位相の双方に対する重みが重みづけ加算器において用いられる。メモリおよび加算器はアナログ技術またはディジタル技術のいずれによって実現されていてもよい。重みづけ加算器は個別の受信信号を重みづけすることによりコンポジット信号を形成して出力する。所定の空間的位置に対するコンポジット信号は立方基本波の応答に関連するサンプルである。
立方基本波情報の取得手段については米国特許第６４９４８４１号明細書に記載されており、その開示内容は本発明に関連するのでこれを参考文献とする。ここでは送信および受信を組み合わせたシーケンスが立方基本波情報の取得に用いられることが記載されている。また、そのほかの送信および受信を組み合わせたシーケンスについて米国特許第６６０２１９５号明細書、米国特許第６６３２１７７号明細書、米国特許第６６３８２２８号明細書および米国特許第６６８２４８２号明細書に記載されており、それらの開示内容も本発明に関連するのでこれを参考文献とする。一般に、種々の振幅および種々の位相を有する複数のパルスから成るパルスシーケンスが送信される。変更された振幅または同相での異なる振幅を用いて立方基本波情報を取得することもできる。受信信号をシーケンスに応じて結合することにより、立方基本波情報を含むサンプルが得られる。立方基本波情報は超音波造影剤に対して高度に特異的である。なぜなら造影剤は立方応答を形成するがトランスデューサおよび組織はきわめて小さな立方応答しか形成しないからである。立方基本波情報から組織のクラッタが除外され、造影剤に対する特異的なイメージングが可能となる。例えば、立方基本波情報を用いると、組織内の毛細血管のイメージングないし識別が容易となる。
イメージプロセッサ１８はＢモード検出器、ドップラー検出器、パルス波ドップラー検出器、造影剤検出器、相関プロセッサ、フーリエ変換プロセッサ、ＡＳＩＣ、汎用プロセッサ、制御プロセッサ、ＦＰＧＡ、ディジタルシグナルプロセッサ、アナログ回路、ディジタル回路またはこれらの組み合わせ、あるいは、これまで既知のあるいは将来開発される超音波サンプルから表示情報を検出する手段である。
造影剤検出器の１つの実施例として、イメージプロセッサ１８により、同じ領域または心門の位置を表す複数のサンプルから高速フーリエ変換を行うことが挙げられる。各サンプルは立方基本波に対応しており、パルス波ドップラーの表示は立方基本波情報から形成される。これについては前掲した各文献の造影剤検出器のいずれも利用可能である。造影剤検出器として他のコンポーネントを用いることもできる。例えばＢモード検出が行われる。他の例として、フィルタにより種々の送信情報を結合し、造影剤の応答（例えば２次高調波または立方基本波）を良好に分離することもできる。フィルタにより、送信された超音波信号の立方基本波または他の周波数帯域での１次情報が得られる。この場合、複数の信号の検出が行われる。
イメージプロセッサ１８は並列に接続されたＢモード検出器を含む。Ｂモード検出器はビームの同じサンプルまたは異なるサンプルに対して動作し、組織、造影剤またはこれらの応答を検出する。例えば、受信ビームのシーケンスから空間的位置ごとに１つずつの受信ビームが立方基本波情報の分離のために１次組織情報をイメージングするＢモード検出器へ供給される。
イメージプロセッサ１８は超音波データフレームを出力する。データフレームは取得フォーマット（例えば極座標系）、表示フォーマットまたはその他のフォーマットでフォーマット化されている。例えばデカルト座標系フォーマットまたは画像フォーマットへの走査変換が行われる。各データフレームは１次元、２次元または３次元の走査領域、例えば患者の運動またはトランスデューサの運動を考慮に入れたイメージング領域全体を表している。データフレームは単一のタイプのデータまたは複数のタイプのデータを含む。例えば１つのデータフレームは造影剤情報のみを含む。他の例として、１つのデータフレームが或る空間的位置での造影剤情報と他の空間的位置での他のタイプの情報、例えばＢモード情報またはドップラー情報とを含んでもよい。同じ空間的位置の同じフレームに異なるタイプのデータが含まれてもよい。さらに別の例として、異なるタイプのデータを異なるデータフレームに形成することもできる。
選択的な実施例として、イメージプロセッサ１８にはデータがネットワークまたはメモリからロードされる。例えば超音波データを取得するために、ＤＩＣＯＭまたは他の画像がロードされる。各画像が１つのデータフレームとなる。１つのフレームが異なるタイプのデータを含み、それぞれが重畳されていてもよい。これに代えて、異なるタイプの種々のデータフレームが唯一のタイプのデータのみを含んでいてもよい。また、他の例として、各フレームが分割されており、そのうちの或る一部が１つのタイプのデータを含み、別の一部が別のタイプのデータを含むこともできる。
相選択プロセッサ２０は、ＡＳＩＣ、相関プロセッサ、フーリエ変換プロセッサ、汎用プロセッサ、制御プロセッサ、ＦＰＧＡ、ディジタルシグナルプロセッサ、アナログ回路、ディジタル回路またはこれらの組み合わせ、あるいは、これまで既知のまたは将来開発されるデータフレームの類似性判別手段および／または変位判別手段である。相選択プロセッサ２０はデータフレームを受信し、フレームがＭＩＰ画像、ＴＩＣ画像または各フレームの情報を結合して形成された他の画像のいずれに含まれるかを判別する。相選択プロセッサ２０は各データフレームに関連する相に基づいて選択を行う。各データフレームが取得されると、心周期に対する位相が当該のデータフレームに対して求められる。造影剤検出器から出力されたデータフレームは１心周期内のほぼ同じ相を表す心周期情報をともなっており、グループごとに選択される。例えばそれぞれ異なる位相に対応する複数のグループが形成される。
相選択プロセッサ２０は造影剤検出器から出力されたデータフレームを空間的にアライメントさせる。各グループのデータフレームまたは位相ごとのデータフレームが相互に位置合わせされる。各データフレームはトランスデューサ位置センサなどの外部のセンサの情報に基づいて空間的にアライメントされる。各データフレームを造影剤データに基づいて空間的にアライメントさせることもできる。別の実施例では、造影剤データフレームはＢモード検出器または他の非造影剤検出器から出力されたデータフレームに基づいて空間的にアライメントされる。Ｂモードフレームは例えば走査線インタリーブ方式またはフレームインタリーブ方式により造影剤フレームと同じ時点または類似の時点で取得される。Ｂモードフレームの空間的アライメントにより、同じ時点または類似の時点で取得された造影剤フレームの空間的アライメントも表される。
相選択プロセッサ２０は残光フィルタ、他のフィルタ、累算器、アルファブレンディングバッファ、他のバッファ、メモリ、プロセッサ、加算器または種々のデータフレームの情報から画像を形成する手段を含む。相選択プロセッサ２０は造影剤検出器から出力され空間的にアライメントされた複数のデータフレームを結合する。心周期情報は心周期内のほぼ同じ時点を有する複数のフレームが結合されることを表している。例えば相選択プロセッサ２０はフレームごとの特定の空間的位置または連続して結合されたフレームの特定の空間的位置のデータを比較する。比較の結果、例えば最大値、平均値、最小値などに基づいて、いずれかの値が選択されるかまたは連続する結合フレームが所望の値を含むように更新される。他の実施例では、相選択プロセッサ２０により、平均値、全体値または時間の関数としての位置または領域を表す他の値が求められる。種々のタイプのイメージングが組み合わされて用いられる。
ディスプレイ２２はＣＲＴ、モニタ、ＬＣＤ、フラットパネルディスプレイ、プロジェクタまたはその他のディスプレイ装置である。ディスプレイ２２は画像を表示するための表示値を受け取る。表示値は１次元画像、２次元画像または３次元画像としてフォーマット化される。或る実施例では、表示値は種々の時点で取得されたデータフレームに基づいて形成されたＴＩＣ画像またはＭＩＰ画像などの画像を表している。造影剤検出器から出力されたデータフレームが結合されて画像が形成される。付加的なデータフレームが取得されて選択されると、画像が更新される。単一データフレームまたはデータフレームの一部から形成された他の画像を表示してもよい。
イメージプロセッサ１８および／または相選択プロセッサ２０は命令にしたがって動作する。メモリ２８はコンピュータで読み取り可能なメモリである。コンピュータで読み取り可能な記憶媒体は、造影エンハンス超音波診断イメージングのためにプログラミングされた１つまたは複数のプロセッサによって実行される複数の命令を表すデータを記憶している。処理、方法および／または技術のステップを実行するための命令は、ここでは、コンピュータで読み取り可能な記憶媒体またはメモリに記憶されている。これは例えばキャッシュ、バッファ、ＲＡＭ、リムーバブルメディア、ハードディスクあるいは他の種々のタイプの揮発性または不揮発性の記憶媒体である。図示の機能、ステップまたはタスクはコンピュータで読み取り可能な媒体に記憶された１つまたは複数の命令のセットに応じて実行される。機能、ステップまたはタスクは命令セット、記憶媒体、プロセッサまたは処理ストラテジなどのタイプに依存せずに実行可能であり、単独または任意に組み合わされたソフトウェア、ハードウェア、集積回路、ファームウェア、マイクロコードその他などによって実行される。同様に、処理ストラテジはマルチプロセシング、マルチタスキング、並列処理その他を含む。或る実施例では、命令はローカル装置またはリモート装置によって読み出されるリムーバブルメディアに記憶される。他の実施例として、命令がリモート位置に記憶されており、コンピュータネットワークまたは電話回線網を介して伝送されてもよい。さらに別の実施例では、命令は所定のコンピュータ、ＣＰＵ、ＧＰＵまたはその他のシステム内に記憶される。
図２には造影エンハンス超音波診断イメージング方法が示されている。本発明の方法は図１の装置１０または他の装置によって実行される。本発明の方法は図２に示されている順序にしたがって行うと有利であるが、もちろん別の順序によって行うこともできる。付加的な別のステップを加えることもできるし、ステップ３０ｂおよび／またはステップ３６など幾つかのステップを省略することもできる。
ステップ３０で、超音波データフレームのシーケンスが取得される。データフレームのシーケンスは、複数の超音波データフレームを取得することにより、または、以前に形成されたデータフレーム（例えばＤＩＣＯＭ画像）を取得することにより、形成される。データフレームはリアルタイムでのライブ走査によって取得したり記憶されたクリップから取得したりすることもできる。シーケンスは、例えば心周期ごとに１回または数回取得され、連続的または周期的である。
データフレームのシーケンスは種々の時点での走査領域を表す複数のデータフレームを含む。各データフレームは同一の領域またはオーバラップした領域を表している。また、患者に対してトランスデューサが動いて領域を外れてしまうアウトオブプレーンなどに起因して、異なる領域を表すフレームも存在する。
当該の領域は造影剤を含む領域または造影剤の投入後にこれを含む可能性のある領域である。造影剤は超音波エネルギに応答する。造影剤データフレームはステップ３０ａで取得される。ここでは複数のデータフレームまたは全てのデータフレームが造影剤情報を含む。造影剤情報は組織または流体からの応答を含む。１つの実施例では、造影剤情報は超音波信号の立方基本波で取得される。例えば、超音波信号は少なくとも２つの異なる振幅レベルおよび異なる位相を有する複数のパルスとして送信される。造影剤の分解を回避するかまたは最小限にするために、小さい振幅の送信、例えばＭＩ＜０．７の送信が用いられ、送信に対応する複数の信号が例えば加算または重みづけ加算により結合される。結合により造影剤に対応した１次データが得られる。データは各データフレームの関心領域ＲＯＩのそれぞれの空間的位置で取得される。
データフレーム、例えば造影剤のみを表すデータフレームまたは造影剤および組織の応答を表すデータフレームには唯一のタイプのデータが存在する。これに代えて、当該のデータフレームが種々のタイプのデータ、例えば同一のフレームのデータまたは異なるセットのデータを表していてもよい。例えば造影剤データフレームとＢモードデータフレームとは別個に取得される。Ｂモードデータフレームはステップ３０ｂで取得される。Ｂモード情報は造影剤情報とは別個に形成される。これに代えて、いずれかのパルス、例えば造影剤情報のための振幅最大値のパルスに応答したエコー信号をＢモード検出のために用いることもできる。Ｂモード情報または組織情報は他の情報を含むこともある。例えばパルスシーケンスの形成および／またはフィルタリングが行われ、基本波、２次高調波またはこれら双方での超音波信号から組織情報が得られる。
造影剤フレームおよびＢモードデータフレームは１心周期の各相に相応する。各フレームは所定の時点で取得される。取得されたフレームには心周期に対する取得のタイミングを表す時点、位相その他のタグが付される。各フレームは心周期の種々の部分に関連している。各フレームは心周期に基づくトリガにより所望の相で取得される。心エコーイメージング法では、ＥＣＧ信号またはＲ波信号が一般的に利用可能である。超音波イメージング法ではトリガシーケンスによるフレームの取得がＲ波に同期される。これに代えて、フレームの取得をトリガなしに行ってもよい。
シーケンスは複数の心周期にわたって続く。各フレームは複数のサイクルにわたって取得されるが、同一の位相に関連している。図３には各フレームが短い黒色のバーで表されており、これらのフレームが６つの心周期のそれぞれで取得されることがわかる。各心周期はＲ波で開始しているが、他の時点で開始していてもよい。データフレームは各サイクルにおけるいずれかの位相ごとに取得される。別の例では、タイミングが正確に一致していないケースも存在する。１つまたは複数の心周期が所定の位相に関連するデータフレームを有さないこともある。超音波データフレームのシーケンスは少なくとも部分的に造影剤情報を表している。
図２のステップ３２では、当該のシーケンスから１つまたは複数のサブセットが選択される。各サブセットには心周期のほぼ同じ位相に相応する複数のデータフレームが存在する。他の心周期や他のグループの隣接相よりも当該の相に近いフレームにおける変化は考慮される。トリガシーケンスにおいても非トリガシーケンスにおいても、種々の心周期にわたる類似した相に関連する複数のフレームはフレームタイムスタンプおよび／またはＲ波タイムスタンプを用いてグループ化される。図３の実施例では、点線の楕円形の囲みによって６つの心周期のシーケンスのフレーム２がグループ化されており、当該のシーケンスのサブセットが表されている。Ｒ波のグループまたは拡張期のグループなど、他の異なるサブセットを形成することもできる。
所定のサブセット内のフレームは他の特性に基づいて放棄されたり使用されなかったりすることもある。例えば、フレーム間の動きにあまり影響されていないデータフレームが選択され、フレーム間の動きに影響されたフレームは選択されない。望ましくない動きをともなうデータフレームは放棄される。このために所望の閾値または他の基準を用いることができる。
図２のステップ３４では各サブセット内のフレームが空間的にアライメントされ、動き補償が行われる。動き補償は心周期のほぼ同じ相に関連する複数のフレームにわたって行われる。所定の相に対して複数のフレームがアライメントされる。つまり各サブセットのフレームが相応のサブセット内で他のフレームに対して位置合わせされるのである。
ステップ３４ではデータフレームに対して動き補償が行われ、フレーム間での平面内の動きであるインプレーンモーションが補償される。インプレーンモーションはトランスデューサの運動、患者の運動および／または関心領域ＲＯＩ内の組織の運動によって生じる。動き補償のために、１次元または数次元での相対的な平行運動および／または回転運動が判別される。１つのフレームからのデータは他のフレームの種々の領域と相関され、最良の一致部分または充分な一致部分が識別される。ここでは相関、相互相関、差の絶対値の最小和および／または他の同一性の尺度が用いられる。大域的な動きも局所的な動きも補償することができる。例えば、複数の異なる領域にわたるフレーム間の動きが求められる。大域的な動きが複数の動き補償から求められるか、または、フレームを歪ませている領域ごとに個別に動き補償が適用される。剛性モデルまたは非剛性モデルが形成される。例えば平行運動および回転運動に加えて歪みなどの非剛性モデルが形成される。データフレーム全体またはデータの一部のウィンドウが最良の一致および相応の動きを判別するために用いられる。
新たなデータフレームが取得されるたびに、以前のフレーム、時間的に隣接するフレームまたは選択されたデータフレームが基準フレームとして用いられる。これに代えて、順次あるいは時間的に間隔を置いたデータフレームの比較に対して同じ基準フレームを用いてもよい。
この場合、フレーム間の変位を用いてフレーム間の空間的位置がアライメントされる。動き補償により、トランスデューサの運動、患者の運動または器官の運動に関連する動きが排除されるかまたは最小化される。これに代えて、フレーム間の動き補償を行わないケースもありうる。
動き補償のためにモーションセンサを用いることもできる。超音波データに基づく動き補償のために、超音波データフレームが用いられる。１つの実施例では、Ｂモードデータフレームが造影剤データフレームのアライメントの判別に用いられる。この場合、動き補償はＢモードデータに基づいて行われる。各フレームが同じ位相に関連している場合、相応のＢモードフレームまたは組織の応答データは、動きの剛性モデルまたは非剛性モデルを用いて動きパラメータを予測することによりアライメントされる。図４には、フレーム間動き予測におけるｎ番目のフレームのポイントＰｉ ｎ、ｎ−１番目のフレームのポイントＰｉ ｎ−１およびｎ−２番目のフレームのポイントＰｉ ｎ−２が示されている。Ｂモードフレーム間または組織フレーム間の類似性によりＢモードフレームまたは組織フレームのアライメントの度合が表される。同じ心周期の同じ位相でアライメントされたフレームが造影剤フレームとして取得されるとき、Ｂモードフレームのアライメントまたは組織フレームのアライメントは造影剤フレームの空間的アライメントを表している。Ｂモードでの動き補償により、造影剤の動きに起因するアライメント時の誤差が回避される。造影剤も運動しているので、造影剤の応答に基づく動き補償は不正確になることがあるからである。
図２のステップ３６で、画像が形成される。画像は１つのサブセットまたは共通の相の動き補償された超音波データフレームに基づいて形成される。画像は少なくとも部分的に心周期の所定の位相での造影剤情報を表している。動き補償により画像の乱れや暈けが制限されるかまたは回避される。画像は同じ位相に関連する複数のフレームが取得されるときに形成されるかまたは更新される。またこのときに他のサブセットの画像または他の位相の画像を形成することもできる。他の画像、例えばＢモード画像を形成することもできる。画像はグレースケール画像、カラー画像またはこれらの組み合わせであってよい。
ステップ３６ａでは、選択されたサブセットの超音波データフレームからの情報が結合される。選択されなかったサブセットの超音波データフレームからの情報は結合されない。新たなデータフレームまたは新たな画像が選択されたフレームのデータに基づいて形成される。選択された超音波データフレームは時間の関数として積算される。積算は数学的積分も含むし、また、複数のソースから画像を形成することも含む。造影剤情報、例えば超音波信号の立方基本波での１次情報を結合することにより、造影剤の循環流および／または毛細血管系の血流をより良く観察することができる。
関心領域ＲＯＩの空間的位置ごとに、または、所定のサブセットのフレームによって表された全ての空間的位置に対して、データが比較されるかまたは値が求められる。各画像のピクセルごとに、共通の相によって選択された残りのフレームのデータに基づいて値が選択される。結合ステップは公知のまたは将来開発されるフレーム間プロセス、例えば、最大強度保持法、最小強度保持法、平均値形成法、時間強度曲線形成法などによって行われる。
例えば、データの平均値、メディアン値または他の統計値が空間的位置ごとのフレームから時間の関数として求められる。１つの実施例では、造影剤データは第１の位相に相応する造影剤データのフレームを平均することにより結合される。選択されたサブセットの超音波データフレームが例えば重みづけ移動平均により平均される。ｎ番目のフレームのポイントＰｉ ｎの強度をＩ（Ｐｉ ｎ）とする。ポイントＰｉ ｎに対する重みづけ移動平均の長さｋ（ｋ個のポイントの重みづけ移動平均）は次のように計算される。
ここで、Ｗｎ，Ｗｎ−ｉ，Ｗｎ−ｋ＋ｉは重み係数である。移動窓は移動平均にかけられる最も近時に取得されたサブセットのフレームの数を定義しており、例えばこの場合には３個のデータフレームが移動平均にかけられる。３個のフレームの移動平均は、心周期の拡張期の終了時に相応するフレームのサブセットである。無限または有限のパルス応答の平均が用いられる。他の実施例として、共通の相を有する全てのサブセットのフレームを重みづけ平均してもよい。
時間平均により一般に画像シーケンスのノイズリダクションが達成される。造影エンハンス超音波診断イメージングでは、所定の位置でのスペックルパターンは造影剤の動きに起因して変化しうる。その結果、時間平均によりスペックルノイズが低減される。
結合ステップの他の実施例として、選択されたフレームのデータに関連する最大値、最小値その他のデータを比較に基づいて選択することができる。選択されたサブセットのフレームは持続的なフレームまたは単一のフレームへと結合される。例えば、最大強度の投影フレームが造影剤データのフレームから形成される。造影剤情報の最大値は画像の空間的位置ごとまたは関心領域ＲＯＩごとに選択される。ｎ番目のフレームのポイントＰｉ ｎの強度をＩ（Ｐｉ ｎ）とする。ポイントｉの時間軸に沿った最大強度の投影は次のように計算される。
最大強度の投影は心周期の同じ相または類似の相を表すサブセットの動き補償されたフレームまたは移動窓にわたって実行される。
心筋造影エコー法では、造影剤は心筋の循環流の時間軌跡を表している。造影剤は所定の時間にわたって移動するので、当該の時間の全体にわたって変化情報を積算し、心筋の循環流を良好に視覚化することが望ましい。
結合ステップの別の例として、強度曲線または他の造影剤応答を表す時間曲線を複数のフレームから求めることもできる。曲線は１つの領域または１つの空間的位置に対して形成される。各フレームはそれぞれ異なる時点に対応しているので、強度曲線は時間の関数となる。
図２のステップ３６ｂでは、結合に基づいて画像が形成される。他のフレームとの結合により形成されたデータフレームに対して、走査変換、色符号化、マッピングおよび／または表示フォーマットへの変換が行われる。例えば、造影剤データはカラーマッピングされる。当該の結合画像を他の超音波情報とさらに結合することもできる。例えば結合画像を組織画像の上に重ねることができる。赤・青・緑（ＲＧＢ）の値または他の表示値は造影剤データから形成される。得られた画像は空間的アライメントを行った後の所定の心周期での造影剤情報を表している。
ステップ３６からステップ３２への戻り線によって示されているように、本発明の方法は反復可能である。図示の実施例では、選択ステップ、動き補償ステップ、結合ステップならびに画像形成ステップが反復される。取得ステップは反復されてもされなくてもよい。例えば、図３に示されているデータフレームが図２のステップ３０で取得される。取得ステップ３０を反復し、選択ステップ３２を別の時点で反復することもできる。動き補償ステップ３４の後、画像は種々の時点に対応する相ごとまたはサブセットごとに形成される。他の実施例として、所定の位相に関連する画像がより多くのデータフレームが取得されたときに更新されるように構成してもよい。
或る実施例では、Ｂモードイメージングは左心室の不透明化をともなって行われる。米国Siemens Medical Solutions社のSequoia c512なる名称の超音波装置および4V1cなる名称のトランスデューサは造影パルスシーケンス技術によりＣＰＳ立方基本波を用いてイメージングを行う。１０回の拍動にわたり先端方向の４室のシネクリップが取得される。当該のＤＩＣＯＭ動画像はソフトウェアによりオフラインで処理される。ソフトウェアにより拍動が時間的に同期され、１０回の拍動のシーケンスに対して３回の拍動ごとに重みづけ移動平均が行われる。最大強度の投影が取得された画像に適用される。左心室の心内膜の血流面が良好に視覚化される。こうして信号雑音比が改善され、また、造影の分解能も向上する。
これまで本発明について様々な実施形態に基づき説明してきたが、本発明の範囲から逸脱することなく種々の変更ないし修正を行えることは自明である。したがって上記の詳細な説明は制限ではなく説明を意図していることに注意されたい。本発明の内容および範囲を規定するのは特許請求の範囲の特徴およびこれと等価の特徴であることは明らかである。
１０ イメージング装置、 １２ 送信ビーム形成装置、 １４ トランスデューサ、 １６ 受信ビーム形成装置、 １８ イメージプロセッサ、 ２０ 相選択プロセッサ、 ２２ ディスプレイ、 ２４ ＥＣＧ検出器、 ２６ ＥＣＧ入力端子、 ２８ メモリ
Claims (20)
- 少なくとも部分的に造影剤情報を表すデータから複数の心周期にわたる超音波データフレームのシーケンスを形成するステップ（３０）、
該シーケンスのうち心周期のほぼ同じ相に相応する複数の超音波データフレームから成る第１のサブセットを選択するステップ（３２）、
前記第１のサブセットの複数の超音波データフレームについて動き補償するステップ（３４）、および、
動き補償された前記第１のサブセットの複数の超音波データフレームに基づいて少なくとも部分的に当該の心周期の相での造影剤情報を表す画像を形成するステップ（３６）
を有する
ことを特徴とする造影エンハンス超音波診断イメージング方法。 - 前記シーケンスを形成するステップ（３０）において、超音波信号の立方基本波の情報としてのデータを取得する（３０ａ）、請求項１記載の造影エンハンス超音波診断イメージング方法。
- 前記データを取得する際に（３０ａ）、少なくとも２つの異なる振幅レベルおよび異なる相を有する複数のパルスの形態の超音波信号を送信し、当該の送信に対応する複数の信号を結合する（３６ａ）、請求項２記載の造影エンハンス超音波診断イメージング方法。
- さらに、複数のパルスのうちの所定のパルスに対応するエコー信号に基づいてＢモード情報を形成し（３０ｂ）、前記動き補償するステップ（３４）において、前記超音波データフレームに相応する前記Ｂモード情報間の類似性に基づいて補償を行う、請求項３記載の造影エンハンス超音波診断イメージング方法。
- 前記動き補償するステップ（３４）においてＢモードデータに基づく補償を行い、前記画像を形成するステップ（３６）において造影剤情報を結合する（３６ａ）、請求項１記載の造影エンハンス超音波診断イメージング方法。
- 前記動き補償するステップ（３４）において、前記第１のサブセットの前記複数の超音波データフレーム間の動きの変位を求め、該動きの変位に基づいて前記第１のサブセットの前記超音波データフレームを空間的にアライメントする、請求項１記載の造影エンハンス超音波診断イメージング方法。
- 前記画像を形成するステップ（３６）において、前記第１のサブセットの前記複数の超音波データフレームを平均する、請求項１記載の造影エンハンス超音波診断イメージング方法。
- 前記画像を形成するステップ（３６）において、前記第１のサブセットの前記超音波データフレームから当該の画像の空間的位置ごとに前記造影剤情報の最大値を選択する、請求項１記載の造影エンハンス超音波診断イメージング方法。
- 前記第１のサブセットを選択するステップ（３２）、前記動き補償するステップ（３４）および前記画像を形成するステップ（３６）を当該の心周期の異なる相に相応する超音波データフレームから成る第２のサブセットについて反復する、請求項１記載の造影エンハンス超音波診断イメージング方法。
- 造影エンハンス超音波診断イメージングのためにプログラミングされたプロセッサによって実行される命令を表すデータを記憶したコンピュータ読み取り可能な記憶媒体（２８）において、
当該の記憶媒体は、複数の心周期にわたって心周期の第１の相に相応する造影剤データフレームおよびＢモードデータフレームを取得する命令（３０）、前記造影剤データフレーム間の動きを前記Ｂモードデータフレームに基づいて動き補償する命令（３４）、ならびに、前記第１の相に相応する造影剤データを結合して造影剤データに基づく画像を形成する命令（３６）を表すデータを記憶している
ことを特徴とするコンピュータ読み取り可能な記憶媒体。 - 前記造影剤データフレームを取得する際に、超音波信号の立方基本波での１次情報が取得され、さらに、超音波信号の基本波、２次高調波またはこれら双方の情報を含むＢモードデータフレームが取得される、請求項１０記載のコンピュータ読み取り可能な記憶媒体。
- 前記立方基本波での情報を取得する際に、少なくとも２つの異なる振幅レベルおよび異なる相を有する複数のパルスの形態の超音波信号が送信され、当該の送信に対応する複数の信号が結合され、さらに、前記Ｂモードデータを取得する際に、複数のパルスのうちの所定のパルスに対応するエコー信号に基づいてＢモード情報が形成される、請求項１１記載のコンピュータ読み取り可能な記憶媒体。
- 前記動き補償を行う際に、前記Ｂモードデータフレーム間の動きの変位が求められ、該動きの変位に基づいて前記造影剤データフレームが空間的にアライメントされる、請求項１０記載のコンピュータ読み取り可能な記憶媒体。
- 前記フレームを結合する際に、前記第１の相に相応する造影剤データのフレームが平均される、請求項１０記載の記コンピュータ読み取り可能な憶媒体。
- 前記フレームを結合する際に、前記造影剤データフレームからのデータにより最大値すなわち最大強度の投影フレームが形成される、請求項１０記載のコンピュータ読み取り可能な記憶媒体。
- 前記フレームの取得、前記動き補償および前記画像の形成が当該の心周期のうち前記第１の相とは異なる第２の相について反復される、請求項１０記載のコンピュータ読み取り可能な記憶媒体。
- 心周期情報を受信するＥＣＧ入力端子（２６）、造影剤検出器（１８）、Ｂモード検出器（１８）、および、前記造影剤検出器から出力されたデータフレームを前記Ｂモード検出器から出力されたデータフレームに基づいて空間的にアライメントするプロセッサ（２０）が設けられており、
前記造影剤検出器から出力されたデータフレームは１心周期内のほぼ同じ時点を表す心周期情報により空間的にアライメントされ、
前記プロセッサにより、空間的にアライメントされた１心周期内のほぼ同じ時点を表す心周期情報を有する複数のデータフレームが結合される
ことを特徴とする造影エンハンス超音波診断イメージング装置。 - 前記造影剤検出器は送信された超音波信号から所定の立方基本波での１次情報を取得するためのフィルタを含み、前記Ｂモード検出器は送信された超音波信号の所定のサブセットに対応する組織情報を検出するように構成されている、請求項１７記載の造影エンハンス超音波診断イメージング装置。
- 前記ＥＣＧ入力端子にＥＣＧ検出器（２４）が接続されており、さらに、前記造影剤検出器から出力されたフレームを結合して画像を形成および表示するディスプレイが設けられている、請求項１７記載の造影エンハンス超音波診断イメージング装置。
- 前記プロセッサは、前記造影剤検出器から出力され空間的にアライメントされたデータフレームを最大強度の投影、平均またはこれらの組み合わせによって結合するように構成されている、請求項１７記載の造影エンハンス超音波診断イメージング装置。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US12/011,178 US20090187106A1 (en) | 2008-01-23 | 2008-01-23 | Synchronized combining for contrast agent enhanced medical diagnostic ultrasound imaging |
| US12/011,178 | 2008-01-23 |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| JP2009172381A true JP2009172381A (ja) | 2009-08-06 |
| JP2009172381A5 JP2009172381A5 (ja) | 2012-03-15 |
| JP5721311B2 JP5721311B2 (ja) | 2015-05-20 |
Family
ID=40303707
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2009012590A Expired - Fee Related JP5721311B2 (ja) | 2008-01-23 | 2009-01-23 | 造影エンハンス超音波診断イメージング方法、コンピュータ読み取り可能な記憶媒体および造影エンハンス超音波診断イメージング装置 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US20090187106A1 (ja) |
| EP (1) | EP2082688A1 (ja) |
| JP (1) | JP5721311B2 (ja) |
| KR (1) | KR101599785B1 (ja) |
| CN (1) | CN101491450B (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2020522335A (ja) * | 2017-06-08 | 2020-07-30 | コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． | 超音波撮像方法 |
Families Citing this family (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2326249B1 (en) * | 2008-08-13 | 2017-01-18 | Koninklijke Philips N.V. | Dynamical visualization of coronary vessels and myocardial perfusion information |
| US20110144495A1 (en) * | 2009-12-14 | 2011-06-16 | Siemens Medical Solutions Usa, Inc. | Perfusion Imaging of a Volume in Medical Diagnostic Ultrasound |
| US9066674B2 (en) * | 2010-03-17 | 2015-06-30 | Siemens Medical Solutions Usa, Inc. | Motion synchronized destruction for three-dimensional reperfusion mapping in medical diagnostic ultrasound imaging |
| KR101120855B1 (ko) * | 2010-08-31 | 2012-03-16 | 삼성메디슨 주식회사 | 대상체의 움직임 주기에 기초하여 초음파 공간 합성 영상을 제공하는 초음파 시스템 및 방법 |
| US20120116218A1 (en) * | 2010-11-10 | 2012-05-10 | Jennifer Martin | Method and system for displaying ultrasound data |
| US20130231655A1 (en) * | 2010-11-18 | 2013-09-05 | Koninklijke Philips Electronics N.V. | Interference reduction and signal to noise ratio improvement for ultrasound cardiac ablation monitoring |
| CN104203108B (zh) * | 2012-03-23 | 2017-05-03 | 皇家飞利浦有限公司 | 用于对周期性移动的对象成像的成像系统 |
| JP2015521876A (ja) * | 2012-06-28 | 2015-08-03 | ビー−ケー メディカル エーピーエス | 超音波イメージング |
| US9974506B2 (en) * | 2013-11-05 | 2018-05-22 | International Business Machines Corporation | Associating coronary angiography image annotations with syntax scores for assessment of coronary artery disease |
| KR101525433B1 (ko) * | 2014-01-09 | 2015-06-09 | 서강대학교산학협력단 | 병소 상태 모니터링 방법 및 의료 영상 시스템 |
| EP3145413A1 (en) * | 2014-05-23 | 2017-03-29 | Koninklijke Philips N.V. | Motion gated-ultrasound thermometry using adaptive frame selection. |
| DE102014009439B4 (de) * | 2014-06-25 | 2018-05-30 | Drägerwerk AG & Co. KGaA | Vorrichtung und Verfahren zur Verarbeitung von tomografischen Daten |
| EP3192053B1 (en) | 2014-09-11 | 2021-03-31 | Koninklijke Philips N.V. | Quality metric for multi-beat echocardiographic acquisitions for immediate user feedback |
| CN104463844B (zh) * | 2014-11-03 | 2017-05-03 | 中国人民解放军第三军医大学第三附属医院 | 基于brox光流法的超声造影图像血管灌注区提取方法 |
| KR102342210B1 (ko) * | 2014-11-26 | 2021-12-28 | 삼성전자주식회사 | 프로브, 초음파 영상 장치 및 초음파 영상 장치의 제어 방법 |
| US20170143311A1 (en) * | 2015-11-12 | 2017-05-25 | Edan Instruments, Inc. | Systems and methods for pulsed wave predictive processing |
| WO2019075621A1 (zh) * | 2017-10-16 | 2019-04-25 | 北京深迈瑞医疗电子技术研究院有限公司 | 超声成像设备、系统及其超声造影成像的图像增强方法 |
| CN108983201B (zh) * | 2018-07-23 | 2022-06-17 | 南京航空航天大学 | 一种基于概率约束的稳健宽带波束形成器设计方法 |
| WO2022049299A1 (en) * | 2020-09-07 | 2022-03-10 | Koninklijke Philips N.V. | Systems and methods for measuring cardiac stiffness |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6659953B1 (en) * | 2002-09-20 | 2003-12-09 | Acuson Corporation | Morphing diagnostic ultrasound images for perfusion assessment |
| JP2004000620A (ja) * | 2003-05-19 | 2004-01-08 | Toshiba Corp | 超音波診断装置 |
| JP2004321688A (ja) * | 2003-04-28 | 2004-11-18 | Toshiba Corp | 超音波診断装置及び画像処理装置 |
| US20050033179A1 (en) * | 2003-07-25 | 2005-02-10 | Gardner Edward A. | Phase selection for cardiac contrast assessment |
| JP2005512650A (ja) * | 2001-12-18 | 2005-05-12 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 組織潅流及び経時的に変化するその他のパラメーターを表示するための超音波画像化システム及び方法 |
| JP2007090075A (ja) * | 2005-09-27 | 2007-04-12 | Siemens Medical Solutions Usa Inc | コントラスト剤を超音波によって特徴付ける方法およびコントラスト剤を超音波によって特徴付けるためにプログラムされたプロセッサによって実行可能な命令を表すデータが記憶されたコンピュータ読み出し可能記憶媒体 |
| JP2007330764A (ja) * | 2006-01-10 | 2007-12-27 | Toshiba Corp | 超音波診断装置及び超音波画像生成方法 |
| JP2009005888A (ja) * | 2007-06-28 | 2009-01-15 | Ge Medical Systems Global Technology Co Llc | 超音波撮像装置 |
Family Cites Families (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5724976A (en) * | 1994-12-28 | 1998-03-10 | Kabushiki Kaisha Toshiba | Ultrasound imaging preferable to ultrasound contrast echography |
| JP4377495B2 (ja) * | 1999-10-29 | 2009-12-02 | 株式会社東芝 | 超音波診断装置 |
| US6447450B1 (en) * | 1999-11-02 | 2002-09-10 | Ge Medical Systems Global Technology Company, Llc | ECG gated ultrasonic image compounding |
| US6494841B1 (en) * | 2000-02-29 | 2002-12-17 | Acuson Corporation | Medical diagnostic ultrasound system using contrast pulse sequence imaging |
| US6602195B1 (en) * | 2000-08-30 | 2003-08-05 | Acuson Corporation | Medical ultrasonic imaging pulse transmission method |
| US6682482B1 (en) * | 2000-08-30 | 2004-01-27 | Acuson Corporation | Medical ultrasonic imaging pulse transmission method |
| US6464643B1 (en) * | 2000-10-06 | 2002-10-15 | Koninklijke Philips Electronics N.V. | Contrast imaging with motion correction |
| US6638228B1 (en) * | 2002-04-26 | 2003-10-28 | Koninklijke Philips Electronics N.V. | Contrast-agent enhanced color-flow imaging |
| US6632177B1 (en) * | 2002-05-01 | 2003-10-14 | Acuson Corporation | Dual process ultrasound contrast agent imaging |
| US6780152B2 (en) * | 2002-06-26 | 2004-08-24 | Acuson Corporation | Method and apparatus for ultrasound imaging of the heart |
| US6673017B1 (en) * | 2002-08-28 | 2004-01-06 | Acuson Corporation | Temporal resolution method and systems for ultrasound imaging |
| US7004906B1 (en) * | 2004-07-26 | 2006-02-28 | Siemens Medical Solutions Usa, Inc. | Contrast agent imaging with agent specific ultrasound detection |
| WO2007015199A2 (en) * | 2005-08-04 | 2007-02-08 | Koninklijke Philips Electronics, N.V. | Motion compensation in functional imaging |
| WO2007080895A1 (ja) * | 2006-01-10 | 2007-07-19 | Kabushiki Kaisha Toshiba | 超音波診断装置及び超音波画像生成方法 |
| US20080081998A1 (en) * | 2006-10-03 | 2008-04-03 | General Electric Company | System and method for three-dimensional and four-dimensional contrast imaging |
| US20080214934A1 (en) * | 2007-03-02 | 2008-09-04 | Siemens Medical Solutions Usa, Inc. | Inter-frame processing for contrast agent enhanced medical diagnostic ultrasound imaging |
| US8057394B2 (en) * | 2007-06-30 | 2011-11-15 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Ultrasound image processing to render three-dimensional images from two-dimensional images |
-
2008
- 2008-01-23 US US12/011,178 patent/US20090187106A1/en not_active Abandoned
- 2008-11-24 EP EP08169745A patent/EP2082688A1/en not_active Ceased
-
2009
- 2009-01-23 CN CN200910006114.6A patent/CN101491450B/zh not_active Expired - Fee Related
- 2009-01-23 KR KR1020090005823A patent/KR101599785B1/ko not_active IP Right Cessation
- 2009-01-23 JP JP2009012590A patent/JP5721311B2/ja not_active Expired - Fee Related
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005512650A (ja) * | 2001-12-18 | 2005-05-12 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 組織潅流及び経時的に変化するその他のパラメーターを表示するための超音波画像化システム及び方法 |
| US6659953B1 (en) * | 2002-09-20 | 2003-12-09 | Acuson Corporation | Morphing diagnostic ultrasound images for perfusion assessment |
| JP2004321688A (ja) * | 2003-04-28 | 2004-11-18 | Toshiba Corp | 超音波診断装置及び画像処理装置 |
| JP2004000620A (ja) * | 2003-05-19 | 2004-01-08 | Toshiba Corp | 超音波診断装置 |
| US20050033179A1 (en) * | 2003-07-25 | 2005-02-10 | Gardner Edward A. | Phase selection for cardiac contrast assessment |
| JP2007090075A (ja) * | 2005-09-27 | 2007-04-12 | Siemens Medical Solutions Usa Inc | コントラスト剤を超音波によって特徴付ける方法およびコントラスト剤を超音波によって特徴付けるためにプログラムされたプロセッサによって実行可能な命令を表すデータが記憶されたコンピュータ読み出し可能記憶媒体 |
| JP2007330764A (ja) * | 2006-01-10 | 2007-12-27 | Toshiba Corp | 超音波診断装置及び超音波画像生成方法 |
| JP2009005888A (ja) * | 2007-06-28 | 2009-01-15 | Ge Medical Systems Global Technology Co Llc | 超音波撮像装置 |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2020522335A (ja) * | 2017-06-08 | 2020-07-30 | コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． | 超音波撮像方法 |
| JP7232199B2 (ja) | 2017-06-08 | 2023-03-02 | コーニンクレッカ フィリップス エヌ ヴェ | 超音波撮像方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101491450A (zh) | 2009-07-29 |
| US20090187106A1 (en) | 2009-07-23 |
| EP2082688A1 (en) | 2009-07-29 |
| CN101491450B (zh) | 2014-07-23 |
| JP5721311B2 (ja) | 2015-05-20 |
| KR101599785B1 (ko) | 2016-03-04 |
| KR20090081344A (ko) | 2009-07-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5721311B2 (ja) | 造影エンハンス超音波診断イメージング方法、コンピュータ読み取り可能な記憶媒体および造影エンハンス超音波診断イメージング装置 | |
| US20080214934A1 (en) | Inter-frame processing for contrast agent enhanced medical diagnostic ultrasound imaging | |
| JP4932984B2 (ja) | 超音波撮像において組織変形の実時間計算および表示を実現する方法 | |
| US8137275B2 (en) | Tissue complex modulus and/or viscosity ultrasound imaging | |
| US10390796B2 (en) | Motion correction in three-dimensional elasticity ultrasound imaging | |
| US20110144495A1 (en) | Perfusion Imaging of a Volume in Medical Diagnostic Ultrasound | |
| US20050288589A1 (en) | Surface model parametric ultrasound imaging | |
| KR101614799B1 (ko) | 의료 진단 초음파 이미징에서의 3차원 변환 매핑을 위한 모션 동기화된 파괴 | |
| JP6385992B2 (ja) | 超音波カラーフローにおけるスパークルアーチファクトの検出 | |
| JP2007518512A (ja) | 心筋灌流を表示するための画像分割 | |
| US20220079551A1 (en) | Methods and systems for valve regurgitation assessment | |
| US8956301B2 (en) | Optimization of lines per second for medical diagnostic ultrasound contrast agent imaging | |
| JP4764209B2 (ja) | 超音波信号解析装置、超音波信号解析方法、超音波解析プログラム、超音波診断装置、及び超音波診断装置の制御方法 | |
| US20090204003A1 (en) | Tracking selection for medical diagnostic ultrasound imaging | |
| US20180256133A1 (en) | Motion artifact suppression for three-dimensional parametric ultrasound imaging | |
| US9320498B2 (en) | Twinkle artifact suppression in ultrasound color flow | |
| US20220079550A1 (en) | Methods and systems for monitoring a function of a heart | |
| EP4132364B1 (en) | Methods and systems for obtaining a 3d vector flow field | |
| US20240032892A1 (en) | Boundary detection in ultrasound data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20101227 |
|
| RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20101228 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20120123 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20120131 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20130530 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20130805 |
|
| A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20131101 |
|
| A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20131107 |
|
| A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20131202 |
|
| A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20131205 |
|
| A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20131227 |
|
| A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20140109 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140205 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20140407 |
|
| A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20140704 |
|
| A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20140709 |
|
| A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20140807 |
|
| A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20140812 |
|
| 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: 20150223 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20150324 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 5721311 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| LAPS | Cancellation because of no payment of annual fees |