JP2022535548A - オフ共鳴アーチファクト補正を伴うスパイラルｍｒ撮像 - Google Patents

Abstract

本発明は、ＭＲデバイス１の検査ボリュームに位置付けられた物体１０のＭＲ撮像の方法に関する。本発明の目的は、強いＢ_０不均一性がある状況においても、効率的で高品質な非デカルトＭＲ撮像を可能とすることである。本発明によると、方法は、少なくとも１つのＲＦ励起パルスと変調された磁場勾配とを含む撮像シーケンスを物体に施すステップと；少なくとも１つの非デカルトｋ空間軌道に沿ってＭＲ信号を取得するステップと；取得されたＭＲ信号からＭＲ画像を再構成するステップと；ＭＲ画像におけるＢ_０不均一性に誘起された不十分なｋ空間サンプリングによって生じる１つ又は複数の不良サンプリングアーチファクトを、深層学習ネットワークを使用して検知するステップとを有する。更には、本発明は、ＭＲデバイス１とコンピュータプログラムとに関する。

Description

本発明は、磁気共鳴（ＭＲ）撮像の分野に関する。本発明は、物体のＭＲ撮像の方法に関連する。本発明は、ＭＲデバイスと、ＭＲデバイス上で実行されるコンピュータプログラムとにも関連する。

２次元的又は３次元的画像を形成するために磁場と核スピンとの間の相互作用を利用する画像形成ＭＲ方法は、今日、広く使用されており、それらは、軟組織の撮像に関して、他の撮像方法よりも多くの観点において優れており、電離放射線を必要とせず、通常は侵襲的でないので、医療診断の分野において顕著である。

ＭＲ方法によると、通常、物体、例えば検査されるべき患者の身体は、強い均一な磁場（Ｂ_０）であって、同時にその方向は測定が基づく座標系の軸（通常はｚ軸）を定める磁場（Ｂ_０）に配置される。磁場は、定められた周波数（いわゆるラーモア周波数又はＭＲ周波数）の交番電磁場（ＲＦ場）の印加によって励起（スピン共鳴）され得た磁場強度に基づいて個々の核スピンについて異なるエネルギーレベルを生む。巨視的な観点から見ると、個々の核スピンの分布は、適切な周波数の電磁パルス（ＲＦパルス）の印加によって平衡状態から外れ得る全体的な磁化を生み、磁化は、ｚ軸の周りで歳差運動を行う。歳差運動は、その開口角度がフリップ角と称されるコーンの面を表す。フリップ角の大きさは、印加された電磁パルスの強度及び持続期間に依存する。いわゆる９０°パルスの場合、スピンはｚ軸から横断面へと外れる（フリップ角９０°）。

ＲＦパルスの停止後に、磁化は元の平衡状態に戻るように緩和され、ｚ方向における磁化は、第１の時間定数Ｔ_１（スピン格子又は長手方向緩和時間）において再び構築され、ｚ方向に直交する方向における磁化は、第２の時間定数Ｔ_２（スピン－スピン又は横断緩和時間）において緩和される。磁化の変化は、ｚ軸に直交する方向において磁化の変化が測定されるようにＭＲデバイスの検査ボリューム内で配置され、方向付けられた受信ＲＦコイルによって検知され得る。例えば９０°パルスの印加後に、横断磁化の減衰は、同一位相の秩序状態から全ての位相角が均一に分布した状態（位相散逸）への（局所的な磁場の不均一性によって誘起された）核スピンの遷移を伴う。位相散逸は、リフォーカスパルス（例えば１８０°パルス）によって補償され得る。このことは、受信コイルにおいてエコー信号（スピンエコー）を生む。

身体における空間解像度を実現するために、３つの主軸に沿って延びる一定の磁場勾配が、均一な磁場に重畳され、これは、スピン共鳴周波数の線形の空間的依存性につながる。すると、受信コイルにおいてピックアップされる信号は、身体におけるの種々の場所に関連付けられ得る種々の周波数の成分を含む。受信コイルを介して得られる信号データは、ｋ空間と呼ばれる空間的周波数ドメインに対応する。ｋ空間データは、通常、異なる位相符号化によって取得された複数のラインを含む。各ラインは、いくつかのサンプルを収集することによってデジタル化される。ｋ空間データのセットは、画像再構成アルゴリズムによってＭＲ画像に変換される。

