JP2005518857A - 弾性撮影法を用いた物体の検査方法 - Google Patents

Abstract

本発明は、弾性撮影法に関し、その診断特性は、機械的振動による例示に対応する偏位時間の変化が分析され非線形性の変形が測定される点で向上されている。かかる変形は、病巣の非線形的な弾性特性の結果であり、情報の明らかな診断項目を構築する。

Description

本発明は、弾性撮影法を用いた物体の検査方法に関するとともに、かかる方法を実行する配置及びかかる配置を制御するコンピュータープログラムに関する。

特許文献1（ＰＨＤ９９−１５２）に開示のＭＲ弾性撮影法（ＭＲ＝磁気共鳴）は、物体のＭＲ画像における位相が物体に作用する機械的振動の影響下で変化するという事実を利用している。かかる変化の程度は、機械的振動の影響下の組織の偏位（ｅｘｃｕｒｓｉｏｎ）に依存する。組織の一定の機械的パラメーターに関連する情報は、従って、ＭＲ位相画像に由来し、つまり、核磁気位相を再生成する画像であって；例えば、かかるパラメーターは弾性に関連する。超音波弾性撮影法において、上記の偏位を直接測定することは可能である。

例えば、乳癌などの悪性病巣は、健常な組織に比較して弾性の面で明らかに異なっている。しかしながら、同様の弾性特性を有する良性の病巣も存在している。変則的な弾性特性は、組織が悪性であることを十分には示さない。
米国特許出願番号第０９／７４３６５９号明細書

本発明の目的は、異なるタイプの組織間で良好な区別を可能とする追加情報が送達される弾性撮影法を提供することである。

この目的は、本発明に従った物体の検査方法により達成され、それは：
ａ）機械的振動により物体を励起するステップ、
ｂ）物体における上記の振動により発生する波長に指向された物体のボクセルの残りの状態に由来する偏位時間の変化を測定するステップ、
ｃ）偏位時間における変化に由来する非線形的な変形を同定するステップ、及び
ｄ）非線形的な変形を評価するステップ、
を有している。

本発明は、以下の考慮及び洞察に基づいている。弾性係数は、（組織がより硬くなるにつれ）増加するので、生体組織に加えられた圧力が変化する際、この組織の変形は、圧力と同程度に変化しない。圧力と変形との間のこの非線形的な関係は、悪性病巣において特に顕著となり、その理由は、この組織が機械的振動により励起された際、偏位時間における変化は、もはや機械的振動の時間変化に対応せず、（非線形的な）変形（ｄｉｓｔｏｒｔｉｏｎ）を示す。従って、ボクセルの偏位の一時的な変化における非線形的な変形は、非線形性の結果である；なぜなら、測定は、組織の分類に最適であるためである。例えば経時的に正弦的に変化する機械的振動による励起の場合、かかる変形は、純粋に正弦的ではなく物体におけるボクセルの偏位を生じ、加えて、機械的振動のより高い調波（ｈａｒｍｏｎｉｃｓ）、特に周波数の二倍の調波を生じる。

各ボクセルは、各ボクセルに関して測定された非線形的な変形（及び可能性のある弾性）を同定すべく個々に測定されてもよい。しかしながら、請求項２に開示したバージョンは、検討者により大きな領域の診断を形成することを簡便にしている。

請求項３は、本発明がＭＲ弾性撮影法手段により実行されることを可能にしており、従って、非線形性に加えて弾性係数の位置依存的な同定を可能にしている。しかしながら、主として、この方法は、超音波により実行されてもよい。

請求項４は、請求項３の方法のバージョンを開示しており、基本周波数及び基本周波数の２倍の周波数におけるボクセルの偏位の同時測定を可能としている。しかしながら、代わりに、基本周波数及び基本周波数の２倍の周波数を有する成分から磁気勾配を構成することも可能である。

請求項５に開示のさらなるバージョンは、磁場の不均一性を位相における変形効果を加えることから阻止する。この場合、勾配磁場の２つの個々の正弦的振動は、同一の位相位置を有していなければならない。

請求項６に開示のバージョンは、（位置依存的な）非線形性及び組織の弾性が、組織の向上された分類を可能にするという事実を考慮している。しなしながら、請求項３に開示の方法は、この目的のための追加の測定を必要としない。

請求項７は、本発明に従った方法を実行するためのＭＲ装置を規定しており、請求項８は、この種のＭＲ装置を制御するためのコンピュータープログラムを開示している。

