JPS63125250A - Ｎｍｒデータの集約方法及び装置 - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 背景技術 （技術分野） 本発明は一般的に核磁気共振（ｎｕｃｌｅａｒ　ｍａｇ
ｎｅｔｉｃｒｅｓｏｎａｎｃｅ　　；　ＮＭＲと略称）
トモグラフィに関するものであり、とくに心臓位相・呼
吸位相面内のＮＭＲｋ−スペース測定値の時間遡及クラ
スタリング（集約化）及び分析に関するものである。本
プロセスを行う装置（システム）及び方法により、心臓
血管の運動及び呼気運動の映像内のアーチファクト（誤
差陰影）を減少でき、映像を用いるかまたは用いずして
心臓構造及び機能の量子化特性判断ができる。

（背景技術） ＮＭＲ分光分析（スペクトロスコピイ）、すなわち均−
磁界内の小サンプルの分析プロセス及び正確なパルス化
無線周波数励振による無線周波データを得る技術はブロ
ック及びパーセル（Ｂｌｏｃｋ及びＰｕｒｃｅｌｌ　）
によって発明された。過去１６年間にわたる研究中に、
分光分析（スペクトロスコピイ）によるＮＭＲ分析は、
物理化学の領域より生化学（バイオロジカルケミストリ
ー）及びバイオロジカル医学応用分野、すなわち健康組
織と病的組織を生化学的にサンプルする領域に移行して
きた。

ローチルバール（Ｌａｕｔｅｒｂｕｒ　）及びダマディ
アン（Ｄａｍａｄ　ｉａｎ　）及び他の人がそれぞれ別
個にＮＭＲ原理を用いて映像を形成することを発明して
いる。

“参考文献” アールΦダマディアン（ＲｌＤａｍａｄｉａｎ）サイエ
ンス１７１．　１１５１．１９７１　；ビー會シー・ロ
ーチルバール（Ｐ、Ｃ１Ｌａｕｔｅｒｂｕｒ）ネーチュ
ア２４２．１９０．　’１９７３　　；ピーψシーｅロ
ーチルバール ピュア・アンド・アプライド・ケミストリー４０、１４
９．１９７４ これらによって得られた装置、すなわちＮＭＲ映像シス
テムは二次元及び三次元データを生じ、その内に表され
るグレースケールは複数のパラメータによる機能を表わ
す。これらのパラメータとは、例えば、とくに解剖学的
映像内の３つのパラメー夕、すなわち核密度にュークラ
イド密度）と、Ｔ、で表される縦方向リラクセーション
時間と、Ｔ２で表わされる横方向リラクセーション時間
を含むものである。

ＮＭＲ映像技術は例えば“プロトンＮＭＲ）モグラフィ
”ビー・アール・ローヘル（Ｐ、Ｒ，Ｌｏｃｈｅｒ）著
フィリップス・テクニカル・レビューＶｏｌ、　４１゜
１９８３／８４　Ｎα３．７３−８８頁に開示され発表
されており、その内容を本明細書中ＮＭＲ映像技術の背
景技術とする。

生物組織のＮＭＲ生態（ｖｉｖｏ）映像では生物組織の
移動によって、かなりの難しさがある。とくに心臓のＮ
ＭＲ映像は、心臓の鼓動によりとくに一層の困難さを生
ずる。心臓のＮＭＲ映像内のぼけ（ｂｌｕｒｒｉｎｇ）
と移動によるアーチファクト（ａｒｔｉｆａｃｔｓ　　
；誤差陰影）をなくすため、測定を患者のＥ、　Ｃ０Ｇ
、に同期させることが知られている。

これは例えば次の如き発表がある。

