JP2023519217A - スピン標識化粒子の拡散伝搬関数又は関連する拡散パラメータを決定するための磁気共鳴法、ソフトウェア製品、及びシステム - Google Patents
スピン標識化粒子の拡散伝搬関数又は関連する拡散パラメータを決定するための磁気共鳴法、ソフトウェア製品、及びシステム Download PDFInfo
- Publication number
- JP2023519217A JP2023519217A JP2022556595A JP2022556595A JP2023519217A JP 2023519217 A JP2023519217 A JP 2023519217A JP 2022556595 A JP2022556595 A JP 2022556595A JP 2022556595 A JP2022556595 A JP 2022556595A JP 2023519217 A JP2023519217 A JP 2023519217A
- Authority
- JP
- Japan
- Prior art keywords
- diffusion
- magnetic resonance
- gradient
- sample
- spin
- 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.)
- Pending
Links
- 238000009792 diffusion process Methods 0.000 title claims abstract description 239
- 239000002245 particle Substances 0.000 title claims abstract description 95
- 238000001646 magnetic resonance method Methods 0.000 title claims abstract description 6
- 238000012546 transfer Methods 0.000 title description 33
- 230000005291 magnetic effect Effects 0.000 claims abstract description 171
- 238000000034 method Methods 0.000 claims abstract description 75
- 238000005259 measurement Methods 0.000 claims abstract description 53
- 230000010363 phase shift Effects 0.000 claims abstract description 21
- 238000002372 labelling Methods 0.000 claims abstract description 3
- 238000002595 magnetic resonance imaging Methods 0.000 claims description 26
- 238000009826 distribution Methods 0.000 claims description 21
- 238000004611 spectroscopical analysis Methods 0.000 claims description 21
- 239000000463 material Substances 0.000 claims description 15
- 238000004590 computer program Methods 0.000 claims description 14
- 239000006185 dispersion Substances 0.000 claims description 12
- 235000013305 food Nutrition 0.000 claims description 6
- 238000004113 cell culture Methods 0.000 claims description 4
- 238000009795 derivation Methods 0.000 claims description 4
- 238000013334 tissue model Methods 0.000 claims description 4
- 229920002994 synthetic fiber Polymers 0.000 claims description 3
- 239000007779 soft material Substances 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 abstract description 4
- 239000000523 sample Substances 0.000 description 90
- 230000006870 function Effects 0.000 description 66
- 210000001519 tissue Anatomy 0.000 description 49
- 239000011148 porous material Substances 0.000 description 27
- 238000005481 NMR spectroscopy Methods 0.000 description 13
- 238000003384 imaging method Methods 0.000 description 11
- 230000033001 locomotion Effects 0.000 description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 7
- 201000010099 disease Diseases 0.000 description 6
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 6
- 239000013598 vector Substances 0.000 description 6
- 230000007246 mechanism Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 239000013074 reference sample Substances 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000001727 in vivo Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000001131 transforming effect Effects 0.000 description 4
- 238000006073 displacement reaction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 210000000056 organ Anatomy 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 230000007480 spreading Effects 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 238000004435 EPR spectroscopy Methods 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000006399 behavior Effects 0.000 description 2
- 210000004027 cell Anatomy 0.000 description 2
- 210000000170 cell membrane Anatomy 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000006698 induction Effects 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 230000005415 magnetization Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000001575 pathological effect Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000011282 treatment Methods 0.000 description 2
- CVOFKRWYWCSDMA-UHFFFAOYSA-N 2-chloro-n-(2,6-diethylphenyl)-n-(methoxymethyl)acetamide;2,6-dinitro-n,n-dipropyl-4-(trifluoromethyl)aniline Chemical compound CCC1=CC=CC(CC)=C1N(COC)C(=O)CCl.CCCN(CCC)C1=C([N+]([O-])=O)C=C(C(F)(F)F)C=C1[N+]([O-])=O CVOFKRWYWCSDMA-UHFFFAOYSA-N 0.000 description 1
- 238000012307 MRI technique Methods 0.000 description 1
- 206010030113 Oedema Diseases 0.000 description 1
- 206010070834 Sensitisation Diseases 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002597 diffusion-weighted imaging Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 239000000839 emulsion Substances 0.000 description 1
- 210000003743 erythrocyte Anatomy 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000000338 in vitro Methods 0.000 description 1
- 239000012497 inhomogeneous sample Substances 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000000302 ischemic effect Effects 0.000 description 1
- 229920002521 macromolecule Polymers 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 210000004165 myocardium Anatomy 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000011896 sensitive detection Methods 0.000 description 1
- 230000008313 sensitization Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 210000000130 stem cell Anatomy 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000008733 trauma Effects 0.000 description 1
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
-
- 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/58—Calibration of imaging systems, e.g. using test probes, Phantoms; Calibration objects or fiducial markers such as active or passive RF coils surrounding an MR active material
-
- 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
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
- G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
-
- 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/46—NMR spectroscopy
-
- 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
-
- 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/58—Calibration of imaging systems, e.g. using test probes, Phantoms; Calibration objects or fiducial markers such as active or passive RF coils surrounding an MR active material
- G01R33/583—Calibration of signal excitation or detection systems, e.g. for optimal RF excitation power or frequency
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Health & Medical Sciences (AREA)
- Signal Processing (AREA)
- Vascular Medicine (AREA)
- Biomedical Technology (AREA)
- Pathology (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)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本開示は、試料中のスピン標識化粒子の拡散パラメータを決定するためのコンピュータ実装拡散磁気共鳴法に関する。この方法(100)は、試料及び磁気共鳴装置を用意するステップ(110)と、磁場勾配パルス列を試料に適用することによって、試料中の一組の粒子をスピン標識化するステップ(120)と、磁気共鳴装置により、スピン標識化粒子について、磁場勾配パルス列に対応する磁気共鳴測定データを取得するステップ(130)と、上記測定データに基づいて、スピン標識化粒子の拡散パラメータを決定するステップ(140)と、を含み、ステップ(140)が、拡散パラメータごとに、上記測定データに基づいて、拡散パラメータを表すフーリエ変換を構成することを含み、各磁場勾配パルス列が、少なくとも3つの勾配パルスを含み、勾配パルスが、試料中のそれぞれの位置に基づいて、スピン標識化粒子に位相シフトを導入するように構成されている。【選択図】 図1
Description
本開示は、磁気共鳴技術による試料の拡散パラメータの決定に関する。
NMRは、周知の技術として、化学的及び構造的環境に関する情報の推測に採用されている。磁気共鳴(MR)イメージング(MRI)は、空間中の異なる位置に異なるNMR信号値を割り当て可能な方法である。局所NMR分光法は、空間中の特定の体積からNMR信号が検出される技術である。
MRIを使用することにより、MR特性の相違に基づいて、試料及び生体組織の領域を識別可能であり、これがMRイメージにおける「コントラスト」となる。また、組織内の分子の運動(流れ及び拡散等)に基づいて、コントラストを生み出すこともできる。
拡散磁気共鳴分光法及び磁気共鳴イメージング(拡散MRI)では、多くの材料及び生体組織中の分子の動きを測定する。拡散MRIは、材料、生体組織、多孔質媒体を含む試料の構造を評価するための重要な技術となっている。
一対の磁場勾配パルスを一般的なMRパルス列に導入することによって、拡散重み付けMR測定を制御可能な方法で実行可能である。これらの拡散勾配がスピンエコーベースの方法に組み込まれた技術をパルス勾配スピンエコー(PGSE)と称する(E.O.Stejskal, and J.E.Tanner,J.Chem.Phys.,42,288,(1965)参照)。より一般的には、拡散勾配は、スピンエコーベースのパルス列並びに刺激エコー及び自由誘導減衰を検出する他のパルス列等に採用される。
本発明の目的は、磁気共鳴技術による試料の拡散パラメータの決定の改良である。
これは、本開示によれば、試料中のスピン標識化粒子の拡散パラメータを決定するためのコンピュータ実装拡散磁気共鳴法によって達成されている。この方法は、試料及び前記試料の磁気共鳴を測定するように構成された磁気共鳴装置を用意するステップと、少なくとも1つの磁場勾配パルス列を前記試料に適用することによって、前記試料に含まれる一組の粒子をスピン標識化するステップと、前記磁気共鳴装置により、前記スピン標識化粒子について、前記少なくとも1つの磁場勾配パルス列に対応する磁気共鳴(MR)測定データを取得するステップと、前記取得した測定データに基づいて、前記スピン標識化粒子の少なくとも1つの拡散パラメータを決定するステップと、を含む。
少なくとも1つの拡散パラメータを決定するステップは、拡散パラメータごとに、前記取得した測定データに基づいて、前記拡散パラメータを表す少なくとも1つのフーリエ変換を構成することを含む。各磁場勾配パルス列は、少なくとも3つの勾配パルスを含み、少なくとも1つの勾配パルスが、前記試料中のそれぞれの位置に基づいて、前記スピン標識化粒子に位相シフトを導入するように構成されている。
これにより、通常推定されるアンサンブル平均拡散伝搬関数(EAP)と比較して、拡散に関する付加的な情報を含む真の拡散伝搬関数等の拡散パラメータを推定できる利点がある。さらに、試料に関する過去未発見の情報にアクセス可能となり、その情報によって、データ及び当該データからのイメージを生成可能であるとともに、様々な数学的パラメータに基づいて新たなコントラストを得られる利点もある。これらの新たなコントラスト機構は、例えばヒト患者の疾患による変化を特徴付ける拡散MRの特異性及び有用性を向上可能な材料又は組織の微細構造に関する付加的な情報を提供するものとする。
拡散パラメータ(diffusion parameter)という用語は、ある体積中の拡散を表す値、一組の値、数式、又はモデルとして了解されるものとする。
いくつかの例において、少なくとも1つの拡散パラメータを決定するステップは、前記スピン標識化粒子の拡散伝搬関数、分散指数、局所的分散、拡散重み付けイメージ、定常状態分布、結合密度分布、及び/又は一組のキュムラントテンソルを決定することを含む。
いくつかの例において、少なくとも1つの磁場勾配パルス列を適用するステップは、2つの異なる時間における粒子の位置に基づいて、試料に含まれるスピン標識化粒子に2つの位相シフトを導入するように構成された勾配パルスを適用することを含む。
いくつかの例において、各勾配パルス列は、それぞれの持続時間が前記勾配パルス列の勾配パルス持続時間の合計の4分の1よりも短い少なくとも2つの勾配パルスを含み、前記少なくとも2つの勾配パルスは、前記試料中のそれぞれの位置に基づいて、前記スピン標識化粒子に位相シフトを導入するように構成され、並びに/又は、各勾配パルス列は、すべての勾配パルスに関する勾配磁場強度の実質的にゼロの時間積分を有する。
これにより、長くて弱い磁場勾配パルスと、勾配磁場強度の時間積分がゼロとなる2つの短くて強い磁場勾配パルスと、に曝露された試料が、2つの短い勾配パルス間で生じた拡散を示すMR測定データを提供できる利点がある。
いくつかの例においては、複数の勾配パルス列が適用され、各勾配パルス列は、勾配パルス間で異なる勾配強度、勾配方向、勾配持続時間、及び/又は遅延時間を有し、少なくとも1つの拡散パラメータを決定するステップは、前記複数の勾配パルス列に対応する測定データに基づいて、各拡散パラメータを表すフーリエ変換を構成することを含む。
これにより、(1つ又は複数の)付加的な勾配パルス列の結果として、前記拡散パラメータを表す前記フーリエ変換の構成の改善に使用し得る付加的な対応するMR測定データを得られる利点がある。
いくつかの例において、この方法は、決定した少なくとも1つの拡散パラメータに基づいて、前記磁気共鳴装置を校正するステップをさらに含む。
これにより、基準試料について決定された拡散伝搬関数等の拡散パラメータを前記基準試料に対応する基準データと比較して、誤った勾配のタイミング又は強度の提供等、磁気共鳴装置の性能を示し得る利点がある。
いくつかの例において、磁気共鳴測定データを取得するステップは、前記用意した磁気共鳴装置によって、磁気共鳴分光法、局所磁気共鳴分光法、及び/又は磁気共鳴イメージングを実行することを含む。
いくつかの例において、少なくとも1つの拡散パラメータを決定するステップは、離散フーリエ変換、高速フーリエ変換、及び/又は解析的導出を利用することを含む。
本開示はさらに、プログラム命令を含むコンピュータプログラムを有する非一時的なコンピュータ可読記憶媒体を含むコンピュータプログラム製品であって、コンピュータプログラムが、プロセッサにロード可能であり、スピン標識化粒子の拡散パラメータを決定するための方法をプロセッサに実行させるように構成された、コンピュータプログラム製品に関する。
本開示はさらに、スピン標識化粒子の拡散パラメータを決定するための磁気共鳴システムに関する。このシステムは、磁気共鳴装置及びコンピュータを備える。コンピュータは、
磁気共鳴装置によって、少なくとも1つの磁場勾配パルス列及び高周波パルス列を試料に与えることと、
磁気共鳴装置によって、複数の異なる勾配持続時間、勾配強度、及び勾配方向で前記試料から磁気共鳴測定データを収集することと、
を行うように前記磁気共鳴装置を制御するように構成されている。
磁気共鳴装置によって、少なくとも1つの磁場勾配パルス列及び高周波パルス列を試料に与えることと、
磁気共鳴装置によって、複数の異なる勾配持続時間、勾配強度、及び勾配方向で前記試料から磁気共鳴測定データを収集することと、
を行うように前記磁気共鳴装置を制御するように構成されている。
コンピュータは、
収集された磁気共鳴測定データに基づいて、少なくとも1つの拡散パラメータごとに、少なくとも1つのフーリエ変換を再構成することと、
対応する再構成した少なくとも1つのフーリエ変換に基づいて、少なくとも1つの拡散パラメータを決定することと、
を行うように構成され、
各勾配パルス列は、少なくとも3つの勾配パルスを含み、少なくとも1つの勾配パルスが、位置に基づいて、試料に含まれるスピン標識化粒子に位相シフトを導入するように構成されている。
収集された磁気共鳴測定データに基づいて、少なくとも1つの拡散パラメータごとに、少なくとも1つのフーリエ変換を再構成することと、
対応する再構成した少なくとも1つのフーリエ変換に基づいて、少なくとも1つの拡散パラメータを決定することと、
を行うように構成され、
各勾配パルス列は、少なくとも3つの勾配パルスを含み、少なくとも1つの勾配パルスが、位置に基づいて、試料に含まれるスピン標識化粒子に位相シフトを導入するように構成されている。
このシステムのいくつかの例において、コンピュータは、2つの異なる時間における粒子の位置に基づいて、試料に含まれるスピン標識化粒子に2つの位相シフトを導入するように構成された勾配パルス列を与えるように構成されている。
このシステムのいくつかの例において、コンピュータは、拡散測定結果に基づいて、磁気共鳴装置を校正するように構成されている。
このシステムのいくつかの例において、前記磁気共鳴装置は、磁気共鳴分光法、局所磁気共鳴分光法、及び/又は磁気共鳴イメージングを利用して、測定データを取得するように構成されている。
図面の全体を通して、同じ参照番号は、同じ部分、概念、及び/又は要素を表す。したがって、ある図面中の参照数字に関する言及は、別段の明示的な定めのない限り、他の図面中の同じ参照数字にも等しく当てはまる。
図1は、スピン標識化粒子の拡散パラメータを決定するための例示的なコンピュータ実装拡散磁気共鳴法を示している。この例示的な方法100は、
試料及び前記試料の磁気共鳴を測定するように構成された磁気共鳴装置を用意するステップ110と、
少なくとも1つの磁場勾配パルス列を前記試料に適用することによって、前記試料に含まれる一組の粒子をスピン標識化するステップ120と、
前記磁気共鳴装置により、前記スピン標識化粒子について、前記適用した勾配パルス列に対応する磁気共鳴測定データを取得するステップ130と、
前記取得した測定データに基づいて、前記スピン標識化粒子の少なくとも1つの拡散パラメータを決定するステップ140と、
を含み、
前記少なくとも1つの拡散パラメータを決定するステップ140が、拡散パラメータごとに、前記取得した測定データに基づいて、前記拡散パラメータを表す少なくとも1つのフーリエ変換を構成することを含み、
各磁場勾配パルス列が、少なくとも3つの勾配パルスを含み、少なくとも1つの勾配パルスが、前記試料中のそれぞれの位置に基づいて、前記スピン標識化粒子に位相シフトを導入するように構成されている。
試料及び前記試料の磁気共鳴を測定するように構成された磁気共鳴装置を用意するステップ110と、
少なくとも1つの磁場勾配パルス列を前記試料に適用することによって、前記試料に含まれる一組の粒子をスピン標識化するステップ120と、
前記磁気共鳴装置により、前記スピン標識化粒子について、前記適用した勾配パルス列に対応する磁気共鳴測定データを取得するステップ130と、
前記取得した測定データに基づいて、前記スピン標識化粒子の少なくとも1つの拡散パラメータを決定するステップ140と、
を含み、
前記少なくとも1つの拡散パラメータを決定するステップ140が、拡散パラメータごとに、前記取得した測定データに基づいて、前記拡散パラメータを表す少なくとも1つのフーリエ変換を構成することを含み、
各磁場勾配パルス列が、少なくとも3つの勾配パルスを含み、少なくとも1つの勾配パルスが、前記試料中のそれぞれの位置に基づいて、前記スピン標識化粒子に位相シフトを導入するように構成されている。
いくつかの例において、少なくとも1つの拡散パラメータを決定するステップ140は、前記スピン標識化粒子の拡散伝搬関数、分散指数、局所的分散、拡散重み付けイメージ、定常状態分布、結合密度分布、及び/又は一組のキュムラントテンソルを決定することを含む。
方法100のいくつかの例において、少なくとも1つの磁場勾配パルス列を適用するステップ120は、2つの異なる時間における粒子の位置に基づいて、試料に含まれるスピン標識化粒子に2つの位相シフトを導入するように構成された勾配パルスを適用することを含む。いくつかの例において、方法100は、少なくとも3つの異なる時間における粒子の位置に基づいて、試料に含まれるスピン標識化粒子に少なくとも2つの位相シフトを導入するように構成された磁場勾配パルス列を適用する(120)。通常、異なる時間は、時間的に分離された複数の磁場勾配パルスに対応する複数の異なる時間間隔等、異なる時間間隔に関する。
いくつかの例において、方法100は、決定した少なくとも1つの拡散パラメータに基づいて、磁気共鳴装置を校正するステップをさらに含む。これらの例のうちのいくつかにおいて、試料は、水が充填された容器等、対応する拡散パラメータデータが予め定められた基準試料である。決定した拡散パラメータと予め定められた対応する拡散パラメータデータとの間の不整合は、勾配タイミング又は勾配強度の所要の校正を示し得る。
いくつかの例において、方法100は、基準試料及び予め定められた対応する拡散パラメータデータを取得するステップと、前記基準試料の決定拡散パラメータ及び予め定められた対応する取得拡散パラメータデータに基づいて、磁気共鳴装置を校正するステップと、を含む。
少なくとも1つの勾配パルス列を適用するステップ120及び磁気共鳴測定データを取得するステップ130は、試料の磁気共鳴測定の一部であり、少なくとも1つの勾配パルス列を適用するステップ120及び磁気共鳴測定データを取得するステップ130は通常、相互依存しており、高周波パルス及び磁気勾配パルスの一致等、時間的に重なる場合もあることが了解されるものとする。磁気共鳴測定の実行方法に関する包括的な詳細については、本明細書の範囲外である。
いくつかの例において、この方法は、用意した磁気共鳴装置によって、試料の磁気共鳴測定を実行するステップ(図示せず)を含み、前記磁気共鳴測定を実行するステップは、少なくとも1つの勾配パルス列を適用すること120と、前記スピン標識化粒子の磁気共鳴測定データを取得すること130と、を含む。
いくつかの例において、各勾配パルス列は、位置に基づいて前記スピン標識化粒子に位相シフトを導入するように構成された少なくとも2つの勾配パルスを含む。
この方法のいくつかの例において、各勾配パルス列は、それぞれの持続時間が前記勾配パルス列の勾配パルス持続時間の合計の4分の1よりも短い少なくとも2つの勾配パルスを含む。これらの例のうちのいくつかにおいて、前記持続時間は、前記勾配パルス列の勾配パルス持続時間の合計の10分の1、30分の1、又は100分の1よりも短い。
いくつかの例においては、第1の勾配パルスの持続時間が長く、勾配強度が低く、第2の勾配パルス及び第3の勾配パルスの持続時間が大幅に短く、勾配強度が高いため、3つの全パルスを含めて時間積分した任意の方向に沿う勾配の成分が実質的にゼロである。これにより、第2の勾配パルス及び第3の勾配パルスによる拡散粒子の位相シフトによって、前記粒子の拡散伝搬関数及び定常状態分布と関連した信号減衰がもたらされ得る(式1参照)。
本例において、q及びq’はそれぞれ、粒子の平均位置により決まる位相シフトを導入するように構成された第2の勾配パルス及び第3の勾配パルスによって決まるベクトルである。EΔ(q,q’)は、信号減衰であって、拡散勾配が適用されない場合の信号で除した測定信号である。本例において、Δは、第2の勾配パルス及び第3の勾配パルスの開始時間差である。拡散伝搬関数P(x’,Δ|x)は、持続時間Δの時間間隔にわたって位置xの粒子がx’まで移動する確率密度を表す。ρ(x)は、定常状態分布である。
本例において、q及びq’はそれぞれ、粒子の平均位置により決まる位相シフトを導入するように構成された第2の勾配パルス及び第3の勾配パルスによって決まるベクトルである。EΔ(q,q’)は、信号減衰であって、拡散勾配が適用されない場合の信号で除した測定信号である。本例において、Δは、第2の勾配パルス及び第3の勾配パルスの開始時間差である。拡散伝搬関数P(x’,Δ|x)は、持続時間Δの時間間隔にわたって位置xの粒子がx’まで移動する確率密度を表す。ρ(x)は、定常状態分布である。
拡散プロセスがd次元で測定される場合、拡散伝搬関数は、信号減衰の2d次元逆フーリエ変換によりρ(x)の推定値に基づいて求められる。
ρ(x)は、以下の式によって、q’=0又はq=0のデータの部分集合を変換することにより推定される。
ρ(x)は、以下の式によって、q’=0又はq=0のデータの部分集合を変換することにより推定される。
いくつかの例において、前記拡散パラメータを決定するステップ140は、前記取得した測定データに基づいて、前記拡散パラメータの少なくとも1つのフーリエ変換を再構成することを含む。いくつかの例において、前記拡散パラメータを決定するステップ140は、前記取得した測定データに基づいて、前記拡散パラメータに対応する少なくとも1つのフーリエ変換を構成することを含む。いくつかの例において、前記拡散パラメータを決定するステップ140は、前記取得した測定データに基づいて、前記拡散パラメータを表す少なくとも1つのフーリエ変換を構成することを含む。いくつかの例において、前記拡散パラメータを決定するステップ140は、前記取得した測定データに基づいて、前記拡散パラメータと関連する少なくとも1つのフーリエ変換を構成することを含む。いくつかの例において、前記拡散パラメータを決定するステップ140は、前記拡散パラメータを表す少なくとも1つのフーリエ変換に対して、前記取得した測定データを適用することを含む。
いくつかの例において、前記少なくとも1つの拡散パラメータを決定するステップ140は、各拡散パラメータを表す数式を構成することを含み、各数式は、前記取得した測定データに基づく少なくとも2つのフーリエ変換を含む。これらの例のうちのいくつかにおいて、前記数式に含まれる少なくとも2つのフーリエ変換の少なくとも一方は、定常状態分布ρ(x)に対応する。いくつかの例において、少なくとも1つの拡散パラメータを決定するステップ140は、前記取得した測定データに基づいて、少なくとも2つのフーリエ変換を含む前記拡散伝搬関数を表す数式を構成することにより前記拡散伝搬関数を決定することを含む。これらの例のうちのいくつかにおいて、前記数式に含まれる前記少なくとも2つのフーリエ変換の一方は、定常状態分布ρ(x)に対応する。
いくつかの例において、少なくとも1つの磁場勾配パルス列を前記試料に適用するステップ120は、一組の勾配コイルを介して磁場勾配を適用することを含む。
この方法のいくつかの例においては、複数の勾配パルス列が適用され(120)、各勾配パルス列が、勾配パルス間で異なる勾配強度、勾配方向、勾配持続時間、及び/又は遅延時間を有する。
この方法のいくつかの例においては、前記取得した測定データに基づいて、少なくとも1つのフーリエ変換を含む前記拡散パラメータを表す数式を構成し、前記取得した測定データは、複数の勾配パルス列に対応する。異なる勾配パルス列に対応する測定データを利用することによって、フーリエ変換がより正確となるため、前記拡散パラメータを表す数式がより正確となり得る。
この方法のいくつかの例において、少なくとも1つの磁場勾配パルス列を適用するステップ120は、粒子の平均位置の分布又は関連するキュムラントテンソルを取得するように構成されたパルス列を適用することを含む。
いくつかの例において、磁気共鳴測定データを取得するステップ130は、前記用意した磁気共鳴装置によって、磁気共鳴分光法、局所磁気共鳴分光法、及び/又は磁気共鳴イメージングを実行することを含む。
いくつかの例において、少なくとも1つの拡散パラメータを決定するステップ140は、前記スピン標識化粒子の定常状態分布を決定することをさらに含む。これらの例のうちのいくつかにおいて、少なくとも1つの拡散パラメータを決定するステップ140は、前記決定した定常状態分布に基づいて、少なくとも1つの別の拡散パラメータを決定することを含む。
いくつかの例において、拡散パラメータを決定するステップ140は、解析的導出、離散フーリエ変換、及び/又は高速フーリエ変換を利用することを含む。
前記スピン標識化粒子を含む用意した試料は、多孔質媒体、軟質物質、食品、培養若しくは死後の生体組織、模擬材料、模擬組織、合成材料、人工組織モデル、組織スライス、及び/又は細胞培養物であってもよい。これらの例のうちのいくつかにおいて、用意した試料は、インビボ試料及び/又はエクスビボ試料である。
図2a及び図2bは、磁場勾配のパルス列に曝露された試料を模式的に示している。図示の例における磁場勾配対時間のプロットは、例示を目的としたものであり、スケールが不正確な場合もある。
図2aは、2つの磁場勾配パルスを含むパルス列に曝露された試料を模式的に示している。本例は、第1の磁場勾配パルス及び第2の磁場勾配パルスへの曝露において、一組の粒子291を含む試料290を示している。第1及び第2のパルスは、勾配矢印と、勾配磁場の強度及び方向を示す勾配磁場強度対時間のプロットと、によって表される。
本例において、試料290に含まれる一組の粒子291は、第1のパルスにおいて、実質的に同じ磁場強度を受けている。第1のパルス後の時間Δtの後、一組の粒子は、それぞれの位置に基づいて、第2のパルスによる磁場強度に曝露される。本例において、第1及び第2のパルスは、持続時間及び勾配強度が同じであり、唯一の相違点として、印加勾配磁場の方向が反対である。
本例において、試料290に含まれる一組の粒子291は、粒子A、B、及びCである。粒子Aは、パルス列のうち、強度及び持続時間が同じで方向が反対の磁場に曝露されている。前記勾配パルス列を用いて正しく設定された磁気共鳴測定の場合、粒子A及び類似の非可動部分は、勾配パルス列が適用されない測定と同じ信号を示し得る。粒子Bは、パルス列のうち、第1のパルスによって、第2のパルスよりも強度が高い磁場に曝露されている。粒子Cは、パルス列のうち、第2のパルスによって、第1のパルスよりも強度が高い磁場に曝露されている。本例において、絶対時間積分磁場強度の差は、粒子Cよりも粒子Bの方が大きく、これは、第2のパルスにおいて、粒子Bが磁場勾配の方向で初期位置から離れた位置にあることに対応する。前記勾配パルス列及び対応する適当な高周波パルスを用いて正しく設定された磁気共鳴測定の場合は、粒子B及び粒子Cの絶対時間積分磁場強度の結果としての差が検出され得る。
試料290全体にわたって磁場勾配が線形の本例の場合、一組の粒子291は、試料290における開始位置に関わらず、磁場勾配の方向における第1及び第2のパルスの位置の差に比例する絶対時間積分磁場強度に曝露されることが了解されるものとする。
図2bは、3つの磁場勾配のパルスを含むパルス列に曝露された試料290を模式的に示している。本例は、第1、第2、及び第3の磁場勾配パルスへの曝露において、一組の粒子291を含む試料290を示している。第1、第2、及び第3のパルスは、勾配矢印と、勾配磁場の強度及び方向を示す勾配磁場強度対時間のプロットと、によって表される。本例においては、図2aの例と同様に、勾配パルス列に対する勾配磁場強度の時間積分が実質的にゼロである。
図2bの例において、試料290に含まれる一組の粒子291は、第1のパルスに曝露されるが、この第1のパルスは、第2の勾配パルス及び第3の勾配パルスよりも持続時間が大幅に長く、勾配磁場強度が大幅に低い。第1のパルスの適用中、粒子291は、ランダムな軌跡を長距離移動する。第1のパルスの効果は、粒子291の平均位置に比例する位相シフトを導入することである。このパルスの持続時間が長いことから、同じ局所環境内で粒子291が移動するすべての軌跡の粒子291ごとの平均は通常、それらが含まれる局所構造292の中心周りに高密度に分布する。このため、同じ局所構造292内のすべての粒子が、試料内のそれぞれの局所環境の中心によって決まる同じ位相シフトを受ける。
図2bにおいて、図示の粒子A 291の軌跡は、長い持続時間の局所構造292の中心近くで経時的な平均位置が終端となり得る様子と、その結果、粒子291がこの持続時間にわたって同様の平均勾配磁場強度を受けることと、を示している。
本例において、第2のパルスの間、一組の粒子291は、当該第2のパルスにおける実質的に同じ磁場強度を受けている。第2のパルス後の時間Δtの後、一組の粒子は、それぞれの位置に基づいて、第3のパルスによる磁場強度に曝露される。本例において、第2のパルス及び第3のパルスは、持続時間及び勾配磁場方向が同じであり、唯一の相違点として、勾配磁場強度は、第3のパルスよりも第2のパルスの方が高い。
いくつかの例において、第1、第2、及び第3の磁場勾配パルスはそれぞれ、異なる方向を有する。いくつかの例において、パルス列は、それぞれが異なる方向を有する少なくとも3つの磁場勾配パルスを含む。これらの例のうちのいくつかにおいて、パルス列に対する経時的な磁場勾配は、実質的にゼロである。
第1のパラメータ例として、図2bのパルス列は、第1、第2、及び第3のパルスの持続時間がそれぞれ100ms、5ms、及び5msであり、第1及び第2のパルス間の時間が0.1msであり、第2及び第3のパルス間の時間が20msであってよい。本例において、第2及び第3のパルスの平均勾配磁場強度が第1のパルスの勾配磁場強度よりも実質的に10倍高いと、勾配パルス列に対する勾配磁場強度の時間積分が実質的にゼロとなる。
第2のパラメータ例として、図2bのパルス列は、第1、第2、及び第3のパルスの持続時間がそれぞれ800ms、0.1ms、及び0.1msであり、第1及び第2のパルス間の時間が0.1msであり、第2及び第3のパルス間の時間が5msであってよい。本例において、第2及び第3のパルスの平均勾配磁場強度が第1のパルスの勾配磁場強度よりも実質的に4000倍高いと、勾配パルス列に対する勾配磁場強度の時間積分が実質的にゼロとなる。
図3は、試料中の拡散を決定するための磁気共鳴ベースのシステムを模式的に示している。システム300は、磁気共鳴イメージングスキャナ等の磁気共鳴装置340及びコンピュータ350を備える。磁気共鳴装置340は、高周波(RF)送信器320、RF受信器330、及び磁場生成器310を備える。RF送信器320は、高周波パルスを試料390に与えるように構成されている。RF受信器330は、試料390中の歳差磁化を検出するように構成されている。磁場生成器310は、調整可能な磁場を試料390に与えるように構成されている。コンピュータは、磁気共鳴装置340から情報を取得するとともに、磁気共鳴装置340を制御するように構成されている。
コンピュータ350は、少なくとも1つの磁場勾配パルス列及び高周波パルス列を試料390に与えることと、複数の異なる勾配持続時間、勾配強度、及び勾配方向で前記試料390から磁気共鳴測定データを収集することと、を行うように磁気共鳴装置340を制御するように構成されている。コンピュータ350は、収集された磁気共鳴測定データに基づいて、少なくとも1つの拡散パラメータごとに、少なくとも1つのフーリエ変換を再構成することと、対応する再構成した少なくとも1つのフーリエ変換に基づいて、少なくとも1つの拡散パラメータを決定することと、を行うようにさらに構成され、
各勾配パルス列は、少なくとも3つの勾配パルスを含み、少なくとも1つの勾配パルスが、位置に基づいて、試料390に含まれるスピン標識化粒子に位相シフトを導入するように構成されている。
各勾配パルス列は、少なくとも3つの勾配パルスを含み、少なくとも1つの勾配パルスが、位置に基づいて、試料390に含まれるスピン標識化粒子に位相シフトを導入するように構成されている。
システム300のいくつかの例において、コンピュータ350は、2つの異なる時間における粒子の位置に基づいて、試料390に含まれるスピン標識化粒子に2つの位相シフトを導入するように構成された勾配パルス列を与えるように構成されている。
図4は、位置及び配向を決定するためのコンピュータプログラム製品を含むデータ処理ユニットを模式的に示している。図4は、非一時的なコンピュータ可読記憶媒体412を含むコンピュータプログラム製品を含むデータ処理ユニット410を示している。非一時的なコンピュータ可読記憶媒体412は、プログラム命令を含むコンピュータプログラムを有する。コンピュータプログラムは、データ処理ユニット410にロード可能であり、図3の説明に記載の位置及び配向を決定するための方法をプロセッサ411に実行させるように構成されている。
データ処理ユニット410は、磁気共鳴ベースのシステム等の装置400に含まれていてもよい。
図5は、磁場勾配パルス列(G)を示している。図5は、MR実験のための提案の3パルス磁場勾配を示している。この例示的なパルス列は、3つのパルス状磁場勾配を含み、q及びq’は、各パルスの持続時間を乗じた勾配強度と関連しており、Δは、第2のパルスの開始と第3のパルスの開始との間の時間である。本例においては、すべてのパルスについて持続時間を乗じた勾配強度の合計がゼロである。図5においては、プロットした勾配列が有効である。すなわち、特定の勾配ベクトルの方向は、実験に基づく設計において存在し得るRFパルスの影響を反映している。例えば、パルス列中に180度の高周波パルスが1つしか存在しない場合、前記高周波パルスの前に適用されるすべての勾配の方向は、当技術の実装において逆となる。
図6は、磁場勾配パルス列(G)及び高周波パルス列(RF)を示している。図6は、提案の技術の特定の一実施形態を示しており、短パルス状磁場勾配が長パルス状磁場勾配の前に適用される。さらに、各勾配が2つに分割されており、一方が180°高周波パルスの前、他方が180°高周波パルスの後である。ここでは、高周波パルスが誘導する位相シフトを考慮するため、180°高周波パルスの前後の部分が反対方向に適用される。本例において、Δは、第1の一対の勾配及び第2の一対の勾配の開始時間差であり、q及びq’は、各パルスの持続時間を乗じた勾配強度と関連する。
図7a~図7cは、試料中の取得信号及び決定拡散伝搬関数を示す等高線プロットを示している。図7aの例示的な等高線は、q及びq’に対してプロットされた信号である。図7bの等高線プロットは、図示のように、x及びx’に対してプロットされた決定拡散伝搬関数を示している。図7cの等高線プロットは、拡散伝搬関数の解析的形態を表しており、フックの復元力を受けた粒子の拡散伝搬関数の解析的表現である。
図8a~図8cは、試料中の決定拡散伝搬関数を示す表面プロットを示している。表面プロットは、Δの値を大きくして推定された拡散伝搬関数を表す。本例において、最後の2つのパルス間の時間は、図8a、図8b、及び図8cそれぞれについて、0.01ms、0.4ms、及び10msである。
図9a~図9cは、1次元細孔の拡散伝搬関数を模式的に示している。図9aには、特定の形状を示している。ここで、破線は、x方向に沿って10単位の細孔内の半透膜の位置を示している。勾配は、+x及び-x方向に適用される。図9bは、提案の方法を用いて推定された伝搬関数の双曲線正接の等高線プロットを示している。図9cは、パルス勾配スピンエコー(PGSE)(E.O.Stejskal, and J.E.Tanner,J.Chem.Phys.,42,288,(1965)参照)法を用いて取得したデータから推定されたアンサンブル平均伝搬関数を示している。
図10a~図10dは、5つの円形細孔を含む試料及び試料の不均質性と関連する数量を模式的に示している。図10aは、半径が12、36、60、84、及び108単位の5つの円形細孔を含む試料の分散指数の推定を示している。図10bは、中心を一致させた同じ細孔を示している。この技術では、図10bに示す系に関する測定を効果的に実行する。図10bの共通中心からの距離をrとする。図10c及び図10dは、提案の方法により得られる2つの数量をrに対してプロットしたものである。図10dに示すように、最初の数量を最大化するr値の採用により、分散指数を推定する。
本発明は、核磁気共鳴(NMR)及び磁気共鳴イメージング(MRI)に関し、より詳細には、試料内の拡散の測定並びにこれと関連するスペクトル及びイメージの生成を行う方法に関する。
複数の磁場勾配を適用することにより、拡散増感勾配強度(G)及び方向、拡散勾配パルス持続時間(δ)、並びに2つのパルスの適用の時間(Δ)に対するMR信号の依存性を調べる実験ツールが与えられる。これらのパラメータの変動に対するMR測定の感度の使用によって、観察サンプルから構造情報が推測されている。磁場勾配を採用する方法を用いて拡散プロセスを特徴付ける方法の1つとして、拡散分子の移動度と関連する拡散率(又は、拡散係数)を計算する。拡散伝搬関数P(x’,t|x)は、持続時間tの時間間隔にわたって位置xの粒子がx’まで移動する確率を表す。この数量は、単純な環境及び複雑な環境での拡散運動を完全に記述しており、制限壁又は半透壁、外力、空間的及び時間的に変化する拡散率等を特徴付け得る。
q空間NMR及びq空間MRIは、変位確率を非侵襲的に測定可能な分光及びイメージングモダリティである。この手法によれば、複数のq値(又は、ベクトル)で信号が観測され、qは、δ及びGに比例する数量である。信号値の結果的なアレイの逆フーリエ変換によって、アンサンブル平均伝搬関数が与えられるが、これは、空間中の粒子の平均的な正味の変位確率である。アンサンブル平均伝搬関数は、q、δ、G、Δ等の正味変位及び実験パラメータの関数として表される。アンサンブル平均伝搬関数は、拡散伝搬関数の位置変数のうちの1つの積分による拡散伝搬関数と関連する。
コンピュータ実装拡散磁気共鳴(拡散MR)システム及び方法では、材料、多孔質媒体、食品、及び生体組織等の試料中の拡散を検出することによって試料を検査する。アレイ状の複数の異なる拡散勾配持続時間、強度、及び方向が採用され、異なる時間に適用された異なる持続時間の勾配を特徴とするパルス列によって、拡散伝搬関数のフーリエ変換が集合的に行われる。そして、MR拡散データが試料から収集される。拡散伝搬関数は、コンピュータ上でデータを変換することにより演算される。このデータの使用によって、定常状態粒子分布及び他の尺度が推定される。計算された伝搬関数及び尺度の使用によって、試料の微細構造の変化の検出並びに臨床用途のMR研究の感度及び特異度の向上が可能なMRイメージングの新たなコントラストが生成され得る。
本発明の目的は、粒子の拡散運動を定量化することである。本発明の別の目的は、測定した拡散に含まれる情報を使用して、その時間的変化を特徴付ける新たなNMR及びMRI技術を提供することである。
本明細書においては、スピン標識化粒子の拡散伝搬関数を測定するための方法及びシステムを記載する。NMR分光計又はイメージャを使用して、NMR装置内のサンプルにパルス列を適用することによりNMR信号を生成し、そこから拡散伝搬関数及び関連する数量を計算することができる。MR拡散データは、拡散勾配持続時間、強度、及び/又は方向が異なるアレイ状の複数のパルス列において、組織から取得されるのが好ましい。コンピュータ実装拡散磁気共鳴(拡散MR)方法では、材料、多孔質媒体、又は生体組織等の試料中の拡散を検出することによって試料を解析する。
計算されたMR特性は、空間座標が異なる軸上に表示される多次元プロットで可視化されるのが好ましい。コンピュータの使用により、MR拡散データから、拡散伝搬関数を表す少なくとも1つの数量が計算される。推定された伝搬関数は、媒体の微細構造と関連する可能性がある。このような基礎となる拡散プロセスを表す数量の空間イメージを生成可能であるとともに、各ボクセル内の他の構造的特性と関連付け可能である。そして、異なる拡散勾配で実行された測定により得られた拡散データがプロットされ、解析される。伝搬関数は、等高線又は等高面で表示される。また、NMR装置自体の校正にも拡散測定結果が採用され得る。
図5、図7、及び図8の簡単な説明
これらの図面には、現時点で好適な実施形態を示すが、本発明は、その思想からも本質的属性からも逸脱することなく、他の形態で具現化可能であることが了解される。
これらの図面には、現時点で好適な実施形態を示すが、本発明は、その思想からも本質的属性からも逸脱することなく、他の形態で具現化可能であることが了解される。
図5は、拡散伝搬関数に対する感度を得るためのMRパルス列に組み込まれる勾配パルス列の一実施形態を示している。水平方向は、時間を示しており、特定の本実施形態においては、異なる時間に3つの勾配パルスが適用される。ボックスの高さは、実験の繰り返しで変化する勾配の強度を示している。高周波パルスが適用される場合は、類似の情報を実現するように勾配の方向が変化する。勾配パルスの一部が高周波(RF)パルスの前後に適用されるようになっていてもよい。勾配パルスのいずれか1つを分割することによって、当該勾配パルスの一部がRFパルスの前に適用され、別の部分がRFパルスの後に適用されるようにすることも可能である。
図7a~図7cは、等高線プロットとして表示される結果の一実施形態を示している。左側のパネルは、(Yolcu et al.,Phys Rev E,93,052602,2016)に記載のような復元力を受ける拡散について、拡散係数=2μm2/ms、閉じ込め=0.5μm-2、第1のパルスの持続時間=300ms、第1及び第2のパルス間の遅延=0.1ms、第2のパルス及び第3のパルスの持続時間=0.0001ms、最後の2つのパルス間の時間=1msといったパラメータで計算される信号を示している。中央のパネルは、本明細書に記載の発明を用いて上記信号プロファイルから演算される伝搬関数を示している。右側のパネルは、フックの復元力を受けた粒子の拡散伝搬関数の解析的表現を示している。
図8a~図8cは、表面プロットとして表示される結果の一実施形態を示している。拡散係数=2μm2/ms、閉じ込め=0.5μm-2、第1のパルスの持続時間=10ms、第1及び第2のパルス間の遅延=0.1ms、第2のパルス及び第3のパルスの持続時間=0.1ms、最後の2つのパルス間の時間=0.01ms(左側)、0.4ms(中央)、及び10ms(右側)で図7a~図7cと同じ方法が採用される。すべてのパネルにおいて、拡散伝搬関数は、本明細書に記載の発明を用いて演算されている。
コンピュータ実装拡散磁気共鳴(拡散MR)方法では、試料又は患者中の拡散を検出することによって、材料、多孔質媒体、及び生体組織を解析する。そして、数学的関係が与えられることにより、取得データから、粒子の動きを表す関数が演算される。また、MR拡散データが試料から収集される。MR拡散データは、拡散勾配持続時間、強度、及び方向が異なるアレイ状の複数の取得物によって、試料から取得されるのが好ましい。
MR拡散データは、磁気共鳴分光法、局所磁気共鳴分光法、又は磁気共鳴イメージング(MRI)を用いて収集可能である。複数の拡散勾配のうちの異なる勾配と関連付けられたデータ信号は、検査対象の試料における拡散プロセスの異なる特性に対応する。
定義
別段の定めのない限り、本明細書において使用するすべての技術的及び科学的用語は、本発明が属する技術分野の当業者が一般に理解する意味を有する。本明細書における使用の通り、以下の用語は、別段の指定のない限り、それぞれに固有の意味を有する。
別段の定めのない限り、本明細書において使用するすべての技術的及び科学的用語は、本発明が属する技術分野の当業者が一般に理解する意味を有する。本明細書における使用の通り、以下の用語は、別段の指定のない限り、それぞれに固有の意味を有する。
本明細書における使用の通り、磁気共鳴(MR)装置(magnetic resonance(MR)device)という用語は、磁気共鳴分光、MRイメージング、電子常磁性共鳴、電子スピン共鳴等の類似技術、又は同等技術が可能なすべての装置を含む。本発明の方法は、このような任意の装置、MR装置若しくは同等装置の変形の使用、又は任意既知のMR方法との併用により実施可能である。磁気共鳴方法及び装置においては、調査対象の試料、組織、又は生体への静磁場の印加によって、関心領域における平衡磁化軸の位置合わせが規定される。そして、静磁場方向と直交する方向での当該領域への高周波磁場の印加により、当該領域において磁気共鳴が励起される。空間情報が必要な場合は、慎重に設計された磁場勾配列の適用によって、局所的なMR分光及びMRイメージングが可能となる。結果としての信号は、関心生体の組織又はエリアに隣接して配置された高周波コイルにより検出される。
本明細書における使用の通り、用語「コンピュータ(computer)」及び「プロセッサ(processor)」は、それぞれの最も広い一般的背景で使用され、このようなすべての装置を含む。本発明の方法は、任意のコンピュータ/プロセッサの使用並びに任意既知のソフトウェア若しくは方法との併用により実施可能である。例えば、コンピュータ/プロセッサとしては、マイクロプロセッサ及びデータ転送バス等の従来の要素を含む従来の汎用デジタルコンピュータ(例えば、個人用「ワークステーション」コンピュータ)が可能である。コンピュータ/プロセッサはさらに、ダイナミックランダムアクセスメモリ、フラッシュメモリ等、又は磁気ディスク、光学ストレージ等の大容量ストレージ等、如何なる形態のメモリ素子をも具備し得る。
本発明の方法は、MRI拡散重み付けイメージング(DWI)を含む。簡潔に、この手法は、分子のランダム運動の測定及び粒子の拡散移動を定量化可能な核磁気共鳴に基づく。拡散イメージングは、この拡散測定をMRIと組み合わせた方法である。この技術では、イメージの各絵素(ピクセル又はボクセル)においてスピン標識化分子の拡散特性を特徴付けることができる。スピン標識化分子の拡散の特性は、化学的及び幾何学的環境と関連する。例えば、拡散イメージングの使用によって、水分子の運動を制限又は妨害する微細構造(例えば、細胞膜及び大型高分子)に関する情報を推測可能である。結果として、拡散イメージングにより、脳等の高次臓器における水の拡散を検出可能である。これらの組織においては、細胞膜による水分子の制限のため、水が全方向に等しく拡散することはなく、異方性拡散と称する特性となる。
特殊な勾配パルス列で得られたMRデータを考察することにより、拡散伝搬関数を推定することによって試料内の拡散を特徴付けられることが分かっている。この得られた拡散伝搬関数は、異なる現象を表し得る。例えば、図5に示す実施形態において、第1のパルスが非常に長く、第2及び第3のパルスが非常に短い場合は、拡散伝搬関数の求めた推定値によって、実際の拡散伝搬関数を正確に表すことができる。パルス持続時間が前記条件と異なる場合も本発明を依然として採用可能であり、求めた数量は、見掛けの拡散伝搬関数を表しており、真の拡散伝搬関数を表さない場合もある。例えば、図5の第2及び第3のパルスの持続時間が有限の場合、推定拡散伝搬関数の独立変数x及びx’は、勾配の持続時間にわたって平均化された粒子位置を表す。
このように拡散を測定する新たな方法によって、試料に関する過去未発見の情報が提供され、その情報によって、データ及び当該データからのイメージを生成可能であるとともに、様々な数学的パラメータに基づいて新たなコントラストを得られる。これらの新たなコントラスト機構は、例えばヒト患者の疾患による変化を特徴付ける拡散MRの特異性及び有用性を向上可能な材料又は組織の微細構造に関する付加的な情報を提供するものとする。この方法は、インビトロ、エクスビボ、又はインビボで実施可能である。
MR磁石及びその操作/解析ソフトウェアを含むMRシステムは、複数の異なる拡散勾配、勾配持続時間、及び拡散時間で試料からMR信号を収集するように構成可能である。このようなMR信号は、ヒト又は動物の特定の材料又は組織内の拡散伝搬関数に変換可能である。このシステム及びコンピュータソフトウェアは、コンピュータを用いた拡散の取得、体系化、解析、及びパラメータ化が可能である。
このシステムは、拡散データの取得、体系化、解析、及びパラメータ化、並びに拡散伝搬関数への変換を行う磁気共鳴イメージング又は分光システム及びソフトウェアとして具現化可能である。拡散MRIへの適用により、現在の拡散MRIコントラスト機構によって、多くの疾患の高感度検出を提供可能ではあるものの、異なる処置が必要な特定の疾患経過(例えば、腫瘍と浮腫)を区別できないことが多い。本発明に係るシステムは、特定の疾患経過に対してより優れた特異性を提供する複数の新たなMRI拡散コントラスト機構を提供し得る。したがって、本発明は、臨床患者を評価する拡散MRIの実用性を向上させ得る。
本発明は、組織中の拡散の変化を検出することにより、過去に不可能であったコントラストを提供して、これらの材料組織及び臓器の微細構造の変化を検出する。また、本発明の方法は、細胞の集合体又は組織の生理状態又は健康状態の変化をモニタリングする方法を提供する。これらの方法は、処置の有効性又は成長及び老化等の過程のモニタリングに特に有用となり得る。本発明に基づくコントラスト機構は、ヒト患者の疾患を特徴付ける拡散MRIの臨床的特異性及び臨床的有用性を向上可能な微細構造に関する付加的な情報を提供するものとする。
この方法は、空間中の異なる点間の移動の尤度を検出することによって、繊維状試料の接続パターンの評価に役立ち得る。特に、この方法によって、交差型、スプレイ型、キス型の繊維組成を区別することが可能と考えられる。
多孔質媒体の構造の変化を検出することによって、媒体の空隙率及び連通性を評価することも可能であり、このような媒体からの水及び炭化水素等の流体の回収し易さを評価するのに役立つと考えられる。
この方法は、懸濁液又は乳化液内の小球の内部及び周囲の拡散を検出することによって、食品の品質の評価に役立つと考えられる。
本発明は、ヒト、動物、又は植物の如何なる生体組織(例えば、ラット、マウス、又はヒトの神経組織、心筋組織、又は筋肉組織)にも適用可能である。データは、インビボ、エクスビボ、培養、又は死後の組織から取得可能である。また、本願は、赤血球のゴースト、組織スライス、(生体組織に由来する細胞を使用した)細胞培養物等、模擬組織又は人工組織モデルにも使用可能である。
また、本発明の使用によって、病理学的に損傷を受けた組織から健常な組織を識別することができる。病理学的組織としては、癌組織、虚血組織、外傷組織、慢性損傷組織、又は変性組織が挙げられるが、これらに限定されない。また、本願は、遺伝子組み換え組織の研究にも使用可能である。また、本願は、組織代替物又は置換物(幹細胞、臓器移植、又は3D骨格により形成された組織等)における拡散及び正常な健常組織における拡散の類似性の認識又は特性化にも使用可能である。
また、本発明の使用によって、MR装置自体を校正することができる。これは、自由水等の単純な系で拡散伝搬関数を測定することにより実現され得る。この場合、拡散プロセスは、ガウス分布になることが知られている。したがって、このようなサンプルから一連の拡散重み付けデータを取得することにより、勾配タイミングの精度に関する正確な情報を与えることができる。
このようなサンプルから一連の拡散重み付けデータを取得することにより、適用方法の精度及びハードウェア性能に関する情報を決定可能である。
拡散を決定する例
以下の項は、例示のための図5のパルス列に基づくが、特許請求の範囲に係る発明を制限するものではない。
以下の項は、例示のための図5のパルス列に基づくが、特許請求の範囲に係る発明を制限するものではない。
波形の積分がゼロとなり、試料を構成するすべてのドメインからの寄与がサンプル内の位置に依存しないように、持続時間の長いパルスが採用される。2つの後続パルスq及びq’は、それぞれx及びx’により示される粒子の平均位置(パルスの持続時間にわたる平均化)によって決まる位相シフトを導入する。qは、d個の成分を有するが、これは、拡散プロセスが測定される次元の数である。そして、信号減衰(拡散勾配が適用されない場合の信号で除した実際の信号)は、以下の関係により、伝搬関数P(x’,Δ|x)及び定常状態分布ρ(x)と関連する。
伝搬関数は、ρ(x)の推定値と併せて、信号の2d次元逆フーリエ変換により求められるが、
これは、長い拡散時間での伝搬関数であって、q’=0のデータ又はq=q’の長いΔで得られたデータの部分集合を以下の式で変換することにより得られる。
これは、長い拡散時間での伝搬関数であって、q’=0のデータ又はq=q’の長いΔで得られたデータの部分集合を以下の式で変換することにより得られる。
N個の分離した細孔を含む構造的に不均質な試料の場合は、EΔ(q,q’)の2d次元逆フーリエ変換によって、以下が得られる。
ここで、fnは、勾配がオフとなった場合にn番目の細孔が信号に寄与する割合を示す。長いΔでは、x=x’の設定によって、数量
が求められる。一方、EΔ(q,0)のd次元逆フーリエ変換は、
である。単分散系の場合、最後の数量の2乗は過去の数量と等しく、それぞれの非類似性の採用によって、すべての細孔の質量中心が同じ点となる場合の空間中の異なる位置での分散を測定することができる。
ここで、fnは、勾配がオフとなった場合にn番目の細孔が信号に寄与する割合を示す。長いΔでは、x=x’の設定によって、数量
が求められる。一方、EΔ(q,0)のd次元逆フーリエ変換は、
である。単分散系の場合、最後の数量の2乗は過去の数量と等しく、それぞれの非類似性の採用によって、すべての細孔の質量中心が同じ点となる場合の空間中の異なる位置での分散を測定することができる。
拡散伝搬関数は、ρ(x)の推定値と併せて、信号
の2d次元逆フーリエ変換により求められるが、これは、長い拡散時間での拡散伝搬関数であって、q’=0又はq=0のデータの部分集合を以下の式で変換することにより得られる。
の2d次元逆フーリエ変換により求められるが、これは、長い拡散時間での拡散伝搬関数であって、q’=0又はq=0のデータの部分集合を以下の式で変換することにより得られる。
こうして、拡散伝搬関数は、以下により求められる。
提案の実験技術の使用により、試料の微視的構造に関連する複数の数量を推定可能である。例えば、数量W(x’,Δ;x,0)は、ある時間にある位置xに存在し、持続時間Δ後に位置x’に存在する粒子の結合分布を示す。この結合密度は、本発明の使用により、以下のフーリエ関係を通じて容易に得られる。
提案の実験技術の使用により、試料の微視的構造に関連する複数の数量を推定可能である。例えば、数量W(x’,Δ;x,0)は、ある時間にある位置xに存在し、持続時間Δ後に位置x’に存在する粒子の結合分布を示す。この結合密度は、本発明の使用により、以下のフーリエ関係を通じて容易に得られる。
qの要素にq’の要素が後続する2d次元ベクトルKを構成可能である。同様に、xの要素にx’の要素が後続する2d次元ベクトルYを構成可能である。信号減衰又はその対数は、Kのマクローリン級数で展開可能である。例えば、以下の方程式によって、信号減衰の低K挙動が決まる。
ここで、Kkl=〈YkYl〉、Kklm=〈YkYlYm〉、及びKklmn=〈YkYlYmYn〉-3〈YkYl〉〈YmYn〉は、時間間隔Δによって分離された2つの時間における粒子の位置のキュムラントを表すテンソルである。これらのキュムラントテンソルは、最小2乗推定等の多くの数値モデルフィッティング技術のうちの1つを採用することにより、本発明を用いて収集されたデータから推定可能である。
ここで、Kkl=〈YkYl〉、Kklm=〈YkYlYm〉、及びKklmn=〈YkYlYmYn〉-3〈YkYl〉〈YmYn〉は、時間間隔Δによって分離された2つの時間における粒子の位置のキュムラントを表すテンソルである。これらのキュムラントテンソルは、最小2乗推定等の多くの数値モデルフィッティング技術のうちの1つを採用することにより、本発明を用いて収集されたデータから推定可能である。
別の例は、N個の細孔を含む構造的に不均質な試料に関する。この場合、n番目の細孔における粒子の定常状態分布は、ρn(x)により示し得る。
さらに、fnは、勾配がオフとなった場合にn番目の細孔が信号に寄与する割合を示す。EΔ(q,0)又はEΔ(0,q’)のd次元逆フーリエ変換は、
であり、積分して1になる。一方、EΔ(q,q’)の2d次元逆フーリエ変換によって、以下が得られる。
であり、積分して1になる。一方、EΔ(q,q’)の2d次元逆フーリエ変換によって、以下が得られる。
「局所的分散」と称し得る数量は、以下により規定され得る。
この数量の平方根(「局所的標準偏差」)は、以下により無次元化され得る。
ここで、x0は、以下をゼロとしない任意の点とすることができる。
x0の特定の選定では、以下が最大化される。
これらのマップは、すべての細孔の質量中心が同じ点となる場合に、試料の空間依存的な不均質性を表すものと考えられる。
この数量の平方根(「局所的標準偏差」)は、以下により無次元化され得る。
ここで、x0は、以下をゼロとしない任意の点とすることができる。
x0の特定の選定では、以下が最大化される。
これらのマップは、すべての細孔の質量中心が同じ点となる場合に、試料の空間依存的な不均質性を表すものと考えられる。
多くのシナリオにおいて、n番目の細孔の定常状態分布ρn(x)は、その値が細孔容積の逆数となる細孔空間の指標関数である。すなわち
であって、細孔空間の外側でゼロとなる。緩和の相違の可能性を無視すると、n番目の細孔の信号割合は、fn=Vn/(ΣnfnVn)である。これの条件下において、μ2(x)を積分すると、〈V〉-1により示される平均細孔容積の逆数となり、また、質量中心を原点とするすべての細孔内にx0が存在する場合は、
に等しくなる。この点は、以下を最大化する任意の点として決定可能である。
一方、μ2(x0)=〈V-1〉〈V〉-1である。特に注目すべきは、数量
であり、これは、すべての細孔が(それぞれの配向に関わらず)同じ容積の場合にゼロとなり、それ以外の場合には大きくなるため、「分散指数」と称する。これにより、紹介の方法で取得されたデータから、以下の式を用いてDIを演算可能である。
であって、細孔空間の外側でゼロとなる。緩和の相違の可能性を無視すると、n番目の細孔の信号割合は、fn=Vn/(ΣnfnVn)である。これの条件下において、μ2(x)を積分すると、〈V〉-1により示される平均細孔容積の逆数となり、また、質量中心を原点とするすべての細孔内にx0が存在する場合は、
に等しくなる。この点は、以下を最大化する任意の点として決定可能である。
一方、μ2(x0)=〈V-1〉〈V〉-1である。特に注目すべきは、数量
であり、これは、すべての細孔が(それぞれの配向に関わらず)同じ容積の場合にゼロとなり、それ以外の場合には大きくなるため、「分散指数」と称する。これにより、紹介の方法で取得されたデータから、以下の式を用いてDIを演算可能である。
データから演算可能なさらに別の数量が、数量の分布である。
これは、2つの短いパルスの適用中の粒子の位置の平均である。uの分布は、以下となるように規定される。
これは、以下の逆フーリエ変換によって、データから求められる。
これは、2つの短いパルスの適用中の粒子の位置の平均である。uの分布は、以下となるように規定される。
これは、以下の逆フーリエ変換によって、データから求められる。
これは、Q/2=q=q’の場合の提案パルス列の特定の実現に関する。信号減衰又はその対数は、qのマクローリン級数で展開可能である。例えば、以下の方程式によって、信号減衰の低Q挙動が決まる。
ここで、Ckl=〈ukul〉、Cklm=〈ukulum〉、及びCklmn=〈ukulumun〉-3〈ukul〉〈umun〉は、粒子の平均位置のキュムラントを表すテンソルである。これらのキュムラントテンソルは、最小2乗推定等の多くの数値モデルフィッティング技術のうちの1つを採用することにより、本発明を用いて収集されたデータから推定可能である。
ここで、Ckl=〈ukul〉、Cklm=〈ukulum〉、及びCklmn=〈ukulumun〉-3〈ukul〉〈umun〉は、粒子の平均位置のキュムラントを表すテンソルである。これらのキュムラントテンソルは、最小2乗推定等の多くの数値モデルフィッティング技術のうちの1つを採用することにより、本発明を用いて収集されたデータから推定可能である。
磁気共鳴データ取得
MRデータ取得の選択肢は、上記実験に限定されない。例えば、分光法、局所分光法、及びイメージング法も採用可能である。長短パルスのオーダを変更することも可能である。各パルスを複数のパルスで置き換えることも可能である。これらのパルスは、それぞれの適用中に一定である必要がなく、強度が時間的に変動し得る。この場合、q及びq’の値は、以下の式により規定される。
q=γ∫G(t)dt (式17)
q’=γ∫G’(t)dt (式18)
ここで、γは、検査対象の粒子のジャイロ磁気(磁気回転)比であり、G(t)及びG’(t)は、それぞれの勾配パルスに対する有効勾配ベクトル(高周波パルスの影響を考慮した後)であり、各勾配の持続時間にわたって積分が計算される。
MRデータ取得の選択肢は、上記実験に限定されない。例えば、分光法、局所分光法、及びイメージング法も採用可能である。長短パルスのオーダを変更することも可能である。各パルスを複数のパルスで置き換えることも可能である。これらのパルスは、それぞれの適用中に一定である必要がなく、強度が時間的に変動し得る。この場合、q及びq’の値は、以下の式により規定される。
q=γ∫G(t)dt (式17)
q’=γ∫G’(t)dt (式18)
ここで、γは、検査対象の粒子のジャイロ磁気(磁気回転)比であり、G(t)及びG’(t)は、それぞれの勾配パルスに対する有効勾配ベクトル(高周波パルスの影響を考慮した後)であり、各勾配の持続時間にわたって積分が計算される。
拡散伝搬関数の演算に必要なMRデータは、(勾配の強度、方向、及び/又は拡散パルス持続時間の変更により得られる)複数のq及びq’の値を有する拡散増感取得により得られる。q空間の特定のサンプリング方式は必要なく、多くの異なる代替手段の採用により、本発明に従って伝搬関数を求め、拡散プロセスを特徴付けることも可能である。データ取得は、複数の勾配分離(Δ等)及び持続時間で実行可能である。
使用するパルス列は、スピンエコー、刺激エコー、勾配エコー、又は自由誘導減衰を利用するもの等、多くの可能性のうちの1つから選択可能である。また、本発明に従って、エコープラナーイメージング等の高速イメージング技術を採用可能である。
拡散伝搬関数の構成
試料中の拡散の伝搬関数を演算するプロセスには、観測された信号減衰の変換の評価を含む。前述の通り、上述の再構成方式には、逆フーリエ変換を含む。フーリエ変換並びに正弦変換及び余弦変換等の関連する変換は、それぞれの逆変換と併せて、コンピュータ上での演算又はいくつかのシステムに対する解析的評価が可能である。コンピュータ上では、離散フーリエ変換(DFT)及び高速フーリエ変換(FFT)等の関連するアルゴリズム並びにそれぞれの逆変換によって、変換を実行可能である。
試料中の拡散の伝搬関数を演算するプロセスには、観測された信号減衰の変換の評価を含む。前述の通り、上述の再構成方式には、逆フーリエ変換を含む。フーリエ変換並びに正弦変換及び余弦変換等の関連する変換は、それぞれの逆変換と併せて、コンピュータ上での演算又はいくつかのシステムに対する解析的評価が可能である。コンピュータ上では、離散フーリエ変換(DFT)及び高速フーリエ変換(FFT)等の関連するアルゴリズム並びにそれぞれの逆変換によって、変換を実行可能である。
例示的な実施形態
第1の例は、拡散を検出することによって、材料、多孔質媒体、食品、及び生体組織等の試料を解析するための拡散磁気共鳴(拡散MR)方法であって、
MRスキャナを制御するコンピュータ上で実行されるパルス列のコンピュータプログラムを用意するステップであり、異なる持続時間、強度、及び/又は方向の勾配パルスの適用によって、異なる時間の粒子の位置を符号化することにより、データが取得される、ステップと、
フーリエ変換によって拡散伝搬関数に関連付ける拡散信号のモデルを用意するステップと、
前記試料からMR拡散データを収集し、コンピュータを用いて、前記取得手順を前記試料に適用することにより取得された前記データを関連付ける少なくとも1つのフーリエ変換、正弦変換、若しくは余弦変換、又はそれぞれの逆変換を計算するステップと、
コンピュータを用いて、試料の拡散特性と関連する尺度を計算するステップと、
を含む、方法である。
第1の例は、拡散を検出することによって、材料、多孔質媒体、食品、及び生体組織等の試料を解析するための拡散磁気共鳴(拡散MR)方法であって、
MRスキャナを制御するコンピュータ上で実行されるパルス列のコンピュータプログラムを用意するステップであり、異なる持続時間、強度、及び/又は方向の勾配パルスの適用によって、異なる時間の粒子の位置を符号化することにより、データが取得される、ステップと、
フーリエ変換によって拡散伝搬関数に関連付ける拡散信号のモデルを用意するステップと、
前記試料からMR拡散データを収集し、コンピュータを用いて、前記取得手順を前記試料に適用することにより取得された前記データを関連付ける少なくとも1つのフーリエ変換、正弦変換、若しくは余弦変換、又はそれぞれの逆変換を計算するステップと、
コンピュータを用いて、試料の拡散特性と関連する尺度を計算するステップと、
を含む、方法である。
第2の例は、前記MR拡散データが、アレイ状の複数の異なる拡散、勾配強度、勾配方向、勾配持続時間、及び拡散時間を特徴とする前記パルス列を用いることにより、前記試料から取得される、上記第1の例である。
第3の例は、前記MR拡散データが、磁気共鳴イメージング、磁気共鳴分光法、又は局所磁気共鳴分光法を用いることにより、前記収集ステップにおいて収集される、上記第1の例である。
第4の例は、前記パルス列が、拡散の測定に用いられる少なくとも3つの勾配パルスを含み、前記勾配パルスのうちの1つが、前記勾配パルスのうちの別の1つよりも長い、上記第1の例である。
第5の例は、前記フーリエ変換が、拡散伝搬関数又は関連する数量を取得するために、解析的導出、離散フーリエ変換、又は高速フーリエ変換を用いて計算される、上記第1の例である。
第6の例は、前記パルス列が、ドメイン形状の分散を含む尺度を取得するために採用される、上記第1の例である。
第7の例は、前記パルス列が、粒子の定常状態分布を取得するために採用される、上記第1の例である。
第8の例は、前記コンピュータを用いて、前記複数の拡散勾配持続時間、勾配強度、勾配方向、及び時間で収集された前記MRデータから、前記組織をイメージングするステップをさらに含み、前記複数の前記拡散勾配は、拡散伝搬関数のフーリエ変換を測定する上記第2の例である。
第9の例は、前記MR拡散データが、材料、多孔質媒体、軟質物質、食品、インビボ、エクスビボ、培養若しくは死後の生体組織、模擬材料、模擬組織、合成材料、人工組織モデル、組織スライス、又は細胞培養物から収集される、上記第1の例である。
第10の例は、前記MR拡散データが、病理学的組織からの健常組織の識別、遺伝子組み換え組織、組織代替物、及び組織置換物の研究に用いられる、上記第1の例である。
第11の例は、前記拡散伝搬関数及び関連する尺度が、前記MR装置を校正するために用いられる、上記第1の例である。
第12の例は、材料試料及び生体組織中の拡散を解析するための拡散MRベースのシステムであって、
磁場勾配パルスと高周波パルスとの列を前記組織に与えて、複数の異なる拡散勾配持続時間、拡散勾配強度、及び拡散勾配方向で前記試料からMR拡散データを収集するMR装置であり、前記複数の拡散時間のうちの異なる時間が、拡散の時間依存性を示す、MR装置と、
コンピュータと、
を備え、
前記MRI装置が、パルス列の実行のため、前記コンピュータに通信可能に接続され、前記MRI装置が、前記MR装置により収集された前記拡散データを解析するとともに前記拡散データをフーリエ変換に適用するため、前記コンピュータに通信可能に接続され、前記拡散伝搬関数が、再構成され、前記MRI装置が、前記MR装置により収集された前記拡散データを解析するとともに前記拡散データを適用して分散の尺度を演算するため、前記コンピュータに通信可能に接続された、システムである。
磁場勾配パルスと高周波パルスとの列を前記組織に与えて、複数の異なる拡散勾配持続時間、拡散勾配強度、及び拡散勾配方向で前記試料からMR拡散データを収集するMR装置であり、前記複数の拡散時間のうちの異なる時間が、拡散の時間依存性を示す、MR装置と、
コンピュータと、
を備え、
前記MRI装置が、パルス列の実行のため、前記コンピュータに通信可能に接続され、前記MRI装置が、前記MR装置により収集された前記拡散データを解析するとともに前記拡散データをフーリエ変換に適用するため、前記コンピュータに通信可能に接続され、前記拡散伝搬関数が、再構成され、前記MRI装置が、前記MR装置により収集された前記拡散データを解析するとともに前記拡散データを適用して分散の尺度を演算するため、前記コンピュータに通信可能に接続された、システムである。
図5は、提案のMR実験の一実施形態である。3つのパルス状磁場勾配(灰色のボックスで示す)が採用され、q及びq’は、各パルスの持続時間を乗じた勾配強度と関連する。
図7は、代表的な解析を示した等高線プロットである。q及びq’に対する信号の等高線プロットを左側に示す。本発明に記載の技術がこのデータに適用されると、中央に示すように、x及びx’に対してプロットされた伝搬関数の推定値が生成される。伝搬関数の解析的形態を使用して、右側の等高線プロットを生成している。
図8は、(左側から右側へと)Δの値を大きくして推定された拡散伝搬関数を表す表面プロットである。
Claims (14)
- 試料中のスピン標識化粒子の拡散パラメータを決定するためのコンピュータ実装拡散磁気共鳴法(100)であって、
試料及び前記試料の磁気共鳴を測定するように構成された磁気共鳴装置を用意するステップ(110)と、
少なくとも1つの磁場勾配パルス列を前記試料に適用することによって、前記試料に含まれる一組の粒子をスピン標識化するステップ(120)と、
前記磁気共鳴装置により、前記スピン標識化粒子について、前記少なくとも1つの磁場勾配パルス列に対応する磁気共鳴測定データを取得するステップ(130)と、
前記取得した測定データに基づいて、前記スピン標識化粒子の少なくとも1つの拡散パラメータを決定するステップ(140)と、
を含み、
前記少なくとも1つの拡散パラメータを決定する前記ステップ(140)が、拡散パラメータごとに、前記取得した測定データに基づいて、前記拡散パラメータを表す少なくとも1つのフーリエ変換を構成することを含み、
各磁場勾配パルス列が、少なくとも3つの勾配パルスを含み、少なくとも1つの勾配パルスが、前記試料中のそれぞれの位置に基づいて、前記スピン標識化粒子に位相シフトを導入するように構成された、方法。 - 少なくとも1つの拡散パラメータを決定する前記ステップ(140)が、前記スピン標識化粒子の拡散伝搬関数、分散指数、局所的分散、拡散重み付けイメージ、定常状態分布、結合密度分布、及び/又は一組のキュムラントテンソルを決定することを含む、請求項1に記載の方法。
- 各勾配パルス列が、それぞれの持続時間が前記勾配パルス列の勾配パルス持続時間の合計の4分の1よりも短い少なくとも2つの勾配パルスを含み、前記少なくとも2つの勾配パルスが、前記試料中のそれぞれの位置に基づいて、前記スピン標識化粒子に位相シフトを導入するように構成された、請求項1又は2に記載の方法。
- 各勾配パルス列が、すべての勾配パルスに関する勾配磁場強度の実質的にゼロの時間積分を有する、請求項1~3のいずれか一項に記載の方法。
- 複数の勾配パルス列が適用され、各勾配パルス列が、勾配パルス間で異なる勾配強度、勾配方向、勾配持続時間、及び/又は遅延時間を有し、前記少なくとも1つの拡散パラメータを決定する前記ステップ(140)が、前記複数の勾配パルス列に対応する測定データに基づいて、各拡散パラメータを表すフーリエ変換を構成することを含む、請求項1~4のいずれか一項に記載の方法。
- 磁気共鳴測定データを取得する前記ステップ(130)が、前記用意した磁気共鳴装置によって、磁気共鳴分光法、局所磁気共鳴分光法、及び/又は磁気共鳴イメージングを実行することを含む、請求項1~5のいずれか一項に記載の方法。
- 前記少なくとも1つの拡散パラメータを決定する前記ステップ(140)が、離散フーリエ変換、高速フーリエ変換、及び/又は解析的導出を利用することを含む、請求項1~6のいずれか一項に記載の方法。
- 前記少なくとも1つの拡散パラメータを決定する前記ステップ(140)が、拡散パラメータごとに、前記拡散パラメータを表す数式を構成するステップを含み、各数式が、前記取得した測定データに基づいた少なくとも2つのフーリエ変換を含む、請求項1~7のいずれか一項に記載の方法。
- 前記スピン標識化粒子を含む前記用意した試料が、多孔質媒体、軟質物質、食品、培養若しくは死後の生体組織、模擬材料、模擬組織、合成材料、人工組織モデル、組織スライス、及び/又は細胞培養物である、請求項1~8のいずれか一項に記載の方法。
- 前記決定した少なくとも1つの拡散パラメータに基づいて、前記磁気共鳴装置を校正するステップをさらに含む、請求項1~9のいずれか一項に記載の方法。
- プログラム命令を含むコンピュータプログラムを有する非一時的なコンピュータ可読記憶媒体(412)を含むコンピュータプログラム製品であって、前記コンピュータプログラムが、プロセッサ(411)にロード可能であり、請求項1~10のいずれか一項に記載のスピン標識化粒子の拡散パラメータを決定するための方法(100)を前記プロセッサ(411)に実行させるように構成された、コンピュータプログラム製品。
- スピン標識化粒子の拡散パラメータを決定するための磁気共鳴システム(300)であって、
磁気共鳴装置(340)と、
前記磁気共鳴装置(340)によって、少なくとも1つの磁場勾配パルスと高周波パルスとの少なくとも1つの列を試料(390)に与えることと、
前記磁気共鳴装置(340)によって、複数の異なる勾配持続時間、勾配強度、及び勾配方向で前記試料(390)から磁気共鳴測定データを収集することと、
を行うように前記磁気共鳴装置(340)を制御するように構成されたコンピュータ(350)と、
を備え、
前記コンピュータ(350)が、
前記収集された磁気共鳴測定データに基づいて、少なくとも1つの拡散パラメータごとに、少なくとも1つのフーリエ変換を再構成することと、
前記対応する再構成した少なくとも1つのフーリエ変換に基づいて、前記少なくとも1つの拡散パラメータを決定することと、
を行うように構成され、
各勾配パルス列が、少なくとも3つの勾配パルスを含み、少なくとも1つの勾配パルスが、位置に基づいて、前記試料(390)に含まれるスピン標識化粒子に位相シフトを導入するように構成された、システム。 - 前記コンピュータ(350)が、拡散測定結果に基づいて、前記磁気共鳴装置を校正するように構成された、請求項12に記載のシステム。
- 前記磁気共鳴装置が、磁気共鳴分光法、局所磁気共鳴分光法、及び/又は磁気共鳴イメージングを利用して、測定データを取得するように構成された、請求項12又は13に記載のシステム。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US202063001266P | 2020-03-28 | 2020-03-28 | |
US63/001,266 | 2020-03-28 | ||
PCT/SE2021/050273 WO2021201753A1 (en) | 2020-03-28 | 2021-03-26 | A magnetic resonance method, software product, and system for determining a diffusion propagator or related diffusion parameters for spin-labelled particles |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2023519217A true JP2023519217A (ja) | 2023-05-10 |
Family
ID=77929295
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2022556595A Pending JP2023519217A (ja) | 2020-03-28 | 2021-03-26 | スピン標識化粒子の拡散伝搬関数又は関連する拡散パラメータを決定するための磁気共鳴法、ソフトウェア製品、及びシステム |
Country Status (5)
Country | Link |
---|---|
US (1) | US11789106B2 (ja) |
EP (1) | EP4127758A4 (ja) |
JP (1) | JP2023519217A (ja) |
CN (1) | CN115769097A (ja) |
WO (1) | WO2021201753A1 (ja) |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006091983A2 (en) * | 2005-02-25 | 2006-08-31 | University Of Florida Research Foundation, Inc. | Variable diffusion-time magnetic resonance-based system and method |
WO2008109783A2 (en) * | 2007-03-06 | 2008-09-12 | The Regents Of The University Of California | Detecting spin perturbations using magnetic resonance imaging |
US9603546B2 (en) * | 2009-01-26 | 2017-03-28 | The United States Of America, As Represented By The Secretary, Department Of Health And Human Services | Phantom for diffusion MRI imaging |
US10031204B2 (en) * | 2012-06-29 | 2018-07-24 | Cr Development Ab | Quantification of the relative amount of water in the tissue microcapillary network |
DE102012217997A1 (de) * | 2012-10-02 | 2014-04-03 | Siemens Aktiengesellschaft | Verfahren zur zeitabhängigen Intensitätskorrektur von diffusionsgewichteten MR-Bildern und Magnetresonanz-Anlage |
EP3146352A1 (en) * | 2014-04-18 | 2017-03-29 | Koninklijke Philips N.V. | Mri involving the acquisition of an angiography weighted image and of a perfusion weighted image |
WO2016132176A1 (en) * | 2015-02-19 | 2016-08-25 | Synaptive Medical (Barbados) Inc. | Systems and methods for measuring global glymphatic flow using magnetic resonance imaging |
US10959642B2 (en) * | 2016-08-19 | 2021-03-30 | University Of Utah Research Foundation | Methods and systems of evaluating axonal loss and demyelination |
USRE49978E1 (en) * | 2016-11-09 | 2024-05-21 | The Brigham And Women's Hospital, Inc. | Method of performing diffusion weighted magnetic resonance measurements on a sample |
EP3493151A1 (en) * | 2017-11-29 | 2019-06-05 | Koninklijke Philips N.V. | Combination of temporally resolved angiographic images with a spatially resolved angiographic image |
-
2021
- 2021-03-26 WO PCT/SE2021/050273 patent/WO2021201753A1/en active Application Filing
- 2021-03-26 CN CN202180024635.3A patent/CN115769097A/zh active Pending
- 2021-03-26 US US17/913,619 patent/US11789106B2/en active Active
- 2021-03-26 JP JP2022556595A patent/JP2023519217A/ja active Pending
- 2021-03-26 EP EP21778992.4A patent/EP4127758A4/en active Pending
Also Published As
Publication number | Publication date |
---|---|
CN115769097A (zh) | 2023-03-07 |
US11789106B2 (en) | 2023-10-17 |
EP4127758A1 (en) | 2023-02-08 |
WO2021201753A1 (en) | 2021-10-07 |
EP4127758A4 (en) | 2024-04-17 |
US20230124954A1 (en) | 2023-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Nilsson et al. | The role of tissue microstructure and water exchange in biophysical modelling of diffusion in white matter | |
Benjamini et al. | White matter microstructure from nonparametric axon diameter distribution mapping | |
Ianuş et al. | Accurate estimation of microscopic diffusion anisotropy and its time dependence in the mouse brain | |
US10649057B2 (en) | Analysis for quantifying microscopic diffusion anisotropy | |
US7643863B2 (en) | Diffusion tensor and q-space MRI specimen characterization | |
Portnoy et al. | Oscillating and pulsed gradient diffusion magnetic resonance microscopy over an extended b‐value range: implications for the characterization of tissue microstructure | |
US8380280B2 (en) | Non-invasive in vivo MRI axon diameter measurement methods | |
US20100033182A1 (en) | Magnetic resonance specimen evaluation using multiple pulsed field gradient sequences | |
Özarslan et al. | Temporal scaling characteristics of diffusion as a new MRI contrast: findings in rat hippocampus | |
US20210270921A1 (en) | Mri tractography based transit time determination for nerve fibers | |
JP7057355B2 (ja) | サンプルに拡散強調磁気共鳴測定を実行する方法 | |
Tax | Estimating chemical and microstructural heterogeneity by correlating relaxation and diffusion | |
US10871539B2 (en) | Determination of a joint probability distribution of radius and length of anisotropic pores from double pulsed field gradient MRI data | |
Chakwizira et al. | Diffusion MRI with free gradient waveforms on a high-performance gradient system: Probing restriction and exchange in the human brain | |
WO2022020102A1 (en) | Time efficient multi-pulsed field gradient (mpfg) mri without concomitant gradient field artifacts | |
WO2006091983A2 (en) | Variable diffusion-time magnetic resonance-based system and method | |
JP2023519217A (ja) | スピン標識化粒子の拡散伝搬関数又は関連する拡散パラメータを決定するための磁気共鳴法、ソフトウェア製品、及びシステム | |
Fick et al. | Diffusion MRI anisotropy: modeling, analysis and interpretation | |
Kaden et al. | Can T 2-spectroscopy resolve submicrometer axon diameters? | |
US11874360B2 (en) | Method and magnetic resonance apparatus for quantitative, highly reproducible tissue differentiation | |
US11768264B2 (en) | System and method for magnetic resonance fingerprinting with relaxation and diffusion data acquisition | |
Assaf et al. | Diffusion as a probe of tissue microstructure | |
Savickas | Advanced diffusion MRI for microstructure imaging: theoretical developments | |
Pieciak et al. | Q-space quantitative diffusion MRI measures using a stretched-exponential representation | |
Ozarslan et al. | Characterization of anomalous diffusion from MR signal may be a new probe to tissue microstructure |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20240321 |
|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20240423 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20240423 |