Ｂ_０不均一性及び組織に誘起される感受性の変化に起因して生じるオフ共鳴効果は、ＭＲ画像において非常に一般的なアーチファクトの主な原因の代表例である。Ｔｈ．ＫｕｅｓｔｎｅｒらによるＩＳＭＲＭ－２０１９のアブストラクト「Ｓｉｍｕｌｔａｎｅｏｕｓ ｄｅｔｅｃｔｉｏｎ ａｎｄ ｉｄｅｎｔｉｆｉｃａｔｉｏｎ ｏｆ ＭＲ ａｒｔｅｆａｃｔ ｔｙｐｅｓ ｉｎ ｗｈｏｌｅ ｂｏｄｙ ｉｍａｇｉｎｇ」、ＩＳＭＲＭ－２０１８（ａｂｓｔｒａｃｔ４３０）は、Ｔ_１重み付き又はＴ_２重み付きＦＳＥシーケンスによって取得された磁気共鳴画像における運動又は場の不均一性から生じるアーチファクトのＣＮＮに基づく検知について言及している。デカルトｋ空間軌道に沿ったＭＲ信号取得では、Ｂ_０不均一性が、画像品質に、一般的には重大な影響を与えない幾何学的歪みを増加させることがよく知られている。しかしながら、スパイラル撮像におけるような非デカルトｋ空間軌道では、オフ共鳴が、通常、理想的な画像のぼけとして現れ、画像の診断上の妥当性に深刻な影響を与え得る。Ｂ_０不均一性によって起こるアーチファクトの範囲は、主磁場の強度とともに拡大又は縮小する。

スパイラル撮像は、効率的なｋ空間カバレッジ及び流動アーチファクトへの低感度から恩恵を受ける迅速なＭＲ撮像技術である。しかしながら、これは、長い取得間隔が使用されるときに（例えば、シングルショットスパイラル撮像において）、Ｂ_０不均一性から特に影響を受けやすい。平行撮像技術と組み合わされたこのような長い取得間隔は、可能な限り高い空間解像度を得るために欠かせないものである。このことは、生理学的な及び患者の運動の問題を緩和するために、機能的ＭＲ撮像（ｆＭＲＩ）において、及び拡散重み付き撮像（ＤＷＩ）においても興味深い。

スパイラルＭＲ撮像のためのぼけ除去方法が当技術分野において知られている。例えば、Ｂ_０マップを取得し、Ｂ_０マップに基づいてＢ_０不均一性効果についてＭＲ信号データを補正することが知られている（例えば、Ａｈｕｎｂａｙらの「Ｒａｐｉｄ ｍｅｔｈｏｄ ｆｏｒ ｄｅ－ｂｌｕｒｒｉｎｇ ｓｐｉｒａｌ ＭＲ ｉｍａｇｅｓ」、Ｍａｇｎ．Ｒｅｓｏｎ．Ｍｅｄ．２０００、ｖｏｌ．４４、４９１～４９４ページ；Ｓｕｔｔｏｎらの「Ｆａｓｔ，ｉｔｅｒａｔｉｖｅ ｉｍａｇｅ ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ ｆｏｒ ＭＲＩ ｉｎ ｔｈｅ ｐｒｅｓｅｎｃｅ ｏｆ ｆｉｅｌｄ ｉｎｈｏｍｏｇｅｎｅｉｔｉｅｓ」、ＩＥＥＥ Ｔｒａｎｓ．Ｍｅｄ．Ｉｍａｇｉｎｇ．２００３、ｖｏｌ．２２、１７８～１８８ページ；Ｎａｙａｋらの「Ｅｆｆｉｃｉｅｎｔ ｏｆｆ－ｒｅｓｏｎａｎｃｅ ｃｏｒｒｅｃｔｉｏｎ ｆｏｒ ｓｐｉｒａｌ ｉｍａｇｉｎｇ」、Ｍａｇｎ．Ｒｅｓｏｎ．Ｍｅｄ．２００１、ｖｏｌ．４５、５２１～５２４ページを参照）。