以下に、図面に対する参照文にて本発明を詳細に述べる。

サンプルにおける圧力Ｐとこの圧力に起因したサンプルの変形とは以下の関係を有している。

ここで、Ｅは、組織の弾性である。変形に関してεは：

である。ここで、Ｌ_０は、非圧縮サンプルの範囲であり、Ｌは、圧縮サンプルの範囲である。既に述べたが、生体組織サンプルの場合、変形εは、圧力Ｐに比例しない。なぜなら、弾性は、変形が増加するにつれ増加するためである。この増加は、簡易に以下の式で近似されてもよい。

ここで、Ｅ_０は、初期弾性（つまり、圧縮しない際の弾性）であり、指数αは非線形性の測定値である。αが０の場合、Ｐとεとの間に線形的な関係が存在する。αが大きくなると、この関係にもある通り、より非線形的になる。

図１は、縦軸に弾性Ｅ_０を、横軸に非線形性を特徴づける指数αをプロットしたグラフを示す。図１の各サンプルに関して、このサンプルに関する弾性Ｅ_０及び指数αを特徴づけるポイントを関連付けられてもよい。このグラフにおいて、複数の組織サンプルが、個々のシンボルにてマークされている。シンボル「＋」は、乳癌由来の組織サンプルを示している。シンボル「○」は線維腺腫の形態の良性病巣を示しており、シンボル「□」は、脂肪組織を示しており、シンボル「△」は良性から悪性状態へと変化する可能性のある乳管の成長組織サンプルを示している。

乳癌の組織サンプルは、より高い弾性を有し（つまり、より硬いということである）、同時に、健康な組織と比較して圧力Ｐと変形εとの間により顕著な非線形的な関係を有する。従って、組織サンプルは、これらパラメーターに基づいて分類されてもよい。このグラフの直線にちょうど乗っている組織サンプルは、悪性病巣を示しており、一方、このグラフの左側に配置された組織サンプルは、良性病巣であることを示す。

図１のグラフは、組織サンプルが増加する圧力Ｐに曝露され、且つ変形εが測定される間の組織サンプルのｉｎ ｖｉｔｒｏでの検討に基づいている。この方式における患者の組織検討は、困難ではない。本発明は、弾性撮影法を利用しながら、同一のパラメーターが患者に関して非侵襲的に同定されるという事実の認識に基づいている。対象が正弦的な機械的振動により励起される際、弾性が一定であり変形に独立している場合、機械的組織におけるボクセルの偏位は、機械的振動につれ経時的に同様の変化を示す。しかしながら、ボクセルが圧力Ｐと変形εとの間で非線形的な関係を示す組織内に存在する場合、ボクセルの偏位は、非線形性の測定で発生する振動及び変形の経時的な変化をもはや正確にフォローし得ない。これら変形及び非線形性の測定を、ＭＲ弾性撮影法を含む例に基づいて以下に詳細に述べる。

図２は、この方法を実行するためのＭＲ装置の概略的な断面図であって、この制御に必要なハードウェア及びソフトウェア成分も示している。このＭＲ装置は、図２の平面に対して垂直及び平行に延びる主磁場Ｂ_０を生成するための主要フィールドマグネット３０を有している。この円筒形の主要フィールドマグネットの内部には、３つの共通の垂直方向に勾配を有する磁場を発生する機能を有する磁場コイルのセット３１が配置されている。これら勾配コイルの内部には、ＲＦ磁場を発生するためのＲＦコイル３２が配置されている。

勾配コイル３１用の電流は、ジェネレーター構造３３により供給され、ＲＦコイル３２用の電流は、ＲＦジェネレーター３４により供給される。ジェネレーター３３及び３４により発生される電流の経時的な変化は、コントロールユニット３５により同定され制御される。乳房にて検査中に発生したＭＲ信号は、レシーバー３６にてさらに処理されるように、ＭＲコイル（図示せず）により受信される。レシーバー３６の出力信号は、受信したＭＲ信号からＭＲ画像を再構築する画像処理ユニット３７に適用される；このＭＲ画像は、モニター３８上に表示される。

