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
Application number
JP2000613339A
Other languages
English (en)
Other versions
JP4679729B2 (ja
JP2002541963A5 (ja
Inventor
ジェリー エル. プリンス,
ナエル エフ. オズマン,
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Johns Hopkins University
Original Assignee
Johns Hopkins University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Johns Hopkins University filed Critical Johns Hopkins University
Publication of JP2002541963A publication Critical patent/JP2002541963A/ja
Publication of JP2002541963A5 publication Critical patent/JP2002541963A5/ja
Application granted granted Critical
Publication of JP4679729B2 publication Critical patent/JP4679729B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7253Details of waveform analysis characterised by using transforms
    • A61B5/7257Details of waveform analysis characterised by using transforms using Fourier transforms
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/563Image 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/45For evaluating or diagnosing the musculoskeletal system or teeth
    • A61B5/4519Muscles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/563Image 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/56308Characterization of motion or flow; Dynamic imaging
    • G01R33/56333Involving spatial modulation of the magnetization within an imaged region, e.g. spatial modulation of magnetization [SPAMM] tagging
    • YGENERAL 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
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10STECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S128/00Surgery
    • Y10S128/92Computer assisted medical diagnostics
    • Y10S128/922Computer 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

(57)【要約】 本発明は、磁気共鳴画像を使用した心臓等の対象物の動きを測定する方法に関する。対象物について関心のある領域を空間的に変調するようにパルスシーケンスが適用され、空間的に変調された対象のフーリエ領域から少なくとも1つの第1のスペクトルピークが獲得される。獲得された第1のスペクトルピークの逆フーリエ変換情報を計算し、計算された第1の高調波位相画像が各スペクトルピークから判定される。上記プロセスが繰り返され、各第2のスペクトルピークから第2の高調波位相画像が作成され、第1および第2の高調波位相画像からひずみが判定される。

Description

【発明の詳細な説明】
【0001】 (関連出願との相互参照) 本出願は、1999年4月22日に出願された米国特許仮出願シリアル番号第
60/130,595号の利益を主張する。
【0002】 (政府の支持) 本発明は、国立衛生研究所授与番号第R29HL47405号および全米科学
財団授与番号第MIP93−50336号によって少なくとも一部は支持された
【0003】 (発明の背景) 1.発明の分野 本発明は、磁気共鳴画像化を使用して心臓の動きを測定することに関し、より
具体的には、心臓における物質点の追跡およびラグランジアンひずみを計算する
ために磁気共鳴画像化を使用して獲得される高調波位相画像(harmonic
phase image)を使用するプロセスに関する。
【0004】 2.従来技術の説明 興味領域内の核を励起し、そして受信機によって受け取られ、場合によっては
コンピュータ手段によって処理され、その後で結果が記録され、結果のハードコ
ピーを表示または生成する磁気エネルギーの応答放出を引き起こすために働くR
Fパルスおよび傾斜パルスを交互に生成する磁場内に、患者または対象が置かれ
る磁気共鳴画像化の有利な使用は、長い間公知であった。概して、Atalar
−McVeighの米国特許第5,512,825号およびConturo−R
obinsonの米国特許第5,281,914号を参照されたい。その両方は
、本発明の所有者に譲渡され、その開示は、参照により本明細書中に明確に援用
される。
【0005】 磁化の空間変調(SPAMM)アプローチを介して磁気共鳴画像化で、2次元
の心臓壁の動きを画像化する場合に、画像面に直交する方向にある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画像化
についての磁化の空間変調の改良方法」を参照されたい。
【0006】 タグ特徴を探知し、そしてその後で、その特徴を、使用されている引き続いて
の補間を伴う変位およびひずみに関する詳細な動きのマップ内へ結合するために
、画像処理技術を使用することは、磁気共鳴のタグ付に関連して公知である。例
えば、Young等のRadiology誌1992年185号241−247
ページの「磁化の空間変調を伴う3次元の動きおよびゆがみ」、およびMcVe
igh等のRadiology誌1991年180号(3)677−683ペー
ジの「MR画像化での心筋のひずみにおける経壁傾斜の非侵襲性測定」を参照さ
れたい。これらのアプローチは、自動化されていない。なぜならいくつかの手動
による変調作業が必要であるから。
【0007】 磁気共鳴タグ付画像シーケンスに関するオプティカルフローの方法を使用する
ことも公知である。概して、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)の代わりに使用される。画像輝度傾
斜は、画像のペアから推定される一時的な微分と共に、「オプティカルフロー」
と一般的に呼ばれる密な動きの推定を提供するために使用され得る特徴である。
そのようなアプローチは、輝度傾斜は傾斜の方向だけの動きについての情報を含
むという事実を補償するための調整を必要とする。
【0008】 米国特許第5,275,163号は、対象の一部の動きをモニターする場合の
磁気共鳴画像化の使用を開示している。パルスおよび傾斜シーケンスは、空間的
に異なるタグ付パターンおよび差し引き像をタグ付画像を形成するために使用、
ペアで適用される。
【0009】 米国特許第5,352,979号は、スライスまたはボリュームにおけるボリ
ューム要素の位相角応答を観察すること、および外部刺激によって引き起こされ
る摂動前および中に生じる画像化を開示している。
【0010】 米国特許第5,379,766号は、タグ付パターンの検出についての高コン
トラストタグ付グリッドを使用することによって、対象の一部の量的な動きの推
定を開示している。
【0011】 米国特許第5,315,248号および第5,545,993号は、動きの追
跡を開示している。
【0012】 磁気共鳴画像化において平面タグ分析を使用することは公知である。心筋の動
き分析と関連してこのようなアプローチを使用することもまた、公知である。そ
のような従来の技術は、通常、変位ベクトルまたはひずみパターンを介して、こ
れらの画像から動きを抽出することを含み、そして後に補間されるタグ特定およ
び位置推定を含む。
【0013】 位相コントラスト磁気共鳴画像化はまた公知である。これは、速度に敏感な属
性を測定すること、および有限差分を使用して計算されたひずみを有する速度場
を再構築することによって、動きを直接測定するための方法を提供する。これら
2つのアプローチに関する課題の1つは、平面タグ付画像は、自動的には正確に
分析され得ないことである。位相コントラスト画像は、自動的に分析され得るが
、信号対雑音比が低くなりがちであり、許容し難い結果になる。
【0014】 心筋におけるひずみ測定は、心臓病の診断および定量化において重要であり得
ることは公知である。タグ付心臓磁気共鳴画像化におけるこの10年間にわたる
発展で、インビボでの心臓の心筋層の詳細なひずみパターンを測定することが可
能となった。MRタグ付は、心筋層の「タグ」と呼ばれる一時的特徴を生成する
ために、被験者の長軸に沿う磁化を空間的に変調するための特別なパルスシーケ
ンスを使用する。タグ付MRIは、正常なおよび異常な心筋の動きのモデルを開
発しおよび洗練するために、心筋の動きの異常と冠状動脈病の相互関係をよりよ
く理解するために、心筋梗塞の治療の効果を理解するためペースメーカーを使用
する心臓の活動化パターンを分析するために、および虚血性心疾患を初期に探知
するためのストレステストと結びついて、使用されてきた。科学的取り組みが成
功しているにもかかわらず、タグ付MRIは、長い画像化ならびに後の処理時間
、画像化中の患者への不適切なアクセスおよび臨床医および同僚によるMR処理
の恩恵の理解の欠如のせいで、普通の臨床使用になることが遅かった。
【0015】 (発明の要旨) 本明細書中に使用される「角度画像」という言葉は、SPAMMタグ付磁気共
鳴画像における孤立したスペクトルピークに対応する画像の位相を言う。この言
葉はまた、高調波位相画像またはHARP画像と本明細書中では呼ばれる。
【0016】 本発明は、磁気共鳴画像化による対象の動きの測定方法を含み、パルスシーケ
ンスが、対象の空間変調された興味領域に適用され、そして空間変調された対象
のフーリエ領域からの少なくとも1つの第一のスペクトルピークが獲得される。
獲得された第一のスペクトルピークの逆フーリエ変換情報が計算され、そして第
一の高調波位相画像が、スペクトルピークの各々から計算される。このプロセス
は、各第二のスペクトルピークから第二の高調波位相画像を生成するために、異
なる時間に関して繰り返され、そしてひずみが、第一および第二の高調波画像か
ら決定される。
【0017】 本方法は、心臓の動きを追跡するのに特に有用であり、そして1つのアプロー
チにおいて、円周方向および半径方向のラグランジアンひずみ両方を決定するこ
とを含む。好適な実施形態において、SPAMMパルスシーケンスが使用される
【0018】 1つのアプローチにおいて、タグ磁気共鳴画像のCINEシーケンスを介して
、明白な動きの追跡が使用される。
【0019】 磁気共鳴画像化を使用して心臓の動きを追跡する改良された手段を提供するこ
とが、本発明の目的である。
【0020】 正確におよび高速で動作するそのようなシステムを提供することが本発明のさ
らなる目的である。
【0021】 点の動きを追跡する高調波位相磁気共鳴画像化を使用する方法を提供し、そし
てそれによって異なる時間に取られた点の画像のシーケンスを比較することによ
って心臓の部分内のひずみの決定を容易にすることが、本発明のさらなる目的で
ある。
【0022】 時間を追って物質点の動きを正確に追跡する高調波位相画像シーケンスを使用
して、心臓の解剖学的組織および心臓と共に動くタグ付特徴の両方の画像化を可
能にし、そしてそれによってラグランジアンひずみの決定を可能にするようなシ
ステムを提供することが、本発明のさらなる目的である。
【0023】 対象のタグ付磁気共鳴画像を使用して、対象の動きの素早いそして正確な視覚
化のための改良された方法を提供することが本発明の目的である。
【0024】 SPAMMタグ付磁気共鳴画像における孤立したスペクトルピークを使用する
ような方法を提供することが、本発明の別の目的である。
【0025】 角度画像が、フーリエ変換情報の2つ以上のスペクトルピークから獲得され、
そして面ひずみまたはテンソルひずみの計算を提供するために使用されるような
システムを提供することが本発明のさらなる目的である。
【0026】 角度画像を使用してそのような計算が自動的におよび素早く行われ得るように
することが、本発明のさらなる目的である。
【0027】 従来の磁気共鳴画像を生成する必要なくフーリエデータから直接角度画像を生
成することが本発明のさらに別の目的である。
【0028】 本発明のこれらおよび他の目的が、本明細書に添付の図面を参照して本発明の
次の説明から、より十分に理解される。
【0029】 (好適な実施形態の説明) 本明細書中に使用される場合、「患者」という言葉は、人類を含む動物界の生
きている一員を意味する。
【0030】 本明細書中に使用される場合、「対象」という言葉は、本発明の方法によって
動きについてモニタされている患者ならびに他のすべての項目およびそれらの部
分を含む。医療使用の中には、心壁、心筋および心腱の動きを測定する場合の使
用がある。モニタされる対象は、例えば血液または脳脊髄液などの流体、あるい
は固体または半固体、あるいはそれらの組合せであり得る。
【0031】 図1−7(b)の次の記載は、本出願において開示されそして請求された本発
明に関する背景的情報を提供し、そして本発明者らの米国特許出願シリアル番号
第09/131,589号の主題である。この従来の出願は、物理的に仮出願シ
リアル番号第60/130,595の一部であって、そしてそこから本出願へ移
行される。図1−7(b)に関連しての図8−18の記載中の「角度画像」の表
現は、「高調波位相画像」または「HARP画像」と表現する。
【0032】 タグ付磁気共鳴画像化(MRI)は、虚血または梗塞形成によって引き起こさ
れる心臓の動きの欠陥を探知およびモニタする標準的アプローチに急速になって
いる。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画像化のための帯域通過オプティカルフロー」もまた参照されたい。
【0033】 平面のゆがみの密な推定は、従来の特徴抽出および光フロー法を使用すること
なく、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パターンについての信号モデルおよび動き
を基本的な搬送周波数の角度変調と解釈することに依存する。本発明は、非常に
小さな変位(エラー変位など)を推定し、タグ線を統合し、およびオプティカル
フローを計算する場合に直接有用であり得る。
【0034】 SPAMMタグ付磁気共鳴画像は、フーリエ領域における明白なスペクトルピ
ークの収束を有する。各スペクトルピークは、特定の方向の動きについての情報
を含む。帯域通過フィルタを使用することによって抽出され、これらのピークの
1つの逆フーリエ変換は、その位相が実際の動きの方向成分に線形に関連する複
雑な画像である。この位相画像は、前に本明細書で定義した「角度画像」である
。この位相画像は、(標準逆アークタンジェント演算子の働きによって)[−π
、π)の範囲にあるよう制限される。これは、角度ラッピングアーティファクト
(angle−wrapping artifact)である。角度ラッピン
グアーティファクトは存在するが、角度画像は、タグパターンを統合するために
使用され得、そして角度画像のペアは、小さな変位フィールド、調整を必要とし
ない画像ペア間のオプティカルフロー、ならびに2次元および3次元のひずみを
測定するために使用され得る。
【0035】 心臓は、2次元または3次元の1−1SPAMMタグパターンを使用して拡張
末期において、繰り返しタグ付される。タグ付パルスシーケンスは、左心室が血
液で一杯であり、心臓が相対的にゆっくり動き、そしてECG信号のQRS群が
存在する心臓周期の一部である拡張末期で課される。推定の目的のために、拡張
末期は、t=0なる時刻とみなされ得、拡張末期における心臓内の点の位置は、
物質座標系として扱われ得る。連続した心臓シーケンス中、k空間は、フーリエ
空間における9の主要なスペクトルローブのうち少なくとも1つを含むフーリエ
変換情報を獲得するために、標準定常状態傾斜エコー画像化パルスシーケンスを
使用してスキャンされる。2つ以上のスペクトルピークが画像化される場合、帯
域通過フィルタは、1つのピークの情報だけを抽出するために適用される。この
データの逆フーリエ変換は、共役対称操作を行うことなく、行われる。共役対称
操作を行うことは、MR画像化においては通例である。その結果生じる複雑な画
像の角度は、角度画像を形成する。
【0036】 傾斜エコー磁気共鳴画像化によってスキャニングを使用する代わりに、スピン
エコー、らせん磁気またはエコー面磁気共鳴画像化などの代わりの公知の手段が
、例えば使用され得ることが理解される。
【0037】 1次元の1−1SPAMMタグパターンは、適用された横断勾配パルス(tr
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
に示されている。
【0038】 これらのスペクトルピークの存在は、搬送高調波を提供するタグプロセスのコ
ンテキスト内で理解され得る。搬送高調波は、画像を空間的に振幅変調させ、そ
れゆえに、対応するスペクトルピークを搬送高調波の位置へと移動させる。
【0039】 数学的コンテキストでこれを置き換えると、時間tに撮られたタグ付けされた
MR画像は、y1が出力方向で、且つy2が位相エンコーディング方向である画像
平面において、任意の点、y=[y12Tにて、輝度値を与える、Ψ(y,t
)により表され得る。スペクトルピークの存在ゆえに、画像Ψ(y,t)は和と
して書かれ得る。
【0040】
【数1】 ここで、各画像Ψk(y,t)は、スペクトルピークに相当する画像である。整
数kは、スペクトルピークを表すIDである。スペクトルピークの位置は、SP
AMMパルスシーケンスによって決定され得るベクトルwk=[w1k2k3 kTによって決定される。スペクトルピークの合計数は2K+1である。この値
は、SPAMMパルスシーケンスの数および特性に依存する。
【0041】 画像Ψkは、複合画像であり、すなわち、マグニチュード(Dk)および位相(
φk)を有する。故に、
【0042】
【数2】 である。タグ分類などの適切な条件下において、Ψkは、帯域通過フィルターを
使用して、Ψから取り出され得る。
【0043】 その角度画像は (3) ak(y,t)≡∠Ψk(y,t), k=−K,...,K を使用して、複合画像から計算され得る。ここで、
【0044】
【数3】 である。ここで、
【0045】
【数4】 は複合画像の想画像構成要素であり、
【0046】
【数5】 は複合画像の実在の要素である。
【0047】 例えば、図2で丸で囲まれたスペクトルピークから計算された角度画像が、図
3に示されている。この画像の鋸歯パターンは、範囲[−π,π)に包み隠され
る複合数量の角度の結果として生じる。
【0048】 角度画像は、いくつかの非常に有用な以後の分析の基礎となり得る。第1に、
画像はタグ付けされたMRI画像における通常の2次元タグに類似した総合タグ
を生成するために使用され得る。この特徴の利点は、データが完全に自動的に生
成され、且ついかなる所望のタグ分類を用いても生成され得ることである。第2
に、画像は、オブジェクトの小さな変位を直接計算するために使用され得る。第
3に、角度画像は、2次元画像面における2次元のゆがみ、または3次元におけ
る完全な3次元ゆがみテンソルを直接計算するために使用され得る。これらのゆ
がみデータは、心筋虚血および心筋梗塞の検出およびそれらを量で表すときに有
用である。第4に、角度画像は、時系列の変位フィールドを表す標準的なオプテ
ィカルフローフィールドを形成するために使用され得る。
【0049】 図4および5は、小さな変位のための、計算された変位および実際の変位をそ
れぞれ示す。計算されたフィールドおよび真の変位フィールドの類似性は、角度
画像の動きが、動きを再構築するに効果的であることを示している。
【0050】 画像の統合タグ付の線にとって、シングル角度画像、ak(y,t)が採用さ
れ得る。タグ付の線は、 (5) ak(y*,t)=a を満たす点{y*}の収束である。
【0051】 図6は、図3の角度画像から生成されるタグ付の線を示す。角度画像を使用し
て総合タグ付の線を生成することには、いくつかの利点がある。第1に、完全な
自動プロセスであることである。第2に、よい等輪郭(isocontour)
アルゴリズムがこの属性を有しているように、タグ付の線がサブピクセル解画像
度を有することである。第3に、全画像が、例えば、左右双方の心室心筋を含む
など、自動的に識別されるこれらのタグを有することである。最後に、範囲[−
π,π)内のN値を選択することによって、Nタグ付の線が、空間的期間2π/
k以上の範囲において、統合されることである。原則的には、これらのタグ付
の線がどのぐらい接近して間隔を開けられ得るかを制限するものは基本的にない
。これは、これらのタグ付の線が互いにどのぐらい接近しているかのという特徴
に対する検出可能性によって制限されないからである。
【0052】 小さな変形では、画像作成面に存在する2本の直線の独立したベクトルwk
よびwlの2つの角度画像(akおよびal)は、
【0053】
【数6】 を使用してtにおいて画像面の変位フィールド(u2)の投影を計算するために
使用され得る。ここで、Wは行列であり、Tは行列の転位であり、Δakは、
【0054】
【数7】 によって計算される。 θkは、パルスシーケンスから決定される既知の角度であり、Wは (8) W(φ)=mod(φ+π,2π)−π によって与えられる非線型ラッピング関数(wrapping functio
n)である。ここで関数x(y)は、 (9) x(y)=y11+y22+x0=Hy+x0 を使用して、画像平面座標システムの任意の点yを、磁気3次元座標システムの
位置x ε IR3にマッピングさせる。ここで、行列H ε IR3×2
[h22]であり、ベクトルh1 ε IR3および、ベクトルh2 ε IR3
、それぞれ画像平面の読み出しおよび位相エンコーディング方向を表し、x0
ε IR3は、画像原点である。行列W ε IR3×2= [w12]である。
【0055】 モニターされている動きの大きさを最小化する方法はたくさんあり、これらは
比較的小さな動きに対し実行される本発明の方法を容易にする。例えば、1つの
方法は、心臓が実質的に動く機会を得る前に、心弛緩終期後、即画像を作成する
ことである。このアプローチは、収縮期の最初の数10ミリ秒において、有用で
あり、且つ潜在的および臨床的に重要である。第2に、低周波タグパターンが使
用される場合、タグパターンの物理期間がだんだん大きくなり、動きは角度のあ
いまい性(ラッピング)を生じない。これを選択した際の潜在的困難性は、低周
波パターンのスペクトルピークが別のスペクトルピークと干渉し、所望されない
動きのアーティファクトが生じることである。第3のアプローチは、心弛緩終期
およびそのすぐ後の画像から、固定の補正値でタグパターンを適用することであ
る。この場合、タグパターンを急速に動く組織へ適用することが必要である。最
後に、この問題を回避する方法は、タグ付け時間からの変位ではなく、画像時間
の間の変位を画像作成することである。これらの変位は小さい。このアプローチ
は、古くは「オプティカルフロー(optical flow)」と呼ばれる。
【0056】 ゆがみは組織に隣接する部分間の変位差に関する。これは、2本の直線的で独
立したベクトルwkおよびwlの、少なくとも2つの角度画像である(akおよび
l)から直接計算され得る。e方向での2次元のゆがみは、
【0057】
【数8】 によって計算される。ここで
【0058】
【数9】 であり、且つ
【0059】
【数10】 であり、且つΔy* lと同様の方程式である。最後の2つの方程式は、角度画像
の派生を計算すると同時に、ラッピングの不連続性を克服するために使用される
。これらの方程式から計算されたゆがみは、オイラー形式(Eulerian
sense)である。完全なゆがみテンソルは、3つのスペクトルピークから来
る3つの角度画像から計算され得る。3つのスペクトルピークの生成は、3次元 SPAMMパターンを使用し、画像平面ではなく画像体積を獲得することによ
って行われる。
【0060】 オプティカルフローは、画像シーケンスにおける輝度パターンの明らかな動き
として定義される。概略は、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年を参照されたい。本発明のコンテキスト内における角
度画像の使用により、規則化を使用することなしに、速度フィールドを直接計算
することが可能である。
【0061】 角度画像をオプティカルフロー方法に適用するために、直線で且つ独立したベ
クトルwkおよびwlを有した、少なくとも4つの角度画像ak(y,t)、al
y,t)、ak(y,t+Δt)、およびal(y,t+Δt)が採用され得る。
2つの画像間の時間間隔、Δt>0は、時間間隔における一定の動きの速度v2
を仮定するのに十分小さいことが好ましい。2次元の動きの速度は
【0062】
【数11】 によって計算され、ここで (14) Δtk(y)≡W[ak(y,t+Δt)−ak(y,t)] である。
【0063】 図7aおよび7bを参照すると、図7aは1−1SPAMMタグを備えた左心
室を示す断面図である。図7bは、本発明の方法により作成される図7aの左心
室の断面図の角度画像を示す。
【0064】 図8から18および関連する開示に例示されるような本発明は、タグ付けされ
た心臓の磁気共鳴画像シーケンスの迅速な分析のための画像処理方法を提供する
。この方法は、心臓の動きに関する情報を含むSPAMMタグの磁気共鳴画像に
おける分離されたスペクトルピークの使用を含む。SPAMMスペクトルピーク
の逆のフーリエ変換は、複合画像であり、その画像の計算された角度は高調波位
相(HARP)画像と呼ばれる。この方法は、2つのHARP画像シーケンスを
使用し、自動的且つ正確に時間的に物質点を追跡する。1実施形態において、迅
速な、半オートメーション化されたプロシージャは、これらの痕跡を使用し、周
囲ゆがみおよび放射状ゆがみ双方を含むラグランジュのゆがみを計算する。これ
らの方法は、2次元のコンテキスト内で開発されたが、3次元で使用され得る。
この新しい計算アプローチにより、スキャン完了後約5〜10分内の迅速な分析
および可視化が可能となる。その性能は、MR画像シーケンスに示され得、正常
な心臓の動きおよび異常な心臓の動きの双方を反映する。
【0065】 タグ付けされた心臓の磁気共鳴画像における過去10年間の主な開発(Zer
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年を
参照されたい。
【0066】 タグ付けされたMRIは、最近の多くの医学研究および科学調査において、著
しく目立ってきた。タグ付けされた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月)。
【0067】 概して、タグ付けされたMR画像の処理および分析は、3つの段階へと分けら
れ得る。その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年)がある。
【0068】 タグ付けされたMR画像処理の第3ステージである、ゆがみ推定は、主に補間
および分化計算であり、有限要素法を含む、文献に記載のいくつかの方法(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処理方法は、完全に自動化されておらず、且つ
緻密なゆがみ推定を形成するために、補間を必要とする。本発明はこれらの関心
の双方を取り扱う。
【0069】 本発明は、高調波位相(HAPR)画像作成と呼ばれる、タグ付けされたMR
画像の分析への新しいアプローチを記載してきた(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画像を採用し得
、画像シーケンスは、好適には、物質点追跡技術およびこれらの追跡された点を
使用して、周辺ゆがみおよび放射状ゆがみを含むラグランジュゆがみを計算する
方法の双方を含む。
【0070】 本発明において提案される方法は、速く且つ完全に自動化されており、迅速に
収集され得るデータを使用する。これらの方法は、2次元画像に直接のみ適用す
る。この結果、推定された動きの量は、3次元の動きから2次元平面への投影を
表し、推定された動きの量は「明らかな」動きとして考えられるべきである。本
方法は、3次元画像へと拡張され得るが、緻密なデータ獲得方法が採用される。
本方法は、即座の臨床効果を有すべきである。これは、その方法が、自動であり
、そして周辺ゆがみが、左心室の動きの分析において特に重要であるからである
【0071】 第1のSPAMMタグ付け方法を考察すると、I(y,t)で、画像座標 y
=[y12T、且つ時間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タグパルスシーケンスによって決定される基本タグ振動数の整数倍である
【0072】 スペクトルピークの2次元パターンは、2次元のSPAMMタグパターンを使
用して生成され得る。本発明の方法はこの場合においても利用され得る。
【0073】 静的なHARP画像法(Osmanらの発表された「Imaging hea
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
において、最適な実施のためのフィルタの選択が、議論される。帯域通過領域の
逆フーリエ変換により、
【0074】
【数12】 により与えられる複雑な高調波画像が生成される。ここで、Dkはマグニチュー
ド画像(the magnitude image)と呼ばれ、φkは位相画像
と呼ばれる。式15におけるIkの使用は、式10のΨkと同じである。
【0075】 マグニチュード画像Dk(y,t)は、心臓のジオメトリの変化とタグフェー
ジングによる画像強度の変化の両方を反映する。図8(c)は、図8(b)のフ
ィルタを使用して図8(a)から抽出された高調波マグニチュード画像を示す。
フィルタリングプロセスによる不鮮明を除いて、その基本画像のように基本的に
は見える。高調波マグニチュード画像にタグパターンがないために、高調波マグ
ニチュードオ画像を使用して、組織を背景と区別するための断片を提供し得る。
単純な閾値を使用して、未加工の断片を提供し得る。ここで、その閾値は、拡張
末期と収縮末期の両方の時点で手動で選択され、これらの時間の間で線形補間さ
れる。
【0076】 位相画像φ(y,t)は、ωkの方向の心筋の動きの詳細図を提供する。原理
的に、Ikの位相は、実部で虚部を割ったアークタンジェントを使用することに
よって計算され得る。Ikの符号を説明すると、この計算の固有の範囲は、例え
ば、C、Fortran、またはMATLABの演算子atan2を使用して、
[−π,−π)に拡張され得る。しかし、これによると、一般的に実際の直線全
体の値を示す実際の位相でなく主値のみが生成される。この主値は、ak(y,
t)によって示され得る;それは、
【0077】
【数13】 によって数学的にIkの真の位相に関連する。非線形ラッピング関数は、
【0078】
【数14】 によって与えられる。
【0079】 akまたはφkのどちらかが高調波位相(HARP)画像と呼ばれ得る。この表
現は、φkでなくakに一般に使用され、データからそれを直接計算し、視覚化し
得る。その2つは混同されるかもしれないが、φkを高調波位相と呼び、akを高
調波位相角度と呼ぶ。図8(b)でアウトラインが示されるスペクトルピークに
対応するHARP角度画像が図8(d)に示される。明瞭さのために、それは、
図8(c)のマグニチュ−ド画像の生の断片を使用して作れられたマスクで表示
される。
【0080】 図8(d)のHARP画像を注意深く調べると、ラッピングされたアーティフ
ァクトによる鋭い転位によって垂直方向に妨げられた強度傾斜がわかる。これら
の転位の位置は、図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トラッキングの基本である。
【0081】 タグ付のパルスシーケンスについての「見かけ上の」動きの問題に移ると、タ
グ付の勾配は、一般に画像平面において適用される。この場合、拡張末期の画像
に生じるタグ付の線は、図9(a)に示されるように、実際に画像平面に直交す
るタグ面の一部である。高調波位相は材料の特性であるので、拡張末期に同じ高
調波位相φを有する一組の点はまた、画像平面に直交する平面であり、単に別の
タイプのタグ面であると考えられ得る。拡張末期で同じHARP角度aを有する
一組の点は、ただ1つの平面ではなく平行平面の集合を含む。これは、次のセク
ションで説明するHARPトラッキングにおいて問題を引き起こす。見かけ上の
動きを説明するために、高調波位相φによって与えられる固有の関係が考えられ
る。
【0082】 水平なタグ画像および垂直なタグ画像の双方を与えられると、拡張末期の同じ
2つの高調波位相を有する一組の点は、単一の点で交差する画像平面に直交する
線を含むことが図9(a)から明らかである。図9(b)に示されるように、タ
グ面が動き下でひずみ、それによってこの線を曲線に歪ませる。動きについての
緩い仮定の下で、この曲線はまだ、1つの点で画像を交差する。次いで、この点
は、拡張末期で対応する点と関連し、画像平面内の見かけ上の動きを表示し得る 見かけ上の動きが、q(y,t)によって示される見かけ上の参照マップを使
用して数学的に説明され得る。この関数は、点を(同じ2つの高調波位相を有す
るという意味において)yが見かけ上の拡張末期にあった画像平面内に与える。
q(y,t)が、拡張末期の真の3次元の質点の位置を画像平面に直交に投影し
たものであることが図9から明らかに示され得る。見かけ上の2次元の動きを計
算するには制限があるが、その計算は、真の3次元の動きと非常によく似た関係
を有する。ひずみ等の見かけ上の動きから導出された動きの量は、等しく厳密な
方法による真の3次元量に関係し得る。
【0083】 CINEタグ付の磁気共鳴画像のHARP処理は、(1)画像平面の質点の見
かけ上の動きのトラッキングと、(2)そのようなトラッキングされた点からの
ラグランジアンひずみの計算を含む。コンパクトな式にするために、ベクトル表
示を使用する。詳細には、ベクトルφ=[φ1 φ2T、およびa=[a12
Tが定義され、それぞれ高調波画像I1およびI2についての一対の高調波位相
画像およびHARP角度画像を説明する。
【0084】 HARPトラッキングにおいて、高調波位相は直接計算できないので、その主
値であるHARP角度が、計算で使用される。多くの点が同じ一対のHARP角
度を有する画像平面にあることが即座に結論される。2つのHARP角度を有す
る所与の質点に対して、同じHARP角度を共有する点の内の1つのみが正しく
マッチングしたものである。すなわち、それはまた、同じ一対の高調波位相を共
有する。見かけ上の動きが1つの画像から次の画像まで小さい場合、最近接のこ
れらの点が正しい点であることは明らかである。本発明の方法は、タグ付のMR
画像のCINEシーケンスを介して見かけ上の動きをトラッキングする。位置が
mで、時間がtmの質点を仮定する。ym+1が、時間tm+1でのこの点の見かけ上
の位置である場合、
【0085】
【数15】 である。この関係は、時間tmからtm+1までymのをトラッキングするための基
準を提供する。
【0086】
【数16】 を満たし、ym+1=yと設定したyと見出すのが望ましい。式(19)は、多次
元且つ非線形で、ルートファインディング問題(root finding p
roblem)である。この問題は、ニュートンラプソン法を使用して反復的に
解かれ得る。単純化の後、ニュートンラプソン反復式(iteration)は
【0087】
【数17】 である。
【0088】 式(20)の直接的な用法にはいくつかの実際的な問題がある。第1の問題は
φが利用できないことであり、その代わりにaを使用しなければならない。幸い
、φを含む式をaを含む式で置き換えるのは比較的簡単である。akの勾配は、
勾配の大きさが理論的に無限−−実際には非常に大きい(図8(d)参照)、ラ
ッピングアーティファクト以外はΦkの勾配と同じであることは式(16)から
明らかである。πをakに加え、再度ラッピングすることによって、1/2の空
間距離だけラッピングアーティファクトをシフトさせ、それによってこの勾配が
結果的にΦkの勾配に等しくなることは、式(16)から明らかである。結果と
して、Φkの勾配は、akおよびW(ak+π)の勾配の(大きさにおいて)小さ
いものに等しい。形式的に、これは、
【0089】
【数18】 として書くことができる。ここで、
【0090】
【数19】 および
【0091】
【数20】 である。
【0092】 式(20)の使用について第2の問題は、式φ(y(n),tm+1)−φ(ym
m)の計算であり、この計算は、それのラッピングされたバージョンである高
調波位相だけでなく位相自体も分からないのでほとんど不可能と思われる。しか
し、k=1,2、|φk(y(n),tm+1)−φk(ym,tm)|<π(微小な動き
を仮定)に関して、
【0093】
【数21】 と示すことが簡単である。式(21)および式(24)を式(20)に代入する
と、フォーマルな解である
【0094】
【数22】 が導かれる。
【0095】 式(25)に基づいてアルゴリズムを実行する際に考慮するいくつかの問題が
ある。第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)反復カウントが超えられるか。
【0096】 これらをまとめて考慮すると、1つの時間フレームに1つの点を次の時間フレ
ームのその見かけ上の点までトラッキングするアルゴリズムが即座に定義され得
る。画像シーケンス全体を介して点をトラッキングするHARPトラッキングア
ルゴリズムを定義するのを容易にするために、より一般的なフレームワークでこ
のアルゴリズムを変換することは有用である。従って、その方法では、yinit
、検索が開始される初期化であると考え、a*が目標HARPベクトルであると
考える。時間tm+1におけるその見かけ上の位置まで時間tmの点ymをトラッキ
ングするために、yinit=ym,a*=a(ym,tm)と設定し、最大反復カウン
トNを求め、次いで以下のアルゴリズムを実行する: (アルゴリズム1(HARP目標設定)n=0とし、y(0)=yinitと設定す
る) 1.
【0097】
【数23】 である場合、アルゴリズムは
【0098】
【数24】 で終了する。 2.適切な補間手順を使用してステップの方向
【0099】
【数25】 を計算する。 3.ステップサイズ
【0100】
【数26】 を計算する。 4.推定値である y(n+1)=y(n)+α(n)(n) を更新する。 5.増分nを増分し、工程1に行く。
【0101】 画像シーケンスを介して点をトラッキングするために、HARP目標設定がシ
ーケンスの各画像にうまく適用される。正しく見かけ上の動きを見出す好適なア
プローチは、全体のシーケンスを介して目標HARPベクトルを同じに、すなわ
ち元の点のHARPベクトルと等しいように維持することであり、シーケンスに
おいて以前に推定された見かけ上の位置にHARP目標設定を初期化することで
ある。これにより、同じHARP角度を有する連続した点を生成し、一般に所望
の解に近いHARP目標設定で使用された初期の点を維持することによって間違
った解に向かうのを避ける。形式的にこのアルゴリズムを説明するために、時間
j+1,tm+2・・・・・の全ての画像を通じて時間tjのyjをトラッキングする
ことを所望すると仮定する。
【0102】 HARPトラッキングは、以下のアルゴリズムにより与えられる。
【0103】 (アルゴリズム2(HARPトラッキング)a*=a(yj,tj),
【0104】
【数27】 ,およびm=jと設定する。) (HARP目標設定に対する)最大反復閾値Nを選択する。 1.
【0105】
【数28】 と設定する。 2.アルゴリズム1(HARP目標設定)を適用し、それによって
【0106】
【数29】 を生成する。 3.mを増分し、工程1に進む。 HARPトラッキングは使用され、それにより時間の順方向と厳密に同じ方法で
時間の逆方向で点をトラッキングし得る。それゆえ、任意の時間の任意の画像の
点を特定し、順方向と逆方向の両方に点をトラッキングし、それによって空間と
時間において任意の点の完全なトラジェクトリを与えることが可能になる。
【0107】 ラグランジアンひずみを判定する際に、参照時間t=0が拡張末期であると仮
定する。HARPトラッキングは、(利用可能な画像倍の集合に対して)時間t
での質点の位置yに対するt=0での質点qをトラッキングし(このセクション
において、全ての位置の量を「見かけ上の」位置と呼ぶ)得る。これにより動き
マップy(q,t)の推定値を提供する。y(q,t)を使用して、有限の違い
を使用して任意の質点qおよび時間tの変形勾配テンソル
【0108】
【数30】 を推定し得る。変形勾配を計算することによって、対象の他の動きの量は容易に
計算される(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トラッキングのさらに強力な適用は、準自動分析内のいくらか
単純な以下の計算で示される。
【0109】 2つの質点qiおよびqjの動きを考える。1単位の伸びまたは単純なひずみが
【0110】
【数31】 により与えられる。点間の距離が変化しない場合、この量はゼロであり、点間の
距離が縮むと負になり、点間の距離が伸びると正になる。左心室壁内の任意の位
置の円形ひずみを測定するために、LV長軸の中心である円に沿って2つの点を
単純に配置する。動径方向ひずみを測定するために、長軸から生じる光線に沿っ
て2つの点を単純に配置する。どちらかの場合に、ひずみは、HARPトラッキ
ングを使用して、2つの点をトラッキングし、eを計算することによって測定さ
れる。HARPトラッキングが基本的に画像内の任意の位置をトラッキングでき
るので、これらの点の位置が「タグ付の線の交差」または画素の位置にさえある
必要はないことを再度強調しなければならない。ひずみのこの測定は、ラグラン
ジアン(またはユークリッドひずみ)テンソルの緻密な計算に対して2つの利点
を有する。第1に、画像全体(または対象の領域)の代わりに2つの点のみのト
ラッキングが必要があるので、それは非常に速い。第2に、点は、1つの画素に
比べ一般にかなり離れて配置されるので、伸びについての計算は、元来ノイズに
敏感ではない。
【0111】 これらの利点を利用するために、心筋内の同心円上の点をトラッキングし、領
域性の動径方向および円形ラグランジアンひずみ(circumferenti
al Lagrangian strain)を計算するアプローチが実行され
た。単純なユーザインターフェースは、図10に示すようにLV壁内の3つの同
心円の配置を可能にする。これらの円は、第1に心室の中心(長軸の位置)上で
クリックし、次いで心外膜に対する1つの円および心内膜に対する別の円をドラ
ッグすることによって手動で配置される。3番目の円は、これら2つの円の間の
中間に自動的に配置される。通常、円は拡張末期の画像で定義されるが、LVの
断面がより円形であり得るので、収縮末期の画像を使用してそれらを定義するの
がより容易である。16個の等間隔の点が各円の外周の周りに自動的に定義され
、アルゴリズム2(HARPトラッキング)を使用して48点全てが(時間の順
方向または逆方向に)トラッキングされる。
【0112】 ひずみは、時間フレームを通じて隣接する点間の距離が変化するのを測定する
ことによって計算される。円は、拡張末期に、収縮末期に、またはその間の任意
の時点に定義されるかどうかに関わらず、拡張末期の画像フレームは、物質標準
として使用される。同じ円にある点間の距離の変化が円形ひずみに対応する。動
径方向の点間の距離の変化は動径方向ひずみに対応する。3つの同心円状がある
ので、心内膜、心外膜、および中間壁の円形周囲ひずみおよび心内膜および心外
膜の動径方向ひずみが計算され得る。ノイズを幾分か低減し、それによって提示
を単純化するために、円は、8個の八分円に分けられ、計算されたひずみは、こ
れらの各八分円で平均化された。慣習により、八分円は、図10に示されるよう
に隔壁の中心から時計方向に番号付けされた。その結果、生じるひずみは、時間
と八分円の両方の関数としてプロットされ、それに断面図内の心臓の動きを空間
時間的に表示する。
【0113】 上記の手順は、点の集合をトラッキングする工程を含む。その工程は、各点に
アルゴリズム2を独立して適用することによってうまく達成され得る。しかし、
場合によっては心筋の大きな動きまたは画像アーティファクトにより、ある時間
フレームで点が誤った目標に収束する(タグジャンプ)ことが起こり得、同様に
、連続したフレームで誤ったトラッキングが生じることになる。本発明は、1つ
以上の正確にトラッキングされた点を使用する改良手順を使用し、それによって
誤ってトラッキングされた点のトラッキングを訂正する。
【0114】 上述したように、アルゴリズム2は、トラッキングされた点の以前に推定され
た位置で初期化される。しかし、2つの時間フレームの間の平面内の動きが非常
に大きい場合、この初期点は正確な解から非常に離れ、誤った点に収束し得る。
改良は、アルゴリズム2に対するよりよい初期化についての体系だった識別に基
づく。所与の円にある1つの点は、全てのフレームを通して正確にトラッキング
されてきたと証明され得ると仮定する。発明者らの実験において、動きが相対的
に小さい隔壁のそのような点が常に見出されてきた。
【0115】 この点は、円にある他の全ての点の初期化および、所望ならば3つの円にある
全ての点の初期化を改善し得、全体にわたって集約されたトラッキング結果を改
良し得る「アンカー」として使用される。
【0116】 正確にトラッキングされたアンカーから始めることで、円に隣接する点の1つ
とアンカーをつなぐ線セグメントにある1画素未満に分けられた点のシーケンス
が定義される。アンカーは、いくつかの特定の時間に点yまでトラッキングされ
ると仮定する。この点が正確にトラッキングされると仮定することによって、そ
のアンカーに近い点の正確なトラッキングの結果がyに近くなる。従って、yは
、アルゴリズム2の初期点として使用され、シーケンスの最初の点をトラッキン
グする。次いで、この結果は、シーケンス上の第2の点等に対する初期点として
使用される。円上のアンカーの第1の隣接する点に達することによって、誤った
結果に収束し、タグをジャンプさせる機会はない。隣接する点は、新しいアンカ
ーとして役に立ち、その手順は、円上の全ての点がトラッキングされるまで円上
の次に隣接する点に対して繰り返される。
【0117】 改良は、動径方向の経路に沿って円の「橋渡し」をし、3つ全ての円をうまく
補正するために使用され得る。一般に、タグジャンプ誤りは、自由壁でのみ起こ
り、改良が計算的に要求されるにしたがって、その演算は一般に一定の時間では
1つの円に抑制される。所望ならば、チェックして、円は、元のアンカーの周囲
をトラッキングすることによって、完全なものになり得る。その結果が異なる場
合、総計としての誤りがあり、円を再定義することが望ましい。タグフェージン
グまたは他の画像アーティファクトによりこのタイプの総計としての誤りが生じ
得るが、平面外の動きによってこのタイプの総計としての誤りがさらに生じ易く
なる。平面外の動きにより、画像化される実際の組織は変化する。この場合、ト
ラッキングされた点に対応する高調波位相角度をもたらす画像平面上に組織が存
在することはあり得ない。要約すると、タグが消え、HARPトラッキングアル
ゴリズムは、正確なHARP角度を有する別の点に簡単に収束する。左心室の動
きとジオメトリとの間の特定の関係のために、これは大した問題ではない。問題
は、左心室の境界近くでより生じ易いので、それにより課される主な制限は、心
外膜または心内膜の非常に近く円を配置してはならないことである。
【0118】 (実施例) 本発明の方法の効率性を確認するために、正常の心臓と異常な心臓の両方につ
いての実験が、犬と人間の両方を含む試験について行われた。異常に拍動する犬
のデータセットが、ドブタミンの応力下にある正常な人間の心臓からのデータと
同様に得られた。HARPトラッキングと従来技術の周知のタグ付トラッキング
法(tag tracking method)の間の質的な比較が為された。
【0119】 (実施例1) HARPトラッキングの性能を示すために、電気的に拍動させた犬の心臓につ
いてのタグ付画像のセットが使用された。これらのデータは、タグ付磁気共鳴画
像法および分析技術を使用した、電気的に拍動させる活性化下での心臓の動きの
研究において、以前に使用された。実験手順およびその結果についての十分な説
明が、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)で生成される。これは、本発明の
方法を使用することにより結果が正確かつ短時間で生成されるという利点を示し
てる。
【0120】 拍動リード(pacing lead)が、犬の心臓の左心室の基底自由壁(
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スライスの厚さ。
【0121】 1つは水平タグであり、残りは垂直タグである20タグ付MR短軸画像の2つ
のシーケンスは、心臓が収縮してる間に14ms間隔で獲得された。画像は、断
片化されたk空間獲得(k−space acquisition)により息を
止めてる間に獲得された。拍動リードの位置に近い基底面(basal pla
ne)において獲得された画像が使用された。図11は、視覚化のためにLVの
周囲の関心ある領域でクロップされた結果として生じる画像を示す。強い初期の
収縮は、約5時(5 o’clock)に拍動リードの近くで見られ得る。フレ
ーム4〜8において、中隔壁(septal wall)が外向きに曲がってる
のが見られ、それは中隔領域への電気的な活性化信号の遅延により引き起こされ
るプリストレッチング(pre−streching)と呼ばれる異常な動きで
ある。この後、LV心筋全体が、収縮期を通じてほとんど正常に連続的に収縮す
ることを経験する。
【0122】 HARPトラッキングにおいて、HARP画像は、図8(b)に示される帯域
通過フィルタを使用して水平タグ付画像シーケンスから計算された。このフィル
タを90°回転させたバージョンが、垂直HARP画像を計算するために使用さ
れた。HARPトラッキング(アルゴリズム2)を示すために、1つは垂直であ
り、残りは水平である2つのタグ付の線が選択され、それぞれのタグ点の収集が
手動で選択された。これらの点の位置、および後に3回にわたる時点で個々にト
ラッキングされたこれらの位置が、図12に示される。1つだけ例外はあるが、
全ての点は、タグを消去した後にさえ点を見るのを期待する箇所にトラッキング
された。詳細には、正常な収縮の内部「反り(bowing)」および異常なプ
リストレッチングの外部「反り」が、HARPトラッキングによって非常によく
取り込まれる。間違ってトラッキングされた点のみが、図12(d)の画像の上
部に見られ得る。画像を注意深く調べることによって、平面の外での動きにより
、第1の時間のフレームのLVの上部にあるタグ付の水平線が長い間に消えるこ
とがわかる。改良によって解決されない間、このタイプの問題は、心筋の境界に
それほど近くない点を選択することによって避けられ得る。
【0123】 HARPトラッキングがタグ付の線の点に限定されないことおよび圧力テンソ
ルを計算するためのポテンシャルを示すために、1画素(1.25mm)によっ
て分けられる点の5×5グリッドが、図13(a)に示されるように、LVの前
横方向側の4つのタグ付の線によって境界を付けられた領域に配置された。これ
らの点は、全ての画像シーケンスを介してそれぞれにトラッキングされた。図1
3(b)〜(c)は、時間フレーム1,5,10,20でそれらの位置を示す拡
大図を示す。トラッキングされた点に対する副画素の解像度がより後の画像にお
いて明瞭に示され、ひずみの基本的なローカルパターン(local patt
ern)を明瞭に見ることができる。四角形からひし形というグリッドの明確な
発展によって、正常な心臓が動きする際に動径方向に厚くなることおよび円形方
向に短縮することの両方が大変明瞭に示される。このデータからひずみテンソル
を計算するために差分が容易に使用され得ることはトラッキングされた点の規則
性から明らかである。これから、主なひずみ(principle strai
nts)の局所領域変化および方向を含む動きに関する様々な量が計算され得る
【0124】 次に、領域ラグランジアンひずみは、アルゴリズム2の後で説明された手順を
使用して計算された。ユーザインターフェースを使用して、心外膜および心内膜
の円(circles)が定義された。収縮末期にLVはほとんど円形なので、
図11の最後の画像が、これらの円を定義するように使用された。結果として生
じる3つの円および定義された八分円が、図14(a)に示される。それぞれの
円の16個の点が、拡張末期まで後方にトラッキングされ、結果として図14(
b)に示される形となった。歪んだ状態のシーケンス全体が、図14(c)に示
される。このシーケンスから、LVの断面積の形状がいくらか伸張し始めるが、
急速に円形になり、次いでほとんど動径方向に収縮するのがわかる。そのような
タグのジャンプが、1つ以上の時間フレームにおいて非常に歪んだ輪郭を描くの
で、図14に誤ってトラッキングされた点はないことを確認にするのは容易であ
る。
【0125】 ラグランジアンひずみプロフィールは、アルゴリズム2において以下で説明さ
れるように図14(c)に示されたトラッキングされた点から計算された。各八
分円における動径方向のひずみについての時間的発展が、図15(a)で示され
る。正の値は、心筋の厚みが増加するのを示し、負の値はその厚みが薄くなるの
を示す。初期の時点で心筋の厚みが増加することは八分円3〜6のみで明らかで
あり、八分円8、1および2は厚みが薄くなるのを示す。これは、初期に拍動リ
ード近くの心筋で生じる大きな収縮と反対の壁(wall)での心筋のプリスト
レッチングの両方の直接的な表現である。時間フレーム5〜10の間で、再収縮
して拡張末期で厚さが最も大きく増加するその前に、八分円5〜7の拍動リード
近くの心筋は弛緩する。心内膜の厚さの増加がより大きい八分円7を除いて、心
外膜と心内膜の動径方向のひずみの間にほとんど違いはない。
【0126】 図15(b)は、各八分円内の円形ひずみの時間的発展を示す。正の値は、円
形方向に伸張するのを示し、負の値は収縮を示す。これらのプロットは、動径方
向のひずみのプロフィールに示されるような同一の一般的な振る舞いを示す。八
分円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)。
【0127】 (実施例2) この実施例は、ドブタミン−誘導の(5μg/kg/min)心臓の応力を経
験してる正常な男性ボランティア(27歳)の心臓の動きのHARPトラッキン
グの使用を含む。この実施例は、ドブタミンの応力下のこの心臓の非常に急速な
動きおよび誤ってトラッキングされた点を補正するためのHARP改良の使用を
含む。
【0128】 その画像は、上述の犬の試験と同じ基礎的な画像手順を使用して同じ磁石で獲
得された。SPAMMタグが、7mmで分けられた画像平面に直交する飽和平面
(saturation planes)を作成するために拡張末期で生成され
た。垂直タグおよび水平タグを備える画像の2つのセットが、個別に息止めして
る際に獲得された。4つのスライスが獲得されるが、この実施例では中間壁の基
本的なスライスのみが使用される。スキャナの設定は以下のようである:視野3
6cm、タグ分離7mm、8mmのスライスの厚さ、TR=6.5ms、E=2
.3ms、15°チップ角度、256×160画像マトリックス、動画フレーム
当たりの5位相−符号化されたビュー。
【0129】 図16(a)は、LVの周囲の関心のある領域に関する、共に加算され、クロ
ップされた結果として生じる水平タグおよび垂直タグ付の画像を示す。これらの
画像に生じる輪郭は、HARPトラッキングを使用する時間内に、心外膜および
心内膜の円を第1の画像に手動で配置し、それらを順方向にトラッキングするこ
とによって生成された。第1および第2の時間フレームの間の大きな動きのため
に、前自由壁(anterior free−wall)のいくつかの点が、第
2のフレームにおいてトラッキングされなかった。基礎的なHARPトラッキン
グアプローチは、現在のフレームにおける初期化として以前にトラッキングされ
た点を使用するので、残りのフレームにおいて、この残りは補正されなかった。
隔壁内の手動で識別された3つのアンカー(各円に1つ)を使用してHARP改
良を適用する結果が、図16(b)に示される。見て調べると、改良された結果
により、各時間フレームで予想される箇所にトラッキングされた点が配置された
のが理解され得る。タグジャンプが排除される。
【0130】 改良されたトラッキングの結果を使用して、各四分円での円形ひずみの時間的
発展が計算された。その結果は、図16(c)に示される。これらのプロットは
、自由壁において最大かつ持続性のある短縮を伴うLV中のかなり画一的な短縮
を示す。これらはまた、正常な筋肉において一般に生じるような心内膜中に生じ
る最大の短縮を示す。
【0131】 公知の従来技術の方法と量的な比較をするために、FindTagsとして知
られる公知の技術の基礎である適切なフィルタリングと比較して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の精
度が、適合したフィルタリングアプローチより同じまたはより優れていることが
示される。
【0132】 次に、正常な人間という対象からのデータを使用してタグ付の線を見出すこと
に注意が向けられた。FindTagsは、垂直なタグ付の短軸のデータセット
からの27個の画像において輪郭(心外膜および心内膜)およびタグ付の線の両
方を推定するために使用された。その27個の画像は、LV内の3つの長軸方向
の位置からの9つの画像シーケンスを含む。HARP画像は、図8(b)に示さ
れたものと同様の第1の高調波および帯域通過フィルタを使用してこれら9つの
画像から生成された。理論によると、タグの底は、πラジアンの位相角度で配置
され、それによって、HARP画像のπ等輪郭(isocontours)が、
FindTagsによって識別されるタグ点に非常に近接することになると予測
される。図17(a)では、タグのかなりの消去が生じる心室の中間についての
7つの時間フレームで撮られた画像におけるこれらのHARP等輪郭に、Fin
dTagsからのタグ点が重ね合わされる。
【0133】 FindTagsを使用して推定されたHAPR等輪郭とタグ輪郭の間はわず
かな違いしかない。HARPは非常に小さく滑らかな結果を生成し、主な違いが
、自由壁(3時)上のHARPによる結果の周りの小さな揺らぎとなる。どれが
視覚的に満足に値するかを決定するのは難しい。FindTagsタグ点の収集
全体を使用して計算された平均HARP角度が、非常にπ(3桁の数(thre
e significant digits))に近く、それによってタグの底
がπというHARP角度を有するという理論的予測を証明する。ローカル位相非
ラッピング(unwrapping)は、HARP画像のラッピング(wrap
ping)アーティファクトの位置に正確にπが対応する場合、この平均角度を
計算するために使用された。
【0134】 精度についての量的な測定を実行するために、同じ画像における各タグ点と最
近接の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の誤り
とほとんど同じであると結論される。
【0135】 HARPトラッキングは同時に2つのHARP画像を使用し、それによってタ
グ付の線自体よりタグ付の線の交差の認識の方がより容易になる。結果として、
この実験では、HARPトラッキングは、FindTagsを使用するタグ付の
線の交差の推定と比較された。本明細書中で説明された拍動する犬の心臓のデー
タを使用して、FindTagsは、20の画像基底短軸画像シーケンス(20
−image basal short−axis image sequen
ce)における垂直と水平の両方の全てのタグ付の線の位置を計算するように使
用された。心内膜および心外膜の輪郭を使用して、またFindTagsの使用
を推定され、心筋内部に落ち込む(falling)タグ付の線の交差点が計算
された。HARPトラッキングは、次いで、ターゲットベクトルa=[ππ]T
ラジアンを求め、第1の時間フレームを含むこれら各々の点で実行された。
【0136】 FindTags交差点とHARPトラッキングされた点の間の平方二乗平均
距離が、図18に時間の関数として示される。これらの誤りは、以前の実験より
いくらか大きい。これは、1つの線の代わりに2つの線を求める場合、誤りが
【0137】
【数32】 だけ増加することを留意することによって、ほとんど説明される。他の小さな違
いは、異なる実験的なセットアップと画像手順によって説明され得る。長時間に
わたって距離が増加するという一般的な傾向はまた、タグが消えるにつれて信号
−ノイズ比が降下するという理由から予測される。
【0138】 図18に示されるプロットにおいて、興味ある特徴の1つは、時間フレーム4
〜7において「こぶ(hump)」が生じることである。HARPで使用される
帯域通過フィルタの設定(placement)により、これを説明することが
可能である。収縮期の間に正常な心臓な動きは収縮するので、基本的なタグ高調
波(the fundamental tag harmonic)のローカル
周波数は一般的に増加する。また、タグが消えるにつれて第1の高調波の周りの
スペクトルとの干渉が生じるので、エネルギー的にDCスペクトルピークが増加
する。それゆえ、わずかに高い周波数に帯域通過フィルタを設定したり、収縮に
関連するその高い周波数スペクトルを同時に取り込んだり、DCスペクトルピー
クからの干渉を避けるために通常は使用される。しかし、フレーム4〜7では、
図11および図15に見られ得るように、プリストレッチング現象が最もはっき
りしており、これによりタグのローカル周波数は増加ではなく減少する。わずか
に高い違いの原因およびおそらく誤りを増加させる原因となるスペクトルのこの
部分は、フィルタでカットオフされ得る。
【0139】 上述の実験では、HARPトラッキング方法およびラグランジアンひずみ計算
法(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秒である。さらに一定の工程の重要な効
率化が可能になる。
【0140】 画像データの編集(organization)および帯域通過フィルタの設
定(definition)もまた、HARP分析の処理時間全体を増加する。
標準のスキャンでは、これらの時間は無視できる程であり、画像シーケンスは、
自動的に構築され得、現行の帯域通過フィルタが使用され得る。特定のスキャン
または実験手順では、いくつかの付加的な時間が、HARP処理のためにデータ
を準備する際に費やされ得る。当業者は、一般的にこれらの付加的な工程を30
分未満で実行し得、便利なユーザインターフェースによりさらにこの時間が低減
する。HARPは、さらにわずかな有効化および最適化を伴う臨床的使用に有用
である。
【0141】 HARP画像を使用して、高密度の増分速度場(incremental v
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トラッキングを初期化すること、または視覚化動きのた
めの非常に高速で高密度の動き場を生成することあるいは高密度のメッシュ上の
他の動き量を計算することについて有用であり得る。
【0142】 本発明の方法では、HARP角度という材料の性質が、速く、正確で、ロバス
トな2次元のHARPトラッキング法を開発するために利用された。点は、左心
室の円形方向および動径方向のラグランジアンひずみを直接に計算するために座
標系においてトラッキングされた。拍動した犬の心臓についての実験は、異常な
動きをトラッキングし、以前に報告された分析と一致するひずみを計算するHA
RPの能力を示した。誤ってトラッキングされた点の補正についての改良技術を
発展させ、ドブタミンの応力を経験する正常な人間の心臓について示された。最
後に、予備的な誤り分析が行われ、FindTags、標準テンプレートマッチ
ング法と比較するとHARPトラッキング法は非常に利点を有することを示した
。HARPトラッキングおよびラグランジアンひずみ分析が、計算的に高速で、
さらに適切な有効化の後では臨床的な使用に用いられることが示された。
【0143】 開示の簡易性の目的のために、本発明の方法を医学的用途に使用することを特
に参考にしてきたが、本発明の方法は、それに制限されず、広範な多様な産業的
用途および他の用途において使用され得る。
【0144】 それゆえ、本発明の特定の実施形態は例示のために本明細書中で説明され、当
業者には、本発明から逸脱することなく詳細部の多数の変更が為され得ることが
明らかである。
【図面の簡単な説明】
【図1】 図1は、2次元の1−1SPAMMタグ付画像のシミュレーション図である。
【図2】 図2は、図1の画像のフーリエ変換の大きさを示す。
【図3】 図3は、複雑な画像の角度を示す。
【図4】 図4は、対象における点の計算された変位を示す。
【図5】 図5は、対象の点を実際の変位とともに示す。
【図6】 図6は、図3に示される角度画像から生成される計算された変位タグ線を示す
【図7a】 図7aは、1−1SPAMMタグを有する左心室の断面図である。
【図7b】 図7bは、図7aの角度画像である。
【図8a】 図8aは、水平SPAMMタグを有する磁気共鳴画像を示す。
【図8b】 図8bは、図8aの画像のフーリエ変換の大きさを示す。
【図8c】 図8cは、図8bの円内のスペクトルピークを抽出することによって生成され
る、大きさを示す複雑な画像を示す。
【図8d】 図8dは、図8cの位相を示す。
【図9a】 図9aは、拡張末期のタグ面を模式的に示す。
【図9b】 図9bは、動きの結果生じるタグ面のゆがみを示す。
【図10】 図10は、左心室壁および八分円の図上に重ねられた同心円の模式図である。
【図11】 図11は、グリッドを提供するために垂直および水平のタグ画像を増やすこと
によって生成される電気刺激で鼓動させた犬の心臓のタグ付磁気共鳴画像のシー
ケンスを示す。画像は、左上に示される拡張末期から右下に示される収縮末期ま
での電気刺激で鼓動させた犬の心臓の動きを示す20個の時間フレームである。
【図12a】 図12aは、時間を追って追跡された時間フレーム1における手動で選択され
た点を示す。
【図12b】 図12bは、時間を追って追跡された時間フレーム5における手動で選択され
た点を示す。
【図12c】 図12cは、時間を追って追跡された時間フレーム10における手動で選択さ
れた点を示す。
【図12d】 図12dは、時間を追って追跡された時間フレーム20における手動で選択さ
れた点を示す。
【図13a】 図13aは、密集して拾われた点の位置を示す拡張末期の心臓を示す。
【図13b】 図13bは、時間フレーム1における追跡物質点の拡大を示す。
【図13c】 図13cは、時間フレーム5における追跡物質点の拡大を示す。
【図13d】 図13dは、時間フレーム10における追跡物質点の拡大を示す。
【図13e】 図13eは、時間フレーム20における追跡物質点の拡大を示す。
【図14a】 図14aは、収縮末期における手動で定義される円を示す。
【図14b】 図14bは、拡張末期へ逆のぼって追跡した後の図14aの円のゆがんだ形を
示す。
【図14c】 図14cは、追跡された円の全シーケンスを示す。
【図15a】 図15aは、各八分円における心外膜の半径方向のひずみ(1点鎖線)および
心内膜の半径方向のひずみ(実線)の時間的進展を示す。
【図15b】 図15bは、各八分円における心外膜の円周方向のひずみ(1点鎖線)、心中
壁の円周方向のひずみ(破線)、および心内膜の円周方向のひずみ(実線)の時
間的進展を示す。
【図16a】 図16aは、HARP改良なしの初めの左上の表示から最後の右下の表示まで
、短軸画像および追跡された円を含むドブタミン誘因ストレスを受ける通常の人
間の心臓を示す。
【図16b】 図16bは、図16aと類似しているが、HARP改良適用後の結果を示す。
【図16c】 図16cは、改良点を使用して計算された円周方向のひずみの時間的進展を示
す。
【図17a】 図17aは、HARPπ輪郭(白いカーブ)が上書きされているFindTa
gs(黒点)からのタグ点を有するタグ付画像を有する通常の人間のデータを使
用してFindTagsと比較してのHARP精度を示す。
【図17b】 図17bは、FindTagsのタグ点とHARPイソ輪郭との間の平均距離
を示す。
【図18】 図18は、FindTagsとHARPとの間のタグ交差推定における平均差
のプロットを示す。
───────────────────────────────────────────────────── フロントページの続き (81)指定国 EP(AT,BE,CH,CY, DE,DK,ES,FI,FR,GB,GR,IE,I T,LU,MC,NL,PT,SE),OA(BF,BJ ,CF,CG,CI,CM,GA,GN,GW,ML, MR,NE,SN,TD,TG),AP(GH,GM,K E,LS,MW,SD,SL,SZ,TZ,UG,ZW ),EA(AM,AZ,BY,KG,KZ,MD,RU, TJ,TM),AE,AG,AL,AM,AT,AU, AZ,BA,BB,BG,BR,BY,CA,CH,C N,CR,CU,CZ,DE,DK,DM,DZ,EE ,ES,FI,GB,GD,GE,GH,GM,HR, HU,ID,IL,IN,IS,JP,KE,KG,K P,KR,KZ,LC,LK,LR,LS,LT,LU ,LV,MA,MD,MG,MK,MN,MW,MX, NO,NZ,PL,PT,RO,RU,SD,SE,S G,SI,SK,SL,TJ,TM,TR,TT,TZ ,UA,UG,US,UZ,VN,YU,ZA,ZW (72)発明者 オズマン, ナエル エフ. アメリカ合衆国 メリーランド 21209, ボルチモア, アパートメント ビー, グリーン メドウ パークウェイ 6220 Fターム(参考) 4C096 AA09 AB50 AC04 AD14 DC25

Claims (27)

    【特許請求の範囲】
  1. 【請求項1】 磁気共鳴画像法により対象物の動きを測定する方法であって
    、該方法は、 該対象物の関心のある領域を空間的に変調するためにパルスシーケンスを適用
    する工程と、 該空間的に変調された対象物のフーリエ領域から少なくとも1つの第1のスペ
    クトルピークを獲得する工程と、 該獲得された第1のスペクトルピークの逆フーリエ変換情報を計算する工程と
    、 該スペクトルピークのそれぞれから第1の高調波位相画像を計算する工程と、 第2のスペクトルピークのそれぞれから第2の高調波位相画像を作成するため
    に異なる時間に上記プロセスを繰り返す工程と、 該第1および第2の高調波位相画像からひずみを判定する工程と、 を包含する方法。
  2. 【請求項2】 心臓の動きをトラッキングするために前記方法を用いる工程
    を包含する、請求項1に記載の方法。
  3. 【請求項3】 円形ラグランジアンひずみを判定するために前記方法を用い
    る工程を包含する、請求項2に記載の方法。
  4. 【請求項4】 動径方向ラグランジアンひずみを判定するために前記方法を
    用いる工程を包含する、請求項2に記載の方法。
  5. 【請求項5】 2次元画像を提供するために前記方法を用いる工程を包含す
    る、請求項1に記載の方法。
  6. 【請求項6】 動いている人間の心臓に前記方法を用いる工程を包含する、
    請求項1に記載の方法。
  7. 【請求項7】 SPAMMパルスシーケンスを前記パルスシーケンスとして
    使用する工程を包含する、請求項1に記載の方法。
  8. 【請求項8】 3次元画像を提供するために前記方法を使用する工程を包含
    する、請求項1に記載の方法。
  9. 【請求項9】 心内膜の、心外膜のおよび中間壁のひずみの少なくとも1つ
    を判定する際に前記方法を使用する工程を包含する、請求項3に記載の方法。
  10. 【請求項10】 心内膜の動径方向の、および心外膜の動径方向のひずみの
    少なくとも1つを判定する際に前記方法を使用する工程を包含する、請求項4に
    記載の方法。
  11. 【請求項11】 心筋におけるひずみを判定するために前記方法を使用する
    工程を包含する、請求項2に記載の方法。
  12. 【請求項12】 2つの異なる時間の点間の前記間隔を判定することにより
    ひずみを判定するために、前記方法を使用する工程を包含する、請求項1に記載
    の方法。
  13. 【請求項13】 2つの異なる時間の点間の距離の増加を測定するために、
    前記方法を使用する工程を包含する、請求項1に記載の方法。
  14. 【請求項14】 2つの異なる時間の点間の距離の減少を測定するために、
    前記方法を使用する工程を包含する、請求項1に記載の方法。
  15. 【請求項15】 前記心臓の左心室内部の動きを測定するために前記方法を
    使用する工程を包含する、請求項6に記載の方法。
  16. 【請求項16】 約拡張末期で前記動きの測定を開始する工程を包含する、
    請求項6に記載の方法。
  17. 【請求項17】 水平方向と垂直方向の両方にタグ付画像を使用する工程を
    包含する、請求項1に記載の方法。
  18. 【請求項18】 画像平面の前記点の見かけ上の動きをトラッキングするこ
    とにより前記ひずみを判定する工程および該トラッキングされた点からのラグラ
    ンジアンひずみを判定する工程を包含する、請求項12に記載の方法。
  19. 【請求項19】 タグ付磁気共鳴画像のCINEシーケンスを介して前記見
    かけ上の動きをトラッキングする工程を包含する、請求項18に記載の方法。
  20. 【請求項20】 前記動きをトラッキングする際に改良手順を使用する工程
    を包含する、請求項2に記載の方法。
  21. 【請求項21】 前記改良手順において、比較的小さな動きをトラッキング
    するための基準として有するアンカー点を使用する工程を包含する、請求項20
    に記載の方法。
  22. 【請求項22】 前記改良手順において、次の隣接する点から1画素未満分
    離されるそれぞれの点の点シーケンスを使用する工程を包含する、請求項21に
    記載の方法。
  23. 【請求項23】 約動径方向経路において、前記改良手順を実行する工程を
    包含する、請求項22に記載の方法。
  24. 【請求項24】 前記高調波位相画像を計算する際に帯域通過フィルタを使
    用する工程を包含する、請求項1に記載の方法。
  25. 【請求項25】 円形伸張の指標として正のひずみ値を使用し、円形収縮の
    指標として負のひずみ値を使用する工程を包含する、請求項3に記載の方法。
  26. 【請求項26】 動径方向の厚みの増加の指標として正のひずみ値を使用し
    、動径方向の厚みの減少の指標として負のひずみ値を使用する工程を包含する、
    請求項4に記載の方法。
  27. 【請求項27】 前記ひずみの判定を実行する際に、実質的には同時に起こ
    る前記第1の高調波位相および第2の高調波位相を使用する工程を包含する、請
    求項1に記載の方法。
JP2000613339A 1999-04-22 2000-04-14 Cine高調波位相(harp)磁気共鳴画像法を用いた心臓動きのトラッキング Expired - Lifetime JP4679729B2 (ja)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1099297A (ja) * 1996-09-18 1998-04-21 Siemens Corporate Res Inc タグ付けされたmrイメージにおける心臓壁の動きを再生するためのハイブリッドモデルおよび方法

Cited By (9)

* Cited by examiner, † Cited by third party
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