しかしながら、前述のタイプのぼけ除去方法を適用した後でも、アーチファクトは、しばしば、非常に強い感受性に誘起された磁場勾配の画像領域に残存する。スパイラルｋ空間軌道の場合、このようなアーチファクトは、特徴的なリンギングアーチファクトとして、再構成及びぼけ除去されたＭＲ画像に現れ、関心対象の解剖学的詳細に重畳し得、又はこれを覆い得る。このような残存するアーチファクトの理由は、強い（典型的には感受性に誘起された）局所的な磁場勾配の場合、スパイラルｋ空間軌道の形状が、それに応じた大きな範囲にわたって、それぞれのボクセルについて理論的なスパイラル形状から外れるからである。これは、図４の２次元的ｋ空間図において示されている。図４ａは、完全に均一な主磁場Ｂ_０が存在する状態においてｋ_ｘ及びｋ_ｙ方向において正弦波的に変調された磁場勾配の適用によって得られる「理想的な」スパイラルｋ空間軌道を図示する。しかしながら、図４ｂにおいては、ｘ方向における強い勾配によってＢ_０は不均一であり、対応するボクセル位置は、理想的なスパイラル形状から著しく外れた歪んだｋ空間軌道になる。この効果の結果、ｋ空間サンプリングが不十分になり、ｋ空間のいくつかの部分においてナイキスト基準に違反する。図４ｂにおいては、ｋ空間の多くの部分から不十分な信号データがサンプリングされる。これらの効果によって、容易には補正できない強い局所的なオフ共鳴の場所において、局所的な勾配に誘起された不十分なｋ空間サンプリングによって生じる不所望の残存アーチファクト（以下において不良サンプリングアーチファクトと称される）がもたらされ得る。

前述のことから、向上されたＭＲ撮像技術に対する需要があることが容易に理解される。本発明の目的は、上述された限界に対処し、強いＢ_０不均一性／強いＢ_０勾配がある状況においても、効率的で高品質な非デカルトＭＲ撮像を可能とすることである。

本発明によると、ＭＲデバイスの検査ボリュームに位置付けられた物体のＭＲ撮像の方法が開示される。方法は、
少なくとも１つのＲＦ励起パルスと変調された磁場勾配とを含む撮像シーケンスを物体に施すステップと；
少なくとも１つの非デカルトｋ空間軌道に沿ってＭＲ信号を取得するステップと；
取得されたＭＲ信号からＭＲ画像を再構成するステップと；
Ｂ_０不均一性に誘起された不十分なｋ空間サンプリングによって生じる１つ又は複数の不良サンプリングアーチファクトを、深層学習ネットワークを使用して検知するステップと
を有する。

再構成された非デカルト（特には、スパイラル）ＭＲ画像において残存する不良サンプリングアーチファクトの原因を、深層学習ネットワークを使用して自動的に特定することが本発明の要旨である。不良サンプリングアーチファクトが検知されると、それらは、目標とされたやり方において、適切なアルゴリズムによって補正され得る。

好ましい実施形態において、再構成されたＭＲ画像は、不良サンプリングアーチファクトを検知するステップの前に、Ｂ_０マップに基づいてぼけ除去され、例えば、これは、ＭＲ画像を再構成するステップに統合される。従来のＢ_０マップに基づくぼけ除去を残余の不良サンプリングアーチファクトを検知する深層学習と組み合わせることによって、強い感受性効果及び鋭い局所的な磁場勾配がある場合でも、特に高品質なＭＲ画像が取得され得る。

