JP6826103B2 - 適応時間処理による微小脈管超音波造影イメージング - Google Patents
適応時間処理による微小脈管超音波造影イメージング Download PDFInfo
- Publication number
- JP6826103B2 JP6826103B2 JP2018507699A JP2018507699A JP6826103B2 JP 6826103 B2 JP6826103 B2 JP 6826103B2 JP 2018507699 A JP2018507699 A JP 2018507699A JP 2018507699 A JP2018507699 A JP 2018507699A JP 6826103 B2 JP6826103 B2 JP 6826103B2
- Authority
- JP
- Japan
- Prior art keywords
- contrast
- time
- image
- diagnostic imaging
- ultrasonic diagnostic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims description 16
- 238000002604 ultrasonography Methods 0.000 title description 26
- 238000012545 processing Methods 0.000 title description 6
- 230000003044 adaptive effect Effects 0.000 title 1
- 239000002872 contrast media Substances 0.000 claims description 38
- 238000012935 Averaging Methods 0.000 claims description 31
- 238000000034 method Methods 0.000 claims description 26
- 238000002059 diagnostic imaging Methods 0.000 claims description 19
- 230000008569 process Effects 0.000 claims description 14
- 238000001514 detection method Methods 0.000 claims description 11
- 239000000523 sample Substances 0.000 claims description 11
- 230000004044 response Effects 0.000 claims description 9
- 230000001427 coherent effect Effects 0.000 claims description 7
- 230000007423 decrease Effects 0.000 claims description 5
- 238000000926 separation method Methods 0.000 claims description 2
- 230000010412 perfusion Effects 0.000 description 14
- 230000033001 locomotion Effects 0.000 description 12
- 230000006870 function Effects 0.000 description 9
- 230000005540 biological transmission Effects 0.000 description 8
- 230000017531 blood circulation Effects 0.000 description 7
- 230000000694 effects Effects 0.000 description 7
- 229940039231 contrast media Drugs 0.000 description 6
- 230000003902 lesion Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 230000002792 vascular Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000002607 contrast-enhanced ultrasound Methods 0.000 description 3
- 230000003247 decreasing effect Effects 0.000 description 3
- 210000004088 microvessel Anatomy 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 210000000481 breast Anatomy 0.000 description 2
- 239000003795 chemical substances by application Substances 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000004941 influx Effects 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000001960 triggered effect Effects 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000000601 blood cell Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000002512 chemotherapy Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 239000002961 echo contrast media Substances 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000010894 electron beam technology Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000001747 exhibiting effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002601 intratumoral effect Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000001991 pathophysiological effect Effects 0.000 description 1
- 230000002688 persistence Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
- 210000005166 vasculature Anatomy 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/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
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
-
- 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/0833—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
- A61B8/085—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures for locating body or organic structures, e.g. tumours, calculi, blood vessels, nodules
-
- 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/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
-
- 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/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
- G01S7/52038—Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target
- G01S7/52039—Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target exploiting the non-linear response of a contrast enhancer, e.g. a contrast agent
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
-
- 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/0825—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the breast, e.g. mammography
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Radiology & Medical Imaging (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Molecular Biology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Hematology (AREA)
- Nonlinear Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Computer Networks & Wireless Communication (AREA)
- Primary Health Care (AREA)
- Epidemiology (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Vascular Medicine (AREA)
- Physiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Description
Nicolas ElifeらによるUS2010 / 060644A1、Marianne Gauthierらによるデコンボリューション法を用いたダイナミックコントラスト強調音響分析による定量的灌流パラメータの評価、Takagi KazuyaによるWO2014 / 061258A1及びThomas Patrice Arsene GauthierらによるUS2012 / 0253190A1は、時間 - 強度曲線などの灌流曲線又はパラメータの定量化及び/又はコントラスト強調超音波画像からの一つ又はそれより多くの特徴の抽出に関連するが、コントラスト強調超音波画像を提供することの改善に関するものはない。
領域の超音波画像は、空間位置に従って構成される領域の超音波データ値を有するものとして知られている。換言すれば、超音波画像の各ピクセル又はボクセルは、対応する空間位置の超音波データ値を表す。 Bモード超音波画像、ドップラー超音波画像、コントラスト超音波画像(コントラスト画像又はコントラスト強調画像とも呼ばれる)など、様々なタイプの超音波画像がある。様々なタイプの超音波画像における超音波データ値は、様々な物理的パラメータを表す。例えば、Bモードの超音波画像におけるピクセル値は、組織によって反射される超音波エコー信号の強度を表す。例えば、ドップラー超音波画像におけるピクセル値は、血流などの組織を移動させることによって引き起こされるドップラー効果を表す。例えば、コントラスト画像のピクセル値は造影剤の量を表す。いくつかの実施形態では、コントラスト時間プロセッサは、高調波エコー信号を処理して異なって処理されるコントラストデータ値を生成する時間平均化計算器及び最大強度検出器を含むことができる。最大強度検出は、コントラスト強度の急激な変化をより良く反映することができ、時間平均化処理は、ピクセル値のノイズ及び動きの効果を低減することができる。いくつかの他の実施形態では、最大強度検出器及び時間平均化計算器のいずれか又は両方は、他の既存又は将来開発される計算機又は検出器によって置き換えられることができる。
いくつかの態様では、MI重み付け回路は、中間段階の間に減少重み付け関数を利用するように構成され、TA重み付け回路は、中間段階の間に増加重み付け関数を利用するように構成される。コントラスト時間プロセッサは、時間 - 強度曲線トリガ回路によって生成されるT-I曲線のパラメータを記憶するように構成されるT-I履歴バッファを含むこともできる。いくつかの態様では、時間 - 強度曲線トリガ回路は、第1のコントラストウォッシュイン期間中に記憶されるパラメータを使用して、第2のコントラストウォッシュイン期間中に重み付け回路をトリガするように構成される。パラメータは、最初のウォッシュイン段階の終わりを画定する第1のパラメータt2と、後のウォッシュアウト段階の始まりを画定する第2のパラメータt3とを含むことができる。
本明細書で更に説明するように、本発明は、超音波システム及び方法を含む。いくつかの態様において、超音波システムは、マイクロプロセッサ、集積回路(例えば、FPGA)、メモリ、ハードドライブなどのコンピュータに標準的な様々な構造を含む。いくつかの実施形態において、本発明は、実行されると、本明細書に記載の様々なステップ及び機能をシステムに実行させる。例えば、本発明は、微小脈管構造のコントラスト強調イメージングのための超音波診断イメージングシステムであって、命令が実行されるとき、システムに、超音波プローブを用いて超音波ビームを送信し、超音波エコー信号を受信するステップと、コヒーレントエコー信号を生成するためにエコー信号をビーム形成するステップと、造影剤から受信される高調波エコー信号を分離するステップと、造影剤の少なくとも一つのウォッシュイン段階を含む時間にわたって異なる方法で造影剤から受信される高調波信号を処理するステップと、所望の画像フォーマットでコントラスト画像を生成するステップとを実行させる命令を有する、超音波診断イメージングシステムと、コントラスト画像を表示するスキャンコンバータに結合されるディスプレイとを含む。
本発明の他の態様によれば、画像領域における微小脈管構造の超音波コントラスト画像を生成する装置であって、造影剤から受信される高調波エコー信号に応答し、造影剤から受信される高調波エコー信号を処理して画像領域における複数の空間位置の各々に対するコントラストデータ値を生成するように構成されるコントラスト時間プロセッサと、生成されるコントラストデータ値を受信するために結合され、期間にわたってコントラスト画像の時間的シーケンスを生成するように構成されるスキャンコンバータとを有し、空間点位置に対するコントラスト値が、造影剤の少なくとも一つのウォッシュイン段階を含む期間にわたって異なる方法で生成される、装置が提供される。この装置は更に、コントラスト画像を表示するためのディスプレイに結合されることができる。これに加えて又はこれに代えて、このようなシステムは、コントラスト画像を表示するためのディスプレイを備えることができる。この装置は、造影剤から受信される高調波エコー信号を提供するための装置又はシステムに結合されることができる。追加的又は代替的に、この装置は、造影剤から受信される高調波エコー信号を提供するための装置又はシステムの一部であってもよい。提供される高調波エコー信号は、高調波エコー信号のフレームのシーケンスであり得る。提供される高調波エコー信号は、コントラストデータフレームが時間に渡って同じ方法で生成される従来の手法に従って生成されるコントラストデータフレームのシーケンスになることができる。造影剤から受信される高調波エコー信号を提供するための装置又はシステムは、超音波ビームを送信し、超音波エコー信号を受信するように構成される超音波プローブを含むことができ、又は病院情報システム(HIS)、画像保管通信システム(PACS)、電子カルテ情報システム(EMR)などのような、コントラスト強調超音波データなどのデータを記憶するように構成される何れかの情報システムになり得る。
いくつかの実施形態では、コントラストピクセルデータは、ウォッシュイン期間、強調安定化又は中間期間、及びウォッシュアウト期間の3つの期間中に異なる方法で適応的に処理される。例えば、コントラストの急速な蓄積をより良く強調するために、ウォッシュインの初期段階の間に最大強度検出が使用される。後のウォッシュアウト段階の間、時間画像データの時間平均化処理が、ノイズ及びモーションアーチファクトを低減するために使用される。中間ピーク強調段階の間に、これらの技術の両方の組み合わせが使用され、この期間中に最大強度検出から時間平均化処理に変化する。
特定の実施形態では、ウォッシュイン期間は、出現段階、増加段階及びピーク段階によって特徴付けられ、コントラストピクセルデータは、これらの期間中に異なる方法で適応的に処理される。ここでは、最大強度検出が、出現段階の間に、コントラストの急速な蓄積をより良好に強調するために使用されることができる。後のピーク段階の間に、時間画像データの時間平均化処理が、ノイズ及びモーションアーチファクトを低減するために使用される。増加段階の間に、これらの技術の両方の組み合わせが使用され、この期間中に最大強度検出から時間平均化処理に変化する。
送信される超音波エネルギーからのエコーは、T / Rスイッチ22を介して結合され、システムがデジタルビームフォーマを使用するとき、ビームフォーマの入力部においてアナログ/デジタル(「A / D」)コンバータによってデジタル化されるエコー信号を生成するアレイ14のトランスデューサ素子によって受信される。代わりにアナログビームフォーマが使用されることもできる。超音波システムの制御及びプローブ選択などのイメージングのための様々な制御設定の制御は、中央システムコントローラ(図示略)に結合され、中央システムコントローラを通じてその制御を適用する制御パネルの制御部のユーザ操作によって行われる。
アレイ14の個々のトランスデューサ素子から受信されるエコー信号は、コヒーレントエコー信号をビーム形成するためにビームフォーマ32によって遅延され、合計される。 2次元アレイによる3Dイメージングでは、米国特許第6,013,032号(Savord)及び米国特許第6,375,617号(Fraser)に記載されるように、プローブに配置されるマイクロビームフォーマとシステムメインフレームのメインビームフォーマとの間にビームフォーマを区画することが好ましい。それから、デジタルコヒーレントエコー信号は、バンドパスフィルタリング、スペックル低減、画像コントラスト強調、組織クラッタ抑圧及び動き補償などの動作を実行する信号プロセッサ34によって処理される。信号プロセッサは、受信される周波数帯域をより低い周波数レンジ又はベースバンド周波数レンジにシフトすることもできる。この実施形態では、送信周波数及び受信周波数は、ビームフォーマ32が、高調波周波数帯域などの送信帯域の周波数帯域とは異なる周波数帯域を自由に受信するように、個別に制御される。
マイクロバブルのような造影剤によるイメージングの間に受信されるエコー信号は、基本波/高調波信号セパレータ36に結合される。基本波/高調波信号セパレータ36は、好ましくは、パルス反転技術によって高調波造影剤から返されるエコーを分離し、画像位置への複数のパルスの送信から生じるエコーは、基本信号成分をキャンセルして高調波成分を強調するように加算的に結合され、基本周波数信号成分を生成するように減算的に結合される。好ましいパルス反転技術は、例えば、米国特許第6,186,950号(Averkiouら)に記載される。基本周波数f及び第2高調波周波数2fを含む信号帯域のような基本周波数及び高調波周波数は、基本周波数fにピークを有する低域通過帯域と、第2高調波周波数2fにピークを有する高域通過帯域とを用いる帯域通過フィルタリングによって分離されることもできる。
空間フォーマットでエコー信号を受信したフェーズドアレイ又は湾曲アレイプローブであってもよく、スキャンコンバータは、これらの信号を、ラスタースキャンされる表示装置のピクセルとしてBモード画像の表示に適したデカルト形式に変換する。 3次元イメージングの場合、スキャン変換される2次元画像データは、ボリュームレンダリングによって3D画像を形成するように処理されてもよい。
本発明の原理によれば、造影剤のマイクロバブルから受信される分離高調波信号は、コントラスト時間プロセッサ50に結合される。コントラスト時間プロセッサは、以下に説明される造影剤ウォッシュインの異なる段階の間に異なる方法でコントラスト信号を検出して処理する。得られたコントラストデータ値は、別個のスキャンコンバータ又は時間多重化される態様で使用されるときにBモード画像スキャン変換のために使用されるものと同じになり得るスキャンコンバータ52によってコントラスト画像を形成するように、Bモード画像と同じフォーマットに変換される。スキャンコンバータ42によって生成されるBモード画像及びスキャンコンバータ52によって生成されるコントラスト画像は、画像プロセッサ24に結合される。2つの画像は両方とも同じ表示フォーマットであるため、コントラスト画像は、造影剤によって描写される血流がBモード画像の組織構造によって囲まれるように、Bモード画像をオーバレイして表示する。画像プロセッサは、例えば、グラフィック情報を追加することによって画像を更に処理することができ、最終画像はディスプレイ装置又はモニタ26に表示される。
図2は、コントラスト時間プロセッサ50の構成及び動作を示す詳細なブロック図である。連続して受信される画像フレームは、フレームバッファ54に記憶される。注入又はインジェクトされる造影剤がイメージングされる体の領域に到達し始めると、連続する画像フレームは、フレームバッファに記憶され、イメージングされる微小脈管構造の漸進的な位置に現れるとき造影剤を捕捉する。プローブ12が患者の体に対して安定して保持されると、同じ空間位置が各画像内の同じピクセル位置に現れる。したがって、コントラストウォッシュイン中のイメージング領域における各固有の位置におけるピクセルデータの時間シーケンスは、連続画像のシーケンスのそれぞれにおける同じピクセル位置をアドレス指定することによって簡単に抽出されることができる。
本発明の原理によれば、各ピクセル位置でのコントラストデータの時間シーケンスは、最大強度検出器72及び時間平均化計算器82によって処理される。最大強度検出器は、ピクセル位置からデータ値のシーケンスを受け取り、それらを順次比較して、コントラストウォッシュインの開始から現在の時点までの最大値、最近受信したコントラスト画像値を検出する。数学的には、これは各ピクセル位置について
Pmaxx、y=P1、P2、P3 ... Pcurrent のうちのPmax
と表されることができる。ここで、x、yは画像内のピクセル位置であり、ウォッシュインの開始から現在の時間までのフレーム1,2,3,現在のピクセル値は、P1、P2、P3 ... Pcurrentである。簡単な例を与えるために、データ値が0,0,3,2,3,5,7及び10であると仮定する。このシーケンスのPmaxx、yの値は、連続的に0,0,3,3,3,5,7及び10になる。この検出プロセスはノイズに対する感度を有することが分かる。最初の2つの0がウォッシュインの開始前のピクセル値であり、これらの最初がノイズ又は動きの影響を受けて3であると仮定する。シーケンスの最大強度検出値は、3,3,3,3,3,5,7及び10になり、所望されないことに、ノイズ効果を反映している。
Ptax、y=P1、P2、P3 ... Pcurrent のうちの平均
と表される。他の実施形態では、時間平均化は、移動ウィンドウ内の加重平均又は非加重平均に等しい移動平均になることができる。
この式では、ノイズの瞬間的な影響が平均化され、それによって値のシーケンスの平均に渡って減少されることがわかる。最大強度検出器及び時間平均化計算器によって実行される動作は繰り返され、それらの出力コントラストデータ値はフレームバッファ内のコントラストデータの各新しいフレームの到達で更新される。
本発明のさらなる態様によれば、コントラスト灌流の全期間は、少なくとも2つの期間に分割される。一実施形態では、少なくとも2つの期間は、第1の期間及び第2の期間を含むことができ、コントラスト時間プロセッサは、第1の期間中の最大強度値に等しく、第2の期間中の時間平均強度値に等しい所与のピクセル位置の瞬時コントラストデータ値を生成する。例えば、コントラスト時間プロセッサは、最大強度値及び時間平均強度値のうちの1つを選択するセレクタを含むように構成することができる。他の実施形態では、少なくとも2つの期間は、第1の期間、第2の期間及び第3の期間を含むことができ、コントラスト時間プロセッサは、第1の期間中の最大強度値に等しく、第3の期間中の時間平均強度値に等しく、前記第1の期間と前記第3の期間との間の前記第2の期間中の前記2つの組み合わせに等しい所与のピクセル位置に対する瞬時コントラストデータ値を生成する。
本発明のさらなる態様によれば、コントラスト時間プロセッサは、ウォッシュインの初期段階中の最大強度値に等しく、ウォッシュインの後の段階の間の時間平均強度値に等しく、コントラスト灌流が最大である中間段階の間の2つの変化する組み合わせに等しい所与のピクセル位置に対する瞬時コントラストデータ値を生成する。これは、一つ又はそれより多くのコントラストピクセル位置に対して時間 - 強度曲線トリガ回路56によって計算される時間 - 強度(T-I)曲線を使用することによって行われる。 T-I曲線は、米国特許第5,833,613号(Averkiou他)に記載されるように、体内のある位置における造影剤の到達及び送達の共通の尺度である。時間 - 強度曲線は、灌流される組織の画像内の各点及びイメージング又は診断に使用するために抽出される各曲線の一つ又はそれより多くのパラメータに対して計算されることができる。例えば、各画像ポイントに対して選択されるパラメータは、米国特許第6,692,438号(Skyba他)に記載されるように、灌流のパラメータ画像を形成するためにグレースケールシェーディング又はカラーコーディングで表示され得る。これらのパラメータは、曲線のピーク及び勾配を含み、それぞれは組織灌流の異なる特性を示す。
によって定義されるような曲線にフィッティングされる。ここでAは曲線のピークであり、t0は造影剤のウォッシュインの開始時間であり、
は曲線の立ち上がりの勾配であり、xは造影剤の量の瞬時測定値である。これらの時間及び強度表現は、訓練される臨床医に組織が灌流される様子を示す。例えば、図3は、造影剤のウォッシュイン中に経時的に受け取ったコントラスト強度値のシーケンスを示す。 x軸は時間軸であり、y軸はピクセル値の振幅(強度)である。これらの値は、それらを結ぶ線分によって示されるように、幾分不規則に変化することがわかる。この不規則なばらつきを克服するために、曲線は、上記のような曲線適合式を使用して、曲線110によって示されるような強度データ値にフィッティングされる。これは、コントラストのウォッシュインのより滑らかな表現を提示し、これは、上に列挙されるパラメータをもたらすように分析されるフィッティング曲線である。
図5は、曲線フィッティングから得られる典型的なT-1ウォッシュイン曲線60を示す。この例では、立ち上がり又は最初のウォッシュイン段階は、63及び時間t1によって示される曲線60のピークAの20%から65及び時間t2によって示される曲線のピークの80%のレベルまでの立ち上がりに後続する持続時間である。造影剤の量がその灌流のピークの周辺にある強調段階は、時間t2における65の80%マークと、67及び時間t3におけるピークの90%への低下との間の持続時間である。コントラスト又はウォッシュアウト段階の低下は、67及び時間t3におけるピークの90%から、69及び時間t4におけるピークの30%までの持続時間である。この例では、t1-t2はウォッシュイン段階、t2-t3はエンハンスメント段階、t3-t4はウォッシュアウト段階である。本発明の原理によれば、最大強度検出器72によって生成される値は、ウォッシュインの初期段階、図の65(t2)までの期間の瞬時コントラストピクセル値として使用される。時間平均化計算器82によって生成される値は、後の段階、67(t3)に続く期間に使用される。ピーク強調段階t2乃至t3の間、両方の値の重み付け組合せが使用される。これは、各段階が達せられるときに時間 - 強度曲線トリガ回路56による重み付け変化をトリガすることによって行われ、2つの重み付け回路、最大強度(MI)重み付け回路74及び時間平均化(TA)重み付け回路84の動作を制御する。これらの重み付け回路は、図6の重み付け関数曲線410及び420によって示されるようなそれぞれの値に重み付けを適用する。この例では、曲線420はTA重み付け回路84によって使用され、曲線410はMI重み付け回路74によって使用される。
ピークが実現される前の第1の時間 - 強度曲線の計算中に時間 - 強度曲線のピークAを正確に予測することは不可能であることが理解されるであろう。この問題を解決するために、本発明の図示される実施形態は、新しいT-I曲線が計算されるまで使用される1つの曲線からのT-1パラメータを記憶するT-1履歴バッファ58も有する。一例として、患者に一定のコントラストの流れが注入されると仮定する。領域内のマイクロバブルを破壊するために、超音波の高出力フラッシュ送信が画像領域に適用される。その後、マイクロバブルの新しい流入の間、時間 - 強度曲線が計算され、生成され、時間t2(最初のウォッシュインの終了)及びt3(ピーク強調の終了)が決定され、T-I履歴バッファに記憶される。他の高出力フラッシュ送信が画像領域に適用され、コントラストの新しい流入がイメージングされる微小脈管構造内にウォッシュインされ始める。マイクロバブルが到達し、画像内のピクセル位置における強度を強調すると、重み付け回路は、先行して決定される時間t2及びt3にそれらの重み付けを変化させるようにトリガされる。同じウォッシュイン期間中、時間 - 強度曲線トリガ回路は新しい時間 - 強度曲線を計算し、現在のウォッシュイン、ウォッシュインサイクルを使用して更新される時間t2及びt3を決定し、更新値をTI履歴バッファに記憶する。したがって、各々の繰り返されるウォッシュイン、ウォッシュインサイクルは、t2及びt3スイッチオーバータイムの最新の更新セットを使用し、重み付け回路機能の変更をトリガする。
当業者であれば、本発明による超音波システムは、ハードウェア、ソフトウェア、又はその両方の組み合わせを使用して構成できることを直ちに認識するであろう。ハードウェア構成では、システムは、記載される発明を実行する回路を含むことができ、又は請求される処理を実行するように構成されるゲートを備えたFPGAのような高度なデジタル回路を使用することができる。今日の超音波システムの大部分の構成であるソフトウェア構成では、フレームバッファは、アドレス可能なメモリ位置にコントラストデータを記憶するデジタルメモリであり、図2に示す他の機能はソフトウェアサブルーチンによって実行される。たとえば、最大強度検出器は、デジタルメモリから連続的にアドレス指定されるコントラストデータ値を減算することによる比較を実行するソフトウェアルーチンによって提供される。時間平均化計算機は、コントラストデータ値のシーケンスを加算し、その合計を値の数で除算するソフトウェアルーチンによって提供される。重み付け回路機能は、時間 - 強度曲線トリガによってトリガされるメモリに記憶される係数を重み付けすることによって、検出器及び計算器ルーチンからの値を掛けることによって提供される。加算器は加算サブルーチンによって提供される。時間 - 強度曲線トリガは、最良の適合が見つかるまで異なる変数値を有する方程式を反復して実行することにより、上記のT-I式を実行するルーチンによって、T-I曲線をコントラストデータ値に適合させる。T-I履歴バッファはデジタルメモリである。加算器90に続いて、必要に応じて、出力データに利得係数値を乗算することによって、出力データは正規化されることができる。ハードウェアとソフトウェアの両方の組み合わせは、当業者には容易に明らかであろう。
Claims (15)
- 画像領域における微小脈管構造の超音波コントラスト画像を生成する装置であって、
高調波エコー信号に応答し、造影剤から受信される高調波信号を処理して、ある期間内のある時点に対する前記画像領域内の複数の空間位置の各々に対してコントラストデータ値を生成するように構成されるコントラスト時間プロセッサであって、前記画像領域内の空間位置に対する前記コントラストデータ値は前記期間にわたって異なる処理で生成され、前記期間は前記造影剤の少なくとも一つのウォッシュイン段階を含む、コントラスト時間プロセッサと、
前記生成されるコントラストデータ値を受信するように結合され、所望の画像フォーマットにおける前記時点に対する前記複数の空間位置に対して前記生成されるコントラストデータ値を構成することによって前記時点に対するコントラスト画像を生成するように構成されるスキャンコンバータと
を有する、装置。 - 画像領域内の微小脈管構造のコントラスト強調イメージングのための超音波診断イメージングシステムであって、
超音波ビームを送信し、超音波エコー信号を受信するように構成される超音波プローブと、
コヒーレントエコー信号を生成するために前記エコー信号をビーム形成するように構成されるビームフォーマと、
前記コヒーレントエコー信号に応答し、造影剤から受信される分離高調波エコー信号を生成するように構成される信号セパレータと、
請求項1に記載の超音波コントラスト画像を生成する装置と、
前記スキャンコンバータに結合され、コントラスト画像を表示するように構成されるディスプレイと
を有する、超音波診断イメージングシステム。 - 前記コントラスト時間プロセッサは、前記高調波エコー信号を処理して異なって処理されるコントラストデータ値を生成するように構成される時間平均化計算器及び最大強度検出器を更に有する、請求項2に記載の超音波診断イメージングシステム。
- 前記コントラスト時間プロセッサは、前記最大強度検出器によって生成されるコントラストデータ値を重み付けするように構成される最大強度重み付け回路と、前記時間平均化計算器によって生成されるコントラストデータ値を重み付けするように構成される時間平均化重み付け回路とを更に有する、請求項3に記載の超音波診断イメージングシステム。
- 前記コントラスト時間プロセッサは、前記重み付けされるコントラストデータ値を合計するように構成される加算器を更に有する、請求項4に記載の超音波診断イメージングシステム。
- 前記信号セパレータから受信される高調波エコー信号を記憶するように構成されるフレームバッファを更に有する、請求項5に記載の超音波診断イメージングシステム。
- 前記コントラスト時間プロセッサは、前記フレームバッファからの高調波エコー信号に応答し、コントラストウォッシュインの連続する段階の時間 - 強度曲線を生成するように構成される、時間 - 強度曲線トリガ回路を更に有する、請求項6に記載の超音波診断イメージングシステム。
- 前記時間 - 強度曲線は、初期ウォッシュイン段階及び後のウォッシュアウト段階を画定し、
前記時間 - 強度曲線は、前記最大強度重み付け回路及び前記時間平均化重み付け回路をトリガして、前記初期段階の間の最大強度コントラスト値及び前記後の段階の間の時間平均化コントラスト値を生成する、
請求項7に記載の超音波診断イメージングシステム。 - 前記時間 - 強度曲線は、前記初期段階に後続する中間段階を更に画定し、
前記時間 - 強度曲線は、前記最大強度重み付け回路及び前記時間平均化重み付け回路をトリガして、前記中間段階の間に最大強度コントラスト値及び時間平均化コントラスト値を生成し、
前記加算器は、前記中間段階の間に前記最大強度コントラスト値及び前記時間平均化コントラスト値の組み合わせを生成するように構成される、
請求項8に記載の超音波診断イメージングシステム。 - 前記最大強度重み付け回路は、前記中間段階の間に減少重み付け関数を利用するように構成され、
前記時間平均化重み付け関数は、前記中間段階の間に増加重み付け関数を利用するように構成される、
請求項9に記載の超音波診断イメージングシステム。 - 前記コントラスト時間プロセッサは、前記時間 - 強度曲線トリガ回路によって生成される時間 - 強度曲線のパラメータを記憶するように構成される時間 - 強度履歴バッファを更に有する、請求項7に記載の超音波診断イメージングシステム。
- 前記時間 - 強度曲線トリガ回路は、第1のコントラストウォッシュイン期間の間に記憶されるパラメータを利用して、第2のコントラストウォッシュイン期間の間に前記重み付け回路をトリガするように構成される、請求項11に記載の超音波診断イメージングシステム。
- 前記パラメータは、初期ウォッシュイン段階の終了を画定する第1のパラメータt2と、後のウォッシュアウト段階の開始を画定する第2のパラメータt3とを有する、請求項12に記載の超音波診断イメージングシステム。
- 前記信号セパレータは、分離基本周波数エコー信号をするように更に構成され、前記システムは、
基本周波数エコー信号に応答し、検出Bモード信号を生成するように構成されるBモード検出器と、
前記検出Bモード信号に応答し、前記所望の画像フォーマットでBモード画像を生成するように構成されるスキャンコンバータと、
前記コントラスト画像及び前記Bモード画像に応答し、前記Bモード画像及び前記コントラスト画像の組み合わせを有する画像を生成するように構成される画像プロセッサと
を更に有する、請求項2に記載の超音波診断イメージングシステム。 - 前記時間 - 強度曲線は、出現段階、増加段階、及びピーク段階を画定し、
前記時間 - 強度曲線は、前記最大強度重み付け回路及び前記時間平均化重み付け回路をトリガして、前記出現段階の間の最大強度コントラスト値、前記ピーク段階の間の時間平均化コントラスト値、及び前記増加段階の間の最大強度コントラスト値と時間平均化コントラスト値との組み合わせを生成する、
請求項7に記載の超音波診断イメージングシステム。
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2015087789 | 2015-08-21 | ||
CNPCT/CN2015/087789 | 2015-08-21 | ||
EP15189930.9 | 2015-10-15 | ||
EP15189930 | 2015-10-15 | ||
PCT/EP2016/069712 WO2017032715A1 (en) | 2015-08-21 | 2016-08-19 | Micro vascular ultrasonic contrast imaging by adaptive temporal processing |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2018523531A JP2018523531A (ja) | 2018-08-23 |
JP2018523531A5 JP2018523531A5 (ja) | 2019-09-26 |
JP6826103B2 true JP6826103B2 (ja) | 2021-02-03 |
Family
ID=56842795
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2018507699A Active JP6826103B2 (ja) | 2015-08-21 | 2016-08-19 | 適応時間処理による微小脈管超音波造影イメージング |
Country Status (7)
Country | Link |
---|---|
US (1) | US11000262B2 (ja) |
EP (1) | EP3337403B1 (ja) |
JP (1) | JP6826103B2 (ja) |
CN (1) | CN107920801B (ja) |
BR (1) | BR112018003058B1 (ja) |
RU (1) | RU2720102C2 (ja) |
WO (1) | WO2017032715A1 (ja) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114245725A (zh) * | 2019-08-05 | 2022-03-25 | 皇家飞利浦有限公司 | 在洗入、洗出期间改变系统操作的对比度增强超声成像 |
US20210228187A1 (en) * | 2020-01-29 | 2021-07-29 | GE Precision Healthcare LLC | System and methods for contrast-enhanced ultrasound imaging |
JP2022018932A (ja) * | 2020-07-16 | 2022-01-27 | コニカミノルタ株式会社 | 超音波診断装置、超音波信号処理方法、及びプログラム |
CN115953401A (zh) * | 2023-03-13 | 2023-04-11 | 深圳鲲为科技有限公司 | 超声造影视频数据分析方法、装置、电子设备及存储介质 |
Family Cites Families (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3720897B2 (ja) | 1996-03-01 | 2005-11-30 | 株式会社東芝 | 動画表示方法およびコンピュータシステム |
US6013032A (en) | 1998-03-13 | 2000-01-11 | Hewlett-Packard Company | Beamforming methods and apparatus for three-dimensional ultrasound imaging using two-dimensional transducer array |
JP4130004B2 (ja) * | 1998-03-13 | 2008-08-06 | ジーイー横河メディカルシステム株式会社 | 超音波撮像装置 |
JP4408988B2 (ja) | 1999-05-31 | 2010-02-03 | 株式会社東芝 | 超音波診断装置 |
US6186950B1 (en) | 1999-11-04 | 2001-02-13 | Atl Ultrasound | Ultrasonic pulse inversion harmonic separation with reduced motional effects |
US6468216B1 (en) | 2000-08-24 | 2002-10-22 | Kininklijke Philips Electronics N.V. | Ultrasonic diagnostic imaging of the coronary arteries |
DE60212917T2 (de) * | 2001-10-16 | 2007-03-01 | Kabushiki Kaisha Toshiba | Vorrichtung zur Berechnung eines Index von örtlichen Blutflüssen |
US6692438B2 (en) | 2001-12-18 | 2004-02-17 | Koninklijke Philips Electronics Nv | Ultrasonic imaging system and method for displaying tissue perfusion and other parameters varying with time |
US6676606B2 (en) | 2002-06-11 | 2004-01-13 | Koninklijke Philips Electronics N.V. | Ultrasonic diagnostic micro-vascular imaging |
US20040077952A1 (en) * | 2002-10-21 | 2004-04-22 | Rafter Patrick G. | System and method for improved diagnostic image displays |
JP3748848B2 (ja) * | 2002-11-11 | 2006-02-22 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 超音波診断装置 |
US20040120559A1 (en) * | 2002-12-20 | 2004-06-24 | Hall Anne Lindsay | Methods and apparatus for contrast agent time intensity curve analyses |
US7591788B2 (en) * | 2003-08-19 | 2009-09-22 | Siemens Medical Solutions Usa, Inc. | Adaptive contrast agent medical imaging |
EP2099364B8 (en) * | 2006-12-21 | 2019-04-03 | Bracco Suisse SA | Detection of the detachment of immobilized contrast agent in medical imaging applications |
WO2008053268A1 (en) * | 2006-12-21 | 2008-05-08 | Institut Gustave Roussy (Igr) | Method and system for quantification of tumoral vascularization |
JP5263867B2 (ja) | 2007-10-15 | 2013-08-14 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 超音波撮像装置 |
US20090204003A1 (en) | 2008-02-07 | 2009-08-13 | Guracar Ismayil M | Tracking selection for medical diagnostic ultrasound imaging |
CN102573647A (zh) * | 2009-10-01 | 2012-07-11 | 皇家飞利浦电子股份有限公司 | 用于监测肝脏治疗的肝脏血流的对比增强超声评估 |
US9131923B2 (en) * | 2010-08-31 | 2015-09-15 | Hitachi Medical Corporation | Ultrasonic diagnostic device and ultrasonic image display method |
JP6295956B2 (ja) * | 2012-10-19 | 2018-03-20 | コニカミノルタ株式会社 | 超音波診断装置、及び超音波診断装置の制御方法 |
CN103876776B (zh) | 2012-12-24 | 2017-09-01 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声造影成像方法及装置 |
-
2016
- 2016-08-19 US US15/753,198 patent/US11000262B2/en active Active
- 2016-08-19 BR BR112018003058-9A patent/BR112018003058B1/pt active IP Right Grant
- 2016-08-19 CN CN201680048022.2A patent/CN107920801B/zh active Active
- 2016-08-19 RU RU2018106243A patent/RU2720102C2/ru active
- 2016-08-19 JP JP2018507699A patent/JP6826103B2/ja active Active
- 2016-08-19 WO PCT/EP2016/069712 patent/WO2017032715A1/en active Application Filing
- 2016-08-19 EP EP16757848.3A patent/EP3337403B1/en active Active
Also Published As
Publication number | Publication date |
---|---|
BR112018003058A2 (ja) | 2018-10-02 |
WO2017032715A1 (en) | 2017-03-02 |
US20190015074A1 (en) | 2019-01-17 |
EP3337403B1 (en) | 2020-11-25 |
CN107920801A (zh) | 2018-04-17 |
CN107920801B (zh) | 2021-05-28 |
RU2018106243A3 (ja) | 2019-09-23 |
RU2018106243A (ru) | 2019-09-23 |
US11000262B2 (en) | 2021-05-11 |
EP3337403A1 (en) | 2018-06-27 |
JP2018523531A (ja) | 2018-08-23 |
BR112018003058B1 (pt) | 2022-11-29 |
RU2720102C2 (ru) | 2020-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6368277B1 (en) | Dynamic measurement of parameters within a sequence of images | |
US8523776B2 (en) | Ultrasonic doppler imaging apparatus and method with blood velocity waveform processing | |
US11969286B2 (en) | Systems and methods for automatic detection and visualization of turbulent blood flow using vector flow data | |
JP6826103B2 (ja) | 適応時間処理による微小脈管超音波造影イメージング | |
JP7232195B2 (ja) | 血管内の壁せん断応力の同時視覚化及び定量化のためのシステム及び方法 | |
EP3871007B1 (en) | Adaptive ultrasound flow imaging | |
US8956301B2 (en) | Optimization of lines per second for medical diagnostic ultrasound contrast agent imaging | |
US20230329669A1 (en) | System and method for concurrent visualization and quantification of blood flow using ultrasound vector flow imaging | |
JP2019534103A (ja) | 造影剤流の肝灌流を特徴付けるシステム及び方法 | |
US10499883B2 (en) | Methods and systems for spatial color flow for diagnostic medical imaging | |
EP3527141B1 (en) | Method of displaying doppler image and ultrasound diagnosis apparatus for performing the method | |
JP2022543539A (ja) | ウォッシュイン、ウォッシュアウト時のシステム動作の変更を伴うコントラスト強調超音波撮像 | |
CN114727807A (zh) | 减少超声图像中的混响伪影以及相关联的设备、系统和方法 | |
US20230368376A1 (en) | Ultrasound time-series data processing device and ultrasound time-series data processing program | |
JP2011024828A (ja) | 超音波診断装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20180711 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20190815 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20190815 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20200526 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20200807 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20201126 |
|
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: 20201222 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210114 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6826103 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |