JP7218288B2 - サンプルに拡散強調磁気共鳴測定を実行する方法 - Google Patents

サンプルに拡散強調磁気共鳴測定を実行する方法 Download PDF

Info

Publication number
JP7218288B2
JP7218288B2 JP2019522714A JP2019522714A JP7218288B2 JP 7218288 B2 JP7218288 B2 JP 7218288B2 JP 2019522714 A JP2019522714 A JP 2019522714A JP 2019522714 A JP2019522714 A JP 2019522714A JP 7218288 B2 JP7218288 B2 JP 7218288B2
Authority
JP
Japan
Prior art keywords
diffusion
measurement
encoding
signal attenuation
measurements
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.)
Active
Application number
JP2019522714A
Other languages
English (en)
Other versions
JP2019534087A (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.)
Random Walk Imaging AB
Original Assignee
Random Walk Imaging AB
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 Random Walk Imaging AB filed Critical Random Walk Imaging AB
Publication of JP2019534087A publication Critical patent/JP2019534087A/ja
Application granted granted Critical
Publication of JP7218288B2 publication Critical patent/JP7218288B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G01R33/56341—Diffusion imaging
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, 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
    • 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/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
    • G01R33/4833—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices
    • G01R33/4835—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices of multiple slices
    • 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/543—Control of the operation of the MR system, e.g. setting of acquisition parameters prior to or during MR data acquisition, dynamic shimming, use of one or more scout images for scan plane prescription
    • 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/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Signal Processing (AREA)
  • Radiology & Medical Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Vascular Medicine (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Biophysics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Optics & Photonics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Description

本発明概念は、サンプルに拡散強調磁気共鳴測定を実行する方法に関する。
磁気共鳴（ＭＲ）または磁気共鳴撮像（ＭＲＩ）実験において、粒子の運動または拡散に関する情報は、運動または拡散エンコーディング磁場勾配を適用することによってエンコードされ得る。
運動または拡散エンコード化信号は、スピン担持粒子の拡散に関する制限／制約を表し得る、構成コンパートメントの組織微細構造、異方性、形状、およびサイズに関する情報を推測するために用いられ得る。これらは、非コヒーレント流または乱流の特性を調査するためにも用いられ得る。拡散テンソル撮像（ＤＴＩ）によって得られる異方性度（ＦＡ）は、コンパートメントの巨視的配向分散によって混乱する。ＦＡにおいて密接に関連し合う、制限を表し得る拡散テンソルの異方性から、配向分散の影響を分離するために、指向性拡散エンコーディング（１Ｄ）は、単一方向を超えて（２Ｄまたは３Ｄまで）展開するエンコーディングスキームと併用またはこれに代替される必要がある。これらのスキームは、複数の非ゼロ固有値を有する拡散エンコーディング／強調テンソル（１）によって説明することができ、様々な程度まで、混乱させる配向分散の影響を緩和または排除し、コンパートメント（拡散テンソル）異方性に特化した感度を提供することができる。
コンパートメント（拡散テンソル）異方性に対する感度を最大化する、Ｌａｓｉｃ他によるアプローチ（２）は、微視的異方性度（μＦＡ）を定量化するために指向性（１Ｄ）および等方性（３Ｄ）エンコーディングを併用する。等方性エンコーディングはたとえば、ｑベクトルのマジック角スピニング（ｑＭＡＳ）によって実現することができ（３）、一方、指向性エンコーディングに関する拡散強調は、拡散時間ｔdおよびｂ値に関してｑＭＡＳの拡散強調と一致する。拡散エンコーディングの異方性をさらに制御することによって、扁円および扁長コンパートメント（拡散テンソル）形状が区別され得る。Ｅｒｉｋｓｓｏｎ他（４）は、拡散エンコーディングテンソルをサイズおよび形状に関してパラメータ化することによって、粉体平均信号に関する単純な式が得られ、コンパートメント（拡散テンソル）形状を定量化することが可能であることを示している。拡散エンコーディング形状を変化させることは、等方性拡散率の非ゼロ分散を有するシステムのために必要である（２）。３Ｄエンコーディングシーケンスの使用により、異方性および等方性拡散寄与の畳み込みを解くことが可能である（５）。
不均質かつ異方性物質の特徴付けにおける上述した進歩にもかかわらず、サンプルの拡散特性および微視的構造に関する更なる情報の抽出が可能になることが望ましい。
本発明概念の目的は、サンプルの拡散特性および微視的構造に関する更なる情報の抽出を可能にする方法を提供することである。追加または代替の目的は、下記により理解され得る。
本発明概念の態様によると、サンプルに拡散強調磁気共鳴測定を実行する方法が提供され、方法は、
サンプルに拡散強調磁気共鳴測定を実行することを備え、上記測定は、
少なくとも２つの非ゼロ固有値を有する第１の拡散強調テンソル表現Ｂ1を有する第１の拡散エンコーディングシーケンスを用いる第１の測定と、
少なくとも２つの非ゼロ固有値を有する第２の拡散強調テンソル表現Ｂ2を有する第２の拡散エンコーディングシーケンスを用いる第２の測定と、
を含み、
第１のテンソル表現Ｂ1および第２のテンソル表現Ｂ2は、同じ数の非ゼロ固有値を有し、第１のテンソル表現Ｂ1の固有値は、第２のテンソル表現Ｂ2の固有値と一致し、
第１および第２の拡散エンコーディングシーケンスは、一致する平均スペクトルコンテントを示し、かつ異なる程度のスペクトル異方性を示す、ように構成される。
本発明概念は、以下の見識に基づくものである。拡散テンソルは、物理的細孔を大まかに表現し得るが、エンコーディング波形のスペクトル特性（同等に時間的特性）を変化させることによって、より正確な形態学的検出が実現され、時間依存非ガウシアン拡散を調査することができる。拡散エンコーディングシーケンス波形のスペクトル特性は、エンコーディング波形の各空間チャネルの平均スペクトルコンテントを含む。したがって拡散エンコーディングシーケンスの「平均スペクトルコンテント」は、概念的に、エンコーディング波形の「色」とみなされ得る。これは、知覚される色を生成する可視光の異なるスペクトル成分の平均または混合を暗示する。「色」に加えて、拡散エンコーディングシーケンス波形の時間的特性は、スペクトル異方性の程度を含む。「スペクトル異方性の程度」は、概念的に、エンコーディング波形の空間チャネルが分光的にどのように異なるかに依存する特性と考えられ得る。
一致するスペクトルコンテントかつ異なる程度のスペクトル異方性を示す第１および第２の拡散エンコーディングシーケンスによって、サンプルにおける測定の組み合わせは、異方性拡散構造からの制限拡散を感知するように構成される。よって、制限拡散を示す等方性拡散構造による、結果として生じる信号減衰への寄与は、最小化または少なくとも低減され得る。
本明細書で用いられる「拡散」は、本明細書において、サンプル内の粒子の運動のランダムまたは確率過程を意味する。拡散は、熱エネルギー、化学エネルギー、および／または濃度差によって生じるランダム分子運動を含んでよい。拡散は、サンプル内でランダムに配向する微細構造内部の分子の分散流またはインコヒーレント流（すなわち、速度分散を有する流動）を含んでよい。本方法において用いられる拡散エンコーディング磁場勾配シーケンスに起因して、サンプル内のインコヒーレント流の影響もまた信号減衰を生じさせ得る。
サンプルの拡散特性（すなわち、等方性拡散の程度、異方性拡散の程度および／または配向など）は一般に、サンプル内の拡散制限構造の幾何学的形状および配向に依存する。異なる拡散特性を有するサンプルの部分体積は、サンプルの異なるコンパートメントと称され得る。よってサンプルのボクセル（すなわち、測定の空間分解能によってもたらされる次元を有するサンプルの部分体積）からの拡散強調測定信号は、ボクセル内の異なるコンパートメント内部の拡散粒子からの信号寄与を含む。
従来技術において、拡散強調磁気共鳴測定に基づくコンパートメント異方性の分析は、拡散過程がガウシアンまたはマルチガウシアンであるものと仮定する。この仮定は、拡散過程が非時間依存性（または同等に非周波数依存性）であるという仮定に等しい。しかし、発明者によって理解されるように、拡散粒子に関する物理的制限の存在は、周波数依存信号減衰によって明らかにされ得る。
拡散過程は、周波数領域において、拡散スペクトラム、すなわち移動する粒子の速度の相関のスペクトラムによって表され得る。本発明の方法に係る測定プロトコルは、意図的に選択されたスペクトル特性、すなわち一致する平均スペクトルコンテントおよび異なる程度のスペクトル異方性を有する第１および第２のエンコーディングシーケンスを含む。その結果、第１および第２の測定による信号減衰の比較によって、拡散スペクトラム、すなわちサンプルの拡散特性の周波数依存性を調査することが可能である。本方法は、拡散スペクトラムを調査することを可能にし、従来の磁気共鳴測定用ハードウェアによって検出され得る周波数範囲内（一般に１ギガヘルツ未満）の拡散スペクトラムの異なるスペクトル成分に対する感度を提供する。
組織サンプルにおいて、コンパートメントは、組織内の細胞によって形成され得る。したがって、組織サンプル内のボクセルは、細胞内部の拡散による（すなわち細胞内コンパートメントからの）信号寄与、および細胞外の拡散（すなわち細胞外コンパートメントからの）信号寄与を含んでよい。拡散スピン担持粒子は、細胞内外の水分子によって形成され得る。細胞膜は、細胞内および細胞外のコンパートメントを分離するコンパートメント壁部を形成し得る。
以下により理解され得るように、上述した方法の態様において言及される拡散エンコーディングシーケンスは、有効勾配シーケンス、すなわち、磁場勾配シーケンスおよび無線周波数（ＲＦ）シーケンスの組み合わせに起因してサンプル内のスピン担持粒子が感じる有効磁場勾配を指してよい。したがって、特に記載されない限り、拡散エンコーディングシーケンスという用語は、（エコー信号減衰の）拡散エンコーディング／強調をもたらすために適用される磁場勾配シーケンスおよびＲＦシーケンスを含むエンコーディングシーケンスを指す。
拡散エンコーディングは、サンプルに、エンコーディング磁場勾配波形およびＲＦパルスのシーケンスを受けさせることによって実現され得る。磁場勾配波形およびＲＦパルスの影響が結合された結果、サンプル（内のスピン担持粒子）は「有効勾配」を受ける。有効勾配の波形は、成分ｇi（ｔ）を有する有効勾配波形ｇ（ｔ）と称されてよく、ここでｉ＝１、２、３であり、たとえばデカルト座標系におけるｘ、ｙ、ｚ軸を表す。
有効勾配は、成分Ｆiを有する時間依存または時間ディフェージングベクトルＦ（ｔ）によって表されてよく、ここでｉ＝１、２、３であり、Ｆ（ｔ）は、
Figure 0007218288000001
によって求められ、式中、γは核ジャイロ磁気比である。
ディフェージングスペクトル、すなわちディフェージングベクトルのスペクトルコンテントは、
Figure 0007218288000002
によって求められ、式中、ωは周波数を示し、τは拡散エンコーディング時間、すなわち有効勾配の持続期間を示す。
ディフェージングスペクトラムＦ（ω）に基づいて、ディフェージングスペクトラムのモーメントは、テンソルＭ(n)として表される、
Figure 0007218288000003
と定義されてよく、式中、Ｍij (n)はテンソル要素である。０番目のモーメントＭ(0)によって拡散強調テンソルＢが与えられる。ガウシアン拡散の場合、Ｂij＝Ｍij (0)のみを考える必要がある。
したがって、拡散エンコーディング磁気勾配シーケンスの拡散強調テンソル表現Ｂは、
Figure 0007218288000004
によって求められる。
ディフェージングベクトルＦ（ｔ）は、ディフェージングベクトルの振幅ｑと正規化ディフェージングベクトル波形Ｆ~（ｔ）とのアダマール積
Figure 0007218288000005
として表されてもよい。
したがってディフェージングスペクトラムは、また／あるいは、
Figure 0007218288000006
によって求められてもよく、式中、Ｆ~（ω）は正規化ディフェージングスペクトラム、すなわち正規化ディフェージングベクトルのスペクトルコンテントとみなされ得る。デカルト座標系におけるｘ、ｙ、ｚ軸を表すｉ＝１、２、３であるスペクトラムＦi（ω）およびＦ~ i（ω）は、それぞれディフェージングベクトルまたは正規化ディフェージングベクトルのスペクトル表現の直交成分または正射影として言及され得る。
式（６）における正規化スペクトラムＦ~（ω）に関して、テンソルＭ(n)の要素は、
Figure 0007218288000007
によって求められる。項ｍij (0)は、一般化拡散時間に対応する。マトリックス形式で表すと、
Figure 0007218288000008
である。
１Ｄエンコーディングにおいて、ディフェージング波形は、３Ｄエンコーディングのディフェージングの大きさによって構成され得るので、１Ｄエンコーディングは、拡散時間に関して３Ｄエンコーディングと一致し、たとえば１Ｄディフェージング振幅は、３Ｄエンコーディングの二乗ディフェージング振幅成分の和、
Figure 0007218288000009
によって求められる。
拡散スペクトラムＤ（ω）は、速度相関テンソル
Figure 0007218288000010
のスペクトルとして定義され得る。
コンパートメントの主軸系において、Ｄ（ω）は対角線であり、マトリックスλ（ω）の対角線に沿った拡散スペクトラムλi（ω）によって求められる。要素λi（ω）は、拡散テンソル固有値とみなされ得る。よって見かけ上の拡散率は、
Figure 0007218288000011
によって求められる。上記から理解され得るように、（有効）拡散エンコーディング磁気勾配シーケンスのスペクトルコンテント、または同等にディフェージングベクトル波形は、ディフェージングスペクトラムＦ（ω）またはその正規化部分Ｆ~（ω）によって求められる。
正規化パワースペクトラムの成分は、
Figure 0007218288000012
によって求められる。正規化パワースペクトラムのｎ番目のモーメントは、テンソルまたはマトリックス形式ｍ(n)で表されてよく、
Figure 0007218288000013
であり、式中、ｍij (n)はｍ(n)の要素を示す。μi (n)によってｍ(n)の固有値が示される。
テンソルＭ(n)のテンソル要素は、
Figure 0007218288000014
によって求められ得る。したがって、拡散強調テンソルＢの要素は、
Figure 0007218288000015
によって求められ得る。（単一の多次元エンコーディングスキームに関する）平均エンコーディングスペクトルコンテントは、
Figure 0007218288000016
によって特徴付けることができ、式中、ｎは任意の正の実数であり、
Figure 0007218288000017
である。鍵括弧〈…〉によって、平均化演算が示される。任意のｎ＞０に関する〈μ(n)〉の差は、２つのエンコーディングスキームを比較するために用いられ得る。異なる拡散エンコーディングシーケンス（または同等に異なる正規化ディフェージングベクトル表現）のスペクトルの一致またはスペクトル同調は、一致する〈μ(n)〉を有することを指してよい。逆に、異なるエンコーディングシーケンスは、それらのスペクトルコンテントが一致しない、すなわち〈μ(n)〉が異なる場合、離調するものとみなされ得る。
（単一の多次元エンコーディングスキームに関する）拡散エンコーディングのスペクトル異方性は、ｍ(n)の異なる固有値、すなわち値μi (n)を比較することによって定量化され得る。テンソル異方性の任意の計測、たとえば異方性度が、テンソルｍ(n)に適用され得る。スペクトル異方性（ＳＡ）を定量化するために、異方性度と同様に、以下の式
Figure 0007218288000018
が（拡散エンコーディングシーケンスの拡散エンコーディングテンソル表現が複数の非ゼロ固有値を有するという条件で）用いられてよく、式中、ｎは任意の正の実数、
Figure 0007218288000019
であり、スペクトル異方性に対する感度を制御するために調整され得る。ただし、ＳＡ(n)は常に０～１の範囲内である。テンソル異方性の定量化に適用可能な、スペクトル異方性の他の定量的計測が用いられてもよいことに留意する。
上記において、第１および第２の拡散エンコーディングシーケンスは、ある実施形態によると、任意のｎ＞０、好適には少なくともｎ＝２の場合、第１の拡散エンコーディングシーケンスに関して計算された〈μ(n)〉が第２の拡散エンコーディングシーケンスに関して計算された〈μ(n)〉と一致するという意味で、一致する平均スペクトルコンテントを示し得る。
さらに上記観点から、第１および第２の拡散エンコーディングシーケンスは、ある実施形態によると、第１の拡散エンコーディングシーケンスに関して計算されたＳＡ(n)が第２の拡散エンコーディングシーケンスに関して計算されたＳＡ(n)と異なるという意味で、異なる程度のスペクトル異方性を示し得る。
よって、第１の拡散エンコーディング勾配シーケンスに関して、
Figure 0007218288000020
であり、式中、μi (n)は、要素
Figure 0007218288000021
を有するテンソルｍ(n)の固有値を示し、
Figure 0007218288000022
はそれぞれ、第１の拡散エンコーディング勾配シーケンスの正規化ディフェージングスペクトラムＦ~（ω）の正規化ディフェージングベクトルのｉ番目およびｊ番目の成分である。
また、第２の拡散エンコーディング勾配シーケンスに関して、
Figure 0007218288000023
であり、式中、μi (n)は、要素
Figure 0007218288000024
を有するテンソルｍ(n)の固有値を示し、Ｆ~ i（ω）およびＦ~* i（ω）はそれぞれ、第２の拡散エンコーディング勾配シーケンスの正規化ディフェージングスペクトラムＦ~（ω）のｉ番目およびｊ番目の固有ベクトルである。
一致するまたは異なる「スペクトルコンテント」の概念は、拡散エンコーディングシーケンス「Ａ」を用いる「第１の試験サンプル」における測定「Ｉ」、および拡散エンコーディングシーケンス「Ｂ」を用いる第１の試験サンプルにおける他の測定「ＩＩ」を考慮することによっても理解され得る。第１の試験サンプルは、直径５μｍの中空球状コンパートメントの集合から成るサンプルである。
シーケンス「Ａ」およびシーケンス「Ｂ」は、シーケンス「Ａ」および「Ｂ」のテンソル表現ＢA、ＢBが互いに等しく（すなわち、それらが同じ拡散エンコーディング強度をもたらし）、同一の固有値を有するようなものである。コンパートメントは球状であるため、主軸系および固有ベクトルは縮退する。よって、シーケンス「Ａ」および「Ｂ」のテンソル表現ＢA、ＢBは、コンパートメントに対する任意の配向を有する３次元デカルト座標系（ｘ，ｙ，ｘ）に関して与えられ得る。
シーケンス「Ａ」およびシーケンス「Ｂ」が試験サンプルに適用されると同じレベルの信号減衰をもたらす場合、シーケンスは、「分光的に一致」しており、すなわち一致する平均スペクトルコンテントを有する。
一方、シーケンス「Ａ」およびシーケンス「Ｂ」が試験サンプルに適用されると異なるレベルの信号減衰をもたらす場合、シーケンスは、「分光的に離調」しており、すなわち異なる平均スペクトルコンテントを有する。
対応して、一致するまたは異なる「スペクトル異方性の程度」の概念は、上記で定義した拡散エンコーディングシーケンス「Ａ」を用いる「第２の試験サンプル」における測定「ＩＩＩ」、および上記で定義した拡散エンコーディングシーケンス「Ｂ」を用いる第２の試験サンプルにおける他の測定「ＩＶ」によって理解され得る。第２の試験サンプルは、均一な配向分散を有する直径５μｍの中空円筒状コンパートメントの集合から成るサンプルである。「均一な配向分散」は、ここでは、円筒状コンパートメントの配向が単位球面にわたり均一に分布しているものとして理解すべきである。
シーケンス「Ａ」および「Ｂ」が分光的に一致していると仮定すると、シーケンス「Ａ」およびシーケンス「Ｂ」が試験サンプルに適用されると同じレベルの信号減衰をもたらす場合、シーケンスは、同じ程度のスペクトル異方性を示す。一方、シーケンス「Ａ」およびシーケンス「Ｂ」が試験サンプルに適用されると異なるレベルの信号減衰をもたらす場合、シーケンスは、異なる程度のスペクトル異方性を示す。
上記観点から、ある実施形態によると、第１および第２の拡散エンコーディングシーケンスは、
第１の拡散エンコーディングシーケンスの第１の正規化ディフェージングベクトル表現Ｆ~ 1と一致する正規化ディフェージングベクトル表現Ｆ~ 3を有し、非ゼロ拡散エンコーディング強度を有する第３の拡散エンコーディングシーケンスが、直径５μｍの球状コンパートメントの集合から成る第１の試験サンプルに適用され、
第２の拡散エンコーディングシーケンスの第２の正規化ディフェージングベクトル表現Ｆ~ 2と一致する正規化ディフェージングベクトル表現Ｆ~ 4を有し、上記非ゼロ拡散エンコーディング強度を有する第４の拡散エンコーディングシーケンスが、上記第１の試験サンプルに適用され、
第３の拡散エンコーディングシーケンスの結果生じる信号減衰が、第４の拡散エンコーディングシーケンスの結果生じる信号減衰と一致するように構成され、かつ、
第３の拡散エンコーディングシーケンスが、均一な配向分散を有する直径５μｍの円筒状コンパートメントの集合から成る第２の試験サンプルに適用され、
第４の拡散エンコーディングシーケンスが、上記第２の試験サンプルに適用され、
第３の拡散エンコーディングシーケンスの結果生じる信号減衰が、第４の拡散エンコーディングシーケンスの結果生じる信号減衰と異なるように構成される。
「第３／第４の拡散エンコーディングシーケンスが…試験サンプルに適用されるように構成される」という表現は、特許請求の範囲に記載の方法において積極的に実行されることが必ず必要なステップのシーケンスとして解釈されてはならない。むしろこの表現は、第１および第２の拡散エンコーディングシーケンスの特性の機能的定義として理解されるべきである。よって、第３および第４の拡散エンコーディングシーケンスが試験サンプルに適用された場合／条件において、上述の信号減衰を伴う。この機能的定義は、当業者が、拡散エンコーディングシーケンスのセットが進歩性のある特性を有するかを決定するための明瞭かつ明確なテストケースとして理解され得る。実際、特性を試験するために現実の「試験サンプル」に実際の測定を実行することは必要ではない。むしろテストケースは、第１および第２の拡散エンコーディングシーケンスのそれぞれのＲＦシーケンスおよび磁場勾配シーケンスを測定することによって、および第１および第２の拡散エンコーディングシーケンスのテンソル表現および正規化ディフェージングベクトル表現を計算する測定から評価されてよく、その後、（上述したエンコーティング強度および正規化ディフェージングベクトル表現を有する）第３および第４の拡散エンコーディングシーケンスが第１および第２の試験用物体に適用された場合、その結果生じる信号減衰を（数値計算またはモンテカルロシミュレーションによって）シミュレートする。
上記および下記において、互いに一致する２つ以上の数量（たとえば、２つ以上の拡散エンコーディング強度、トレース、固有値、〈μ(n)〉またはＳＡ(n)）が述べられる。そのような文脈において、「一致する」は、数量が等しいこと、または少なくとも実質的に等しいものとして理解され得る。すなわち、数量は厳密に等しい必要はないが、コンパートメント異方性、スペクトルコンテント、およびスペクトル異方性の影響が測定から識別可能であるような量のみ異なるべきである。数量は、好適には１０％未満、より好適には５％未満、さらに好適には１％未満異なるべきである。理解され得るように、同等性の実現可能な最大レベルは、実際には、とりわけ機器の性能に依存する。
拡散エンコーディングシーケンスの「拡散エンコーディング強度」は、本明細書において理解されるように、拡散エンコーディングシーケンスに関する拡散エンコーディングテンソルのトレースによって求められるｂ値を指す。
本方法によると、拡散強調磁気共鳴測定が実行される。測定は、少なくとも２つの（すなわち、「第１」および「第２」の）測定を含み、各測定は、サンプルに拡散エンコーディングシーケンスを受けさせることを含む。本明細書で用いられる場合、測定／エンコーディングシーケンスの「第１」および「第２」というラベルは、測定がその特定の順序で実行されることを暗示するものではなく、測定は任意の順序で実行されてよい。
上記拡散強調磁気共鳴測定の各々は一般に、エンコーディングブロックおよび後続する検出ブロックを含んでよい。エンコーディングブロックにおいて、サンプルは、拡散エンコーディングシーケンスを受ける。検出ブロックにおいて、拡散エンコーディングに起因して減衰された信号が検出および取得され得る。検出された信号または測定信号は、減衰エコー信号であってよい。
エンコーディングブロックはさらに、サンプル内の磁化を促すために適合された無線周波数（ＲＦ）パルスシーケンスを含んでよい。ＲＦパルスシーケンスは、縦方向緩和のみ、横方向緩和のみ、または縦方向緩和および横方向緩和の両方に起因する減衰をエンコード化してよい。よって、検出ブロックにおいて検出された信号の減衰は、拡散エンコーディング磁気勾配パルスシーケンスおよびＲＦパルスシーケンスの両方に起因する減衰、すなわち、有効勾配シーケンスｇ（ｔ）または対応するディフェージングベクトルＦ（ｔ）に起因する減衰の結果である。
好適には、上記複数の測定は、同一のタイミングを有するＲＦパルスシーケンスを用いて実行される。すなわち、核緩和に起因する同じレベルの信号減衰が、各測定においてエンコード化される。これによって、測定に影響を及ぼす変化するパラメータの数が低減され得るので、データ分析が単純化され得る。
（減衰された）測定信号を取得するために、各拡散エンコーディングシーケンスは、当該技術において周知であるように、１または複数のイメージング磁気勾配および任意選択的に磁気勾配補正勾配によって補足され得る。イメージング磁気勾配シーケンスおよび補正磁気勾配シーケンスは、エンコーディングブロック中にサンプルに適用され得る。場合によっては、これらのシーケンスは、拡散エンコーディング磁気勾配パルスシーケンスと少なくとも部分的に重なってよい。ただしそのような場合、複合勾配パルスシーケンスの少なくとも一部は、上述したように説明または特徴付けされ得る拡散エンコーディングシーケンスを含む。
ある実施形態によると、上記第１の測定は、上記第１の拡散エンコーディングシーケンスの結果生じる第１の信号減衰を取得することを含み、
上記第２の測定は、上記第２の拡散エンコーディングシーケンスの結果生じる第２の信号減衰を取得することを含む。
ある実施形態によると、上記測定は、
上記第１の測定と、参照測定フレームに対して異なる回転によるサンプルに適用される上記第１の拡散エンコーディングシーケンスを用いて実行される複数の追加の測定と、を含む第１の測定セットと、
上記第２の測定と、参照測定フレームに対して異なる回転によるサンプルに適用される上記第２の拡散エンコーディングシーケンスを用いて実行される複数の追加の測定と、を含む第２の測定セットと、
を含む。
これによって、（「粉体平均化」としても知られる）方向性平均化信号減衰を決定することが可能である。
したがって、ある実施形態によると、
上記第１の測定セットの各測定は、上記拡散エンコーディングシーケンスの結果生じるそれぞれの信号減衰を取得することを含み、方法は、上記それぞれの信号減衰に基づいて第１の平均信号減衰を決定することを備え、
上記第２の測定セットの各測定は、上記拡散エンコーディングシーケンスの結果生じるそれぞれの信号減衰を取得することを含み、方法は、上記それぞれの信号減衰に基づいて第２の平均信号減衰を決定することを備える。
ある実施形態によると、方法はさらに、第１および第２の信号減衰の差または第１および第２の平均信号減衰の差を示す出力を生成することを備える。よって出力は、時間依存異方性拡散、すなわち異方性領域内の時間依存拡散に起因する信号減衰を感知可能であってよい。
ある実施形態によると、上記第１のテンソル表現は、３つの一致する非ゼロ固有値を有する。
ある実施形態によると、上記第２のテンソル表現は、３つの一致する非ゼロ固有値を有する。
ある実施形態によると、上記第１のテンソル表現は、厳密に２つの一致する非ゼロ固有値を有する。
ある実施形態によると、上記第２のテンソル表現は、厳密に２つの一致する非ゼロ固有値を有する。
ある実施形態によると、上記測定を実行することは、サンプルの関心領域内の複数のボクセルの各々から、上記測定の各々の結果生じるそれぞれの信号減衰を測定することを含む。
ある実施形態によると、方法はさらに、第１の測定において取得された信号減衰が第２の測定において取得された信号減衰と異なるボクセルの指標を含む出力を生成することを備える。
本発明概念の第２の態様によると、サンプルに拡散強調磁気共鳴測定を実行する方法が提供され、方法は、
サンプルに拡散強調磁気共鳴測定を実行することを備え、
上記測定は、
異なる拡散エンコーディング強度および一致する平均スペクトルコンテントおよび一致する程度のスペクトル異方性を有する拡散エンコーディングシーケンスを用いて実行される第１の測定セットと、
異なる拡散エンコーディング強度および一致する平均スペクトルコンテントおよび一致する程度のスペクトル異方性を有する拡散エンコーディングシーケンスを用いて実行される第２の測定セットと、
を含み、
第１のセットの拡散エンコーディングシーケンスのスペクトル異方性の程度は、第２のセットの拡散エンコーディングシーケンスのスペクトル異方性の程度と異なる。
第１の態様に関連する利点および詳細な説明は、対応して第２の態様にも適用される。
一致する平均スペクトルコンテントおよび一致する程度のスペクトル異方性の上述した数学的かつテストベースの定義は、対応して、第１および第２のセットのエンコーディングシーケンスの各々にも適用される。
ある実施形態によると、方法はさらに、
第１の信号減衰曲線を推定するために上記第１の測定セットを表す第１のデータセットに第１の関数をフィッティングすることと、
第２の信号減衰曲線を推定するために上記第２の測定セットを表す第２のデータセットに第２の関数をフィッティングすることと、
を備える。
第１および第２の測定セットの各々に対応するそれぞれの信号減衰曲線を推定することによって、拡散特性およびそれらの周波数依存性の更なる分析が可能である。同じフィッティング関数が、第１および第２のデータセットへのフィッティングに用いられ得る。
ある実施形態によると、方法はさらに、第１の関数の少なくとも１つのパラメータ、第２の関数の少なくとも１つのパラメータ、および第３の関数の少なくとも１つのパラメータに基づいて出力を生成することを備える。
ある実施形態によると、第１のセットおよび第２のセットの各々の拡散エンコーディングシーケンスはそれぞれのテンソル表現を有し、テンソル表現の非ゼロ固有値の数は等しい。
ある実施形態によると、第１のセットの拡散エンコーディングシーケンスの各々の平均スペクトルコンテントは、第２のセットの拡散エンコーディングシーケンスの各々の平均スペクトルコンテントと一致する。
本発明概念の上記ならびに追加の目的、特徴、および利点は、添付図面を参照して、以下に示す本発明概念の好適な実施形態の例示的かつ非限定的な詳細な説明によってより良く理解される。図面において、特に記載されない限り、同様の要素に関して同様の参照番号が用いられる。
２つの拡散エンコーディングシーケンスに関するディフェージング波形およびパワースペクトラム、および対応する減衰曲線のシミュレーションを示す。 図２は、サンプルに拡散強調磁気共鳴測定を実行する方法のフローチャートを示す。 図３ａおよび図３ｂは、パルスシーケンスの例を示す。
詳細な説明
本発明概念の理解を容易にするために、以下、いくつかの理論的概念の説明が提供される。
理論
拡散エンコーディング時間τにおける任意の有効な勾配波形ｇ（ｔ）を考えると、信号は、アンサンブル平均
Figure 0007218288000025
によって求められ、式中、γは核ジャイロ磁気比であり、Δｒ（ｔ）は変位である。複数のコンパートメントシステムが考慮される場合、アンサンブル平均は、異なる拡散特性、たとえば孔径、形状、および配向を有するコンパートメントであるサブアンサンブルにわたって行われる。低い拡散エンコーディング勾配の極限において、信号減衰は、ガウシアン位相近似（ＧＰＡ）に対応して、二次キュムラント展開（６）によって近似され得る。
拡散エンコーディング中の変位相関性が無視できる、見かけ上の拡散テンソルＤを特徴とするガウシアン拡散過程を有する単一コンパートメントまたはサブアンサンブルの場合、信号減衰は、Ｅ＝ｅｘｐ（－β）によって求められ、式中、
Figure 0007218288000026
であり、Ｆ（ｔ）は、時間ディフェージングベクトル
Figure 0007218288000027
であり、式中、γは核ジャイロ磁気比であり、ｇ（ｔ）は拡散エンコーディング勾配シーケンスの波形である。
拡散テンソルＤは、コンパートメント異方性に関する情報を伝達する。ディフェージングベクトルＦ（ｔ）は、振幅ｑと正規化波形Ｆ~（ｔ）とのアダマール積、
Figure 0007218288000028
として表すことができる。
非ガウシアン拡散の場合、および信号キュムラント展開（２次キュムラント展開）のガウシアン近似の場合、減衰は、
Figure 0007218288000029
によって求められ、ここで
Figure 0007218288000030
であり、Ｄ（ω）は拡散スペクトラム、すなわち速度相関テンソル
Figure 0007218288000031
のスペクトラムである。周波数領域分析は、任意の勾配波形に容易に適用され得る。最も重要な点として、これによって、３Ｄ拡散エンコーディングにおける制限拡散効果の直感的な理解が可能であり、異なるエンコーディング波形の設計が容易になり得る。
アインシュタインの縮約記法（すなわち、繰り返す指数を総和すること）によると、
Figure 0007218288000032
である。コンパートメント／制限の主軸系（ＰＡＳＣ）において（７）、拡散スペクトラムは、
Figure 0007218288000033
によって求められ、式中、λ（ω）は対角行列であり、Ｒは回転行列である。またディフェージングスペクトラムは、
Figure 0007218288000034
としても求められる。表記を簡略化するために、積分演算子
Figure 0007218288000035
を
Figure 0007218288000036
と定義し、ここで、ｉ，ｊ，ｋ∈１，２，３であり、上付き記号＊は複素共役化を示す。ただし、
Figure 0007218288000037
である。アインシュタインの縮約記法を取り入れ、式（２２）における定義を用いると、式（１８）は、
Figure 0007218288000038
となる。減衰（２３）は、見かけ上の拡散率
Figure 0007218288000039
、および、Ｂテンソル
Figure 0007218288000040
を
Figure 0007218288000041
として、表され得る。
ガウシアン拡散、すなわち非時間／周波数依存性の場合、
Figure 0007218288000042
、かつ
Figure 0007218288000043
である。
Ｄij＝ＲkiＲkjλkを示すことによって、式（２７）は、テンソルＢおよびＤの内積〈・，・〉
Figure 0007218288000044
と書き換えることができる。λi（ω）に関する異なるモデルが適用され得る。低周波数において、λk（ω）は、
Figure 0007218288000045
としてテイラー展開されてよく、式中、ｎは整数であり、λk (n)（０）は、ゼロ周波数におけるλk（ω）のｎ番目の導関数である。制限拡散の場合、
Figure 0007218288000046
である。この場合、式（２９）は、制限サイズに関して表される等しいパワーのωおよびλi (n)（０）に関して書き表され得る（７）。
式（２９）を式（２３）に挿入することにより、
Figure 0007218288000047
が導かれ、ここで、
Figure 0007218288000048
は、テンソルＭ(n)として表されるディフェージングスペクトラムのモーメントであり、Ｍij (n)はテンソル要素である。ガウシアン拡散の場合、Ｂij＝Ｍij (0)のみを考慮する必要がある。
正規化スペクトラムＦ~（ω）に関して、テンソルＭ(n)は、
Figure 0007218288000049
によって求められる。
スペクトルコンテント
拡散エンコーディングシーケンスのスペクトルコンテントは、ディフェージングスペクトラムまたは正規化部分に関して特徴付けられ得る。正規化パワースペクトルの成分は、
Figure 0007218288000050
によって求められ、正規化パワースペクトラムのモーメントは、
Figure 0007218288000051
である。ただし、Ｍij (0)＝０である場合、ｍij (0)＝０であり、Ｍij (0)＞０である場合、ｍij (0)＝１である。μi (n)によって、ｍ(n)の固有値が示される。
（単一の多次元エンコーディングスキームに関する）平均エンコーディングスペクトルコンテントは、
Figure 0007218288000052
によって特徴付けられてよく、式中、ｎは任意の正の実数であり、
Figure 0007218288000053
である。鍵括弧〈…〉によって、平均化演算が示される。任意のｎ＞０に関する〈μ(n)〉の差は、２つのエンコーディングスキームを比較するために用いられ得る。式（３６）における分母の値は、拡散エンコーディングテンソルＢの非ゼロ固有値の総数を提供する。
スペクトル異方性
（単一の多次元エンコーディングスキームに関する）拡散エンコーディングのスペクトル異方性は、ｍ(n)の異なる固有値、すなわち値μi (n)を比較することによって定量化され得る。テンソル異方性の様々な計測、たとえば異方性度が、テンソルｍ(n)に適用され得る。ただし、スペクトル異方性は、拡散エンコーディングテンソルの複数の非ゼロ固有値を有するエンコーディングスキームに関してのみ定義され得る。スペクトル異方性（ＳＡ）を定量化するために、異方性度と同様に、以下の式
Figure 0007218288000054
が用いられてよく、式中、ｎは任意の正の実数、
Figure 0007218288000055
であり、スペクトル異方性に対する感度を制御するために調整され得る。ただし、ＳＡ(n)は常に０～１の範囲内である。
一致するスペクトルコンテントまたは多次元拡散エンコーディングスキームおよび指向性エンコーディングスキームの「同調」は、数値的に行われ得る。以下の手順は、多次元エンコーディング勾配波形が、何らかの先行する最適化に基づいて得られている場合、有用である。１．入力波形の（たとえばｘ、ｙ、ｚ軸に沿った）回転の範囲に関して、積
Figure 0007218288000056
を計算する、ここでｎの選択は、「同調」、すなわち一致するスペクトルパワーにおいて優先することが選択され得る周波数範囲に依存するものとし、いくつかのモーメントが計算され得る。２．上記ステップにおいて計算されたｐの最大値、およびｐの最大値を生じる回転Ｒを求める。３．回転Ｒを用いて入力波形を変換する。４．変換された波形の各々について、モーメントｍ(n)を計算し、モーメントｍ(n)が平均モーメントに最も近くなる波形を選択する（２つまたは３つの波形および対応するモーメントを考慮する）。５．ステップ４において選択された波形の形状を指向性エンコーディングのために用いる。
粉体平均信号
式（３１）における
Figure 0007218288000057
を定義すると、これは、テンソルＭ(n)とＤ(n)との内積
Figure 0007218288000058
として書き換えることもでき、ここで、
Figure 0007218288000059
である。２次（Ｍ(n)のパワー）に至るまで、粉体（方向性）平均信号のキュムラント展開によって
Figure 0007218288000060
が導かれ、ここで、
Figure 0007218288000061
、および
Figure 0007218288000062
である。内積表記によって表すと、
Figure 0007218288000063
および
Figure 0007218288000064
であり、ここで、
Figure 0007218288000065
および
Figure 0007218288000066
である。
減衰項（４３）は、拡散尖度または見かけ上の拡散係数の分布の分散に関連する、信号減衰曲線の曲率を定量化するために用いられ得る。ただし、定義（４４）および（４５）は、指数ｉおよびｊの順列の下で不変である。ｉ＝０および２の項（制限拡散に関する近似）のみを考えると、
Figure 0007218288000067
および
Figure 0007218288000068
である。
ただし、参照文献（１０）における式３および６を参照すると、
Figure 0007218288000069
は、ガウシアン項に対応する。もし、Ｄkl (0)およびＤkl (2)が相関しない場合、すなわち（長期の）屈曲の極限における拡散率が時間的変化（低周波数における傾斜）と無関係である場合、
Figure 0007218288000070
である。
有効制限サイズとの関連性：制限拡散の場合
式（２９）における制限拡散スペクトラムの第２の拡散係数は、
Figure 0007218288000071
によって求められ（７）、式中、Ｃiは幾何学的係数であり、γiは、主軸ｉに沿った制限のサイズである。有効制限テンソルを定義すると、
Figure 0007218288000072
である。
よって式（４６）は、
Figure 0007218288000073
となる。
λk (0)（０）が方向と無関係かつＤ0に等しい場合、上の第１項は、
Figure 0007218288000074
によって求められる。κ2は、以下の
Figure 0007218288000075
Figure 0007218288000076
および
Figure 0007218288000077
を考慮して評価され得る。
Ｄkl (0)およびＤkl (2)が相関しない場合、
Figure 0007218288000078
および
Figure 0007218288000079
である。
等方性拡散強調の場合
定義
Figure 0007218288000080
を考えると、
Figure 0007218288000081
である。
Figure 0007218288000082
Figure 0007218288000083
式中、
Figure 0007218288000084
はテンソル積を示す。式（５６）における第１項は、ゼロ周波数（長時間体制）における等方性拡散率の分散、すなわちテンソルトレースの分散に起因し、その他の項は、コンパートメントのサイズおよび異方性に起因する分散を表す。
スペクトル異方性を有さない等方性拡散強調の場合
エンコーディングが分光的に等方性である場合、
Figure 0007218288000085
Figure 0007218288000086
および
Figure 0007218288000087
である。
式（５６）と式（５８）とを比較すると、留意すべき点として、分光的に等方性のエンコーディングの場合、コンパートメント異方性に起因する追加の分散が排除され、コンパートメントサイズに起因する分散のみが残る。スペクトル異方性の変更は、異方性コンパートメントに関連する時間依存拡散効果を分離するために用いられ得る。
軸対称性
単一コンパートメントの場合
ここで、簡略な表記Ｆ＝Ｆ（ω）およびλ＝λ（ω）を取り入れる。制限／コンパートメントの主軸系（ＰＡＳＣ）における拡散スペクトラム
Figure 0007218288000088
および拡散スペクトラムを有する軸対称拡散エンコーディングの場合、
Figure 0007218288000089
であり、式（２３）における減衰によって、
Figure 0007218288000090
が生じ、ここで、Ｐ2（ｘ）＝（３ｘ2－１）／２は第２のルジャンドル多項式であり、
Figure 0007218288000091
は、拡散エンコーディングおよび拡散テンソルの主対称軸間の角度である。軸対称拡散エンコーディングは、ここでは厳密に、
Figure 0007218288000092
の対称性を指す。ガウシアン（非時間／周波数依存性）拡散の場合、ディフェージングスペクトラムのクロス積を含む式（６０）における最後の項は、
Figure 0007218288000093
である場合、ゼロになる。プランシュレルの定理によると、上記条件は、
Figure 0007218288000094
に変換され、これは、ベクトルＦ（ｔ）が直円錐面に常に平行であり、ｑ軌跡が少なくとも三回対称を有する場合（３、４）、または単純に常に積
Figure 0007218288000095
である場合、満たされる。低いコンパートメント異方性に関して式（６０）における最後の項はゼロになることは明らかである。ｑ－ＭＡＳ波形（２、１１）を用いる数値計算は、たとえば円筒から成る大部分が異方性のコンパートメントに関して、この項が比較的小さいことを示す。
クロス積の項がない場合、式（６０）は、
Figure 0007218288000096
のように簡略化することができ、ここで、
Figure 0007218288000097
および
Figure 0007218288000098
である。ここで、表記
Figure 0007218288000099
および
Figure 0007218288000100
が用いられ、
Figure 0007218288000101
である。Λijは、拡散エンコーディング波形ｉに起因するＰＡＳＣにおける軸ｊに沿った見かけ上の拡散係数である。
サンプルにおける全てのコンパートメントが同一である場合、粉体／方向性平均〈ｅ-β〉は、
Figure 0007218288000102
によって求められ、式中、
Figure 0007218288000103
である。単一コンパートメントの場合に平均拡散率（ＭＤ）に対応する等方性拡散率（Ｄ）は、
Figure 0007218288000104
によって求められる。ただし、２つの波形が同一のスペクトルコンテントを有する、すなわち
Figure 0007218288000105
であり、式中、
Figure 0007218288000106
である。したがって、Δｂによって拡散エンコーディングの形状を変化させることにより、ゼロ周波数における拡散異方性Ｄ(0)を定量化することが可能である（４）。
コンパートメント異方性は、見かけ上の拡散係数の分散ももたらし、これは、式（６３）におけるβ-による時間依存拡散の影響も受ける。ここでＶAと示される、コンパートメント異方性に起因する見かけ上の拡散分布の第２の中心モーメント（２、４）は、
Figure 0007218288000107
によって求められ、ここで、
Figure 0007218288000108
である。
よって粉体平均信号は、
Figure 0007218288000109
によって求められる。
Figure 0007218288000110
を有する等方性エンコーディングの場合、
Figure 0007218288000111
であり、
Figure 0007218288000112
である指向性エンコーディングの場合、
Figure 0007218288000113
である。等方性拡散率の分散がない、
Figure 0007218288000114
である純粋に等方性の拡散の場合、
Figure 0007218288000115
である。ただし、異方性制限／コンパートメントがある場合のみ、ＶA＞０である。最も重要な点として、拡散エンコーディングが分光的に等方性である場合のみ、異方性制限に関して
Figure 0007218288000116
であるが、分光的に異方性のエンコーディングの場合、
Figure 0007218288000117
である。したがって、変化するスペクトル異方性により、異方性コンパートメントにおける時間／周波数依存拡散効果を検出することが可能であり得る。一方、コンパートメント異方性に関するバイアスのない情報、すなわち時間／周波数依存拡散効果によって混乱しない情報は、異なる形状のＢ＝Ｍ(0)を有する拡散エンコーディング波形が分光的に同調する場合、すなわち十分な程度類似した平均スペクトルコンテントを有する場合、得られ得る。Ｂを生じるために用いられる異なる波形は、自己同調であってよく、すなわち低いスペクトル異方性を有してもよい。
複数のコンパートメント
各々が見かけ上の等方性拡散率Ｄおよび形状異方性に起因する見かけ上の拡散分散（異方性分散）ＶAによって特徴付けられる、複数のコンパートメントからの粉体平均化信号を考える。「見かけ上の」とは、Ｂの異方性およびエンコーディング波形の時間またはスペクトル特性に依存する値を意味する。合計粉体平均信号は、
Figure 0007218288000118
によって求められ、式中、
Figure 0007218288000119
である。
合計分散Ｖは、等方性分散ＶI＝〈Ｄ2〉－〈Ｄ〉2と平均異方性分散〈ＶA〉との和によって求められる。ただし、ＶIは、等方性コンパートメントおよび異方性コンパートメントの両方からのサイズ分散効果を備える。分光的等方性でもある等方性エンコーディングテンソルＢの場合、〈ＶA〉＝０である。上述したように、等方性Ｂを有する分光的異方性エンコーディングの場合、〈ＶA〉＞０であるように残差分散が予想される。
異方性コンパートメントにおける時間／周波数依存拡散に特化した感度
スペクトル異方性がどのように異方性コンパートメントにおける時間／周波数依存拡散効果の検出に特化した感度を提供し得るかを説明するために、ο［ｂ3］と示された式（７０）における高次項が考慮される必要がある。議論のために、剰余οにおける高次キュムラントは低次キュムラントから独立したものではないということは、認めるに足るものである。
全ての波形が同一の（互いに同調かつ自己同調し、すなわちスペクトル異方性のない）スペクトルコンテントを有する、指向性および等方性拡散エンコーディングによる信号減衰ΔｌｎＥの差を考える。この場合、減衰差は、
Figure 0007218288000120
によって求められ、ここで、剰余における低次キュムラントへの依存性が留意される。すなわち、剰余は、〈Ｄ〉、Ｖ1、および〈ＶA〉に関する情報を密接に関連させる。次に、互いに同調され分光的等方性であり、かつ第１のエンコーディングペアと比べて異なるタイミングパラメータ（τ）を有する、すなわち異なる平均スペクトルコンテントを有する、指向性および等方性エンコーディングスキームの類似ペアを考える。ここで、第１および第２のエンコーディングペアに関する減衰差の差、
Figure 0007218288000121
を考えることができる。剰余項に起因して、差は、〈Ｄ〉およびＶ1における時間依存性（時間依存拡散）によって変調される。したがって差は、異方性コンパートメントにおける時間／周波数依存効果を分離するものではない。この問題を回避し、異方性構造の時間依存性（時間依存拡散）を完全に分離するために、スペクトル異方性を変化させることが提案される。
異なる程度のスペクトル異方性を有する２つの等方性エンコーディングスキーム、たとえば分光的異方性および等方性エンコーディングスキームを考える。これら２つのエンコーディングスキームに関する減衰差は、
Figure 0007218288000122
によって求められ、ここで
Figure 0007218288000123
は、スペクトル異方性に起因する残差分散を示すために用いられる。この場合、平均スペクトルコンテントは変更されず、異方性コンパートメントのみに影響を及ぼすと予想されるスペクトル異方性のみが変更されるので、
Figure 0007218288000124
のみが影響を受け、〈Ｄ〉またはＶ1は影響を受けない。その結果、剰余項も、異方性コンパートメントにおける時間／周波数依存拡散効果によってのみ影響を受ける。
例
図１は、２つの拡散エンコーディングシーケンスに関するディフェージング波形およびパワースペクトラム、および対応する減衰曲線のシミュレーションを示す。挿入図Ａ１は、２つの直交する方向である点線および破線に沿った（すなわち対応するＢテンソルの固有ベクトルに沿った）第１の拡散エンコーディングシーケンスに関するスケーリングされたディフェージングベクトル波形を示す。挿入図Ｂ１は、波形に関する対応するパワースペクトラムを示す。挿入図Ａ２は、２つの直交する方向である点線および破線に沿った（すなわち。対応するＢテンソルの固有ベクトルに沿った）第２の拡散エンコーディングシーケンスに関するスケーリングされたディフェージングベクトル波形を示す。挿入図Ｂ１は、波形に関する対応するパワースペクトラムを示す。
第１の拡散エンコーディングシーケンスは、面内分光的等方性（ＳＡ＝０）である。これは、挿入図Ｂ１における重なり合うパワースペクトラムから分かる。第２の拡散エンコーディングシーケンスは、面内分光的異方性（ＳＡ＝０．６８）である。挿入図Ｂ２において、パワースペクトラムは重なり合わない。ただし第１および第２の拡散シーケンスは、分光的に同調している。両方のスキームは、同一のＢテンソルならびにｍ(2)テンソルのトレースを生じる。
挿入図Ｃは、τ＝９５ｍｓにおける正規化信号減衰を示す。実線および破線はそれぞれ、分光的等方性波形（Ａ１、Ｂ１）および分光的異方性波形（Ａ２、Ｂ２）を用いた計算から得られる。マーカは、分光的等方性波形に関するモンテカルロシミュレーションの結果（丸印、Ａ１、Ｂ１）および分光的異方性波形に関するモンテカルロシミュレーションの結果（バツ印、Ａ２、Ｂ２）を示す。信号減衰差は、時間依存異方性拡散のみに関して観測され、時間依存等方性拡散に関して観測されない。
計算における粉体平均化は、１０００の均一に分布した回転に適用された。軸対称楕円細孔における制限拡散に関する拡散スペクトラムは、２つの軸に沿って１μｍ、１つの軸に沿って５μｍの半径を有する球形構造に関してＤ（ω）を用いて近似された。固有の拡散率Ｄ0＝１０-9ｍ2／ｓが用いられた。エンコーディング時間ｔe＝２τ＋ｔmは、０．０５～６２５ｍｓの範囲内で対数的に間隔をあけられた。混合時間は、ｔm＝０．２ｔeに調整された。加えて、モンテカルロシミュレーションは、５・１０5のウォーカおよび時間ステップを用いて同じ幾何学的形状および固有拡散率によって実行された。シミュレーションにおける粉体平均化は、ｂ＝０～１８・１０9ｓ／ｍ2の範囲内で１０の直線的に間隔をおいた拡散強調ステップを有する２４回転に適用された。
結果が示すように、一般の多次元拡散エンコーディング（ＭＤＥ）実験における時間依存拡散の影響は、周波数領域において分析され得る。（ディフェージングモーメントテンソルと称され得る）ｍまたはＭテンソルは、スペクトルコンテントに関してエンコーディング波形の追加の特徴付けを提供する。Ｂテンソルは、長期またはガウシアン拡散率に対する感度を提供するが、ディフェージングモーメントテンソルは、時間依存拡散に対する感度を提供する。Ｂテンソルの異方性を変化させることによって長期拡散異方性に特化したコントラストがもたらされるが、スペクトル異方性を変化させることによって、時間依存拡散異方性に特化したコントラストが生じる。平面拡散エンコーディングを用いるこの例において、スペクトル異方性は、平均スペクトルコンテントから独立して変化してよい。その結果生じる信号減衰差は、異方性領域における時間依存拡散を分離する、特異性の高い情報を伝達する。
実施形態の説明
図２は、サンプルに拡散強調磁気共鳴測定を実行する方法の一般的なフローチャートを示す。サンプルはたとえば、脳組織または任意の器官細胞の（懸濁液の）生体標本といった、水を含む生体サンプルであってよい。ただし方法は、汎用性を有し、たとえば岩石など他の種類のサンプルの分析にも用いられ得る。より一般的には、サンプルは、核磁気共鳴技術によって測定され得る特性を有する核スピン系を含む。
方法の理解を容易にするために、以下において、単一ボクセル、すなわち（ＭＲＩ法の場合）単一空間チャネルまたは（ＮＭＲ法の場合）単一周波数チャネルからの減衰エコー信号が参照される。当該技術において周知であるように、この分解能は、エンコーディングシーケンス中にサンプルに更なる磁気勾配（たとえば、ＭＲＩ法の場合、イメージング勾配）を適用することによって実現され得る。ボクセルに対応するサンプルの部分体積からエコー信号成分を識別／隔離するために、サンプルからの測定信号は、当該技術において周知であるように高速フーリエ変換を受けてよく、それによって、各エコー信号のスペクトル成分は、サンプルからサンプルの複数の空間または周波数領域に変換される。
当該技術において周知であるように、ＮＭＲスペクトロメータまたはＭＲＩデバイスの空間分解能は、とりわけ磁界強度、サンプルに適用される勾配パルスシーケンスの大きさ、およびスルーレートによって制限される。したがって、ボクセルに関するエコー信号は一般に、ボクセルに対応するサンプルの部分体積内の複数の微視的コンパートメントからの寄与を含む。
方法は、最先端のＮＭＲスペクトロメータまたはＭＲＩデバイスを用いて実行され得る。当該技術において周知であるように、そのようなデバイスは、デバイスの動作、とりわけ磁気勾配パルスシーケンスの生成、信号の取得ならびに取得した信号（すなわち測定信号）を表すデータを形成するための取得した信号のサンプリングおよびデジタル化を制御するためのコントローラを含んでよい。コントローラは、ＭＲＩデバイスの１または複数のプロセッサに実装されてよく、緩和エンコーディングシーケンスおよび磁気勾配パルスシーケンスの生成は、コンピュータ可読媒体（たとえば非一時的コンピュータ可読記憶媒体）に格納されデバイスの１または複数のプロセッサによって実行され得るソフトウェア命令を用いて実装され得る。ソフトウェア命令はたとえば、デバイスの１または複数のプロセッサがアクセス権を有する、デバイスのメモリのプログラム／制御部に格納され得る。ただし、ソフトウェア命令を用いる代わりに、コントローラ方法は、たとえば数例を挙げると１または複数の集積回路、１または複数の特定用途向け集積回路（ＡＳＩＣ）またはフィールドプログラマブルゲートアレイ（ＦＰＧＡ）など、デバイス／コンピュータの専用回路の形式で実装されてもよい。
取得した信号を表すデータは、デバイスの、またはデバイスに接続され得るコンピュータなどのデータメモリに格納され得る。
データ処理は、処理デバイスによって実行され得る。動作は、非一時的コンピュータ可読媒体において格納または具体化され処理デバイスによって実行され得るソフトウェア命令のセットに実装され得る。たとえばソフトウェア命令は、ＮＭＲスペクトロメータ／ＭＲＩデバイスのメモリのプログラム／制御部に格納され、スペクトロメータ／デバイスの１または複数のプロセッサユニットによって実行され得る。ただし、ＮＭＲスペクトロメータまたはＭＲＩデバイスから独立したデバイス、たとえばコンピュータ上で計算を実行することも同様に可能である。デバイスおよびコンピュータはたとえば、ＬＡＮ／ＷＬＡＮなどの通信ネットワークを介して、または他の何らかの直列または並列通信インタフェースを介して通信するように構成され得る。さらに留意すべき点として、ソフトウェア命令を用いる代わりに、データ処理は、たとえば数例を挙げると１または複数の集積回路、１または複数のＡＳＩＣまたはＦＰＧＡなど、デバイス／コンピュータの専用回路の形式で処理デバイスに実装されてもよい。
図２を参照すると、方法は、サンプルに複数の拡散強調磁気共鳴測定を実行すること（ステップ７０２－１～７０２ｎ）を備える。ＮＭＲ／ＭＲＩデバイスは、拡散エンコーディングシーケンスを生成し、その結果生じる信号減衰を取得してよい。各測定は、エンコーディングブロック、および後続の検出ブロックを含んでよい。エンコーディングブロックは、ＲＦシーケンスおよび拡散エンコーディング磁気勾配シーケンスを含んでよい。当該技術において周知であるように、エコー信号を取得するために、各拡散エンコーディング磁気勾配パルスシーケンスは、１または複数のイメージング磁気勾配および任意選択的に磁気勾配補正勾配によって補足され得る。
信号検出のブロックに先行して緩和および拡散エンコーディングのブロックを備えるパルスシーケンスの例が図３ａに示され、特定の実装が図３ｂに示される。したがって、図３ａは、緩和速度および拡散エンコーディング磁気勾配シーケンスの値に従ってエコー信号を変調するエンコーディングブロック、およびエコー信号が（たとえばスペクトラムまたは画像として）読み出される検出ブロックを示す。図３ｂは、９０°および１８０°のＲＦパルス（細い縦線および太い縦線）を有するパルスシーケンス、３つの直交する方向（実線、破線、および点線）における変調勾配、および検出信号（太線）を示す。信号は、縦方向の復元、横方向の緩和、および拡散によって変調される。
ゼロに等しい複合横方向磁化ｍｘｙを有する初期状態から開始し、第１の９０°ＲＦパルスは、縦方向磁化ｍｚを横断面にフリップする。持続期間τ1を有する時間遅延の間、縦方向磁化は、縦方向緩和速度Ｒ１で熱平衡値ｍ0に向かって復元する。第２の９０°パルスは、復元磁化を横断面にフリップし、ここで、検出されるまでの期間τ2の間、横方向緩和速度Ｒ１でゼロに向かって衰退する。期間τ2中、時間依存磁場勾配が適用される。
一般に、スピンエコーエンコーディングおよび励起エコーエンコーディングの両方が用いられ得る。いずれの場合も、ＲＦ信号シーケンスは、縦方向緩和のみ、横方向緩和のみ、または縦方向緩和および横方向緩和の両方に起因する減衰をエンコードしてよい。シーケンスの一例は、単一９０°パルスおよび単一１８０°パルスを含んでよい。１８０°パルスに関する勾配パルスシーケンスのタイミングは変化してよい。たとえば、勾配パルスシーケンスは、１８０°パルスに先行または後続して実行され得る。そのようないくつかのシーケンスは、取得／検出まで反復され得る。励起エコーシーケンスの例は、第１の９０°パルス、第２の９０°パルス、および第３の９０°パルスを含んでよい。勾配パルスシーケンスは、第１および第２の９０°パルスの間に、および／または第３の９０°パルスに続いて（すなわち、検出ブロックより前に）実行され得る。ただし、これらのシーケンス例は典型例として提供されたものにすぎず、他のシーケンスも可能である。好適には、同じＲＦ信号シーケンスが、複数の測定の全て、すなわち第１、第２、および第３のセットの全ての測定において用いられる。
方法は、少なくとも２つの非ゼロ固有値を有する第１の拡散強調テンソル表現Ｂ1を有する第１の拡散エンコーディングシーケンスを用いて第１の測定７０２－１を実行することを備えてよい。方法は、少なくとも２つの非ゼロ固有値を有する第２の拡散強調テンソル表現Ｂ2を有する第２の拡散エンコーディングシーケンスを用いて第２の測定７０２－２を実行することを備えてよい。第１のテンソル表現Ｂ1および第２のテンソル表現Ｂ2は、同じ数の非ゼロ固有値を有してよい。また、第１のテンソル表現Ｂ1の固有値は、第２のテンソル表現Ｂ2の固有値と一致してよい（たとえば、１０％未満異なり、あるいは等しく、たとえば１％未満異なるなど少なくとも実質的に等しい）。また、第１および第２の拡散エンコーディングシーケンスは、上で定義されたように、一致する平均スペクトルコンテントを示し、異なる程度のスペクトル異方性を示すように構成され得る。
２つのそのような拡散エンコーディングシーケンスの例が図１に示される。図１において、シーケンスは、平面エンコーディングシーケンスであり、すなわち、厳密に２つの非ゼロ固有値を有するテンソル表現を有する。ただし、３Ｄ拡散エンコーディングシーケンス、すなわち厳密に３つの非ゼロ固有値を有するテンソル表現を有する拡散エンコーディングシーケンスを生成することも可能である。拡散エンコーディングシーケンスは、上述した数値的アプローチを用いて設計され得る。追加または代替として、拡散エンコーディングシーケンスは、上述したテストケースを利用して設計され得る。テストケースが拡散エンコーディングシーケンスの初期設計に合わない場合、拡散エンコーディングシーケンスのうち１または複数の１または複数のパラメータが変更され得る反復アプローチが適用され得る。その後新たな試験が、修正された拡散エンコーディングシーケンスを用いて実行され得る。反復は、テストケースが合うまで行われてよい。試験を目的として、第１および第２の拡散エンコーディングシーケンスのｂ値は、１＊１０9～４＊１０9ｓ／ｍ2の範囲内の任意の単一の値であってよい。上述したようなテストケースは、直径５μｍの球状コンパートメントおよび円筒状コンパートメントを参照するが、たとえば５～１５μｍの範囲内のいずれかなど他の直径を用いることも可能である。この範囲内のコンパートメント径は、２つの拡散エンコーディングシーケンスがスペクトルの一致および異なる程度のスペクトル異方性という要件を満たすかを試験することも可能にする。
各測定の検出ブロックは、エンコーディングブロックに後続してエコー信号を検出することを含んでよい。複数の測定の結果生じる信号は、データとして取得および記録され得る。エコー信号は、データを形成するためにサンプリングおよびデジタル化され得る。データは、更なるデータ処理のために格納され得る。データは例えば、デバイスの、またはデバイスに接続され得るコンピュータなどのデータメモリに格納され得る。
測定は粉体平均化されてよく、第１および第２の測定の各々は、異なる測定方向、好適には複数の方向の数だけ反復され得る。第１の平均信号減衰は、各測定方向に関して、第１の拡散エンコーディングシーケンスの適用の結果生じる信号減衰の平均として計算され得る。対応して、第２の平均信号減衰は、各測定方向に関して、第２の拡散エンコーディングシーケンスの適用の結果生じる信号減衰の平均として計算され得る。
方法のステップ７０４において、処理デバイスは、磁気共鳴測定７０２－１および７０２－２の結果生じる、測定された第１および第２の（任意選択的に平均化された）信号減衰に基づいて出力を生成してよい。出力は例えば、第１および第２の信号減衰の差、または第１および第２の信号減衰の比に基づいてよい。処理デバイスは、サンプルの関心領域内の複数のボクセルの各々に関して測定された信号減衰に基づいて対応する出力を生成してよい。処理デバイスは、第１の測定において取得された信号減衰が第２の測定において取得された信号減衰と閾値より大きく異なるボクセルの指標を含む出力をデジタル画像の形式で生成してよい。たとえば、第１の測定が第２の測定において取得された信号減衰と異なるボクセルは、輝度の増加、色の逸脱、境界ボックス、および／または他の何らかのグラフィック要素を用いて強調され得る。
上記において、第１および第２の拡散エンコーディングシーケンスをそれぞれ用いる第１および第２の測定が参照されたが、方法は、たとえば第３の拡散エンコーディングシーケンスを用いる第３の測定など、更なる測定を含んでよいことに留意すべきである。第３の拡散エンコーディングシーケンスは、少なくとも２つの非ゼロ固有値を有する第２の拡散強調テンソル表現Ｂ3を有してよい。テンソル表現Ｂ3は、テンソル表現Ｂ1およびＢ2と同じ数の非ゼロ固有値を有してよい。また、第３のテンソル表現Ｂ３の固有値は、テンソル表現Ｂ1およびＢ2の固有値と一致してよい。さらに、第３の拡散エンコーディングシーケンスは、第１の拡散エンコーディングシーケンスの平均スペクトルコンテントおよび第２の拡散エンコーディングシーケンスの平均スペクトルコンテントと一致する平均スペクトルコンテントを示し、第１の拡散エンコーディングシーケンスのスペクトル異方性の程度および第２の拡散エンコーディングシーケンスのスペクトル異方性の程度と異なる程度のスペクトル異方性を示すように構成され得る。
図２を参照すると、方法のステップ７０２－１は、異なる拡散エンコーディング強度（ｂ値）および一致する平均スペクトルコンテントおよび一致する程度のスペクトル異方性を有する拡散エンコーディングシーケンスを用いて実行される第１の測定セットを実行することを含んでもよい。対応して、ステップ７０２－２は、異なる拡散エンコーディング強度および一致する平均スペクトルコンテントおよび一致する程度のスペクトル異方性を有する拡散エンコーディングシーケンスを用いて実行される第２の測定セットを実行することを含んでよい。第１および第２のセットの各々の拡散エンコーディングシーケンスは、それぞれのテンソル表現を有してよく、第１および第２のセットのそれぞれの拡散エンコーディングシーケンスに関するテンソル表現の非ゼロ固有値の数は等しい。
検出ブロックは、信号減衰の第１のセットを取得すること、および取得した信号減衰の第１のセットを表す第１のデータセットを記録することを含んでよい。対応して検出ブロックは、信号減衰の第２のセットを取得すること、および取得した信号減衰の第２のセットを表す第２のデータセットを記録することを含んでよい。任意選択的に、各セットの各測定は、上述したような対応する方法で粉体平均化され得る。
方法のステップ７０４において、処理デバイスは結果的に、第１の信号減衰曲線Ｅ（ｂ）を推定するために、（任意選択的に平均化された）信号減衰の第１のセットを表す第１のデータセットに関数をフィッティングしてよい。処理デバイスはさらに、第２の信号減衰曲線Ｅ（ｂ）を推定するために、（任意選択的に平均化された）信号減衰の第２のセットを表す第２のデータセットに関数をフィッティングしてよい。任意のフィッティング関数が、拡散スペクトラムＤ（ω）の成分または固有モードλi（ω）またはλk (n)（０）、すなわちゼロ周波数におけるλk（ω）のｎ番目の導関数（式４０）または有効制限テンソルρ4（式５２～５６）に関するパラメータを含む。処理デバイスはその結果、第１のフィッティングおよび第２のフィッティングの少なくとも１つのパラメータに基づいて出力を生成してよい。上記と同様、処理デバイスは、サンプルの関心領域内の複数のボクセルの各々に関して測定された信号減衰に基づいて対応する出力を生成してよい。
上記において、本発明概念は主に、限られた数の例を参照して説明された。しかし、当業者によって容易に理解されるように、開示した例ではない他の例が、添付の特許請求の範囲によって定義される本発明概念の範囲内で同等に可能である。
参照文献リスト
上記開示において、丸括弧または鍵括弧内の１または複数の数字は、以下の参照文献リストにおける、対応した番号の参照文献を参照するものである。
1. Westin C-F, Szczepankiewicz F, Pasternak O, 他: Measurement Tensors in Diffusion MRI: Generalizing the Concept of Diffusion Encoding. Med Image Comput Comput Assist Interv 2014; 17:209-216.
2. Lasic S, Szczepankiewicz F, Eriksson S, Nilsson M, Topgaard D: Microanisotropy imaging: quantification of microscopic diffusion anisotropy and orientational order parameter by diffusion MRI with magic-angle spinning of the q-vector. Front Phys 2014; 2:1-14.
3. Eriksson S, Lasic S, Topgaard D: Isotropic diffusion weighting in PGSE NMR by magic-angle spinning of the q-vector. J Magn Reson 2013; 226:13-8.
4. Eriksson S, Lasic S, Nilsson M, Westin C-F, Topgaard D: NMR diffusion-encoding with axial symmetry and variable anisotropy: Distinguishing between prolate and oblate microscopic diffusion tensors with unknown orientation distribution. J Chem Phys 2015; 142:104201.
5. Szczepankiewicz F, Lasic S, van Westen D,他: Quantification of microscopic diffusion anisotropy disentangles effects of orientation dispersion from microstructure: applications in healthy volunteers and in brain tumors. Neuroimage 2015; 104:241-52.
6. Stepisnik J: Validity limits of Gaussian approximation in cumulant expansion for diffusion attenuation of spin echo. Phys B 1999; 270:110-117.
7. Stepisnik J: Time-dependent self-diffusion by NMR spin-echo. Phys B 1993; 183:343-350.
以下に、本願出願の当初の特許請求の範囲に記載された発明を付記する。
［１］
サンプルに拡散強調磁気共鳴測定を実行する方法であって、
前記サンプルに拡散強調磁気共鳴測定を実行することを備え、
前記測定は、
少なくとも２つの非ゼロ固有値を有する第１の拡散強調テンソル表現Ｂ 1 を有する第１の拡散エンコーディングシーケンスを用いる第１の測定と、
少なくとも２つの非ゼロ固有値を有する第２の拡散強調テンソル表現Ｂ 2 を有する第２の拡散エンコーディングシーケンスを用いる第２の測定と
を含み、
前記第１の拡張強調テンソル表現Ｂ 1 および前記第２の拡張強調テンソル表現Ｂ 2 は、同じ数の非ゼロ固有値を有し、前記第１の拡張強調テンソル表現Ｂ 1 の前記固有値は、前記第２の拡張強調テンソル表現Ｂ 2 の前記固有値と一致し、
前記第１および前記第２の拡散エンコーディングシーケンスは、一致する平均スペクトルコンテントを示し、かつ、異なる程度のスペクトル異方性を示すように構成される、方法。
［２］
前記第１および前記第２の拡散エンコーディングシーケンスは、
前記第１の拡散エンコーディングシーケンスの第１の正規化ディフェージングベクトル表現Ｆ ~ 1 と一致する正規化ディフェージングベクトル表現Ｆ ~ 3 を有し、非ゼロ拡散エンコーディング強度を有する第３の拡散エンコーディングシーケンスが、直径５μｍの球状コンパートメントの集合から成る第１の試験サンプルに適用され、
前記第２の拡散エンコーディングシーケンスの第２の正規化ディフェージングベクトル表現Ｆ ~ 2 と一致する正規化ディフェージングベクトル表現Ｆ ~ 4 を有し、前記非ゼロ拡散エンコーディング強度を有する第４の拡散エンコーディングシーケンスが、前記第１の試験サンプルに適用され、
前記第３の拡散エンコーディングシーケンスの結果生じる信号減衰が、前記第４の拡散エンコーディングシーケンスの結果生じる信号減衰と一致する、
ように構成され、
前記第３の拡散エンコーディングシーケンスが、均一な配向分散を有する直径５μｍの円筒状コンパートメントの集合から成る第２の試験サンプルに適用され、
前記第４の拡散エンコーディングシーケンスが、前記第２の試験サンプルに適用され、
前記第３の拡散エンコーディングシーケンスの結果生じる信号減衰が、前記第４の拡散エンコーディングシーケンスの結果生じる信号減衰と異なる、
ように構成される、［１］の方法。
［３］
前記第１の測定は、第１の信号減衰を取得することを含み、前記第２の測定は、第２の信号減衰を取得することを含む、［１］または［２］の方法。
［４］
前記測定は、
前記第１の測定と、参照測定フレームに対して異なる回転による前記サンプルに適用される前記第１の拡散エンコーディングシーケンスを用いて実行される複数の追加の測定と、を含む第１の測定セットと、
前記第２の測定と、前記参照測定フレームに対して異なる回転による前記サンプルに適用される前記第２の拡散エンコーディングシーケンスを用いて実行される複数の追加の測定と、を含む第２の測定セットと、
を含む、［１］～［３］のいずれかの方法。
［５］
前記第１の測定セットの各測定は、それぞれの信号減衰を取得することを含み、前記方法は、前記それぞれの信号減衰に基づいて第１の平均信号減衰を決定することを備え、
前記第２の測定セットの各測定は、それぞれの信号減衰を取得することを含み、前記方法は、前記それぞれの信号減衰に基づいて第２の平均信号減衰を決定することを備える、［４］の方法。
［６］
前記第１および第２の信号減衰の差または前記第１および第２の平均信号減衰の差を示す出力を生成することをさらに備える、［１］または［２］の方法。
［７］
前記第１の拡張強調テンソル表現は、３つの一致する非ゼロ固有値を有する、［１］～［６］のいずれかの方法。
［８］
前記第２の拡張強調テンソル表現は、３つの一致する非ゼロ固有値を有する、［１］～［７］のいずれかの方法。
［９］
前記測定を実行することは、前記サンプルの関心領域内の複数のボクセルの各々から、前記測定の各々の結果生じるそれぞれの信号減衰を測定することを含む、［１］～［８］のいずれかの方法。
［１０］
前記第１の測定において取得された信号減衰が前記第２の測定において取得された信号減衰と異なるボクセルの指標を含む出力を生成することをさらに備える、［９］の方法。
［１１］
サンプルに拡散強調磁気共鳴測定を実行する方法であって、
前記サンプルに、拡散強調磁気共鳴測定を実行することを備え、
前記測定は、
異なる拡散エンコーディング強度および一致する平均スペクトルコンテントおよび一致する程度のスペクトル異方性を有する拡散エンコーディングシーケンスを用いて実行される第１の測定セットと、
異なる拡散エンコーディング強度および一致する平均スペクトルコンテントおよび一致する程度のスペクトル異方性を有する拡散エンコーディングシーケンスを用いて実行される第２の測定セットと、
を含み、
前記第１の測定セットの前記拡散エンコーディングシーケンスの前記スペクトル異方性の程度は、前記第２の測定セットの前記拡散エンコーディングシーケンスの前記スペクトル異方性の程度と異なる、方法。
［１２］
第１の信号減衰曲線を推定するために前記第１の測定セットを表す第１のデータセットに第１の関数をフィッティングすることと、
第２の信号減衰曲線を推定するために前記第２の測定セットを表す第２のデータセットに第２の関数をフィッティングすることと、
をさらに備える、［１１］の方法。
［１３］
前記第１の関数の少なくとも１つのパラメータおよび前記第２の関数の少なくとも１つのパラメータに基づいて出力を生成することをさらに備える、［１２］の方法。
［１４］
前記第１の測定セットおよび前記第２の測定セットの各々の前記拡散エンコーディングシーケンスはそれぞれのテンソル表現を有し、前記テンソル表現の非ゼロ固有値の数は等しい、［１１］～［１３］のいずれか１項の方法。
［１５］
前記第１の測定セットの前記拡散エンコーディングシーケンスの各々の前記平均スペクトルコンテントは、前記第２の測定セットの前記拡散エンコーディングシーケンスの各々の前記平均スペクトルコンテントと一致する、［１１］～［１４］のいずれか１項の方法。

Claims (16)

  1. サンプルに拡散強調磁気共鳴測定を実行する方法であって、
    前記サンプルに拡散強調磁気共鳴測定を実行することを備え、
    前記測定は、
    少なくとも２つの非ゼロ固有値を有する第１の拡散強調テンソル表現Ｂ1を有する第１の拡散エンコーディングシーケンスを用いる第１の測定と、
    少なくとも２つの非ゼロ固有値を有する第２の拡散強調テンソル表現Ｂ2を有する第２の拡散エンコーディングシーケンスを用いる第２の測定と
    を含み、
    前記第１の拡散強調テンソル表現Ｂ1および前記第２の拡散強調テンソル表現Ｂ2は、同じ数の非ゼロ固有値を有し、前記第１の拡散強調テンソル表現Ｂ1の前記固有値は、前記第２の拡散強調テンソル表現Ｂ2の前記固有値と一致し、
    前記第１および前記第２の拡散エンコーディングシーケンスは、一致する平均スペクトルコンテントを示し、かつ、異なる程度のスペクトル異方性を示すように構成される、方法。
  2. 前記第１および前記第２の拡散エンコーディングシーケンスは、
    前記第１の拡散エンコーディングシーケンスの第１の正規化ディフェージングベクトル表現Ｆ~ 1と一致する正規化ディフェージングベクトル表現Ｆ~ 3を有し、非ゼロ拡散エンコーディング強度を有する第３の拡散エンコーディングシーケンスが、直径５μｍの球状コンパートメントの集合から成る第１の試験サンプルに適用され、
    前記第２の拡散エンコーディングシーケンスの第２の正規化ディフェージングベクトル表現Ｆ~ 2と一致する正規化ディフェージングベクトル表現Ｆ~ 4を有し、前記非ゼロ拡散エンコーディング強度を有する第４の拡散エンコーディングシーケンスが、前記第１の試験サンプルに適用され、
    前記第３の拡散エンコーディングシーケンスの結果生じる信号減衰が、前記第４の拡散エンコーディングシーケンスの結果生じる信号減衰と一致する、
    ように構成され、
    前記第３の拡散エンコーディングシーケンスが、均一な配向分散を有する直径５μｍの円筒状コンパートメントの集合から成る第２の試験サンプルに適用され、
    前記第４の拡散エンコーディングシーケンスが、前記第２の試験サンプルに適用され、
    前記第３の拡散エンコーディングシーケンスの結果生じる信号減衰が、前記第４の拡散エンコーディングシーケンスの結果生じる信号減衰と異なる、
    ように構成される、請求項１の方法。
  3. 前記第１の測定は、第１の信号減衰を取得することを含み、前記第２の測定は、第２の信号減衰を取得することを含む、請求項１または２の方法。
  4. 前記測定は、
    前記第１の測定と、参照測定フレームに対して異なる回転による前記サンプルに適用される前記第１の拡散エンコーディングシーケンスを用いて実行される複数の追加の測定と、を含む第１の測定セットと、
    前記第２の測定と、前記参照測定フレームに対して異なる回転による前記サンプルに適用される前記第２の拡散エンコーディングシーケンスを用いて実行される複数の追加の測定と、を含む第２の測定セットと、
    を含む、請求項１～３のいずれかの方法。
  5. 前記第１の測定セットの各測定は、それぞれの信号減衰を取得することを含み、前記方法は、前記それぞれの信号減衰に基づいて第１の平均信号減衰を決定することを備え、
    前記第２の測定セットの各測定は、それぞれの信号減衰を取得することを含み、前記方法は、前記それぞれの信号減衰に基づいて第２の平均信号減衰を決定することを備える、請求項４の方法。
  6. 前記第１の測定は、第１の信号減衰を取得することを含み、前記第２の測定は、第２の信号減衰を取得することを含み、
    前記方法は、前記第１および第２の信号減衰の差を示す出力を生成することをさらに備える、請求項１または２の方法。
  7. 前記第１の拡散強調テンソル表現は、３つの一致する非ゼロ固有値を有する、請求項１～６のいずれかの方法。
  8. 前記第２の拡散強調テンソル表現は、３つの一致する非ゼロ固有値を有する、請求項１～７のいずれかの方法。
  9. 前記測定を実行することは、前記サンプルの関心領域内の複数のボクセルの各々から、前記測定の各々の結果生じるそれぞれの信号減衰を測定することを含む、請求項１～８のいずれかの方法。
  10. 前記第１の測定において取得された信号減衰が前記第２の測定において取得された信号減衰と異なるボクセルの指標を含む出力を生成することをさらに備える、請求項９の方法。
  11. サンプルに拡散強調磁気共鳴測定を実行する方法であって、
    前記サンプルに、拡散強調磁気共鳴測定を実行することを備え、
    前記測定は、
    異なる拡散エンコーディング強度および一致する平均スペクトルコンテントおよび一致する程度のスペクトル異方性を有する拡散エンコーディングシーケンスを用いて実行される第１の測定セットと、
    異なる拡散エンコーディング強度および一致する平均スペクトルコンテントおよび一致する程度のスペクトル異方性を有する拡散エンコーディングシーケンスを用いて実行される第２の測定セットと、
    を含み、
    前記第１の測定セットの前記拡散エンコーディングシーケンスの前記スペクトル異方性の程度は、前記第２の測定セットの前記拡散エンコーディングシーケンスの前記スペクトル異方性の程度と異なる、方法。
  12. 第１の信号減衰曲線を推定するために前記第１の測定セットを表す第１のデータセットに第１の関数をフィッティングすることと、
    第２の信号減衰曲線を推定するために前記第２の測定セットを表す第２のデータセットに第２の関数をフィッティングすることと、
    をさらに備える、請求項１１の方法。
  13. 前記第１の関数の少なくとも１つのパラメータおよび前記第２の関数の少なくとも１つのパラメータに基づいて出力を生成することをさらに備える、請求項１２の方法。
  14. 前記第１の測定セットおよび前記第２の測定セットの各々の前記拡散エンコーディングシーケンスはそれぞれのテンソル表現を有し、前記テンソル表現の非ゼロ固有値の数は等しい、請求項１１～１３のいずれか１項の方法。
  15. 前記第１の測定セットの前記拡散エンコーディングシーケンスの各々の前記平均スペクトルコンテントは、前記第２の測定セットの前記拡散エンコーディングシーケンスの各々の前記平均スペクトルコンテントと一致する、請求項１１～１４のいずれか１項の方法。
  16. 前記測定は、
    前記第１の測定と、参照測定フレームに対して異なる回転による前記サンプルに適用される前記第１の拡散エンコーディングシーケンスを用いて実行される複数の追加の測定と、を含む第１の測定セットと、
    前記第２の測定と、前記参照測定フレームに対して異なる回転による前記サンプルに適用される前記第２の拡散エンコーディングシーケンスを用いて実行される複数の追加の測定と、を含む第２の測定セットと、
    を含み、
    前記第１の測定セットの各測定は、それぞれの信号減衰を取得することを含み、前記方法は、前記それぞれの信号減衰に基づいて第１の平均信号減衰を決定することを備え、
    前記第２の測定セットの各測定は、それぞれの信号減衰を取得することを含み、前記方法は、前記それぞれの信号減衰に基づいて第２の平均信号減衰を決定することを備え、
    前記方法は、前記第１および第２の平均信号減衰の差を示す出力を生成することをさらに備える、請求項１または２の方法。
JP2019522714A 2016-11-09 2017-11-09 サンプルに拡散強調磁気共鳴測定を実行する方法 Active JP7218288B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
SE1651469-7 2016-11-09
SE1651469 2016-11-09
PCT/SE2017/051126 WO2018088955A1 (en) 2016-11-09 2017-11-09 A method of performing diffusion weighted magnetic resonance measurements on a sample

Publications (2)

Publication Number Publication Date
JP2019534087A JP2019534087A (ja) 2019-11-28
JP7218288B2 true JP7218288B2 (ja) 2023-02-06

Family

ID=62111280

Family Applications (2)

Application Number Title Priority Date Filing Date
JP2019522646A Active JP7057355B2 (ja) 2016-11-09 2017-11-09 サンプルに拡散強調磁気共鳴測定を実行する方法
JP2019522714A Active JP7218288B2 (ja) 2016-11-09 2017-11-09 サンプルに拡散強調磁気共鳴測定を実行する方法

Family Applications Before (1)

Application Number Title Priority Date Filing Date
JP2019522646A Active JP7057355B2 (ja) 2016-11-09 2017-11-09 サンプルに拡散強調磁気共鳴測定を実行する方法

Country Status (8)

Country Link
US (3) US11061096B2 (ja)
EP (2) EP3538910B1 (ja)
JP (2) JP7057355B2 (ja)
KR (2) KR102574708B1 (ja)
CN (2) CN109923426B (ja)
AU (2) AU2017356785B2 (ja)
CA (2) CA3041785A1 (ja)
WO (2) WO2018088955A1 (ja)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3538910B1 (en) * 2016-11-09 2023-06-07 Random Walk Imaging AB A method of performing diffusion weighted magnetic resonance measurements on a sample
SE543292C2 (en) * 2019-04-26 2020-11-17 Cr Dev Ab A method of performing diffusion weighted magnetic resonance measurements
CN110680322B (zh) * 2019-10-17 2020-07-31 河海大学 一种描述磁共振成像信号非指数衰减的方法及其应用
WO2021201753A1 (en) * 2020-03-28 2021-10-07 Oezarslan Evren A magnetic resonance method, software product, and system for determining a diffusion propagator or related diffusion parameters for spin-labelled particles
CN114646913B (zh) * 2020-12-21 2023-06-02 北京大学第三医院（北京大学第三临床医学院） 生物组织微观结构的无创测量方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070238969A1 (en) 2006-01-24 2007-10-11 Schlumberger Technology Corporation Diffusion-based magnetic resonance methods for characterizing bone structure
US20150115957A1 (en) 2012-05-04 2015-04-30 Cr Development Ab Pulse sequence method for mri
US20150130458A1 (en) 2012-05-04 2015-05-14 Cr Development Ab Analysis for quantifying microscopic diffusion anisotropy
WO2015119569A1 (en) 2014-02-10 2015-08-13 Cr Development Ab Method for quantifying isotropic diffusion and/or anisotropic diffusion in a sample
US20150253410A1 (en) 2012-09-28 2015-09-10 Children's Medical Center Corporation Diffusion-weighted mri using multiple b-values and constant echo time

Family Cites Families (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3514337B2 (ja) * 1995-04-19 2004-03-31 株式会社日立メディコ ３次元データ位置合わせ方法
US5969524A (en) * 1997-04-14 1999-10-19 The United States Of America As Represented By The Department Of Health And Human Services Method to significantly reduce bias and variance of diffusion anisotrophy measurements
US6288540B1 (en) * 1999-05-21 2001-09-11 University Of Rochester Optimized orthogonal gradient technique for fast quantitative diffusion MRI on a clinical scanner
JP4837218B2 (ja) 2000-03-31 2011-12-14 ザ ジェネラル ホスピタル コーポレーション 組織の拡散イメージング
WO2002082376A2 (en) 2001-04-06 2002-10-17 Regents Of The University Of California Method for analyzing mri diffusion data
US6806705B2 (en) 2002-05-15 2004-10-19 Koninklijke Philips Electronics N.V. Diffusion tensor magnetic resonance imaging including local weighted interpolation
EP1611452A1 (en) * 2003-03-31 2006-01-04 Koninklijke Philips Electronics N.V. A method of magnetic resonance perfusion imaging
DE102005053269B3 (de) 2005-11-08 2007-04-12 Siemens Ag Verfahren und Vorrichtung zur Ermittlung eines diffusionsgewichteten Bildes
CN101143093B (zh) 2006-09-11 2010-09-29 西门子（中国）有限公司 磁共振扩散成像方法
US7355407B1 (en) * 2006-12-03 2008-04-08 Toshiba Medical Systems Corp. Methods and apparatus for single-shot magnetic resonance imaging with optimized isotropic diffusion weighting
US20090091322A1 (en) * 2007-04-25 2009-04-09 Stefan Posse Single-shot magnetic resonance spectroscopic imaging with partial parallel imaging
US7906965B2 (en) * 2007-11-05 2011-03-15 The United States Of America As Represented By The Department Of Health And Human Services Methods and apparatuses for estimating the elliptical cone of uncertainty
CN102077108B (zh) 2008-04-28 2015-02-25 康奈尔大学 分子mri中的磁敏度精确量化
JP5591493B2 (ja) * 2008-07-17 2014-09-17 株式会社東芝 磁気共鳴イメージング装置
JP5686729B2 (ja) 2009-04-22 2015-03-18 株式会社日立メディコ 磁気共鳴イメージング装置及び繊維状組織の走行方向表示方法
DE102009019895B4 (de) 2009-05-04 2011-05-12 Siemens Aktiengesellschaft Verfahren und Magnetresonanzanlage zur diffusionsgewichteten Aufnahme von MR-Signalen
SE0950363A1 (sv) 2009-05-22 2010-06-29 Cr Dev Ab Metod och system för magnetisk resonanstomografi, samt användning därav.
DE102010001577B4 (de) 2010-02-04 2012-03-08 Siemens Aktiengesellschaft Verfahren zur Reduktion von Verzerrungen in der Diffusionsbildgebung und Magnetresonanzanlage
US8742754B2 (en) * 2010-02-16 2014-06-03 Board Of Regents Of The University Of Texas System Method and system for diffusion tensor imaging
DE102011006851B4 (de) 2011-04-06 2012-12-06 Siemens Aktiengesellschaft Gruppierung von diffusionsgewichteten MR-Bildern nach den verwendeten Diffusionsgradienten zur Erstellung einer Diffusionsinformation
US9720124B2 (en) * 2011-08-10 2017-08-01 Schlumberger Technology Corporation Logging in gas shale and other unconventional reservoirs
US10228335B2 (en) * 2013-01-03 2019-03-12 Schlumberger Technology Corporation Method for nuclear magnetic resonance diffusion measurements
KR20150110685A (ko) 2013-01-25 2015-10-02 리전츠 오브 더 유니버스티 오브 미네소타 다중 채널 송신기를 위한 다중 대역 rf/mri 펄스 디자인
US10048345B2 (en) * 2013-03-07 2018-08-14 The Board Of Trustees Of The University Of Illinois Fractional order and entropy bio-markers for biological tissue in diffusion weighted magnetic resonance imaging
WO2015042416A1 (en) 2013-09-20 2015-03-26 Children's Medical Center Corporation Methods and apparatus for modeling diffusion-weighted mr data acquired at multiple non-zero b-values
US10817500B2 (en) * 2014-03-13 2020-10-27 Sybase, Inc. Reduce log contention by batching log record transfers to the log
EP2955536B1 (en) * 2014-06-12 2020-08-12 Commissariat A L'energie Atomique Et Aux Energies Alternatives MRI method to quantify iron amount in tissues using diffusion magnetic resonance imaging
US9233763B1 (en) 2014-08-19 2016-01-12 Gulfstream Aerospace Corporation Methods and systems for aircraft systems health trend monitoring
EP3081955A1 (en) 2015-04-13 2016-10-19 Commissariat A L'energie Atomique Et Aux Energies Alternatives Mri method for determining signature indices of an observed tissue from signal patterns obtained by motion-probing pulsed gradient mri
SE538834C2 (sv) * 2015-12-29 2016-12-20 Cr Dev Ab Method of extracting information about a sample by nuclear magnetic resonance measurements
DE102016202254B4 (de) 2016-02-15 2017-11-30 Siemens Healthcare Gmbh Modellfreies Ermitteln von Bildbereichen mit anomaler Diffusion anhand von diffusionsgewichteten Magnetresonanzbilddaten
WO2017190029A1 (en) * 2016-04-29 2017-11-02 Washington University System and method for multi-modality quantification of neuroinflammation in central nervous system diseases
EP3538910B1 (en) * 2016-11-09 2023-06-07 Random Walk Imaging AB A method of performing diffusion weighted magnetic resonance measurements on a sample

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070238969A1 (en) 2006-01-24 2007-10-11 Schlumberger Technology Corporation Diffusion-based magnetic resonance methods for characterizing bone structure
JP2009524830A (ja) 2006-01-24 2009-07-02 シュルンベルジェ ホールディングス リミテッド 骨構造をキャラクタライズするための拡散型磁気共鳴法
US20150115957A1 (en) 2012-05-04 2015-04-30 Cr Development Ab Pulse sequence method for mri
US20150130458A1 (en) 2012-05-04 2015-05-14 Cr Development Ab Analysis for quantifying microscopic diffusion anisotropy
JP2015518408A (ja) 2012-05-04 2015-07-02 セェーアール ディベロップメント アクチエボラグＣｒ Ｄｅｖｅｌｏｐｍｅｎｔ Ａｂ Ｍｒｉのパルス・シーケンス方法
JP2015518568A (ja) 2012-05-04 2015-07-02 セェーアール ディベロップメント アクチエボラグＣｒ Ｄｅｖｅｌｏｐｍｅｎｔ Ａｂ 微視的な拡散異方性を定量化するための解析
US20150253410A1 (en) 2012-09-28 2015-09-10 Children's Medical Center Corporation Diffusion-weighted mri using multiple b-values and constant echo time
WO2015119569A1 (en) 2014-02-10 2015-08-13 Cr Development Ab Method for quantifying isotropic diffusion and/or anisotropic diffusion in a sample
JP2017507698A (ja) 2014-02-10 2017-03-23 シーアール ディベロップメント アーベー サンプルにおける等方性拡散及び／又は異方性拡散を定量化するための方法

Also Published As

Publication number Publication date
KR20190077448A (ko) 2019-07-03
JP7057355B2 (ja) 2022-04-19
AU2017356784B2 (en) 2023-10-19
AU2017356785B2 (en) 2023-10-12
AU2017356785A1 (en) 2019-05-02
JP2019534087A (ja) 2019-11-28
KR20190073500A (ko) 2019-06-26
US20190265323A1 (en) 2019-08-29
CA3041784A1 (en) 2018-05-17
CN109923426A (zh) 2019-06-21
EP3538910B1 (en) 2023-06-07
JP2019534086A (ja) 2019-11-28
US10948560B2 (en) 2021-03-16
CN109923427B (zh) 2021-09-07
CN109923427A (zh) 2019-06-21
KR102620599B1 (ko) 2024-01-02
CN109923426B (zh) 2021-09-03
KR102574708B1 (ko) 2023-09-04
EP3538909A1 (en) 2019-09-18
CA3041785A1 (en) 2018-05-17
AU2017356784A1 (en) 2019-05-02
US11061096B2 (en) 2021-07-13
WO2018088955A1 (en) 2018-05-17
EP3538910A4 (en) 2020-07-15
USRE49978E1 (en) 2024-05-21
EP3538909A4 (en) 2020-07-15
EP3538910A1 (en) 2019-09-18
US20190361083A1 (en) 2019-11-28
WO2018088954A1 (en) 2018-05-17

Similar Documents

Publication Publication Date Title
JP7218288B2 (ja) サンプルに拡散強調磁気共鳴測定を実行する方法
Westin et al. Q-space trajectory imaging for multidimensional diffusion MRI of the human brain
AU2015214638B2 (en) Method for quantifying isotropic diffusion and/or anisotropic diffusion in a sample
Lundell et al. Diffusion encoding with general gradient waveforms
Cai et al. Rapid detection of the presence of diffusion exchange
JP7227438B2 (ja) 核磁気共鳴測定により試料に関する情報を抽出する方法
US10871539B2 (en) Determination of a joint probability distribution of radius and length of anisotropic pores from double pulsed field gradient MRI data
Zong et al. Determining mean fractional anisotropy using DDCOSY: preliminary results in biological tissues
Ianuş et al. Measuring microscopic anisotropy with diffusion magnetic resonance: from material science to biomedical imaging
Topgaard Probing biological tissue microstructure with magnetic resonance diffusion techniques
JP2023519217A (ja) スピン標識化粒子の拡散伝搬関数又は関連する拡散パラメータを決定するための磁気共鳴法、ソフトウェア製品、及びシステム
Zhang et al. Diffusion configuration of water molecules in diffusion weighted imaging
Mesri Investigating Microstructure in the Brain Using Diffusion Weighted Magnetic Resonance Imaging

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20201106

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20210728

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210803

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20211101

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220222

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20220422

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220817

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20221124

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20221124

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: 20230110

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230125

R150 Certificate of patent or registration of utility model

Ref document number: 7218288

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150