JPH06142079A - 心筋速度の位相差磁気共鳴画像診断法による非侵襲的心筋運動分析法 - Google Patents
心筋速度の位相差磁気共鳴画像診断法による非侵襲的心筋運動分析法Info
- Publication number
- JPH06142079A JPH06142079A JP3236526A JP23652691A JPH06142079A JP H06142079 A JPH06142079 A JP H06142079A JP 3236526 A JP3236526 A JP 3236526A JP 23652691 A JP23652691 A JP 23652691A JP H06142079 A JPH06142079 A JP H06142079A
- Authority
- JP
- Japan
- Prior art keywords
- region
- interest
- magnetic resonance
- velocity
- frame
- 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
Classifications
-
- 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
-
- 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/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
- G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
- G01R33/5617—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using RF refocusing, e.g. RARE
-
- 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/56316—Characterization of motion or flow; Dynamic imaging involving phase contrast techniques
-
- 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/567—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution gated by physiological signals, i.e. synchronization of acquired MR data with periodical motion of an object of interest, e.g. monitoring or triggering system for cardiac or respiratory gating
- G01R33/5676—Gating or triggering based on an MR signal, e.g. involving one or more navigator echoes for motion monitoring and correction
Landscapes
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Health & Medical Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- General Health & Medical Sciences (AREA)
- Vascular Medicine (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Image Processing (AREA)
- Closed-Circuit Television Systems (AREA)
- Image Analysis (AREA)
Abstract
(57)【要約】 (修正有)
【目的】器官の運動を非侵襲的に評価し映像化する、改
良された装置と方法を提供する。位相差磁気共鳴法によ
り心筋運動を評価する、改良された装置と方法を提供す
る。 【構成】運動と変形が、心臓周期などの関心期間におい
て、複数のフレームにおける複数の領域の速度を測定す
る位相差磁気共鳴画像法により、評価され、映像化され
る。関心領域は、フレームからフレームへ追跡され、既
知の最後の位置と既知の位置における速度とにもとずい
ている。張力と回転は、選択された軸に沿った速度変化
から決定される。
良された装置と方法を提供する。位相差磁気共鳴法によ
り心筋運動を評価する、改良された装置と方法を提供す
る。 【構成】運動と変形が、心臓周期などの関心期間におい
て、複数のフレームにおける複数の領域の速度を測定す
る位相差磁気共鳴画像法により、評価され、映像化され
る。関心領域は、フレームからフレームへ追跡され、既
知の最後の位置と既知の位置における速度とにもとずい
ている。張力と回転は、選択された軸に沿った速度変化
から決定される。
Description
【0001】
【産業上の利用分野】本発明は、一般的には核磁気共鳴
画像診断法(MRI)に関し、具体的には、本発明は、
器官内の組織の速度の位相差核磁気共鳴画像診断マップ
を使用して、心臓など器官の運動分析と画像診断法に関
する。
画像診断法(MRI)に関し、具体的には、本発明は、
器官内の組織の速度の位相差核磁気共鳴画像診断マップ
を使用して、心臓など器官の運動分析と画像診断法に関
する。
【0002】
【従来の技術】心筋運動の評価は、いくつかの心臓病理
学の特徴づけと、好結果の干渉の発展と評価とにとって
絶対に必要である。エコー心電図検査は、心筋力学の臨
床研究には好適な方法になっており、内壁の運動と肥大
化の局部測定を行う。エコー心電図検査法と、CT及び
核磁気共鳴画像診断法などの断層撮影法の基本的欠点
は、これらの方法が、一定の外部基準方式を採用してお
り、心筋の個々の部分を追跡出来ないことである。これ
らの方法により、一定の画像面に現われる心筋が映像化
される。従って、映像面内に見られる心筋組織が、心臓
周期を通じて、組織の異なるサンプルを実際に示してい
る。その結果、内壁の運動と肥大化の量が、間違って受
けとられる。心筋内に定められた特定の点の正確な識別
によって、この制約は適切に解決される。この情報を与
えることが出来る一つの方法は、注入された放射線不透
過標識を使用しているが、非常に侵襲的で健康的組織を
冒す。
学の特徴づけと、好結果の干渉の発展と評価とにとって
絶対に必要である。エコー心電図検査は、心筋力学の臨
床研究には好適な方法になっており、内壁の運動と肥大
化の局部測定を行う。エコー心電図検査法と、CT及び
核磁気共鳴画像診断法などの断層撮影法の基本的欠点
は、これらの方法が、一定の外部基準方式を採用してお
り、心筋の個々の部分を追跡出来ないことである。これ
らの方法により、一定の画像面に現われる心筋が映像化
される。従って、映像面内に見られる心筋組織が、心臓
周期を通じて、組織の異なるサンプルを実際に示してい
る。その結果、内壁の運動と肥大化の量が、間違って受
けとられる。心筋内に定められた特定の点の正確な識別
によって、この制約は適切に解決される。この情報を与
えることが出来る一つの方法は、注入された放射線不透
過標識を使用しているが、非常に侵襲的で健康的組織を
冒す。
【0003】磁気共鳴画像診断法は非侵襲的であり、全
体的心筋機能、心室、及び局部内壁肥大化の正確な測定
を行うことが明らかにされた。磁気共鳴画像診断法、聴
覚的窓を必要とする制約を受けないので、すべての実質
的映像面によるエコー心電図検査法と同じ情報を提供す
ることが出来る。その上、磁気共鳴画像診断法は、全領
域の解剖学的画像を提供し、操作者の技量に余り左右さ
れない。心臓の研究に非常に有用であることが発見され
た一つの磁気共鳴画像診断法は、映写磁気共鳴画像診断
法(CINE MRI)と呼ばれており、米国特許No.
4,710,717 に記載されている。簡単に言えば、この方法
により、データは高速で得られ、位相符号化勾配の大き
さの増分は、生理学的制御(例えば、心電図検査)によ
り制御される。これが行われている間、各エコーが得ら
れる心臓周期内の時間的位置も測定される。この時間情
報と補間法とにより、心臓周期中の対象の状態を描写す
る画像が形成される。心臓の研究には有用ではあるが、
これらの画像には、上述の欠点がある。
体的心筋機能、心室、及び局部内壁肥大化の正確な測定
を行うことが明らかにされた。磁気共鳴画像診断法、聴
覚的窓を必要とする制約を受けないので、すべての実質
的映像面によるエコー心電図検査法と同じ情報を提供す
ることが出来る。その上、磁気共鳴画像診断法は、全領
域の解剖学的画像を提供し、操作者の技量に余り左右さ
れない。心臓の研究に非常に有用であることが発見され
た一つの磁気共鳴画像診断法は、映写磁気共鳴画像診断
法(CINE MRI)と呼ばれており、米国特許No.
4,710,717 に記載されている。簡単に言えば、この方法
により、データは高速で得られ、位相符号化勾配の大き
さの増分は、生理学的制御(例えば、心電図検査)によ
り制御される。これが行われている間、各エコーが得ら
れる心臓周期内の時間的位置も測定される。この時間情
報と補間法とにより、心臓周期中の対象の状態を描写す
る画像が形成される。心臓の研究には有用ではあるが、
これらの画像には、上述の欠点がある。
【0004】磁気共鳴画像診断法の画像は、一般に、ス
ピン密度と緩和効果とに依存する画像強度を有する。し
かし、この強度が速度に比例する画像を形成する方法
も、また、例証されている。例えば、オー.ドンネル
著、医療物理学、12、59〜64、1985年;スピ
ーツァー著、放射線学、176、255〜262、19
90年、及びネーラー他著、コンピュータ断層写真法、
10、715〜722、1986年を参照。一般に、こ
れらの方法は、位相差磁気共鳴画像診法と呼ばれる方法
のクラスに属する。ノールバート・ジェイ・ペール他の
“3方向位相差磁気共鳴画像診断の符号化戦略”MRI
出版誌、Vol.4 405−413、7月−8月、199
1年は、速度の3成分を同時に測定する、特に有用にし
て効率的な方法と、その方法を行う装置を開示してい
る。位相差の原理は、上述の映写画像診断法と組合せら
れており、心臓周期中の多数点における速度分布を描写
する画像の形成を可能にしている。
ピン密度と緩和効果とに依存する画像強度を有する。し
かし、この強度が速度に比例する画像を形成する方法
も、また、例証されている。例えば、オー.ドンネル
著、医療物理学、12、59〜64、1985年;スピ
ーツァー著、放射線学、176、255〜262、19
90年、及びネーラー他著、コンピュータ断層写真法、
10、715〜722、1986年を参照。一般に、こ
れらの方法は、位相差磁気共鳴画像診法と呼ばれる方法
のクラスに属する。ノールバート・ジェイ・ペール他の
“3方向位相差磁気共鳴画像診断の符号化戦略”MRI
出版誌、Vol.4 405−413、7月−8月、199
1年は、速度の3成分を同時に測定する、特に有用にし
て効率的な方法と、その方法を行う装置を開示してい
る。位相差の原理は、上述の映写画像診断法と組合せら
れており、心臓周期中の多数点における速度分布を描写
する画像の形成を可能にしている。
【0005】しかし、最近導入された心筋標識法は、注
入された標識から引き出されたと同じ様な個々の心筋部
位の運動に関する情報を得る非侵襲的方法を可能にして
いる。ツエールホーニィ他著、“人体の心臓:磁気共鳴
画像による標識づけ−心筋運動の非侵襲的評価のための
方法”、放射線学、1988年、169、59〜63、
及びエクセル著“磁気化の空間変調による運動の磁気共
鳴画像”放射線学、1989年、171、841−84
5を参照。
入された標識から引き出されたと同じ様な個々の心筋部
位の運動に関する情報を得る非侵襲的方法を可能にして
いる。ツエールホーニィ他著、“人体の心臓:磁気共鳴
画像による標識づけ−心筋運動の非侵襲的評価のための
方法”、放射線学、1988年、169、59〜63、
及びエクセル著“磁気化の空間変調による運動の磁気共
鳴画像”放射線学、1989年、171、841−84
5を参照。
【0006】しかし、心筋標識法にもとずく磁気共鳴画
像診断が現在行われているが、それにはいくつかの欠点
がある。スピン・エコー画像を使用すると、3次元を求
めることが困難であるように、時間的分解能が制約さ
れ、試験時間が延長する。“位置した帯域(sat band
s)”の厚さ、数、及び間隔は、心筋運動分析の空間分解
能を限定する。運動分析の分解能は、さらに、画像の空
間分解能を制約する。運動をすべての3次元で可視化す
る方法は、今なお、例証されなければならない。最終的
に、標識の対比は、緩和時間T′により衰退し、これに
より、心臓周期全体の試験は困難になる。
像診断が現在行われているが、それにはいくつかの欠点
がある。スピン・エコー画像を使用すると、3次元を求
めることが困難であるように、時間的分解能が制約さ
れ、試験時間が延長する。“位置した帯域(sat band
s)”の厚さ、数、及び間隔は、心筋運動分析の空間分解
能を限定する。運動分析の分解能は、さらに、画像の空
間分解能を制約する。運動をすべての3次元で可視化す
る方法は、今なお、例証されなければならない。最終的
に、標識の対比は、緩和時間T′により衰退し、これに
より、心臓周期全体の試験は困難になる。
【0007】ウエディーン(医療の磁気共鳴協会の会誌
1990年8月、462頁)は、歪に関する変数が、速
度の磁気共鳴画像のマップが得られることを明らかにし
た。しかし、単一の時間フレームだけが扱われており、
多くの時間フレームがどのように組合せられるかについ
ては、教示されていない。さらに、報告書は、一定の回
転と変換に関する開発情報の存在を認識していない。
1990年8月、462頁)は、歪に関する変数が、速
度の磁気共鳴画像のマップが得られることを明らかにし
た。しかし、単一の時間フレームだけが扱われており、
多くの時間フレームがどのように組合せられるかについ
ては、教示されていない。さらに、報告書は、一定の回
転と変換に関する開発情報の存在を認識していない。
【0008】
【課題を解決するための手段】本発明は、心筋運動を評
価する改良された方法を意図するもので、この方法は、
音によってのみ制約される運動の分解能により、磁気共
鳴標識化の減衰する問題もなく、また、精密なパターン
認識のアルゴリズムの必要もなく、選択された領域の運
動、速度、回転、及び剪断を決定する反復的方法を採用
している。従来技術の標識化磁気共鳴画像法と標識のい
ずれにも比較して、かなりの大きい心筋領域が追跡さ
れ、反復して選択することが出来る。人体の心筋の比較
的高い3次元運動分析を使用出来ることにより、トル
ク、捩れ、剪断及び運動に関するデータが、正常な病理
学的状態の下で非侵襲的に得られる。初めて、右心室の
非侵襲的研究が行うことが出来る。患者の前後の介入を
非侵襲的に研究することが出来ることにより、各種の治
療と、個々の治療から得られるこれら患者の母集団を最
もよく予測する重要な機会が得られる。非侵襲的に入手
出来る、複雑な生来の心臓病の正常な履歴に関する基本
的生理の情報によって、治療の選択と時機が改善され
る。
価する改良された方法を意図するもので、この方法は、
音によってのみ制約される運動の分解能により、磁気共
鳴標識化の減衰する問題もなく、また、精密なパターン
認識のアルゴリズムの必要もなく、選択された領域の運
動、速度、回転、及び剪断を決定する反復的方法を採用
している。従来技術の標識化磁気共鳴画像法と標識のい
ずれにも比較して、かなりの大きい心筋領域が追跡さ
れ、反復して選択することが出来る。人体の心筋の比較
的高い3次元運動分析を使用出来ることにより、トル
ク、捩れ、剪断及び運動に関するデータが、正常な病理
学的状態の下で非侵襲的に得られる。初めて、右心室の
非侵襲的研究が行うことが出来る。患者の前後の介入を
非侵襲的に研究することが出来ることにより、各種の治
療と、個々の治療から得られるこれら患者の母集団を最
もよく予測する重要な機会が得られる。非侵襲的に入手
出来る、複雑な生来の心臓病の正常な履歴に関する基本
的生理の情報によって、治療の選択と時機が改善され
る。
【0009】従って、本発明の目的は、器官の運動を非
侵襲的に評価し映像化する、改良された装置と方法を提
供することである。本発明のほかの目的は、位相差磁気
共鳴法により心筋運動を評価する、改良された装置と方
法を提供することである。本発明の特徴は、器官の小部
分の位置、速度及び変形を時間の関数として測定するこ
とにより、器官の運動を評価することである。
侵襲的に評価し映像化する、改良された装置と方法を提
供することである。本発明のほかの目的は、位相差磁気
共鳴法により心筋運動を評価する、改良された装置と方
法を提供することである。本発明の特徴は、器官の小部
分の位置、速度及び変形を時間の関数として測定するこ
とにより、器官の運動を評価することである。
【0010】本発明と目的、及びその特徴は、図面を参
照して、次の詳細な説明と添付請求の範囲から、さらに
容易に明らかになるであろう。
照して、次の詳細な説明と添付請求の範囲から、さらに
容易に明らかになるであろう。
【0011】
【実施例】磁気共鳴データが取得されると、各画素が大
きさと位相とを有する複合された量である画像が形成す
る。大きは、多くの要素内の正味の横断磁化の大きさに
相当する。普通の磁気共鳴画像法の場合、位相は重要で
なく、放棄される。核磁気共鳴による映像化の開発のか
なり前に、磁界の勾配は、スピンが動くとき、スピンは
変化している磁界を受けるので、運動情報を核磁気共鳴
信号の位相へ符号化するために、使用出来ることが知ら
れていた。対象がx方向に動いているとすると、その位
置は、時間の関数として、テイラー級数に展開される。
すなわち、 x(t)=xo +Vt+at2 /2+…… (1) ここで、xo 、V、及びaは、それぞれ、時間t=0に
おける位置、速度、加速度、である。時間t′における
スピンの位相は、x磁界の勾配Gx の存在により、次式
で表される。
きさと位相とを有する複合された量である画像が形成す
る。大きは、多くの要素内の正味の横断磁化の大きさに
相当する。普通の磁気共鳴画像法の場合、位相は重要で
なく、放棄される。核磁気共鳴による映像化の開発のか
なり前に、磁界の勾配は、スピンが動くとき、スピンは
変化している磁界を受けるので、運動情報を核磁気共鳴
信号の位相へ符号化するために、使用出来ることが知ら
れていた。対象がx方向に動いているとすると、その位
置は、時間の関数として、テイラー級数に展開される。
すなわち、 x(t)=xo +Vt+at2 /2+…… (1) ここで、xo 、V、及びaは、それぞれ、時間t=0に
おける位置、速度、加速度、である。時間t′における
スピンの位相は、x磁界の勾配Gx の存在により、次式
で表される。
【0012】Δφ(t′)=γMo xo +γM1 V +γM2 a/2+…… (2) ここで、γは磁気回転比であり、Mo、M1 及びM
2 は、それぞれ、勾配波形の0次、1次、2次のモーメ
ントである。すなわち、 M0= ∫0 t' Gx(t)dt, M1=∫0 t' Gx(t)tdt, M2= ∫0 t' Gx(t)t2dt,(3) これは、例えば、磁界の非同次性による静止位相移動の
潜在性を無視している。式(1)のMo 項は、空間的局
地化に関する項である。Mo の各種レベルにおける測定
は、フーリェ変換により、核磁気共鳴信号を各種の値の
x(すなわち、各種位置)における成分に分離する。従
って、Mo 項は画像の構造に取り入れられ、画像におい
て観測された位相のずれは、残りの項による。すなわ
ち、 Δφ(t′)=γM0 V +γM2 a/2+…… (4) 静止非同次性の影響を消去するため、各種の1次モーメ
ントにより行われた測定値の引算が使用される。加速度
とより高次の運動、あるいは、より高次のモーメントが
小さいとするならば、観測された位相のずれは、速度に
比例する。比定定数はγΔM1 であり、ここで、ΔM1
は1次モーメントの変化量である。γΔM1 は、単位速
度当り位相の大きさを有しており、流れあるいは運動の
符号化と呼ばれている。比例定数γΔM1 は、核磁気共
鳴実験の変数により完全に決定され、高精度で知られて
いることについて、留意する必要がある。従って、この
現象が、速度の測定に利用されるならば、現象は、人為
的較正、あるいは、外部の基準を必要としない速度の絶
対測定値生成する。心筋追跡 単一方向(例えば、X)運動の位相差映写画像は、その
方向で心臓周期において一様な間隔で、速度の画素マッ
プにより画素を生成する。これらのデータにより、一つ
の画素程度に小さい心筋サンプルの空間軌道の反復的誘
導が可能になる。フレーム1に選択された小さい関心領
域(ROI)は、そのフレーム内にサンプルの位置を設
定する。そのフレーム内の画像の画素値から誘導された
サンプルの速度とフレーム間の時間Δtとの積に加え
て、この関心領域は、次のフレームのサンプルの期待位
置の予測値を設定する。この過程は、繰り返されて、心
臓周期全体を通してサンプルの経路を反復して形成す
る。
2 は、それぞれ、勾配波形の0次、1次、2次のモーメ
ントである。すなわち、 M0= ∫0 t' Gx(t)dt, M1=∫0 t' Gx(t)tdt, M2= ∫0 t' Gx(t)t2dt,(3) これは、例えば、磁界の非同次性による静止位相移動の
潜在性を無視している。式(1)のMo 項は、空間的局
地化に関する項である。Mo の各種レベルにおける測定
は、フーリェ変換により、核磁気共鳴信号を各種の値の
x(すなわち、各種位置)における成分に分離する。従
って、Mo 項は画像の構造に取り入れられ、画像におい
て観測された位相のずれは、残りの項による。すなわ
ち、 Δφ(t′)=γM0 V +γM2 a/2+…… (4) 静止非同次性の影響を消去するため、各種の1次モーメ
ントにより行われた測定値の引算が使用される。加速度
とより高次の運動、あるいは、より高次のモーメントが
小さいとするならば、観測された位相のずれは、速度に
比例する。比定定数はγΔM1 であり、ここで、ΔM1
は1次モーメントの変化量である。γΔM1 は、単位速
度当り位相の大きさを有しており、流れあるいは運動の
符号化と呼ばれている。比例定数γΔM1 は、核磁気共
鳴実験の変数により完全に決定され、高精度で知られて
いることについて、留意する必要がある。従って、この
現象が、速度の測定に利用されるならば、現象は、人為
的較正、あるいは、外部の基準を必要としない速度の絶
対測定値生成する。心筋追跡 単一方向(例えば、X)運動の位相差映写画像は、その
方向で心臓周期において一様な間隔で、速度の画素マッ
プにより画素を生成する。これらのデータにより、一つ
の画素程度に小さい心筋サンプルの空間軌道の反復的誘
導が可能になる。フレーム1に選択された小さい関心領
域(ROI)は、そのフレーム内にサンプルの位置を設
定する。そのフレーム内の画像の画素値から誘導された
サンプルの速度とフレーム間の時間Δtとの積に加え
て、この関心領域は、次のフレームのサンプルの期待位
置の予測値を設定する。この過程は、繰り返されて、心
臓周期全体を通してサンプルの経路を反復して形成す
る。
【0013】時間を含む点Vvc(x、y、z、t)にお
ける各画素の速度ベクトルVvc=(Vx、Vy、Vz)
のマップを形成するように、心筋運動の反復追跡は、心
臓周期において単一時間の三つの速度成分の画像を処理
することにより、3次元に普遍化される。γVC=(x、
y、z)とすると、これはVvc(γvc、t)と書くこと
も出来る。なお、ここで添え字VCはベクトルを意味す
る。i番目のフレーム内のサンプルの時間ti における
位置γi vc=(xi 、yi 、zi )から、i+1番目の
フレームの位置γi vc+1は次のように計算される。
ける各画素の速度ベクトルVvc=(Vx、Vy、Vz)
のマップを形成するように、心筋運動の反復追跡は、心
臓周期において単一時間の三つの速度成分の画像を処理
することにより、3次元に普遍化される。γVC=(x、
y、z)とすると、これはVvc(γvc、t)と書くこと
も出来る。なお、ここで添え字VCはベクトルを意味す
る。i番目のフレーム内のサンプルの時間ti における
位置γi vc=(xi 、yi 、zi )から、i+1番目の
フレームの位置γi vc+1は次のように計算される。
【0014】 γi vc+1=γi vc+Vvc(γi vc,ti )Δt 上式の反復積分は、反対方向に行うことが出来ること
は、本技術の専門家には明らかであろう。 γi −1vc=γi vc+Vvc(γi vc,ti )Δt 形成される軌道の精度は、速度の推定値にの精度に依存
する。従って、単一画素の軌道は使用することが出来る
が、平均値が使用されるならば、領域の大きさは増大す
るので、向上した精度が得られる。この平均値は、領域
の平均的運動であり、従って、この関心領域は、全体的
に異なる挙動の構造体を含むほどに大きくてはならな
い。容積は選択可能であり、下限の値で0.1ccよりはる
かに小さく制限されるが、0.25ccより小さい容積の領
域が、心筋運動の研究には好適である。空間分解能の考察 本発明の位相差法により現在使用可能な速度マップの空
間分解能は、切断厚さ3から10mmで、平均で1から2
mmである。この寸法は、所望の情報と比較して粗いが、
本位相差法は、準画素の変位とそれに対応する歪みを、
空間分解能により制約されない程度まで測定することが
出来る。各フレームの各画素は、画素の中心点の速度に
関する情報を有する。この速度とフレームの持続時間と
の積は、一つのフレーム中に起った移動量である。予測
された移動量の精度と正確さは、速度推定値の精度と正
確さとによってのみ、限定される。速度の推定値が、フ
レーム時間における点がうける速度を表すならば、速度
マップの空間分解能は重要ではない。従って、次のフレ
ームのそのときのサンプルの位置は、微小な画素として
推定される。次のフレーム上の移動を計算するために、
次のフレームのサンプルの速度の推定が必要である。こ
の画素データから、これは、補間法(例えば、3次元、
2次元、あるいは、最も近い隣接点)により得られ、こ
の推定された速度は、増加移動量の計算に使用される。
この推定された速度の妥当性は、この推定速度がサンプ
ルに対し適切である程度に依存する。画像の分解能が、
心筋の速度の空間変動を特性づけるに十分である限り、
これは許容される。変形 個々の心筋サンプルを追跡する能力は、非常に有用であ
るが、しばしば、病理学的情報が、局部的変形に含まれ
ている。例えば、心臓が収縮すると、筋肉は、いくつか
の方向に縮小し、ほかの方向に肥大化し、また、剪断変
形を受けると考えられている。虚血性筋肉は期待されて
いるように挙動せず、また、変形しないか、あるいは、
健康な筋肉の場合期待されものと反対の変形を呈するこ
とが、侵襲的研究から知られている。
は、本技術の専門家には明らかであろう。 γi −1vc=γi vc+Vvc(γi vc,ti )Δt 形成される軌道の精度は、速度の推定値にの精度に依存
する。従って、単一画素の軌道は使用することが出来る
が、平均値が使用されるならば、領域の大きさは増大す
るので、向上した精度が得られる。この平均値は、領域
の平均的運動であり、従って、この関心領域は、全体的
に異なる挙動の構造体を含むほどに大きくてはならな
い。容積は選択可能であり、下限の値で0.1ccよりはる
かに小さく制限されるが、0.25ccより小さい容積の領
域が、心筋運動の研究には好適である。空間分解能の考察 本発明の位相差法により現在使用可能な速度マップの空
間分解能は、切断厚さ3から10mmで、平均で1から2
mmである。この寸法は、所望の情報と比較して粗いが、
本位相差法は、準画素の変位とそれに対応する歪みを、
空間分解能により制約されない程度まで測定することが
出来る。各フレームの各画素は、画素の中心点の速度に
関する情報を有する。この速度とフレームの持続時間と
の積は、一つのフレーム中に起った移動量である。予測
された移動量の精度と正確さは、速度推定値の精度と正
確さとによってのみ、限定される。速度の推定値が、フ
レーム時間における点がうける速度を表すならば、速度
マップの空間分解能は重要ではない。従って、次のフレ
ームのそのときのサンプルの位置は、微小な画素として
推定される。次のフレーム上の移動を計算するために、
次のフレームのサンプルの速度の推定が必要である。こ
の画素データから、これは、補間法(例えば、3次元、
2次元、あるいは、最も近い隣接点)により得られ、こ
の推定された速度は、増加移動量の計算に使用される。
この推定された速度の妥当性は、この推定速度がサンプ
ルに対し適切である程度に依存する。画像の分解能が、
心筋の速度の空間変動を特性づけるに十分である限り、
これは許容される。変形 個々の心筋サンプルを追跡する能力は、非常に有用であ
るが、しばしば、病理学的情報が、局部的変形に含まれ
ている。例えば、心臓が収縮すると、筋肉は、いくつか
の方向に縮小し、ほかの方向に肥大化し、また、剪断変
形を受けると考えられている。虚血性筋肉は期待されて
いるように挙動せず、また、変形しないか、あるいは、
健康な筋肉の場合期待されものと反対の変形を呈するこ
とが、侵襲的研究から知られている。
【0015】非常に多くの点に関する位置のデータは、
生成され、各データは放射線不透過マーカーと等価であ
るので、部分的な収縮、内壁の肥大化、剪断張力、ねじ
れの測定値を生成することは、可能である。これは、心
臓拡張期フレームにおける多数の関心領域を選択し、心
臓周期を通してその位置を追跡し、マーカーにより処理
されるように、これらのデータを分析することにより達
成される。しかし、位相差磁気共鳴画像法がデータによ
るいくつかの変数(例えば、ひずみ)、計算のより直接
的方法が発見された。
生成され、各データは放射線不透過マーカーと等価であ
るので、部分的な収縮、内壁の肥大化、剪断張力、ねじ
れの測定値を生成することは、可能である。これは、心
臓拡張期フレームにおける多数の関心領域を選択し、心
臓周期を通してその位置を追跡し、マーカーにより処理
されるように、これらのデータを分析することにより達
成される。しかし、位相差磁気共鳴画像法がデータによ
るいくつかの変数(例えば、ひずみ)、計算のより直接
的方法が発見された。
【0016】一方向(x)の運動のような単純な場合、
x方向に指向し、心臓周期において厳密に動く長さLの
小要素は、小要素に沿って一様(しかし、ゼロでなく、
時間が変化する)である、測定された速度分布(Vx )
を有することが、位相差磁気共鳴画像法により立証され
た。変形を受けている小要素の場合、小要素に沿った速
度は、空間的には一様でない。変形(張力)に関する情
報は、速度の空間勾配あるいは傾斜dVx /dxに含ま
れている。i番目時間フレームを考察することして、小
要素をxo から(xo +L)へ延長する。xo における
速度がVo ならば、(xo +L)における速度は、{V
o +L(dVx /dx)}である。もし、これらの点
が、これらの速度でフレーム時Δtの間続くならば、そ
の距離は、ΔL=L(dVx /dx)Δtだけ増大す
る。このフレームの間に発生した張力εxxは、従って、
次式で表される。
x方向に指向し、心臓周期において厳密に動く長さLの
小要素は、小要素に沿って一様(しかし、ゼロでなく、
時間が変化する)である、測定された速度分布(Vx )
を有することが、位相差磁気共鳴画像法により立証され
た。変形を受けている小要素の場合、小要素に沿った速
度は、空間的には一様でない。変形(張力)に関する情
報は、速度の空間勾配あるいは傾斜dVx /dxに含ま
れている。i番目時間フレームを考察することして、小
要素をxo から(xo +L)へ延長する。xo における
速度がVo ならば、(xo +L)における速度は、{V
o +L(dVx /dx)}である。もし、これらの点
が、これらの速度でフレーム時Δtの間続くならば、そ
の距離は、ΔL=L(dVx /dx)Δtだけ増大す
る。このフレームの間に発生した張力εxxは、従って、
次式で表される。
【0017】εxx=ΔL/L=(dVx /dx)Δt 正のεxxは、小要素がこのフレームの間、伸長されたこ
とを示し、負のεxxは収縮されたことを表す。連続的変
数の場合、dVx /dx=dεxx/dtである。従っ
て、各時間フレームの、x方向のVx の導関数は、時間
の導関数、あるいは、x方向縦の張力の成長速度であ
る。
とを示し、負のεxxは収縮されたことを表す。連続的変
数の場合、dVx /dx=dεxx/dtである。従っ
て、各時間フレームの、x方向のVx の導関数は、時間
の導関数、あるいは、x方向縦の張力の成長速度であ
る。
【0018】ほかの方向の成長速度も、次式により計算
することが出来る。 dVy /dy=dεyy/dt 及び dVz /dz=dεzz/dt 剪断張力の成長速度も、推定することが出来る。例え
ば、2次元速度の2次元マップ、すわなち、Vx (x,
y)とVy (x,y)に関して、二つのタイプの剪断変
形が可能であり、θxyとθyxとで示されて、図1Aの正
の方向で変形角として定義される。
することが出来る。 dVy /dy=dεyy/dt 及び dVz /dz=dεzz/dt 剪断張力の成長速度も、推定することが出来る。例え
ば、2次元速度の2次元マップ、すわなち、Vx (x,
y)とVy (x,y)に関して、二つのタイプの剪断変
形が可能であり、θxyとθyxとで示されて、図1Aの正
の方向で変形角として定義される。
【0019】上記と同じ誘導段階により、これらの剪断
張力の成長速度は次式で表されることが明らかになる。 dθyx/dt=dVx /dz及び dθxy/dt=dVx /dy 一般に、どちらのタイプの剪断運動は、同時に存在す
る。すなわちdVx /dyとdVy /dxは、同時には
ゼロでない。
張力の成長速度は次式で表されることが明らかになる。 dθyx/dt=dVx /dz及び dθxy/dt=dVx /dy 一般に、どちらのタイプの剪断運動は、同時に存在す
る。すなわちdVx /dyとdVy /dxは、同時には
ゼロでない。
【0020】このタイプの運動は、実際には、二つの、
一層容易に理解される運動より成っている。純粋な剪断
の場合、θxyはθyxに等しくなければならず、θxy=−
θyxならば、対象物は厳密な回転を呈しており、少しも
変形しない。この運動は、図1Bに示されている。従っ
て、二つの導関数dVx /dyとdVy /dxから、回
転速度dα/dtと純粋な剪断の成長速度dθ/dtと
を計算することが出来る。すなわち、 dα/dt=1/2(dVy /dx−dVx /dy) dθ/dt=(dVy /dx−dVx /dy) 軸からの角度はθ/2であるので、図1Bに定義された
剪断角は、異常であるように見える。しかし、この定義
は、例えば、θyx=0ならばθ=θxyであることを保証
している。
一層容易に理解される運動より成っている。純粋な剪断
の場合、θxyはθyxに等しくなければならず、θxy=−
θyxならば、対象物は厳密な回転を呈しており、少しも
変形しない。この運動は、図1Bに示されている。従っ
て、二つの導関数dVx /dyとdVy /dxから、回
転速度dα/dtと純粋な剪断の成長速度dθ/dtと
を計算することが出来る。すなわち、 dα/dt=1/2(dVy /dx−dVx /dy) dθ/dt=(dVy /dx−dVx /dy) 軸からの角度はθ/2であるので、図1Bに定義された
剪断角は、異常であるように見える。しかし、この定義
は、例えば、θyx=0ならばθ=θxyであることを保証
している。
【0021】2次元データのこれらの分析結果を要約す
ると、領域内のxとyの速度の平均値は、この領域と関
係がある。組合せ(dVy /dx−dVx /dy)は、
この領域の厳密な回転と関係がある。dVx /dxとd
Vy /dyは、縦の張力の成長速度と関係があり、(d
Vy /dx+dVx /dy)は、純粋な剪断張力の成長
速度である。
ると、領域内のxとyの速度の平均値は、この領域と関
係がある。組合せ(dVy /dx−dVx /dy)は、
この領域の厳密な回転と関係がある。dVx /dxとd
Vy /dyは、縦の張力の成長速度と関係があり、(d
Vy /dx+dVx /dy)は、純粋な剪断張力の成長
速度である。
【0022】この式は、3次元の運動、回転、変形へ容
易に一般化される。単一の時間フレームにおける三つの
速度成分の画像から、速度の空間導関数の画像が、作成
され、これらの変形の成長速度の空間分布を示す。例え
ば、dVx/dxの画像は、x方向の縦張力の成長速度
の空間分布を示す。明るい画素値は、この時間フレーム
でx方向に伸長している領域を示し、暗い(負)画素値
は、収縮の速さを示す。同様に、ほかの縦方向の剪断張
力成長速度の画像は、形成される。動的フォーマットで
表示されると、これらの画像は、心筋の収縮が空間と時
間にどのように伝播するかを示す。領域張力 張力成長速度は、心筋追跡法と組合せられて、全心臓周
期あるいはその一部において、成長した全張力を推定す
る。端末心臓拡張期フレームで識別された心筋領域に関
して、各フレームの速度マップは、反復して使用され、
領域をフレームからフレームへ追跡する。各フレームの
(運動中の)領域内の速度勾配は、各フレームの各張力
成分の成長速度を計算するために使用される。このプロ
セスは、心臓周期を通して続行され、小さい容積の心筋
内の縦方向の剪断張力を推定する。これらから、主要な
張力と最も大きい収縮の方向とが、計算される。心筋領
域の形状と位置は、導関数の計算を可能にする重要な各
方向の、最も小さい少なくとも二つの画素により、任意
に選択することが出来る。
易に一般化される。単一の時間フレームにおける三つの
速度成分の画像から、速度の空間導関数の画像が、作成
され、これらの変形の成長速度の空間分布を示す。例え
ば、dVx/dxの画像は、x方向の縦張力の成長速度
の空間分布を示す。明るい画素値は、この時間フレーム
でx方向に伸長している領域を示し、暗い(負)画素値
は、収縮の速さを示す。同様に、ほかの縦方向の剪断張
力成長速度の画像は、形成される。動的フォーマットで
表示されると、これらの画像は、心筋の収縮が空間と時
間にどのように伝播するかを示す。領域張力 張力成長速度は、心筋追跡法と組合せられて、全心臓周
期あるいはその一部において、成長した全張力を推定す
る。端末心臓拡張期フレームで識別された心筋領域に関
して、各フレームの速度マップは、反復して使用され、
領域をフレームからフレームへ追跡する。各フレームの
(運動中の)領域内の速度勾配は、各フレームの各張力
成分の成長速度を計算するために使用される。このプロ
セスは、心臓周期を通して続行され、小さい容積の心筋
内の縦方向の剪断張力を推定する。これらから、主要な
張力と最も大きい収縮の方向とが、計算される。心筋領
域の形状と位置は、導関数の計算を可能にする重要な各
方向の、最も小さい少なくとも二つの画素により、任意
に選択することが出来る。
【0023】多くの時間フレームから変形寄与を組合せ
る場合、恒久的な合理的な座標系を維持することが重要
である。例えば、一つの時間フレームにおいて、筋肉繊
維がx方向に指向し、局部的収縮がx方向の張力に反映
されている仮定する。また、周囲の筋肉の働きにより、
検査中の領域が回転していると仮定する。次の時間フレ
ームにおいて、張力の成長速度は、再度計算される。し
かし、最初のフレームの結果と組合せる前に、この第2
の時間フレームの結果は、第1の結果との整合へ回帰さ
れなければならない。
る場合、恒久的な合理的な座標系を維持することが重要
である。例えば、一つの時間フレームにおいて、筋肉繊
維がx方向に指向し、局部的収縮がx方向の張力に反映
されている仮定する。また、周囲の筋肉の働きにより、
検査中の領域が回転していると仮定する。次の時間フレ
ームにおいて、張力の成長速度は、再度計算される。し
かし、最初のフレームの結果と組合せる前に、この第2
の時間フレームの結果は、第1の結果との整合へ回帰さ
れなければならない。
【0024】図2は、時間間隔における各種フレームの
領域の速度測定値にもとずいた時間間隔の間、関心領域
の動きと変形の全体を計算する段階を示す流れ図であ
る。説明を明確にするために、2次元の運動に限定され
ているが、3次元の運動への拡張は、本技術の専門家に
は明らかである。最初に、フレームカウンタiとして、
回転角α、xとyの変位、ΔxとΔy、ずれ(θ)、及
び縦方向の張力(εxx、εyy)は、ゼロに設定される。
領域の速度測定値にもとずいた時間間隔の間、関心領域
の動きと変形の全体を計算する段階を示す流れ図であ
る。説明を明確にするために、2次元の運動に限定され
ているが、3次元の運動への拡張は、本技術の専門家に
は明らかである。最初に、フレームカウンタiとして、
回転角α、xとyの変位、ΔxとΔy、ずれ(θ)、及
び縦方向の張力(εxx、εyy)は、ゼロに設定される。
【0025】第1フレーム(i=1)の場合、関心領域
(ROI)は選択されなければならない。各時間フレー
ムにおいて、関心領域内の速度は次のようにモデル化さ
れる。 Vx =a+b(x−xav)+c(y−yav) Vy =d+e(x−xav)+f(y−yav) xavとyavは、関心領域内のxとyの平均座標値であ
る。この時間フレームの関心領域内の画素値を使用し
て、変数a、b、c、d、e、fは、例えば最小二乗法
を使用して推定される。従って、このモデルから、xの
平均速度、dVx /dx、dVx /dy、yの平均速
度、dVx /dx、dVy /dyの適合値は、それぞ
れ、a、b、c、d、e、fである。
(ROI)は選択されなければならない。各時間フレー
ムにおいて、関心領域内の速度は次のようにモデル化さ
れる。 Vx =a+b(x−xav)+c(y−yav) Vy =d+e(x−xav)+f(y−yav) xavとyavは、関心領域内のxとyの平均座標値であ
る。この時間フレームの関心領域内の画素値を使用し
て、変数a、b、c、d、e、fは、例えば最小二乗法
を使用して推定される。従って、このモデルから、xの
平均速度、dVx /dx、dVx /dy、yの平均速
度、dVx /dx、dVy /dyの適合値は、それぞ
れ、a、b、c、d、e、fである。
【0026】このようにして、各フレームiに関して
は、増分の縦方向張力Δεxx,iとΔε yy,i、変位Δxi
とΔyi、回転Δαi 、剪断張力Δθi が計算される。
張力テンソルは計算される。すなわち、εは、
は、増分の縦方向張力Δεxx,iとΔε yy,i、変位Δxi
とΔyi、回転Δαi 、剪断張力Δθi が計算される。
張力テンソルは計算される。すなわち、εは、
【0027】
【化1】
【0028】となり、実験室的座標系に関して、最初の
フレームと一致するように、−αだけ回転しなければな
らない。この回転により、回転張力ε′、このフレーム
の、回転した縦方向の剪断張力Δε′xx,iとΔ
ε′yy,i、及びΔθ′i が形成し、次式が成立する。
フレームと一致するように、−αだけ回転しなければな
らない。この回転により、回転張力ε′、このフレーム
の、回転した縦方向の剪断張力Δε′xx,iとΔ
ε′yy,i、及びΔθ′i が形成し、次式が成立する。
【0029】
【化2】
【0030】次に、連続している、全縦方向張力εxxと
εyy及び全剪張力Δθは、それぞれ、Δε′xx,i、Δ
ε′yy,i及びΔθ′i だけ増分する。第1フレームに関
して、α=0であり、回転が影響しないことに留意する
ことが必要である。しかし、この後で明らかになるが、
ほかのフレームに関しては、αはゼロではない。
εyy及び全剪張力Δθは、それぞれ、Δε′xx,i、Δ
ε′yy,i及びΔθ′i だけ増分する。第1フレームに関
して、α=0であり、回転が影響しないことに留意する
ことが必要である。しかし、この後で明らかになるが、
ほかのフレームに関しては、αはゼロではない。
【0031】これに従って、次のフレームの用意とし
て、正味の変換量ΔxとΔy、及び正味の回転αは、そ
れぞれ、Δxi、Δyi、及びΔαiだけ増分する。フ
レーム・カウンタが、すべてのフレームが処理されたこ
とを示すならば、計算は完了している。さもなければ、
初期の関心領域の形状は、ΔxとΔyだけ変換され、α
だけ回転し、εxx、εyy、及びθだけ変形して、次のフ
レームが処理される。
て、正味の変換量ΔxとΔy、及び正味の回転αは、そ
れぞれ、Δxi、Δyi、及びΔαiだけ増分する。フ
レーム・カウンタが、すべてのフレームが処理されたこ
とを示すならば、計算は完了している。さもなければ、
初期の関心領域の形状は、ΔxとΔyだけ変換され、α
だけ回転し、εxx、εyy、及びθだけ変形して、次のフ
レームが処理される。
【0032】次に、画像は構成されて、計算された変
位、回転、張力にもとづいて表示している時間の間、サ
ンプルの関心領域を示す。図3は、フレームiとフレー
ムi+1の幅Wと高さHの関心領域10の測定・計算さ
れた変数を示す。フレームの間には、関心領域が、Δx
とΔyだけ変換され、αだけ回転し、ΔLだけ伸長し、
−ΔHだけ収縮している。示されていないが、関心領域
は剪断も受けている。図2から計算された多様な変数
は、運動の映写評価に使用される。
位、回転、張力にもとづいて表示している時間の間、サ
ンプルの関心領域を示す。図3は、フレームiとフレー
ムi+1の幅Wと高さHの関心領域10の測定・計算さ
れた変数を示す。フレームの間には、関心領域が、Δx
とΔyだけ変換され、αだけ回転し、ΔLだけ伸長し、
−ΔHだけ収縮している。示されていないが、関心領域
は剪断も受けている。図2から計算された多様な変数
は、運動の映写評価に使用される。
【0033】3軸(x、y、z)方向の測定値と計算値
を得ることにより、心臓周期中の動き、ひずみ、ねじれ
が評価され、可視化される。前述のように、心臓サンプ
ルが本発明により追跡される空間分解能は画像により制
約され、画像の大きさに対しあまり敏感でない。心筋追
跡を十分に行うために画像に必要な空間分解能は、測定
された変位の大きさは依存しないが、運動速度が著しく
変化する大きさに依存する。同様に、張力を適確に測定
するために、空間の変化を特性づけることが出来る画像
分解能が、これらの変数に必要とされる。これは、従来
技術の標識法による放射線不透過マーカー、あるいは、
磁気共鳴画像法による測定法と対照的であり、この場
合、測定された変位の精度は、点を局地化する能力に直
接に依存し、これは、装置の空間分解能とマーカー(放
射線不透過あるいは磁気の)とにより順次制御される。
本発明の方法も、マーカーの位置を自動的に検出する精
密なアルゴリズムを必要としないほかの利点を有してい
る。
を得ることにより、心臓周期中の動き、ひずみ、ねじれ
が評価され、可視化される。前述のように、心臓サンプ
ルが本発明により追跡される空間分解能は画像により制
約され、画像の大きさに対しあまり敏感でない。心筋追
跡を十分に行うために画像に必要な空間分解能は、測定
された変位の大きさは依存しないが、運動速度が著しく
変化する大きさに依存する。同様に、張力を適確に測定
するために、空間の変化を特性づけることが出来る画像
分解能が、これらの変数に必要とされる。これは、従来
技術の標識法による放射線不透過マーカー、あるいは、
磁気共鳴画像法による測定法と対照的であり、この場
合、測定された変位の精度は、点を局地化する能力に直
接に依存し、これは、装置の空間分解能とマーカー(放
射線不透過あるいは磁気の)とにより順次制御される。
本発明の方法も、マーカーの位置を自動的に検出する精
密なアルゴリズムを必要としないほかの利点を有してい
る。
【0034】図4に関し、本発明も取り入れており、ジ
ェネラル・エレクトリック社により商標“シグナ(SI
GNA)”で販売されている、好適な核磁気共鳴装置の
主要な構成要素が構成図の形で示されている。装置の全
動作は、主コンピュータ101(データ・ジェネラル
エム・ヴイ4000)より成る、ホスト・コンピュータ
装置、一般表示100により制御されている。コンピュ
ータ100は、複数のコンピュータ周辺装置と核磁気共
鳴装置とが主コンピュータ101へ接続しているインタ
フェース102より成っている。コンピュータ周辺装置
の間に、主コンピュータ101の命令により、患者デー
タと画像データをテープにとるために、使用される磁気
テープ駆動装置がある。また、処理済の患者データは、
画像ディスク記憶装置110に格納される。アレイ・プ
ロセッサ106は、得られた核磁気共鳴データの事前処
理と画像再構成のために使用される。画像プロセッサ1
08の機能は、相互作用画像表示操作を行うことであ
り、この操作には、拡大、画像比較、グレイ・スケール
調整、及びリアルタイムデータ表示などがある。コンピ
ュータ100も、ディスク・データ記憶装置112を使
用している、核磁気共鳴の生データ(すなわち、画像構
成の前に)を格納する装置より成っている。操作コンソ
ール116も、インタフェース102により主コンピュ
ータ101へ接続しており、患者の研究に適切なデータ
と、較正、スキャンの始動と終了などの、適切な核磁気
共鳴装置に必要な、入力装置によりデータを操作者に提
供する。操作コンソールは、また、ディスクあるいは磁
気テープに記憶された画像を表示するために使用され
る。
ェネラル・エレクトリック社により商標“シグナ(SI
GNA)”で販売されている、好適な核磁気共鳴装置の
主要な構成要素が構成図の形で示されている。装置の全
動作は、主コンピュータ101(データ・ジェネラル
エム・ヴイ4000)より成る、ホスト・コンピュータ
装置、一般表示100により制御されている。コンピュ
ータ100は、複数のコンピュータ周辺装置と核磁気共
鳴装置とが主コンピュータ101へ接続しているインタ
フェース102より成っている。コンピュータ周辺装置
の間に、主コンピュータ101の命令により、患者デー
タと画像データをテープにとるために、使用される磁気
テープ駆動装置がある。また、処理済の患者データは、
画像ディスク記憶装置110に格納される。アレイ・プ
ロセッサ106は、得られた核磁気共鳴データの事前処
理と画像再構成のために使用される。画像プロセッサ1
08の機能は、相互作用画像表示操作を行うことであ
り、この操作には、拡大、画像比較、グレイ・スケール
調整、及びリアルタイムデータ表示などがある。コンピ
ュータ100も、ディスク・データ記憶装置112を使
用している、核磁気共鳴の生データ(すなわち、画像構
成の前に)を格納する装置より成っている。操作コンソ
ール116も、インタフェース102により主コンピュ
ータ101へ接続しており、患者の研究に適切なデータ
と、較正、スキャンの始動と終了などの、適切な核磁気
共鳴装置に必要な、入力装置によりデータを操作者に提
供する。操作コンソールは、また、ディスクあるいは磁
気テープに記憶された画像を表示するために使用され
る。
【0035】コンピュータ装置100は、システム制御
器118と勾配増幅装置128により、核磁気共鳴装置
に対して制御を行う。記憶プログラムの命令により、コ
ンピュータ100は、本技術分野の熟練者には周知の方
法で、順次通信ネットワーク103(例えば、イサーネ
ット(Ethernet) ネットワーク)によりシステム制御器
118と通信する。システム制御器118は、数台の付
属装置より成っており、これらの装置には、パルス制御
モジュール(PCM)120、ラジオ周波数トランシー
バー122、ステータス・コントロール・モジュール
(SCM)124、電源126がある。パルス制御モジ
ュール120は、主コンピュータ装置101によりプロ
グラム制御で発生した制御信号を使用して、勾配コイル
励起を制御するデジタル波形と、無線周波励起パルスを
変調する。トランシーバー122に使用される無線周波
包絡波形とを発生する。勾配波形は、Gx、Gy、Gz
の増幅器130、132、134より構成されている勾
配増幅装置128へ送られる。各増幅器130、13
2、134は、磁石組立体146の一部である組立体1
36内の関連の勾配コイルを励起するために使用され
る。励起されると、勾配コイルは、磁界の勾配Gx、G
y、Gzを形成する。勾配磁界は、研究中の患者の領域
から発している核磁気共鳴信号へ空間情報を符号化する
ために、トランシーバー122、無線周波増幅器12
3、及び無線周波コイル138によって発生した無線周
波パルスと組合せて使用される。パルス制御モジュール
120により形成された波形と制御信号は、無線周波数
の搬送波変調とモード制御のためのに、トランシーバー
122により使用される。送信モードにおいて、トラン
シーバーは、無線周波信号を無線周波電力増幅器123
へ送り、次に、増幅器123は、主磁石組立体146内
にある無線周波コイル138を励起する。患者の励起さ
れたスピンにより発信された核磁気共鳴信号は、送信に
使用される、同一か、あるいは、異なる無線周波コイル
により感知される。信号は、トランシーバー122の受
信部分で、検出、増幅、復調、濾波、及び符号化が行わ
れる。処理された信号は、インターフェース102とト
ランシーバー122を連結している、専用で、単一方向
の高速度デジタル・リンクにより、主コンピュータ10
1へ送られる。
器118と勾配増幅装置128により、核磁気共鳴装置
に対して制御を行う。記憶プログラムの命令により、コ
ンピュータ100は、本技術分野の熟練者には周知の方
法で、順次通信ネットワーク103(例えば、イサーネ
ット(Ethernet) ネットワーク)によりシステム制御器
118と通信する。システム制御器118は、数台の付
属装置より成っており、これらの装置には、パルス制御
モジュール(PCM)120、ラジオ周波数トランシー
バー122、ステータス・コントロール・モジュール
(SCM)124、電源126がある。パルス制御モジ
ュール120は、主コンピュータ装置101によりプロ
グラム制御で発生した制御信号を使用して、勾配コイル
励起を制御するデジタル波形と、無線周波励起パルスを
変調する。トランシーバー122に使用される無線周波
包絡波形とを発生する。勾配波形は、Gx、Gy、Gz
の増幅器130、132、134より構成されている勾
配増幅装置128へ送られる。各増幅器130、13
2、134は、磁石組立体146の一部である組立体1
36内の関連の勾配コイルを励起するために使用され
る。励起されると、勾配コイルは、磁界の勾配Gx、G
y、Gzを形成する。勾配磁界は、研究中の患者の領域
から発している核磁気共鳴信号へ空間情報を符号化する
ために、トランシーバー122、無線周波増幅器12
3、及び無線周波コイル138によって発生した無線周
波パルスと組合せて使用される。パルス制御モジュール
120により形成された波形と制御信号は、無線周波数
の搬送波変調とモード制御のためのに、トランシーバー
122により使用される。送信モードにおいて、トラン
シーバーは、無線周波信号を無線周波電力増幅器123
へ送り、次に、増幅器123は、主磁石組立体146内
にある無線周波コイル138を励起する。患者の励起さ
れたスピンにより発信された核磁気共鳴信号は、送信に
使用される、同一か、あるいは、異なる無線周波コイル
により感知される。信号は、トランシーバー122の受
信部分で、検出、増幅、復調、濾波、及び符号化が行わ
れる。処理された信号は、インターフェース102とト
ランシーバー122を連結している、専用で、単一方向
の高速度デジタル・リンクにより、主コンピュータ10
1へ送られる。
【0036】パルス制御モジュール120とステイタス
制御モジュール124は、独立したサブシステムであ
り、そのどらちも、主コンピュータ101及び患者位置
づけ装置152などの周辺位置と接続し、また、直列通
信リンク103により相互に接続している。パルス制御
モジュール120とステイタス制御モジュール124
は、いずれも、主コンピュータ101からの命令を処理
する16ビットのマイクロプロセッサ(インテル808
6など)より成っている。ステイタス制御モジュール1
24は、患者の初期の位置と移動可能な患者位置合せ光
線ファンビーム(Light fan beam) (図示せず)とに関
する情報を得る手段を有している。この情報は、主コン
ピュータ101に使用されて、画像ディスプレイと復元
変数を修正する。また、ステイタス制御モジュール12
4は、患者の転送と位置合せ装置の作動などの機能を始
動する。勾配コイル組立体136と無線周波送受信コイ
ル138とは、分極化している磁界の形成に使用され磁
石の穴のなかに取り付けられている。磁石は、患者位置
合せ装置148、粗調整コイル電源140、及び主磁石
電源142とより成る主磁石組立体の一部である。主電
源142は、磁石により形成された、分極化している磁
界を1.5テスラの適切な動作強度へもっていくために使
用され、次に、接続を切断される。
制御モジュール124は、独立したサブシステムであ
り、そのどらちも、主コンピュータ101及び患者位置
づけ装置152などの周辺位置と接続し、また、直列通
信リンク103により相互に接続している。パルス制御
モジュール120とステイタス制御モジュール124
は、いずれも、主コンピュータ101からの命令を処理
する16ビットのマイクロプロセッサ(インテル808
6など)より成っている。ステイタス制御モジュール1
24は、患者の初期の位置と移動可能な患者位置合せ光
線ファンビーム(Light fan beam) (図示せず)とに関
する情報を得る手段を有している。この情報は、主コン
ピュータ101に使用されて、画像ディスプレイと復元
変数を修正する。また、ステイタス制御モジュール12
4は、患者の転送と位置合せ装置の作動などの機能を始
動する。勾配コイル組立体136と無線周波送受信コイ
ル138とは、分極化している磁界の形成に使用され磁
石の穴のなかに取り付けられている。磁石は、患者位置
合せ装置148、粗調整コイル電源140、及び主磁石
電源142とより成る主磁石組立体の一部である。主電
源142は、磁石により形成された、分極化している磁
界を1.5テスラの適切な動作強度へもっていくために使
用され、次に、接続を切断される。
【0037】外部からの干渉を最小にするため、磁石、
勾配コイル組立体、及び無線送受信コイル、また、患者
取扱い装置より成る核磁気共鳴装置の構成要素は、無線
周波シールド室一般呼称144内に封じ込められてい
る。この遮蔽は、室全体を取り囲む銅あるいはアルミニ
ウムのスクリーン網により形成されている。スクリーン
網は、装置により発生した無線周波信号を封じ込める働
きをしし、室外で発生した無線周波信号から装置を遮蔽
する。
勾配コイル組立体、及び無線送受信コイル、また、患者
取扱い装置より成る核磁気共鳴装置の構成要素は、無線
周波シールド室一般呼称144内に封じ込められてい
る。この遮蔽は、室全体を取り囲む銅あるいはアルミニ
ウムのスクリーン網により形成されている。スクリーン
網は、装置により発生した無線周波信号を封じ込める働
きをしし、室外で発生した無線周波信号から装置を遮蔽
する。
【0038】特に図4と図5に関して、トランシーバー
122は、無線周波励起磁界B1 を電力増幅器123に
よりコイル138Aに形成する構成要素とコイル138
Bに誘起した結果の核磁気共鳴信号を受信する構成要素
とより成っている。基本、すなわち搬送波である無線周
波励起磁界の周波数は、通信リンク103を経て主コン
ピュータ101から、一組のデジタル信号を受信する周
波数シンセタイザーにより生成する。これらのデジタル
信号は、出力201において1ヘルツの分解能で生成さ
れる周波数を示している。この指令された無線周波搬送
波は、変調器202へ送られる。この場合、送られるの
は、線路203を経て受信した信号に応答した周波数と
振幅であり、また、結果としての無線励起信号は、パル
ス制御モジュール120から線路204を経て受信され
制御信号に応答して、オンとオフに切換えられる。線路
205を経て送られる無線周波励起パルス出力の大きさ
は、主コンピュータ101から通信リンク103を経て
デジタル信号を受信する送信減衰回路206により減衰
する。減衰した無線周波励起パルスは、無線周波送信コ
イル138Aを作動する電力増幅器へ送られる。対象体
内の励起されたスピンにより発生した核磁気共鳴信号
は、受信器コイル138Bにより捕捉され、受信器20
7の入力へ送られる。受信器207は、核磁気共鳴信号
を増幅し、次に、主コンピュータ101からリンク10
3を経て受信されたデジタル減衰信号により決定された
量だけ減衰する。また、核磁気共鳴信号が、行われてい
る特定の取得により必要される時間間隔の間だけ取得さ
れるように、受信器207は、線路208を経てパルス
制御モジュールから送られる信号により、オンとオフに
切換えられる。受信された核磁気共鳴信号は、直角検出
器209により変調され、疑似信号防止(anti-aliasin
g)フィルタ216と217を経て、215で一括して示
されている一対のアナログ−デジタル変換器へ接続して
いる二つの信号IとQを発生する。また、直角検出器2
09は、無線基準信号をもう一つの周波数シンセタイザ
ー210から受信し、この基準信号は、直角検査器20
9により使用され、無線周波基準信号と同じ位相の核磁
気共鳴信号の成分(I信号)の振幅と、基準信号と直角
位相である核磁気共鳴信号の成分(Q信号)の振幅とを
感知する。
122は、無線周波励起磁界B1 を電力増幅器123に
よりコイル138Aに形成する構成要素とコイル138
Bに誘起した結果の核磁気共鳴信号を受信する構成要素
とより成っている。基本、すなわち搬送波である無線周
波励起磁界の周波数は、通信リンク103を経て主コン
ピュータ101から、一組のデジタル信号を受信する周
波数シンセタイザーにより生成する。これらのデジタル
信号は、出力201において1ヘルツの分解能で生成さ
れる周波数を示している。この指令された無線周波搬送
波は、変調器202へ送られる。この場合、送られるの
は、線路203を経て受信した信号に応答した周波数と
振幅であり、また、結果としての無線励起信号は、パル
ス制御モジュール120から線路204を経て受信され
制御信号に応答して、オンとオフに切換えられる。線路
205を経て送られる無線周波励起パルス出力の大きさ
は、主コンピュータ101から通信リンク103を経て
デジタル信号を受信する送信減衰回路206により減衰
する。減衰した無線周波励起パルスは、無線周波送信コ
イル138Aを作動する電力増幅器へ送られる。対象体
内の励起されたスピンにより発生した核磁気共鳴信号
は、受信器コイル138Bにより捕捉され、受信器20
7の入力へ送られる。受信器207は、核磁気共鳴信号
を増幅し、次に、主コンピュータ101からリンク10
3を経て受信されたデジタル減衰信号により決定された
量だけ減衰する。また、核磁気共鳴信号が、行われてい
る特定の取得により必要される時間間隔の間だけ取得さ
れるように、受信器207は、線路208を経てパルス
制御モジュールから送られる信号により、オンとオフに
切換えられる。受信された核磁気共鳴信号は、直角検出
器209により変調され、疑似信号防止(anti-aliasin
g)フィルタ216と217を経て、215で一括して示
されている一対のアナログ−デジタル変換器へ接続して
いる二つの信号IとQを発生する。また、直角検出器2
09は、無線基準信号をもう一つの周波数シンセタイザ
ー210から受信し、この基準信号は、直角検査器20
9により使用され、無線周波基準信号と同じ位相の核磁
気共鳴信号の成分(I信号)の振幅と、基準信号と直角
位相である核磁気共鳴信号の成分(Q信号)の振幅とを
感知する。
【0039】受信された核磁気共鳴信号のIとQの成分
は、取得期間において64KHzの抽出速度で、アナログ
−デジタル変換器218によって、連続的に抽出され
て、デジタル化される。一組の256個のデジタル数字
が、核磁気共鳴信号の各IとQの成分に関して、同時に
得られ、これらのデジタル数字は、直列リンク105を
経て主コンピュータ101へ送られる。各一対の数字
は、I+iQに等しい複素数値と考えられる。図4の核
磁気共鳴装置は、一連のパルス・シーケンスを行い、上
述のように、所望の画像復元に十分な核磁気共鳴データ
を収集する。個々の心筋サンプルと、心臓周期における
複数フレームのこれらサンプル内のひずみとを反復して
追跡することによる位相差磁気共鳴画像法にもとずい
て、心筋運動を評価する、改良された方法と装置が説明
された。また、張力成長の速度に依存する画像形成法も
説明された。本発明は、ほかの器官及び筋肉の運動の研
究に、また、非生物学的サンプルの研究に適用すること
が出来る。
は、取得期間において64KHzの抽出速度で、アナログ
−デジタル変換器218によって、連続的に抽出され
て、デジタル化される。一組の256個のデジタル数字
が、核磁気共鳴信号の各IとQの成分に関して、同時に
得られ、これらのデジタル数字は、直列リンク105を
経て主コンピュータ101へ送られる。各一対の数字
は、I+iQに等しい複素数値と考えられる。図4の核
磁気共鳴装置は、一連のパルス・シーケンスを行い、上
述のように、所望の画像復元に十分な核磁気共鳴データ
を収集する。個々の心筋サンプルと、心臓周期における
複数フレームのこれらサンプル内のひずみとを反復して
追跡することによる位相差磁気共鳴画像法にもとずい
て、心筋運動を評価する、改良された方法と装置が説明
された。また、張力成長の速度に依存する画像形成法も
説明された。本発明は、ほかの器官及び筋肉の運動の研
究に、また、非生物学的サンプルの研究に適用すること
が出来る。
【0040】このように、本発明は、特定の実施例に関
して説明されたが、説明は本発明を示しており、本発明
を限定するものとして解釈されるべきではない。多様な
修正と適用は、添付請求の範囲により定義されているよ
うに、本発明の本来の精神と範囲から逸脱することな
く、本技術分野の精通者には、想起されるものである。
して説明されたが、説明は本発明を示しており、本発明
を限定するものとして解釈されるべきではない。多様な
修正と適用は、添付請求の範囲により定義されているよ
うに、本発明の本来の精神と範囲から逸脱することな
く、本技術分野の精通者には、想起されるものである。
【図1】図1Aは速度導関数の関数として剪断張力を示
す。図1Bは速度導関数の関数として剪断張力を示す。
す。図1Bは速度導関数の関数として剪断張力を示す。
【図2】本発明の実施例による、器官の組織の運動を計
算する方法の流れ図である。
算する方法の流れ図である。
【図3】本発明による、フレーム間の関心領域の移動
と、移動の、測定と計算の変数を示す。
と、移動の、測定と計算の変数を示す。
【図4】本発明を使用している核磁気共鳴装置の構成図
である。
である。
【図5】図4の核磁気共鳴装置の一部を形成しているト
ランシーバーの電気的構成図である。
ランシーバーの電気的構成図である。
100 コンピュータ装置 101 主コンピュータ 102 インタフェース 118 システム制御器 120 パルス制御モジュール 122 トランシーバー 124 ステイタス制御モジュール 128 勾配増幅装置 136 勾配コイル組立体 138 無線周波送受信コイル 142 主磁石電源 144 主磁石組立体 152 患者位置づけ装置
─────────────────────────────────────────────────────
【手続補正書】
【提出日】平成5年8月30日
【手続補正1】
【補正対象書類名】図面
【補正対象項目名】図2
【補正方法】変更
【補正内容】
【図2】
【手続補正2】
【補正対象書類名】図面
【補正対象項目名】図4
【補正方法】変更
【補正内容】
【図4】
【手続補正3】
【補正対象書類名】図面
【補正対象項目名】図5
【補正方法】変更
【補正内容】
【図5】
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.5 識別記号 庁内整理番号 FI 技術表示箇所 G06F 15/70 405 8837−5L H04N 7/18 U
Claims (7)
- 【請求項1】 対象物内の関心領域の運動と変形とを評
価する装置にして、 (a) 前記対象物内に磁界を形成する第1の手段と、 (b) 複数の時間フレームにおいて前記領域の少なく
とも一つの速度成分を表す磁気共鳴画像信号を得る第2
の手段と、 (c) 一つの時間フレームにおいて前記画像内の関心
領域の位置と形状とを識別する第3の手段と、 (d) 前記一つの時間フレームにおいて前記関心領域
内の速度成分から誘導された変数を計算する第4の手段
と、 (e) 前記速度変数にもとずいた前記の一つの時間フ
レーム以外の時間フレームにおいて前記関心領域の新し
い位置と新しい形状を推定する第5の手段と、 (f) 前記関心領域を前記の新しい位置へ変換する第
6の手段と、 (g) 前記の新しい位置における前記関心領域を変形
する第7の手段と、 (h) 前記関心領域の運動と変形を決定するために、
次の時間フレームに関して段階(c)から段階(g)ま
でを反復する第8の手段とより成ることを特徴とする前
記の装置。 - 【請求項2】 前記第4の手段が速度Vx と速度Vy を Vx =a+b(x−xav)+C(y−yav) Vy =d+e(x−xav)+f(y−yav) としてモデル化し、この場合、a、b、c、d、e、f
は、各フレームの前記関心領域に関して最小二乗法によ
り計算された定数であることを特徴とする請求項1に記
載の装置。 - 【請求項3】 前記第4の手段が各フレームにおいて縦
方向の張力の変化を Δεx,xi=bΔt Δεy,yi=fΔt とし、変位の変化を Δxi=aΔt Δyi=dΔt とし、剪断の変化を Δθi =〔(eΔt)〕+(cΔt)〕とし、回転の変
化を Δα=1/2〔(eΔt)〕−(cΔt)〕として、計
算することを特徴とする請求項2に記載の装置。 - 【請求項4】 前記第5の手段が、前記関心領域の組合
せられた新しい位置と新しい形状とを推定する場合に、
前記関心領域の座標系をフレームの間で回転することを
特徴とする請求項3に記載の装置。 - 【請求項5】 前記関心領域の前記運動と変形の映写画
像を形成する第9の手段よりなることを特徴とする請求
項4に記載の装置。 - 【請求項6】 対象物内の関心領域の変形を評価する装
置にして、 (a) 磁界を前記対象物の全体に形成する第1の手段
と、 (b) 複数のフレームにおける前記領域の複数の速度
成分を表す磁気共鳴画像信号を得る第2の手段と、 (c) 前記領域内の複数の点における前記速度成分に
関する変数を計算する第3の手段と、 (d) 速度の空間導関数にもとずいた前記複数の時間
フレームにおける前記関心領域の変形を評価する第4の
手段とより成ることを特徴とする前記の装置。 - 【請求項7】 前記第4の手段が、前記領域内の張力の
成長速度に関する画像を形成する手段を有することを特
徴とする請求項6に記載の装置。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US07/617,904 US5195525A (en) | 1990-11-26 | 1990-11-26 | Noninvasive myocardial motion analysis using phase contrast mri maps of myocardial velocity |
US617904 | 1990-11-26 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH06142079A true JPH06142079A (ja) | 1994-05-24 |
JP3118278B2 JP3118278B2 (ja) | 2000-12-18 |
Family
ID=24475527
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP03236526A Expired - Fee Related JP3118278B2 (ja) | 1990-11-26 | 1991-09-17 | 心筋速度の位相差磁気共鳴画像診断法による非侵襲的心筋運動分析法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US5195525A (ja) |
EP (1) | EP0488496B1 (ja) |
JP (1) | JP3118278B2 (ja) |
DE (1) | DE69130975T2 (ja) |
ES (1) | ES2131505T3 (ja) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002522146A (ja) * | 1998-08-10 | 2002-07-23 | ザ ジョンズ ホプキンス ユニバーシティ | 角度画像を用いてタグ付き磁気共鳴イメージングにおける対象物の動きを測定する方法 |
JP2003506174A (ja) * | 1999-08-05 | 2003-02-18 | アメリカ合衆国 | 磁気共鳴映像法において位相ラベル付けにより目標の内部の運動および全体的な運動の写像を得る方法および装置 |
JP2013223726A (ja) * | 2012-04-19 | 2013-10-31 | General Electric Co <Ge> | 磁気共鳴血管造影のシステムおよび方法 |
JP2015511527A (ja) * | 2012-03-30 | 2015-04-20 | 日立アロカメディカル株式会社 | 超音波イメージング方法及び超音波イメージング装置 |
JP2015511526A (ja) * | 2012-03-30 | 2015-04-20 | 日立アロカメディカル株式会社 | 超音波イメージング方法及び超音波イメージング装置 |
Families Citing this family (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3920997C1 (ja) * | 1989-06-27 | 1990-08-23 | Bruker Analytische Messtechnik Gmbh, 7512 Rheinstetten, De | |
US5348011A (en) * | 1991-11-14 | 1994-09-20 | Picker International, Inc. | Shared excitation phase encode grouping for improved throughput cardiac gated MRI cine imaging |
US5329925A (en) * | 1991-11-14 | 1994-07-19 | Picker International, Inc. | Reduced scan time cardiac gated magnetic resonance cine and flow imaging |
US5352979A (en) * | 1992-08-07 | 1994-10-04 | Conturo Thomas E | Magnetic resonance imaging with contrast enhanced phase angle reconstruction |
DE4310993A1 (de) * | 1993-04-03 | 1994-10-06 | Philips Patentverwaltung | MR-Abbildungsverfahren und Anordnung zur Durchführung des Verfahrens |
US5318026A (en) * | 1993-08-11 | 1994-06-07 | Board Of Trustees Of The Leland Stanford Junior University | Method and apparatus for tracking of deformable regions by phase contrast MRI |
DE4329922A1 (de) * | 1993-09-04 | 1995-03-09 | Philips Patentverwaltung | MR-Abbildungsverfahren und Anordnung zur Durchführung des Verfahrens |
US5517992A (en) * | 1994-12-27 | 1996-05-21 | General Electric Company | Magnetic resonance imaging of shear-rate within mechanically vibrated materials |
US5680862A (en) * | 1995-02-01 | 1997-10-28 | The Board Of Trustees Of The Leland Stanford Junior University | Iterative method of determining trajectory of a moving region in a moving material using velocity measurements in a fixed frame of reference |
AU7837898A (en) | 1997-06-12 | 1998-12-30 | Johns Hopkins University School Of Medicine, The | Method for measuring myocardial motion and the like |
US7572263B2 (en) * | 1998-04-01 | 2009-08-11 | Arthrocare Corporation | High pressure applicator |
WO1999049819A1 (en) * | 1998-04-01 | 1999-10-07 | Parallax Medical, Inc. | Pressure applicator for hard tissue implant placement |
US6240310B1 (en) * | 1998-06-26 | 2001-05-29 | Siemens Medical Systems, Inc. | Method of acquiring MR data especially for cardiac cine data sets |
US20050240445A1 (en) * | 1998-09-29 | 2005-10-27 | Michael Sutherland | Medical archive library and method |
US6954802B2 (en) | 1998-09-29 | 2005-10-11 | Tdk Electronics Corporation | Removable media recording station for the medical industry |
US7050610B2 (en) * | 2001-04-04 | 2006-05-23 | Siemens Medical Solutions Usa, Inc. | Method and system for improving the spatial resolution for strain imaging |
US6438196B1 (en) | 2001-06-28 | 2002-08-20 | General Electric Company | EKG driven CT image reconstruction for cardiac imaging |
US7047060B1 (en) | 2001-11-26 | 2006-05-16 | Koninklijke Philips Electronics N.V. | Multiple preparatory excitations and readouts distributed over the cardiac cycle |
WO2003071950A1 (en) * | 2002-02-27 | 2003-09-04 | Amid Srl | M-tracking for space-time imaging |
JP4060615B2 (ja) * | 2002-03-05 | 2008-03-12 | 株式会社東芝 | 画像処理装置及び超音波診断装置 |
US6994673B2 (en) * | 2003-01-16 | 2006-02-07 | Ge Ultrasound Israel, Ltd | Method and apparatus for quantitative myocardial assessment |
US8064979B2 (en) * | 2003-06-09 | 2011-11-22 | General Electric Company | Tempero-spatial physiological signal detection method and apparatus |
AU2004251360A1 (en) * | 2003-06-25 | 2005-01-06 | Siemens Medical Solutions Usa, Inc. | Automated regional myocardial assessment for cardiac imaging |
US7653227B2 (en) * | 2004-02-09 | 2010-01-26 | Siemens Medical Solutions Usa, Inc. | Hierarchical modeling in medical abnormality detection |
US7338436B1 (en) | 2005-06-27 | 2008-03-04 | Pacesetter, Inc. | Implantable cardiac physical motion mapping device |
EP1960809A2 (en) * | 2005-11-29 | 2008-08-27 | Ullevål Universitetssykehus HF | Choosing variables in tissue velocity imaging |
US8170307B2 (en) | 2008-09-23 | 2012-05-01 | The Methodist Hospital | Automated wall motion quantification in aortic dissections |
DE102014209803B4 (de) * | 2014-05-22 | 2016-03-03 | Siemens Aktiengesellschaft | 4D-Geschwindigkeitsprofil-Bildaufnahme mit einer Magnetresonanzanlage |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4710717A (en) * | 1986-12-29 | 1987-12-01 | General Electric Company | Method for fast scan cine NMR imaging |
JPH064066B2 (ja) * | 1987-10-15 | 1994-01-19 | 株式会社東芝 | 磁気共鳴イメージング装置 |
JP2646663B2 (ja) * | 1988-06-07 | 1997-08-27 | 株式会社日立製作所 | 動体イメージング方法およびその装置 |
US4995394A (en) * | 1989-08-11 | 1991-02-26 | General Electric Company | Fast NMR cardiac profile imaging |
-
1990
- 1990-11-26 US US07/617,904 patent/US5195525A/en not_active Expired - Lifetime
-
1991
- 1991-08-01 ES ES91307075T patent/ES2131505T3/es not_active Expired - Lifetime
- 1991-08-01 EP EP91307075A patent/EP0488496B1/en not_active Expired - Lifetime
- 1991-08-01 DE DE69130975T patent/DE69130975T2/de not_active Expired - Lifetime
- 1991-09-17 JP JP03236526A patent/JP3118278B2/ja not_active Expired - Fee Related
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002522146A (ja) * | 1998-08-10 | 2002-07-23 | ザ ジョンズ ホプキンス ユニバーシティ | 角度画像を用いてタグ付き磁気共鳴イメージングにおける対象物の動きを測定する方法 |
JP2003506174A (ja) * | 1999-08-05 | 2003-02-18 | アメリカ合衆国 | 磁気共鳴映像法において位相ラベル付けにより目標の内部の運動および全体的な運動の写像を得る方法および装置 |
JP2015511527A (ja) * | 2012-03-30 | 2015-04-20 | 日立アロカメディカル株式会社 | 超音波イメージング方法及び超音波イメージング装置 |
JP2015511526A (ja) * | 2012-03-30 | 2015-04-20 | 日立アロカメディカル株式会社 | 超音波イメージング方法及び超音波イメージング装置 |
JP2013223726A (ja) * | 2012-04-19 | 2013-10-31 | General Electric Co <Ge> | 磁気共鳴血管造影のシステムおよび方法 |
Also Published As
Publication number | Publication date |
---|---|
DE69130975D1 (de) | 1999-04-15 |
US5195525A (en) | 1993-03-23 |
EP0488496B1 (en) | 1999-03-10 |
EP0488496A3 (en) | 1992-10-14 |
DE69130975T2 (de) | 1999-07-08 |
EP0488496A2 (en) | 1992-06-03 |
JP3118278B2 (ja) | 2000-12-18 |
ES2131505T3 (es) | 1999-08-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3118278B2 (ja) | 心筋速度の位相差磁気共鳴画像診断法による非侵襲的心筋運動分析法 | |
US7949170B2 (en) | Image processing method, image processing device, computer aided detection, and method for filtering along the time axis | |
JP5198859B2 (ja) | 位相マッピングと、位相基準として用いる基準媒体が関係するmri温度測定 | |
US5594849A (en) | Biomedical magnetism imaging apparatus and method | |
JP4379559B2 (ja) | 磁気共鳴イメージング装置 | |
JPH11321A (ja) | 速度符号化核磁気共鳴画像から変形を示す画像を発生する方法 | |
JP2003528663A (ja) | マイクロコイルを用いた磁気共鳴撮像方法及びシステム | |
JPWO2002022012A1 (ja) | 磁気共鳴イメージング装置 | |
JPWO2014185521A1 (ja) | 磁気共鳴イメージング装置、画像処理装置、画像診断装置、画像解析装置、mri画像作成方法およびプログラム | |
US6317620B1 (en) | Method and apparatus for rapid assessment of stenosis severity | |
JPS6188132A (ja) | 核磁気共鳴の手段による流動的画像形成 | |
US10362961B2 (en) | System and method for neutral contrast magnetic resonance imaging of calcifications | |
JPH03184530A (ja) | Nmr装置 | |
JP5465565B2 (ja) | 磁気共鳴イメージング装置 | |
US6507749B1 (en) | Method and apparatus for tracking the motion of fluid and determining a velocity spectrum thereof from MR data acquired in a single cycle | |
JP3501182B2 (ja) | 流速画像を算出できる磁気共鳴イメージング装置 | |
US20210311146A1 (en) | Magnetic resonance fingerprinting thermometry | |
JP3869056B2 (ja) | 医用画像診断装置における画像表示方法 | |
US11754654B2 (en) | Method and device for determining a motion field from k-space data | |
HU213190B (en) | Method and magnetic resonance based apparatus for determining movement inside of a part of examined person without surgical operation | |
JPH08317915A (ja) | 磁気共鳴イメージング装置 | |
JP4118119B2 (ja) | 磁気共鳴イメージング装置 | |
KR100222143B1 (ko) | 공간의 자화변조를 이용한 심장벽운동 분석장치 | |
JPH05123314A (ja) | 磁気共鳴イメージング装置におけるマルチスライス撮像方法 | |
JPH01274749A (ja) | 核磁気共鳴イメージング装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
LAPS | Cancellation because of no payment of annual fees |