ホラセル（ＦＯ８Ｓｅｌ）他米国特許（口、Ｓ、　Ｐ、
　）４．４０９．５５０ホツセル（ＦＯ８Ｓｅｌ）他米
国特許（Ｕ、Ｓ、　Ｐ、　）４．４１３．２３３核磁気
共振（ＮＭＲ）映像の非移動的特性と、骨格構造による
妨害の無いことより、この映像再現方法は心臓の映像化
に極めて好ましい技術である。しかし比較的に長い走査
時間を必要とすることから動きによるアーチファクトが
生ずる。心臓の鼓動サイクルに映像のシーケンスを同期
させることにより、このようなアーチファクトを大幅に
減少させ得る。

ＮＭＲ法では、心臓の鼓動サイクルの任意の位相（フェ
ース）を映像化することが可能であるため、収縮期末よ
り拡張期末までの量的測定を可能にする。さらにある映
像のループにおける連続する位相の映像を表示すること
によって動作を評価することも可能である。心臓と大動
脈の解剖学的状況を表示するのに、非トリがまたは電気
カージオグラム（ＢＣＧ））’Ｊガの両者によりＮＭＲ
映像を使用すること、及びその利点は良く知られている
。この技術において、心臓内の解剖上の細部が示される
他、ＮＭＲ心臓映像法は、梗塞の検出と寸法（大きさ）
測定に重要な、組織の特性検査を可能にする。

速度（ベロシティ）映像を用いて、心臓壁の運動及び血
流速度がＮＭＲ映像法によって量子化される。

二次元フーリエ変換（２ＤＦＴ）映像法を用いると、複
合２ＤＦＴによる時系列信号よりＮＭＲ映像を再構成（
リコンストラクション）できる。１２８　Ｘ１２８画素
（ピクセル）の映像マトリックスに対し、１２８の連続
映像シーケンスの信号が得られる。心臓の映像化には、
患者の６ＣＧのＲ波により導出したパルスによってこれ
らの各シーケンスをトリガする。

このパルスの遅延時間により心臓鼓動サイクル内の映像
化位相が定まる。各シーケンスは、一連の無線周波数（
ＲＦ）及び磁界グレーディエンド（傾斜）パルスにより
それぞれ形成され、信号の識別と信号内の空間情報を提
供する。ＲＦ励振パルス（９０゜パルス）の後に供給さ
れるグラディエンド磁界は、異なる周波数で励振スライ
スのブレセス（歳差運動前進）内の異なる位置のスピン
を与える。このようなスピンは共振周波数位相に対し位
相外れを開始させる。所定の位置に対して生ずる位相外
れは、グラディエンドの振幅とグラディエンドの加わっ
ている時間に比例する。さらにこれに加えて、固定スピ
ンと移動スピンの間に位相シフトの差が存する。グレー
ディエンドの方向内に均一に移動するスピンに対し、グ
レーディエンド方向の速度に正比例する特別項が加算さ
れる。

（発明の目的） 本発明は、従来より心臓トリガリング（カーブイアツク
トリガリング）及び呼吸ゲーティング（レスピラトリー
ゲーティング）に一般に使用されている心臓／呼吸モニ
タ装置よりのＮＭＲデータを相関させる改良した方法と
システムを得ることを目的とする。

（従来技術） 従来もこの問題を解決するため多くの試みがなされてき
ていた。ＢＣＧ　Ｒ・波形の後に一定の遅延を設けてＮ
ＭＲデータ収集を開始すると異常なＲ−Ｒ間隔を排除で
きる。典型的な既知の回路では、Ｒ・ピークの後にある
一定の固定した遅延を設けて各データ収集を開始させる
。多スライスによる研究では、データ収集は次のＲ・ピ
ークに接近して継続される。心臓が正常に規則正しく鼓
動している場合には、各測定瞬時は心臓サイクル内の同
じ位相に位置する。

しかしながら、不整脈が存する場合には、いくつかのデ
ータ収集に対し心臓位置のずれが生じ、これによってぼ
けと移動によるアーチファクトが生ずる。またこのよう
な不整脈は反復時間ＴＲに変化を生じる結果をもたらし
易く、このためＴ、測定の正確度の低下と雑音レベルの
増加を由来し、また−次映像の全体的映像品質を低下さ
せる。これについては次の発表がある。