更なる好ましい実施形態において、深層学習ネットワークは、ＭＲ画像からアーチファクトマップを導出するように訓練され、アーチファクトマップは、検知された不良サンプリングアーチファクトのみの図的表現である。これを達成するために、深層学習ネットワークは、好ましくは、その出力（最後のネットワーク層）においてモデル化されたアーチファクトマップのセットによって、及びその入力（最初のネットワーク層）においてそれぞれのモデル化されたアーチファクトマップと訓練ＭＲ画像との重畳によって訓練される。例えば、モデル化されたアーチファクトマップは、種々のオフ共鳴周波数において計算された単一又は複数のボクセルのオフ共鳴の点広がり関数を含み得る。換言すれば、本発明の手法は、不良サンプリングアーチファクトを、それらの最も単純な形態において、強い局所的なオフ共鳴が存在する状態における点広がり関数（使用された非デカルトｋ空間サンプリングスキームに対応する）として模倣し、これらを、深層学習ネットワークを訓練するために使用することである。より現実的な外観のために、不良サンプリングアーチファクトは、種々の点広がり関数の小さな、畳み込まれた局所的な集合体としてモデル化され、その各々は僅かに異なるオフ共鳴周波数におけるものであるが空間的に接近している。アーチファクトモデルのなおも更なる改善のために、オフ共鳴ボクセルのより複雑な局所的集合体も考慮するように小さな幾何学的歪曲が点広がり関数に付与され得る。こうして、深層学習ネットワークは、（ネットワークの出力において）予想される不良サンプリングアーチファクトのみを示す画像としてのモデル化されたアーチファクトマップによって訓練され、それに加えて、（ネットワークの入力において）それぞれの不良サンプリングアーチファクト及び（アーチファクトのない）訓練ＭＲ画像の対応する重畳によって訓練される。この手法は、「残余学習」の概念を採用する。これは、複雑で大きなＭＲ画像データに適用され得る。

深層学習ネットワークの出力としてアーチファクトマップが使用できるようになると、検知された不良サンプリングアーチファクトは、例えば、再構成された（及び、任意選択的には、ぼけ除去された）ＭＲ画像からアーチファクトマップを減算することによって、容易に補正され得る。

更なる好ましい実施形態において、補正は、予め定められた画像領域及び／又は不均一性又は主磁場の局所的な変化（勾配）の程度が所与の閾値を超えることをＢ_０マップが示す画像領域に制限される。画像品質を最適化するために、Ｂ_０マップにおける情報が、アーチファクト減算処理を、アーチファクトの潜在的な原因であることが明白である画像領域に適切に制限又は制約するために使用され得る。このようにして、不良サンプリングアーチファクトに対して構造的類似性を示す解剖学構造が除去されて、ＭＲ画像における重要な診断情報を潜在的に損なうことが自動的に回避される。可能な実施形態において、Ｂ_０マップ情報は、アーチファクト検知をガイドするために深層学習ネットワークを訓練する処理において直接的に使用される。ネットワークの入力として残余の不良サンプリングアーチファクトとともにＭＲ画像のみを使用する代わりに、Ｂ_０マップもネットワークに対してその入力として供給される。このようにして、深層学習ネットワークは、Ｂ_０マップ情報を自動的に考慮に入れ、それによって、有用な画像構造をアーチファクトとして間違って特定することを回避する。

更なる実施形態において、解剖学的地図が撮像データ（好ましくは、ボリュメトリックデータ、例えば、３Ｄ又はマルチスライスデータ）に適用され得、不良サンプリングアーチファクトの外観（例えば、頭部においては、眼の前頭皮質近くの領域又は内耳空洞など）に本質的な予め定められた画像領域の特定を可能とする。これは、例えばＢ_０マップが使用できない場合に、１つ又は複数の不良サンプリングアーチファクトの検知を制限するために使用され得る。このようなイベントにおいて、特定されたアーチファクトの補正も、これらの予め定められた領域に制限される。

