JP2010538718A - 超複素数を使用したmri空間符号化 - Google Patents

超複素数を使用したmri空間符号化 Download PDF

Info

Publication number
JP2010538718A
JP2010538718A JP2010524529A JP2010524529A JP2010538718A JP 2010538718 A JP2010538718 A JP 2010538718A JP 2010524529 A JP2010524529 A JP 2010524529A JP 2010524529 A JP2010524529 A JP 2010524529A JP 2010538718 A JP2010538718 A JP 2010538718A
Authority
JP
Japan
Prior art keywords
signal
space
complex
acquired
point
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
Application number
JP2010524529A
Other languages
English (en)
Inventor
グレニール,デニス
Original Assignee
セントレ ナショナル デ ラ リシェルシェ サイエンティフィック(セ・エン・エル・エス)
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by セントレ ナショナル デ ラ リシェルシェ サイエンティフィック(セ・エン・エル・エス) filed Critical セントレ ナショナル デ ラ リシェルシェ サイエンティフィック(セ・エン・エル・エス)
Publication of JP2010538718A publication Critical patent/JP2010538718A/ja
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/5608Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本発明は、複素信号を処理する方法であって、信号を複素数の形式で取得する工程と、関連付けられた超複素成分を、取得された複素信号に基づいて求める工程であって、上記成分が、少なくとも、取得された複素信号の位相の時間に対する導関数に対応する工程と、それによって求められた超複素信号を処理する工程であって、上記処理から結果として生じる信号が、取得信号の成分の数よりも大きな成分数を含む工程とを含む方法に関する。

Description