コイル３１及び３２により規定された検査領域において、患者テーブル２に配置している勾配コイル３１上に配置された女性患者４７が配置される。勾配コイル３１の前面では、クッション４１と同一の高さを有するマンモグラフィーアクセサリー１が配置されている。患者の乳房は、若干圧縮されつつ機械的振動が暴露されているマンモグラフィーアクセサリーの開口部にぶら下がっており、この振動は、経時的に正弦的に変化しつつ、（脚から頭部へと）水平方向に伝播する。マンモグラフィーアクセサリー１により発生される機械的振動の経時的な変化は、ジェネレーター３３及び３４により発生された電流が経時的に変化しつつ、コントロールユニット３５により同期される。上記のＭＲ装置は、独国特許出願第１０１５６１７８．４号（ＰＨＤＥ０１０３０２）に詳細に述べられており、重複を避けるために参考文とする。

図３は、図２に示したＭＲ装置によりＭＲ弾性撮影法検査のためのＭＲデータの取得を経時的な実行を示している。初期化（１００）が完了し患者４０が図１と同様に配置されると、患者の乳房は、固定され、弾性的に変形可能であるように若干圧縮される。続いて、経時的に正弦的に変化する機械的振動は、マンモグラフィーアクセサリー１において発生され、この振動は、患者の乳房に作用し、ＭＲ取得全体に関して持続する（ブロック１０１）。

この振動が静止した後、ブロック１０２では、図４に示した経時的に変化する第１に生成した第１ＭＲシーケンスが存在する。このシーケンスは、９０°のＲＦパルスＲＦ_１と続いて再集束した（ｒｅｆｏｃｕｓｉｎｇ）１８０°のＲＦパルスＲＦ_２（図４の４番目のライン）とを有している。二つのＲＦパルスＲＦ_１とＲＦ_２との間に、単一で一時的で正弦的に変化する勾配磁場の振動Ｇ_１が発生する（図４の三番目のライン）；この振動の持続時間は、機械的振動の持続時間のちょうど半分（例えば、５０Ｈｚの機械的振動の周波数の場合、１０ｍｓ）に等しい。第２ＲＦパルスの後、勾配磁場Ｇの第２の振動Ｇ_２が発生される：この第２振動は、振動Ｇ_１とちょうど同じ経時的変化を有しており、その開始は、第１振動の終端に相対して、機械的振動の１／４の時間（５ｍｓ）だけシフトされている。

２つのＲＦパルスＲＦ_１及びＲＦ_２の間、主勾配磁場は、ｘ方向に適用され、この勾配磁場は、ＲＦパルスが、患者の長手方向に対して垂直に延びるスライスにおいて排他的に活性化することを可能としている；勾配及び読取勾配（ｒｅａｄ ｇｒａｄｉｅｎｔ）をコード化する位相と同様に、この勾配磁場は、図４に示していない。

図４の第１のラインは、機械的振動の経時的な変化の一部を示しており、第２のラインは、２倍の周波数におけるこれら振動のより高い調波を示している。かかるより高い調波は、機械的振動により励起された波が乳房に伝播し且つかかる領域に到達する際、変形と弾性との間で非線形的な関係が存在する領域において連続的に発生してもよい。さらに、これら２つの第１のラインの破線は、特定の周波数において屈折されたスピンがパルスＲＦ_２により再集束されたことを示している。

勾配振動Ｇ_１及びＧ_２は、機械的振動により屈折された核スピンの位相を変化させる。この変化の程度は、特定の勾配の方向における偏位の程度に依存するが、この変化の兆候（ｓｉｇｎ）は、勾配振動と核スピンの偏位との間の位相位置に依存する。

図４から明らかなのは、勾配振動Ｇ_１の兆候が変化する一方で、機械的な５０Ｈｚの振動の兆候（第１のライン）は同様に保持されているということである。従って、位相の変化は起こらない。しかしながら、第２の勾配振動Ｇ_２は、位相の変化を引き起こす。なぜなら、５０Ｈｚの振動は、勾配振動Ｇ_２と同様の程度において交差するゼロを有しているためである。この２つの勾配振動Ｇ_１及びＧ_２は、基本周波数の調波（１００Ｈｚ）において屈折されたスピンに関する位相の同一の大きな変化を生じ、この変化は、示した変化に起因して加算される。

図４に示したシーケンスは、基本周波数における偏位と共に基本周波数の２倍の周波数における偏位に対して感受的である。この場合でないとすると、偏位の基本周波数とこの調波は、測定時間が遅延するように分離して測定されなければならない。