なおも別の好ましい実施形態において、深層学習ネットワークは畳み込みネットワークである。いわゆるＦ－Ｎｅｔネットワークアーキテクチャ（Ｂｒｏｓｃｈ Ｔ、Ｓａａｌｂａｃｈ Ａの「Ｆｏｖｅａｌ Ｆｕｌｌｙ Ｃｏｎｖｏｌｕｔｉｏｎａｌ Ｎｅｔｓ ｆｏｒ Ｍｕｌｔｉ－Ｏｒｇａｎ Ｓｅｇｍｅｎｔａｔｉｏｎ」、Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ｔｈｅ ＳＰＩＥ、ｖｏｌｕｍｅ１０５７４、２０１８を参照）が使用され得、これは複数の解像度レベルに依存して高レベル及び低レベルの特徴を抽出する。代替的に、他のネットワークアーキテクチャが用いられてよい。例えば、いわゆるＵ－Ｎｅｔアーキテクチャ（Ｒｏｎｎｅｂｅｒｇｅｒ，Ｏ、Ｆｉｓｃｈｅｒ，Ｐ、Ｂｒｏｘ，Ｔの「Ｕ－ｎｅｔ：Ｃｏｎｖｏｌｕｔｉｏｎａｌ ｎｅｔｗｏｒｋｓ ｆｏｒ ｂｉｏｍｅｄｉｃａｌ ｉｍａｇｅ ｓｅｇｍｅｎｔａｔｉｏｎ」、Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｃｏｎｆｅｒｅｎｃｅ ｏｎ Ｍｅｄｉｃａｌ ｉｍａｇｅ ｃｏｍｐｕｔｉｎｇ ａｎｄ ｃｏｍｐｕｔｅｒ－ａｓｓｉｓｔｅｄ ｉｎｔｅｒｖｅｎｔｉｏｎ、２３４～２４１ページ、２０１５を参照）が使用されてもよい。もしもＢ_０マップ及び入力ＭＲ画像が、非常に異なる解像度を有するならば、調整されたＹ－Ｎｅｔアーキテクチャ（Ｍｏｈａｍｍｅｄ，Ａらの「Ｙ－Ｎｅｔ：Ａ ｄｅｅｐ Ｃｏｎｖｏｌｕｔｉｏｎａｌ Ｎｅｕｒａｌ Ｎｅｔｗｏｒｋ ｆｏｒ Ｐｏｌｙｐ Ｄｅｔｅｃｔｉｏｎ」、ａｒＸｉｖ ｐｒｅｐｒｉｎｔ ａｒＸｉｖ：１８０６．０１９０７、２０１８を参照）が適切である。このアーキテクチャにおいて、しばしばより小さな解像度を有するＢ_０マップの符号化パスが、より少ない層及び／又はチャンネルを使用して実現され得、それによって、ネットワークの全体的なサイズが低減される。

ここまでに説明された本発明の方法は、検査ボリューム内に均一な静的磁場を生成するための少なくとも１つの主電磁コイルと、検査ボリューム内で種々の空間的方向におけるスイッチされた磁場勾配を生成するためのいくつかの勾配コイルと、検査ボリューム内にＲＦパルスを生成するため及び／又は検査ボリュームに位置付けられた物体からのＭＲ信号を受信するための少なくとも１つのＲＦコイルと、時間的に連続するＲＦパルス及びスイッチされた磁場勾配を制御するための制御ユニットと、受信されたＭＲ信号からＭＲ画像を再構成するための再構成ユニットとを含むＭＲデバイスによって実行され得る。本発明の方法は、例えば、ＭＲデバイスの再構成ユニット及び／又は制御ユニットの対応するプログラミングによって実現され得る。

本発明の方法は、現在の臨床的使用においてほとんどのＭＲデバイスにおいて有利に実行され得る。このためには、ただ単に、ＭＲデバイスが上に説明された発明の方法のステップを実施するようにＭＲデバイスを制御するコンピュータプログラムを利用することが必要なだけである。コンピュータプログラムは、データ記憶媒体に存在してよく、又はＭＲデバイスの制御ユニットへのインストールのためにダウンロードされるようにデータネットワークに存在してよい。

更には、本発明の不良サンプリングアーチファクトの検知及び補正の手法は、診断ワークステーションにインストールされる遡及的なアーチファクト除去ソフトウェアツールとして実現され得る。

添付の図面は、本発明の好ましい実施形態を開示する。しかしながら、図面は例示のみを目的とするものと企図されており、本発明の限界を定めるものではないことが理解されるべきである。

本発明の方法を実行するためのＭＲデバイスを図示する。 Ｂ_０マップ並びに再構成及びぼけ除去されたＭＲ頭部画像を種々の局所的なオフ共鳴のためにモデル化された点広がり関数のセットとともに図示する。 スパイラルＭＲ頭部画像における不良サンプリングアーチファクト補正のための本発明の手法を示す。 スパイラルＭＲ撮像の例における不良サンプリングアーチファクトの原因を示すｋ空間の図を図示する。

図１を参照すると、ＭＲデバイス１が概略的に図示されている。デバイスは、略均一で時間的に一定な主磁場が検査ボリュームを通ってｚ軸に沿って生まれるように、超伝導又は抵抗性主電磁コイル２を備える。

