JP2002541963A - Cine高調波位相(harp)磁気共鳴画像法を用いた心臓動きのトラッキング - Google Patents
Cine高調波位相(harp)磁気共鳴画像法を用いた心臓動きのトラッキングInfo
- Publication number
- JP2002541963A JP2002541963A JP2000613339A JP2000613339A JP2002541963A JP 2002541963 A JP2002541963 A JP 2002541963A JP 2000613339 A JP2000613339 A JP 2000613339A JP 2000613339 A JP2000613339 A JP 2000613339A JP 2002541963 A JP2002541963 A JP 2002541963A
- Authority
- JP
- Japan
- Prior art keywords
- image
- harp
- tagged
- tracking
- distortion
- 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
- 230000033001 locomotion Effects 0.000 title claims abstract description 131
- 210000002216 heart Anatomy 0.000 title claims abstract description 53
- 238000002595 magnetic resonance imaging Methods 0.000 title claims description 43
- 238000000034 method Methods 0.000 claims abstract description 133
- 230000003595 spectral effect Effects 0.000 claims abstract description 43
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 claims abstract description 21
- 230000008569 process Effects 0.000 claims abstract description 12
- 210000004165 myocardium Anatomy 0.000 claims description 27
- 210000005240 left ventricle Anatomy 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 7
- 230000006872 improvement Effects 0.000 claims description 6
- 230000007423 decrease Effects 0.000 claims description 5
- 230000000977 initiatory effect Effects 0.000 claims 1
- 238000003384 imaging method Methods 0.000 description 34
- 230000002107 myocardial effect Effects 0.000 description 31
- 238000013459 approach Methods 0.000 description 30
- 238000004422 calculation algorithm Methods 0.000 description 23
- 238000006073 displacement reaction Methods 0.000 description 22
- 238000004458 analytical method Methods 0.000 description 19
- 230000003287 optical effect Effects 0.000 description 18
- 238000004364 calculation method Methods 0.000 description 16
- 230000000747 cardiac effect Effects 0.000 description 16
- 230000006870 function Effects 0.000 description 15
- 239000000463 material Substances 0.000 description 15
- 239000013598 vector Substances 0.000 description 13
- 210000001519 tissue Anatomy 0.000 description 12
- 230000002861 ventricular Effects 0.000 description 11
- 241000282472 Canis lupus familiaris Species 0.000 description 10
- 230000008901 benefit Effects 0.000 description 10
- 238000012545 processing Methods 0.000 description 10
- 230000008602 contraction Effects 0.000 description 9
- 210000001174 endocardium Anatomy 0.000 description 9
- 230000005415 magnetization Effects 0.000 description 9
- 238000004904 shortening Methods 0.000 description 9
- 238000010009 beating Methods 0.000 description 8
- 238000002474 experimental method Methods 0.000 description 8
- 239000011159 matrix material Substances 0.000 description 8
- 230000002159 abnormal effect Effects 0.000 description 6
- 230000004913 activation Effects 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 6
- 230000002123 temporal effect Effects 0.000 description 6
- JRWZLRBJNMZMFE-UHFFFAOYSA-N Dobutamine Chemical compound C=1C=C(O)C(O)=CC=1CCNC(C)CCC1=CC=C(O)C=C1 JRWZLRBJNMZMFE-UHFFFAOYSA-N 0.000 description 5
- 239000002131 composite material Substances 0.000 description 5
- 229960001089 dobutamine Drugs 0.000 description 5
- 230000014509 gene expression Effects 0.000 description 5
- 238000001514 detection method Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000012800 visualization Methods 0.000 description 4
- 208000012661 Dyskinesia Diseases 0.000 description 3
- 210000004369 blood Anatomy 0.000 description 3
- 239000008280 blood Substances 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000002059 diagnostic imaging Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 238000010191 image analysis Methods 0.000 description 3
- 208000031225 myocardial ischemia Diseases 0.000 description 3
- 238000003672 processing method Methods 0.000 description 3
- 230000000541 pulsatile effect Effects 0.000 description 3
- 241000282412 Homo Species 0.000 description 2
- 206010061216 Infarction Diseases 0.000 description 2
- 238000013184 cardiac magnetic resonance imaging Methods 0.000 description 2
- 208000029078 coronary artery disease Diseases 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000005562 fading Methods 0.000 description 2
- 239000012634 fragment Substances 0.000 description 2
- 230000007574 infarction Effects 0.000 description 2
- 208000028867 ischemia Diseases 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 208000010125 myocardial infarction Diseases 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000010200 validation analysis Methods 0.000 description 2
- 206010006322 Breath holding Diseases 0.000 description 1
- 101100402621 Homo sapiens MSANTD4 gene Proteins 0.000 description 1
- AYFVYJQAPQTCCC-GBXIJSLDSA-N L-threonine Chemical compound C[C@@H](O)[C@H](N)C(O)=O AYFVYJQAPQTCCC-GBXIJSLDSA-N 0.000 description 1
- 101150104791 MYOC gene Proteins 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 102100031642 Myb/SANT-like DNA-binding domain-containing protein 4 Human genes 0.000 description 1
- 238000002940 Newton-Raphson method Methods 0.000 description 1
- 102100026827 Protein associated with UVRAG as autophagy enhancer Human genes 0.000 description 1
- 101710102978 Protein associated with UVRAG as autophagy enhancer Proteins 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 210000001175 cerebrospinal fluid Anatomy 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 230000007012 clinical effect Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 229910003460 diamond Inorganic materials 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 230000003205 diastolic effect Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000010247 heart contraction Effects 0.000 description 1
- 208000019622 heart disease Diseases 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 238000004137 mechanical activation Methods 0.000 description 1
- 230000008450 motivation Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 238000000718 qrs complex Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000010410 reperfusion Effects 0.000 description 1
- 230000003362 replicative effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
- 210000002435 tendon Anatomy 0.000 description 1
- 230000008719 thickening Effects 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
- 230000009012 visual motion Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/45—For evaluating or diagnosing the musculoskeletal system or teeth
- A61B5/4519—Muscles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
- G01R33/56308—Characterization of motion or flow; Dynamic imaging
- G01R33/56333—Involving spatial modulation of the magnetization within an imaged region, e.g. spatial modulation of magnetization [SPAMM] tagging
-
- 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/92—Computer assisted medical diagnostics
- Y10S128/922—Computer assisted medical diagnostics including image analysis
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Signal Processing (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Psychiatry (AREA)
- Physiology (AREA)
- Vascular Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
Description
60/130,595号の利益を主張する。
財団授与番号第MIP93−50336号によって少なくとも一部は支持された
。
具体的には、心臓における物質点の追跡およびラグランジアンひずみを計算する
ために磁気共鳴画像化を使用して獲得される高調波位相画像(harmonic
phase image)を使用するプロセスに関する。
コンピュータ手段によって処理され、その後で結果が記録され、結果のハードコ
ピーを表示または生成する磁気エネルギーの応答放出を引き起こすために働くR
Fパルスおよび傾斜パルスを交互に生成する磁場内に、患者または対象が置かれ
る磁気共鳴画像化の有利な使用は、長い間公知であった。概して、Atalar
−McVeighの米国特許第5,512,825号およびConturo−R
obinsonの米国特許第5,281,914号を参照されたい。その両方は
、本発明の所有者に譲渡され、その開示は、参照により本明細書中に明確に援用
される。
の心臓壁の動きを画像化する場合に、画像面に直交する方向にある2組のタグ付
面を使用することは、公知である。米国特許第5,054,489号、ならびに
第5,111,820号、および5,217,016号を参照されたい。また、
Axel等のRadiology誌1989年171号841−845ページの
「磁化の空間変調を伴うMR画像化」、およびAxel等のRadiology
誌1989年172号(1)349−350ページの「心壁の動き:MR画像化
についての磁化の空間変調の改良方法」を参照されたい。
の補間を伴う変位およびひずみに関する詳細な動きのマップ内へ結合するために
、画像処理技術を使用することは、磁気共鳴のタグ付に関連して公知である。例
えば、Young等のRadiology誌1992年185号241−247
ページの「磁化の空間変調を伴う3次元の動きおよびゆがみ」、およびMcVe
igh等のRadiology誌1991年180号(3)677−683ペー
ジの「MR画像化での心筋のひずみにおける経壁傾斜の非侵襲性測定」を参照さ
れたい。これらのアプローチは、自動化されていない。なぜならいくつかの手動
による変調作業が必要であるから。
ことも公知である。概して、Prince等のIEEE Trans. on
Medical Imaging 1992年6月11(2)238−24
9ページの「タグ付MR画像シーケンスから動きの推定」、Gupta等のTe
chnical Report誌1995年JHU/ECE95−13の「タグ
付MRIのための可変輝度オプティカルフローについて」、およびGupta等
の1997年Santa Barbaraでの画像処理に関するIEEE国際会
議のプロシーディング「タグ付MR画像化のための帯域通過オプティカルフロー
」を参照されたい。そのようなアプローチにおいて、正弦波状タグパターンが、
飽和面(saturated plane)の代わりに使用される。画像輝度傾
斜は、画像のペアから推定される一時的な微分と共に、「オプティカルフロー」
と一般的に呼ばれる密な動きの推定を提供するために使用され得る特徴である。
そのようなアプローチは、輝度傾斜は傾斜の方向だけの動きについての情報を含
むという事実を補償するための調整を必要とする。
磁気共鳴画像化の使用を開示している。パルスおよび傾斜シーケンスは、空間的
に異なるタグ付パターンおよび差し引き像をタグ付画像を形成するために使用、
ペアで適用される。
ューム要素の位相角応答を観察すること、および外部刺激によって引き起こされ
る摂動前および中に生じる画像化を開示している。
トラストタグ付グリッドを使用することによって、対象の一部の量的な動きの推
定を開示している。
跡を開示している。
き分析と関連してこのようなアプローチを使用することもまた、公知である。そ
のような従来の技術は、通常、変位ベクトルまたはひずみパターンを介して、こ
れらの画像から動きを抽出することを含み、そして後に補間されるタグ特定およ
び位置推定を含む。
性を測定すること、および有限差分を使用して計算されたひずみを有する速度場
を再構築することによって、動きを直接測定するための方法を提供する。これら
2つのアプローチに関する課題の1つは、平面タグ付画像は、自動的には正確に
分析され得ないことである。位相コントラスト画像は、自動的に分析され得るが
、信号対雑音比が低くなりがちであり、許容し難い結果になる。
ることは公知である。タグ付心臓磁気共鳴画像化におけるこの10年間にわたる
発展で、インビボでの心臓の心筋層の詳細なひずみパターンを測定することが可
能となった。MRタグ付は、心筋層の「タグ」と呼ばれる一時的特徴を生成する
ために、被験者の長軸に沿う磁化を空間的に変調するための特別なパルスシーケ
ンスを使用する。タグ付MRIは、正常なおよび異常な心筋の動きのモデルを開
発しおよび洗練するために、心筋の動きの異常と冠状動脈病の相互関係をよりよ
く理解するために、心筋梗塞の治療の効果を理解するためペースメーカーを使用
する心臓の活動化パターンを分析するために、および虚血性心疾患を初期に探知
するためのストレステストと結びついて、使用されてきた。科学的取り組みが成
功しているにもかかわらず、タグ付MRIは、長い画像化ならびに後の処理時間
、画像化中の患者への不適切なアクセスおよび臨床医および同僚によるMR処理
の恩恵の理解の欠如のせいで、普通の臨床使用になることが遅かった。
鳴画像における孤立したスペクトルピークに対応する画像の位相を言う。この言
葉はまた、高調波位相画像またはHARP画像と本明細書中では呼ばれる。
ンスが、対象の空間変調された興味領域に適用され、そして空間変調された対象
のフーリエ領域からの少なくとも1つの第一のスペクトルピークが獲得される。
獲得された第一のスペクトルピークの逆フーリエ変換情報が計算され、そして第
一の高調波位相画像が、スペクトルピークの各々から計算される。このプロセス
は、各第二のスペクトルピークから第二の高調波位相画像を生成するために、異
なる時間に関して繰り返され、そしてひずみが、第一および第二の高調波画像か
ら決定される。
チにおいて、円周方向および半径方向のラグランジアンひずみ両方を決定するこ
とを含む。好適な実施形態において、SPAMMパルスシーケンスが使用される
。
、明白な動きの追跡が使用される。
とが、本発明の目的である。
らなる目的である。
てそれによって異なる時間に取られた点の画像のシーケンスを比較することによ
って心臓の部分内のひずみの決定を容易にすることが、本発明のさらなる目的で
ある。
して、心臓の解剖学的組織および心臓と共に動くタグ付特徴の両方の画像化を可
能にし、そしてそれによってラグランジアンひずみの決定を可能にするようなシ
ステムを提供することが、本発明のさらなる目的である。
化のための改良された方法を提供することが本発明の目的である。
ような方法を提供することが、本発明の別の目的である。
そして面ひずみまたはテンソルひずみの計算を提供するために使用されるような
システムを提供することが本発明のさらなる目的である。
することが、本発明のさらなる目的である。
成することが本発明のさらに別の目的である。
次の説明から、より十分に理解される。
きている一員を意味する。
動きについてモニタされている患者ならびに他のすべての項目およびそれらの部
分を含む。医療使用の中には、心壁、心筋および心腱の動きを測定する場合の使
用がある。モニタされる対象は、例えば血液または脳脊髄液などの流体、あるい
は固体または半固体、あるいはそれらの組合せであり得る。
明に関する背景的情報を提供し、そして本発明者らの米国特許出願シリアル番号
第09/131,589号の主題である。この従来の出願は、物理的に仮出願シ
リアル番号第60/130,595の一部であって、そしてそこから本出願へ移
行される。図1−7(b)に関連しての図8−18の記載中の「角度画像」の表
現は、「高調波位相画像」または「HARP画像」と表現する。
れる心臓の動きの欠陥を探知およびモニタする標準的アプローチに急速になって
いる。Zerhouni等のRadiology誌1988年169巻1号59
−63ページの「人間の心臓:MR画像化におけるタグ付心筋の動きの非侵襲性
アセスメントのための方法」、McVeigh等のRadiology誌199
1年180巻3号677−683ページの「MR画像化での心筋のひずみにおけ
る経壁傾斜の非侵襲性測定」、Moore等のJournal of Mag.
Res. Imaging1992年3月/4月号2巻165−175ページ
の「2平面タグ付MR画像からの3次元的左心室ひずみの計算」、を参照された
い。タグ付MRIは、所定のパターン中の心臓内および周囲の組織の磁気特性を
一時的に変化させるためにMRスキャナを使用し、組織は、次に心臓およびパタ
ーンの両方がゆがみを受ける場合に、画像化され得る。そうして生じる画像のシ
ーケンスにおけるゆがんだパターンの分析は、心壁内の心筋の動きについての情
報を生成する。画像分析は、通常特徴抽出方法を使用して行われるが、光フロー
法も検討されてきた。Fischer等のMag.Res.in Medici
ne1994年31巻401−413ページの「真の心筋の動きの追跡」、De
nney等のIEEE.Trans.Med.Imag.1995年14巻4号
625−635ページの「平面タグ付MR画像から3次元の左心室の動きの再構
築:推定理論的アプローチ」、Prince等のIEEE Trans.on
Medical Imaging 1992年6月11巻238−249ページ
の「タグ付MR画像シーケンスからの動き推定」、Amartur等のMag.
Res.Med.,1993年 29巻1号59−67ページの「心臓の動きを
研究するための新アプローチ:シネMR画像のオプティカルフロー」、およびG
uptaの1997年(Santa Barbaraにおける)画像処理に関す
るIEEE国際会議のプロシーディング3巻364−367ページの「タグ付M
R画像化のための帯域通過オプティカルフロー」もまた参照されたい。
なく、SPAMMタグ付画像から直接生成され得る。Osman等の1988年
(米国San Diegoにおける)医療画像化に関するSPIEの国際シンポ
ジウムのプロシーディングの「正弦波状MRタグ付を使用して心筋のひずみの2
次元成分の直接計算」、Axel等のRadiology誌1989年171巻
841−845ページの「磁化の空間変調を伴う動きのMR画像化」、およびA
xel等のRadiology誌1989年172巻1号349−350ページ
の「心壁の動き:MR画像化についての磁化の空間変調の改良方法」を参照され
たい。このアプローチは、SPAMMパターンについての信号モデルおよび動き
を基本的な搬送周波数の角度変調と解釈することに依存する。本発明は、非常に
小さな変位(エラー変位など)を推定し、タグ線を統合し、およびオプティカル
フローを計算する場合に直接有用であり得る。
ークの収束を有する。各スペクトルピークは、特定の方向の動きについての情報
を含む。帯域通過フィルタを使用することによって抽出され、これらのピークの
1つの逆フーリエ変換は、その位相が実際の動きの方向成分に線形に関連する複
雑な画像である。この位相画像は、前に本明細書で定義した「角度画像」である
。この位相画像は、(標準逆アークタンジェント演算子の働きによって)[−π
、π)の範囲にあるよう制限される。これは、角度ラッピングアーティファクト
(angle−wrapping artifact)である。角度ラッピン
グアーティファクトは存在するが、角度画像は、タグパターンを統合するために
使用され得、そして角度画像のペアは、小さな変位フィールド、調整を必要とし
ない画像ペア間のオプティカルフロー、ならびに2次元および3次元のひずみを
測定するために使用され得る。
末期において、繰り返しタグ付される。タグ付パルスシーケンスは、左心室が血
液で一杯であり、心臓が相対的にゆっくり動き、そしてECG信号のQRS群が
存在する心臓周期の一部である拡張末期で課される。推定の目的のために、拡張
末期は、t=0なる時刻とみなされ得、拡張末期における心臓内の点の位置は、
物質座標系として扱われ得る。連続した心臓シーケンス中、k空間は、フーリエ
空間における9の主要なスペクトルローブのうち少なくとも1つを含むフーリエ
変換情報を獲得するために、標準定常状態傾斜エコー画像化パルスシーケンスを
使用してスキャンされる。2つ以上のスペクトルピークが画像化される場合、帯
域通過フィルタは、1つのピークの情報だけを抽出するために適用される。この
データの逆フーリエ変換は、共役対称操作を行うことなく、行われる。共役対称
操作を行うことは、MR画像化においては通例である。その結果生じる複雑な画
像の角度は、角度画像を形成する。
エコー、らせん磁気またはエコー面磁気共鳴画像化などの代わりの公知の手段が
、例えば使用され得ることが理解される。
ansverse gradient pulse)が後に続くα度パルスを適
用することによって生成され得、これは、タグパターンの完成がクラッシャ傾斜
(crasher gradient)である別のα度パルスが後に続く画像面
内であり、これにより、コヒーレント水平磁化が不要になる。1D1−1SPA
MMパターンは、各々異なる周波数で生じ結果として1D1−1SPAMMタグ
付画像のフーリエ変換における3つのスペクトルピークが存在することになる3
つの複雑な画像の合計である。2次元の1−1SPAMMパターンは、2つの1
D1−1SPAMMシーケンスを急速に連続して適用することによって生成され
得る。この結果は、そのような画像中に9つのスペクトルピークに結果的になる
9つの複雑な画像の合計である2次元 1−1SPAMMタグ付画像となる。3
次元の1−1SPAMMパターンは、3つの1D1−1SPAMMパルスシーケ
ンスを急速に連続して適用することによって生成され得る。この結果は、そのよ
うな画像のフーリエ領域において27個のスペクトルピークの結果になる27個
の複雑な画像の合計である3次元の1−1SPAMMタグ付画像となる。一般的
に、複雑な画像およびタグ付画像中のスペクトルピークの数は、SPAMMパル
スシーケンスの数および特性に依存する。合成された、またはシミュレートされ
た2次元1−1SPAMMパターンは、リング型の対象に適用され、図1に示さ
れており、そしてそのフーリエ領域の大きさは9のスペクトルピークを示す図2
に示されている。
ンテキスト内で理解され得る。搬送高調波は、画像を空間的に振幅変調させ、そ
れゆえに、対応するスペクトルピークを搬送高調波の位置へと移動させる。
MR画像は、y1が出力方向で、且つy2が位相エンコーディング方向である画像
平面において、任意の点、y=[y1y2]Tにて、輝度値を与える、Ψ(y,t
)により表され得る。スペクトルピークの存在ゆえに、画像Ψ(y,t)は和と
して書かれ得る。
数kは、スペクトルピークを表すIDである。スペクトルピークの位置は、SP
AMMパルスシーケンスによって決定され得るベクトルwk=[w1k w2k w3 k ]Tによって決定される。スペクトルピークの合計数は2K+1である。この値
は、SPAMMパルスシーケンスの数および特性に依存する。
φk)を有する。故に、
使用して、Ψから取り出され得る。
3に示されている。この画像の鋸歯パターンは、範囲[−π,π)に包み隠され
る複合数量の角度の結果として生じる。
画像はタグ付けされたMRI画像における通常の2次元タグに類似した総合タグ
を生成するために使用され得る。この特徴の利点は、データが完全に自動的に生
成され、且ついかなる所望のタグ分類を用いても生成され得ることである。第2
に、画像は、オブジェクトの小さな変位を直接計算するために使用され得る。第
3に、角度画像は、2次元画像面における2次元のゆがみ、または3次元におけ
る完全な3次元ゆがみテンソルを直接計算するために使用され得る。これらのゆ
がみデータは、心筋虚血および心筋梗塞の検出およびそれらを量で表すときに有
用である。第4に、角度画像は、時系列の変位フィールドを表す標準的なオプテ
ィカルフローフィールドを形成するために使用され得る。
れぞれ示す。計算されたフィールドおよび真の変位フィールドの類似性は、角度
画像の動きが、動きを再構築するに効果的であることを示している。
れ得る。タグ付の線は、 (5) ak(y*,t)=a を満たす点{y*}の収束である。
て総合タグ付の線を生成することには、いくつかの利点がある。第1に、完全な
自動プロセスであることである。第2に、よい等輪郭(isocontour)
アルゴリズムがこの属性を有しているように、タグ付の線がサブピクセル解画像
度を有することである。第3に、全画像が、例えば、左右双方の心室心筋を含む
など、自動的に識別されるこれらのタグを有することである。最後に、範囲[−
π,π)内のN値を選択することによって、Nタグ付の線が、空間的期間2π/
wk以上の範囲において、統合されることである。原則的には、これらのタグ付
の線がどのぐらい接近して間隔を開けられ得るかを制限するものは基本的にない
。これは、これらのタグ付の線が互いにどのぐらい接近しているかのという特徴
に対する検出可能性によって制限されないからである。
よびwlの2つの角度画像(akおよびal)は、
使用され得る。ここで、Wは行列であり、Tは行列の転位であり、Δakは、
n)である。ここで関数x(y)は、 (9) x(y)=y1h1+y2h2+x0=Hy+x0 を使用して、画像平面座標システムの任意の点yを、磁気3次元座標システムの
位置x ε IR3にマッピングさせる。ここで、行列H ε IR3×2 =
[h2h2]であり、ベクトルh1 ε IR3および、ベクトルh2 ε IR3は
、それぞれ画像平面の読み出しおよび位相エンコーディング方向を表し、x0
ε IR3は、画像原点である。行列W ε IR3×2= [w1w2]である。
比較的小さな動きに対し実行される本発明の方法を容易にする。例えば、1つの
方法は、心臓が実質的に動く機会を得る前に、心弛緩終期後、即画像を作成する
ことである。このアプローチは、収縮期の最初の数10ミリ秒において、有用で
あり、且つ潜在的および臨床的に重要である。第2に、低周波タグパターンが使
用される場合、タグパターンの物理期間がだんだん大きくなり、動きは角度のあ
いまい性(ラッピング)を生じない。これを選択した際の潜在的困難性は、低周
波パターンのスペクトルピークが別のスペクトルピークと干渉し、所望されない
動きのアーティファクトが生じることである。第3のアプローチは、心弛緩終期
およびそのすぐ後の画像から、固定の補正値でタグパターンを適用することであ
る。この場合、タグパターンを急速に動く組織へ適用することが必要である。最
後に、この問題を回避する方法は、タグ付け時間からの変位ではなく、画像時間
の間の変位を画像作成することである。これらの変位は小さい。このアプローチ
は、古くは「オプティカルフロー(optical flow)」と呼ばれる。
立したベクトルwkおよびwlの、少なくとも2つの角度画像である(akおよび
al)から直接計算され得る。e方向での2次元のゆがみは、
の派生を計算すると同時に、ラッピングの不連続性を克服するために使用される
。これらの方程式から計算されたゆがみは、オイラー形式(Eulerian
sense)である。完全なゆがみテンソルは、3つのスペクトルピークから来
る3つの角度画像から計算され得る。3つのスペクトルピークの生成は、3次元 SPAMMパターンを使用し、画像平面ではなく画像体積を獲得することによ
って行われる。
として定義される。概略は、E.C.Hildrethの「Measureme
nt of Visual Motion」、MIT Press、Cambr
idge、1984年を参照されたい。現在のコンテキストにおいては、「明ら
かな」という語は、真の3次元の動きではなく、画像平面の動きを示唆する。従
来技術のコンテキストにおいて、オプティカルフローの定義は、速度フィールド
を含み、概してこの速度フィールドの緻密な推定を得るために、ある種の規則化
が必要である。概略は、Hornらの「Determining Optica
l Flow」、Artificial Intelligence、17:1
85−203、1981年を参照されたい。本発明のコンテキスト内における角
度画像の使用により、規則化を使用することなしに、速度フィールドを直接計算
することが可能である。
クトルwkおよびwlを有した、少なくとも4つの角度画像ak(y,t)、al(
y,t)、ak(y,t+Δt)、およびal(y,t+Δt)が採用され得る。
2つの画像間の時間間隔、Δt>0は、時間間隔における一定の動きの速度v2
を仮定するのに十分小さいことが好ましい。2次元の動きの速度は
室を示す断面図である。図7bは、本発明の方法により作成される図7aの左心
室の断面図の角度画像を示す。
た心臓の磁気共鳴画像シーケンスの迅速な分析のための画像処理方法を提供する
。この方法は、心臓の動きに関する情報を含むSPAMMタグの磁気共鳴画像に
おける分離されたスペクトルピークの使用を含む。SPAMMスペクトルピーク
の逆のフーリエ変換は、複合画像であり、その画像の計算された角度は高調波位
相(HARP)画像と呼ばれる。この方法は、2つのHARP画像シーケンスを
使用し、自動的且つ正確に時間的に物質点を追跡する。1実施形態において、迅
速な、半オートメーション化されたプロシージャは、これらの痕跡を使用し、周
囲ゆがみおよび放射状ゆがみ双方を含むラグランジュのゆがみを計算する。これ
らの方法は、2次元のコンテキスト内で開発されたが、3次元で使用され得る。
この新しい計算アプローチにより、スキャン完了後約5〜10分内の迅速な分析
および可視化が可能となる。その性能は、MR画像シーケンスに示され得、正常
な心臓の動きおよび異常な心臓の動きの双方を反映する。
houniらの「Human heart:Tagging with MR
imaging−−a method for noninvasive as
sessment of myocardial motion」、Radio
logy、169(1):59−−63、1988年、Axelらの「MR i
maging of motion with spatial modula
tion of magnetization」、Radiology、171
:841−−845、1989年、McVeighらの「Cardiac ta
gging with breath−hold cine MRI」、Mag
n.Res.Med.、28:318−−327、1992年、Fischer
らの「Improved myocardial tagging contr
ast」、Magn.Res.Med.、30:191−−200、1993年
、Atalarらの「Minimization of dead−perio
ds in MRI pulse sequences for imagin
g oblique planes」、Magn.Res.Med.、32(6
):773−−777、1994年12月、およびFischerらの「Tru
e myocardial motion tracking」、Magn.R
es.Med.、31:401−−413、1994年)により、ビボ心臓にお
ける心筋の詳細なゆがみパターンを測ることが可能になった。(Youngらの
「Three−dimensional motion and deform
ation of the heart wall:Estimation w
ith spatial modulation of magnetizat
ion −− a model−based approach」、Radio
logy、185:241−−247、1992年、Mooreらの「Calc
ulation of three−dimensional left ve
ntricular strains from biplanar tagg
ed MR images」、J.Magn.Res.Imag.、2(2):
165−−175、1992年3月/4月、Parkらの「Analysis
of left ventricular wall motion base
d on volumetric deformable models an
d MRI−SPAMM」、Med.Image Anal.、1(1):53
−−71、1996年、Denney,Jr.らの「Model−free r
econstruction of three−dimensional m
yocardial strain from planar tagged
MR images」、J.Magn.Res.Imag.、7:799−−8
10、1997年、およびE.R.McVeighの「Regional my
ocardial function」、Cardiology Clinic
s、16(2):189−−206、1998年)MRタグ付けは、特別なパル
スシーケンスを使用し、被験体の縦方向の磁化を空間的に変調し、心筋内にタグ
と呼ばれる一時的な機能を作り出す。すぐだめになる勾配エコー画像作成技術は
、CINEシーケンスを作成するために使用され、心臓の解剖、および心臓と共
に動くタグ機能の双方の動きを示す。異なる方向且つ異なる時間から撮られた多
くの画像におけるタグ機能の動き分析は、3次元における物質点を追跡するため
に使用され得、心筋内のゆがみパターンの詳細なマップを導く。概略は、E.R
.McVeighの「Regional myocardial functi
on」、Cardiology Clinics,16(2):189−−20
6、1998年およびE.R.McVeighの「MRI of myocar
dial function:motion tracking techni
ques」、Mag.Res.Imag.、14(2):137、1996年を
参照されたい。
しく目立ってきた。タグ付けされたMRIは、正常および異常な心筋の動きのモ
デルを発展および精錬するために(Mooreらの「Calculation
of three−dimensional left ventricula
r strains from biplanar tagged MR im
ages」、J.Magn.Res.Imag.、2(2):165−−175
、1992年3月/4月、E.R.McVeighの「MRI of myoc
ardial function:motion tracking tech
niques」、Mag.Res.Imag.、14(2):137、1996
年、Clarkらの「Circumferential myocardial
shortening in the normal human left
ventricle」、Circ.、84:67−−74、1991年、Mc
Veighらの「Noninvasive measurements of
transmural gradients in myocardial s
train with MR imaging」、Radiology、180
(3):677−−683、1991年、およびLugo−Olivieriら
の「The effects of ischemia on the tem
poral evolution of radial myocardial
deformation in humans」、Radiology、19
3:161、1994年)、異常な心筋の動きがある、冠状動脈疾患の相関関係
をよりよく理解するために(McVeighらの「Imaging async
hronous mechanical activation of the
paced heart with tagged MRI」、Magn.R
es.Med.、39:507−−513、1998年)、ペースメーカーを使
用して心臓の活性化パターンを分析するために(Limaらの「Segment
al motion and deformation of transmu
rally infarcted myocardium in acute
postinfarct period」、Am.J.Physiol.、26
8(3):H1304−−12、1995年)、心筋梗塞後の処置の効果を理解
するために(Croisilleらの「Combined dobutamin
e stress 3D tagged and contrast enha
nced MRI differentiate viable from n
on−viable myocardium after acute inf
arction and reperfusion」、Circ.、92(8)
:I−−508、1995年)、心筋虚血の早期発見のためのストレステスト(
Budingerらの「Cardiac MR imaging:Report
of a working group sponsored by the
national heart,lung,and blood insti
tute」、Radiology、208(3):573−−576、1998
年)心臓のMRIのための動機)と共に使用されてきた。これらの成功した使用
例にもかかわらず、タグ付けされたMRIは日常的臨床使用に参入するのにかな
り時間がかかっている。これは一部には、長い画像作成および後処理時間、画像
作成の間の患者への不十分なアクセス、MR処理および臨床家およびその提携者
による利益に関する理解の欠如のためである(Youngらの「Trackin
g and finite element analysis of str
ipe deformation in magnetic resonanc
e tagging」、IEEE Trans.Med.Imag.、14(3
):413−−421、1995年9月)。
れ得る。その3つの段階は、(1)2次元画像で左心室(LV)心筋を見つけ、
(2)LV壁内のタグ機能の位置を推定し、(3)これらの測定からゆがみフィ
ールドを推定することである。準自動化されたアプローチが同様に提案されてき
てはいるが(T.S.Denneyの「Identification of
myocardial tags in tagged MR images
without prior knowledge of myocardia
l contours」、In J.DuncanおよびG.Gindi、編集
、Proc.Inf.Proc.Med.Imag.、ページ327−−340
、1997年)、多くの既知のアプローチが、心内膜および心外膜の完全な手作
業の輪郭取りに頼っている(Aminiらの「Coupled B−snake
grids and constrained thin−plate sp
lines for analysis of 2次元 tissue def
ormations from tagged MRI」、IEEE Tran
s.Med.Imag.、17(3):344−356、1998年6月、およ
びGuttmanらの「Tag and contour detection
in tagged MR images of the left ven
tricle」、IEEE Trans.Med.Imag.、13(1):7
4−−88、1994年)。最近の研究はまた、完全に自動化された輪郭取りも
示唆してきている(Kerwinらの「Tracking MR tag su
rfaces using a spatiotemporal filter
and interpolator」、Int.J.Imag.Sys.Te
ch.、10(2):128−−142、1999年)。ほとんどの場合、輪郭
取りの結果は、タグ機能の推定段階において必要であり、この推定段階には、い
くつかの利用可能な準自動化された方法(Aminiらの「Coupled B
−snake grids and constrained thin−pl
ate splines for analysis of 2次元 tiss
ue deformations from tagged MRI」、IEE
E Trans.Med.Imag.、17(3):344−−356、199
8年6月、およびT.S.Denneyの「Identification o
f myocardial tags in tagged MR image
s without prior knowledge of myocard
ial contours」、In J.DuncanおよびG.Gindi、
編集、Proc.Inf.Proc.Med.Imag.、ページ327−−3
40、1997年)、ならびに完全な自動化にとって恐らく非常に有望な、新し
いアルゴリズム(Kerwinらの「Tracking MR tag sur
faces using a spatiotemporal filter
and interpolator」、Int.J.Imag.Sys.Tec
h.、10(2):128−−142、1999年、およびMoultonらの
「Spline surface interpolation for ca
lculating 3D ventricular strains fro
m MRI tissue tagging」、Am.J.Physiol.(
Heart Circ.Physiol.)、270:H281−−H297、
1996年)がある。
および分化計算であり、有限要素法を含む、文献に記載のいくつかの方法(Mc
Veighらの「Noninvasive measurements of
transmural gradients in myocardial s
train with MR imaging」、Radiology、180
(3):677−−683、1991年、Aminiらの「Coupled B
−snake grids and constrained thin−pl
ate splines for analysis of 2次元 tiss
ue deformations from tagged MRI」、IEE
E Trans.Med.Imag.、17(3):344−356、1998
年6月、およびO’Dellらの「Three−dimensional my
ocardial deformations:Calculations w
ith displacement field fitting of ta
gged MR images」、Radiology、195:829−−8
35、1995年)、グローバルポリノミアルフィッティングアプローチ(gl
obal polynomial fitting approach)(De
nneyらの「Reconstruction of 3D left ven
tricular motion from planar tagged c
ardiac MR images:An estimation theor
etic approach」、IEEE.Trans.Med.Imag.、
14(4):625−−635、1995年)、ならびにいわゆるモデル無し確
率論的推定アプローチ(model−free stochastic est
imation approach)(Denney,Jrらの「Model−
free reconstruction of three−dimensi
onal myocardial strain from planar t
agged MR images」、J.Magn.Res.Imag.、7:
799−−810、1997年、およびKerwinらの「Cardiac m
aterial markers from tagged MR image
s」、Med.Image Anal.、2(4):339−−353、199
8年)がある。タグ表面を推定する方法(O’Dellらの「Three−di
mensional myocardial deformations:Ca
lculations with displacement field f
itting of tagged MR images」、Radiolog
y、195:829−−835、1995年、Aminiらの「Flexibl
e shapes for segmentation and tracki
ng of cardiovascular data」、In Proc.I
EEE Int.Conf.Image Proc.、ページ5−−9、IEE
E Comp.Soc.Press、1998年、およびOsmanらの「Di
rect calculation of 2次元 components o
f myocardial strain using sinusoidal
MR tagging」、In Proc.SPIE Med.Imag.C
onf.、1998年2月」San Diego)は、タグ識別およびゆがみ推
定の間にある中間段階を表す。これらのタグ付けされたMRI処理方法の間にあ
る明らかな差にもかかわらず、それらはすべて2つの要となる制限を共有する。
すなわち、タグつけされたMRI処理方法は、完全に自動化されておらず、且つ
緻密なゆがみ推定を形成するために、補間を必要とする。本発明はこれらの関心
の双方を取り扱う。
画像の分析への新しいアプローチを記載してきた(Osmanらの「Motio
n estimation from tagged MR images u
sing angle images」、In Proc.Int.Conf.
Imag.Proc.、Volume 1、ページ704−−708.Comp
.Soc.Press、1998年、Chicago、Osmanらの「Ima
ging heart motion using harmonic pha
se MRI」、1998年10月」提出、およびShinnarらの「Inv
ersion of the Bloch equation」、J.Chem
.Phys.、98(8):6121−−6128、1993年4月)。このア
プローチは、SPAMMタグパターンに基づいており(Axelらの「MR i
maging of motion with spatial modula
tion of magnetization」、Radiology、171
:841−−845、1989年)、SPAMMタグパターンは、基礎を成す画
像を振幅変調し、フーリエドメインにずらりと並んだスペクトルピークを生成す
る。これらのスペクトルピークのそれぞれは、組織の動きの特定の構成要素に関
する情報を伝える。この情報は、位相変調方法を使用して抽出され得る。Shi
nnarらの「Inversion of the Bloch equati
on」、J.Chem.Phys.、98(8):6121−−6128、19
93年4月には、統合タグ付の線を再構築し、小さな変位フィールドを計算し、
およびオイラーゆがみ画像を計算する、シングルショットHARP画像分析技術
と呼ばれ得るものが記載されている。これらの方法は、心周期内でのシングル位
相(時間フレーム)のみからのデータを必要とするが、限られている。なぜなら
ば、これらの方法は、動きの物質的性質を計算し得ないからである。本発明にお
いて、この方法は画像シーケンス−CINEタグ付けされたMR画像を採用し得
、画像シーケンスは、好適には、物質点追跡技術およびこれらの追跡された点を
使用して、周辺ゆがみおよび放射状ゆがみを含むラグランジュゆがみを計算する
方法の双方を含む。
収集され得るデータを使用する。これらの方法は、2次元画像に直接のみ適用す
る。この結果、推定された動きの量は、3次元の動きから2次元平面への投影を
表し、推定された動きの量は「明らかな」動きとして考えられるべきである。本
方法は、3次元画像へと拡張され得るが、緻密なデータ獲得方法が採用される。
本方法は、即座の臨床効果を有すべきである。これは、その方法が、自動であり
、そして周辺ゆがみが、左心室の動きの分析において特に重要であるからである
。
=[y1 y2]T、且つ時間tの際のタグ付けされた心臓MR画像の強度を示す
。犬の心臓の異常な動きを示す典型的な画像が、図8(a)に示されている。左
心室(LV)が画像中心で輪のように見えている。タグ付けの効果は、タグパタ
ーンによって、基礎を成す画像の増加として記載され得る。図8aにあるパター
ンは、1次元SPAMMタグパターン(グリッド)であり(Axelらの「MR
imaging of motion with spatial modu
lation of magnetization」、Radiology、1
71:841−−845、1989年)、これはある基本振動数を有する有限の
余弦級数として書かれ得る(M.E.Gurtin、「An Introduc
tion to Continuum Mechanics」、Academi
c Press、Inc.、1981年)。このパターンを増加することによっ
て、基礎を成す画像の振幅変調が生じ、フーリエ変換を図8bに示されるような
パターンへと複製する。フーリエ空間におけるスペクトルピークの位置は、SP
AMMタグパルスシーケンスによって決定される基本タグ振動数の整数倍である
。
用して生成され得る。本発明の方法はこの場合においても利用され得る。
rt motion using harmonic phase MRI,O
ctober 1998.」, Shinnar et al.,「Inver
sion of the Bloch equation,」J.Chem.P
hys.,98(8):6121−−6128,April 1993)は、帯
域通過フィルタを使用し、一定のタグ方向において通常最低次の高調波周波数で
ある周波数ωkに集中しているk番目のスペクトルピークを切りはなす(iso
late)。帯域通過フィルタは通常、リンギング(ringing)を減少さ
せるために滑らかにロールオフする端(edges)付の楕円形サポート(el
liptical support)を有する。図8(b)に描かれた輪郭は、
この場合は円であるが、このデータを処理するために使用される帯域通過フィル
タの−3dBの等輪郭を示す。一旦フィルタが選択されると、回転されたバージ
ョンを使用して垂直タグ付の画像を処理することを除いて、全ての画像のシーケ
ンスにおいて同じフィルタが使用される。Shinnar et al.,「I
nversion of the Bloch equation,」J.Ch
em.Phys.,98(8):6121−−6128,April 1993
において、最適な実施のためのフィルタの選択が、議論される。帯域通過領域の
逆フーリエ変換により、
ド画像(the magnitude image)と呼ばれ、φkは位相画像
と呼ばれる。式15におけるIkの使用は、式10のΨkと同じである。
ジングによる画像強度の変化の両方を反映する。図8(c)は、図8(b)のフ
ィルタを使用して図8(a)から抽出された高調波マグニチュード画像を示す。
フィルタリングプロセスによる不鮮明を除いて、その基本画像のように基本的に
は見える。高調波マグニチュード画像にタグパターンがないために、高調波マグ
ニチュードオ画像を使用して、組織を背景と区別するための断片を提供し得る。
単純な閾値を使用して、未加工の断片を提供し得る。ここで、その閾値は、拡張
末期と収縮末期の両方の時点で手動で選択され、これらの時間の間で線形補間さ
れる。
的に、Ikの位相は、実部で虚部を割ったアークタンジェントを使用することに
よって計算され得る。Ikの符号を説明すると、この計算の固有の範囲は、例え
ば、C、Fortran、またはMATLABの演算子atan2を使用して、
[−π,−π)に拡張され得る。しかし、これによると、一般的に実際の直線全
体の値を示す実際の位相でなく主値のみが生成される。この主値は、ak(y,
t)によって示され得る;それは、
現は、φkでなくakに一般に使用され、データからそれを直接計算し、視覚化し
得る。その2つは混同されるかもしれないが、φkを高調波位相と呼び、akを高
調波位相角度と呼ぶ。図8(b)でアウトラインが示されるスペクトルピークに
対応するHARP角度画像が図8(d)に示される。明瞭さのために、それは、
図8(c)のマグニチュ−ド画像の生の断片を使用して作れられたマスクで表示
される。
ァクトによる鋭い転位によって垂直方向に妨げられた強度傾斜がわかる。これら
の転位の位置は、図8(a)のタグ付の線とほとんどかなり一致し、両方は、収
縮期の間に生じる心筋の動きを反映している。HARP画像の強度傾斜は、実際
、元の画像において十分明白な情報以上の濃密な動きの情報を含む。例えば、H
ARP角度画像の計算された等輪郭は、任意に分離した心筋中のタグ付の線を生
成し得る(Osman et al.,「Imaging heart mot
ion using harmonic phase MRI」,Octobe
r 1998 および Shinnar et al,「Inversion
of the Bloch equation.」J.Chem.Phys., 98(8):6121−−6128,April 1993)。基本的な原理
は、高調波位相およびHARP角度の両方は、タグ付の組織という材料の性質と
いうことである;それゆえ、質点は、その動きを通じてそのHARP角度を保持
する。これは、動きについてHARPトラッキングの基本である。
グ付の勾配は、一般に画像平面において適用される。この場合、拡張末期の画像
に生じるタグ付の線は、図9(a)に示されるように、実際に画像平面に直交す
るタグ面の一部である。高調波位相は材料の特性であるので、拡張末期に同じ高
調波位相φを有する一組の点はまた、画像平面に直交する平面であり、単に別の
タイプのタグ面であると考えられ得る。拡張末期で同じHARP角度aを有する
一組の点は、ただ1つの平面ではなく平行平面の集合を含む。これは、次のセク
ションで説明するHARPトラッキングにおいて問題を引き起こす。見かけ上の
動きを説明するために、高調波位相φによって与えられる固有の関係が考えられ
る。
2つの高調波位相を有する一組の点は、単一の点で交差する画像平面に直交する
線を含むことが図9(a)から明らかである。図9(b)に示されるように、タ
グ面が動き下でひずみ、それによってこの線を曲線に歪ませる。動きについての
緩い仮定の下で、この曲線はまだ、1つの点で画像を交差する。次いで、この点
は、拡張末期で対応する点と関連し、画像平面内の見かけ上の動きを表示し得る 見かけ上の動きが、q(y,t)によって示される見かけ上の参照マップを使
用して数学的に説明され得る。この関数は、点を(同じ2つの高調波位相を有す
るという意味において)yが見かけ上の拡張末期にあった画像平面内に与える。
q(y,t)が、拡張末期の真の3次元の質点の位置を画像平面に直交に投影し
たものであることが図9から明らかに示され得る。見かけ上の2次元の動きを計
算するには制限があるが、その計算は、真の3次元の動きと非常によく似た関係
を有する。ひずみ等の見かけ上の動きから導出された動きの量は、等しく厳密な
方法による真の3次元量に関係し得る。
かけ上の動きのトラッキングと、(2)そのようなトラッキングされた点からの
ラグランジアンひずみの計算を含む。コンパクトな式にするために、ベクトル表
示を使用する。詳細には、ベクトルφ=[φ1 φ2]T、およびa=[a1 a2
]Tが定義され、それぞれ高調波画像I1およびI2についての一対の高調波位相
画像およびHARP角度画像を説明する。
値であるHARP角度が、計算で使用される。多くの点が同じ一対のHARP角
度を有する画像平面にあることが即座に結論される。2つのHARP角度を有す
る所与の質点に対して、同じHARP角度を共有する点の内の1つのみが正しく
マッチングしたものである。すなわち、それはまた、同じ一対の高調波位相を共
有する。見かけ上の動きが1つの画像から次の画像まで小さい場合、最近接のこ
れらの点が正しい点であることは明らかである。本発明の方法は、タグ付のMR
画像のCINEシーケンスを介して見かけ上の動きをトラッキングする。位置が
ymで、時間がtmの質点を仮定する。ym+1が、時間tm+1でのこの点の見かけ上
の位置である場合、
準を提供する。
元且つ非線形で、ルートファインディング問題(root finding p
roblem)である。この問題は、ニュートンラプソン法を使用して反復的に
解かれ得る。単純化の後、ニュートンラプソン反復式(iteration)は
、
φが利用できないことであり、その代わりにaを使用しなければならない。幸い
、φを含む式をaを含む式で置き換えるのは比較的簡単である。akの勾配は、
勾配の大きさが理論的に無限−−実際には非常に大きい(図8(d)参照)、ラ
ッピングアーティファクト以外はΦkの勾配と同じであることは式(16)から
明らかである。πをakに加え、再度ラッピングすることによって、1/2の空
間距離だけラッピングアーティファクトをシフトさせ、それによってこの勾配が
結果的にΦkの勾配に等しくなることは、式(16)から明らかである。結果と
して、Φkの勾配は、akおよびW(ak+π)の勾配の(大きさにおいて)小さ
いものに等しい。形式的に、これは、
tm)の計算であり、この計算は、それのラッピングされたバージョンである高
調波位相だけでなく位相自体も分からないのでほとんど不可能と思われる。しか
し、k=1,2、|φk(y(n),tm+1)−φk(ym,tm)|<π(微小な動き
を仮定)に関して、
と、フォーマルな解である
ある。第1に、位相ラッピングのために、その解はもはや固有ではない;実際、
両方の方向におけるほぼタグ期間(tag peirod)毎にa(y,tm+1
)=a(ym,tm)を満たす解が期待され得る。それゆえ、「よい」初期の点で
開始することおよび間違った解に向かうのを防ぐためのステップサイズを保持す
ることの両方が望ましい。従って、y(0)=ymでアルゴリズムを初期化すること
且つ1ステップを1画素の距離に制限することが望ましい。第2の問題は、任意
のyに対してa(y,tm+1)の推定を行わなければならない。直線の双一次補
間は通常、働く;しかしこの場合、ラッピングアーティファクトが誤った結果を
引き起こす。これらの誤りを防ぐために、yの近傍のaのローカル位相非ラッピ
ング(local phase unwapping)を実行し、非ラッピング
角度を双一次補間し、次いでその結果をラッピングし、補完されたHARP角度
を作成する。最後に考慮すべきことは、判定基準を止めることである。以下の2
つの判定基準が使用される:(1)計算された角度が望ましい目標HARPベク
トルに十分近いかまたは(2)反復カウントが超えられるか。
ームのその見かけ上の点までトラッキングするアルゴリズムが即座に定義され得
る。画像シーケンス全体を介して点をトラッキングするHARPトラッキングア
ルゴリズムを定義するのを容易にするために、より一般的なフレームワークでこ
のアルゴリズムを変換することは有用である。従って、その方法では、yinitは
、検索が開始される初期化であると考え、a*が目標HARPベクトルであると
考える。時間tm+1におけるその見かけ上の位置まで時間tmの点ymをトラッキ
ングするために、yinit=ym,a*=a(ym,tm)と設定し、最大反復カウン
トNを求め、次いで以下のアルゴリズムを実行する: (アルゴリズム1(HARP目標設定)n=0とし、y(0)=yinitと設定す
る) 1.
ーケンスの各画像にうまく適用される。正しく見かけ上の動きを見出す好適なア
プローチは、全体のシーケンスを介して目標HARPベクトルを同じに、すなわ
ち元の点のHARPベクトルと等しいように維持することであり、シーケンスに
おいて以前に推定された見かけ上の位置にHARP目標設定を初期化することで
ある。これにより、同じHARP角度を有する連続した点を生成し、一般に所望
の解に近いHARP目標設定で使用された初期の点を維持することによって間違
った解に向かうのを避ける。形式的にこのアルゴリズムを説明するために、時間
tj+1,tm+2・・・・・の全ての画像を通じて時間tjのyjをトラッキングする
ことを所望すると仮定する。
時間の逆方向で点をトラッキングし得る。それゆえ、任意の時間の任意の画像の
点を特定し、順方向と逆方向の両方に点をトラッキングし、それによって空間と
時間において任意の点の完全なトラジェクトリを与えることが可能になる。
定する。HARPトラッキングは、(利用可能な画像倍の集合に対して)時間t
での質点の位置yに対するt=0での質点qをトラッキングし(このセクション
において、全ての位置の量を「見かけ上の」位置と呼ぶ)得る。これにより動き
マップy(q,t)の推定値を提供する。y(q,t)を使用して、有限の違い
を使用して任意の質点qおよび時間tの変形勾配テンソル
計算される(Atalar et al.,「Optimization of
tag thickness for measuring positio
n with magnetic resonace imaging,」IE
EE Trans.Med.Imag.,13(I):152−−160,19
94)。HARPトラッキングのさらに強力な適用は、準自動分析内のいくらか
単純な以下の計算で示される。
距離が縮むと負になり、点間の距離が伸びると正になる。左心室壁内の任意の位
置の円形ひずみを測定するために、LV長軸の中心である円に沿って2つの点を
単純に配置する。動径方向ひずみを測定するために、長軸から生じる光線に沿っ
て2つの点を単純に配置する。どちらかの場合に、ひずみは、HARPトラッキ
ングを使用して、2つの点をトラッキングし、eを計算することによって測定さ
れる。HARPトラッキングが基本的に画像内の任意の位置をトラッキングでき
るので、これらの点の位置が「タグ付の線の交差」または画素の位置にさえある
必要はないことを再度強調しなければならない。ひずみのこの測定は、ラグラン
ジアン(またはユークリッドひずみ)テンソルの緻密な計算に対して2つの利点
を有する。第1に、画像全体(または対象の領域)の代わりに2つの点のみのト
ラッキングが必要があるので、それは非常に速い。第2に、点は、1つの画素に
比べ一般にかなり離れて配置されるので、伸びについての計算は、元来ノイズに
敏感ではない。
域性の動径方向および円形ラグランジアンひずみ(circumferenti
al Lagrangian strain)を計算するアプローチが実行され
た。単純なユーザインターフェースは、図10に示すようにLV壁内の3つの同
心円の配置を可能にする。これらの円は、第1に心室の中心(長軸の位置)上で
クリックし、次いで心外膜に対する1つの円および心内膜に対する別の円をドラ
ッグすることによって手動で配置される。3番目の円は、これら2つの円の間の
中間に自動的に配置される。通常、円は拡張末期の画像で定義されるが、LVの
断面がより円形であり得るので、収縮末期の画像を使用してそれらを定義するの
がより容易である。16個の等間隔の点が各円の外周の周りに自動的に定義され
、アルゴリズム2(HARPトラッキング)を使用して48点全てが(時間の順
方向または逆方向に)トラッキングされる。
ことによって計算される。円は、拡張末期に、収縮末期に、またはその間の任意
の時点に定義されるかどうかに関わらず、拡張末期の画像フレームは、物質標準
として使用される。同じ円にある点間の距離の変化が円形ひずみに対応する。動
径方向の点間の距離の変化は動径方向ひずみに対応する。3つの同心円状がある
ので、心内膜、心外膜、および中間壁の円形周囲ひずみおよび心内膜および心外
膜の動径方向ひずみが計算され得る。ノイズを幾分か低減し、それによって提示
を単純化するために、円は、8個の八分円に分けられ、計算されたひずみは、こ
れらの各八分円で平均化された。慣習により、八分円は、図10に示されるよう
に隔壁の中心から時計方向に番号付けされた。その結果、生じるひずみは、時間
と八分円の両方の関数としてプロットされ、それに断面図内の心臓の動きを空間
時間的に表示する。
アルゴリズム2を独立して適用することによってうまく達成され得る。しかし、
場合によっては心筋の大きな動きまたは画像アーティファクトにより、ある時間
フレームで点が誤った目標に収束する(タグジャンプ)ことが起こり得、同様に
、連続したフレームで誤ったトラッキングが生じることになる。本発明は、1つ
以上の正確にトラッキングされた点を使用する改良手順を使用し、それによって
誤ってトラッキングされた点のトラッキングを訂正する。
た位置で初期化される。しかし、2つの時間フレームの間の平面内の動きが非常
に大きい場合、この初期点は正確な解から非常に離れ、誤った点に収束し得る。
改良は、アルゴリズム2に対するよりよい初期化についての体系だった識別に基
づく。所与の円にある1つの点は、全てのフレームを通して正確にトラッキング
されてきたと証明され得ると仮定する。発明者らの実験において、動きが相対的
に小さい隔壁のそのような点が常に見出されてきた。
全ての点の初期化を改善し得、全体にわたって集約されたトラッキング結果を改
良し得る「アンカー」として使用される。
とアンカーをつなぐ線セグメントにある1画素未満に分けられた点のシーケンス
が定義される。アンカーは、いくつかの特定の時間に点yまでトラッキングされ
ると仮定する。この点が正確にトラッキングされると仮定することによって、そ
のアンカーに近い点の正確なトラッキングの結果がyに近くなる。従って、yは
、アルゴリズム2の初期点として使用され、シーケンスの最初の点をトラッキン
グする。次いで、この結果は、シーケンス上の第2の点等に対する初期点として
使用される。円上のアンカーの第1の隣接する点に達することによって、誤った
結果に収束し、タグをジャンプさせる機会はない。隣接する点は、新しいアンカ
ーとして役に立ち、その手順は、円上の全ての点がトラッキングされるまで円上
の次に隣接する点に対して繰り返される。
補正するために使用され得る。一般に、タグジャンプ誤りは、自由壁でのみ起こ
り、改良が計算的に要求されるにしたがって、その演算は一般に一定の時間では
1つの円に抑制される。所望ならば、チェックして、円は、元のアンカーの周囲
をトラッキングすることによって、完全なものになり得る。その結果が異なる場
合、総計としての誤りがあり、円を再定義することが望ましい。タグフェージン
グまたは他の画像アーティファクトによりこのタイプの総計としての誤りが生じ
得るが、平面外の動きによってこのタイプの総計としての誤りがさらに生じ易く
なる。平面外の動きにより、画像化される実際の組織は変化する。この場合、ト
ラッキングされた点に対応する高調波位相角度をもたらす画像平面上に組織が存
在することはあり得ない。要約すると、タグが消え、HARPトラッキングアル
ゴリズムは、正確なHARP角度を有する別の点に簡単に収束する。左心室の動
きとジオメトリとの間の特定の関係のために、これは大した問題ではない。問題
は、左心室の境界近くでより生じ易いので、それにより課される主な制限は、心
外膜または心内膜の非常に近く円を配置してはならないことである。
いての実験が、犬と人間の両方を含む試験について行われた。異常に拍動する犬
のデータセットが、ドブタミンの応力下にある正常な人間の心臓からのデータと
同様に得られた。HARPトラッキングと従来技術の周知のタグ付トラッキング
法(tag tracking method)の間の質的な比較が為された。
いてのタグ付画像のセットが使用された。これらのデータは、タグ付磁気共鳴画
像法および分析技術を使用した、電気的に拍動させる活性化下での心臓の動きの
研究において、以前に使用された。実験手順およびその結果についての十分な説
明が、Lima等の「Segmental motion and defor
mation of transmurally infarcted myo
cardium in acute postinfarct period.
」Am.J.Physiol.,268(3):H1304−−12,1995
において与えらている。この結果は、完全な3次元において見かけ上の動きおよ
びひずみを説明するのでなく、1つの断面積上の見かけ上の動きおよびひずみの
みを説明する。しかし、本発明で得られた結果はLima等による方法の結果と
ほとんど同じであるが、Lima等による方法にかかる時間のほんの一部(a
fraction of the time)で生成される。これは、本発明の
方法を使用することにより結果が正確かつ短時間で生成されるという利点を示し
てる。
basal free wall)内に配置された。ソフトウェアリリース4.
7(General Electric Medical Systems,M
ilwaukee,WI)を備えた標準的な1.5Tスキャナで、磁気共鳴画像
法が実行された。6msのSPAMMパルスシーケンスが心筋についてのタグ付
パターン(tag pattern)を生成するように使用され、そのタグ付パ
ターンは画像平面の中で5.5mmで分けられた平行面飽和帯(paralle
l plane saturation bands)を含む。パルスシーケン
スにタグ付けすることがペーサー(pacer)からの信号でトリガーされ、パ
ルスにタグ付けした後3msでパルスシーケンスの画像化が開始する。画像スキ
ャン変数は、以下の通りである:TR=6.5ms、TE=2.1、読み出し帯
域=±3.2kHz、320mmの視野、256×96の獲得マトリックス(a
cquisition matrix)、部分的エコー、動画フレーム当たり2
つの読み出しおよび6mmスライスの厚さ。
のシーケンスは、心臓が収縮してる間に14ms間隔で獲得された。画像は、断
片化されたk空間獲得(k−space acquisition)により息を
止めてる間に獲得された。拍動リードの位置に近い基底面(basal pla
ne)において獲得された画像が使用された。図11は、視覚化のためにLVの
周囲の関心ある領域でクロップされた結果として生じる画像を示す。強い初期の
収縮は、約5時(5 o’clock)に拍動リードの近くで見られ得る。フレ
ーム4〜8において、中隔壁(septal wall)が外向きに曲がってる
のが見られ、それは中隔領域への電気的な活性化信号の遅延により引き起こされ
るプリストレッチング(pre−streching)と呼ばれる異常な動きで
ある。この後、LV心筋全体が、収縮期を通じてほとんど正常に連続的に収縮す
ることを経験する。
通過フィルタを使用して水平タグ付画像シーケンスから計算された。このフィル
タを90°回転させたバージョンが、垂直HARP画像を計算するために使用さ
れた。HARPトラッキング(アルゴリズム2)を示すために、1つは垂直であ
り、残りは水平である2つのタグ付の線が選択され、それぞれのタグ点の収集が
手動で選択された。これらの点の位置、および後に3回にわたる時点で個々にト
ラッキングされたこれらの位置が、図12に示される。1つだけ例外はあるが、
全ての点は、タグを消去した後にさえ点を見るのを期待する箇所にトラッキング
された。詳細には、正常な収縮の内部「反り(bowing)」および異常なプ
リストレッチングの外部「反り」が、HARPトラッキングによって非常によく
取り込まれる。間違ってトラッキングされた点のみが、図12(d)の画像の上
部に見られ得る。画像を注意深く調べることによって、平面の外での動きにより
、第1の時間のフレームのLVの上部にあるタグ付の水平線が長い間に消えるこ
とがわかる。改良によって解決されない間、このタイプの問題は、心筋の境界に
それほど近くない点を選択することによって避けられ得る。
ルを計算するためのポテンシャルを示すために、1画素(1.25mm)によっ
て分けられる点の5×5グリッドが、図13(a)に示されるように、LVの前
横方向側の4つのタグ付の線によって境界を付けられた領域に配置された。これ
らの点は、全ての画像シーケンスを介してそれぞれにトラッキングされた。図1
3(b)〜(c)は、時間フレーム1,5,10,20でそれらの位置を示す拡
大図を示す。トラッキングされた点に対する副画素の解像度がより後の画像にお
いて明瞭に示され、ひずみの基本的なローカルパターン(local patt
ern)を明瞭に見ることができる。四角形からひし形というグリッドの明確な
発展によって、正常な心臓が動きする際に動径方向に厚くなることおよび円形方
向に短縮することの両方が大変明瞭に示される。このデータからひずみテンソル
を計算するために差分が容易に使用され得ることはトラッキングされた点の規則
性から明らかである。これから、主なひずみ(principle strai
nts)の局所領域変化および方向を含む動きに関する様々な量が計算され得る
。
使用して計算された。ユーザインターフェースを使用して、心外膜および心内膜
の円(circles)が定義された。収縮末期にLVはほとんど円形なので、
図11の最後の画像が、これらの円を定義するように使用された。結果として生
じる3つの円および定義された八分円が、図14(a)に示される。それぞれの
円の16個の点が、拡張末期まで後方にトラッキングされ、結果として図14(
b)に示される形となった。歪んだ状態のシーケンス全体が、図14(c)に示
される。このシーケンスから、LVの断面積の形状がいくらか伸張し始めるが、
急速に円形になり、次いでほとんど動径方向に収縮するのがわかる。そのような
タグのジャンプが、1つ以上の時間フレームにおいて非常に歪んだ輪郭を描くの
で、図14に誤ってトラッキングされた点はないことを確認にするのは容易であ
る。
れるように図14(c)に示されたトラッキングされた点から計算された。各八
分円における動径方向のひずみについての時間的発展が、図15(a)で示され
る。正の値は、心筋の厚みが増加するのを示し、負の値はその厚みが薄くなるの
を示す。初期の時点で心筋の厚みが増加することは八分円3〜6のみで明らかで
あり、八分円8、1および2は厚みが薄くなるのを示す。これは、初期に拍動リ
ード近くの心筋で生じる大きな収縮と反対の壁(wall)での心筋のプリスト
レッチングの両方の直接的な表現である。時間フレーム5〜10の間で、再収縮
して拡張末期で厚さが最も大きく増加するその前に、八分円5〜7の拍動リード
近くの心筋は弛緩する。心内膜の厚さの増加がより大きい八分円7を除いて、心
外膜と心内膜の動径方向のひずみの間にほとんど違いはない。
形方向に伸張するのを示し、負の値は収縮を示す。これらのプロットは、動径方
向のひずみのプロフィールに示されるような同一の一般的な振る舞いを示す。八
分円4〜6における初期の時点での収縮が短縮(shortening)として
示され、これと同じ期間内に八分円8、1および2は大きく伸張してるのを示す
。いくらかの期間の後に、全ての心筋組織は収縮性の短縮(contracti
le shortening)を示す。これらのプロットはまた、心外膜より大
きく短縮する中間壁に比べて一貫してより大きく心内膜が短縮を有するのを示す
。これは、収縮する間の左心室の心筋の公知の振る舞いと一致する(Clark
et al., 「Circumferential myocardial
shortening in the normal human left
ventricle.」Cric.,84:67−−74,1991)。
験してる正常な男性ボランティア(27歳)の心臓の動きのHARPトラッキン
グの使用を含む。この実施例は、ドブタミンの応力下のこの心臓の非常に急速な
動きおよび誤ってトラッキングされた点を補正するためのHARP改良の使用を
含む。
得された。SPAMMタグが、7mmで分けられた画像平面に直交する飽和平面
(saturation planes)を作成するために拡張末期で生成され
た。垂直タグおよび水平タグを備える画像の2つのセットが、個別に息止めして
る際に獲得された。4つのスライスが獲得されるが、この実施例では中間壁の基
本的なスライスのみが使用される。スキャナの設定は以下のようである:視野3
6cm、タグ分離7mm、8mmのスライスの厚さ、TR=6.5ms、E=2
.3ms、15°チップ角度、256×160画像マトリックス、動画フレーム
当たりの5位相−符号化されたビュー。
ップされた結果として生じる水平タグおよび垂直タグ付の画像を示す。これらの
画像に生じる輪郭は、HARPトラッキングを使用する時間内に、心外膜および
心内膜の円を第1の画像に手動で配置し、それらを順方向にトラッキングするこ
とによって生成された。第1および第2の時間フレームの間の大きな動きのため
に、前自由壁(anterior free−wall)のいくつかの点が、第
2のフレームにおいてトラッキングされなかった。基礎的なHARPトラッキン
グアプローチは、現在のフレームにおける初期化として以前にトラッキングされ
た点を使用するので、残りのフレームにおいて、この残りは補正されなかった。
隔壁内の手動で識別された3つのアンカー(各円に1つ)を使用してHARP改
良を適用する結果が、図16(b)に示される。見て調べると、改良された結果
により、各時間フレームで予想される箇所にトラッキングされた点が配置された
のが理解され得る。タグジャンプが排除される。
発展が計算された。その結果は、図16(c)に示される。これらのプロットは
、自由壁において最大かつ持続性のある短縮を伴うLV中のかなり画一的な短縮
を示す。これらはまた、正常な筋肉において一般に生じるような心内膜中に生じ
る最大の短縮を示す。
られる公知の技術の基礎である適切なフィルタリングと比較してHARPトラッ
キングの精度についての2つの予備調査が行われた(T.S. Denney,
「Identification of myocardial tags i
n tagged MR images without prior kno
wledge of myocardial contours,」In J.
Duncan and G.Gindi,editors,Proc.Inf
.Proc.Med.Imag.,pages 327−−340,1997)
。画像のノイズ比(CNR)との対比に依存して、FindTagsの精度が0
.1〜0.2画素の範囲であることが、理論およびファントム確証(phant
om validation)の両方によって示される(Mooreら.,「T
agged MR imaging in a deforming phan
tom:photographic validation,」Radiolo
gy,190:765−−769,1994)。その結果により、HARPの精
度が、適合したフィルタリングアプローチより同じまたはより優れていることが
示される。
に注意が向けられた。FindTagsは、垂直なタグ付の短軸のデータセット
からの27個の画像において輪郭(心外膜および心内膜)およびタグ付の線の両
方を推定するために使用された。その27個の画像は、LV内の3つの長軸方向
の位置からの9つの画像シーケンスを含む。HARP画像は、図8(b)に示さ
れたものと同様の第1の高調波および帯域通過フィルタを使用してこれら9つの
画像から生成された。理論によると、タグの底は、πラジアンの位相角度で配置
され、それによって、HARP画像のπ等輪郭(isocontours)が、
FindTagsによって識別されるタグ点に非常に近接することになると予測
される。図17(a)では、タグのかなりの消去が生じる心室の中間についての
7つの時間フレームで撮られた画像におけるこれらのHARP等輪郭に、Fin
dTagsからのタグ点が重ね合わされる。
かな違いしかない。HARPは非常に小さく滑らかな結果を生成し、主な違いが
、自由壁(3時)上のHARPによる結果の周りの小さな揺らぎとなる。どれが
視覚的に満足に値するかを決定するのは難しい。FindTagsタグ点の収集
全体を使用して計算された平均HARP角度が、非常にπ(3桁の数(thre
e significant digits))に近く、それによってタグの底
がπというHARP角度を有するという理論的予測を証明する。ローカル位相非
ラッピング(unwrapping)は、HARP画像のラッピング(wrap
ping)アーティファクトの位置に正確にπが対応する場合、この平均角度を
計算するために使用された。
近接のHARPπラジアン等輪郭(π radian isocontour)
との間の距離が計算された。全時間と空間的位置にわたる平方二乗平均(rms
)距離が、約0.15画素である。基底の、中間の心室の、短軸の先端の画像に
対して、時間の関数として平均距離がプロットされて、図17(b)に示される
。FindTagsは現実を表す場合、HARPは、約0.15〜0.25画素
誤りを有すると結論され得る。これらの平均距離はFindTagsの固有の誤
りと同じオーダーであるので、この結論は引き出すことは正しい。実際、HAR
Pは大きく誤りを減少させ得、その場合には、これらの距離はFindTags
の固有の誤りを簡潔に表示してる。しかし、FindTagsの誤りに加えて、
HARPがこの誤りを有し、その場合にはHARPは約0.3画素平均誤りを有
する可能性もまたある。HARPトラッキング誤りは、FindTagsの誤り
とほとんど同じであると結論される。
グ付の線自体よりタグ付の線の交差の認識の方がより容易になる。結果として、
この実験では、HARPトラッキングは、FindTagsを使用するタグ付の
線の交差の推定と比較された。本明細書中で説明された拍動する犬の心臓のデー
タを使用して、FindTagsは、20の画像基底短軸画像シーケンス(20
−image basal short−axis image sequen
ce)における垂直と水平の両方の全てのタグ付の線の位置を計算するように使
用された。心内膜および心外膜の輪郭を使用して、またFindTagsの使用
を推定され、心筋内部に落ち込む(falling)タグ付の線の交差点が計算
された。HARPトラッキングは、次いで、ターゲットベクトルa=[ππ]T
ラジアンを求め、第1の時間フレームを含むこれら各々の点で実行された。
距離が、図18に時間の関数として示される。これらの誤りは、以前の実験より
いくらか大きい。これは、1つの線の代わりに2つの線を求める場合、誤りが
いは、異なる実験的なセットアップと画像手順によって説明され得る。長時間に
わたって距離が増加するという一般的な傾向はまた、タグが消えるにつれて信号
−ノイズ比が降下するという理由から予測される。
〜7において「こぶ(hump)」が生じることである。HARPで使用される
帯域通過フィルタの設定(placement)により、これを説明することが
可能である。収縮期の間に正常な心臓な動きは収縮するので、基本的なタグ高調
波(the fundamental tag harmonic)のローカル
周波数は一般的に増加する。また、タグが消えるにつれて第1の高調波の周りの
スペクトルとの干渉が生じるので、エネルギー的にDCスペクトルピークが増加
する。それゆえ、わずかに高い周波数に帯域通過フィルタを設定したり、収縮に
関連するその高い周波数スペクトルを同時に取り込んだり、DCスペクトルピー
クからの干渉を避けるために通常は使用される。しかし、フレーム4〜7では、
図11および図15に見られ得るように、プリストレッチング現象が最もはっき
りしており、これによりタグのローカル周波数は増加ではなく減少する。わずか
に高い違いの原因およびおそらく誤りを増加させる原因となるスペクトルのこの
部分は、フィルタでカットオフされ得る。
法(Lagrangian strain computations)を包含
する全ての計算が、MATLAB(The Mathworks、Natick
MA)を使用して400MHzのIntel Pentium(R) IIプ
ロセッサで実行された。MATLABコードを最適化(例えば、ループをなくす
こと)しない事実にも関わらず、HARP処理は、知られてる他の方法と比較し
て非常に速い。拍動する犬の心臓に対して、20の垂直HARP画像および20
の水平HARP画像の計算が約30秒で終了した。LVの心筋内への心内膜およ
び心外膜の円の配置によって、人間に相互作用があるのは約20秒である。20
の時間フレーム全てを介するこのプロセスにより定義された48の点のトラッキ
ングは、わずか約5秒で終了し、ラグランジアンひずみの計算もまた、わずか約
5秒で終了する。画像からひずみを計算する時間全体は、ボタンにクリックし、
画像を調整する時間を含んでわずか約2秒である。さらに一定の工程の重要な効
率化が可能になる。
定(definition)もまた、HARP分析の処理時間全体を増加する。
標準のスキャンでは、これらの時間は無視できる程であり、画像シーケンスは、
自動的に構築され得、現行の帯域通過フィルタが使用され得る。特定のスキャン
または実験手順では、いくつかの付加的な時間が、HARP処理のためにデータ
を準備する際に費やされ得る。当業者は、一般的にこれらの付加的な工程を30
分未満で実行し得、便利なユーザインターフェースによりさらにこの時間が低減
する。HARPは、さらにわずかな有効化および最適化を伴う臨床的使用に有用
である。
elocity field)であるオプティカルフローを計算するための方法
を以前に提案した(Osman et al .,「Imaging hear
t motion using harmonic phase MRI,」O
ctober 1998)。このHARPオプティカルフロー方法では、材料の
性質としてHARP角度のアイデアおよび見かけ上の動きの観念を用いたが、2
つのHARP値を共有して点を繰り返して求めはしなかった。代わりに、各画素
で単純な2×2のマトリックスインバースを使用して次の画像フレームにおいて
直接におよび急速に各画素の大体の位置を計算するための複数の制約光フローの
アイデアを用いる。画像における全ての画素に適用され、次の時間フレームにの
みトラッキングされる場合、HARPトラッキングは、より優れた精度のみにつ
いて本質的に同じ結果を与える。HARPオプティカルフロー方法はまた、さら
に多くの時間おそらく4〜5倍の長い時間が必要となる。HARPオプティカル
フロー方法は、HARPトラッキングを初期化すること、または視覚化動きのた
めの非常に高速で高密度の動き場を生成することあるいは高密度のメッシュ上の
他の動き量を計算することについて有用であり得る。
トな2次元のHARPトラッキング法を開発するために利用された。点は、左心
室の円形方向および動径方向のラグランジアンひずみを直接に計算するために座
標系においてトラッキングされた。拍動した犬の心臓についての実験は、異常な
動きをトラッキングし、以前に報告された分析と一致するひずみを計算するHA
RPの能力を示した。誤ってトラッキングされた点の補正についての改良技術を
発展させ、ドブタミンの応力を経験する正常な人間の心臓について示された。最
後に、予備的な誤り分析が行われ、FindTags、標準テンプレートマッチ
ング法と比較するとHARPトラッキング法は非常に利点を有することを示した
。HARPトラッキングおよびラグランジアンひずみ分析が、計算的に高速で、
さらに適切な有効化の後では臨床的な使用に用いられることが示された。
に参考にしてきたが、本発明の方法は、それに制限されず、広範な多様な産業的
用途および他の用途において使用され得る。
業者には、本発明から逸脱することなく詳細部の多数の変更が為され得ることが
明らかである。
。
る、大きさを示す複雑な画像を示す。
によって生成される電気刺激で鼓動させた犬の心臓のタグ付磁気共鳴画像のシー
ケンスを示す。画像は、左上に示される拡張末期から右下に示される収縮末期ま
での電気刺激で鼓動させた犬の心臓の動きを示す20個の時間フレームである。
た点を示す。
た点を示す。
れた点を示す。
れた点を示す。
示す。
心内膜の半径方向のひずみ(実線)の時間的進展を示す。
壁の円周方向のひずみ(破線)、および心内膜の円周方向のひずみ(実線)の時
間的進展を示す。
、短軸画像および追跡された円を含むドブタミン誘因ストレスを受ける通常の人
間の心臓を示す。
す。
gs(黒点)からのタグ点を有するタグ付画像を有する通常の人間のデータを使
用してFindTagsと比較してのHARP精度を示す。
を示す。
のプロットを示す。
Claims (27)
- 【請求項1】 磁気共鳴画像法により対象物の動きを測定する方法であって
、該方法は、 該対象物の関心のある領域を空間的に変調するためにパルスシーケンスを適用
する工程と、 該空間的に変調された対象物のフーリエ領域から少なくとも1つの第1のスペ
クトルピークを獲得する工程と、 該獲得された第1のスペクトルピークの逆フーリエ変換情報を計算する工程と
、 該スペクトルピークのそれぞれから第1の高調波位相画像を計算する工程と、 第2のスペクトルピークのそれぞれから第2の高調波位相画像を作成するため
に異なる時間に上記プロセスを繰り返す工程と、 該第1および第2の高調波位相画像からひずみを判定する工程と、 を包含する方法。 - 【請求項2】 心臓の動きをトラッキングするために前記方法を用いる工程
を包含する、請求項1に記載の方法。 - 【請求項3】 円形ラグランジアンひずみを判定するために前記方法を用い
る工程を包含する、請求項2に記載の方法。 - 【請求項4】 動径方向ラグランジアンひずみを判定するために前記方法を
用いる工程を包含する、請求項2に記載の方法。 - 【請求項5】 2次元画像を提供するために前記方法を用いる工程を包含す
る、請求項1に記載の方法。 - 【請求項6】 動いている人間の心臓に前記方法を用いる工程を包含する、
請求項1に記載の方法。 - 【請求項7】 SPAMMパルスシーケンスを前記パルスシーケンスとして
使用する工程を包含する、請求項1に記載の方法。 - 【請求項8】 3次元画像を提供するために前記方法を使用する工程を包含
する、請求項1に記載の方法。 - 【請求項9】 心内膜の、心外膜のおよび中間壁のひずみの少なくとも1つ
を判定する際に前記方法を使用する工程を包含する、請求項3に記載の方法。 - 【請求項10】 心内膜の動径方向の、および心外膜の動径方向のひずみの
少なくとも1つを判定する際に前記方法を使用する工程を包含する、請求項4に
記載の方法。 - 【請求項11】 心筋におけるひずみを判定するために前記方法を使用する
工程を包含する、請求項2に記載の方法。 - 【請求項12】 2つの異なる時間の点間の前記間隔を判定することにより
ひずみを判定するために、前記方法を使用する工程を包含する、請求項1に記載
の方法。 - 【請求項13】 2つの異なる時間の点間の距離の増加を測定するために、
前記方法を使用する工程を包含する、請求項1に記載の方法。 - 【請求項14】 2つの異なる時間の点間の距離の減少を測定するために、
前記方法を使用する工程を包含する、請求項1に記載の方法。 - 【請求項15】 前記心臓の左心室内部の動きを測定するために前記方法を
使用する工程を包含する、請求項6に記載の方法。 - 【請求項16】 約拡張末期で前記動きの測定を開始する工程を包含する、
請求項6に記載の方法。 - 【請求項17】 水平方向と垂直方向の両方にタグ付画像を使用する工程を
包含する、請求項1に記載の方法。 - 【請求項18】 画像平面の前記点の見かけ上の動きをトラッキングするこ
とにより前記ひずみを判定する工程および該トラッキングされた点からのラグラ
ンジアンひずみを判定する工程を包含する、請求項12に記載の方法。 - 【請求項19】 タグ付磁気共鳴画像のCINEシーケンスを介して前記見
かけ上の動きをトラッキングする工程を包含する、請求項18に記載の方法。 - 【請求項20】 前記動きをトラッキングする際に改良手順を使用する工程
を包含する、請求項2に記載の方法。 - 【請求項21】 前記改良手順において、比較的小さな動きをトラッキング
するための基準として有するアンカー点を使用する工程を包含する、請求項20
に記載の方法。 - 【請求項22】 前記改良手順において、次の隣接する点から1画素未満分
離されるそれぞれの点の点シーケンスを使用する工程を包含する、請求項21に
記載の方法。 - 【請求項23】 約動径方向経路において、前記改良手順を実行する工程を
包含する、請求項22に記載の方法。 - 【請求項24】 前記高調波位相画像を計算する際に帯域通過フィルタを使
用する工程を包含する、請求項1に記載の方法。 - 【請求項25】 円形伸張の指標として正のひずみ値を使用し、円形収縮の
指標として負のひずみ値を使用する工程を包含する、請求項3に記載の方法。 - 【請求項26】 動径方向の厚みの増加の指標として正のひずみ値を使用し
、動径方向の厚みの減少の指標として負のひずみ値を使用する工程を包含する、
請求項4に記載の方法。 - 【請求項27】 前記ひずみの判定を実行する際に、実質的には同時に起こ
る前記第1の高調波位相および第2の高調波位相を使用する工程を包含する、請
求項1に記載の方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13059599P | 1999-04-22 | 1999-04-22 | |
US60/130,595 | 1999-04-22 | ||
PCT/US2000/010232 WO2000064344A1 (en) | 1999-04-22 | 2000-04-14 | Cardiac motion tracking using cine harmonic phase (harp) magnetic resonance imaging |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2002541963A true JP2002541963A (ja) | 2002-12-10 |
JP2002541963A5 JP2002541963A5 (ja) | 2007-05-31 |
JP4679729B2 JP4679729B2 (ja) | 2011-04-27 |
Family
ID=22445423
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2000613339A Expired - Lifetime JP4679729B2 (ja) | 1999-04-22 | 2000-04-14 | Cine高調波位相(harp)磁気共鳴画像法を用いた心臓動きのトラッキング |
Country Status (8)
Country | Link |
---|---|
US (1) | US6892089B1 (ja) |
EP (1) | EP1171031B1 (ja) |
JP (1) | JP4679729B2 (ja) |
AT (1) | ATE545873T1 (ja) |
AU (1) | AU780349B2 (ja) |
CA (1) | CA2368237A1 (ja) |
MX (1) | MXPA01010656A (ja) |
WO (1) | WO2000064344A1 (ja) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010104614A (ja) * | 2008-10-30 | 2010-05-13 | Toshiba Corp | 画像処理装置、磁気共鳴イメージング装置および画像管理システム |
JP2010178927A (ja) * | 2009-02-05 | 2010-08-19 | Toshiba Corp | 磁気共鳴イメージング装置 |
JP2012516179A (ja) * | 2009-01-29 | 2012-07-19 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 生物学的構造における時変パラメータの視覚化 |
JP5237957B2 (ja) * | 2007-10-18 | 2013-07-17 | 株式会社日立メディコ | 磁気共鳴撮影装置 |
JP2016101470A (ja) * | 2014-11-27 | 2016-06-02 | 株式会社東芝 | 医用画像処理装置、医用画像処理方法および医用イメージング装置 |
US9779510B2 (en) | 2014-08-26 | 2017-10-03 | Toshiba Medical Systems Corporation | Medical image processing apparatus, medical image processing method and medical image device |
KR101784502B1 (ko) | 2016-09-05 | 2017-11-06 | 성균관대학교산학협력단 | 동적 태그 자기 공명 영상 장치 및 방법 |
Families Citing this family (53)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6597935B2 (en) * | 2000-02-10 | 2003-07-22 | The Johns Hopkins University | Method for harmonic phase magnetic resonance imaging |
US6723537B2 (en) | 2001-05-18 | 2004-04-20 | Rigel Pharmaceuticals, Incorporated | Directed evolution of protein in mammalian cells |
US6961454B2 (en) | 2001-10-04 | 2005-11-01 | Siemens Corporation Research, Inc. | System and method for segmenting the left ventricle in a cardiac MR image |
DE10160205A1 (de) * | 2001-12-07 | 2003-06-18 | Philips Intellectual Property | Verfahren zur Rekonstruktion eines 3D-Bilddatensatzes eines Untersuchungsbereiches |
JP4176989B2 (ja) * | 2001-12-12 | 2008-11-05 | 株式会社東芝 | 磁気共鳴診断装置 |
JP4060615B2 (ja) * | 2002-03-05 | 2008-03-12 | 株式会社東芝 | 画像処理装置及び超音波診断装置 |
US7835783B1 (en) | 2002-04-22 | 2010-11-16 | The United States Of America As Represented By The Department Of Health And Human Services | Magnetic resonance imaging methods and apparatus for time-series motion tracking with inversion recovery compensation |
AU2003256925A1 (en) | 2002-07-29 | 2004-02-16 | William Gregory Hundley | Cardiac diagnostics using wall motion and perfusion cardiac mri imaging and systems for cardiac diagnostics |
US7081749B2 (en) * | 2004-08-02 | 2006-07-25 | Albert Macovski | Magnetic resonance imaging using demodulated k-space segments |
SE526438C2 (sv) * | 2003-04-01 | 2005-09-13 | Sectra Imtec Ab | Metod och system för mätning i en dynamisk sekvens av medicinska bilder |
US8064979B2 (en) * | 2003-06-09 | 2011-11-22 | General Electric Company | Tempero-spatial physiological signal detection method and apparatus |
US7495438B2 (en) * | 2005-05-04 | 2009-02-24 | John Hopkins University | Three dimensional magnetic resonance motion estimation on a single image plane |
US7352370B2 (en) * | 2005-06-02 | 2008-04-01 | Accuray Incorporated | Four-dimensional volume of interest |
DE102005051323B4 (de) * | 2005-10-26 | 2008-01-31 | Siemens Ag | Verfahren und Gerät zur Magnet-Resonanz-Bildgebung auf Basis einer Gradientenechosequenz |
US7911207B2 (en) | 2005-11-16 | 2011-03-22 | Board Of Regents, The University Of Texas System | Method for determining location and movement of a moving object |
US7498811B2 (en) * | 2005-11-16 | 2009-03-03 | Macfarlane Duncan L | Apparatus and method for patient movement tracking |
US7977942B2 (en) * | 2005-11-16 | 2011-07-12 | Board Of Regents, The University Of Texas System | Apparatus and method for tracking movement of a target |
WO2007087702A1 (en) * | 2006-01-31 | 2007-08-09 | Canadian Space Agency | Method and system for increasing signal-to-noise ratio |
JP5105848B2 (ja) * | 2006-02-06 | 2012-12-26 | 株式会社東芝 | 磁気共鳴イメージング装置および磁気共鳴イメージング装置における撮影条件設定方法 |
CN100570393C (zh) * | 2006-02-06 | 2009-12-16 | 株式会社东芝 | 磁共振成像装置及磁共振成像方法 |
WO2007099525A2 (en) * | 2006-03-03 | 2007-09-07 | Medic Vision - Brain Technologies Ltd. | System and method of automatic prioritization and analysis of medical images |
US8081806B2 (en) * | 2006-05-05 | 2011-12-20 | General Electric Company | User interface and method for displaying information in an ultrasound system |
US7813537B2 (en) * | 2006-05-15 | 2010-10-12 | Siemens Medical Solutions Usa, Inc. | Motion-guided segmentation for cine DENSE images |
EP2023812B1 (en) | 2006-05-19 | 2016-01-27 | The Queen's Medical Center | Motion tracking system for real time adaptive imaging and spectroscopy |
US8073523B2 (en) * | 2006-06-06 | 2011-12-06 | California Institute Of Technology | Method for obtaining strain from radially-tagged magnetic resonance imaging (MRI) |
US8290266B2 (en) * | 2006-08-24 | 2012-10-16 | Agency For Science, Technology And Research | Model-based method and system for image segmentation and modelling |
WO2008105690A1 (en) * | 2007-02-28 | 2008-09-04 | St. Jude Medical Ab | An implantable heart analysing device system and method |
TW200930850A (en) * | 2008-01-03 | 2009-07-16 | Green Energy Technology Inc | Cooling structure for body of crystal growth furnace |
US8940315B2 (en) * | 2008-04-18 | 2015-01-27 | Medtronic, Inc. | Benzodiazepine formulation in a polyorthoester carrier |
US8390291B2 (en) * | 2008-05-19 | 2013-03-05 | The Board Of Regents, The University Of Texas System | Apparatus and method for tracking movement of a target |
WO2010084446A1 (en) * | 2009-01-23 | 2010-07-29 | Koninklijke Philips Electronics N.V. | Cardiac image processing and analysis |
JP5491836B2 (ja) * | 2009-01-30 | 2014-05-14 | 株式会社東芝 | 超音波診断装置、超音波画像処理装置、医用画像診断装置及び医用画像処理装置 |
US8538104B2 (en) * | 2009-04-30 | 2013-09-17 | Kyoto University | Method of analyzing motion of imaging target by means of tagged MR images, and MRI device |
DE102009021492B4 (de) * | 2009-05-15 | 2012-01-12 | Siemens Aktiengesellschaft | Verfahren zur Bestimmung einer MR-Relaxationszeit im Herzmuskel bei einer Magnetresonanzuntersuchung, Magnetresonanzanlage, Computerprogrammprodukt und elektronisch lesbarer Datenträger |
US8478012B2 (en) * | 2009-09-14 | 2013-07-02 | General Electric Company | Methods, apparatus and articles of manufacture to process cardiac images to detect heart motion abnormalities |
WO2012006578A2 (en) | 2010-07-08 | 2012-01-12 | The Regents Of The University Of California | End-to-end visual recognition system and methods |
EP2670302A4 (en) | 2011-02-03 | 2015-04-15 | Univ California | CURVILIGNE MARKING IN CARDIAC MRI |
US20130035588A1 (en) * | 2011-08-03 | 2013-02-07 | Siemens Corporation | Magnetic resonance imaging for therapy planning |
EP2747641A4 (en) | 2011-08-26 | 2015-04-01 | Kineticor Inc | METHOD, SYSTEMS AND DEVICES FOR SCAN INTERNAL MOTION CORRECTION |
US20140371574A1 (en) * | 2011-09-16 | 2014-12-18 | University Of Pittsburgh - Of The Commonwealth System Of Higher Education | Noncontact electrophysiological measurement and imaging of the heart |
CN111154375A (zh) | 2012-08-29 | 2020-05-15 | 赫普有限公司 | 包含空心玻璃球和导电颜料的防腐蚀锌底漆涂料组合物 |
US9305365B2 (en) | 2013-01-24 | 2016-04-05 | Kineticor, Inc. | Systems, devices, and methods for tracking moving targets |
US10327708B2 (en) | 2013-01-24 | 2019-06-25 | Kineticor, Inc. | Systems, devices, and methods for tracking and compensating for patient motion during a medical imaging scan |
US9717461B2 (en) | 2013-01-24 | 2017-08-01 | Kineticor, Inc. | Systems, devices, and methods for tracking and compensating for patient motion during a medical imaging scan |
US9782141B2 (en) | 2013-02-01 | 2017-10-10 | Kineticor, Inc. | Motion tracking system for real time adaptive motion compensation in biomedical imaging |
US10004462B2 (en) | 2014-03-24 | 2018-06-26 | Kineticor, Inc. | Systems, methods, and devices for removing prospective motion correction from medical imaging scans |
DE102014209803B4 (de) * | 2014-05-22 | 2016-03-03 | Siemens Aktiengesellschaft | 4D-Geschwindigkeitsprofil-Bildaufnahme mit einer Magnetresonanzanlage |
WO2016014718A1 (en) | 2014-07-23 | 2016-01-28 | Kineticor, Inc. | Systems, devices, and methods for tracking and compensating for patient motion during a medical imaging scan |
US9943247B2 (en) | 2015-07-28 | 2018-04-17 | The University Of Hawai'i | Systems, devices, and methods for detecting false movements for motion correction during a medical imaging scan |
US10716515B2 (en) | 2015-11-23 | 2020-07-21 | Kineticor, Inc. | Systems, devices, and methods for tracking and compensating for patient motion during a medical imaging scan |
CA3020896A1 (en) | 2016-04-27 | 2017-11-02 | Myocardial Solutions, Inc. | Rapid quantitative evaluations of heart function with strain measurements from mri |
US10109050B2 (en) * | 2016-06-01 | 2018-10-23 | Siemens Healthcare Gmbh | Spatiotemporal background phase correction for phase contrast velocity encoded MRI |
US11872019B2 (en) | 2019-06-05 | 2024-01-16 | Myocardial Solutions, Inc. | MRI-derived strain-based measurements and related image data acquisitions, image data processing, patient evaluations and monitoring methods and systems |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1099297A (ja) * | 1996-09-18 | 1998-04-21 | Siemens Corporate Res Inc | タグ付けされたmrイメージにおける心臓壁の動きを再生するためのハイブリッドモデルおよび方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4486708A (en) | 1981-12-21 | 1984-12-04 | Albert Macovski | Selective material projection imaging system using nuclear magnetic resonance |
US5363044A (en) * | 1993-02-09 | 1994-11-08 | Sunnybrook Health Science Center | Motion ghost manipulation in magnetic resonance imaging |
US5910728A (en) | 1996-11-12 | 1999-06-08 | Beth Israel Deaconess Medical Center | Simultaneous acquisition of spatial harmonics (SMASH): ultra-fast imaging with radiofrequency coil arrays |
AU8472398A (en) * | 1997-06-24 | 1999-01-04 | California Institute Of Technology | A method for suppressing noise in measurements |
US6453187B1 (en) * | 1998-08-10 | 2002-09-17 | The Johns Hopkins University | Method of employing angle images for measuring object motion in tagged magnetic resonance imaging |
-
2000
- 2000-04-14 CA CA002368237A patent/CA2368237A1/en not_active Abandoned
- 2000-04-14 AT AT00922265T patent/ATE545873T1/de active
- 2000-04-14 WO PCT/US2000/010232 patent/WO2000064344A1/en active Application Filing
- 2000-04-14 JP JP2000613339A patent/JP4679729B2/ja not_active Expired - Lifetime
- 2000-04-14 AU AU42474/00A patent/AU780349B2/en not_active Expired
- 2000-04-14 MX MXPA01010656A patent/MXPA01010656A/es active IP Right Grant
- 2000-04-14 US US09/959,311 patent/US6892089B1/en not_active Expired - Lifetime
- 2000-04-14 EP EP00922265A patent/EP1171031B1/en not_active Expired - Lifetime
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1099297A (ja) * | 1996-09-18 | 1998-04-21 | Siemens Corporate Res Inc | タグ付けされたmrイメージにおける心臓壁の動きを再生するためのハイブリッドモデルおよび方法 |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5237957B2 (ja) * | 2007-10-18 | 2013-07-17 | 株式会社日立メディコ | 磁気共鳴撮影装置 |
JP2010104614A (ja) * | 2008-10-30 | 2010-05-13 | Toshiba Corp | 画像処理装置、磁気共鳴イメージング装置および画像管理システム |
JP2012516179A (ja) * | 2009-01-29 | 2012-07-19 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 生物学的構造における時変パラメータの視覚化 |
JP2010178927A (ja) * | 2009-02-05 | 2010-08-19 | Toshiba Corp | 磁気共鳴イメージング装置 |
US9779510B2 (en) | 2014-08-26 | 2017-10-03 | Toshiba Medical Systems Corporation | Medical image processing apparatus, medical image processing method and medical image device |
US10489910B2 (en) | 2014-08-26 | 2019-11-26 | Canon Medical Systems Corporation | Medical image processing apparatus, medical image processing method and medical image device |
JP2016101470A (ja) * | 2014-11-27 | 2016-06-02 | 株式会社東芝 | 医用画像処理装置、医用画像処理方法および医用イメージング装置 |
KR101784502B1 (ko) | 2016-09-05 | 2017-11-06 | 성균관대학교산학협력단 | 동적 태그 자기 공명 영상 장치 및 방법 |
WO2018044135A1 (ko) * | 2016-09-05 | 2018-03-08 | 성균관대학교산학협력단 | 동적 태그 자기 공명 영상 장치 및 방법 |
Also Published As
Publication number | Publication date |
---|---|
CA2368237A1 (en) | 2000-11-02 |
US6892089B1 (en) | 2005-05-10 |
MXPA01010656A (es) | 2002-06-04 |
EP1171031A4 (en) | 2006-05-10 |
WO2000064344A1 (en) | 2000-11-02 |
ATE545873T1 (de) | 2012-03-15 |
EP1171031B1 (en) | 2012-02-15 |
AU4247400A (en) | 2000-11-10 |
JP4679729B2 (ja) | 2011-04-27 |
EP1171031A1 (en) | 2002-01-16 |
AU780349B2 (en) | 2005-03-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4679729B2 (ja) | Cine高調波位相(harp)磁気共鳴画像法を用いた心臓動きのトラッキング | |
Osman et al. | Cardiac motion tracking using CINE harmonic phase (HARP) magnetic resonance imaging | |
Osman et al. | Visualizing myocardial function using HARP MRI | |
Osman et al. | Imaging heart motion using harmonic phase MRI | |
Ibrahim | Myocardial tagging by cardiovascular magnetic resonance: evolution of techniques–pulse sequences, analysis algorithms, and applications | |
US9176211B2 (en) | Method for quantitatively mapping myocardial contractile function with magnetic resonance based multiparametric strain analysis | |
US9395431B2 (en) | Multi-contrast delayed enhancement cardiac magnetic resonance imaging | |
US6597935B2 (en) | Method for harmonic phase magnetic resonance imaging | |
US7495438B2 (en) | Three dimensional magnetic resonance motion estimation on a single image plane | |
JP4527881B2 (ja) | 磁気共鳴イメージング方法 | |
CN105283774B (zh) | 在不利的心脏条件下改进受试者心脏成像的系统与方法 | |
US10776998B1 (en) | Method and system for analysis of 3D deformations and regional function of a heart with 3D SinMod | |
JP2003515404A (ja) | 素材の変形及びひずみを検査するための高速3次元磁気共鳴タッギング | |
Maier et al. | Acquisition and evaluation of tagged magnetic resonance images of the human left ventricle | |
Wang et al. | Orthogonal CSPAMM (OCSPAMM) MR tagging for imaging ventricular wall motion | |
Osman | Measuring regional cardiac function using harmonic phase magnetic resonance imaging | |
Augenstein et al. | Magnetic resonance imaging and ventricle mechanics | |
Yang et al. | A multi-scale phase-based optical flow method for motion tracking of left ventricle | |
Ibrahim | Harmonic Phase (HARP) Analysis | |
Lee et al. | Validation of myocardial elastography using MR tagging in normal and abnormal human hearts in vivo | |
CN115137333A (zh) | 一种无外加门控的五维血流磁共振成像方法 | |
WO2004042656A1 (en) | A method for processing medical images | |
Simpson | Development of a high temporal resolution spiral 3D phase velocity mapping technique for analysis of regional myocardial mechanics | |
Sampath et al. | Imaging the evolution of three-dimensional myocardial strains using a fast MR imaging technique | |
Auger | 3D cine DENSE MRI: ventricular segmentation and myocardial stratin analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070330 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20070330 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100126 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20100420 |
|
A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20100427 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20100526 |
|
A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20100602 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100628 |
|
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: 20110106 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20110202 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4679729 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140210 Year of fee payment: 3 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
EXPY | Cancellation because of completion of term |