シー・ガロナト デー・ジェー・トロスト ニス・ニス・プラトー −、；−・ｙイｔｆンハークＳＭＲＭＶｏｌ、２．１９
８５゜核医療分野においては、不整脈中、およびその直
後において、収集されたデータを排除（リジェクト）す
ることにより映像品質と正確さとが改善されることが良
く知られている。しかしながらこれ迄においては、この
ような方法のデータ収集を簡単かつ効果的に行うことは
困難であった。

呼吸サイクルのタイミングを用いたＮＭＲｋ・スペース
サンプリングも既に利用されている。しかしながらこれ
らでは、異常な生理学的サイクルを排除（ディスカード
）したり、生理学的サイクルに連係したｋ・スペースサ
ンプルを行うと、実際上の異常性を見落とすこととなり
、ＮＭＲ走査による診療の目的に対しては欠陥となって
いた。冠状動脈結滞のようなときはＲ−Ｒ間隔の変化に
よって、心臓サイクルの終わりにおいては、心臓トリガ
は有効性が少なくなり、心臓の前（プレ）・Ｒ・波形が
心臓収縮異常に対する検出の端緒となることがある。

従来考えられていたＮＭＲデータ収集の解決手段は、ス
ライス、位相、トリガ、遅延、ウィンドウ（時間枠）、
データ排除、記録等を組合せた複雑なものであり、心臓
位相の不正確さ、解像時間の限定、スピン飽和の不均一
等の問題を生じていた。

（発明の要旨） 本発明はこれらの各難点を解決するためになされたもの
である。

本発明によるタイムクラスタ心臓・呼吸エンコータハ、
マイクロコンビ二一夕をベースとスル装置（システム）
であり、心臓（カーブイアツク；Ｃ）Ｒ・波時間、呼吸
（レスビラトリ；Ｒ）ダイヤフラム位置、及びＮＭＲデ
ータ取得（アクイジション；Ａ）タイミングデータを並
列に収集するもので、標準磁気共振イメジャー（映像器
）上に自走（フリーランニング）像として収集する。多
数の映像に等価の原データが収集された後、Ｃ−Ａ−Ｒ
位相タイミング　データをマイクロ　コンピュータより
そのスキャナのミニコンピユータにアップロードし、こ
れらの原データを新規なデータセットとして使用できる
ように集約化（クラスタ）を行い、所望の各心臓・呼吸
（Ｃ−Ｒ）位相の組合せに１つの映像が等価となるよう
にする。量子化したＣ−Ｒ位相平面内のＮＭＲプロフィ
ルの位置が、クラスター映像を決定し、これは例えば、
何れかの心臓位相と呼吸位相の組合せに属することとな
る。次でかくクラスタされたデータを記入（フィル）し
、かつふるい分け（フィルタ）して、非均一ｋ・スペー
スサンプリングを補償し、最終的には一般の二次元（２
Ｄ）または三次元（３０）フーリエ変換によって再構成
する。

２つの複式（マルチプル）プロフィル映像を得るため、
すなわち心臓のサイクルの２つの異なった位相における
心臓の状態を表示するため、本発明を使用すると、例え
ば心臓サイクルの最初の５分の１の周期期間を示すもの
では、この期間中のそれぞれのＲ−Ｒ間隔にわたって求
めたすべてのプロフィルよりのタイム・クラスタ映像を
得、また心臓サイクルの第４番目の５分の１周期期間を
示すものではこの第４５分の１周期期間中の対応のＲ−
Ｒ間隔中に求めたプロフィルによる心臓サイクルのタイ
ム・クラスタ映像を得ることができる。