磁気共鳴生成及び操作システムは、ＭＲ撮像を実施するために、一連のＲＦパルス及びスイッチされた磁場勾配を印加して、核磁気スピンを反転若しくは励起すること、磁気共鳴を誘起すること、磁気共鳴をリフォーカスさせること、磁気共鳴を操作すること、空間的に若しくは他のやり方で磁気共鳴を符号化すること、スピンを飽和させることなどを行う。

より具体的には、勾配パルス増幅器３は、検査ボリュームのｘ、ｙ及びｚ軸に沿った全身勾配コイル４、５、及び６のうちの選択された１つに電流パルスを印加する。デジタルＲＦ周波数送信器７は、送信／受信スイッチ８を介して、ＲＦパルス又はパルスパケットを全身ボリュームＲＦコイル９に送信して、検査ボリューム内にＲＦパルスを送信する。典型的なＭＲ撮像シーケンスは、持続時間の短いＲＦパルスセグメントのパケットから構成され、これらは、互いとともに及び任意の印加された磁場勾配とともに、核磁気共鳴の選択された操作を達成する。ＲＦパルスは、共鳴を飽和、励起し、磁化を反転し、共鳴をリフォーカスさせ、又は共鳴を操作し、検査ボリューム内に位置付けられた身体１０の一部分を選択するために使用される。ＭＲ信号は全身ボリュームＲＦコイル９によってもピックアップされる。

身体１０の限定的な領域のＭＲ画像の生成のために、局所的アレイＲＦコイル１１、１２、１３のセットが、撮像のために選択された領域に隣接して配置される。アレイコイル１１、１２、１３は、身体コイルＲＦ送信によって誘起されたＭＲ信号を受信するために使用され得る。

結果的なＭＲ信号は、全身ボリュームＲＦコイル９及び／又はアレイＲＦコイル１１、１２、１３によってピックアップされ、好ましくはプレ増幅器（図示せず）を含む受信器１４によって復調される。受信器１４は、送信／受信スイッチ８を介してＲＦコイル９、１１、１２及び１３に接続される。

ホストコンピュータ１５は、本発明によるスパイラルｋ空間軌道に沿ったＭＲ信号を取得するために、勾配パルス増幅器３及び送信器７を制御して、エコー平面撮像（ＥＰＩ）、エコーボリューム撮像、勾配及びスピンエコー撮像、高速スピンエコー（ＴＳＥ）撮像などの複数のＭＲ撮像シーケンスのうちの任意のものを生成する。選択されたシーケンスについて、受信器１４は、各ＲＦ励起パルスに続いて、単一の又は複数のＭＲデータをそれぞれのｋ空間軌道に沿って迅速に連続して受信する。データ取得システム１６は、受信された信号のアナログ－デジタル変換を実施し、更なる処理に適したデジタルフォーマットに各ＭＲ信号を変換する。最新のＭＲデバイスにおいては、データ取得システム１６は、未加工画像データの取得に特化された別個のコンピュータである。

最終的に、デジタル未加工画像データは、フーリエ変換又は他の適切な再構成アルゴリズムを適用する再構成プロセッサ１７によって画像表現に再構成される。ＭＲ画像は、患者を通る平面スライス、平行平面スライスのアレイ、３次元的ボリュームなどを表現する。次いで、画像は、画像メモリに記憶され、画像メモリにおいて、画像は、画像表現のスライス、投影又は他の部分を、例えば結果的なＭＲ画像の人間が読み取り可能な表示を提供するビデオモニタ１８を介した視覚化のための適切なフォーマットに変換するためにアクセスされる。

ＭＲデバイス１は、例えばホストコンピュータ１５及び再構成プロセッサ１７の適切なプログラミングによって、本明細書において上記及び下記に説明される本発明の撮像方法を実施するように構成される。

引き続き図１を参照するとともに更に図２及び図３を参照すると、本発明の撮像手法の実施形態が説明されている。