ステップ１０２において得たＭＲ信号の受信の後、このＭＲシーケンスは、（図２において水平方向に延びている）ｚ軸方向にシフトされる、異なるスライスに関して、ステップ１０３において繰り返される（ブロック１０３）。まれに、ＲＦパルスＲＦ_１及びＲＦ_２の搬送周波数は、勾配磁場の経時的な変化が同一に保持されつつ他のスライスが励起される様式にて変化される。ステップ１０２及び１０３を有するループは、検査される三次元領域におけるスライスと同一の回数、例えば、２０回、繰り返され、各時間において、異なるスライスが励起される。

続いて、勾配をコード化する位相は、ステップ１０４において変化され、全てのスライスは、そこで生成されるスピンエコー信号を取得するように再び励起される。個々のスライスは、機械的振動の関係においてＲＦパルス及び勾配磁場の経時的に同一の位置で励起される。勾配をコード化する位相は、要求される位相コード化ステップと同一の回数、例えば１２８又は２５６回、繰り返される。従って、ＭＲ画像は、そこから再構築可能なように全てのスライスから取得される；この画像は、スライスにより形成された三次元領域における核磁気分布の位相を示している。

三次元領域の第１のＭＲ位相画像に関するＭＲデータの取得の後に、ステップ１０５では、経時的なシフトがＭＲシーケンスと機械的振動との間で行われ、つまり、ＭＲ実験の開始が時間の画分ΔＴ、例えばΔＴ＝Ｔ／８（ステップ１０５）だけシフトされるようにシフトが起こる。続いて、ステップ１０２、１０３及び１０４を有するループは、検査される三次元領域のさらなるＭＲ位相画像が取得されるように再度完了される。この２つのＭＲ位相画像は、そこで励起される組織又は核スピンがｘ軸方向に屈折される領域に関してのみ異なっている。その後、機械的振動とシーケンスとの間の経時的なシフトが変更される（例えば、ΔＴ＝２Ｔ／８、３Ｔ／８・・・７Ｔ／８）さらなる位相画像が形成される。ユーザーに入力され或いはユーザーによるエントリーにより入力されたデータは、モニター８上に表示されてもよく、或いは、異なるディスプレイ上に表示されてもよい。

最後に、複数のＭＲ位相画像のＭＲデータ（この例では８個）が取得され、このＭＲデータは、勾配磁場Ｇ_１、Ｇ_２の勾配がｘ軸方向に延びることを仮定すると、ｘ軸方向におけるスピンの偏位に依存している。

ｘ軸方向以外の方向にも偏位の測定を可能とするために、正弦的勾配Ｇ_１及びＧ_２は、ステップ１０２、１０３、１０４及び１０５を有するループが完了した後、ｙ軸方向におけるスピンの偏位により影響される三次元領域のＭＲ射創画像のセットを取得するように、ｘ軸とは異なる方向、例えば、ｙ軸方向、に形成される（ブロック１０６）。最終的に、この工程は、最終的に、ｘ軸方向、ｙ軸方向又はｚ軸方向におけるスピンの偏位にそれぞれ依存する８個のＭＲ位相画像の３つのセットが利用可能となるように、ｘ軸方向に関して繰り返される（つまり、正弦勾配Ｇ１及びＧ２はｚ軸方向に適用される）。この取得方法は、その後終了される（ブロック１０７）。

取得したＭＲデータの処理は、図６に対して参照しながら以下に詳細に述べる。既に示しているように、ブロック２００における初期化の後、ＭＲ位相画像は、ステップ１において、取得したＭＲデータから形成される。ステップ２０２において、偏位ｕ_ｘ、ｕ_ｙ及びｕ_ｚの経時的な変化は、３つの空間的方向、つまり、三次元検査ゾーンの個々のボクセルに関するそれぞれで同定される。

示すことを目的として、図５に概略的に表示した。上部ラインは、それぞれ、ｘ軸方向、ｙ軸方向及びｚ軸方向における偏位に依存するＭＲ位相画像の３つのセットを概略的に示している。セットに関連付けられた各ＭＲ画像において、その位相は、同一のボクセルにおいて考慮される。８つの画像は、特定のボクセルにおけるｘ軸方向、ｙ軸方向又はｚ軸方向における位相又は偏位の経時的な変化に関する８個のベースポイントを生成する。このボクセルが非線形的な弾性特性を示す領域に配置されている場合、このボクセルにおける偏位は、経時的に純粋に正弦的ではない。なぜなら、このボクセルにおける屈折は、基本的な振動の周波数を有する成分に加えて、より高い調波、特に基本周波数の２倍の成分を有しているからである。