本発明の装置（システム）は任意のＮＭＲ装置と組合せ
て構成でき、かつ任意のパルスシーケンスで使用できる
。本装置は、従来の心臓トリガ及び呼吸ゲートに共通に
使用されている標準の心臓／呼吸モニタ装置より容易に
得られる出力信号を必要とするのみで、またこれに伴っ
てスキャナよりのタイミング信号と異なるパルスシーケ
ンスに適合するようにソフトウェアのパラメータを変化
させることを必要とするのみである。

本発明は心臓位相・呼吸位相面におけるＮＭＲｋ・スペ
ース測定プロフィルの時間遡及ならびに集約（クラスタ
リング）を行う簡単かつフレキシブルなシステムを提供
しようとするもので、とくに次の点を考慮しているもの
である。

１）　心筋及び呼吸フロー／運動映像アーチファクトの
減少、 ２）　映像化を行うか、行わないかの双方に対し、心臓
の構造／機能の量子化特性づけ。

（実施例） 本発明は心臓位相、呼吸位相平面内のＮＭＲｋ・スペー
ス測定プロフィルの時間遡及（レトロスペクティブ）集
約化（クラスタ）に対するエンコーダ方法及び装置に関
するものであり、この方法と装置とは次の如くの機能を
行う。

１）　心筋及び呼吸フロー運動の映像アーチファクトを
減少させること、及び ２）　映像を映写するか、或いは映写を行わずして心臓
の構造／機能の量子化特性測定を行うこと。

第１図は本発明によるタイムクラスタ心臓、呼吸（カー
ディオレスピラトリー）エンコーダ（ＴＣＣＲＢ　）の
全体図を示し、第２図はこれに使用するハードウェアで
あり、第３図はこのエンコーダ及びその機能的使用方法
の説明図である。これら各図面において心ＩＩ（カーデ
ィアク（Ｃ））時間、呼吸（レスピラト’Ｊ−（Ｒ））
時間、またはダイアフラム位置、及びＮＭＲデータ収集
（アクイジション（Ａ））サイクルマーカデータをマイ
クロコンピュータ（１）によって収集し、これと並列に
標準磁気共振映像器（２）用にフリーランニングの特に
安定状態とするを可とする映像収集を行う。複数の映像
に対応する原データを収集した後、Ｃ−Ａ−Ｒ位相マー
カデータをマイクロコンピュータ（１）よりＮＭＲスキ
ャナーのミニコンピユータ（３）ニ。

−ドし、映像の原データのＣ−Ｒ平面クラスタに使用す
る。量子化したＣ−Ｒ位相平面内のＮＭＲプロフィルの
位置を決定し、次いで新しいデータのセットとしてクラ
スタを行う。この１つの映像はそれぞれ各々の所望のＣ
−Ｒ位相組合せに等価なものである。次いでこれらの原
データを記入し、またふるい分け（フィルタ）を加え、
均一でないｋ・スペースサンプリングに対する補償を行
い、最終的に再構成を行う。

第１図及び第３図に示した本発明の方法は次の如くであ
る。

（１）フリーランニングの、また可能な場合は速度符号
化パルスシーケンスよりＮＭＲ７’ロフイルを収集し、
同時にディジタル化εＣＧ信号を収集する。

Ｐ（ｋ　　ｘ、　　ｋ　　ｙ、　　ｋ　　ｚ、　　ｔ）
ＢＣＧ　（ｔ）→ＢＣＧ−ｊ（ｔ）　；心臓サイクルの
ｊ分の１、（２）直列核医療データに使用すると同じよ
うな技術を用い、利用可能な信号対雑音に適当な量子化
された心臓サイクル時間枠（ビン）を検索する。

即ち、プロフィルの時間を使用し、これより前位のＲ・
波の時間をマイナスし、これらを後続するＲ・波時間で
割り、さらに前位のＲ・波時間をマイナスし、プロフィ
ルを設けるべき時間枠を決定する。呼吸サイクル内のプ
ロフィル位置は同様にして次で定まる。