例示的な実施形態において、スピンエコーＴ_２重み付きスパイラル撮像が実施される。６０ｍｓのエコー時間における５０ｍｓの取得窓中に、ただ１つのスパイラル軌道が取得される。減衰因子Ｒ＝３の平行撮像が使用される。反復的なスパイラルＳＥＮＳＥ画像再構成及び適切なＢ_０マップに基づくぼけ除去を使用して、スパイラルスピンエコー（ＳＥ）ＭＲ画像（図２における上段右の画像を参照）が得られる。また、ぼけ除去のために使用される（別個に取得された）Ｂ_０磁場マップが図２に図示されている（上段左の画像）。スキャン効率と、Ｂ_０マップにおいて現れ得る高い、変化の激しい局所的な場の不均一性とを最適化するために選ばれた長い取得窓に起因して、Ｂ_０マップに基づいたオフ共鳴補正（ぼけ除去）の後にも、再構成されたＭＲ画像には未補正の不良サンプリングアーチファクトが残存している。残存している不良サンプリングアーチファクトは、矢印によって示されるリンギングタイプの構造として現れる。図２の下段の行において、計算された点広がり関数（単一ボクセル信号表現）が、Ｈｚの単位で与えられた種々のオフ共鳴周波数について、及びＭＲ信号取得において使用される所与のスパイラルｋ空間軌道について図示される。上段右のＭＲ画像における内耳に近い場所について、Ｂ_０マップは、約２５０Ｈｚの逸脱を予測している。画像において、点広がり関数の対応する内側リング寸法が重ね合されている。この点広がり関数は、マークされた場所（小さい円）の中央のボクセル又はボクセルの集合体から生じる信号の寄与を模倣し、これは、強い局所的な勾配の存在によって認識できないほどに歪んでおり、Ｂ_０マップに基づいてこれ以上の補正を行うことはできない。不良サンプリングアーチファクトは、１つのボクセルからもたらされるだけではない。異なる強度を有するいくつかの互いに接近したボクセルの小集合体から不良サンプリングアーチファクトがもたらされ得、これは実際の局所的な不均一性の程度が異なり、より複雑な、歪んだ、非回転的な又は捩れたアーチファクトパターンをもたらす。上述されたように、計算された点広がり関数が、深層学習ネットワークを訓練するためのモデルアーチファクトマップとして使用される。図示される実施形態において、複数の解像度レベルに依存して高レベル及び低レベルの両方の特徴を抽出するＦ－Ｎｅｔアーキテクチャが使用される。それぞれが２つの畳み込み層を有する３つの異なる解像度レベルが用いられる。

図３において示されるように、訓練された深層学習ネットワークは、不良サンプリングアーチファクトを含むスパイラルＭＲ画像（上段左の画像）を分析する。ネットワークは、ＭＲ画像からアーチファクトマップを導出し、このアーチファクトマップは、ＭＲ画像における不良サンプリングアーチファクトの推定である。検知されたアーチファクトは、ＭＲ画像からアーチファクトマップを減算することで、アーチファクトマップに基づいて補正される。補正されたＭＲ画像は、図３における下段の画像である。アーチファクトマップにおける矢印は、深層学習ネットワークの解釈の誤りに基づいて間違って検知された不良サンプリングアーチファクトを示す。Ｂ_０マップを使用して、強い不均一性及び鋭い勾配を有する領域が特定され得る。残余の不良サンプリングアーチファクトは、ＭＲ画像のこれらの領域にのみ生じるものと予想され得る。従って、導出されたアーチファクトマップに基づくＭＲ画像の補正は、これらの領域にのみ許可される。Ｂ_０マップからの情報は、対応する減算マスクを導出するために使用され得、この減算マスクは、強い不均一性及び鋭い勾配がある場所において「１」であり、疑わしさがない場所、すなわち局所的な不均一性変化が臨界閾値よりも低く、オフ共鳴勾配に誘起された不十分なサンプリングに起因するアーチファクトが予想され得ない場所において、滑らかに「０」まで下がる。

更なる実施形態において、解剖学的地図が撮像データ（好ましくは、ボリュメトリックデータ、例えば、３Ｄ又はマルチスライスデータ）に適用され得、このようなアーチファクトの外観（頭部においては、眼の前頭皮質近くの領域又は内耳空洞など）に本質的な、強い局所的な感受性勾配のある領域の特定を可能とし、磁場マップ情報が使用できない場合に、この重み付きアーチファクト減算を支援する。

このようにして、アーチファクトに対して類似性を示す解剖学的特徴が誤って除去されて、画像の臨床的価値に影響を与えることが自動的に回避される。

Claims (12)

