JP5827989B2 - 放射状mriデータに対するオシログラム・フィルタリング方法 - Google Patents

放射状mriデータに対するオシログラム・フィルタリング方法 Download PDF

Info

Publication number
JP5827989B2
JP5827989B2 JP2013502845A JP2013502845A JP5827989B2 JP 5827989 B2 JP5827989 B2 JP 5827989B2 JP 2013502845 A JP2013502845 A JP 2013502845A JP 2013502845 A JP2013502845 A JP 2013502845A JP 5827989 B2 JP5827989 B2 JP 5827989B2
Authority
JP
Japan
Prior art keywords
oscillogram
sinogram
projection
data
axis
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2013502845A
Other languages
English (en)
Other versions
JP2013523280A (ja
Inventor
エー. コルム,カーティス
エー. コルム,カーティス
イディヤチューリン,ジャウダットエス.
エス. イディヤチューリン,ジャウダット
ムーラー,ステーン
ジー. ガーウッド,マイケル
ジー. ガーウッド,マイケル
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Minnesota
Original Assignee
University of Minnesota
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 University of Minnesota filed Critical University of Minnesota
Publication of JP2013523280A publication Critical patent/JP2013523280A/ja
Application granted granted Critical
Publication of JP5827989B2 publication Critical patent/JP5827989B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/4816NMR imaging of samples with ultrashort relaxation times such as solid samples, e.g. MRI using ultrashort TE [UTE], single point imaging, constant time imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7225Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/725Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7253Details of waveform analysis characterised by using transforms
    • A61B5/7257Details of waveform analysis characterised by using transforms using Fourier transforms
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271Specific aspects of physiological measurement analysis
    • A61B5/7278Artificial waveform generation or derivation, e.g. synthesising signals from measured signals
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/74Details of notification to user or communication with user or patient ; user input means
    • A61B5/742Details of notification to user or communication with user or patient ; user input means using visual displays
    • 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/4818MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • G01R33/4824MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory
    • 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/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • 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/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56509Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling
    • 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/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/5659Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the RF magnetic field, e.g. spatial inhomogeneities of the RF magnetic field

Landscapes

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

Description

優先権の主張
本特許出願は、米国特許法第119条(e)に基づいて、出願日2010年3月31日、Curtis A. Corumの米国特許仮出願第61/319,703号、名称「METHOD OF OSCILLOGRAM FILTERING FOR RADIAL MRI DATA」(代理人整理番号、600.746PRV)に対する優先権の利益を主張し、また、同仮出願全体を参照により本明細書に組み込む。
米国政府支援による研究または開発についての申告
本発明は、米国国立衛生研究所、NIH BTRC Core3からの授与番号5P41RR008079に基づく政府支援を受けてなされた。米国政府は、本発明の一定の権利を有する。
磁気共鳴画像中のアーチファクトには、断続的で広帯域のノイズが含まれることがある。現在の技術は、こうしたアーチファクトを軽減するのには不十分である。
図面は、必ずしも原寸に比例して描かれておらず、異なる図の中の同様のコンポーネントを同様の数字で説明していることがある。異なる添え字を持つ同様の数字が、同様のコンポーネントの異なる例を表すことがある。図面は、本文書で議論する様々な実施形態を限定ではなく例として概ね例示している。
ソノグラムおよび対応するオシログラムの図である。 オシログラムおよび対応する投影の図である。 オシログラムおよび対応する投影の図である。 1つの例による方法の図である。 1つの例によるシステムの図である。
以下の例は、任意の順序または組合せで組み合わせることができる。
本発明による主題には、放射状MRIデータ(radial MRI data)から断続的で広帯域のノイズをフィルタリングする方法が含まれる。オシログラム・ドメイン中に示す際の放射状MRIデータのまばらさ、凸性により、画像データへの影響を最小限に抑えながら、フィルタリングを行って、ノイズ(または他の人為現象的信号)を分離および抑制することが可能になる。
RF干渉およびジッパーアーチファクト
無線周波数(RF)アーチファクトが、デカルト・サンプルMRI画像中に「ジッパー(zipper)」として頻繁に現れる。結果として得られる画像内で、ジッパーアーチファクトは一般に、読出し方向には狭くなり、位相符号化方向に沿って延びる。
放射状(投影(プロジェクション))イメージングにおいては、ジッパーは、最終画像中の特定の線(または面)に局在せず、中間のサイノグラム内に現れる。図1は、視野角について一定であるアーチファクトAおよび視野角について正弦的に現れたアーチファクトBを有するサイノグラム(スパイラログラム(spiralogram)またはスパイログラムと称することもある)を示す。
2次元または3次元放射状データセット内には、ジッパーは、視野順序(view ordering)ならびにジッパー干渉の周波数、帯域、位相に依存した構造である「ファン(fan)」アーチファクトまたは「タービン」アーチファクトとして現れることがある。
概ね垂直の(または変動しない)ジッパーアーチファクトに加えて、ジッパーが移動中であること、または時間的に変動する周波数オフセットを有することもある。このことは、放射状取得の視界がアイソセンタに対してオフセットである場合に起こる。この場合、RF干渉は、実際には固定周波数であってもよいが、MRIデータは、非アイソセンタ取得の周波数オフセット・テーブルに起因してオフセットである。
移動ジッパーは、一般に、アイソセンタからの距離および視野順序に依存した正弦的経路を描くこととなる。図示の例では、ジッパーのベースは、r-インデックスBで始まる。
他のフィルタ法
メジアン・フィルタまたは他のフィルタを、(サーチ・アルゴリズムが時間を浪費することがあるが)垂直ジッパーまたは移動ジッパーに適用することもできる。
オシログラム
オシログラムは、サイノグラムの図を変換したものを表す。フーリエ変換が、視野方向に沿って、また、サイノグラムの読出し(r-インデックス)方向に沿って行われる。サイノグラム表現では、ジッパーは、移動するまたは変動しない概ね垂直のフィーチャとして現れ、オシログラム表現では、ジッパーは、水平線として現れる。
本発明の主題の一例によれば、オシログラムの振幅をフィルタリングに使用することができる。
オシログラム・フィルタ
本発明の主題の一例は、オシログラム内のRFジッパーアーチファクトに対応する水平線を自動的に認識し、排除するように構成される。1つの方法は、水平軸全体にわたって積分し、1次元投影を作成することを伴う。
次いで、オシログラムの水平線は、投影上で背景の上方に延びるピークに変換される。
実施例
差分画像を作成して、タービン・アーチファクトの除去および所望の画像データの維持を示すことができる。
最終画像内では、タービン・アーチファクトが、復元プロセスによりタービンと同様の放射状の対称形パターンに散らばるジッパーとして現れる。この構造は、ジッパー干渉の特性ならびに放射状の視野順序に依存する。
他の応用例
周期的またはほぼ周期的な動きもオシログラムから検出可能である。これには、拍動流、心臓の動き、呼吸運動が含まれる。オシログラムまたはオシログラムの部分集合から、ナビゲータ型の情報を作成することもできる。
図1は、SWIFTシーケンスを使用することにより作成したサイノグラム内のジッパーアーチファクトを示す。
SWIFTは、空間的情報を符号化するのに使用される磁場勾配の新たな値(または向き)を少なくとも設定するのに必要な時間だけパルス長Tを上回る短期の繰返し時間Tを有する、周波数変調されたパルスのシーケンスを伴う。画像は、3次元バックプロジェクション・リコンストラクションを使用することにより処理する。一例では、双曲正割ファミリー(hyperbolic secant family)からの周波数変調されたパルス(HSnパルス)を使用する。SWIFTでは、形成されたパルスが、時間依存の振幅および位相を有するN個の異なるサブパルス要素を含むことがある。FMパルス中は、瞬時の共鳴に至るまで、等色線は、有効RFフィールド・ベクトルに追従する。共鳴時には、等色線は、RFパルスの「ハグ(hug)」から解放され、次いで、衰退していく小さな変調を伴ったほぼ自由な歳差運動を示し、これによりスペクトル・コンタミネーションがもたらされる。したがって、こうしたスピン系反応からスペクトル情報を抽出するために、ストカスティック(確率論的)NMRにおいて位相情報を回復する方法と同様の相互相関手法を使用する処理を行う。TR<<T1のSWIFTについて単位時間毎に理論的に達成可能な信号対雑音比(SNR)は、パルスFT(pulsed FT)についてのものと同一である。SWIFT取得中、適用される画像化勾配は、磁化率または不均質性に起因して、本来備わる固有の勾配の全てを通常は超えるものである。この条件の場合、取得される画像は、横緩和に全く影響されず、信号の強さは、T1およびスピン密度にのみ依存する。最大Tコントラストは、有効フリップ角に依存し、感度とコントラストとの最良の妥協点は、エルンスト角の2倍を越えるフリップ角を有することとなる。フリップ角が非常に小さい場合、Tコントラストは、さほど重要でなく、コントラストは、スピン濃度に完全に由来する。画像取得に先立つ、または、画像取得の途中に行われる適当な準備シーケンスによって、他の種類のコントラストを得ることもできる。
SWIFTは、放射状イメージング技法の一例にすぎず、他の放射状イメージング技術を、本発明の主題とともに使用することもできる。
図1では、r-インデックスAにおいて、高強度のジッパーが、垂直になり、中心の左側に位置する。また、第2のジッパーをr-インデックスBに示す。
図1の左側区画に示したサイノグラム・データは、垂直軸上に視野角を示し、水平軸上に読出し方向を示す。視野角データは順次作り出されるので、視野角を時間軸とみなすこともできる。サイノグラム・データは、2次元または3次元の放射状イメージング・プロトコルで作り出すことができる。
図1に示すサイノグラムは、SWIFTイメージング・プロトコルによって作られる自由誘導減衰(FID)として見ることができる。正弦波帯域は、検査中の対象のマージンを表す。
図示のジッパーアーチファクトは、様々な周波源を有することがある。たとえば、共鳴周波数における干渉が、こうしたアーチファクトを引き起こすことがある。さらに、エイリアシングされた周波数またはその周波数の調波からの干渉が、ジッパーアーチファクトにつながることもある。さらに、ジッパーは、システム内のRF源またはそれに近い周波源により引き起こされることもある。ジッパーアーチファクトは、オフセット・テーブルに関連することもある。
本発明の主題は、フィルタを使用することによりサイノグラムを調整する方法およびシステムを含む。フィルタには、メジアン・フィルタ、ノッチ・フィルタ、または、アーチファクトを除去するようになされた他の種類のフィルタが含まれることがある。メジアン・フィルタは、1次元または複数次元で機能することができ、入力の中央値の出力表示を提供する。メジアン・フィルタは、歪みを低減し、平滑化機能を提供することができる。
SWIFTプロトコルを使用するイメージングは、相対的に高速の切換えおよび1〜4μ秒のオーダの短い持続期間を伴う。こうした短い持続期間により、本明細書で記載するジッパーアーチファクトを作り出すことのある無線周波数干渉源がもたらされることがある。
図1に示すように、サイノグラムをフーリエ変換して、オシログラムを作成する。図1のオシログラムでは、垂直軸(ここでは角周波数を指す)を遅軸(時間)とみなすことができるか、または、これを任意値nでインデックス付けすることができる。オシログラムの水平軸も時間軸であり、速軸を指し、ここでは解放前部分と解放後部分とに分けられている。オシログラムのデータは、視野ごとの順序に生のk空間データを表す。
図1に示すように、サイノグラム・ドメイン中のAにおける干渉が、図のオシログラム中のA’およびA’に変換する。さらに、サイノグラム・ドメイン中のBにおける干渉が、オシログラム中のB’およびB’に変換する。RF干渉Aは、狭帯域もしくは広帯域になること、または、一定周波数を有することがある。RF干渉Bは、移動(moving)周波数に対応する。
画像情報(非移動)は、図に示すようにハッシュ・マークを対角線が横切る三角形領域に限定される。垂直のハッシュ・マークを有する三角形領域で示すように、周期的な動きにより、位置の変わったコンポーネントがもたらされることがある。画像情報の三角形の頂点と周期運動の三角形との垂直距離は、干渉の周波数に相当する。さらに、対象の幅が広いほど、オシログラム表現中の三角形の幅は広くなる。
対象データが視野順序とともに変化するので、本発明の主題により、対象の動きとアーチファクトとを識別することができる。干渉は、位相がコヒーレントであることもある。
サイノグラム内のデータを2つのフーリエ変換によって変換して、オシログラムを作成する。上で述べたように、垂直ジッパーアーチファクトは、時間において相対的に狭く、n軸において、または周波数軸に沿って相対的に広い。
図2は、tfast軸に沿ったオシログラム投影の1次元ノッチ・フィルタを示す。
投影I(A1’)は、オシログラム投影内のアーチファクト・ピークまたは干渉ピークを示し、以下のように表すことができる。
Figure 0005827989
この投影を振幅について積分し、一例では、虚数成分または位相成分について積分することができる。
アーチファクトは、画像情報を越えて延び、自動ピーク発見(検出)アルゴリズムによって識別することができる。図2では、オシログラム投影内の矩形ボックスは位置を示す。
オシログラム内のマスキング
RF干渉の角周波数座標を使用することにより、2次元オシログラム内にマスク-アウト・ストリップ(mask-out strip)を提供する。さらに、オシログラムの解放前部分(左側)を捨てることもできる。マスク-アウトすることは、フィルタリングを行うことまたは0に設定することを含む。
オシログラム投影内でのマスキングに続いて、変換および復元により画像サイノグラムを作成する。
一例では、オシログラム内の中央の水平線に「キープアウト」・ゾーンを設定することができる。このキープアウト・ゾーンは、ユーザ・インタフェースを使用することによりユーザが選択可能な幅を有することがある。一例では、キープアウト・ゾーンは、投影の中心の周囲約10%または5〜20%である。
様々なフィルタを使用することもできる。たとえば、こうしたフィルタを視界の外の画像データに適用する(帯域通過)ことができ、あるいは、こうしたフィルタには、RFノッチ・フィルタが含まれることがある。
さらに、SWIFTを使用する脂肪抑制は、磁化プレパレーション(magnetization preparation)・パルスを含むことがある。磁化プレパレーション・プロトコルでは、複数のノッチ・フィルタを投影上で利用することがある。たとえば、複数のフィルタ、すなわち櫛形フィルタを本発明の主題とともに使用することができる。対象の範囲を超える線形領域を減衰させるために、マスクを使用することもできる。
周期信号に加えて、本発明の主題を使用することにより、(心血流や呼吸等の)動きを修正することもできる。
図3は、オシログラム・ドメイン内で構成され、先の図中に示したようなアーチファクトを除去するようになされたマスク・フィルタを有するオシログラム投影を示す。
図4は、一例による方法400の流れ図を示す。この方法は、放射状に規則化された磁気共鳴データに対応する第1のサイノグラムを受け取ることを410で含むが、この第1のサイノグラムは、磁気共鳴データについての読出し方向の関数として視野角を表す。さらに、本方法は、第1のサイノグラムを変換して、角周波数軸を有するオシログラムを形成することを420で含む。また430で、本方法は、オシログラムの選択された軸に沿って投影を形成することを含むが、この選択された軸は、角周波数軸に直交する。本方法は、投影内のピークをフィルタリングすることを440で含み、このフィルタリングされた投影を変換して第2のサイノグラムを形成することを450で含む。
図5は、一例によるシステム500を示す。このシステムは、データ・レシーバ510、サイノグラム・ジェネレータ520、プロセッサ530、フィルタ・モジュール540、出力モジュール550を含むことができる。これらの要素は、別々のコンポーネントまたはモジュールとすること、あるいは、任意の組合せで一体化することまたは組み合わせることができる。さらに、これらの要素の2つ以上を単一のモジュールとして組み合わせることができる。データ・レシーバ510は、放射状に規則化された磁気共鳴データを受け取るように構成される。サイノグラム・ジェネレータ520は、磁気共鳴データに対する読出し方向の関数としての視野角に対応する第1のサイノグラムを作り出すように構成される。プロセッサ530は、角周波数軸を有するオシログラムを作り出すように構成される。オシログラムは、第1のサイノグラムを変換したものに相当する。フィルタ・モジュール540は、オシログラムの選択された軸に沿って形成された投影を選択的にフィルタリングするように構成される。選択された軸は、角周波数軸に直交する。出力モジュール550は、フィルタリングされた投影を変換したものに相当する第2のサイノグラムを形成するように構成される。システム500は、コンピュータと図示のように構成されたハードウェアとの任意の組合せとして実装することができる。
追記
上の詳細な説明は、その一部分を形成する添付の図面への言及を含んでいる。図面には、本発明を実施することのできる具体的な実施形態が例証として示されている。本明細書では、これらの実施形態を「例」とも称している。こうした例は、図示または説明した要素に加えて他の要素を含むことがある。ただし、本発明者らは、図示または説明した要素のみ提供される例もまた企図している。さらに、本発明者らは、特定の例(もしくはその1つまたは複数の側面)または図示もしくは本明細書で説明した他の例(もしくはその1つまたは複数の側面)について、図示または説明した要素を任意に組み合わせたまたは置き換えた例(またはその1つまたは複数の側面)を使用する例も企図している。
本文書中で参照した全ての刊行物、特許、および特許文献は、個別的に参照により本明細書に組み込むかのごとく、その全体を参照により本明細書に組み込む。本文書と上記のように参照により組み込んだ文献とに矛盾した用法がある場合、組み込んだ参照文献における用法は本文書における用法の補助となるものと考えられたい。すなわち、融和することのできない矛盾の場合、本文書における用法が優先する。
本文書中では、特許文献では一般的であるように、「at least one(少なくとも1)」もしくは「one or more(1または複数)」の使用または何か他の例とは無関係に、「a」または「an」の語を1または複数を含めるように使用している。本文書では、「or」の語を非排他的であることを指すのに用いており、すなわち、特に明記しない限り、「A or B」が「A but not B(AだがBでない)」、「B but not A(BだがAでない)」および「A and B(AおよびB)」を含むように「or」の語を使用している。添付の特許請求の範囲では、「including」および「in which」の語が、ぞれぞれ「comprising」および「wherein」の語の平易な英語の等価物として使用されている。また、添付の特許請求の範囲では、「including」および「comprising」の語は限定的ではなく、すなわち、請求項中でこうした語の後に挙げる要素とは他の要素を含むシステム、デバイス、物品、または方法が、依然としてその請求項の範囲内のものであるとみなされる。さらに、添付の特許請求の範囲では、「first(第1)」、「second(第2)」、「third(第3)」等の語を単なるラベルとして使用しており、これらは、その対象に数的条件を課すことを意図しない。
本明細書に記載する方法の例は、少なくとも部分的には機械またはコンピュータで実施することが可能である。いくつかの例には、上記の例で説明した方法を実施するための電子デバイスを構成するのに使用可能な命令を用いて符号化したコンピュータ読取り可能な媒体または機械読取り可能な媒体が含まれることがある。こうした方法の一実装例には、マイクロ符号、アセンブリ言語符号、高水準言語符号等の符号が含まれることがある。こうした符号には、様々な方法を実行するためのコンピュータ読取り可能な命令が含まれることがある。この符号は、コンピュータ・プログラム製品の一部分を形成することができる。さらに、この符号は、実行中または他の時間に、1つまたは複数の揮発性または不揮発性の有形的コンピュータ読取り可能媒体に有形的に記憶することができる。これらの有形的コンピュータ読取り可能媒体の例には、ハードディスク、リムーバブル磁気ディスク、リムーバブル光ディスク(たとえば、コンパクト・ディスクおよびデジタル・ビデオ・ディスク)、磁気カセット、メモリカードまたはメモリスティック、ランダム・アクセス・メモリ(RAM)、リード・オンリ・メモリ(ROM)等が含まれることがあるが、これらに限定されない。
上の説明は、例証となることを意図しており、限定することを意図していない。たとえば、上述の例(またはそれら例の1つまたは複数の側面)を、互いに組み合わせて用いることができる。上の説明を検討することに応じて当業者等により、他の実施形態が用いられる可能性がある。読者が技術的開示内容の特性を迅速に確かめることが可能になるように、米国特許法施行規則(37C.F.R)1.72(b)に従って要約書を提供している。この要約書は、それが特許請求の範囲または請求項の意味を解釈または制限するのに用いられないという理解の基に提出したものである。また、上の詳細な説明では、様々な特徴をまとめて、開示内容を合理化することができる。このことは、特許請求されていない開示された特徴が請求項には必須であることを意図するものと解釈されるべきではない。反対に、発明主題は、開示された特定の実施形態の特徴全てには存在しないことがある。したがって、ここで、添付の特許請求の範囲を詳細な説明に組み込み、各請求項は、単独の実施形態として独立し、また、これら実施形態を様々な組合せまたは順序で互いに組み合わせることが可能であることが企図される。本発明の範囲は、添付の特許請求の範囲に関連して、こうした特許請求の範囲の権利が与えられる均等物の完全な範囲とともに確定するべきである。

Claims (22)

  1. 放射状に規則化された磁気共鳴データを受け取るように構成され、空間的情報を符合化するのに使用される磁場勾配の新たな値又は向きを少なくとも設定するのに必要な時間だけパルス長を上回る短期の繰り返し時間を有する周波数変調されたパルスを伴うデータ・レシーバと、
    前記磁気共鳴データに対する読出し方向の関数としての視野角に対応する第1のサイノグラムを作り出すように構成されたサイノグラム・ジェネレータと、
    角周波数軸を有し、前記第1のサイノグラムを変換したものに相当するオシログラムを作り出すように構成されたプロセッサと、
    前記オシログラムの選択された軸に沿って形成された投影を選択的にフィルタリングするように構成されたフィルタ・モジュールであって、前記選択された軸が、前記角周波数軸に直交する、フィルタ・モジュールと、
    前記フィルタリングされた投影を変換したものに相当する第2のサイノグラムを形成するように構成された出力モジュールと、
    を備えるシステム。
  2. 前記データ・レシーバが、球状に規則化されたデータを受け取るように構成される、
    請求項1に記載のシステム。
  3. 前記データ・レシーバが、2次元データまたは3次元データのうちの少なくとも1つを受け取るように構成される、
    請求項1または2のうちのいずれか1項に記載のシステム。
  4. 前記プロセッサが、フーリエ変換または逆フーリエ変換のうちの少なくとも1つを実装するように構成される、
    請求項1〜3のうちのいずれか1項に記載のシステム。
  5. 前記プロセッサが、前記視野角に沿って、かつ、前記読出し方向に沿ったフーリエ変換を実装するように構成される、
    請求項1〜4のうちのいずれか1項に記載のシステム。
  6. 前記フィルタ・モジュールが、前記オシログラムの振幅成分を使用することにより前記投影を形成するように構成される、
    請求項1〜5のうちのいずれか1項に記載のシステム。
  7. 前記フィルタ・モジュールが、前記オシログラムの時間軸と一直線にそろった方向に1次元投影を形成するように構成される、
    請求項1〜6のいずれか1項に記載のシステム。
  8. 前記フィルタ・モジュールが、ノッチ・フィルタおよびマスク・フィルタのうちの少なくとも1つを実装するように構成される、
    請求項1〜7のいずれか1項に記載のシステム。
  9. 前記フィルタ・モジュールが、メジアン・フィルタを実装するように構成される、
    請求項1〜8のいずれか1項に記載のシステム。
  10. 前記フィルタ・モジュールが、ピーク検出器を実装することを含む、
    請求項1〜9のいずれか1項に記載のシステム。
  11. 前記出力モジュールが、前記第2のサイノグラムを使用することにより画像を作り出すように構成される、
    請求項1〜10のいずれか1項に記載のシステム。
  12. 空間的情報を符合化するのに使用される磁場勾配の新たな値又は向きを少なくとも設定するのに必要な時間だけパルス長を上回る短期の繰り返し時間を有する周波数変調されたパルスを伴い、放射状に規則化された磁気共鳴データに対応し、前記磁気共鳴データに対する読出し方向の関数としての視野角を表す第1のサイノグラムを受け取ることと、
    角周波数軸を有するオシログラムを形成するために、前記第1のサイノグラムを変換することと、
    前記オシログラムの選択された軸に沿って投影を形成することであって、前記選択された軸が、前記角周波数軸に直交することと、
    前記投影内のピークをフィルタリングすることと、
    第2のサイノグラムを形成するために、前記フィルタリングされた投影を変換することと、
    を含む方法。
  13. 放射状に規則化された磁気共鳴データに対応する前記第1のサイノグラムを受け取ることが、球状に規則化されたデータを受け取ることを含む、
    請求項12に記載の方法。
  14. 放射状に規則化された磁気共鳴データに対応する前記第1のサイノグラムを受け取ることが、2次元データまたは3次元データのうちの少なくとも1つを受け取ることを含む、請求項12または13のうちのいずれか1項に記載の方法。
  15. 変換することが、フーリエ変換または逆フーリエ変換のうちの少なくとも1つを実行することを含む、
    請求項12〜14のうちのいずれか1項に記載の方法。
  16. オシログラムを形成するために前記第1のサイノグラムを変換することが、前記視野角に沿ってフーリエ変換を実行することと、前記読出し方向に沿ってフーリエ変換を実行することとを含む、
    請求項12〜15のうちのいずれか1項に記載の方法。
  17. 前記投影を形成することが、前記オシログラムの振幅成分を使用することにより前記投影を形成することを含む、
    請求項12〜16のうちのいずれか1項に記載の方法。
  18. 前記選択された軸に沿って前記投影を形成することが、前記オシログラムの時間軸と一直線にそろった方向に1次元投影を形成することを含む、
    請求項12〜17のうちのいずれか1項に記載の方法。
  19. フィルタリングすることが、ノッチ・フィルタリングを行うことまたはマスク・フィルタリングを行うことのうちの少なくとも1つを含む、
    請求項12〜18のうちのいずれか1項に記載の方法。
  20. フィルタリングすることが、メジアン・フィルタ・アルゴリズムを実装することを含む、請求項12〜19のうちのいずれか1項に記載の方法。
  21. フィルタリングすることが、ピーク検出アルゴリズムを実装することを含む、
    請求項12〜20のうちのいずれか1項に記載の方法。
  22. 前記第2のサイノグラムを使用することにより画像を生成することを更に含む、
    請求項12〜21のうちのいずれか1項に記載の方法。
JP2013502845A 2010-03-31 2011-03-31 放射状mriデータに対するオシログラム・フィルタリング方法 Active JP5827989B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US31970310P 2010-03-31 2010-03-31
US61/319,703 2010-03-31
PCT/US2011/030728 WO2011123647A1 (en) 2010-03-31 2011-03-31 Method of oscillogram filtering for radial mri data

Publications (2)

Publication Number Publication Date
JP2013523280A JP2013523280A (ja) 2013-06-17
JP5827989B2 true JP5827989B2 (ja) 2015-12-02

Family

ID=44072663

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013502845A Active JP5827989B2 (ja) 2010-03-31 2011-03-31 放射状mriデータに対するオシログラム・フィルタリング方法

Country Status (4)

Country Link
US (1) US10537265B2 (ja)
JP (1) JP5827989B2 (ja)
DE (1) DE112011101171B4 (ja)
WO (1) WO2011123647A1 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5827989B2 (ja) 2010-03-31 2015-12-02 リージェンツ オブ ザ ユニバーシティ オブ ミネソタ 放射状mriデータに対するオシログラム・フィルタリング方法
US9000183B2 (en) * 2012-06-20 2015-04-07 Astrazeneca Ab Cyclohexane-1,2′-indene-1′,2″-imidazol compounds and their use as BACE inhibitors
US10605882B2 (en) * 2014-05-30 2020-03-31 The General Hospital Corporation Systems and methods for removing background phase variations in diffusion-weighted magnetic resonance imaging
CN105748074A (zh) * 2016-04-28 2016-07-13 董印 核磁共振诊断装置

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5865143A (ja) * 1981-10-14 1983-04-18 株式会社東芝 コンピユ−タ・トモグラフイ装置
US5525906A (en) 1994-07-11 1996-06-11 General Electric Detection and elimination of wide bandwidth noise in MRI signals
US6204663B1 (en) 1997-03-26 2001-03-20 Numar Corporation Pulse sequence and method for suppression of magneto-acoustic artifacts in NMR data
US6324242B1 (en) * 1999-04-09 2001-11-27 Arch Development Corporation Fast reconstruction with uniform noise properties in half-scan tomography
US6420873B1 (en) 2001-02-15 2002-07-16 Netcom, Inc. MRI noise detector and suppressor
US7457655B2 (en) * 2004-10-08 2008-11-25 Mayo Foundation For Medical Education And Research Motion correction of magnetic resonance images using moments of spatial projections
US7916144B2 (en) * 2005-07-13 2011-03-29 Siemens Medical Solutions Usa, Inc. High speed image reconstruction for k-space trajectory data using graphic processing unit (GPU)
US7403006B2 (en) 2005-10-11 2008-07-22 Regents Of The University Of Minnesota Frequency swept excitation for magnetic resonance
US7676073B2 (en) * 2006-08-29 2010-03-09 Siemens Medical Solutions Usa, Inc. System and method for reducing circular artifacts in tomographic imaging
CN101573629B (zh) * 2007-01-02 2013-07-10 威斯康星校友研究基金会 利用相衬合成图像进行高度约束反投影重建的对比增强磁共振血管造影
CN101711126B (zh) * 2007-06-14 2012-03-28 株式会社日立医药 磁共振成像装置及由倾斜磁场引起的误差修正方法
CN101470180B (zh) * 2007-12-29 2016-01-20 西门子(中国)有限公司 磁共振成像中失真校准的方法和装置
WO2009126285A1 (en) * 2008-04-10 2009-10-15 Regents Of The University Of Minnesota Rf pulse distortion correction in mri
US8472688B2 (en) * 2008-04-17 2013-06-25 Wisconsin Alumni Research Foundation Method for image reconstruction employing sparsity-constrained iterative correction
US8611981B2 (en) * 2008-08-26 2013-12-17 Bracco Imaging S.P.A. MRI-CEST diagnostic technique based on non-punctual analysis
US9841480B2 (en) 2009-04-02 2017-12-12 Regents Of The University Of Minnesota Localization of magnetic particles by means of swift-MRI
JP5827989B2 (ja) 2010-03-31 2015-12-02 リージェンツ オブ ザ ユニバーシティ オブ ミネソタ 放射状mriデータに対するオシログラム・フィルタリング方法

Also Published As

Publication number Publication date
DE112011101171B4 (de) 2015-05-21
DE112011101171T5 (de) 2013-01-10
JP2013523280A (ja) 2013-06-17
WO2011123647A1 (en) 2011-10-06
US20130085374A1 (en) 2013-04-04
US10537265B2 (en) 2020-01-21

Similar Documents

Publication Publication Date Title
Godenschweger et al. Motion correction in MRI of the brain
US8693761B2 (en) Method of generating nuclear magnetic resonance images using susceptibility weighted imaging and susceptibility mapping (SWIM)
US9797974B2 (en) Nonrigid motion correction in 3D using autofocusing with localized linear translations
US9275294B2 (en) Compressed sensing using regional sparsity
US9830687B2 (en) Image processing device, magnetic resonance imaging apparatus and image processing method
KR101939642B1 (ko) 모션 분해 mri를 위한 시스템 및 방법
US10204412B2 (en) Magnetic resonance imaging apparatus, medical image processing device, and image processing method
KR101359134B1 (ko) 도전율 및 자화율 복원 장치 및 방법과 이에 관한 기록매체
KR20160044079A (ko) 세그먼트 주성분 분석 기법을 이용한 뇌파 노이즈 제거 장치 및 뇌파 노이즈 제거 방법
JP5827989B2 (ja) 放射状mriデータに対するオシログラム・フィルタリング方法
US10605882B2 (en) Systems and methods for removing background phase variations in diffusion-weighted magnetic resonance imaging
Chen et al. T2* phase imaging and processing for brain functional magnetic susceptibility (χ) mapping
CN111065328B (zh) 用于有策略地获取的梯度回波成像的系统和方法
Droigk et al. Direct multi-dimensional Chebyshev polynomial based reconstruction for magnetic particle imaging
Huang et al. Continuous MR bone density measurement using water-and fat-suppressed projection imaging (WASPI) for PET attenuation correction in PET-MR
Krug et al. Filtering the magnetohydrodynamic effect from 12-lead ECG signals using independent component analysis
Wang et al. Artifacts reduction in strain maps of tagged magnetic resonance imaging using harmonic phase
KR20140104630A (ko) 자기 공명(MR; Magnetic Resonance) 영상의 복셀(Voxel)에 대한 MR 스펙트럼을 획득하는 방법 및 장치
KR20190081656A (ko) 자기공명영상의 왜곡 보정을 위한 방법 및 이를 적용한 자기공명영상 장치
JP4163526B2 (ja) Mrタギング画像による変位測定システムおよび方法
Otto et al. Angular smoothing and radial regularization of ODF fields: application on deterministic crossing fiber tractography
da S Senra Filho et al. Anomalous diffusion process applied to magnetic resonance image enhancement
Bartusek et al. Optimization of wavelet-based de-noising in MRI
Callaghan et al. Padé methods for reconstruction and feature extraction in magnetic resonance imaging
US11112480B2 (en) Compressed sensing with regularization parameter

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20140317

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20140317

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20150213

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150310

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20150608

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150619

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20150709

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20150709

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150805

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20150805

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20151006

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20151019

R150 Certificate of patent or registration of utility model

Ref document number: 5827989

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250