Ｐ　　ｉ（ｋ　　ｘ、ｋ　　ｙ、ｋ　　ｚ、　　ｔ）；
心臓サイクルの１分の１の位相 ｔｉ−１＜ｔ＜ｔｉ これらのデータは（速度）映像のセットを形成し、これ
らデータの１つは移動している心臓組織に対し選択した
それぞれの心臓サイクル位相の時間枠に対応するが、こ
のセットの核映像はｋ・スペースを不完全にカバーし、
また少し異なるようにこれを表わす。

（３）　　周波数スペースを保ちながら、データをパラ
メータ付のダイナミックな心臓モデルにフィツトさせる
。

（４）ディスクリートな、均一でないフーリエ変換を用
い、非均一なｋ・スペースカバーを克服するため、この
データを記入（充填）し、ふるい分けし、変換する。

（５）所望に応じ、心臓のダイナミックシュミレーショ
ン（Ｄ）をそのパラメータと共に集積グラフィックワー
クステーション（図表表示装置）　（Ｗ）にディスプレ
イ　（表示）を行う。

（６）所望に応じ、新規の再構成データを用い、ダイナ
ミック心臓モデルをアップディトする。

ＴＣＣＲＢ映像はより少ないフロー／モーションアーテ
ィファクトを示し、より良好な心臓の表示を示す。これ
は、比較位相トリガ映像に比べて良好であり、特に冠状
動脈閉塞が生じている箇所の心臓サイクルの限界的（ク
リティカル）な後側の部分においてこれが顕著である。

第１図及び第３図に示す如＜　、ＴＣＣＲＢはディジタ
ル化した心臓Ｒ・波信号と、アナログまたはディジタル
呼吸入力及びＮＭＲデータ収集時間データをそれぞれ並
列に受信する。

第３図に示すように、このＴＣＣＲＥ　　（タイムクラ
スタートカーディオレスビラトリーエンコーダ）は１つ
の走査中に収集できるか、または１つの観測期間で集め
られる全てのｌｔＭＲ及び生理学的（フィジオロジカル
（ＰＤ））データ、（例えば心臓εＣＧ及び呼吸信号等
）を数学的心臓構造／機能モデル（Ｍ）並びに同じ患者
の心臓の以前に調べたタイムクラスタート心臓、呼吸エ
ンコーダスキャン”を用いて集約化（インチグレート）
する。この結果は映像化を行うか、或いは行わない双方
に対し、可能な限り完全な心臓の特性化したものとなる
。

全てのデータは可能な限り急速に、かつ均一となるよう
にこれらを収集し、磁界不均一の修正とデータの修正を
容易に行わせるためのスピン飽和及び渦流の平衡化を生
じさせる。この集約化されたデータは、映像品質を改良
するために利用でき、また所望の構造／機能バラメーク
を計算するために利用できる。このシステムの出力は、
心臓のダイナミック三次元シュミレーション（Ｄ）、或
いは最大の可能性をもってフィツトされた心臓構造機°
能（ＣＦ）特性化パラメータ（Ｐ）または映像（ＩＭ）
であり、これらは呼吸フロー運動アーチファクトを減少
させ、改良したものである。

第１図及び第３図にふいて使用したＮＭＲ及び生理学的
タイミング信号（ＰＤ）はＮＭＲパルスプログラマタイ
ミング信号（Ａ）、心臓ＥＣＧ　Ｒ・波（Ｃ）及びアナ
ログダイアフラム位置呼吸信号（Ｒ）である。

これらの全ての信号は測定中の全ての時間においてマイ
クロプロセッサ（１）によって収集される。