1. ＭＲデバイスの検査ボリューム内に位置付けられた物体のＭＲ撮像の方法であって、前記方法は、
   少なくとも１つのＲＦ励起パルスと変調された磁場勾配とを含む撮像シーケンスを前記物体に施すステップと、
   少なくとも１つの非デカルトｋ空間軌道に沿ってＭＲ信号を取得するステップと、
   取得された前記ＭＲ信号からＭＲ画像を再構成するステップと、
   前記ＭＲ画像における不均一性に誘起された不十分なｋ空間サンプリングによって生じる１つ又は複数の不良サンプリングアーチファクトを、深層学習ネットワークを使用して検知するステップと
   を有する、方法。
2. 前記非デカルトｋ空間軌道は、スパイラルｋ空間軌道である、請求項１に記載の方法。
3. 再構成された前記ＭＲ画像は、残存する不良サンプリングアーチファクトを前記検知するステップの前に、Ｂ_０マップに基づいてぼけ除去される、請求項１又は２に記載の方法。
4. 前記深層学習ネットワークは、前記ＭＲ画像からアーチファクトマップを導出するように訓練され、前記アーチファクトマップは、少なくとも１つの検知された前記不良サンプリングアーチファクトのみの図的表現である、請求項１から３のいずれか一項に記載の方法。
5. 前記深層学習ネットワークは、その出力においてモデル化されたアーチファクトマップのセットによって、及びその入力においてそれぞれのモデル化された前記アーチファクトマップと訓練ＭＲ画像との重畳によって訓練される、請求項４に記載の方法。
6. モデル化された前記アーチファクトマップは、使用された前記撮像シーケンスに関して計算された単一又は複数のボクセルのオフ共鳴の点広がり関数を含む、請求項５に記載の方法。
7. 検知された前記不良サンプリングアーチファクトは、前記深層学習ネットワークによって、再構成された前記ＭＲ画像から導出された前記アーチファクトマップに基づいて補正される、請求項４から６のいずれか一項に記載の方法。
8. 前記１つ又は複数の不良サンプリングアーチファクトを検知するステップは、予め定められた画像領域及び／又は前記不均一性又は主磁場の局所的な変化の程度が所与の閾値を超えることをＢ_０マップが示す画像領域に制限される、請求項１から７のいずれか一項に記載の方法。
9. Ｂ_０マップが、前記不良サンプリングアーチファクトの検知中に前記深層学習ネットワークの更なる入力として使用される、請求項１から８のいずれか一項に記載の方法。
10. 前記深層学習ネットワークは畳み込みネットワークである、請求項１から９のいずれか一項に記載の方法。
11. 検査ボリューム内に均一な静的磁場を生成するための少なくとも１つの主電磁コイルと、前記検査ボリューム内で種々の空間的方向におけるスイッチされた磁場勾配を生成するためのいくつかの勾配コイルと、前記検査ボリューム内にＲＦパルスを生成するため及び／又は前記検査ボリュームに位置付けられた物体からのＭＲ信号を受信するための少なくとも１つのＲＦコイルと、時間的に連続するＲＦパルス及びスイッチされた磁場勾配を制御するための制御ユニットと、受信された前記ＭＲ信号からＭＲ画像を再構成するための再構成ユニットとを含むＭＲデバイスであって、前記ＭＲデバイスは、
    少なくとも１つのＲＦ励起パルスと変調された磁場勾配とを含む撮像シーケンスを前記物体に施すことと、
    少なくとも１つの非デカルトｋ空間軌道に沿ってＭＲ信号を取得することと、
    取得された前記ＭＲ信号からＭＲ画像を再構成することと、
    前記ＭＲ画像における不均一性に誘起された不十分なｋ空間サンプリングによって生じる１つ又は複数の不良サンプリングアーチファクトを、深層学習ネットワークを使用して検知することと
    を実施する、ＭＲデバイス。
12. 非デカルトｋ空間サンプリングを使用して取得されたＭＲ信号からＭＲ画像を再構成することと、
    前記ＭＲ画像における不均一性に誘起された不十分なｋ空間サンプリングによって生じる１つ又は複数の不良サンプリングアーチファクトを、深層学習ネットワークを使用して検知することと
    のための命令を含む、コンピュータプログラム。

Images