これら成分の振幅（ａｍｐｌｉｔｕｄｅ）を同定するために、ステップ２０３において、上記の８つのベースポイントにより規定されたボクセルの位相又は偏位に対してフーリエ変換が実行される。このフーリエ変換は、基本周波数におけるｘ軸方向、ｙ軸方向及びｚ軸方向における偏位の振幅を示す３つの量ａ_１ｘ、ａ_１ｙ及びａ_１ｚを生成し、且つ、基本周波数の２倍におけるｘ軸方向、ｙ軸方向及びｚ軸方向における偏位の振幅を示す３つの量ａ_２ｘ、ａ_２ｙ及びａ_２ｚも生成する。基本周波数を有する成分の振幅ａ_１は、以下の方程式に一致して算出されてもよい：

基本周波数の２倍を有する成分の振幅に関して、

に従う。

上記の２つの成分は、以下の方程式に一致して、他の一つと関連している：

非線形的な測定である値が取得される。振幅ａ_２自体、既に非線形の測定であるが、ａ_２は、非線形性のみならず物体における励起により発生する波の減衰にも依存する。パラメーターａは、高い程度で非依存的である。これは、方程式３に一致した非線形の値αに対応している。この対応は、より高い項目（ｔｅｒｍ）を考慮しつつ、変形張筋（ｄｉｓｔｏｒｔｉｏｎ ｔｅｎｓｏｒ）とともに自由エネルギーの正確な適用を介して与えられる。αとａとの関係は、厳格に単調に増加する。従って、図１と同様に、等価な診断的見解は、Ｅ０及びａにより与えられてもよい。

個々の値が全てのボクセルに関して同定された後、ステップ２０４において、検査ゾーンにおける弾性の非線形性を示す画像が形成される。

非線形性自体が、一般に、図１に対する参照文にて既述したように、病巣の良性又は悪性特性として表現されることを可能としないので、ステップ２０５において、検査ゾーンにおける弾性係数の空間的分布を示す画像がさらに生成される。この目的に必要な測定値は、図３と一致した取得の後に既に利用可能である。これら値から弾性係数を算出することは、上述した米国特許出願番号第０９／７４３６５９号（ＰＨＤ９９．１５２）明細書にて詳細に述べられている。

その後、ユーザーは、例えばカーソルによりマークされた領域など、画像における個々の領域を選択してもよく、且つ、Ｅ値を示してもよく、且つ、これらの領域に関して、例えば、図１に概略的に一致した領域に関して示してもよい（ステップ２０６）。この方法はその後終了する。

既述したように、超音波弾性撮影法手段もまた、物体における機械的振動により励起される波の伝播を追跡するのに適しており、かかる追跡は、ＭＲ弾性撮影法よりも有意に速い。組織の非線形的な特性をもたらす経時的な非線形的な変形は、上述と同様の様式にてその後同定されてもよい。乳房に加えて、他の領域を検討してもよく、例えば、前立腺などが挙げられる。

異なるタイプの組織の分類に関するグラフを示す。 本発明に従った方法を実行するためのＭＲ装置を示す。 ＭＲ位相画像の取得を示すフローチャートを示す。 ＭＲ取得を意図したシーケンスの一つを示す。 評価工程の一部の概略的表示である。 評価工程のフローチャートを示す。

Claims (8)

1. 物体の検査方法であって：
   ａ）機械的振動により前記物体を励起するステップ；
   ｂ）前記物体において前記振動により生じた波に指向される前記物体のボクセルの静止状態から、偏位の経時的な変化を測定するステップ；
   ｃ）前記の偏位の経時的な変化から、非線形的変形を同定するステップ；及び
   ｄ）前記非線形的変形を評価するステップ；
   を有することを特徴とする方法。