これらのデータはミニコンピユータ（３）に移され、こ
のミニコンピユータ内でソフトウェアにより量子化した
心臓及び呼吸サイクルの基準フレーム、即ちＣ−Ｒ平面
内にＮＭＲＫ・スペースプロフィルデータをクラスタす
る。必要に応じ、再構成前または再構成中に数学的モデ
ル（Ｍ）をこのデータに適応させることができ、このモ
デル（Ｍ）は予め存在しているデータにより変形するこ
とができる。

このようにプロセスを加えたデータは映像（ＩＭ）に再
構成でき、これをディスプレイ表示しうるか、この構造
／機能パラメータ（Ｐ）は脈動している心臓のシュミレ
ーションを駆動するに使用することができる。これらを
第１図及び第３図に図表的に示しである。

このシステムは連続的に集められたＮＭＲデータとＥＣ
Ｇ　Ｒ・波信号とディジタル化した呼吸信号を映像品質
改良の目的によりノルマル化（量子化）したＣ−Ｒ平面
内に時間遡及的クラスタを行うために使用できる。

第４Ａ図に示した図表は２つの６４プロフイル映像の収
集を示すものであり、心臓のサイクルの２つの異なる位
相を示し、呼吸クラスタリングを伴っていないものを示
している。この図表を参照するとＴＣＣＲＢによりクラ
スタされた映像Ａはサイクルの最初の５分の１の期間に
対応し、これらはプロフィル；ａ２．　ａ７．　・・・
、ｅ５８により形成されるものである。ＴＣＣＲＢクラ
スタ映像Ｂはサイクルの４番目の５分の１に対応し、プ
ロフィルａ５．　ａ９．　ｅ６３及びｅ６４により形成
される。短いＲ２−Ｒ３間隔及び長いＲ６４−Ｒ６５間
隔は種々の心臓サイクル内で固定トリガによる遅延によ
る問題を示すものであり、プロフィルＢ２は、Ｒ−Ｒの
間隔終わりよりくるものであり、またプロフィル８６４
はその中央部よりくるものである。

呼吸のクラスタリングが望まれる場合には、各信号クラ
スタを同様にブレークダウンして量子化した呼吸サイク
ル軸に沿って細かいクラスタ表示とする。例えば第４Ｂ
図に示すような４×３心臓呼吸クラスタとする。

本発明方法を行うための装置を第２図に示しである。第
２図の装置（システム）は次の如き構成部分を有してい
る。

Ｍ６８０２０単一ボードコンピュータ（ＶＭＥブス／電
力供給ラック搭載ケース内）； ４−８チヤネルＡ［）Ｃカード（ＡＤＣ）　　；８チヤ
ネル光学的に分離されたディジタルＩ１０カード（Ｉｌ
ｏ）； 調整可能なラッチングを設けたディジタル化した入力ラ
ッチボー／　クス（ＬＡＴＣＨＢＯＸ）　；ディジタル
　エクイップメント　コーポレーション製　ＶＡＸ／１
１／７５０（３）　　’／　７　トウ！　７ヘースー　
７レイプロセツサ付； 映像再構成用フィリップス　フレキシブル　レコンスト
ラクション　パッケージ（ＦＬＥＸＩＲＩｌｉＣ）　　
；Ｍ６８０２０ベースグラフィック　ワークステーショ
ン（Ｗ）　　（３−Ｄ分析及び表示ソフトウェア付）。

【図面の簡単な説明】

第１図は、心臓診察システム内に設けた本発明によるタ
イムクラスタ心臓一呼吸エンコーダを示す機能的ブロッ
ク図、 第２図は、ＮＭＲ心臓診断システム内の本発明によるタ
イムクラスタ心臓一呼吸エンコーダのハードウェアブロ
ック図、 第３図は、本発明によるタイムクラスタ心臓一呼吸エン
コーダの使用する方法ステップの機能的ブロック図、 第４Ａ図は、タイムクラスタ映像と心臓トリガ映像の関
係を示すプロット、 第４Ｂ図は、平常な心臓及び呼吸サイクル平面内の２次
元的クラスタアレイを示す図である。 ■・・・マイクロコンピュータ ２・・・核磁気共振イメージヤ− ３・・・ミニコンピユータ ４・・・アレイプロセッサ ム− Φ ムー

