JP2020531076A - パルス磁気粒子イメージングシステム及び方法 - Google Patents
パルス磁気粒子イメージングシステム及び方法 Download PDFInfo
- Publication number
- JP2020531076A JP2020531076A JP2020508019A JP2020508019A JP2020531076A JP 2020531076 A JP2020531076 A JP 2020531076A JP 2020508019 A JP2020508019 A JP 2020508019A JP 2020508019 A JP2020508019 A JP 2020508019A JP 2020531076 A JP2020531076 A JP 2020531076A
- Authority
- JP
- Japan
- Prior art keywords
- pulse
- excitation
- ffr
- magnetic
- magnetic field
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 168
- 239000006249 magnetic particle Substances 0.000 title claims abstract description 69
- 238000000034 method Methods 0.000 title claims description 152
- 230000005291 magnetic effect Effects 0.000 claims abstract description 362
- 230000005284 excitation Effects 0.000 claims abstract description 327
- 239000000700 radioactive tracer Substances 0.000 claims abstract description 253
- 238000001514 detection method Methods 0.000 claims abstract description 92
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 claims abstract description 69
- 239000002122 magnetic nanoparticle Substances 0.000 claims abstract description 54
- 238000009826 distribution Methods 0.000 claims abstract description 50
- 230000001052 transient effect Effects 0.000 claims abstract description 39
- 230000005415 magnetization Effects 0.000 claims description 110
- 230000008569 process Effects 0.000 claims description 68
- 238000012545 processing Methods 0.000 claims description 43
- 238000002360 preparation method Methods 0.000 claims description 26
- 230000008859 change Effects 0.000 claims description 25
- 239000002245 particle Substances 0.000 claims description 19
- 238000009877 rendering Methods 0.000 claims description 18
- 238000005259 measurement Methods 0.000 claims description 13
- 239000000696 magnetic material Substances 0.000 claims description 5
- 238000011109 contamination Methods 0.000 claims description 4
- 230000004069 differentiation Effects 0.000 claims description 3
- 230000002093 peripheral effect Effects 0.000 claims description 2
- 230000005389 magnetism Effects 0.000 claims 1
- 238000010586 diagram Methods 0.000 abstract description 10
- 230000004044 response Effects 0.000 description 74
- 239000000523 sample Substances 0.000 description 71
- 230000006870 function Effects 0.000 description 38
- 239000011159 matrix material Substances 0.000 description 27
- 238000013459 approach Methods 0.000 description 24
- 230000000694 effects Effects 0.000 description 22
- 230000007704 transition Effects 0.000 description 22
- 230000000116 mitigating effect Effects 0.000 description 21
- 238000005070 sampling Methods 0.000 description 21
- 238000013461 design Methods 0.000 description 18
- 230000008929 regeneration Effects 0.000 description 17
- 238000011069 regeneration method Methods 0.000 description 17
- 230000001976 improved effect Effects 0.000 description 14
- 230000000630 rising effect Effects 0.000 description 11
- 238000002474 experimental method Methods 0.000 description 10
- 239000000203 mixture Substances 0.000 description 10
- 238000004422 calculation algorithm Methods 0.000 description 9
- 238000001914 filtration Methods 0.000 description 9
- 239000002105 nanoparticle Substances 0.000 description 8
- 230000002829 reductive effect Effects 0.000 description 8
- 230000005540 biological transmission Effects 0.000 description 7
- 230000000737 periodic effect Effects 0.000 description 7
- 230000003595 spectral effect Effects 0.000 description 7
- 230000008901 benefit Effects 0.000 description 6
- 238000000354 decomposition reaction Methods 0.000 description 6
- 230000001939 inductive effect Effects 0.000 description 6
- 230000010354 integration Effects 0.000 description 6
- 230000033001 locomotion Effects 0.000 description 6
- 230000014759 maintenance of location Effects 0.000 description 6
- 230000035945 sensitivity Effects 0.000 description 6
- 238000005284 basis set Methods 0.000 description 5
- 239000003990 capacitor Substances 0.000 description 5
- 238000009472 formulation Methods 0.000 description 5
- 230000006872 improvement Effects 0.000 description 5
- 230000006698 induction Effects 0.000 description 5
- 230000000670 limiting effect Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000007423 decrease Effects 0.000 description 4
- 230000001965 increasing effect Effects 0.000 description 4
- 238000002347 injection Methods 0.000 description 4
- 239000007924 injection Substances 0.000 description 4
- 230000003993 interaction Effects 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 238000011002 quantification Methods 0.000 description 4
- 229920006395 saturated elastomer Polymers 0.000 description 4
- 230000000638 stimulation Effects 0.000 description 4
- 230000002123 temporal effect Effects 0.000 description 4
- 238000013519 translation Methods 0.000 description 4
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 3
- 230000003321 amplification Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 229910052802 copper Inorganic materials 0.000 description 3
- 239000010949 copper Substances 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000001419 dependent effect Effects 0.000 description 3
- 238000009795 derivation Methods 0.000 description 3
- 230000004907 flux Effects 0.000 description 3
- WTFXARWRTYJXII-UHFFFAOYSA-N iron(2+);iron(3+);oxygen(2-) Chemical compound [O-2].[O-2].[O-2].[O-2].[Fe+2].[Fe+3].[Fe+3] WTFXARWRTYJXII-UHFFFAOYSA-N 0.000 description 3
- 239000003446 ligand Substances 0.000 description 3
- 238000013507 mapping Methods 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 238000003199 nucleic acid amplification method Methods 0.000 description 3
- 230000010412 perfusion Effects 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 238000000926 separation method Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 230000003851 biochemical process Effects 0.000 description 2
- 210000004027 cell Anatomy 0.000 description 2
- 230000003727 cerebral blood flow Effects 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000001143 conditioned effect Effects 0.000 description 2
- 239000012809 cooling fluid Substances 0.000 description 2
- 230000003111 delayed effect Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 230000012010 growth Effects 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 229910052751 metal Inorganic materials 0.000 description 2
- 230000036961 partial effect Effects 0.000 description 2
- 230000010287 polarization Effects 0.000 description 2
- 229920013655 poly(bisphenol-A sulfone) Polymers 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 239000003642 reactive oxygen metabolite Substances 0.000 description 2
- 230000001172 regenerating effect Effects 0.000 description 2
- 230000003252 repetitive effect Effects 0.000 description 2
- 230000009885 systemic effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- OJQSISYVGFJJBY-UHFFFAOYSA-N 1-(4-isocyanatophenyl)pyrrole-2,5-dione Chemical compound C1=CC(N=C=O)=CC=C1N1C(=O)C=CC1=O OJQSISYVGFJJBY-UHFFFAOYSA-N 0.000 description 1
- 229910000881 Cu alloy Inorganic materials 0.000 description 1
- 101001099381 Homo sapiens Peroxisomal biogenesis factor 19 Proteins 0.000 description 1
- 206010061218 Inflammation Diseases 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 102100038883 Peroxisomal biogenesis factor 19 Human genes 0.000 description 1
- 208000010378 Pulmonary Embolism Diseases 0.000 description 1
- 208000006011 Stroke Diseases 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000013475 authorization Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 239000012472 biological sample Substances 0.000 description 1
- 230000000740 bleeding effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 230000002490 cerebral effect Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000003750 conditioning effect Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000005292 diamagnetic effect Effects 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 208000035475 disorder Diseases 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000003631 expected effect Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000005294 ferromagnetic effect Effects 0.000 description 1
- NBVXSUQYWXRMNV-UHFFFAOYSA-N fluoromethane Chemical compound FC NBVXSUQYWXRMNV-UHFFFAOYSA-N 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000003306 harvesting Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 208000015181 infectious disease Diseases 0.000 description 1
- 230000004054 inflammatory process Effects 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000007620 mathematical function Methods 0.000 description 1
- 150000002739 metals Chemical class 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- RVZRBWKZFJCCIB-UHFFFAOYSA-N perfluorotributylamine Chemical compound FC(F)(F)C(F)(F)C(F)(F)C(F)(F)N(C(F)(F)C(F)(F)C(F)(F)C(F)(F)F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)F RVZRBWKZFJCCIB-UHFFFAOYSA-N 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 229920000889 poly(m-phenylene isophthalamide) Polymers 0.000 description 1
- 230000002685 pulmonary effect Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000011270 sentinel node biopsy Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 210000000130 stem cell Anatomy 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000012385 systemic delivery Methods 0.000 description 1
- 238000009423 ventilation Methods 0.000 description 1
- 230000035899 viability Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/12—Measuring magnetic properties of articles or specimens of solids or fluids
- G01R33/1276—Measuring magnetic properties of articles or specimens of solids or fluids of magnetic particles, e.g. imaging of magnetic nanoparticles
-
- 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/0515—Magnetic particle imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux
- G01R33/10—Plotting field distribution ; Measuring field distribution
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Nanotechnology (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
Abstract
Description
本出願は、米国仮出願第62/546,395号(2017年8月16日に出願)に対する優先権を主張する。なおこの文献の全体の内容は、参照により本明細書に組み込まれている。
パルス波形を用いてMPI信号をエンコードすることが、正弦波的な連続波形を用いた場合には一般に実現できない方法で可能である。たとえば、トレーサー信号から送信フィードスルー干渉を時間的に分離することができる。またパルス波形によって、ナノ粒子トレーサー分布の定常状態磁気物理過程の測定を可能にすることができ、トレーサーに付随する磁気応答時間が有限であるために生じる連続的な正弦波MPIにおけるぼかし効果を潜在的に回避することができる。同時に、パルス波形によって、これらの磁気応答時間を定量化できる方法及びエンコーディング領域を得ることができ、画像コントラストの増加及び付加情報を測定する可能性をもたらすことができる。たとえば、トレーサー密度画像に加えて、トレーサー緩和物理過程の測定を定量化して空間分解して報告する緩和画像も生成され得る。他の実施形態では、完全な4次元イメージングデータセットが与えられ得る。このデータは、同じ励起パルスシーケンスからまたは異なる励起パルスシーケンスの連続印加を介して同時に取得され得る。
図1Bに示すように、MPIにおける基本的なパルスエンコーディングに対する要件には、無磁場領域(FFR)として知られる空間的に変化する磁場パターンまたは構造、時間的に制御可能で空間的に均一な励起磁場源、及び1つ以上の磁場パルスが含まれ、各磁場パルスには、励起磁場源によって生成される1つ以上の実質的に一定の磁場波形が含まれる。いくつかの実施形態では、較正する場合、または空間的に局在化された画像を生成する必要がない検知応用例では、FFRを省略してもよい。いくつかの実施形態では、パルス励起磁場に加えて1つ以上の空間的に均一なバイアシング磁場を印加して、さらなる側面のデータ(たとえば、緩和データ)を、FFRを用いてまたは用いずに取得してもよい。
磁気飽和
MPIは、磁気的に飽和可能な磁気トレーサーをイメージングする。印加磁場とナノ粒子磁化との間の関係は磁化曲線(M−H曲線として知られる)によって記述することができる。磁化曲線は、印加磁場に対するナノ粒子集合の磁化応答を与える(図7を参照)。印加磁場が小さいときには、磁気ナノ粒子はその磁化を急速に変化させる。印加磁場が大きくなると、M−H曲線は磁化のプラトーまたは一定値に漸近的に近づく。そして、全印加磁場の大きさが大きいところでは、これらのトレーサーは、印加磁場をさらに増加させてもその磁化を実質的に変化させず、そのため粒子は「飽和した」と言われる。超常磁性酸化鉄(SPIO)トレーサーの場合、印加磁場に対するSPIOの磁化応答はランジュバン曲線に従うと理解されている。
典型的なMPIスキャニングシステムには、特定の空間的な磁場パターン(本明細書ではFFRと言う)を発生させるために用いられる1つ以上のアクティブまたはパッシブな磁場源が含まれている。FFRは、印加磁場が磁気飽和値を下回る空間領域として規定される。これは、印加磁場が磁気飽和値を上回る近くの空間領域とは異なる状態とすることができる。FFRには多くの場合、必ずしもそうではないが、線形磁場勾配が含まれているため、磁場がFFRと近くの飽和磁場領域との間で滑らかに移行する。FFRは、無磁場点、無磁場線、またはマクスウェルの方程式によって課される制限を受ける図2に示したようなより一般的な形状であってもよい。
典型的なMPIシステムでは、FFRを形成する勾配磁場に重ね合わせて、空間的に均一だが時間変化する磁場を印加することによって、トレーサー信号を発生させる。このように印加された均一磁場は一般的に、FFRの形状を維持しながら、FFRを空間及び時間においてシフトまたは平行移動させるように作用する。いくつかの実施形態では、たとえば、FFL内の線と同軸の印加磁場による均一な励起の場合、均一磁場を付加してもFFRの平均場所はシフトせず、その代わりに、たとえば、磁気的バイアスを加えてFFRを低磁場領域(LFR)に変換する。
典型的なMPI軌跡には、正弦波的に変わり連続的に送信する励起磁場(たとえば、45kHzでの正弦曲線)が含まれるが、本発明のいくつかの実施形態では、非正弦波スキャニング軌跡を伴う励起波形が提供される。本明細書ではこれらの波形をパルス波形と言う。これらのパルス励起波形の重要な特徴は、磁場の大きさがある時間の間実質的に一定に保持されるいくつかの周期的または非周期的成分が含まれることである。実質的に一定の意味は、磁場の大きさが、何らかの最初の過渡的または立ち上がり時間の後及び/または以後の過渡的または立ち上がり時間の前に、何らかの望ましいまたは理想的な磁場の大きさから、規定されたエラーを超えて外れることはないということである。エラーの大きさ例としては、これらに限定されないが、望ましい値の10%、望ましい値の5%、及び望ましい値の1%が挙げられる。この実質的に一定の値に磁場が維持される時間は500ナノ秒〜500ミリ秒の範囲であってもよいが、この範囲に限定されない。図4に、種々の振幅、極性、及び継続時間(t1で示す)を伴う実質的に一定の波形成分を例示する。パルスシーケンス全体において、図4に示すような種々の振幅、極性、及びt1継続時間を伴う多くの実質的に一定の波形成分を用いてもよい。
励起波形における実質的に一定の値を、種々の基本形状のパルス波形成分によって分離してもよい。概略的に、これらの基本形状を正確なタイミングで互いにつないで、より複雑なパルス波形を形成してもよい。このようなパルス波形基本形状の1つの重要なタイプは、立ち上がり及び/または立ち下がりエッジを通って急速に移行するものである。このような急速に磁場が移行する間、時間変化する磁場波形は任意の形状であってもよい。たとえば、波形は、図5及び図6に例示するように、直線的に、指数関数的に、正弦波的に、またはその重ね合わせで立ち上がってもよい。図5に示すのは、立ち上がり及び/または立ち下がりエッジを介した急速な移行を例示する典型的な励起波形であり、最大の継続時間t2を伴う成分で示す。この継続時間t2は100ナノ秒〜10マイクロセカンドの範囲であってもよい。また図5には、より大きいパルスMPI波形またはパルスシーケンスを形成するために用いてもよい他の典型的な過渡パルス及びパルス状波形成分も示す。
MPIで用いる磁気トレーサーの特徴は、その定常状態応答(M−H曲線(図7)で示されることが多い)と印加磁場に対してその時間変化する応答との両方である。発明者らは、MPIにおいて、磁化曲線の印加磁場に関する微分が、理想的な、ネイティブの、または定常状態の点広がり関数(PSF)に比例していることを、理論的及び実験的に見いだした。この定常状態のナノ粒子応答によって、逆畳み込みなどのツールを用いなければ、MPIにおいて可能な分解能に対する限定が生じるということが、従来理解されていた。実際には、トレーサーの動的な磁気緩和挙動があるために、このような分解能を標準的な正弦波MPIにおいて実現するのは難しい可能性がある。さらに、本明細書で説明するように、パルスMPI法によって、これまで考えられていたこの限定を超えて、定常状態ランジュバン理論によって予測されるものより良好な分解能を実現することが、逆畳み込みを行わずに、緩和ダイナミクスを利用することによって可能となる。
波形基本形状(たとえば、実質的に一定である過渡パルス及び時間)を組み合わせて、より複雑なパルス波形及びMPI励起パルスシーケンスを形成してもよい。MPIパルスシーケンスの2つの基本成分には、1つ以上の別個の磁化準備基本形状と1つ以上の別個の信号読み出し時間とが含まれていてもよい。これについて、図8及び図9に概念的に例示する。いくつかの実施形態では、一連の磁化準備パルスは継続時間が以後のより長い読み出し時間より短くてもよい。いくつかの実施形態では、準備及び読み出し時間の相対的な継続時間が連続シーケンスにおいて変化してもよい。いくつかの実施形態では、実質的に一定の時間の間の移行が、励起及び/または準備として機能する。
磁場が実質的に一定に保持される1つ以上の時間を、磁化準備の間及び信号読み出し時間の間の両方で用いてもよい。磁場が実質的に一定に保持される時間の間、トレーサー分布の磁化は新しい磁化定常状態に向かって進展する。この定常状態は、空間的に変化する可能性があり(なぜならば、たとえば、FFR及びそれが形成する空間的な磁場パターンと磁場依存のトレーサー緩和特性との間の相互作用が存在するため)、またトレーサーを支配する磁気緩和物理過程によって進展する。概略的に、磁化の進展はほぼステップ応答のようであるため、誘導的に受信された信号はインパルスのような応答を記録する。図11に、勾配が存在しない状態でサンプルに印加されるパルス励起波形を形成することができるテーブルトップリラクソメータシステムを用いて取られた典型的なデータを示す。本明細書では、このようなデバイスを任意波形リラクソメータ(AWR)という。図11では、方形波状の励起パルスシーケンスによって、誘導受信部コイルに減衰インパルス応答信号が発生する。これは、磁化の成長及び減衰に対応する。なぜならば、サンプル分布の正味の磁化は、方形波の各実質的な一定値に対応する全印加磁場状態に付随する定常状態値に向かって進展するからである。実質的に一定の時間が十分に長い継続時間ならば、定常状態が実現され、トレーサー分布の定常状態の2つの別個な値が、トレーサーに対応付けられる定常状態M−Hまたはランジュバン曲線に従ってサンプリングされている。図57に、可変の実質的に一定の保持時間(または方形波半周期時間)を伴う方形波パルスMPIを用いた実験データを示す。保持時間がトレーサー磁気緩和プロセスに対して短すぎると、トレーサー磁化は、次の過渡パルスの前に定常状態を実現することができない。この結果、図57に示すように、信号強度及び分解能の両方が低下する。保持時間が、長さが不十分の保持時間から長くなると、定常状態に対応付けられる最大信号及び最適な分解能への漸近アプローチを観察することができる。
勾配磁場またはFFR構造がパルス波形の実質的に一定の時間の間に存在する場合、これは、空間的に異なる緩和時間に起因するトレーサー分布に関係する空間情報の第2のエンコーディングとなる。図12に、この空間情報が、空間的に変化する磁場(たとえば、FFR)の存在下で、任意の実質的に一定の時間に対応付けられる受信信号において、時間的にどのようにエンコードされるかを例示する。実質的に一定の時間に十分に長い(たとえば、観察される最長の磁気緩和現象よりもはるかに長い)継続時間が適用された場合、トレーサー分布全体は、ある定常状態分布から別の定常状態分布へ進展する。しかし、実質的に一定の値の正確な継続時間及び定常状態が実現されたか否かには関係なく、時間の間の時間ドメイン信号によって空間情報がエンコードされる。
読み出し時間の前に、または可能性として読み出し時間と同時に、1つ以上のパルス成分によって、トレーサー分布の磁化状態を準備するか、セットアップするか、そうでなければ動的に成形してもよい。図9及び図10に、このような準備時間で用いてもよいいくつかの典型的なパルス及びパルスを連続して集めたものを示す。これらのパルスは、トレーサー分布内の空間依存する磁気緩和ダイナミクスを利用するのに十分に短い継続時間であってもよい。たとえば、準備時間における一連のパルスは、非飽和領域内のどこででも定常状態磁化が達成されるのを防ぐのに十分に短い継続時間であってもよい。しかし、このパルシング時間における磁化の変化または成長の度合いは、FFR構造の正確な軌跡に空間的に依存する。概略的に、磁気緩和時間がFFRからの距離とともに単調減少する(緩和時間と全印加磁場強度との間の単調な関係に起因する)ために、緩和ダイナミクスはFFRアイソセンタに近いほど最も遅く、FFRアイソセンタから遠いほど最も速い。このように、読み出し時間における信号を、準備時間における分布を事前に操作することによって予め加重することができる。
パルス波形の概略的な特徴は、たとえば、磁化準備時間の後に別個の信号読み出し時間を規定できることである。いくつかの実施形態では、読み出し時間は、トレーサー分布の磁化の進展が定常状態位置に近づくことが観察される実質的に一定の時間から構成される。いくつかの実施形態では、この読み出しは、急速に移行する立ち上がり及び/または立ち下がりエッジに即時に続いてもよい。いくつかの実施形態では、実質的に一定の読み出し時間の継続時間を、対象とするサンプル内でのトレーサー分布の磁気的定常状態の実現が確実になるように故意にデザインする。他の実施形態では、実質的に一定の読み出し時間を、定常状態条件が実現されないように故意にデザインする。いくつかの実施形態では、読み出し時間には非定常成分が含まれていてもよい。図9、図10、図11、図12、及び図13に、パルスMPI励起波形における典型的な読み出し時間を示す。
パルス波形を用いてMPIにおける均一な励起コイルを駆動してもよい。これらのパルス波形はFFRの状態及び/または場所に影響し、磁気トレーサー分布が存在する場合に信号を発生させる。
いくつかの実施形態では、パルスの効果がFFR構造の場所を平行移動させることとなるように、パルス波形を均質コイルによって印加する。これらの波形の実質的に一定の時間の間、FOV内でFFR構造は静止している。図14に、2つの励起コイルの主軸によって規定される平面内のどこへでもFFRの場所を平行移動させるように、直交する均質コイル内のパルスがどのように作用するかを例示する。FFR構造のアイソセンタが3つのすべての寸法において局在化している場合(たとえば、FFP)、均一なコイルを用いてFFR構造を任意の方向に平行移動させることができる。FFRアイソセンタの連続的な場所を接続する曲線(たとえば、図14の線1、2、及び3)に沿って位置する任意のトレーサーは、パルスによって励起される。軌跡は単純な線より一般的な形状であってもよく、トレーサー信号を支配する磁気緩和は、この軌跡と、FFRの場所とトレーサー分布との間の瞬間的な相対ジオメトリとに概ね依存する。単純な方形状パルスの場合、各ステップ状の移行によって、FFRアイソセンタの最初の場所とステップ後の場所とを接続するx空間内の線が規定される。パルスの時間に対して、最初のFFR場所に近い励起線の付近に配置されたトレーサーの応答は、最終的なFFR場所に近い励起線の付近に配置されたトレーサーよりも速い。
いくつかの実施形態では、励起パルスを無磁場線(FFL)の線に沿った方向に印加してもよい。図15に例示するように、線方向にパルスが印加される間、勾配磁場構造の状態は変更されるが、構造の場所は線に直交する平面内で変化しない。励起パルスを印加する直前の磁場の値に応じて、FFLは、FFLから低磁場線(LFL)へ、LFLからFFLへ、またはある方向に分極されたLFLから反対方向に分極されたLFLへと移行する。LFLを、同じかまたは類似の空間的に変化する構造によって特徴付けてもよいが、ゼロ磁場大きさを通って横断する点(複数可)がないように大きさをバイアスしてもよい。FFLの場合と同様に、磁場のベクトル方向は、LFLに直交する平面の一方の側から他方へ回転する。この励起領域では、FFLまたはLFLを規定する線に沿ったすべてのトレーサーは同様に励起されて、全印加磁場への依存性の点で同じ緩和物理過程を受ける。これは、線に沿っているが直交面内での投影励起を伴わない完全なFOV投影エンコーディングを表す。線に沿った励起パルスを、線に直交する平面内の励起パルスと平行または直列に行ってもよい。
いくつかの実施形態では、FFLの方向にある均一なバイアス磁場を通常、たとえば、直交面内で、線に沿って、または両方で別個に印加される別個の励起パルスとともに印加してもよい。このようなバイアス磁場によって、FFLの代わりにLFLが設定される。図15に、バイアシング波形を伴うLFLを例示する。バイアシング磁場の付加は、LFLが直交面内で動いたときにトレーサー集合体の回転の対称性を壊す働きをしてもよい。またトレーサーは任意の点においてゼロ磁場領域によって励起されることはないため、より急速な緩和ダイナミクスが通常、よりSNR効率的な信号エンコーディングに対して誘起されてもよい。
いくつかの実施形態では、勾配磁場または空間的に変化するFFR構造が存在しない状態で、サンプルに均一パルス励起波形を印加する。この領域では、線形バイアシング磁場をゆっくりと立ち上げてもよく、直線勾配磁場が用いられたときにスキャニングシステム内の点状源の経験に直線的に比例する印加磁場大きさ空間での軌跡を通してサンプル全体が取られる。このようなシステムを用いて、MPIパルスシーケンスの効果を迅速にテストすることができ、また集合体内のトレーサーサンプルまたはサンプル物体に対応付けられる1次元(1D)のPSFを測定することができる。なお、このようなアプローチによって、すべてのトレーサーが、場所とは関係なく同様に励起される。こうして、1DのPSFによって、観察体積内のすべてのトレーサーの集合体または平均挙動が記述される。このタイプのスキャンは、FFRを発生させる磁気源をターンオフできるならば、小さいテーブルトップリラクソメータシステム(たとえば、任意波形リラクソメータ(AWR))を用いて、またはフルイメージングスキャナを用いて、行うことができる。
パルスMPIにおいて利用可能なエンコーディング機械(種々のパルスタイプ、概念上のパルス波形成分、及びパルス波形によって駆動されるコイル配置を含む)を組み合わせて、完全なパルスMPIパルスシーケンスを形成し、イメージングFOVを十分にサンプリングすることができる。図17及び55に、典型的なパルスMPIパルスシーケンス図を例示する。典型的なパルス励起波形が他のシステムシーケンスと関連付けて示されている。これらの実施形態では、1つ以上のパルスシーケンスを用いて、2Dスライスを、FFLスキャニングシステムを用いて十分にサンプリングする。サンプルとFFLとの間で回転するジオメトリを伴う複数のスライスをサンプリングして、投影再構成技術を用いたトモグラフィ3D画像形成を可能にする。比較的速いパルス励起の他に、もっと遅いシフティング成分によって、2DFOV全体のサンプリングが可能になる。これを図18に示す典型的なラスタースキャニング軌跡に例示する。いくつかの実施形態では、1つ以上の磁化準備成分と1つ以上の別個の読み出し時間とを含む励起パルス波形を何度も繰り返してもよい。各繰り返しは、もっとゆっくりと変化するシフト波形によって規定されるFFR構造の別個の平均場所に対応付けられる。概略的に、遅いシフト成分に対応付けられるはるかに遅いスルーレートは、励起波形における実質的に一定の磁場時間に対して課された実質的に一定であるという基準(たとえば、望ましい値の1%または10%など)に違反しない。いくつかの実施形態では、この違反は許容してもよいし望ましくてもよい。
パルス波形を伴うエンコーディングMPI信号の利点の多くを例示するパルスMPIパルスシーケンスの1つは、方形波状励起波形を用いることである。これらの波形には、立ち上がり及び立ち下がりエッジを介した急速な移行によって分離された磁場の実質的に一定の値の交互に変わる時間が含まれる。図4、図5、及び図6にこれらの成分を示す。用語方形波状は、急速な移行及び実質的に一定の値の時間が前述の特性を有する限り、繰り返し励起パルス成分の全体形状が、方形に、台形に、またはより複雑なもの(たとえば、指数関数的または正弦波的に変化する成分など)に類似してもよいことを強調するために用いられる。このように規定された方形波状励起パルスシーケンスは、最小の磁化準備時間を有すると考えられ得る。前の実質的に一定の読み出し時間で磁化を特定の状態で準備し、次の実質的に一定の読み出し時間への急速な移行で、勢いを与える信号を与え、第2の実質的に一定の読み出し時間で、読み出し中に磁化分布が進展する先の定常状態ターゲットが設定される。
定常状態再生シーケンスによって、他のエンコーディング領域で、または応答時間が大変異なる多重応答ダイナミクスが全体磁化進展(したがって受信信号)に寄与する場合、はっきりしないトレーサー応答ダイナミクスをエンコードすることができる。
より複雑なエンコーディング領域及びパルスシーケンスが、たとえば、より大きなSNR効率を得るためには望ましい場合がある。図21に、直交配置で駆動された直交コイルに一連のパルスが印加されると、FFRアイソセンタがFOV内の特定の場所の周りに、特定の励起パルス列の間にどのように回転するかを例示する。説明を目的として、図示するように、取り囲むパターンの半径またはサイズは、FFRジオメトリ、イメージングFOV、及びスローシフトラスターパターンに対して必ずしも一定の比率ではない。たとえば、取り囲むパターンは実際には比較的はるかに小さくてもよい。さらに、各場所で印加されるパルスは単純な実質的に一定の時間でなくてもよい。図5に例示するように、多くのパルシング形状の何れかが急速な移行によって分離されてもよい(たとえば、平均FFRの場所に急速に移行した後の短時間の正弦波パルス)。中心点を囲む番号付けされた場所のそれぞれに対応付けられるパルスの継続時間は、定常状態条件を設定するのに必要なものよりはるかに短くてもよいが、急速な移行及び/または通常の励起に対する応答に対応付けられる著しい時間ピークを引き出す。いくつかの実施形態では、定常状態を誘起する実質的に一定の時間を、このような速い励起及び読み出し波形と組み合わせて用いてもよい。概略的に、回転励起があると、短い時間にわたってイメージングFOV内の場所または小さい領域に対応付けられる比較的大きな総合信号が得られる場合があり、より単純なラスタースキャンと比べたときにSNR効率が増加する。この取り囲むパルス波形をスローシフティング磁場と共に繰り返して、時間とともに全体イメージングFOVをサンプリングすることができる。
図17、図18、及び図19に例示するように、パルス励起波形を通常、他の重要なシーケンスとともに用いて、完全なトモグラフィまたは投影イメージングを行う。概略的に、pMPIエンコーディングをシステム行列法とともに用いてもよいことに注意されたい。システム行列法では、多くの較正スキャンを、イメージング視野(FOV)内の異なる場所に配置された点状源を用いて行って、システム行列を構成する。システム行列を反転させて、可能性として同じシステム及びスキャニング軌跡を用いて取られるもっと複雑なサンプルの画像が生成され得る。
大きなイメージング視野(FOV)(たとえば、励起波形に対応付けられるズレよりも大きい)をカバーするために、ゆっくりと変化する均一波形を繰り返し励起波形と重ねてもよい。この結果、図17、18、及び55に示すように、FFRアイソセンタの平均場所が時間とともにゆっくりと平行移動することができる。スローシフトティング波形は、周期的な励起波形の個々の励起成分の時間スケールと比べて十分に遅くてもよく、個々の励起時間の間にFFR構造を移動させることに対するスローシフトティング波形の影響は無視することができ、またたとえば、パルス波形の実質的に一定の時間における望ましい大きからの最大偏位さに対する時間保証が維持される。他の実施形態として、特に速い取得が望まれる場合には、スローシフトティング波形の効果によって、トレーサーの励起を無視できない仕方で増大させてもよい。
FFR構造を発生させる勾配磁場の供給源として電磁石を用いる場合、図22に示す典型的な勾配波形に例示するように、応用例に応じて異なるタイプの勾配波形が望ましい場合がある。たとえば、一定の勾配強度をイメージングスキャンの全体にわたって印加してもよい。この強度を複数のスキャンの間で変更して、異なる勾配強度において複数のイメージングデータセットを取るようにしてもよい。たとえば、低分解能において取る速いスカウトスキャンを、勾配強度が高くて固有の空間分解能エンコーディングが向上されたより長いスキャンの前に行ってもよい。スカウトスキャンから得られる情報を用いて、以後のスキャンにおけるイメージングパラメータを設定してもよい。
パルスMPIエンコーディング方式では、受信データから画像を形成する複数の方法が認められる。図23に、パルスMPIx空間デジタル信号処理及び再構成における主要なコンポーネントを記載する。受信部コイルから得たアナログ受信信号をサンプリングして、受信エレクトロニクスチェーンの終わりにデジタル化する。これらの信号を、スキャンにおいて印加される既知のパルスシーケンスとともに用いて、種々の信号処理及び調整タスクを実行する。これらには、タスク、たとえばデジタルフィルタリング、区分的な減少ステップ、たとえば積分または閾値化、及び位相補正が含まれていてもよい。これらのタスクを、データをたとえばフーリエ空間内に変換することなく、時間ドメインにおいて完全に行ってもよい。他の実施形態では、データを好適な基底系上に変換または投影してもよい。
トレーサー密度の画像を、種々のグリッディング方法を用いて再構成してもよい。実質的に一定の読み出し時間またはより複雑な繰り返しパルス励起波形に対応付けられる時間ドメインデータを積分して、結果として得られる値を、ピクセルまたはボクセル場所にグリッディングしてもよい。グリッディング手順では、FFRアイソセンタの既知のx空間軌跡を用いて、時間ドメイン信号を画像空間場所にマッピングしてもよい。いくつかの実施形態では、繰り返しパルス励起の間のFFRアイソセンタの平均場所、または実質的に一定の読み出し時間に対応付けられるFFRのアイソセンタの単一の場所によって、このマッピングに対するルールが得られる。
読み出し時間に対応付けられる時間ドメインデータの積分を、読み出し中のFFRの場所に基づいて、規則的なピクセルまたはボクセルグリッド上に直接内挿またはグリッディングしてもよい。またグリッド場所をFFRの場所によって部分的に決定することを、読み出す前に、たとえば磁化準備または前の読み出し時間の間に、行ってもよい。
図28に例示する一実施形態では、これによって、1つ以上の受信コイルに対する受信電圧信号の曲線の下の積分領域であって、各方形波状半周期ステップに対応付けられる積分領域を、平均FFRアイソセンタ場所にグリッディングする。勾配磁場がない状態でのスキャンの場合またはAWRを用いたときは、このプロセスは、バイアス磁場を用いたときの全印加磁場の平均にグリッディングすることに対応する。図28によって記載される実施形態では、受信信号を積分することから得た単一値を、前の半周期と現在の半周期とにおけるFFRアイソセンタ場所の平均値に、言い換えると、各ステップの前後のFFRの場所を接続する線に沿った中間点にグリッディングする。勾配がない状態でのスキャンの場合、印加磁場空間における対応するグリッド場所は、前及び現在のステップに対応付けられる2つの全印加磁場値の平均である。方形波励起におけるステップジャンプは投影励起を表すため、隣接する方形波半周期を接続する線に沿ったすべてのトレーサーをある程度まで励起すると(この実施形態の場合と同様に直接積分グリッディングアプローチ)、方形波振幅のサイズに基づく投影ボケを被る。またボケは方形波状半周期の継続時間に、各半周期の間に十分な定常状態条件が実現されるまで依存する。定常状態が実現されると、ボケは最大に達して中間点の周りに対称的になり、それ自身を矩形関数として数学的に特徴付けられ得る。
いくつかの実施形態では、直接積分グリッディングの代わりに投影再構成技術を適用してもよい。これは、大きい投影ボケペナルティを被ることなくもっと大きいパルス励起振幅を用いるためには、望ましいアプローチであり得る。パルス励起を励起方向に沿った局所投影と見なしてもよい。直交励起コイル、またはサンプルとパルスの励起方向との間の相対ジオメトリを回転させる他のメカニズムを用いて、異なる励起角度における局所的な部分視野(pFOV)のサンプリングを行ってもよく、図29に例示するような局所投影再構成技術が可能になる。
パルスMPI励起シーケンスとして、各連続的な読み出し時間の終わりまでにトレーサー分布の磁化が定常状態を実現するものを用いるとき、磁化画像を中間または最終画像として直接再構成することを考えることができる。
磁化準備シーケンスを伴ういくつかのパルス波形は本来的に、図10及び図12に例示するような勾配またはFFRの状況で存在する空間的に異なる磁気緩和状態を用いて信号を加重する。任意の磁化準備に加えて、少なくとも1つの実質的に一定の領域を含むすべての読み出し時間は、良好に特徴付けられた固有の空間エンコーディングを、図12に例示するような一定の読み出し時間に対応付けられる時間ドメインデータに含む。
概略的に、パルスMPIにおける受信時間ドメイン信号内で磁気緩和情報を直接エンコードする。この磁気緩和情報を測定し、フィッティングさせ、定量化し、そうでなければ特徴付けてもよい。この定量化は、非イメージング、センサ、またはスペクトルフォーマットの勾配磁場を用いずに行うことがきるか、または勾配磁場の状況で行って緩和画像及び/または4Dイメージングデータセットを形成してもよい。緩和画像を形成するときに、磁気緩和現象の測定値を、FFRアイソセンタの既知の軌跡に基づいて画像空間場所にグリッディングし、トレーサー密度画像とは別個でこれに大きくまたは完全に直交する情報を用いて緩和画像を形成する。いくつかのエンコーディング方式を用いて、緩和マップを、印加磁場の関数として、像空間内の各ピクセルまたはボクセルに対応付けることができる。勾配磁場またはFFR構造が存在しない場合、全サンプル体積を表す印加磁場の関数としての緩和マップが急速に形成され得る。他の実施形態では、空間及び緩和関連変数の関数としてトレーサー密度を記述する4Dデータセットを得てもよい。いくつかの実施形態では、システム行列法を用いて緩和画像を構成することができる。好適なシステム行列を構成するときに、動的パラメータの部分的なx空間グリッディング及びフィッティング、測定、または定量化が望ましい場合がある。
いくつかの実施形態では、各読み出し時間に対応付けられる未処理の時間ドメイン信号によって磁化インパルス応答を直接エンコードする。たとえば、急速な移行の後に実質的に一定の時間が存在する場合、磁気トレーサー分布はステップ応答のような方法で新しい定常状態磁化分布に進展する。誘導受信コイルを用いることによって適用される時間微分によって、未処理の時間ドメイン信号において直接観察できるインパルスのような応答が実現する。これらの状況では、1つ以上の特性緩和パラメータを、各読み出し時間に対して測定するかまたはフィッティングさせて、規則的なピクセルまたはボクセルグリッド上にグリッディングするかまたは内挿して、緩和画像を形成してもよい。いくつかの実施形態では、各ステップに対応付けられる画像ドメイングリッド点は、そのステップ中の既知のFFRアイソセンタ場所である。
図35に、方形波状パルスシーケンスエンコーディングを使用し、図24及び図25記載するように再構成された緩和イメージングの一実施形態を例示する。この実施形態では、方形波励起は、交互に変わる時間ドメインインパルス応答をもたらす。これらのインパルス応答はそれぞれ、FFRアイソセンタの付近のトレーサーに対応付けられる緩和現象の加重重ね合わせの結果である。概略的に、個々の構成する緩和応答は、FFRアイソセンタからの距離及び局所的な微環境状態によって決定され、ランジュバン物理過程に従ってFFRアイソセンタからの距離によって加重される。これに関連して、各インパルス応答を、局所的に集合した(FFRに対して)トレーサー応答と見なしてもよい。各インパルス応答を、たとえば、単一の平均、特性、または総合の指数関数的な緩和時定数を、各方形波ステップまたは半周期に対応付けられる時間ドメインデータにフィッティングすることによって定量化してもよい。各フィッティングされた値を次に、ステップ中のFFRアイソセンタの平均位置、または勾配磁場が存在しない状態での励起における全印加磁場に対応付けてもよい。トレーサー密度イメージングの場合と同様に、結果は、概略的に、サンプル点と最終的なイメージンググリッドとの間の多対1の関係である。内挿グリッディングステップを行って、測定した緩和時間のセットをピクセルまたはボクセルの最終的な規則的グリッドにマッピングしてもよい。緩和ダイナミクスを測定し、フィッティングし、または定量化する前に、前述したように緩和加重またはフィルタリングを適用して、たとえば、緩和フィッティングの測定をさらに局在化してもよい。いくつかの実施形態では、動的な緩和値をフィッティングして、フィッティングされた値をx空間にグリッディングまたは内挿する前に、FFLベースまたは他の投影FFRベースのパルスMPIシステムを用いて、投影再構成を最初に、FFL/FFRと画像ボリュームとの間の異なる相対角度で収集した未処理の時間ドメインデータセットに適用してもよい。このような方法では、時間ドメインデータを画像ドメイン値に低減する前に、時間ドメインデータを、3つのすべての空間の寸法に分解した場所に対応付けてもよく、その結果、動的なフィッティングの良条件化が改善されて、3D空間内での真の緩和挙動がより良好に局在化され得る。
いくつかの実施形態では、パルスシーケンスによって、所与のFFRの平均場所の付近のトレーサーに対する定常状態を実現するための時間を、定常状態再生パルス励起波形を用いて調べてもよい。図20に、このタイプのエンコーディング及び関連する生データの実施形態を示す。この場合、定常状態再生または緩和時間の測定値を定量化またはフィッティングしてもよい。定常状態再生励起によって、はるかに長いまたは総合時定数を、はるかに速い時定数及び同時に作用する関連する物理的現象の存在下で、測定可能または検出可能とすることができる。ある状態下のいくつかのMPIトレーサーを、特徴時間が一桁以上異なる同時の緩和プロセスによって特徴付けてもよい。さらに、遅いプロセスがトレーサー応答の全エネルギーのかなりの割合に関与する場合があるが、このプロセスのダイナミクスを未処理の時間ドメイン信号に取り込むことは、速いプロセス及びノイズフロアと並べられると難しい場合がある。この場合、定常状態再生エンコーディングによって、より遅いプロセス、及び/またはプロセスに関連する総合または結合緩和時間を定量化するための堅固な方法が得られる。
いくつかの実施形態では、パルス励起をFFLの線に沿った方向に行ってもよい。この方式では、線は直交面内では平行移動しないが、その代わりに、FFLからLFLへ、またはある強度もしくは分極におけるLFLから別の強度もしくは分極におけるLFLへ移行する。線方向に印加された励起パルスとは関係なく、直交面内でのFFLまたはLFLの平均場所は一定である。いくつかの実施形態では、異なる印加磁場強度において線に沿ってトレーサーの緩和を調べることが可能である。AWRにおいて緩和マップをパルスエンコーディングによって構成するときまたは勾配がないときと同様に、励起パルスシーケンスを用いてトレーサー分布の緩和ダイナミクスを全印加磁場強度に対して調べてもよい。しかしこの情報は、ここでは励起線によって空間的に局在化される。この励起をFFR及び/またはFFLの異なる平均場所で適用してもよい。投影モードにおいてまたは完全なトモグラフィ投影再構成に続いて、1D緩和マップまたは他の非スカラ値を個々のピクセル及び/またはボクセルに対応付けてもよい。概略的に、他の緩和画像実施形態の場合と同様に、これらの多値データセットをメリットのスカラ値に低減してスカラ緩和画像を得てもよい。またこれらのデータを、空間及び時間緩和に関連する変数の関数としてトレーサー強度を記述する4Dデータセットに低減または変換してもよい。
トレーサー密度イメージングの場合と同様に、時間ドメイン信号自体に固有の空間時間エンコーディングを緩和画像を構成するときに用いてもよい。緩和加重トレーサー密度イメージングで適用されたものと同じかまたは類似の加重、ウィンドウ処理、またはフィルタリング手順を、緩和画像を形成する前に時間ドメイン信号に適用してもよい。いくつかの実施形態では、メリットの緩和パラメータ(たとえば、指数関数的時定数)を、時間ドメイン信号のウィンドウを開けられた多くの別個の分解またはフィルタリングされた別個のバージョンのそれぞれにフィッティングしてもよい。早い態様の読み出し信号ほど速い緩和プロセス(たとえば、FFRアイソセンタからより遠い粒子に起因するもの)によって多く加重され、遅い態様ほど遅い緩和プロセス(たとえば、FFRアイソセンタにより近い粒子に起因するもの)によるため、異なるウィンドウまたはフィルタに対応付けられる緩和パラメータの測定値によって、特定の現象に対応付けられる緩和時間に対して除去または選択するメカニズムが得られる。いくつかの実施形態では、この手順によって緩和画像の分解能が向上する。いくつかの実施形態では、空間加重とは関係なく、この手順を用いて別個の緩和状態(たとえば、結合対未結合)に対して除去または選択することができる。
別個のトレーサー密度画像及び緩和状態画像を形成するある実施形態について説明してきた。数学的に、これらは、形式ロー(x、y、z)及びr(x、y、z)の別個の画像をもたらすと記述することができる。ここで、ローはスカラトレーサー密度画像に対応し、rはスカラ緩和画像に対応する。概略的に、より明示的に情報を組み合わせて形式I(x、y、z、タウ)の4D画像データセットを形成することができてもよい。この形式では、再構成されたデータセットにはトレーサーについての強度情報が空間及び緩和時間の関数として含まれる。このような再構成によって、点広がり関数(PSF)及び画像(MPIプロセスをより十分に特性評価してトレーサーダイナミクスについてのスペクトル情報を与える)が得られる。いくつかの実施形態では、4D画像がタウ方向に沿って投影された場合、スカラ強度変数Iは崩壊してトレーサー密度ローになる。いくつかの実施形態では、未処理の時間ドメインデータに適用されるウィンドウを開けられた分解技術を用いて、このような4D画像データセットを形成してもよい。他の実施形態では、ラプラス変換などの数学的ツールを利用して、緩和時間または緩和スペクトルデータを再生し、分解し、または再構成してもよい。他の実施形態では、未処理の時間ドメインデータを好適な基底系上に投影してもよいし、または変換した後に好適な基底系上に投影してもよい。基底系を、たとえば、トレーサーを特性評価する既知の物理的緩和プロセスに対応付けてもよい。
パルスMPIの状況で種々の後処理技術を利用してもよい。概略的に、各スキャンに対して信号プロセッサによって1つ以上の画像を再構成してもよく、また単一のサンプルに1つ以上のスキャンを対応付けてもよい。このように得た画像のセットを種々の仕方で処理して、ユーザに対する最終出力を形成してもよい。たとえば、単一画像を処理し、フィルタリングし、または閾値処理してもよい。複数の画像を、何らかの方法で組み合わせる前に、同様に処理し、フィルタリングし、または閾値処理してもよい。
いくつかの実施形態では、複数の画像を、異なるネイティブ分解能及び/またはSNRレベルをもたらすエンコーディングまたは再構成から得てもよい。たとえば、動的に変化する勾配、異なる固定値の勾配大きさで取った複数のスキャンを用いること、動的に変化するパルス励起振幅、異なる固定値の励起振幅で取った複数のスキャン、及び/または異なるパルスMPIエンコーディングを用いて(たとえば、異なる磁化準備デザインを通して)取ったスキャンを用いることをすべて組み合わせて、より小さいセットまたは単一の画像を形成してもよい。最終画像(複数可)は、分解能、SNR、及び画像コントラストのうちの1つ以上を最大にするようにデザインしてもよく、また所望の応用例に依存してもよい。
いくつかの実施形態では、1つ以上の別個の密度及び/または1つ以上の別個の緩和画像を組み合わせて1つ以上の新しい画像にしてもよい。いくつかの実施形態では、緩和情報を用いてトレーサー密度画像をカラー化して、画像コントラストを向上させ、緩和情報をトレーサー密度画像上にマッピングする。緩和情報を、連続的に変化するカラーマップを伴う別の様式から、密度画像及び/または解剖学的参考画像上に直接重ねてもよい。いくつかの実施形態では、他の手段によって、連続的に変化する緩和情報を閾値処理するかまたは別個のカラー上にマッピングしてもよい。たとえば、順モデルベースのアルゴリズムを用いて、対象とする緩和種を分離してもよい(たとえば、結合または未結合状態のMPIトレーサーを、またはサンプル体積に複数が適用される場合の別個のトレーサーの分離において)。これらの結果を用いてスタンドアローンの画像を形成してもよいし、またはカラーマッピングとして密度画像上に適用して重ねてもよい。
いくつかの実施形態では、信号処理及び再構成によって、MPI強度値が3つの空間座標変数及び時間緩和変数に対してマッピングされる4DMPIデータセットが得られる。種々の寸法に沿って投影することによって、情報が豊富なデータセットをシームレスに照合することができてもよい。たとえば、1つ以上の空間寸法に沿って投影すると、残りの空間寸法に対して緩和スペクトル情報を視覚化することができる。緩和寸法に沿って投影することでトレーサー密度画像が再生され得る。4Dデータセットを通るスライスがあると、容易に視覚化された3Dデータセット(たとえば、空間内に局在化されたMPI強度)が、緩和スペクトル変数の列挙を通ってステップするときに与えられ得る。
より大きいコアトレーサーを用いることによるイメージングの向上
パルスエンコーディング及び再構成技術を用いれば、MPIの現在の応用例がすべて、イメージングパラメータ(たとえば、分解能、SNR、及び画像コントラスト)の点に関して改善され得る。たとえば、平均コアサイズが25nmよりも大きいトレーサーをパルスMPI技術とともに用いれば、アクティブなMPI開発の現在領域(たとえば、幹細胞追跡、癌診断イメージング、肺換気及び灌流イメージング、心血管造影応用例、脳卒中診断イメージング、及び灌流イメージング、たとえば脳血流量(CBF)及び脳血液量(CBV))において分解能を向上させることができる。
パルスMPIの固有の態様(たとえば、緩和情報の定量化、4Dイメージングの可能性、及び磁化準備を通して磁化成形ができること)によって、MPIのこれまでに予測されていない新しい応用例が得られ得る。緩和情報、特に環境に敏感な物理過程(たとえば、ブラウン緩和)に対応付けられるものを頑強に定量化できることを用いて、これらのパルスシーケンスを、新しい分子イメージング応用例の広範囲にわたって用いてもよい。パルスMPIパルスシーケンスにより、任意的に平均磁気コアサイズが25nm以上のトレーサーを用いて、以下のすべてを可能にすることができる。インビボ粘度測定イメージング、pHイメージング、炎症イメージング、活性酸素種イメージング、機能化されたトレーサーを用いた結合事象のイメージング、MPI磁気スイッチの設計及び使用、機能化されたまたは専用のトレーサーを介した生化学的動力学の観察、身体内の出血を明らかにすること、感染イメージング、肺塞栓症及び肺の灌流障害の診断、センチネル節生検。
実際にはパルスMPI励起を実現するには複数のハードウェアシステムが必要である。図42に、パルスMPIスキャニングシステムにおける構成サブシステムのコンポーネント及びそれらの間の関係を記載する。図43にはさらに、典型的なパルスMPIスキャニングシステムの態様を例示する。
図2に例示するように、すべてのMPIシステムの最も重要な点はあるタイプのFFR構造を出現させることを通して信号の空間位置を確認することである。FFR構造を形成するために、典型的なMPIスキャニングシステムでは、1つ以上のアクティブまたはパッシブな磁場源を用いて、FFRを構成する特定の空間磁場パターンを発生させる。FFRの特徴は、印加磁場がトレーサー固有の飽和値を下回る少なくとも1つの領域と、印加磁場がトレーサー固有の飽和値を上回る1つの領域とを含むことである。FFRは多くの場合に、必ずしもそうではないが、磁場に線形空間勾配を含むため、これらの領域の間を滑らかで直線的に移行する。FFRは、無磁場点、無磁場線、または図2に示したようなより一般的な形状を有していてもよい。FFRは、図44に例示するように、1つ以上のアクティブまたはパッシブな磁気源(たとえば、電磁石コイル、超伝導磁石、または永久磁石のアレイ)によって生成してもよい。また、図44に示すようなパッシブな軟磁性材料(たとえば鉄)を磁束集中ガイドとしてデザインで用いて、FFRパターンの望ましい態様を改善してもよい(たとえば、電力を低減した状態での直線勾配強度)。
図17に示すように、典型的なMPIスキャンは、対象とするイメージングFOVに対してFFRパターンの場所をシフトさせて、対象とするイメージングFOVをサンプリングすることと考えてもよい。FFLを用いるこの典型的なx空間パルスシーケンス実施形態では、線に直交する平面内でFFL線アイソセンタの平均場所をゆっくりとラスター化して、投影画像を取得する。FFLサンプル相対ジオメトリが回転すると、トモグラフィ投影再構成に対して複数の投影を得ることができる。図のクロスハッチ領域は、線に直交する平面内で励起パルスシーケンスによってサンプリングまたはカバーされたpFOV領域を例示する。代替的に、いくつかの実施形態では、励起は全体的に線の方向であってもよい。より大きいFOVを十分にサンプリングするためにはシフティング波形が必要である。これらの領域は説明用であり、及びまたは必ずしも一定の比率ではない。たとえば、いくつかの実施形態では、pFOV領域幅は相対的なサイズがはるかに小さくてもよい。さらに、図示した特定のラスター軌跡は限定ではない。他のラスターまたはデカルト軌跡、ならびに半径方向、螺旋状、または他のタイプの軌跡を用いて、全体FOVをカバーしてもよい。またpMPI励起を、システム行列再構成法用の較正またはシステム行列を作成するときに用いることができる。
パルスMPIでは、速い励起または送信ハードウェアシステムによって、パルス波形によって記述される時間変化する磁場を実現できなくてはならない。典型的なMPI軌跡には、速くて時間的に正弦波的に変わる励起または駆動磁場が含まれるが、本発明のいくつかの実施形態では、非正弦波スキャニング軌跡をもたらす非正弦波的な時間変化を伴うパルス励起波形を用いることを考える。これらの非正弦波励起波形には、磁場の特定値がある時間の間だけ実質的に一定に保持される周期的または非周期的な成分、ならびに立ち上がり及び/または立ち下がりエッジを通る急速な移行、ならびにより多くの任意の過渡的な波形成分が含まれていてもよい。
いくつかの実施形態では、図45に例示するように、送信、駆動、または励起エレクトロニクスチェーンにはLR回路が含まれ、このLR回路は、このインダクタンスを実質的に変えない非共振フィルターチェーン内の線形増幅器によってパワー供給される。電力増幅器を、制御された電圧または制御された電流モードで実行してもよい。図45に例示するように、制御された電圧モードで動作するとき、線形電力増幅器に供給される入力電圧信号にプリエンファシスを適用してもよい。所望のパルス状電流波形(急速な過渡事象及び実質的に一定の時間が含まれ得る)を実現するために不可欠な電圧波形を、電力増幅器に入力する前に事前に計算する。これは、回路のシミュレーション、数学的な最適化技術(たとえば、凸最適化定式化)、及び/またはスキャンを実行する前の測定による回路特性評価によって行ってもよい。概略的に、線形増幅器を伴うLR回路にパワー供給するとき、最大のスルーレートはシステムのLR時定数によって決定される。したがって、適切なスルーレートを得るために低インダクタンスの送信/励起/駆動コイルをデザインすることが望ましいかまたは必須であり得る。
いくつかの実施形態では、切替器回路を用いてパルス励起波形を実現してもよい。このアプローチでは、送信回路システムが、図46に示すように正確なタイミング図に従って種々の状態の間で切り替わる。送信コイルにキャパシタを共振周期の一部の間接続して、送信コイルの電流を、キャパシタの接続を切る前に、高速な切り替え時間による制御可能な切り替えにより急速に変えて、コイル内の実質的に一定の電流をグライディング時間にわたって実現する。所望のスルーレートを実現するようにキャパシタのキャパシタンスを選択して、可能性としてコイルのインダクタンスにもっとフレキシビリティが与えられるようにしてもよい。実質的に一定の領域間で電流が移行する必要がある度に、キャパシタを短時間で再接続する。切り替える正確なタイミングが、回路動作にとって非常に重要であり、適切にデザインされた切り替えによって可能になる。このような切り替えは、たとえば10マイクロセカンドよりも良好であるか、または1マイクロセカンドよりも良好であるか、または0.1マイクロセカンドよりも良好である時間分解能によって制御することができる。図46に、所望の特徴、結果として生じる励起コイルを時間とともに流れる電流、及び各別個のステージにおける回路トポロジを形成するために用いる典型的な切り替えタイミング図を示す。これらのデザインは教育的で典型的であることが意図されており、決して限定ではない。
パルス波形を送信、駆動、または励起コイルに印加してもよい。これらの励起コイルは全般的に、空間的に均一だが時間変化する磁場をイメージングFOVに与えるようにデザインされる。パルスMPIにおけるある望ましい信号エンコーディングストラテジ及びMPI軌跡を容易にするために、本開示では複数の送信コイル配置を考える。ある回転または投影エンコーディングストラテジをサポートするために、互いに直交する主軸を伴う励起コイルを有することが望ましい場合があり得る。たとえば、FFLMPIシステムにおいて線の方向をy軸として選ぶ場合、図3に示すように、励起コイルをx軸及びz軸と同心で配置してもよい。図14に例示するように、これらのコイルを直交モードで駆動して、励起またはシフトパルスをxz平面の任意の方向に沿って与えてもよい。FFPシステムの場合、少なくとも3つの直交励起コイル(設定された3つの主軸のそれぞれと位置合わせされている)のシステムであれば、励起方向の点に関して類似のフレキシビリティを得ることができる。
図17及び図18に例示するように、一般的なパルスMPI取得にはパルス励起波形を繰り返し印加することが伴い、一方で、シフト波形によって、時間とともに大きなイメージングFOVのサンプリングを可能にする全体x空間軌跡が得られる。これらのパルスシーケンス図に示すシフト波形は、イメージングFOV及びサンプルに対するFFRアイソセンタ及び/または全体FFR構造場所の時間を通じての動きを指す。この相対的な動きは複数の方法で行ってもよい。いくつかの実施形態では、電磁石の中心軸に沿って全体FFR構造をシフトさせるように配設された高出力の均一な電磁石にシフト波形を印加する。いくつかの実施形態では、サンプル自体を機械的手段によって(たとえば、サンプルホルダー、可動ステージ、及びモータの配置を用いて)動かしてもよい。いくつかの実施形態では、サンプルの機械的な動き及びFFRの電磁石シフトティングの両方を全体x空間イメージングスキャンにおいて用いる。図43に、これに関連して、機械装置及びシフト磁石をどのように配設し得るかを例示する。
MPIスキャニングシステムの多くの態様において動作中に冷却する必要があり得る。いくつかの実施形態では、送信または励起コイルならびにシフト磁石及びFFR発生電磁石を積極的に冷却する。いくつかの実施形態では、渦電流シールディングシステムの態様もまた冷却してよい。いくつかの実施形態では、水または代替的な冷却流体(たとえば、フルオロカーボンベースのフロリナート)を用いて、中空の電磁石ワイヤ及びコイルを冷却する。いくつかの実施形態では、要素を熱流体槽内に配置することによって冷却する。
このような大きいスルーレートを用いるときには、対象物に対する比吸収率(SAR)考察の点について注意しなければならない。場合によっては、磁気刺激効果もまた適用され得る。パルスMPI内の高スルー時間は定義上は短いため、単一の励起時間内の小さい負荷サイクルを表す場合がある。パルス励起波形の周期性がより大きい状況では、パルスMPI波形に、たとえば、2乗平均平方根(RMS)という意味で計算されたものと同じかまたは類似のSAR制限測定を施す。
MPIでは、誘導受信部コイルを用いて磁気ナノ粒子トレーサーから磁化信号を検出する。従来、信号を受信する間、連続的な正弦波励起波形が印加される。送受信コイルは強く結合されているため、全般的に、この結果、著しい励起フィードスルーが受信部コイル内に入る。緩和がない場合は、このフィードスルーはトレーサー信号よりも何桁も大きい。幾何学的デカップリングストラテジ(たとえば、使用する受信部コイルが、送信コイルによって規定されるジオメトリに対して正味面積がゼロの磁気傾度測定用で巻かれている)が、典型的に用いられる。
パルスMPIは、励起及び受信チェーンにおいて従来のMPI方法から大きくずれることが要求される。ハードウェアの他の態様は維持してもよい。たとえば、磁場源コンポーネントとしては、FFRを発生させ、FFRをゆっくりと(励起時間と比べて)シフトさせて大きいイメージングFOVをカバーするものを、正弦波MPIシステムで用いたものと同じかまたは類似の方法で構成して使用してもよい。図51に例示するように、モジュラパルスMPI励起及び受信システムがモジュラ正弦波励起及び受信システムとスワップ可能または切り替え可能であるMPIスキャニングシステムをデザインすることが望ましい場合がある。このようなデザインであれば、特性の異なるパルスMPI励起/受信システムのモジュラ交換または切り替えも容易になる。これは、このようなシステムによって基本パラメータ(たとえば、システム応答時間及びイメージングSNR)がトレードオフされた場合、望ましい可能性がある。主要なデザイン関係はシステム応答時間とトレーサー磁気緩和時間との間にあるため、あるパルスMPI励起/受信システムを、特定のトレーサーまたはトレーサーの特定の応用例に最適にマッチングさせてもよい。
ここで、理想的な方形波励起を用いた1次元のパルスMPI法の基本的な導出を示す。この分析は少数の理想的なパルスMPI法に対するものであり、基本的考え方のいくつかをさらに例示することが意図されており、決して限定ではない。
Claims (56)
- パルス磁気粒子イメージングシステムであって、
少なくとも1つの磁石を含む磁場発生システムであって、前記磁気粒子イメージングシステムの観察領域内に空間的に構造化された磁場を設けて、前記空間的に構造化された磁場が、内部に磁気ナノ粒子トレーサー分布を有する観察物体に対する無磁場領域(FFR)を有するようにする前記磁場発生システムと、
前記観察領域に隣接して配設されたパルス励起システムであって、電磁石とパルスシーケンス発生器とを含み、前記パルスシーケンス発生器は前記電磁石に電気的に接続されて前記電磁石に励起波形を与え、前記電磁石は、前記励起波形が与えられると、前記FFRの場所または状態をシフトさせることの少なくとも一方によって、前記観察領域内に励起磁場を発生させてそこから励起信号を誘起する、前記パルス励起システムと、
前記観察領域に隣接して配設された検出システムであって、前記励起信号を検出して検出信号を得るように構成された前記検出システムと、を含み、
前記励起波形は過渡的部分と実質的に一定の部分とを含む、前記パルス磁気粒子イメージングシステム。 - 少なくとも前記磁場発生システムと前記パルス磁気粒子イメージングシステムの周囲の環境とから前記観察領域を電磁的に分離するように前記観察領域を内部に囲むように配設された電磁シールドをさらに含む請求項1に記載のパルス磁気粒子イメージングシステム。
- 少なくとも前記パルス励起システムの送信部分と前記検出システムの受信部分とは前記電磁シールド内に囲まれる請求項2に記載のパルス磁気粒子イメージングシステム。
- 前記磁場発生システムはさらに、軟磁性材料を含むパッシブな磁場収束素子であって、内部で磁力線を所望の形状に収束または形成するような形状に構成され前記少なくとも1つの磁石に対して配設された前記磁場収束素子含む、請求項1に記載のパルス磁気粒子イメージングシステム。
- 前記検出システムは、実質的に前記励起波形の前記実質的に一定の部分の間にのみ前記励起信号を検出して前記検出信号を得るように構成されている請求項1に記載のパルス磁気粒子イメージングシステム。
- 前記励起波形の前記実質的に一定の部分は少なくとも500ナノ秒で500ミリ秒未満である請求項1に記載のパルス磁気粒子イメージングシステム。
- 前記励起波形の前記実質的に一定の部分は、前記励起波形の目標振幅の約10%以内で一定である請求項1〜6のいずれか1項に記載のパルス磁気粒子イメージングシステム。
- 前記励起波形の前記過渡的部分は継続時間が少なくとも100ナノ秒で100マイクロセカンド未満である請求項1に記載のパルス磁気粒子イメージングシステム。
- 前記励起波形は磁化準備部分と読み出し部分とを含み、前記磁化準備部分は少なくとも前記過渡的部分の一部を含み、前記読み出し部分は少なくとも前記一定部分の一部を含む請求項1〜8のいずれか1項に記載のパルス磁気粒子イメージングシステム。
- 前記磁化準備部分は、前記FFRの付近にあるトレーサー磁化の状態を、前記読み出し部分の前の前記トレーサーの磁気緩和特性に基づいて動的に構成する請求項9に記載のパルス磁気粒子イメージングシステム。
- 前記磁化準備は、前記読み出し部分の間の特定の緩和状態に対応付けられるトレーサーからの信号を選択的にゼロにする請求項10に記載のパルス磁気粒子イメージングシステム。
- 前記磁化準備は前記トレーサー磁化をエンコードして、磁気緩和特性に基づいて信号を選択的に減衰させるかまたはゼロにする請求項10に記載のパルス磁気粒子イメージングシステム。
- 前記読み出し部分は、前記FFRにおける前記磁気ナノ粒子トレーサーに対する緩和時間よりも長い請求項9に記載のパルス磁気粒子イメージングシステム。
- 前記読み出し部分は、前記FFRにおける前記磁気ナノ粒子トレーサー内に定常状態磁化を設定するのに十分に長い請求項9に記載のパルス磁気粒子イメージングシステム。
- 前記読み出し部分は、前記FFRにおける前記磁気ナノ粒子トレーサーに対する緩和時間よりも短い請求項9に記載のパルス磁気粒子イメージングシステム。
- 前記励起波形は複数のパルスを含み、前記複数のパルスの各パルスは過渡的部分を含む、請求項1に記載のパルス磁気粒子イメージングシステム。
- 前記励起波形は複数の一定部分を含み、
前記励起波形は複数の磁化準備部分と複数の読み出し部分とを含み、前記複数の磁化準備部分の各磁化準備部分が、前記複数のパルスの少なくとも1つにある少なくとも1つの過渡的部分を含み、前記複数の読み出し部分の各前記読み出し部分が、前記複数の一定部分の少なくとも1つの一定部分の少なくとも一部を含む、請求項16に記載のパルス磁気粒子イメージングシステム。 - 前記FFRは長手方向を規定する無磁場線であり、前記励起波形は前記無磁場線の前記長手方向に少なくとも部分的に印加されて、磁場空間内の対応する無磁場線構造を移動させるが前記無磁場線構造の形状及び空間的位置を維持することによって、前記無磁場線の状態を変える請求項1〜16のいずれか1項に記載のパルス磁気粒子イメージングシステム。
- 前記FFRは長手方向を規定する無磁場線であり、前記励起波形は前記無磁場線の前記長手方向に直交する平面において少なくとも部分的に印加されて、前記無磁場線の前記位置を変える請求項1〜16のいずれか1項に記載のパルス磁気粒子イメージングシステム。
- 前記パルスシーケンス発生器は、前記複数のパルスを、それぞれが、予め選択された形状、大きさ、幅、またはパルス間時間のうちの少なくとも1つを有して特定のパルスシーケンスエンコーディングをもたらすように、与えるように構成されている請求項16に記載のパルス磁気粒子イメージングシステム。
- 前記パルスシーケンス発生器は、パルス間部分によって分離された複数のパルスを含む励起波形を与えるように構成され、
前記複数のパルスの各パルスは過渡的部分を含む請求項1〜20のいずれか1項に記載のパルス磁気粒子イメージングシステム。 - 前記パルスシーケンス発生器は、一定のパルス間部分によって分離された複数のパルスを含む励起波形を与えるように構成され、
前記複数のパルスの各パルスは過渡的部分の間に実質的に一定の部分を有する請求項1〜20のいずれか1項に記載のパルス磁気粒子イメージングシステム。 - 前記励起波形は有限時間方形波に近似する請求項20に記載のパルス磁気粒子イメージングシステム。
- 前記パルス波形の少なくとも一部に対して、前記複数のパルスの連続パルス間の各一定部分は先行する一定部分よりも長い請求項20に記載のパルス磁気粒子イメージングシステム。
- 前記パルス励起システムは、非共振フィルターチェーン内の線形増幅器によってパワー供給されるLR回路を含む請求項1〜24のいずれか1項に記載のパルス磁気粒子イメージングシステム。
- 前記LR回路はインダクタンスが1マイクロヘンリー〜50マイクロヘンリーである請求項25に記載のパルス磁気粒子イメージングシステム。
- 前記LR回路はインダクタンスが1マイクロヘンリー〜30マイクロヘンリーである請求項25に記載のパルス磁気粒子イメージングシステム。
- 前記パルス励起システムは共振切替器回路を含む請求項1〜24のいずれか1項に記載のパルス磁気粒子イメージングシステム。
- 前記検出システムは、前記検出信号のフィードスルー汚染が多い部分と比べて前記検出信号のフィードスルー汚染が少ない部分を増幅する利得制御回路を含む請求項5に記載のパルス磁気粒子イメージングシステム。
- 前記検出システムと通信してそこから検出信号を受け取るように構成された信号プロセッサをさらに含み、
前記信号プロセッサはさらに、前記FFRが横断した前記観察物体の領域に対応する画像をレンダリングするためのイメージング信号を発生させるように構成されている請求項1〜24のいずれか1項に記載のパルス磁気粒子イメージングシステム。 - 前記磁気粒子イメージング装置の前記観察領域内に配設されるように構成された物体ホルダーをさらに含む請求項1〜30のいずれか1項に記載のパルス磁気粒子イメージングシステム。
- 前記物体ホルダーまたは前記磁場発生システムの少なくとも一方に動作可能に接続されて、前記FFRの相対位置を平行移動させるかまたは回転させることの少なくとも一方を行う機械アセンブリをさらに含む請求項31に記載のパルス磁気粒子イメージングシステム。
- スローシフト電磁石を含むスローシフト電磁システムと、前記磁気粒子イメージング装置の前記観察領域に隣接して配置されたスローシフト波形発生器とをさらに含み、
前記スローシフト波形発生器は、前記励起パルスの時間スケールと比べて遅い時間スケールで前記FFRの位置をシフトさせるために前記スローシフト電磁石に波形を与える請求項1〜32のいずれか1項に記載のパルス磁気粒子イメージングシステム。 - 前記磁場発生システムと通信するように構成されたスローシフト波形発生器を含むスローシフト電磁システムをさらに含み、
前記スローシフト波形発生器は、前記励起パルスの時間スケールと比べて遅い時間スケールで前記FFRの位置をシフトさせるために前記磁場発生システムに波形を与える請求項1〜32のいずれか1項に記載のパルス磁気粒子イメージングシステム。 - 前記磁場発生システムは動的に構成可能である請求項1〜32のいずれか1項に記載のパルス磁気粒子イメージングシステム
- 前記FFRを動的に変えて、前記検出信号においてエンコードされた信号対雑音比及び分解能を動的に変えるように、前記磁場発生システムは動的に構成可能である請求項35に記載のパルス磁気粒子イメージングシステム。
- 前記検出システムと通信して前記検出信号を受信するように構成された信号処理及び画像レンダリングシステムをさらに含み、
前記信号処理及び画像レンダリングシステムは、前記検出信号を処理するように、また前記磁気ナノ粒子トレーサーを含み前記FFRによってアドレスされた前記観察物体の部分に対応する画像をレンダリングするように構成されている請求項1〜34のいずれか1項に記載のパルス磁気粒子イメージングシステム。 - 前記信号処理及び画像レンダリングシステムは、前記検出信号を処理するように、また少なくとも1.5mmの空間分解能で前記画像をレンダリングするように構成されている請求項37に記載のパルス磁気粒子イメージングシステム。
- 前記信号処理及び画像レンダリングシステムは、前記検出信号を処理するように、また少なくとも1000μm〜100μmの空間分解能で前記画像をレンダリングするように構成されている請求項37に記載のパルス磁気粒子イメージングシステム。
- 前記信号処理及び画像レンダリングシステムは、前記検出信号を処理するように、また前記画像をレンダリングして、対応する画像場所における前記トレーサーの密度、質量、濃度、またはその微分のうちの少なくとも1つを表すように構成されている請求項37に記載のパルス磁気粒子イメージングシステム。
- 前記信号処理及び画像レンダリングシステムは、前記検出信号を処理するように、また前記画像をレンダリングして、対応する画像場所における前記トレーサーの磁気緩和ダイナミクスの局所測定を表すように構成されている請求項37または40に記載のパルス磁気粒子イメージングシステム。
- 前記信号処理及び画像レンダリングシステムは、前記検出信号を処理するように、また前記画像をレンダリングして、対応する画像場所における局所粘度を表すように構成されている請求項37または41に記載のパルス磁気粒子イメージングシステム。
- 磁気ナノ粒子トレーサーを用いて物体をイメージングする方法であって、
前記物体に前記磁気ナノ粒子トレーサーを与えることと、
FFRを有する空間的に構造化された磁場を印加して、前記FFRと前記空間的に構造化された磁場の周囲領域とが、前記磁気ナノ粒子トレーサーの少なくとも一部を含む領域において前記観察物体を遮断するようにすることと、
前記FFRの特性または前記FFRの位置を変えることの少なくとも一方によって前記磁気ナノ粒子トレーサーの一部を励起することと、
前記FFRの前記特性及び前記FFRの前記位置が実質的に一定である間に、前記励起に起因する前記磁気ナノ粒子トレーサーの磁化の変化を検出して、検出信号を取得することと、
前記物体内の前記FFRの複数の異なる場所に対して前記励起及び検出を繰り返して、複数の検出信号を取得することと、
前記複数の検出信号を処理して前記物体の領域の画像をレンダリングすることと、を含む前記物体をイメージングする方法。 - 前記観察物体に前記磁気ナノ粒子トレーサーを投与することをさらに含み、
投与された前記磁気ナノ粒子トレーサーに含まれる磁気ナノ粒子は、アンサンブル平均直径が少なくとも10nmで100nm未満である請求項43に記載の物体をイメージングする方法。 - 前記磁気ナノ粒子トレーサーの前記磁気ナノ粒子は、少なくとも25nmで50nm未満である請求項43または44に記載の物体をイメージングする方法。
- 前記磁気ナノ粒子トレーサーの前記磁気ナノ粒子は、5nmの変動内で一様である請求項44または45に記載の物体をイメージングする方法。
- 前記検出中の少なくとも500ナノ秒で500ミリ秒未満の間、前記FFRの前記特性及び前記FFRの前記位置が実質的に一定である請求項43〜46のいずれか1項に記載の物体をイメージングする方法。
- 前記退出中の前記FFRの変化量及び前記FFRの前記位置の約10%以内で、前記FFRの前記特性及び前記FFRの前記位置が実質的に一定である請求項43〜47のいずれか1項に記載の物体をイメージングする方法。
- 前記FFRの前記特性または前記FFRの前記位置を変えることの継続時間は、少なくとも100ナノ秒で10マイクロセカンド未満である請求項43〜48のいずれか1項に記載の物体をイメージングする方法。
- 前記複数の検出信号を処理して前記物体の前記領域の前記画像をレンダリングすることによって、磁気粒子密度画像または磁気緩和動的パラメータ画像の少なくとも一方をレンダリングする請求項43〜49のいずれか1項に記載の物体をイメージングする方法。
- 前記複数の検出信号を処理して前記物体の前記領域の前記画像をレンダリングすることによって、局所粘度をレンダリングする請求項43〜50のいずれか1項に記載の物体をイメージングする方法。
- パルス磁気粒子イメージングシステムとともにまたはその一部として用いるデバイスであって、
サンプル観察領域に隣接して配設されたパルス励起システムであって、電磁石とパルスシーケンス発生器とを含み、前記パルスシーケンス発生器は前記電磁石に電気的に接続されて前記電磁石に励起波形を与え、前記電磁石は、前記サンプル観察領域内に磁場を与えて、サンプルが前記サンプル観察領域内にサンプルホルダーによって保持されたときに前記サンプルから励起信号を発生させる、前記パルス励起システムと、
前記サンプル観察領域に隣接して配設された検出システムであって、前記サンプルから前記励起信号を検出して検出信号を得るように構成された前記検出システムと、を含み、
前記励起波形は過渡的部分と実質的に一定の部分とを含む、前記デバイス。 - 前記パルス励起システムと前記検出システムとは、非パルスMPIシステムをパルスMPIシステムに変換するように適合されたモジュラ構造の一部である請求項52に記載のデバイス。
- 前記検出システムと通信して前記検出信号を受信して処理するように構成された信号プロセッサをさらに含む請求項52に記載のデバイス。
- 前記信号プロセッサはさらに、前記検出信号を処理して前記サンプル内の磁気粒子に対する磁化緩和時間を決定するように構成されている請求項52に記載のデバイス。
- 前記サンプル観察領域を画定するサンプルホルダーをさらに含む請求項52〜55のいずれか1項に記載のデバイス。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201762546395P | 2017-08-16 | 2017-08-16 | |
US62/546,395 | 2017-08-16 | ||
PCT/US2018/000231 WO2019035956A1 (en) | 2017-08-16 | 2018-08-16 | SYSTEMS AND METHODS FOR PULSED MAGNETIC PARTICLE IMAGING |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2020531076A true JP2020531076A (ja) | 2020-11-05 |
Family
ID=65362549
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020508019A Pending JP2020531076A (ja) | 2017-08-16 | 2018-08-16 | パルス磁気粒子イメージングシステム及び方法 |
Country Status (5)
Country | Link |
---|---|
US (3) | US11231469B2 (ja) |
EP (2) | EP3669201B1 (ja) |
JP (1) | JP2020531076A (ja) |
CN (2) | CN111183364B (ja) |
WO (1) | WO2019035956A1 (ja) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11890488B2 (en) * | 2019-03-13 | 2024-02-06 | Magnetic Insight, Inc. | Magnetic particle actuation |
MX2021011111A (es) * | 2019-03-25 | 2021-10-22 | Promaxo Inc | Bobinas de campo de gradiente de formacion de imagenes por resonancia magnetica (mri) de un solo lado y sus aplicaciones. |
CN110975152A (zh) * | 2019-12-17 | 2020-04-10 | 华中科技大学 | 一种可连续工作的磁刺激装置及方法 |
CN111726925B (zh) * | 2020-06-12 | 2022-05-31 | 北方工业大学 | 一种脉冲模式下x射线脉冲剂量波形的优化方法 |
KR102545062B1 (ko) | 2020-06-25 | 2023-06-20 | 한국전자통신연구원 | 나노 자성 입자 영상화 장치 및 방법 |
DE102020118408A1 (de) | 2020-07-13 | 2022-01-13 | RUHR-UNIVERSITäT BOCHUM | Verfahren zum Bestimmen von Gewebeeigenschaften eines biologischen Gewebes |
CN112014639B (zh) * | 2020-09-02 | 2022-07-05 | 安徽一天电能质量技术有限公司 | 一种交流电力谐波向量测量方法 |
CN113129403B (zh) * | 2021-04-19 | 2022-06-10 | 中国科学院自动化研究所 | 基于前向模型的磁粒子成像系统矩阵图像重建方法及系统 |
CN114089426B (zh) * | 2021-05-26 | 2023-11-10 | 华北科技学院(中国煤矿安全技术培训中心) | 一种改进的u形螺线源瞬变电磁全空间定向探测方法 |
CN113331812A (zh) * | 2021-06-03 | 2021-09-03 | 中国科学院自动化研究所 | 面向全身的三维磁粒子成像方法、系统及设备 |
CN113499052A (zh) * | 2021-07-08 | 2021-10-15 | 中国科学院自动化研究所 | 磁纳米粒子成像系统矩阵测量的栅格状探测板及测量方法 |
CN113558597B (zh) * | 2021-07-21 | 2022-05-03 | 北京航空航天大学 | 基于磁滞效应的无磁场点磁纳米粒子成像方法 |
US11733324B2 (en) | 2021-09-24 | 2023-08-22 | Mitsubishi Electric Corporation | Magnetic particle imaging system and magnetic particle imaging method |
CN114533016B (zh) * | 2022-04-22 | 2022-07-19 | 北京航空航天大学 | 基于梯形波激励的多色定量磁纳米粒子成像方法及系统 |
CN115120221B (zh) * | 2022-06-27 | 2024-06-18 | 北京航空航天大学 | 基于布朗弛豫编码的磁纳米粒子成像系统及方法 |
CN115568841B (zh) * | 2022-09-27 | 2023-04-14 | 北京航空航天大学 | 一种基于尼尔弛豫的磁纳米粒子检测与成像方法 |
CN115797493B (zh) * | 2023-02-03 | 2023-04-25 | 北京航空航天大学 | 基于一维系统矩阵稀疏采样的磁场自由线磁粒子成像方法 |
CN116965795B (zh) * | 2023-09-20 | 2023-12-12 | 暨南大学附属第一医院(广州华侨医院) | 一种高灵敏度的磁粒子成像系统及方法 |
CN117572309B (zh) * | 2023-10-20 | 2024-04-09 | 暨南大学附属第一医院(广州华侨医院) | 基于谐波相位偏转的磁粒子谱仪、定量方法和存储介质 |
CN117100243B (zh) * | 2023-10-23 | 2024-02-20 | 中国科学院自动化研究所 | 基于系统矩阵像素压缩的磁粒子成像系统、方法及设备 |
CN117530675B (zh) * | 2024-01-03 | 2024-04-02 | 中国科学院自动化研究所 | 基于弛豫矫正的磁粒子成像分辨率提升系统、方法及设备 |
CN117598680B (zh) * | 2024-01-23 | 2024-05-07 | 辽宁嘉玉科技有限公司 | 磁粒子磁化感知距离测量装置与方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006525506A (ja) * | 2003-04-15 | 2006-11-09 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 状態変数及び状態変数の変化を決定する方法 |
US20100066363A1 (en) * | 2008-08-15 | 2010-03-18 | The University Of Houston System | Imaging method for obtaining spatial distribution of nanoparticles in the body |
JP2012525193A (ja) * | 2009-04-30 | 2012-10-22 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 磁性粒子に影響し及び/又は磁性粒子を検出し、磁気共鳴撮像のための装置及び方法 |
US20150015247A1 (en) * | 2008-06-23 | 2015-01-15 | The Regents Of The University Of California | Magnetic Particle Imaging Devices and Methods |
Family Cites Families (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2613076B1 (fr) * | 1987-03-25 | 1990-05-18 | Thomson Cgr | Procede d'imagerie rapide par resonance magnetique nucleaire |
GB9506909D0 (en) | 1995-04-04 | 1995-05-24 | Scient Generics Ltd | Spatial magnetic interrogation system |
GB9608329D0 (en) | 1996-04-23 | 1996-06-26 | Scient Genarics Ltd | Improved methods for coding magnetic tags |
BR9709864A (pt) | 1996-06-19 | 2000-01-11 | Flying Null Ltd | Detetor para sentir a presença de uma etiqueta magnética e leitora de nota bancária |
GB9619896D0 (en) | 1996-09-24 | 1996-11-06 | Flying Null Ltd | Improvements in or relating to magnetic sensors |
GB9620190D0 (en) | 1996-09-27 | 1996-11-13 | Flying Null Ltd | Improved methods for coding magnetic tags |
ES2196372T3 (es) | 1996-10-09 | 2003-12-16 | Flying Null Ltd | Tecnicas de interrogacion magnetica. |
GB9717574D0 (en) | 1997-08-19 | 1997-10-22 | Flying Null Ltd | Catheter location |
AU8815298A (en) | 1997-08-19 | 1999-03-08 | Flying Null Limited | Improvements relating to surgical devices and their location |
GB9800064D0 (en) | 1998-01-05 | 1998-03-04 | Sentec Ltd | Uni-directional magnetic tag |
GB9805824D0 (en) | 1998-03-18 | 1998-05-13 | Flying Null Ltd | Magnetic patterns |
GB9806923D0 (en) | 1998-03-31 | 1998-05-27 | Flying Null Ltd | Position sensing |
GB9817803D0 (en) | 1998-08-14 | 1998-10-14 | Flying Null Ltd | Magnetic information carriers |
DE69909842T2 (de) | 1998-12-23 | 2004-01-29 | Flying Null Ltd | Lesevorrichtungen für magnetische etiketten |
US6812842B2 (en) * | 2001-12-20 | 2004-11-02 | Calypso Medical Technologies, Inc. | System for excitation of a leadless miniature marker |
WO2004091390A2 (en) * | 2003-04-15 | 2004-10-28 | Philips Intellectual Property & Standards Gmbh | Device and method for examination and use of an electrical field in an object under examination containing magnetic particles |
WO2008099331A1 (en) * | 2007-02-15 | 2008-08-21 | Philips Intellectual Property & Standards Gmbh | Arrangement for magnetic particle imaging, method for influencing and/or detecting magnetic particles and magnetic particle |
EP2375976B1 (en) * | 2008-12-10 | 2016-05-11 | Philips Intellectual Property & Standards GmbH | Arrangement with variable selection field orientation for magnetic particle imaging |
WO2010070558A1 (en) * | 2008-12-19 | 2010-06-24 | Philips Intellectual Property & Standards Gmbh | Arrangement for imaging an object including a vessel using magnetic particle imaging |
CN102469952B (zh) | 2009-08-07 | 2014-11-12 | 皇家飞利浦电子股份有限公司 | 用于确定至少一个电磁量的设备和方法 |
CN102481111B (zh) * | 2009-08-21 | 2014-12-17 | 皇家飞利浦电子股份有限公司 | 用于生成和移动具有无场线的磁场的设备和方法 |
WO2011030247A1 (en) * | 2009-09-11 | 2011-03-17 | Koninklijke Philips Electronics N.V. | Apparatus and method for influencing and/or detecting magnetic particles in a field of view |
WO2011116229A2 (en) * | 2010-03-17 | 2011-09-22 | The Regents Of The University Of California | Magnetic particle imaging devices and methods |
RU2594815C2 (ru) * | 2010-12-10 | 2016-08-20 | Конинклейке Филипс Электроникс Н.В. | Устройство и способ для воздействия на и/или детектирования магнитных частиц |
CN103892809B (zh) * | 2014-04-01 | 2016-02-24 | 华中科技大学 | 一种磁纳米温度成像方法及系统 |
KR101623116B1 (ko) | 2015-06-25 | 2016-05-23 | 서울대학교 산학협력단 | 영상 획득 방법 |
CN105054932B (zh) * | 2015-07-03 | 2017-07-18 | 华中科技大学 | 一种基于磁纳米粒子一次谐波幅值的成像方法 |
DE102015214071B3 (de) | 2015-07-24 | 2016-09-22 | Bruker Biospin Mri Gmbh | MPI-Verfahren |
-
2018
- 2018-08-16 JP JP2020508019A patent/JP2020531076A/ja active Pending
- 2018-08-16 EP EP18846668.4A patent/EP3669201B1/en active Active
- 2018-08-16 CN CN201880064884.3A patent/CN111183364B/zh active Active
- 2018-08-16 CN CN202310859434.6A patent/CN116879818A/zh active Pending
- 2018-08-16 EP EP23199044.1A patent/EP4309575A3/en active Pending
- 2018-08-16 US US15/998,525 patent/US11231469B2/en active Active
- 2018-08-16 WO PCT/US2018/000231 patent/WO2019035956A1/en unknown
-
2022
- 2022-01-21 US US17/581,684 patent/US11709212B2/en active Active
-
2023
- 2023-06-05 US US18/205,937 patent/US20240103103A1/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006525506A (ja) * | 2003-04-15 | 2006-11-09 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 状態変数及び状態変数の変化を決定する方法 |
US20150015247A1 (en) * | 2008-06-23 | 2015-01-15 | The Regents Of The University Of California | Magnetic Particle Imaging Devices and Methods |
US20100066363A1 (en) * | 2008-08-15 | 2010-03-18 | The University Of Houston System | Imaging method for obtaining spatial distribution of nanoparticles in the body |
JP2012525193A (ja) * | 2009-04-30 | 2012-10-22 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 磁性粒子に影響し及び/又は磁性粒子を検出し、磁気共鳴撮像のための装置及び方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2019035956A1 (en) | 2019-02-21 |
EP3669201A1 (en) | 2020-06-24 |
US11709212B2 (en) | 2023-07-25 |
US20190079149A1 (en) | 2019-03-14 |
CN111183364B (zh) | 2023-08-04 |
WO2019035956A9 (en) | 2019-04-18 |
EP4309575A2 (en) | 2024-01-24 |
EP4309575A3 (en) | 2024-04-10 |
US20220221537A1 (en) | 2022-07-14 |
CN111183364A (zh) | 2020-05-19 |
CN116879818A (zh) | 2023-10-13 |
EP3669201B1 (en) | 2023-10-25 |
US11231469B2 (en) | 2022-01-25 |
US20240103103A1 (en) | 2024-03-28 |
EP3669201A4 (en) | 2021-05-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2020531076A (ja) | パルス磁気粒子イメージングシステム及び方法 | |
Tay et al. | Pulsed excitation in magnetic particle imaging | |
KR100329467B1 (ko) | 자기공명단층촬상장치 | |
KR101301490B1 (ko) | 자기공명영상장치 및 확산강조영상획득방법 | |
AU2015229128A1 (en) | System and method for free radical imaging | |
US20150115956A1 (en) | System and method for quiet magnetic resonance imaging | |
JP6162142B2 (ja) | 磁気共鳴イメージング装置、及び、sarの予測方法 | |
CN103997958A (zh) | Mpi中的背景去除 | |
JP2017513627A (ja) | 縮小視野磁気共鳴イメージングのシステムおよび方法 | |
US10698053B2 (en) | System and method for gradient-modulated sweep imaging with fourier transformation magnetic resonance imaging | |
Galante et al. | Fast room temperature very low field-magnetic resonance imaging system compatible with magnetoencephalography environment | |
Jia et al. | Gradient-based pulsed excitation and relaxation encoding in magnetic particle imaging | |
CN104797189A (zh) | Mpi中的动态背景校正 | |
KR900007543B1 (ko) | 기본 자장을 펄스식으로한 핵자기공명 영상형성 장치 | |
RU2063702C1 (ru) | Способ магнитно-резонансной томографии и устройство для его осуществления | |
Bazzi | Quantitative magnetic resonance imaging methodology development | |
US10436864B1 (en) | Magnetic resonance imaging | |
JP5562010B2 (ja) | 磁気共鳴イメージング装置 | |
McDonough | Design and Characterization of a Single Sided Field Free Linescanner for Magnetic Particle Imaging | |
JP5037956B2 (ja) | 磁気共鳴イメージング装置 | |
CN118021275A (zh) | 基于焦斑调制的磁粒子成像方法及系统、图像处理装置 | |
JP2023046466A (ja) | 磁気共鳴イメージング装置、および、画像処理方法 | |
Goodwill | Narrowband and x-space magnetic particle imaging | |
JPH06114030A (ja) | Mrイメージング装置 | |
JP2002000580A (ja) | Mr画像撮影装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20210810 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20220630 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20220705 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20220930 |
|
RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20221104 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20221116 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20221205 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20230328 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20230628 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20230926 |