2. 前記の評価するステップは、前記非線形的変形の程度を示す画像を形成するステップを有することを特徴とする請求項１に記載の方法。
3. ａ）一時的で周期的な、好ましくは正弦的に変化する機械的振動により前記物体を励起するステップ；
   ｂ）前記機械的振動に同期する勾配磁場と組み合わせて前記物体に核磁気を励起し、且つ、ＭＲ位相画像を形成するように前記物体に生じるＭＲ信号を受信するステップであって、前記勾配磁場の経時的な変化が、前記ＭＲ画像が基本周波数における偏位、及び前記基本周波数の少なくとも一つのより高い調波における偏位、により同定されるように選択される、ステップ；
   ｃ）さらなる位相画像を形成するように、前記勾配磁場の勾配の方向を変化し、及び／又は前記機械的振動と前記勾配磁場との位相差を変化しつつ、複数回前記ステップｂ）を繰り返すステップ；
   ｄ）前記ＭＲ位相画像に基づいて前記基本周波数におけるスピンの偏位の振幅を同定するステップ；
   ｅ）前記ＭＲ位相画像に基づいて前記基本周波数のより高い調波におけるスピンの偏位の振幅を同定するステップ；並びに
   ｆ）前記振幅の比率に依存する画像を形成するステップ；
   を有することを特徴とする請求項２に記載の方法。
4. 前記勾配磁場は、前記機械的振動の持続時間の半分の、２つの異なる正弦的振動を有しており、
   前記振動間の時間における距離は、前記機械的振動の前記持続時間の１／４であることを特徴とする請求項３に記載の方法。
5. 再集束ＲＦパルスは、前記の２つの正弦的振動間の前記検査ゾーンに作用することを特徴とする請求項４に記載の方法。
6. ａ）前記検査ゾーンにおける弾性を追加的に算出するステップ；並びに
   ｂ）同一のボクセルに関して同定された変形及び弾性の値を評価するステップ；
   を有することを特徴とする請求項３に記載の方法。
7. ＭＲ装置、機械的振動ジェネレーター、評価ユニット、勾配磁場の経時的な変化を同定するジェネレーター並びに前記ＭＲ装置、前記ジェネレーター及び前記評価ユニットを制御するコントロールユニットにより請求項１に記載の方法を実行する装置であって：
   ａ）一時的で周期的な、好ましくは正弦的に変化する機械的振動により前記物体を励起するステップ；
   ｂ）前記機械的振動に同期する勾配磁場と組み合わせて前記物体に核磁気を励起し、且つ、ＭＲ位相画像を形成するように前記物体に生じるＭＲ信号を受信するステップであって、前記勾配磁場の経時的な変化が、前記ＭＲ画像が基本周波数における偏位、及び前記基本周波数の少なくとも一つのより高い調波における偏位、により同定されるように選択される、ステップ；
   ｃ）さらなる位相画像を形成するように、前記勾配磁場の勾配の方向を変化し、及び／又は前記機械的振動と前記勾配磁場との位相差を変化しつつ、複数回前記ステップｂ）を繰り返すステップ；
   ｄ）前記ＭＲ位相画像に基づいて前記基本周波数におけるスピンの偏位の振幅を同定するステップ；
   ｅ）前記ＭＲ位相画像に基づいて前記基本周波数のより高い調波におけるスピンの偏位の振幅を同定するステップ；並びに
   ｆ）前記振幅の比率に依存する画像を形成するステップ；
   を実行するようにプログラムされた装置。
8. 請求項１に記載された：
   ａ）一時的で周期的な、好ましくは正弦的に変化する機械的振動により前記物体を励起するステップ；
   ｂ）前記機械的振動に同期する勾配磁場と組み合わせて前記物体に核磁気を励起し、且つ、ＭＲ位相画像を形成するように前記物体に生じるＭＲ信号を受信するステップであって、前記勾配磁場の経時的な変化が、前記ＭＲ画像が基本周波数における偏位、及び前記基本周波数の少なくとも一つのより高い調波における偏位、により同定されるように選択される、ステップ；
   ｃ）さらなる位相画像を形成するように、前記勾配磁場の勾配の方向を変化し、及び／又は前記機械的振動と前記勾配磁場との位相差を変化しつつ、複数回前記ステップｂ）を繰り返すステップ；
   ｄ）前記ＭＲ位相画像に基づいて前記基本周波数におけるスピンの偏位の振幅を同定するステップ；
   ｅ）前記ＭＲ位相画像に基づいて前記基本周波数のより高い調波におけるスピンの屈折の振幅を同定するステップ；並びに
   ｆ）前記振幅の比率に依存する画像を形成するステップ；
   を有する方法を実行するように、ＭＲ装置、振動ジェネレーター及び評価ユニットに作用するコントロールユニット用コンピュータープログラム。

Images