Claims (1)

1. 【特許請求の範囲】 １、心臓及び呼吸サイクルのＮＭＲデータの集約（遡及
   タイム・クラスタ）装置において、 エレクロカーデイオグラム（ＢＣＧ）装置よりの心臓Ｒ
   ・波データ、呼吸ダイヤフラム位置データ、ＮＭＲゲー
   タ取得タイミングデータを含み、かつこれらに付随する
   フリーランニングの、安定状態を可とするＮ−Ｒイメー
   ジャーより収集映像データを受信し、互に関連させるデ
   ータ収集工程を行うデータ収集・プロセス手段と、 前記の収集データを新しいデータセットに 集約（クラスタ）し、各心臓・呼吸位相平面領域に等価
   である１つの映像と、量子化したＣ−Ｒ位相平面内のＮ
   ＭＲプロフィルの位置をこれらが存するクラスタ映像心
   臓位相・呼吸位相によって決定し、かつ既知の心臓位相
   、無制約時間分析及びスピン飽和に対するデータを形成
   する計算手段と、 前記収集データを非均一ｋ・スペースサン プリングを行うための記入及びふるい分けを行う手段と
   、 二次元（２０）または三次元（３０）フーリエ変換技術
   により前記データを再構成する手段とを具えてなること
   を特徴とするＮＭＲデータの集約装置。 ２、数学的構造／機能モデルを蓄積する手段と、前記集
   約（タイム・クラスタ）ＮＭＲデータを前記数学的構造
   ／機能モデルと共にプロセスし、ただし依然として周波
   数スペース内で心臓機能のダイナミックシュミレーショ
   ンまたは再構成映像を改善するためのパラメータを形成
   する手段をさらに含むことを特徴とする特許請求の範囲
   第１項記載の装置。 ３、前記再構成データにもとづいて、数学的構造／機能
   モデルをアップデート（更新）する特許請求の範囲第２
   項記載の装置。 ４、前記装置に接続され、前記再構成データに基づく映
   像をディスプレイ表示する図表ワークステーションを含
   む第１項、第２項または第３項に記載の装置。 ５、静止状態を可とするフリーランニング符号化パルス
   　シーケンス（列）よりＮＭＲプロフィルを収集する工
   程と、 同時にディジタル化データ信号を収集する 工程と、 前記信号プロフィルを量子化した心臓サイ クル時間枠内に分類（ソート）し、信号対雑音係数を補
   償する工程と、 運動している心臓組織に対し選択した各心 臓サイクル時間枠に１つづつの映像の新しいセットを形
   成し、前記セットの各映像がｋ・スペースを不完全かつ
   僅か異なった態様でカバーする如くする工程と、 前記データの記入とふるい分けを行う工程 と、 ディスクリートな不均一フーリエ変換を用 いて、前記不均一ｋ・スペースカバーを打消す工程とを
   含んでなることを特徴とする心臓位相・呼吸位相内のＮ
   ＭＲｋ・スペースプロフィルを集約（時間遡及クラスタ
   ）する方法。 ６、多くのパラメータを有するダイナミック心臓モデル
   を選択する工程と、 前記変換データを該ダイナミック心臓モデ ルに適合させて、再生映像を改良するか、あるいは心臓
   のダイナミックシュミレーションを得る工程を含んでな
   る特許請求の範囲第５項記載の方法。 ７、前記変換再構成データを用いて、前記ダイナミック
   心臓モデルをアップデートする特許請求の範囲第６項記
   載の方法。 ８、前記変換再構成データを用いて、心臓のダイナミッ
   ク三次元シュミレーションを表示する特許請求の範囲第
   ５、６または７項記載の方法。