(57)【特許請求の範囲】
【請求項１】
対象内における拡散プロセスを記述するための磁気共鳴による拡散テンソルの係数を決定するために、対象のボリューム要素から、種々の拡散エンコーディングのための制御データの考慮のもとに時間順に空間分解され種々に拡散エンコーディングされた拡散データが発生される磁気共鳴による拡散テンソルの係数の決定方法において、
各ボリューム要素について、最初の拡散エンコーディングにより測定された拡散データから拡散中間データが発生され、
各ボリューム要素について、先行して発生された拡散データの拡散中間データと、他の拡散エンコーディングにより測定された現在の拡散データとから、更新拡散中間データが決定され、
種々の拡散エンコーディングによる複数回の更新後に、最後に決定された更新拡散中間データから、拡散テンソルの成分が決定される
ことを特徴とする磁気共鳴による拡散テンソルの係数の決定方法。
【請求項２】
拡散中間データおよび更新拡散中間データは、ムーア・ペンローズの擬似逆の使用のもとに正規方程式により決定されることを特徴とする請求項１記載の方法。
【請求項３】
更新拡散中間データは、先行して決定された拡散中間データに、現在の測定から求められた増分の拡散中間データを加算することによって決定されることを特徴とする請求項１又は２記載の方法。
【請求項４】
増分の拡散中間データを決定するために、拡散データの自然対数が使用されることを特徴とする請求項３記載の方法。
【請求項５】
増分の拡散中間データを決定するために、個々の互いに独立な方向における拡散傾斜磁場を表すベクトルが、対数化された拡散データにより重み付けされることを特徴とする請求項３又は４記載の方法。
【請求項６】
拡散中間データおよび更新拡散中間データが特異値分解法により決定されることを特徴とする請求項１又は３乃至５の１つに記載の方法。
【請求項７】
拡散テンソルの成分を決定するために、付加的に画像化傾斜磁場の拡散エンコーディングされる作用が考慮されることを特徴とする請求項１乃至６の１つに記載の方法。
【請求項８】
制御データは拡散データの発生にともなって更新されることを特徴とする請求項１又は５乃至７の１つに記載の方法。
【請求項９】
制御データは測定経過の決定後に拡散データの発生前に決定されることを特徴とする請求項１乃至８の１つに記載の方法。
【請求項１０】
拡散データは、制御データの状態数に基づいて決定される時間順で、先ず拡散テンソルの成分を強く決定する拡散データが発生され、その後に拡散テンソルの成分を弱く決定する拡散データが発生されるように発生されることを特徴とする請求項１乃至９の１つに記載の方法。
【請求項１１】
拡散データは、少なくとも７つの種々の拡散エンコーディングにより発生されることを特徴とする請求項１乃至１０の１つに記載の方法。
【請求項１２】
拡散エンコーディングの１つは零であることを特徴とする請求項１乃至１１の１つに記載の方法。
【請求項１３】
拡散テンソルの成分が予め与えられた精度で決定され得る際に、信号が出力されることを特徴とする請求項１乃至１２の１つに記載の方法。
【請求項１４】
拡散テンソルの成分のための予め定められた精度基準が到達された際に、測定が終了されることを特徴とする請求項１乃至５の１つ又は７乃至１３の１つに記載の方法。
【請求項１５】
拡散テンソルの成分の良さが決定され、品質尺度として出力されることを特徴とする請求項１乃至１４の１つに記載の方法。
【請求項１６】
制御データから状態数が拡散テンソルの成分を決定するための尺度として求められることを特徴とする請求項１乃至１５の１つに記載の方法。
【請求項１７】
拡散テンソルの成分が画像表示されることを特徴とする請求項１乃至１６の１つに記載の方法。
【請求項１８】
請求項１乃至１６の１つに記載の方法を実施するための装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、対象内における拡散プロセスを記述するための磁気共鳴による拡散テンソルの係数を決定するために、対象のボリューム要素から種々の拡散エンコーディングのための制御データの考慮のもとに時間順に空間分解され種々に拡散エンコーディングされた拡散データが発生される磁気共鳴による拡散テンソルの係数方法に関する。
【０００２】
本発明は、同様にこの方法を実施するための装置に関する。
【背景技術】
【０００３】
上述の如き方法および装置は公知である（例えば、特許文献１参照）。これに示された方法は１９６５年にＳｔｅｊｓｋａｌおよびＴａｎｎｅｒによって提案された磁気共鳴測定シーケンスに基づいている。この特許文献１に記載された方法は、ＳｔｅｊｓｋａｌおよびＴａｎｎｅｒによって提案された測定シーケンスを、検査領域から各ボクセルについて拡散テンソルの係数または要素が求められて画像表示されるように拡張している。具体的には、最も強い拡散プロセスの方向を再現する主軸を有する拡散楕円面による表示が行なわれる。個々の方向における広がりは、対応方向における拡散プロセスの数値を表している。医学用途においては、拡散の磁気共鳴画像化により、水分子から体内組織への相対移動度が測定される。組織への拡散は例えば繊維方向の如き組織構造に依存し得ることから、拡散テンソルから医学上重要な情報を導き出すことができる。
【０００４】
拡散データは、磁気共鳴測定の場合、拡散エンコーディングに使用された拡散傾斜磁場の大きさおよび方向によって決定される。この場合に、高周波１８０°リフォーカスパルスに対して対称的に配置されている強い傾斜磁場パルスが使用される。１８０°リフォーカスパルスの前の第１傾斜磁場パルスは全てのスピンに対して位相シフトを生じさせ、第２傾斜磁場パルスはこの位相シフトを反転させる。それにより、固定の分子（医用画像化においてはプロトン）の場合には、位相シフトが再び打ち消される。しかしながら、ブラウン運動により第２傾斜磁場パルスの作用時に第１傾斜磁場パルスの作用時とは異なる位置にある分子の場合、位相シフトが完全に相殺されない。信号減弱をもたらす残留位相シフトがとどまる。このような傾斜磁場パルスの大きさおよび方向を介して拡散エンコーディングが制御される。
【０００５】
磁気共鳴技術およびこれに続く画像表示による医用拡散テンソル測定（ＤＴＩ＝Ｄｉｆｆｕｓｉｏｎ  Ｔｅｎｓｏｒ  Ｉｍａｇｉｎｇ）においては、大量の測定データが生じ、これらの測定データからボクセルごとに、得ようとする６つのテンソルパラメータまたはテンソル係数、すなわち３×３の対称テンソル行列の独立要素または成分が算出される。典型的には、引続いてテンソルパラメータは診断に重要なパラメータマップを算出するために使用される。例えば、拡散テンソルの等方性成分または拡散テンソルの異方性成分が相応のパラメータマップに示される（平均見かけ拡散係数マップまたはＡＤＣａｖマップもしくは部分異方性マップまたはＦＡマップ）。多数の測定データのために、ボクセルごとの拡散テンソルの算出は平均化による未知パラメータの決定を必要とする。多変量の線形回帰から知られた方法、例えば、擬似逆を形成するかまたは特異値分析を行なう方法が使用される。これらの方法は、もちろんメモリ占有率が非常に高く、計算費用がかかる。なぜならば、公知の方法の場合、全データが評価のために予め保持され、求めようとするパラメータが測定終了後に平均化方法によって決定されるからである。代替として、拡散エンコーディングの特定のセットに対して強度値の換算規則が適用される。しかしながら、この種の方法は非常に柔軟性がない。
【特許文献１】米国特許第５５３９３１０号明細書
【発明の開示】
【発明が解決しようとする課題】
【０００６】
本発明の課題は、柔軟性のある評価が可能であり、所要メモリ占有容量が低減されている方法および装置を提供することにある。
【課題を解決するための手段】
【０００７】
この課題は、本発明によれば、請求項１に記載された方法および請求項１８に記載された装置によって解決される。この場合に、各ボリューム要素について、最初の拡散エンコーディングにより測定された拡散データから拡散中間データが発生され、各ボリューム要素について、先行して発生された拡散データの拡散中間データと他の拡散エンコーディングにより測定された現在の拡散データとから更新された拡散中間データすなわち更新拡散中間データが決定され、種々の拡散エンコーディングによる複数回の更新後に、最後に決定された更新拡散中間データから、拡散テンソルの成分が決定される。
【０００８】
この増分方法により、到来する測定値を即座に処理して、拡散中間データの更新に使用することができる。最後に更新された中間データから、計算の割合が非常に多くはない少ない処理ステップにて、テンソル成分を決定することができる。それゆえ、従来技術のように、評価が開始できる前に、先ず全ての測定値を求めてメモリに格納することは必要でない。所要メモリが本発明による方法の場合には著しく減少する。計算の一部は既にデータ取得と並行して行なうことができ、測定後に結果計算のために必要な時間も同様に著しく減少する。測定終了後になおも残っている計算はもはや直接に測定データを用いない。計算は測定データへのアクセスに要する部分と、中間結果に積み重ね得る部分とに分けられている。
【０００９】
第１の特に有利な実施態様においては、拡散中間データおよび更新拡散中間データが、ムーア・ペンローズの擬似逆の使用のもとに正規方程式により決定される。この方法は、各部分測定にとって、つまり定められた拡散エンコーディングを用いる各測定にとって、この部分測定のための情報およびそれからの情報だけが必要である。したがって、原則的には予測できない測定中断後にも拡散成分を決定することができる。その際に、結果の精度はもちろん測定データの品質に依存する。
【００１０】
他の有利な実施態様においては、拡散中間データおよび更新拡散中間データが特異値分解法により決定される。この場合、確かに、"（測定個数）×（６つの決定すべき拡散成分）"のように全ての拡散エンコーディングに対する情報が必要であるが、それにもかかわらず、この場合にも測定データの増分処理が相応に低減された所要メモリ量および計算量で可能である。
【００１１】
上記の両実施態様の場合、高い最適化された標準ライブラリに利用できるかまたは現代のプロセッサのハードウェア拡張も利用する線形代数学からの行列計算のための標準方法が使われる。
【００１２】
他の有利な実施態様においては、拡散テンソルの成分を決定するために、付加的に画像化傾斜磁場の拡散エンコーディングされる作用が考慮される。それにともなって、そのように決定された拡散テンソル成分の品質改善がもたらされる。なぜならば、拡散エンコーディングに影響を及ぼす全ての傾斜磁場が一緒に考慮されるからである。
【００１３】
他の有利な実施態様において、拡散データは、制御データの状態数に基づいて決定される時間順で、先ず拡散テンソルの成分を強く決定する拡散データが発生され、その後に拡散テンソルの成分を弱く決定する拡散データが発生されるように発生される。短い測定時間後に既に拡散テンソルの成分のための一次近似が算出される。さらにこの種の測定順序は予定より早い測定中断に対して寛大である。
なお、別の有利な実施態様において、拡散データは、少なくとも７つの種々の拡散エンコーディングにより発生される。さらに別の有利な実施態様において、拡散エンコーディングの１つは零である。
【００１４】
他の有利な実施態様においては、拡散テンソルの成分が予め与えられた精度で決定され得る際に信号が出力される。別の有利な実施態様において、拡散テンソルの成分のための予め定められた精度基準が到達された際に、測定が終了される。他の有利な実施態様においては、拡散テンソルの成分の良さが決定され、品質尺度として出力される。
【００１５】
他の有利な実施態様においては、制御データから状態数が拡散テンソルの成分を決定するための尺度として求められる。
なお、別の有利な実施態様において、更新拡散中間データは、先行して決定された拡散中間データに、現在の測定から求められた増分の拡散中間データを加算することによって決定される。
増分の拡散中間データを決定するために、拡散データの自然対数が使用される。増分の拡散中間データを決定するために、個々の互いに独立な方向における拡散傾斜磁場を表すベクトルが、対数化された拡散データにより重み付けされる。
他の有利な実施態様において、制御データは拡散データの発生にともなって更新される。別の有利な実施態様において、制御データは測定経過の決定後に拡散データの発生前に決定される。
別の有利な実施態様において、拡散テンソルの成分が画像表示される。
【発明を実施するための最良の形態】
【００１６】
以下において、本発明の実施例を２つの図に基づいて説明する。
図１は拡散テンソルを決定するための主要な構成要素を有する診断用磁気共鳴装置の概略図を示し、
図２は拡散テンソルの成分を決定する方法の基本構成をブロック図で示す。
【００１７】
生体組織の拡散特性はそのミクロ構造に対する重要な情報を与える。この種の拡散過程は拡散テンソルによって記述され、その成分は磁気共鳴技術によって求められる。拡散エンコーディングされたまたは拡散強調された磁気共鳴信号は、今日において公知の診断用磁気共鳴装置により発生させることができる。
【００１８】
拡散強調された磁気共鳴信号は次式にしたがって記述される。
【数１】
0005235283.tif 000002
ただし、
Ｓ0      拡散強調された傾斜磁場パルスなしの初期信号強度
ｂ        拡散強調係数
Ｄ        拡散係数（ＡＤＣ、「見かけ上の拡散係数」とも呼ばれる。）
【００１９】
式（１）による拡散係数Ｄにより、等方性の拡散特性を有する組織のみを記述することができる。異方性の拡散の場合も把握することができるようにするために、式（１）は次のように拡張される。
【数２】
0005235283.tif 000003
Ｄは対称の３×３の拡散テンソルを記述する。
【数３】
0005235283.tif 000004
そして、
【数４】
0005235283.tif 000005
は、正規化され拡散エンコーディングされる傾斜磁場ベクトルを記述する。Ｄxx，Ｄxy，Ｄxz，Ｄyy等は拡散テンソルの個々の成分または要素である。時々、これはテンソルパラメータとも呼ばれる。
【００２０】
使用される傾斜磁場の作用は、拡散行列Ｂの導入によって記述することができる。これは次のとおり展開される。
【数５】
0005235283.tif 000006
ここでは簡略表示の理由から、座標方向ｘ，ｙ，ｚのインデックスの代わりに、式においては類似させた数値インデックス１，２，３を使用した。Ｂ行列により、全ての拡散画像化傾斜磁場を測定シーケンスにおけるそれらの交差項と同様に考慮することができる。
【００２１】
拡散テンソルを決定するために、種々の拡散エンコーディング方向による少なくとも６つの測定が必要である。拡散傾斜磁場（ベースライン画像）の作用なしの追加測定により、拡散テンソルの要素および各ボクセルのためのＴ２強調信号Ｓ0が測定範囲または検査範囲から決定される。
【００２２】
７よりも多い測定が行なわれる場合、これは、良好な信号雑音比の理由から好ましく、係数はその際に生じる誤差が最小になるように決定される。この場合に次の関係が生じる。
【００２３】
測定または観察：
測定回数：    ｎ                              （６）
測定ベクトル：
【数６】
0005235283.tif 000007
測定方程式：
【数７】
0005235283.tif 000008
【００２４】
未知数：
未知数の個数：    ｕ＝７                        （９）
未知数のベクトル：
【数８】
0005235283.tif 000009
【００２５】
方程式系の冗長性：    ｒ＝ｎ−ｕ                  （１１）
【００２６】
解モデル：
モデル関数：
【数９】
0005235283.tif 000010
改善のための目標関数ｖ：
【数１０】
0005235283.tif 000011
解：
【数１１】
0005235283.tif 000012
計画行列：
【数１２】
0005235283.tif 000013
これは次の偏微分をともなう。
【数１３】
0005235283.tif 000014
【００２７】
拡散テンソルは通常は患者座標系において指定される。より広範囲な解析のために、テンソルの固有ベクトルの方向に対応する主軸を有する座標系への座標変換を行なうと有利である。このために、拡散行列Ｄから固有値−固有ベクトル分解を行なうことが必要である。
【数１４】
0005235283.tif 000015
ただし、
ｅi:方向１，２，３における単位固有ベクトル、
ｅ1：主拡散方向、
λi：ｅiの方向における拡散率（λ1≧λ2≧λ3）。
【００２８】
逆の正規方程式行列Ｎ-1は、未知のテンソル要素および評価されたＳ0値の精度に関する情報を含んでいる。
【数１５】
0005235283.tif 000016
であり、
【数１６】
0005235283.tif 000017
である。
【００２９】
行列Σxxは拡散エンコーディングのための傾斜磁場パルスによって完全に定められている。これは、不確実の形状または不確実の空間的姿に関する情報を供給する。
【００３０】
スカラーσ02は、ピクセルからピクセルへと変化し、測定データとテンソルモデルとの整合性の良さに関する情報を提供する。
【００３１】
正規方程式行列の状態数は、方程式（１４）が当面のデータにより安定した解をもたらすかどうかについての尺度として使用可能である。状態数の算出は拡散エンコーディングのための制御データだけを必要とする。行列は、その状態数が無限である場合には、特異あり、したがって解くことができない。状態数が大きい場合には、行列は不良状態にある。それは、状態数の逆数が使用された計算法の浮動小数点精度に近づくことと同意である。それにともなって、状態数の高い値は、測定方法の不良な計画品質に関する情報を与える。状態数を決定するのに公知の方法を使用することができる。
【００３２】
行列Ｎおよびベクトル
【数１７】
0005235283.tif 000018
の増分決定の第１の実現が、以下において導き出される。方程式（１６）の偏微分の使用のもとに、まずベクトルａi（ｉ＝１…ｍ）を導入する。
【数１８】
0005235283.tif 000019
【００３３】
検査領域の各測定後に行列Ｎiが更新される。なぜならば、
【数１９】
0005235283.tif 000020
なる関係が存在するからである。これに加えて、各ボクセルについてベクトル
【数２０】
0005235283.tif 000021
が、
【数２１】
0005235283.tif 000022
なる関係にしたがって更新される。
【００３４】
測定された各画像またはボリューム後、すなわち各拡散エンコーディング後に必要な更新は、ボクセルごとに７つの乗算および加算を含むベクトルａiおよび行列Ｎiに関する幾つかの少ない簡単な操作からなるだけである。この理由は、測定された画像もしくはボリュームごとのベクトルａiは変化しないままとどまることにある。これらの計算操作は測定中（インライン）も行なうことができる。
【００３５】
全ての画像もしくはボリュームが個々の拡散エンコーディングで測定され、測定値が上述の方程式（２１），（２２）にしたがって継続処理された後には、方程式（１４）において、行列Ｎiは行列Ｎに一致し、ベクトル
【数２２】
0005235283.tif 000023
はベクトル
【数２３】
0005235283.tif 000024
に一致する。行列Ｎの反転後に未知のテンソル成分が各ボクセルについて方程式（１４）にしたがって決定される。
【００３６】
同様に、現在の拡散中間データの精度に対する値は次のようにして決定することができる。重みユニットの分散がピクセルごとに、先に式（１９）にて示されているように算出される。ベクトルｖの成分が必要でない場合には、そのための効率のよい増分解法を指定することができる。式（１２），（１９）に基づいて、
【数２４】
0005235283.tif 000025
を導き出すことができる。式（２３）の議論が次の４つの成分を提供する。
【００３７】
【数２５】
0005235283.tif 000026
この成分は正規方程式行列Ｎ＝ＡTＡおよび解ベクトルを含む。その行列も解ベクトルも、式（２０）〜（２２）にて示されている増分解法との関連で既に算出されている。
【数２６】
0005235283.tif 000027
この成分は、上述と同様に増分解法の一部として既に決定された解ベクトルｘを含む。
【数２７】
0005235283.tif 000028
この成分は上述と同様に既に算出されている。
【数２８】
0005235283.tif 000029
この成分だけが付加的に算出され記憶されなければならない。
【００３８】
既に述べたように、σ^02の解は増分計算することができる。なおも決定されなければならない個々の付加的なデータはｌｎＳi（ｘ，ｙ，ｚ）の２乗の和である。
【００３９】
他の実施形態においては、拡散テンソルの要素の増分計算のために、特異値分解（ＳＶＤ；Ｓｉｎｇｕｌａｒ  Ｖａｌｕｅ  Ｄｅｃｏｍｐｏｓｉｔｉｏｎ）が使用される。特異値分解も過剰決定された連立方程式の最小二乗誤差解法を見いだす。原則的には、列数Ｎに等しいか又は大きい行数Ｍを有する各Ｍ・Ｎ行列Ａは、Ｍ・Ｎ列直交行列Ｕと正または零要素（特異値）を有するＭ・Ｎ対角線行列Ｗおよび転置されたＭ・Ｎ直交行列Ｖとの積として指定される。
【数２９】
0005235283.tif 000030
【００４０】
方程式（２５）により、方程式（１２）のための最小二乗誤差解法が次のようにもたらされる。ただし、行列Ａ+はいわゆるムーア・ペンローズ行列または擬似逆行列である。
【数３０】
0005235283.tif 000031
【００４１】
先に既に述べたように、行列Ａの要素が拡散測定のための制御データつまり測定プロトコルから決定可能である。これらのデータは測定値に依存しない。その結果、擬似逆行列Ａ+が本来のデータ取得前に決定される。解ベクトル
【数３１】
0005235283.tif 000032
はピクセルごとに測定データから次のように増分計算される。
【数３２】
0005235283.tif 000033
ただし、ａi+は擬似逆行列Ａ+のｉ番目の列を示す。最後の拡散エンコーディングによる測定が終了した後には、同時に解ベクトルが存在する。
【００４２】
２つの実施形態の上記の基礎説明に続いて、図１の概略図に基づいて、拡散テンソルを決定するための主要構成要素を有する診断用磁気共鳴装置１を説明する。適切な診断用磁気共鳴装置は、例えばシーメンス社から「Ｍａｇｎｅｔｏｎ  Ｓｙｍｐｈｏｎｙ」、「Ｍａｇｎｅｔｏｎ  Ａｖａｎｔｏ」、「Ｍａｇｎｅｔｏｎ  Ｔｒｉｏ」、「Ｍａｇｎｅｔｏｎ  Ｅｓｐｒｅｅ」という製品名で市販されている。診断用磁気共鳴装置１は、高い磁場強度の均一の磁場を発生する静磁場磁石２を含む。磁気共鳴信号の励起および受信のために、静磁場の強度および画像化のための原子核の磁気回転比によって決定される運転周波数を有する高周波システム４が設けられている。医用画像化においては、主に、生体組織内に多数存在する水素原子核またはプロトンが磁気共鳴画像化に使用され、また拡散テンソルを決定するために使用される。傾斜磁場システム６により核スピンが励起中、励起後および受信時に、投入された傾斜磁場により空間エンコーディングされる。相応のソフトウェアを有するコンピュータアーキテクチャによって十分に実現されている制御部８は、検査領域における磁気共鳴信号の発生から、検査領域から再び放射された磁気共鳴信号の受信までの測定経過全体を制御する。
【００４３】
冒頭に既に述べたように、測定シーケンスは、拡散強調されたまたは拡散エンコーディングされた磁気共鳴信号が発生されるように構成されている。磁気共鳴装置の制御部は相応に構成されている。拡散強調に使用される傾斜磁場のための制御データ１０は、拡散エンコーディングの方向および強さに応じて予め与えられる。検査領域から受信され相応の拡散エンコーディングを有する磁気共鳴信号は、ここでは拡散データ１２として評価ユニットへ導かれる。評価ユニットは、以下に更に説明するようにして、測定された拡散データ１２の完全なセットから拡散テンソル１４の成分を決定する。
【００４４】
図２に基づいて、ブロック図にて拡散テンソルの成分を決定する方法の基本構成を説明する。拡散テンソルの成分を決定する方法の主要な特徴は、特定の拡散エンコーディングのための各部分測定から、先行測定の結果にさらに加えられなければならない増分成分が算出されることにある。最後の測定の増分成分の考慮に基づいて直接的に、反転により拡散テンソルの成分を提供する行列が使用できる。
【００４５】
測定方法の開始２０の後に、まず特定の個数の拡散エンコーディング２２の制御データが予め定められる。この方法においては並行して、拡散テンソルの決定すべき成分の安定な解を得るために、拡散エンコーディングに使用された制御データの状態数２４が十分であるかどうかがチェックされる。制御データは、拡散エンコーディングに使用された傾斜磁場パルスの強さ、持続時間および時間的間隔のための設定を含み、これらの設定は磁気共鳴装置によって相応の拡散強調されたパルスシーケンスに変換される。拡散信号の空間エンコーディングは、一般に、励起時、受信時およびその中間における傾斜磁場により行なわれる。各拡散エンコーディングのために求められ空間エンコーディングされた測定データ１２が処理されて拡散データ２６を形成し、この拡散データ２６を用いて拡散中間データ２８を更新し、したがってこれまでに行なわれた全ての拡散エンコーディングを考慮する現在の拡散中間データ３０を作成する。先に基礎において既に説明したように、実際に現在の拡散中間データの精度のための値も算出される。
【００４６】
精度３１が十分である場合、拡散テンソル１４の成分を決定するために、拡散中間データの行列が反転される。
【００４７】
拡散テンソルの成分は、画像表示装置において、例えば拡散楕円面による公知の表示にて、検査領域における位置に応じて表示される。
【図面の簡単な説明】
【００４８】
【図１】拡散テンソルを決定するための主要な構成要素を有する診断用磁気共鳴装置の概略図
【図２】拡散テンソルの成分を決定する方法の基本構成を示すブロック図
【符号の説明】
【００４９】
１        診断用磁気共鳴装置
２        静磁場磁石
４        高周波システム
６        傾斜磁場システム
８        制御部
１０        制御データ
１２        拡散データ
１４        拡散テンソル
２０        方法の開始
２２        拡散エンコーディング
２４        状態数
２６        拡散データ
２８        拡散中間データ
３０        現在の拡散中間データ
３１        現在の拡散中間データの精度
【図１】
0005235283.tif 000034
【図２】
0005235283.tif 000035