本発明は、超複素数による信号の取得及び処理に関し、効果的には、核磁気共鳴の分野に適用される。
核磁気共鳴における信号処理の一般的な理論的事項
核磁気共鳴(NMR)は、複素系の縮退を削減する(すなわち、複素系について、存在する分子、その量、及びその位置を空間中において求める)ことを可能にする手法である。
NMR実験では、それがどのようなものであっても、信号が、陰時間関数又は陽時間関数として取得される。
測定される量は、磁化(回転磁気モーメントの和)であり、小さな差分要素によってもたらされるdSで表すNMR信号は、以下の関係によって表される。
Figure 2010538718
ここで、ρ(x,y,z,v,t)は、求められる量(位置(x,y,z)において観測される核のスピン密度)を表し、vは、磁気モーメントの回転の周波数が、化学的環境に応じた遮蔽係数で修正され得るということ(NMR分光法)を考慮に入れるために導入され、f(x,y,z,v,t)は、既知であり、陰時間関数又は陽時間関数として定義される。
信号dS(x,y,z,v,t)は、回転磁気モーメントの位相及び強度を表すので、複素(複素数フィールドの要素)である。
サンプルを全体とみなした場合、合計取得信号は、
Figure 2010538718
に比例し、ここで、
Δx、Δy、Δz、Δvは、x方向、y方向、z方向及び周波数方向それぞれに沿った信号の積分範囲であり、
信号S(t)は複素信号であり、複素信号の和とみなされる。
取得された信号の各点は、2つの成分(すなわち、強度成分及び位相成分)を有する。
前述の成分は何れも、取得にわたって測定され、測定された位相は、ラジアンで表され、よって、無次元である。
NMRでは、時間の関数としての取得された信号から、生成した関数ρ(x,y,z,v,t)を求めることが可能である。
特に、NMRでは、取得された信号から画像を再構成することが可能であり、取得された信号の各点は、(患者の体などの)画像オブジェクトの点に対応する。
信号から画像を再構成するということは、
S(t)=Mρ(x,y,z,v,t)
というタイプの等式系を解くということになる。ここで、Sは既知のベクトルであり、その要素は複素数であり、その強度及び位相は、既知であり、信号の経時的なサンプリングによって表され、ρは、求める対象のベクトルであり、Mは行列であり、その要素はそれ自体がやはり既知である。実験条件によって規定されるからである。
行列Mの各行の各要素は、特定の取得時間において再構成される対象のボリュームの各要素の位相によって表される。
この目的で、NMR画像の従来の再構成は、関数f(x,y,z,v,t)を、
Figure 2010538718
として書き換えることができるということを用いることによって実現される。
実際に、前述の条件下では、指数が無次元であるため、値k(t)、k(t)、k(t)、k(t)は、x、y、z、vの共役変数であり、S(t)は、k(t)、k(t)、k(t)、k(t)として書き換えることができる。S及びρはフーリエ変換によって関係付けられる。
よって、NMR取得の目標は、取得された信号S(t)を、S(k,k,k,k)の形式にすることができるということを確実にすることである。
この条件は、磁場グラディエントにより、k個の空間(又は、更にフーリエ空間)を「走査する」ことにより、NMR撮像について満たされる。それにより、十分な数の座標点(k,k,k,k)の取得が可能になり、ボリュームをx,y,z,v空間において再構成することができ、それにより、取得された信号の縮退を削減することができる。
よって、高空間分解能で画像をすばやく取得することができることは、k空間を、広い範囲の周波数にわたってすばやく包含する(高強度の場のグラディエントで、すばやく切り換える)ことによって調節される。
図1a及び図1bは、2D画像の場合における「エコー・プラナー・イメージング」シーケンスの範囲内のk空間の走査を示す。
前述のシーケンス中、横磁化が無線周波数RFによってもたらされる(図1aを参照)。
マーク1とマークNとの間の、経時的な、取得信号のK個の点が信号S(t)を表す。
グラディエントG(t)及びG(t)の印加により、信号S(t)が、k空間内のデカルト座標をたどることになる(図1a及び図1bを参照されたい。)。
前述のシーケンスでは、何れの現在のNMR実験とも同様に、アクティブに信号に注入される情報は、経時的な磁化の位相に制限される。
磁化は、経時的に、位相の変動を受け、更に、速度及び加速度の変動も受けるが、測定される唯一の成分は位相である。測定された磁化それぞれは、(患者などの)撮像する対象のオブジェクトの点に対応する。
S(k,k)として事前に書き換えられる信号に対して行われるデュアル・フーリエ変換は、非常に容易に前述の系を解くことができる数学的なトリックであり(、又は上記トリックに過ぎず)、よって、S(k,k)からρ(x,y)に転送することが可能になるからである。
取得された信号の縮退は、例えば、(画素のサイズの程度の)非常に小さな残留エラー・レートでかなり削減される。極めてすばやい変動が、正の最大振幅と負の最大振幅との間のxに沿ったグラディエントに課される(図1)。(1kHzに近い周波数で+1テスラ/mという値と−1テスラ/mという値との間で切り換わらなければならないことがあり得る。)。
上記動作手順の限界は、ハイゼンベルグ方程式(不確定性原理)によって支配される。
空間(x,y,z,v)における空間解決は、共役空間(k,k,k,k)における空間分解能に反比例する。
空間分解能及び時間分解能は、(患者の安全基準によって課される制約、ハードウェアによってもたらされる制限、使用されるマシンに固有のパラメータである)信号を符号化するために印加され得る場のグラディエントの強度、及び(各組織に固有のパラメータである)信号が消滅する前の利用可能な時間に依存する。
よって、k空間を走査するために、極めてすばやい変動が、xに沿ったグラディエントに対して課される。グラディエントは、1kHz程度の周波数における(メートル毎の1テスラの程度の)その限度値間で振動しなければならないことがあり得る。
前述の現象は、強度が非常に大きな音響雑音をもたらし、これは、NMR実験中に、フィ―ルド内に配置された患者の鼓膜の回復不能な損傷をもたらし得る。
更に、グラディエント強度が大きい場合、前述の振動は、有害な直接的な神経刺激(指に入った、痛ましいピンや針)をもたらす。
更に、信号を取得することができる時間は一般に、2百又は3百ミリ秒に制限され、よって、観測可能な信号を生成するために使用される無線周波数パルスは、画像を取得するために反復されなければならない。前述の反復は、組織の局所的な加熱(体にとって危険であり得る比吸収率(SAR))を引き起こす。
最後に、体の特定の組織(軟骨、腱、骨等)は、短い緩和時間によって特徴付けられる。
前述の組織の場合、信号は、(骨の場合、10ミリ秒程度の)非常に短い期間にわたって取得することが可能であるに過ぎない。前述の限定性は、達成することが可能な最大空間分解能を大きく制限するか、又は、組織(すなわち、骨)を見ることを妨げる。
多くの欠点を犠牲にして、NMRでは、k空間を走査することにより、取得された信号の縮退を削減することが可能であり、前述の走査は、空間を走査する軌跡及び複素数を施し、前述の欠点をもたらす。
取得された信号は、複素数の系列であり、信号の各点は振幅及び位相を表し、それ以上を表すことができない。実際に、2つの情報(その実数部分において1つ及びその虚数部分において1つ、又は、更に、関心がある、より特殊な場合では、その振幅において1つ、及びその位相において別のもののみ)を収容することができるに過ぎない。
よって、例えば、同じ位相によって特徴付けられる2点は、再構成後、同じ点となる。
適切な点が、大きく異なる速度又は加速度によって更に特徴付けられるものであっても、再構成された信号は次いで、縮退される。
本発明は、信号の縮退を削減するための新たな手法を提案している。
この目的で、本発明は、複素数の系列のものよりも「優れた」形態を使用することにより、取得された信号の縮退を削減することが可能である。
信号を規定するために、位相情報のみよりもずっと多くの情報を注入することができる。
本発明では、超複素数の項は、四元数及び八元数、並びに、クリフォード代数によって定義される要素、及び拡張されるか、又は複素代数を超える代数の要素を規定するために使用される。
第1の局面によれば、本発明は、複素信号を処理する方法であって、複素数として信号を取得する工程と、関連付けられた超複素成分を取得複素信号から求める工程であって、上記成分が、取得複素信号の位相の少なくとも時間導関数に対応する工程と、処理から生じる信号が取得信号の成分の数よりも大きな成分数を含むように、超複素信号に対する処理動作が求められる工程とを含む方法に関する。
方法の他の局面は以下の通りである。
上記方法は、求められた空間において少なくとも2つの次元でオブジェクトを、取得前に符号化する工程を含み、
空間を符号化するために使用される別々の導関数の線形結合によって形成される多項式が、符号化空間次元に関連付けられ、多項式は、直交である2つの直交次元に関連付けられる。
符号化は、信号の各点の時間導関数の形式で求められる空間の信号の各点の情報の特徴的情報を作成し、よって、符号化空間の各点には、超複素数が関連付けられ、その成分は、この点を空間的に特徴付ける別々の時間導関数の値である。
取得信号の各点には、超複素数が関連付けられ、その成分は、空間を符号化するために使用される信号の別々の時間導関数である。
更に、超複素空間において解くことにより、点の空間密度を信号から求める工程を含む復号化工程を含み、一方で、線形等式系が、取得信号によって形成される超複素ベクトルによって形成され、他方で、上述の空間の点を符号化するための超複素行列によって形成される。
方法は、核磁気共鳴システムにおいて適用され、処理する対象の信号は、核磁気共鳴信号である。
符号化は、磁場グラディエントによって行われる。
第2の局面によれば、本発明は、核磁気共鳴撮像システムに関する。
本発明の装置は、本発明の第1の局面による方法を適用する手段を含む。
システムの縮退を削減するための超複素数の使用は、既に、核磁気共鳴の範囲において開示されている(例えば、
E. Kupceらによる「Projection reconstruction
technique for speeding up multidimensional NMR spectroscopy, (J. Am. Chem. Soc., Vol. 126, 2004,p. 6429−6440)」、
Kazimierczukらによる「Two dimensional Fourier transform of arbitrary sample, NMR data sets、JMR 179, 2006, p. 323−328」、
Chen Jらによる「Ultra High resolution 3D NMR
spectra from limited−size data sets, JMR, No. 169, 2004, p. 215−224」、
Granwherらによる「Sensitivity quantification of remote detection NMR and MRI」)
上記ケース全てにおいて、同じ情報(位相、超複素数の各成分は、特定の時点におけるか(、又は、上記Granwehrの場合、特定の「飛行時間」後の)信号の位相を含む)を記憶するために超複素数が使用される。
本願の方法では、取得された超複素信号の各成分が、固有に、符号化空間次元に沿った別の性質(位相、速度、加速度、角加速度等)のものであること(、及び、上記刊行物に開示されている、何れかの次元のタイプ(位相、位相、位相等)でないこと)であるということにより、磁化の空間位置の縮退の問題の、他の手法全てと比較した、より大きな削減を可能にする。
超複素数により、信号を処理することにより、k空間という概念を捨てることが可能である。
従前の研究の全てにおいて、かつ、現在の研究において、フーリエ変換をもう使用していなくても、逆格子空間(k空間)は、しかしながら使用される。従前の研究では全て、軌跡により、逆格子空間において、関心の平面又はボリュームが走査されるということが課される。このボリュームは、符号化する(すなわち、3次元オブジェクト(x,y,z)を符号化する)対象のボリュームと同じ次元数を有するので、3次元逆格子空間(kx,ky,kz)のボリュームを現時に取得しなければならない。
オブジェクトを符号化するための超複素数及びモーメントが使用された場合、空間情報の符号化は、各点で実現され、符号化する対象のオブジェクトと同じ次元数のボリュームを走査する。前述の新たな手法は、よって、オブジェクトを符号化するために使用される軌跡の単純化を大いに可能にする。
従来の手法では、獲得された軌跡の各点は、複素数によって規定され、したがって、2次元(実数次元及び虚数次元(一方は複素平面を表す)を有する。
前述の軌跡の各点は、特定の時点における信号の振幅及び位相を特徴付ける。
ボリュームを再構成することができるためには、取得信号の経時的な軌跡は、符号化する対象のオブジェクトと同じ次元数Nの空間においてボリュームをたどるべきである。
したがって、N個の次元でオブジェクトを符号化するために使用される独立した次元の合計数は、2Nであり、これは、N個の次元で、空間に配置された複素数の次元数を表す。
本願提案の方法では、取得点はそれぞれ、次元数(2N)の超複素数であり、その成分対はもう、磁化の振幅及び位相のみを表すものでなく、振幅及び角速度、振幅及び角加速度等も表す。
取得点の次元数が既に2Nであり、各成分に含まれる情報が別々の次元数のものであるので、次元数Nのボリュームを再構成することができる軌跡の最小次元数は再び1に低下する。
よって、本発明では、
(超複素形式において処理される)信号に符号化中に情報がどのようにして注入されるかにおける柔軟性を大きく増加させ、
音響雑音を削減するか、又は完全に抑制し、NMR画像取得中に、特定の場合において、120dB超から20dB未満の値まで通過させ、
患者における直接的な神経刺激を削減するか、又は抑制し、
患者の組織の加熱を削減し、
NMR画像の時間分解能を改善し、
NMR画像の空間分解能を改善する。
前述の範囲内で、全体的に、強度の大きなグラディエントを使用することなく、非常にすばやく切り換える必要なしで、画像又はボリュームの全体を符号化することが可能である。
よって、本発明では、ハイゼンブルグ方程式によってもう制限されない範囲内でこれを、超複素数を使用して、特に、画像又はボリュームを取得し、再構成することが可能である。
特に、取得時間と空間分解能との比は、既知の手法と比較して効果的である。
更に、本発明の方法は、既知のタイプのNMR装置に特にうまく一体化される。
2D画像の場合における「エコー・プラナー・イメージング」シーケンスの範囲内のk空間の走査を示す図である。 2D画像の場合における「エコー・プラナー・イメージング」シーケンスの範囲内のk空間の走査を示す図である。 本発明の方法を示すフロー図である。 本発明の方法において使用される、k空間内の軌跡(k空間内の放物線軌跡)を示す図である。 本発明の方法において使用される、k空間内の軌跡(k空間内の放物線軌跡)を示す図である。 (システムである)符号化される対象の画像を示す図である。 本発明の方法中に使用される軌跡をk空間において示す図である。 既知であり、本発明の方法によって得られるものとして、フーリエ変換を使用することによる、復号化後の再構成された画像を示す図である。 既知であり、本発明の方法によって得られるものとして、フーリエ変換を使用することによる、復号化後の再構成された画像を示す図である。
本発明の他の構成及び効果は、単に例証的であり、非限定的であり、かつ、添付図面を参照して読まれるものとする以下の説明から更に明らかになるものである。
図2は、4次元の超複素数(すなわち、四元数)を使用してNMR信号を取得する方法のフロー図を示す。
NMR適用例の範囲内で、符号化する対象のオブジェクトが、磁場内に配置される。
このケースは、2次元により、オブジェクトを符号化したい場合である。
このオブジェクトを時間信号で符号化するために、無線周波数信号(RF)E1が印加され、これは、観測可能な磁化をもたらし、2次元空間を走査するために使用され、E2が2次元に沿って印加され、磁場グラディエントは次いで、その位置に応じて小さな差分要素それぞれの磁化の特性(位相、回転速度、加速度)を調節する。
図3bは、図3aのグラディエントを印加することによる、k空間の走査を示す。方向xに沿って、一定の磁場グラディエントが印加される。前述の次元に沿って、磁化の回転周波数が、その位置の線形関数になる。方向yに沿って、磁場グラディエントがランプ(勾配)に沿って印加される。この次元に沿って、磁化の回転速度の加速度がその位置の線形関数になる。したがって、k空間の概念を使用した従来の表現では、符号化する対象の空間を包含する軌跡は放物線である(図3b参照)。
特に、一定の磁場グラディエントは、第1に方向に沿ったその位置に応じて磁化の回転速度に関する情報を信号に注入することを可能にし、ランプ形状の磁場グラディエントは、第2の方向に沿ったその位置に応じて磁化の回転速度の加速度に関する情報を信号に注入することを可能にする。画像の再構成を可能にするうえでは、前述の単一の放物線軌跡を取得することで十分である。
従来の手法と比較すれば、符号化する対象の空間を包含するために使用される軌跡はずっと単純である。
更に、使用される磁場グラディエント、強度の変動が非常に低く、それにより、信号取得処理によって生成される雑音、及びその生成に必要なエネルギが非常に大きく削減される。
更に、取得時間は従来の手法よりも短いことがあり得、それにより、体に有害な神経刺激及び組織の加熱を避けることが可能である。
取得信号E4は、超複素形式で使用される。
実際には、差分回路と結合されたアンテナを使用することにより、又はさもなければ、差分信号を直接供給するアンテナを使用することによって取得することができる。
この例では、4つの成分を備えた超複素信号が使用され、この場合、これは四元数である。
四元数は4つの成分(rと表す実数成分、iと表す虚数成分、並びに、j及びkと表す他の2つの成分)を含む。四元数信号の4つの成分は、
=j=k=i.j.k=−1
という関係を介して互いに関係付けられる。
四元数としての取得信号は、
Figure 2010538718
という式で表される。ここで、S(t)は、(既知の手法におけるように)信号の複素成分を表し、S(t)は、S(t)の時間導関数を表す。演算子R及びIはそれぞれ、信号の実数部分及び虚数部分を呈する。前述の取得信号の式では、関数f(t)は、符号化する対象の空間を包含するために使用される軌跡を表す。
この目的で、使用される多項式軌跡は例えば、
Figure 2010538718
として表される。ここで、指数nは、符号化される対象の次元を表す。
他の軌跡を使用することができ、使用される軌跡は、実際に、信号の縮退を「十分に」削減することができる(すなわち、包含された軌跡は、解く対象の等式系に、十分な情報をもたらし、よって、等式系を解くことにより、所望の局所化に有用な情報を得ることができる(関心のあるオブジェクトの外部で上記情報が完全に縮退されているか否かは当然、意味を持たない)。
直交方向に対応する、対応する多項式自体が直交であるという唯一の条件下で、軌跡は例えば、それ自体が多項式であり得る。
特に、前述の例では、取得信号は、印加された磁場グラディエントに応じて、
Figure 2010538718
として表され、因数2π及びγは、Hz、Hz/T、T/m、T/m/sの間の次元数の理由で導入される。
(t)をGtで置き換えることにより、前述の式は、
Figure 2010538718
又は、更に、
Figure 2010538718
になる。
このことを図3bに示しているように、k空間を包含するために使用される放物線軌跡は、k空間の一様な走査を可能にせず、達しない点はこの場合、ゼロである。
上述のように、取得信号は、S(t)=Mρ(x,y)として表すことができ、ここで、S(t)及びρ(x,y)はベクトルであり、Mは行列である。k空間の概念を使用して取得信号Sから密度関数ρに渡そうとした場合、信号に注入した情報は回復することも使用することも可能でなく、再構成により、ρについて完全に「異常な」値がもたらされる。
特に、示された例では、信号の位相のみがk空間/フーリエ変換手法によって使用され、画像を再構成するための唯一のやり方は、信号内に注意深く符号化した、速度に関する情報を使用するということである。補償として、前述の手法は、超複素線形等式系の直接の解決を必要とする。
図4aは、符号化されることが意図されたオブジェクトを示す。
図4bは、(放物線軌跡である)方向yに沿ったランプ及び方向xに沿った一定のグラディエントを印加することによって取得される信号の実数部分を示す。
(非ゼロ)情報を含む唯一の点は、放物線軌跡上に配置される点である。
図5aは、フーリエ変換を使用することによって得られる画像を示す。この例では、フーリエ変換を使用するための条件により、情報を平面上に分布させることが必要となる。取得情報をこの形式にすることが可能でないということにより、適切に使用されることが妨げられ、再構成画像は、決して、元の画像に匹敵するものでない。この場合、縮退状態を超える状態の再構成画像は、符号化画像と全く異なる。
図5bは、超複素線形等式系を解いた上述の四元数として取得信号を事前に記憶した後の復号化画像を示す。
この手法では、従来の手法の欠点なしで、縮退を削減し、元のものと同一の画像を再構成することが可能である。
なお多くの情報を記憶するためには、より高次のモーメント(加速度の導関数等)を使用することで十分である。

Claims (9)

  1. 複素信号を処理する方法であって、
    信号を複素数として取得する工程と、
    関連付けられた超複素成分を、取得された複素信号から求める工程であって、前記成分が、少なくとも、前記取得された複素信号の位相の時間導関数に対応する工程と、
    それによって求められた超複素信号を処理する工程であって、前記処理から結果として生じる信号が、前記取得信号の成分の数よりも大きな成分数を含む工程とを含む方法。
  2. 請求項1記載の方法であって、前記取得前に、求められた空間における少なくとも2つの次元でオブジェクトを符号化する工程を含む方法。
  3. 請求項2記載の方法であって、符号化された空間次元それぞれには、前記空間を使用するために使用される別々の導関数の線形結合によって形成される多項式が関連付けられ、前記多項式は、直交である2つの直交次元に関連付けられる方法。
  4. 請求項2乃至3のうちの何れか一項に記載の方法であって、前記符号化する工程は、前記信号の各点の時間導関数として求められる、前記空間の前記信号の点毎に、特徴付ける情報を生成し、よって、前記符号化空間の各点には超複素数が関連付けられ、その成分は、前記点を空間的に特徴付ける別の時間導関数の値である方法。
  5. 請求項1乃至4のうちの何れか一項に記載の方法であって、前記取得信号の各点には、超複素数が関連付けられ、その成分は、空間を符号化するために使用される、別々の時間導関数の値である方法。
  6. 請求項1乃至5のうちの何れか一項に記載の方法であって、他方で、請求項4記載の空間の点を符号化するための超複素行列により、一方で、請求項5記載の取得信号に形成された超複素ベクトルによって形成される線形等式系を超複素空間において解くことにより、点の空間密度を信号から求める工程を含む復号化工程を更に含む方法。
  7. 請求項1乃至6のうちの何れか一項に記載の方法であって、核磁気共鳴システムにおいて適用され、処理する対象の信号が核磁気共鳴信号である方法。
  8. 請求項4記載の方法であって、前記符号化が、磁場グラディエントによって行われる方法。
  9. 核磁気共鳴撮像システムであって、請求項1乃至8のうちの何れか一項に記載の方法を適用する手段を備える核磁気共鳴撮像システム。











JP2010524529A 2007-09-18 2008-09-18 超複素数を使用したmri空間符号化 Pending JP2010538718A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR0757654A FR2921160B1 (fr) 2007-09-18 2007-09-18 Procede de traitement de signaux hypercomplexes
PCT/EP2008/062480 WO2009037326A1 (fr) 2007-09-18 2008-09-18 Encodage spatial en irm au moyen de nombres hypercomplexes

Publications (1)

Publication Number Publication Date
JP2010538718A true JP2010538718A (ja) 2010-12-16

Family

ID=39311077

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2010524529A Pending JP2010538718A (ja) 2007-09-18 2008-09-18 超複素数を使用したmri空間符号化

Country Status (6)

Country Link
US (1) US8558544B2 (ja)
EP (1) EP2191288A1 (ja)
JP (1) JP2010538718A (ja)
CA (1) CA2699749C (ja)
FR (1) FR2921160B1 (ja)
WO (1) WO2009037326A1 (ja)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9557397B2 (en) * 2013-11-04 2017-01-31 Aspect Imaging Ltd. Method for manipulating the MRI's protocol of pulse-sequences
CN104200421B (zh) * 2014-08-01 2017-07-28 南京信息工程大学 基于四元数正交变换的彩色图像加密方法及解密方法
CN108122233B (zh) * 2017-12-18 2021-11-19 辽宁师范大学 基于局部像素综合特征的彩色图像分割方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02232034A (ja) * 1979-08-20 1990-09-14 General Electric Co <Ge> 核磁気共鳴装置

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5998996A (en) * 1997-03-27 1999-12-07 General Electric Company Correction of artifacts caused by Maxwell terms in phase contrast angiography
US20040217760A1 (en) * 2000-02-11 2004-11-04 Madarasz Frank L. Bayesian methods for flow parameter estimates in magnetic resonance imaging
WO2005004703A2 (en) * 2003-06-30 2005-01-20 Board Of Regents, The University Of Texas System Methods and apparatuses for fast chemical shift magnetic resonance imaging
DE602005021326D1 (de) * 2005-06-04 2010-07-01 Bruker Biospin Ag Automatische Projektions-Spektroskopie
US8200025B2 (en) * 2007-12-07 2012-06-12 University Of Ottawa Image classification and search

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02232034A (ja) * 1979-08-20 1990-09-14 General Electric Co <Ge> 核磁気共鳴装置

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
JPN5010013293; GRANWEHR: '"Sensitivity quantification of remote detection NMR and MRI"' JOURNAL OF MAGNETIC RESONANCE V179 N2, 200604, P280-P289, ACADEMIC PRESS *
JPN5010013294; CHEN J: '"Ultra-high resolution 3D NMR spectra from limited-size data sets"' JOURNAL OF MAGNETIC RESONANCE V169 N2, 200408, P215-P224, ACADEMIC PRESS *
JPN5010013295; KAZIMIERCZUK: '"Two-dimensional Fourier transform of arbitrarily sampled NMR data sets"' JOURNAL OF MAGNETIC RESONANCE V179 N2, 200604, P323-P328, ACADEMIC PRESS *
JPN5010013297; E. KUPCE: '"Projection-Reconstruction Technique for Speeding up Multidimensional NMR Spectroscopy"' J.AM.CHEM.SOC. V126 N20, 2004, P6429-P6440 *

Also Published As

Publication number Publication date
WO2009037326A1 (fr) 2009-03-26
EP2191288A1 (fr) 2010-06-02
CA2699749C (fr) 2017-05-23
CA2699749A1 (fr) 2009-03-26
FR2921160B1 (fr) 2009-11-20
US20100213936A1 (en) 2010-08-26
FR2921160A1 (fr) 2009-03-20
US8558544B2 (en) 2013-10-15

Similar Documents

Publication Publication Date Title
JP6513398B2 (ja) 事前情報に制約される正則化を使用するmr画像再構成
Mickevicius et al. Investigation of undersampling and reconstruction algorithm dependence on respiratory correlated 4D-MRI for online MR-guided radiation therapy
US8680861B1 (en) System and method for multislice fast magnetic resonance imaging
US9478051B2 (en) Method for accelerating magnetic resonance imaging using varying k-space sampling density and phase-constrained reconstruction
WO2018010134A1 (en) System and method for magnetic resonance imaging
US9103898B2 (en) Motion correction in accelerated T1-weighted magnetic resonance imaging
US20090009167A1 (en) Magnetic Resonance Spectroscopic Imaging with Short Echo Times
JP6332891B2 (ja) サイドバンドアーチファクトを抑制する並列マルチスライスmr撮像
JP2016512780A (ja) 圧縮センシングを使用したmr画像の再構成
Witschey et al. Localization by nonlinear phase preparation and k‐space trajectory design
CN106574954A (zh) 针对epi的具有奈奎斯特伪影校正的并行mr成像
US10203387B2 (en) MR imaging with enhanced susceptibility contrast
Gao et al. PCLR: phase‐constrained low‐rank model for compressive diffusion‐weighted MRI
US10698053B2 (en) System and method for gradient-modulated sweep imaging with fourier transformation magnetic resonance imaging
JPWO2012077543A1 (ja) 磁気共鳴イメージング装置及びコントラスト強調画像取得方法
US9274197B2 (en) Magnetic resonance imaging data sampling methods and systems
JP2020527089A (ja) ディクソン型水/脂肪分離mrイメージング
Ilbey et al. Single point imaging with radial acquisition and compressed sensing
US9535148B2 (en) Dynamic contrast enhanced magnetic resonance imaging with high spatial-temporal resolution
JP2010538718A (ja) 超複素数を使用したmri空間符号化
JP2019503750A (ja) 動きに起因する拡散勾配不整合の補正を備えるdti
KR101165434B1 (ko) 키홀 방법을 이용한 자기공명 온도 감시 영상 획득 방법
CN112334786A (zh) 用于螺旋采集的mr体模
US20230400545A1 (en) Spin echo mr imaging with spiral acquisition
CN113544526B (zh) 使用波编码的并行mr成像

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20110822

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130820

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20130821

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20140212