特許 7187715 加速された磁気共鳴温度測定

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2022-12-02

(45)【発行日】2022-12-12

(54)【発明の名称】加速された磁気共鳴温度測定

(51)【国際特許分類】

   A61B 5/055 20060101AFI20221205BHJP

【ＦＩ】

A61B5/055 376

【請求項の数】 23

(21)【出願番号】P 2021561953

(86)(22)【出願日】2020-04-21

(65)【公表番号】

(43)【公表日】2022-06-23

(86)【国際出願番号】 IB2020000294

(87)【国際公開番号】W WO2020217098

(87)【国際公開日】2020-10-29

【審査請求日】2021-10-18

(31)【優先権主張番号】62/931,525

(32)【優先日】2019-11-06

(33)【優先権主張国・地域又は機関】US

(31)【優先権主張番号】62/838,409

(32)【優先日】2019-04-25

(33)【優先権主張国・地域又は機関】US

(31)【優先権主張番号】62/971,449

(32)【優先日】2020-02-07

(33)【優先権主張国・地域又は機関】US

(73)【特許権者】

【識別番号】508154863

【氏名又は名称】インサイテック・リミテッド

(74)【代理人】

【識別番号】100078282

【弁理士】

【氏名又は名称】山本 秀策

(74)【代理人】

【識別番号】100113413

【弁理士】

【氏名又は名称】森下 夏樹

(74)【代理人】

【識別番号】100181674

【弁理士】

【氏名又は名称】飯田 貴敏

(74)【代理人】

【識別番号】100181641

【弁理士】

【氏名又は名称】石川 大輔

(74)【代理人】

【識別番号】230113332

【弁護士】

【氏名又は名称】山本 健策

(72)【発明者】

【氏名】ズール， ユヴァル

(72)【発明者】

【氏名】シャピラ， ボアズ

(72)【発明者】

【氏名】レヴィ， ヨアフ

【審査官】下村 一石

(56)【参考文献】

【文献】特表２０１４－５０６８２２（ＪＰ，Ａ）

【文献】米国特許出願公開第２０１２／０２５７８０６（ＵＳ，Ａ１）

【文献】国際公開第２０１８／１２６３９６（ＷＯ，Ａ１）

【文献】YUAN, Jing et al.，Towards fast and accurate temperature mapping with proton resonance frequency-based MR thermometry，Quantitative Imaging in Medicine and Surgery，AME Publishing company，2012年01月26日，2012;2(1)，pp.21-32

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ａ６１Ｂ ５／０５５

(57)【特許請求の範囲】

【請求項1】

撮像システムであって、
磁気共鳴映像法（ＭＲＩ）装置であって、前記磁気共鳴映像法（ＭＲＩ）装置は、
（ｉ）第１のマルチエコーパルスシーケンスを実行することにより、標的領域を励起することと、
（ｉｉ）前記標的領域の複数のベースラインｋ空間画像を入手することであって、各ベースライン画像は、前記第１のマルチエコーパルスシーケンス内のエコーに関連付けられている、ことと、
（ｉｉｉ）第２のマルチエコーパルスシーケンスを実行することにより、前記標的領域を励起することと、
（ｉｖ）第２の複数のｋ空間画像を入手することであって、前記第２の複数のｋ空間画像のそれぞれは、前記第２のマルチエコーパルスシーケンス内のエコーに関連付けられており、前記第２の複数のｋ空間画像のうちの少なくとも１つは、前記ベースラインｋ空間画像と比較して情報が欠落するように、アンダーサンプリングされる、ことと
を行うように構成されている、磁気共鳴映像法（ＭＲＩ）装置と、
少なくとも１つのベースラインｋ空間画像および前記第２の複数のｋ空間画像のうちの少なくとも１つに少なくとも部分的に基づいて、標的画像を再構築するように構成されている算出ユニットと
を備え、
前記算出ユニットは、前記第２の複数のｋ空間画像のうちの少なくとも１つに関して、（ｉ）それからの欠落情報を算出的に推定することと、（ｉｉ）それに関連付けられているエコー時間に少なくとも部分的に基づいて、前記推定される欠落情報を算出的に更新することと、（ｉｉｉ）前記更新された推定される欠落情報と、それに対応する前記第２の複数のｋ空間画像のうちの前記少なくとも１つと、前記第２の複数のｋ空間画像のうちの前記少なくとも１つに対応する前記ベースラインｋ空間画像とに少なくとも部分的に基づいて、マップを生成することとを行うように構成されている、撮像システム。

【請求項2】

前記標的画像は、熱マップ、磁化率強調画像、Ｔ_２ ^＊画像、微小出血画像、または、ＢＯＬＤｆ－ＭＲＩ画像のうちの１つである、請求項１に記載のシステム。

【請求項3】

前記算出ユニットは、ニューラルネットワークを備える、請求項１に記載のシステム。

【請求項4】

前記ニューラルネットワークは、畳み込みニューラルネットワークである、請求項３に記載のシステム。

【請求項5】

前記マップは、前記標的領域の解剖学的画像または熱マップのうちの少なくとも一方を含む、請求項１に記載のシステム。

【請求項6】

前記算出ユニットは、デフォルト値、前記対応するベースラインｋ空間画像のうちの１つの一部、または、前記第２の複数のｋ空間画像のうちの異なるものに少なくとも部分的に基づいて、前記欠落情報を算出的に推定するようにさらに構成されている、請求項１に記載のシステム。

【請求項7】

前記複数のベースラインｋ空間画像のうちのそれぞれ１つは、複数のピクセルを備え、前記算出ユニットは、前記複数のベースラインｋ空間画像のうちの１つの中のピクセルに関連付けられている大きさに少なくとも部分的に基づいて、前記推定される欠落情報を更新するようにさらに構成されている、請求項１に記載のシステム。

【請求項8】

前記システムは、第１のＭＲＩコイルと第２のＭＲＩコイルとをさらに備え、前記複数のベースラインｋ空間画像のそれぞれおよび前記第２の複数のｋ空間画像のそれぞれは、前記第１のＭＲＩコイルまたは前記第２のＭＲＩコイルのうちの少なくとも一方に関連付けられている、請求項１に記載のシステム。

【請求項9】

前記算出ユニットは、（ａ）前記複数のベースラインｋ空間画像から複数のベースライン画像を算出的に再構築することと、（ｂ）前記第２の複数のｋ空間画像のうちのそれぞれ１つに関して、それに関連付けられている前記算出的に推定される欠落情報に少なくとも部分的に基づいて、第３の画像セット内の画像を算出的に再構築することとを行うようにさらに構成されている、請求項１に記載のシステム。

【請求項10】

前記複数のベースラインｋ空間画像のうちのそれぞれ１つは、前記第１のマルチエコーパルスシーケンスに応答してエコーから入手され、前記第２の複数のｋ空間画像のうちのそれぞれ１つは、前記第２のマルチエコーパルスシーケンスに応答してエコーから入手され、前記算出ユニットは、それに関連付けられている前記エコーのエコー時間に少なくとも部分的に基づいて、前記第３の画像セット内の各再構築された画像を前記再構築されたベースライン画像のうちの１つに対応させるようにさらに構成されている、請求項９に記載のシステム。

【請求項11】

前記システムは、複数のＭＲＩコイルをさらに備え、前記複数のベースラインｋ空間画像のそれぞれおよび前記第２の複数のｋ空間画像のそれぞれは、前記ＭＲＩコイルのうちの少なくとも１つに関連付けられており、前記算出ユニットは、前記関連付けられている少なくとも１つのＭＲＩコイルに少なくとも部分的に基づいて、前記第３の画像セット内の各再構築された画像を前記再構築されたベースライン画像のうちの１つに関連させるようにさらに構成されている、請求項１０に記載のシステム。

【請求項12】

前記第３の画像セット内の前記再構築された画像に関連付けられている前記エコー時間は、前記対応する再構築されたベースライン画像に関連付けられている前記エコー時間と同一である、請求項１０に記載のシステム。

【請求項13】

前記算出ユニットは、（ｃ）前記再構築されたベースライン画像のうちの１つの中のピクセルと前記第３の画像セット内の前記対応する再構築された画像の中の対応するピクセルとの間の位相差を決定するようにさらに構成されている、請求項１０に記載のシステム。

【請求項14】

前記算出ユニットは、（ｄ）それから前記第３の画像セット内の前記対応する再構築された画像が再構築される前記ｋ空間画像のエコー時間に少なくとも部分的に基づいて、前記位相差を更新するようにさらに構成されている、請求項１３に記載のシステム。

【請求項15】

前記第３の画像セット内の前記再構築された画像に関連付けられている前記位相差は、前記第２の複数のｋ空間画像のエコー時間と正に相関する、請求項１４に記載のシステム。

【請求項16】

前記算出ユニットは、前記第２の複数のｋ空間画像内の少なくとも２つの異なる画像に関連付けられている前記エコー時間に少なくとも部分的に基づいて、前記位相差を更新するようにさらに構成されている、請求項１４に記載のシステム。

【請求項17】

前記算出ユニットは、（ｅ）前記更新された位相差に少なくとも部分的に基づいて、前記第３の画像セット内の前記対応する再構築された画像を更新することと、（ｆ）前記更新された画像を第４の画像セット内のｋ空間画像に変換することとを行うようにさらに構成されている、請求項１４に記載のシステム。

【請求項18】

前記算出ユニットは、（ｇ）同一のエコー時間において入手される前記第２の複数のｋ空間画像内の前記対応する画像に少なくとも部分的に基づいて、前記第４の画像セット内の前記ｋ空間画像を更新するようにさらに構成されている、請求項１７に記載のシステム。

【請求項19】

前記算出ユニットは、終了条件が満たされるまで、ステップ（ｂ）～（ｇ）を反復的に実行するようにさらに構成されている、請求項１８に記載のシステム。

【請求項20】

前記終了条件は、
所定の限界値を超過する反復の回数、または、
２回の反復の間の前記第４の画像セット内の前記更新されたｋ空間画像または前記第３の画像セット内の前記再構築された画像の所定の最小値を下回る変化
のうちの１つ以上のものに対応する、請求項１９に記載のシステム。

【請求項21】

前記算出ユニットは、ステップ（ｂ）～（ｇ）を少なくとも２回実行するために、人工ニューラルネットワークを実行するようにさらに構成されている、請求項１８に記載のシステム。

【請求項22】

前記算出ユニットは、前記第４の画像セット内の前記更新されたｋ空間画像に少なくとも部分的に基づいて、第５の画像セット内の画像を算出的に再構築することと、前記第５の画像セット内の前記再構築された画像からノイズを算出的に低減させることとを行うようにさらに構成されている、請求項１８に記載のシステム。

【請求項23】

前記算出ユニットは、前記第３の画像セットまたは前記第５の画像セットのうちの少なくとも一方の中の前記画像のうちの少なくとも２つに、その中の前記ノイズを低減させるために、局所的に低ランクの正則化を適用するようにさらに構成されている、請求項２２に記載のシステム。

【発明の詳細な説明】

【技術分野】

【0001】

本願は、２０１９年４月２５日に出願された、米国仮出願第６２／８３８，４０９号、２０１９年１１月６日に出願された、第６２／９３１，５２５号、および２０２０年２月７日に出願された、第６２／９７１，４４９号の優先権および利益を主張する。前述の文書の全ての開示全体は、参照することによって本明細書に組み込まれる。

【0002】

本発明は、概して、磁気共鳴（ＭＲ）温度測定、特に、ＭＲ温度測定を加速させるためのアプローチに関する。

【背景技術】

【0003】

患者の頭蓋骨または他の身体領域内の良性または悪性腫瘍または血餅等の組織が、組織を外科手術的に除去することによって侵襲的に、または、例えば、熱アブレーションを使用することによって非侵襲的に治療され得る。両方のアプローチとも、身体内のある局所的条件を効果的に治療するが、そうでなければ健康な組織を破壊または損傷させることを回避するための繊細な手技を伴い得る。健康な組織が救われ得ない、またはその破壊が生理学的機能に悪影響を及ぼす可能性がある限り、外科手術は、罹患組織が健康な組織内に統合される条件に関して適切ではあり得ない。

【0004】

集束された超音波を使用して遂行され得るような熱アブレーションは、超音波エネルギーの影響が、明確に画定される標的領域に限局され得るため、健康な組織または器官によって囲繞される、またはそれに隣接する罹患組織を治療するための特定の魅力を有する。超音波エネルギーが、超音波の比較的に短い波長（例えば、１メガヘルツ（１ＭＨｚ）における断面が１．５ミリメートル（ｍｍ）と同程度に小さい）に起因して数ミリメートルのみの断面を有する区域に集束され得る。また、音響エネルギーは、概して、軟組織を通して良好に透過するため、介在する解剖学的構造は、多くの場合、所望の焦点区域を画定することに対する障害物を課さない。したがって、超音波エネルギーは、囲繞健康組織を有意に損傷させることなく罹患組織をアブレートするために、小さい標的に集束され得る。

【0005】

超音波集束システムは、概して、音響トランスデューサ表面またはトランスデューサ表面のアレイを利用し、超音波ビームを発生させる。トランスデューサは、幾何学的に成形され、患者内の標的組織腫瘤に対応する「焦点区域」に超音波エネルギーを集束させるように位置付けられ得る。組織を通した波伝搬の間、超音波エネルギーの一部が、吸収され、上昇された温度、最終的には、好ましくは、焦点区域の中の標的組織腫瘤における、細胞壊死につながる。トランスデューサアレイの個々の表面または「要素」は、典型的には、個々に制御可能であり、すなわち、それらの位相および／または振幅が、（例えば、連続波の場合には好適な遅延または位相偏移を伴う「ビーム形成器」および要素のための増幅器回路網を使用して）相互から独立して設定され、ビームが、所望の方向に操向され、所望の距離に集束されることを可能にし、焦点区域性質が、必要に応じて成形されることを可能にすることができる。したがって、焦点区域は、トランスデューサ要素の中に入力される電気信号の振幅および位相を独立して調節することによって、迅速に変位および／または再成形されることができる。

【0006】

しかしながら、人体は、可撓性であり、患者が静止した状態を保つように位置付けられているときでも（例えば、呼吸、またはわずかな不随意移動に起因して）移動するため、複数の超音波処理として経時的に送達される治療は、超音波処理が相互から数秒以内に送達されるときでも、標的化への、および／または１つ以上の治療パラメータへの暫定的調節を要求し得る。運動のための補償が、したがって、超音波ビームが、標的上に集束されたままであり、囲繞健康組織を損傷させないことを確実にするために必要である。

【0007】

故に、磁気共鳴映像法（ＭＲＩ）等の撮像モダリティが、非侵襲性療法の間に集束する超音波と連動して使用され、標的組織および超音波焦点の両方の場所を監視し得る。概して、図１に描写されるように、ＭＲＩシステム１００は、静磁場磁石１０２と、１つ以上の勾配磁場コイル１０４と、高周波（ＲＦ）伝送機１０６と、ＲＦ受信機（図示せず）とを含む。（いくつかの実施形態では、同一のデバイスが、ＲＦ伝送機または受信機として交互に使用される。）磁石は、その中に患者１１０を受容するための領域１０８を含み、患者にわたって静的かつ比較的に均質な磁場を提供する。勾配磁場コイル１０４によって発生される時変磁場勾配が、静磁場と重畳される。ＲＦ伝送機１０６は、患者１１０にわたってＲＦパルスシーケンスを伝送し、患者の組織に（時変）ＲＦ応答信号を放出させ、これは、全体的な（２または３次元）撮像領域にわたって統合され、ＲＦ受信機によってサンプリングされ、未加工画像データを構成する、応答信号の時系列を生産する。本未加工データは、算出ユニット１１２に伝達される。時系列内の各データ点は、ｋ空間（すなわち、波ベクトル空間）内の特定の点における位置依存の局所的磁化のフーリエ変換の値として解釈されることができる（波ベクトルｋは、勾配磁場の時間展開の関数である）。したがって、応答信号の時系列を逆フーリエ変換することによって、算出ユニット１１２は、未加工データから組織の画像（すなわち、空間座標の関数として、測定された磁化に影響を及ぼす組織の性質を示す画像）を再構築することができる。磁気共鳴（ＭＲ）画像が、次いで、ユーザに表示され得る。ＭＲＩシステム１００は、医療手技を計画し、かつ手技の間に治療の進行を監視するために使用され得る。例えば、ＭＲＩは、解剖学的領域を撮像する、領域内の標的組織（例えば、腫瘍）を位置特定する、超音波トランスデューサ１１４によって発生されるビームを標的組織に誘導する、および／または標的組織の中およびそれを囲繞する温度を監視するために使用され得る。

【0008】

ＭＲＩは、種々の異なる治療シナリオにおいて、ＭＲＩ誘導下の集束超音波（ＭＲｇＦＵＳ）システム等の画像誘導下システムの有利な使用を可能にするが、多くの場合、撮像率（すなわち、連続ＭＲＩ画像が入手され得る率）は、標的の１つ以上の特性が変化する率に遅れを取る。例えば、集束超音波治療または他の療法エネルギーへの暴露の間に、標的の場所または温度は、急速に変化し得、１つの全体ＭＲＩ走査から別のものへの大規模な不連続ジャンプをもたらす。十分に細かいタイムスケールにおいて標的における変化を追跡することの本実行不可能性は、（治療を停止し、再設定することが、変化を考慮するために必要であり得るため）治療プロセスにおける有意な非効率、または標的の治療エネルギーへの、または標的化されていない組織への危険な暴露をももたらし得る。

【0009】

本問題を無効にするための１つの従来のアプローチが、走査するステップ（すなわち、ｋ空間データを「アンダーサンプリングする」ステップ）の間に部分的な未加工のｋ空間データのみを入手することによって、ＭＲ撮像率を加速させ、コスト関数を伴うアルゴリズムが、次いで、反復的に実施され、アンダーサンプリングされたｋ空間データを使用して画像を発生させ得る。コスト関数は、概して、最初に入手される疎データへの忠実性を確実にする、データ忠実度項と、経時的なデータの平滑性を最大化するモデルとして実現され得る、制約項とを含む。ｋ空間データ点の入手を低減させることは、ＭＲ撮像率を上昇させ得るが、撮像再構築プロセスの間の「欠落している」（または入手されていない）ｋ空間データについての情報を補完するためのプロセスは、再構築されるべき画像についての以前の知識の欠如に起因して時間がかかり得、結果として、上昇されたＭＲ撮像率が、損なわれる。概して、取得された以前の知識が多いほど、より少ないサンプリングされたｋ空間データが、要求され、これは、次いで、ＭＲ撮像率のより大きい加速を達成することができる。故に、再構築されるべき画像についての以前の付加的な情報を実装し、それによって、画像再構成のためにＭＲ走査の間に要求されるｋ空間データの量をさらに低減させ、かつＭＲ撮像率をさらに加速させる、改良されたＭＲＩ入手アプローチの必要性が、存在する。

【発明の概要】

【課題を解決するための手段】

【0010】

本発明の実施形態は、不完全なｋ空間データから再構築されるべき画像についての以前の知識を利用し、それによって、正確な再構築を促進する、加速されたＭＲ温度測定を提供する。種々の実施形態では、欠落データは、ニューラルネットワーク等の機械学習アルゴリズムを使用して算出的に推定され、画像が、反復的に更新された推定される欠落情報に基づいて発生される。

【0011】

故に、一側面では、本発明は、撮像システムであって、種々の実施形態では、標的領域を励起するように構成される、ＭＲＩ装置を備える、撮像システムに関する。標的領域は、特徴（例えば、ヒトまたは動物の解剖学的特徴、治療標的、非治療標的、撮像ファントム等のファントム等）を含む、それから本質的に成る、またはそれから成ってもよい。種々の実施形態では、ＭＲＩ装置は、（ｉ）第１のマルチエコーパルスシーケンスを実行し、標的領域を励起し、（ｉｉ）各ベースライン画像が、第１のマルチエコーパルスシーケンス内のエコーと関連付けられる、標的領域の複数のベースラインｋ空間画像を入手し、（ｉｉｉ）第２のマルチエコーパルスシーケンスを実行し、標的領域を励起し、（ｉｖ）第２の複数のｋ空間画像を入手するように構成される。第２の複数のｋ空間画像はそれぞれ、第２のマルチエコーパルスシーケンス内のエコーと関連付けられてもよく、第２の複数のｋ空間画像のうちの少なくとも１つは、ベースラインｋ空間画像と比較して情報が欠落するように、アンダーサンプリングされてもよい。いくつかの実施形態では、第２の複数のｋ空間画像の全てが、アンダーサンプリングされる。一実施形態では、第２の複数のｋ空間画像のうちの少なくとも１つは、異なるパターンにおいてアンダーサンプリングされる。加えて、画像システムは、少なくとも部分的に、少なくとも１つのベースラインｋ空間画像および第２の複数のｋ空間画像のうちの少なくとも１つに基づいて、標的画像を再構築するように構成される、算出ユニットを含んでもよい。

【0012】

いくつかの実施形態では、算出ユニットは、ニューラルネットワーク、例えば、畳み込みニューラルネットワーク等の機械学習アルゴリズムを備える、または実行する。標的画像は、熱マップ、磁化率強調画像、Ｔ_２ ^＊画像、微小出血画像、またはＢＯＬＤ ｆ－ＭＲＩ画像であってもよい。

【0013】

種々の実施形態では、本システムは、第２の複数のｋ空間画像のうちの少なくとも１つに関して、欠落情報が、ニューラルネットワークを使用して算出的に推定され、標的画像が、少なくとも部分的に、更新された推定される欠落情報に基づいて発生されるように、構成される。算出ユニットは、第２の複数のｋ空間画像のうちの少なくとも１つに関して、（ｉ）ニューラルネットワークを使用して、対応する解剖学的画像を算出的に推定し、（ｉｉ）少なくとも部分的に、更新された推定される解剖学的画像に基づいて、標的画像を発生させるように構成されてもよい。代替として、または加えて、算出ユニットは、第２の複数のｋ空間画像のうちの少なくとも１つに関して、ｉ）それからの欠落情報を推定し、対応する解剖学的画像を再構築し、（ｉｉ）ニューラルネットワークを使用して、解剖学的画像を算出的に補正し、（ｉｉｉ）少なくとも部分的に、更新された推定される解剖学的画像に基づいて、標的画像を発生させるように構成されてもよい。さらに別の代替物では、算出ユニットは、第２の複数のｋ空間画像のうちの少なくとも１つに関して、（ｉ）それからの欠落情報を算出的に推定し、（ｉｉ）少なくとも部分的に、それと関連付けられるにエコー時間に基づいて、推定される欠落情報を算出的に更新し、（ｉｉｉ）少なくとも部分的に、更新された推定される欠落情報、それに対応する第２の複数のｋ空間画像のうちの少なくとも１つ、および第２の複数のｋ空間画像のうちの少なくとも１つに対応するベースラインｋ空間画像に基づいて、マップを発生させるように構成される。算出ユニットはさらに、少なくとも部分的に、デフォルト値、対応するベースラインｋ空間画像のうちの１つの一部、または第２の複数のｋ空間画像のうちの異なるものに基づいて、欠落情報を算出的に推定するように構成されてもよい。

【0014】

いくつかの実施形態では、マップは、標的領域の解剖学的画像または熱マップのうちの少なくとも一方を含む。ベースラインｋ空間画像はそれぞれ、複数のピクセルを備えてもよく、算出ユニットはさらに、少なくとも部分的に、ベースラインｋ空間画像のうちの１つの中のピクセルと関連付けられる大きさに基づいて、推定される欠落情報を更新するように構成されてもよい。

【0015】

本システムは、第１および第２のＭＲＩコイルを含んでもよく、ベースラインｋ空間画像のそれぞれおよび第２の複数のｋ空間画像のそれぞれは、第１および／または第２のＭＲＩコイルと関連付けられる。いくつかの実施形態では、算出ユニットはさらに、（ａ）複数のベースラインｋ空間画像から複数のベースライン画像を算出的に再構築し、（ｂ）第２の複数のｋ空間画像のうちのそれぞれ１つに関して、少なくとも部分的に、それと関連付けられる算出的に推定される欠落情報に基づいて、第３の画像セット内の画像を算出的に再構築するように構成される。

【0016】

種々の実施形態では、ベースラインｋ空間画像はそれぞれ、第１のマルチエコーパルスシーケンスに応答してエコーから入手され、第２の複数のｋ空間画像はそれぞれ、第２のマルチエコーパルスシーケンスに応答してエコーから入手される。算出ユニットはさらに、少なくとも部分的に、それと関連付けられるエコーのエコー時間に基づいて、第３の画像セット内の各再構築された画像を再構築されたベースライン画像のうちの１つに対応させるように構成されてもよい。

【0017】

本システムはさらに、複数のＭＲＩコイルを備えてもよく、ベースラインｋ空間画像のそれぞれおよび第２の複数のｋ空間画像のそれぞれは、ＭＲＩコイルのうちの少なくとも１つと関連付けられる。算出ユニットはさらに、少なくとも部分的に、関連付けられる少なくとも１つのＭＲＩコイルに基づいて、第３の画像セット内の各再構築された画像を再構築されたベースライン画像のうちの１つに関連させるように構成されてもよい。第３の画像セット内の再構築された画像と関連付けられる、エコー時間は、対応する再構築されたベースライン画像と関連付けられる、エコー時間と同一であってもよい。

【0018】

いくつかの実施形態では、算出ユニットはさらに、再構築されたベースライン画像のうちの１つの中のピクセルと第３の画像セット内の対応する再構築された画像の中の対応するピクセルとの間の位相差を決定するように構成される。算出ユニットはさらに、少なくとも部分的に、それから第３の画像セット内の対応する再構築された画像が再構築される、ｋ空間画像のエコー時間に基づいて、位相差を更新するように構成されてもよい。第３の画像セット内の再構築された画像と関連付けられる、位相差は、第２の複数のｋ空間画像のエコー時間と正に相関し得る。

【0019】

種々の実施形態では、算出ユニットはさらに、少なくとも部分的に、第２の複数のｋ空間画像内の少なくとも２つの異なる画像と関連付けられる、エコー時間に基づいて、位相差を更新するように構成される。代替として、または加えて、算出ユニットはさらに、少なくとも部分的に、更新された位相差に基づいて、第３の画像セット内の対応する再構築された画像を更新し、更新された画像を第４の画像セット内のｋ空間画像に変換するように構成されてもよい。なおもさらに、算出ユニットは、少なくとも部分的に、同一のエコー時間において入手される第２の複数のｋ空間画像内の対応する画像に基づいて、第４の画像セット内のｋ空間画像を更新するように構成されてもよい。

【0020】

いくつかの実施形態では、算出ユニットは、終了条件が満たされるまで、種々のステップを反復的に実施するように構成される。終了条件は、所定の限界値を超過する、反復の回数、または２回の反復の間の、第４の画像セット内の更新されたｋ空間画像または第３の画像セット内の再構築された画像の、所定の最小値を下回る変化のうちの１つ以上のものであってもよい。算出ユニットは、種々のステップを実施するために、人工ニューラルネットワーク（例えば、畳み込みまたは再帰型ニューラルネットワーク）等の機械学習アルゴリズムを実行するように構成されてもよい。

【0021】

いくつかの実施形態では、算出ユニットは、少なくとも部分的に、第４の画像セット内の更新されたｋ空間画像に基づいて、第５の画像セット内の画像を算出的に再構築し、第５の画像セット内の再構築された画像からノイズを算出的に低減させるように構成される。算出ユニットは、第３の画像セットまたは第５の画像セットのうちの少なくとも一方の中の画像のうちの少なくとも２つに、その中のノイズを低減させるために、局所的に低ランクの正則化を適用するように構成されてもよい。

【0022】

上記に説明される動作を実装する方法もまた、本発明の範囲内である。

【0023】

本明細書で使用されるように、用語「ほぼ」、「概ね」、および「実質的に」は、±１０％、いくつかの実施形態では、±５％を意味する。本明細書の全体を通して、「一実施例」、「ある実施例」、「一実施形態」、または「ある実施形態」の言及は、実施例と関連して説明される特定の特徴、構造、または特性が、本技術の少なくとも１つの実施例に含まれることを意味する。したがって、本明細書内全体を通した種々の場所における語句「一実施例では」、「ある実施例では」、「一実施形態」、または「ある実施形態」の表出は、必ずしも全てが同一の実施例を参照しているわけではない。さらに、特定の特徴、構造、ルーチン、ステップ、または特性は、本技術の１つ以上の実施例において、任意の好適な様式で組み合わせられてもよい。本明細書に提供される見出しは、便宜のためにすぎず、請求される技術の範囲または趣意を限定または解釈することを意図していない。
（項目１）
撮像システムであって、
磁気共鳴映像法（ＭＲＩ）装置であって、上記磁気共鳴映像法（ＭＲＩ）装置は、
（ｉ）第１のマルチエコーパルスシーケンスを実行し、標的領域を励起することと、
（ｉｉ）上記標的領域の複数のベースラインｋ空間画像を入手することであって、各ベースライン画像は、上記第１のマルチエコーパルスシーケンス内のエコーと関連付けられる、ことと、
（ｉｉｉ）第２のマルチエコーパルスシーケンスを実行し、上記標的領域を励起することと、
（ｉｖ）第２の複数のｋ空間画像を入手することであって、第２の複数のｋ空間画像はそれぞれ、上記第２のマルチエコーパルスシーケンス内のエコーと関連付けられ、上記第２の複数のｋ空間画像のうちの少なくとも１つは、上記ベースラインｋ空間画像と比較して情報が欠落するように、アンダーサンプリングされる、ことと
を行うように構成される、磁気共鳴映像法（ＭＲＩ）装置と、
算出ユニットであって、上記算出ユニットは、少なくとも部分的に、少なくとも１つのベースラインｋ空間画像および上記第２の複数のｋ空間画像のうちの少なくとも１つに基づいて、標的画像を再構築するように構成される、算出ユニットと
を備える、システム。
（項目２）
上記標的画像は、熱マップ、磁化率強調画像、Ｔ _２ ^＊ 画像、微小出血画像、またはＢＯＬＤｆ－ＭＲＩ画像のうちの１つである、項目１に記載のシステム。
（項目３）
上記算出ユニットは、ニューラルネットワークを備える、項目１に記載のシステム。
（項目４）
上記ニューラルネットワークは、畳み込みニューラルネットワークである、項目３に記載のシステム。
（項目５）
上記算出ユニットは、
上記第２の複数のｋ空間画像のうちの少なくとも１つに関して、（ｉ）上記ニューラルネットワークを使用して、それからの欠落情報を算出的に推定することと、（ｉｉ）少なくとも部分的に、上記推定される欠落情報に基づいて、標的画像を発生させることと
を行うように構成される、項目３に記載のシステム。
（項目６）
上記算出ユニットは、
上記第２の複数のｋ空間画像のうちの少なくとも１つに関して、（ｉ）上記ニューラルネットワークを使用して、対応する解剖学的画像を算出的に推定することと、（ｉｉ）少なくとも部分的に、上記推定される解剖学的画像に基づいて、標的画像を発生させることと
を行うように構成される、項目３に記載のシステム。
（項目７）
上記算出ユニットは、
上記第２の複数のｋ空間画像のうちの少なくとも１つに関して、（ｉ）それからの欠落情報を推定し、対応する解剖学的画像を再構築することと、（ｉｉ）上記ニューラルネットワークを使用して、上記解剖学的画像を算出的に補正することと、（ｉｉｉ）少なくとも部分的に、上記補正された解剖学的画像に基づいて、標的画像を発生させることと
を行うように構成される、項目３に記載のシステム。
（項目８）
上記算出ユニットは、
上記第２の複数のｋ空間画像のうちの少なくとも１つに関して、（ｉ）それからの欠落情報を算出的に推定することと、（ｉｉ）少なくとも部分的に、それと関連付けられるにエコー時間に基づいて、上記推定される欠落情報を算出的に更新することと、（ｉｉｉ）少なくとも部分的に、上記更新された推定される欠落情報、それに対応する上記第２の複数のｋ空間画像のうちの少なくとも１つ、および上記第２の複数のｋ空間画像のうちの少なくとも１つに対応する上記ベースラインｋ空間画像に基づいて、マップを発生させることと
を行うように構成される、項目１に記載のシステム。
（項目９）
上記マップは、上記標的領域の解剖学的画像または熱マップのうちの少なくとも一方を含む、項目８に記載のシステム。
（項目１０）
上記算出ユニットはさらに、少なくとも部分的に、デフォルト値、上記対応するベースラインｋ空間画像のうちの１つの一部、または上記第２の複数のｋ空間画像のうちの異なるものに基づいて、上記欠落情報を算出的に推定するように構成される、項目８に記載のシステム。
（項目１１）
上記ベースラインｋ空間画像のうちのそれぞれ１つは、複数のピクセルを備え、上記算出ユニットはさらに、少なくとも部分的に、上記ベースラインｋ空間画像のうちの１つの中のピクセルと関連付けられる大きさに基づいて、上記推定される欠落情報を更新するように構成される、項目８に記載のシステム。
（項目１２）
第１および第２のＭＲＩコイルをさらに備え、上記ベースラインｋ空間画像のそれぞれおよび上記第２の複数のｋ空間画像のそれぞれは、上記第１または第２のＭＲＩコイルのうちの少なくとも一方と関連付けられる、項目１に記載のシステム。
（項目１３）
上記算出ユニットはさらに、（ａ）上記複数のベースラインｋ空間画像から複数のベースライン画像を算出的に再構築することと、（ｂ）上記第２の複数のｋ空間画像のうちのそれぞれ１つに関して、少なくとも部分的に、それと関連付けられる上記算出的に推定される欠落情報に基づいて、第３の画像セット内の画像を算出的に再構築することとを行うように構成される、項目８に記載のシステム。
（項目１４）
上記ベースラインｋ空間画像のうちのそれぞれ１つは、上記第１のマルチエコーパルスシーケンスに応答してエコーから入手され、上記第２の複数のｋ空間画像のうちのそれぞれ１つは、上記第２のマルチエコーパルスシーケンスに応答してエコーから入手され、上記算出ユニットはさらに、少なくとも部分的に、それと関連付けられる上記エコーのエコー時間に基づいて、上記第３の画像セット内の各再構築された画像を上記再構築されたベースライン画像のうちの１つに対応させるように構成される、項目１３に記載のシステム。
（項目１５）
複数のＭＲＩコイルをさらに備え、上記ベースラインｋ空間画像のそれぞれおよび上記第２の複数のｋ空間画像のそれぞれは、上記ＭＲＩコイルのうちの少なくとも１つと関連付けられ、上記算出ユニットはさらに、少なくとも部分的に、上記関連付けられる少なくとも１つのＭＲＩコイルに基づいて、上記第３の画像セット内の各再構築された画像を上記再構築されたベースライン画像のうちの１つに関連させるように構成される、項目１４に記載のシステム。
（項目１６）
上記第３の画像セット内の上記再構築された画像と関連付けられる上記エコー時間は、上記対応する再構築されたベースライン画像と関連付けられる上記エコー時間と同一である、項目１４に記載のシステム。
（項目１７）
上記算出ユニットはさらに、（ｃ）上記再構築されたベースライン画像のうちの１つの中のピクセルと上記第３の画像セット内の上記対応する再構築された画像の中の対応するピクセルとの間の位相差を決定するように構成される、項目１４に記載のシステム。
（項目１８）
上記算出ユニットはさらに、（ｄ）少なくとも部分的に、それから上記第３の画像セット内の上記対応する再構築された画像が再構築される上記ｋ空間画像のエコー時間に基づいて、上記位相差を更新するように構成される、項目１７に記載のシステム。
（項目１９）
上記第３の画像セット内の上記再構築された画像と関連付けられる上記位相差は、上記第２の複数のｋ空間画像のエコー時間と正に相関する、項目１８に記載のシステム。
（項目２０）
上記算出ユニットはさらに、少なくとも部分的に、上記第２の複数のｋ空間画像内の少なくとも２つの異なる画像と関連付けられる上記エコー時間に基づいて、上記位相差を更新するように構成される、項目１８に記載のシステム。
（項目２１）
上記算出ユニットはさらに、（ｅ）少なくとも部分的に、上記更新された位相差に基づいて、上記第３の画像セット内の上記対応する再構築された画像を更新することと、（ｆ）上記更新された画像を第４の画像セット内のｋ空間画像に変換することとを行うように構成される、項目１８に記載のシステム。
（項目２２）
上記算出ユニットはさらに、（ｇ）少なくとも部分的に、同一のエコー時間において入手される上記第２の複数のｋ空間画像内の上記対応する画像に基づいて、上記第４の画像セット内の上記ｋ空間画像を更新するように構成される、項目１７に記載のシステム。
（項目２３）
上記算出ユニットはさらに、終了条件が満たされるまで、ステップ（ｂ）－（ｇ）を反復的に実施するように構成される、項目２２に記載のシステム。
（項目２４）
上記終了条件は、
所定の限界値を超過する反復の回数、または
２回の反復の間の上記第４の画像セット内の上記更新されたｋ空間画像または上記第３の画像セット内の上記再構築された画像の所定の最小値を下回る変化
のうちの１つ以上のものに対応する、項目２３に記載のシステム。
（項目２５）
上記算出ユニットはさらに、ステップ（ｂ）－（ｇ）を少なくとも２回実施するために、人工ニューラルネットワークを実行するように構成される、項目２２に記載のシステム。
（項目２６）
上記算出ユニットはさらに、少なくとも部分的に、上記第４の画像セット内の上記更新されたｋ空間画像に基づいて、第５の画像セット内の画像を算出的に再構築し、上記第５の画像セット内の上記再構築された画像からノイズを算出的に低減させるように構成される、項目２２に記載のシステム。
（項目２７）
上記算出ユニットはさらに、上記第３の画像セットまたは上記第５の画像セットのうちの少なくとも一方の中の上記画像のうちの少なくとも２つに、その中の上記ノイズを低減させるために、局所的に低ランクの正則化を適用するように構成される、項目２６に記載のシステム。
（項目２８）
磁気共鳴映像法の方法であって、上記方法は、
第１のマルチエコーパルスシーケンスを実行し、標的領域を励起することと、
上記標的領域の複数のベースラインｋ空間画像を入手することであって、各ベースライン画像は、上記第１のマルチエコーパルスシーケンス内のエコーと関連付けられる、ことと、
第２のマルチエコーパルスシーケンスを実行し、上記標的領域を励起することと、
第２の複数のｋ空間画像を入手することであって、上記第２の複数のｋ空間画像はそれぞれ、上記第２のマルチエコーパルスシーケンス内のエコーと関連付けられ、上記第２の複数のｋ空間画像のうちの少なくとも１つは、上記ベースラインｋ空間画像と比較して情報が欠落するように、アンダーサンプリングされる、ことと、
少なくとも部分的に、少なくとも１つのベースラインｋ空間画像および上記第２の複数のｋ空間画像のうちの少なくとも１つに基づいて、標的画像を再構築することと
を含む、方法。
（項目２９）
上記標的画像は、熱マップ、磁化率強調画像、Ｔ _２ ^＊ 画像、微小出血画像、またはＢＯＬＤ ｆ－ＭＲＩ画像のうちの１つである、項目２８に記載の方法。
（項目３０）
上記第２の複数のｋ空間画像のうちの少なくとも１つに関して、（ｉ）それからの欠落情報を算出的に推定することと、（ｉｉ）少なくとも部分的に、上記推定される欠落情報に基づいて、標的画像を発生させることとをさらに含む、項目２８に記載の方法。
（項目３１）
上記第２の複数のｋ空間画像のうちの少なくとも１つに関して、（ｉ）対応する解剖学的画像を算出的に推定することと、（ｉｉ）少なくとも部分的に、上記推定される解剖学的画像に基づいて、標的画像を発生させることとをさらに含む、項目２８に記載の方法。
（項目３２）
上記第２の複数のｋ空間画像のうちの少なくとも１つに関して、（ｉ）それからの欠落情報を推定し、対応する解剖学的画像を再構築することと、（ｉｉ）上記解剖学的画像を算出的に補正することと、（ｉｉｉ）少なくとも部分的に、上記補正された解剖学的画像に基づいて、標的画像を発生させることとをさらに含む、項目２８に記載の方法。
（項目３３）
上記第２の複数のｋ空間画像のうちの少なくとも１つに関して、（ｉ）それからの欠落情報を算出的に推定することと、（ｉｉ）少なくとも部分的に、それと関連付けられるエコー時間に基づいて、上記推定される欠落情報を算出的に更新することと、（ｉｉｉ）少なくとも部分的に、上記更新された推定される欠落情報、それに対応する上記第２の複数のｋ空間画像のうちの少なくとも１つ、および上記第２の複数のｋ空間画像のうちの少なくとも１つに対応する上記ベースラインｋ空間画像に基づいて、マップを発生させることとをさらに含む、項目２８に記載の方法。
（項目３４）
上記マップは、上記標的領域の解剖学的画像または熱マップのうちの少なくとも一方を含む、項目３３に記載の方法。
（項目３５）
上記欠落情報は、少なくとも部分的に、デフォルト値、上記対応するベースラインｋ空間画像のうちの１つの一部、または上記第２の複数のｋ空間画像のうちの異なるものに基づいて、算出的に推定される、項目３３に記載の方法。
（項目３６）
上記ベースラインｋ空間画像のうちのそれぞれ１つは、複数のピクセルを備え、上記推定される欠落情報は、少なくとも部分的に、上記ベースラインｋ空間画像のうちの１つの中のピクセルと関連付けられる大きさに基づいて更新される、項目３３に記載の方法。
（項目３７）
（ａ）上記複数のベースラインｋ空間画像から複数のベースライン画像を算出的に再構築することと、（ｂ）上記第２の複数のｋ空間画像のうちのそれぞれ１つに関して、少なくとも部分的に、それと関連付けられる上記算出的に推定される欠落情報に基づいて、第３の画像セット内の画像を算出的に再構築することとをさらに含む、項目３３に記載の方法。
（項目３８）
上記ベースラインｋ空間画像のうちのそれぞれ１つは、上記第１のマルチエコーパルスシーケンスに応答してエコーから入手され、上記第２の複数のｋ空間画像のうちのそれぞれ１つは、上記第２のマルチエコーパルスシーケンスに応答してエコーから入手され、少なくとも部分的に、それと関連付けられる上記エコーのエコー時間に基づいて、上記第３の画像セット内の各再構築された画像を上記再構築されたベースライン画像のうちの１つに対応させることをさらに含む、項目３７に記載の方法。
（項目３９）
上記ベースラインｋ空間画像のそれぞれおよび上記第２の複数のｋ空間画像のそれぞれは、複数のＭＲＩコイルのうちの少なくとも１つと関連付けられ、少なくとも部分的に、上記関連付けられる少なくとも１つのＭＲＩコイルに基づいて、上記第３の画像セット内の各再構築された画像を上記再構築されたベースライン画像のうちの１つに関連させることをさらに含む、項目３８に記載の方法。
（項目４０）
上記第３の画像セット内の上記再構築された画像と関連付けられる上記エコー時間は、上記対応する再構築されたベースライン画像と関連付けられる上記エコー時間と同一である、項目３８に記載の方法。
（項目４１）
（ｃ）上記再構築されたベースライン画像のうちの１つの中のピクセルと上記第３の画像セット内の上記対応する再構築された画像の中の対応するピクセルとの間の位相差を決定することをさらに含む、項目３８に記載の方法。
（項目４２）
（ｄ）少なくとも部分的に、それから上記第３の画像セット内の上記対応する再構築された画像が再構築される上記ｋ空間画像のエコー時間に基づいて、上記位相差を更新することをさらに含む、項目４１に記載の方法。
（項目４３）
上記第３の画像セット内の上記再構築された画像と関連付けられる上記位相差は、上記第２の複数のｋ空間画像のエコー時間と正に相関する、項目４２に記載の方法。
（項目４４）
上記位相差は、少なくとも部分的に、上記第２の複数のｋ空間画像内の少なくとも２つの異なる画像と関連付けられる上記エコー時間に基づいて更新される、項目４２に記載の方法。
（項目４５）
（ｅ）少なくとも部分的に、上記更新された位相差に基づいて、上記第３の画像セット内の上記対応する再構築された画像を更新することと、（ｆ）上記更新された画像を第４の画像セット内のｋ空間画像に変換することとをさらに含む、項目４２に記載の方法。
（項目４６）
上記算出ユニットはさらに、（ｇ）少なくとも部分的に、同一のエコー時間において入手される上記第２の複数のｋ空間画像内の上記対応する画像に基づいて、上記第４の画像セット内の上記ｋ空間画像を更新するように構成される、項目４１に記載の方法。
（項目４７）
ステップ（ｂ）－（ｇ）が、終了条件が満たされるまで反復的に実施される、項目４６に記載の方法。
（項目４８）
上記終了条件は、
所定の限界値を超過する反復の回数、または
２回の反復の間の上記第４の画像セット内の上記更新されたｋ空間画像または上記第３の画像セット内の上記再構築された画像の所定の最小値を下回る変化
のうちの１つ以上のものに対応する、項目４７に記載の方法。
（項目４９）
ステップ（ｂ）－（ｇ）が、少なくとも２回実施される、項目４６に記載の方法。
（項目５０）
少なくとも部分的に、上記第４の画像セット内の上記更新されたｋ空間画像に基づいて、第５の画像セット内の画像を算出的に再構築することと、上記第５の画像セット内の上記再構築された画像からノイズを算出的に低減させることとをさらに含む、項目４６に記載の方法。
（項目５１）
上記第３の画像セットまたは上記第５の画像セットのうちの少なくとも一方の中の上記画像のうちの少なくとも２つに、その中の上記ノイズを低減させるために、局所的に低ランクの正則化を適用することをさらに含む、項目５０に記載の方法。

【図面の簡単な説明】

【0024】

図面では、同様の参照記号は、概して、異なる図全体を通して同一の部分を指す。また、図面は、必ずしも縮尺通りではなく、代わりに、概して、本発明の原理を図示することに重点が置かれている。以下の説明では、本発明の種々の実施形態が、以下の図面を参照して説明される。

【0025】

【図1】図１は、代表的なＭＲＩ装置を断面的に描写する。

【0026】

【図2A】図２Ａは、ＰＲＦ温度測定のための例示的な単一エコーＧＲＥパルスシーケンスを描写する。

【0027】

【図2B】図２Ｂは、エコートレイン内のエコーが入手される、マルチエコーＧＲＥパルスシーケンスを描写する。

【0028】

【図2C】図２Ｃは、マルチエコーＧＲＥパルスシーケンスを使用して入手される、画像スライスの完全にサンプリングされたデータセットを描写する。

【0029】

【図3】図３は、種々の実施形態による、加速されたＭＲ温度測定への例示的アプローチを図示する、フローチャートである。

【0030】

【図4A】図４Ａは、あるデータが欠落している、マルチエコーＧＲＥパルスシーケンスを使用して入手される、画像スライスのデータセットを描写する。

【0031】

【図4B】図４Ｂは、種々の実施形態による、ｋ空間データをアンダーサンプリングするための例示的アプローチを描写する。

【0032】

【図4C】図４Ｃは、図４Ｂに示されるような１つのＴＲ内に複数の位相エンコード値を有する、エコーを入手するために実装される、Ｇ_ｙ勾配ブリップのセットを描写する。

【0033】

【図4D】図４Ｄは、いくつかの実施形態において有用である、別のｋ空間データ軌道を図示する。

【0034】

【図5A】図５Ａは、種々の実施形態による、反復的推定プロシージャを実施するための例示的アプローチを表す、フローチャートである。

【0035】

【図5B】図５Ｂは、図５Ａに記載される推定アプローチを使用した、反復間での画質の改良を描写する。

【0036】

【図6】図６は、それらの間に周波数差を有する、２つのＲＦパルスを印加し、異なる位相を有する２つのスライスを同時に入手するための例示的アプローチをグラフで描写する。

【0037】

【図7】図７Ａ－７Ｃは、種々の実施形態による、ｋ空間データが渦巻状の様式において入手され得る方法を図示する。

【発明を実施するための形態】

【0038】

本発明は、画像誘導下の治療／診断手技の間にリアルタイムで標的（例えば、治療標的）または着目領域の中の他の着目物体を監視するためのシステムおよび方法を提供する。手技は、例えば、癌性である場合には、組織を壊死させる、アブレートする、または別様に破壊することのいずれかを行うように加熱する目的のために、または疼痛の改善または温熱療法の制御された誘引等の非破壊的治療のために、例えば、物質、組織、または器官（すなわち、その超音波処理）への集束された超音波の印加を伴い得る。超音波はまた、例えば、神経変調等の他の非熱的タイプの治療のために使用され得る。代替として、手技は、例えば、高周波（ＲＦ）放射線、Ｘ線またはガンマ線、または荷電粒子等の異なる形態の療法エネルギーを使用する、または凍結アブレーション等の他の治療モダリティを伴い得る。種々の治療手技における標的の温度および／または場所を監視するステップが、標的への、および／または他の非標的組織および器官の周囲に療法エネルギービームを誘導する役割を果たし、すなわち、影響を及ぼされる解剖学的領域の画像に基づいて、ビームの焦点、プロファイル、および／または方向を調節し得、これはまた、いくつかの実施形態では、ビーム焦点を可視化させ得る。ＭＲＩは、そのような画像ベースの熱的および／または運動追跡のために広く使用される技法である。しかしながら、例えば、Ｘ線撮像、Ｘ線コンピュータ断層撮影（ＣＴ）、または超音波撮像を含む、他の撮像技法もまた、使用され得、本発明の範囲内である。加えて、温度および／または運動監視は、１つ以上の２次元画像および／または３次元画像を使用して達成され得る。本明細書に説明される技法が実装され得るＭＲＩシステムは、当技術分野において周知であり、例示的システムが、図１に示される。

【0039】

いくつかの実施形態では、手技の間の撮像が、生体内温度を定量的に監視するために同時に使用される。これは、特に、ＭＲ誘導下の熱療法（例えば、ＭＲｇＦＵＳ治療）において有用であり、標的治療面積（例えば、熱によって破壊されるべき腫瘍）の温度は、治療の進行を査定し、熱伝導およびエネルギー吸収の局所的な差異を補正し、治療面積を囲繞する組織への損傷を回避するために、持続的に監視されるべきである。ＭＲ撮像を用いた温度の監視（例えば、測定および／またはマッピング）は、概して、ＭＲ温度測定またはＭＲ熱撮像と称される。

【0040】

ＭＲ温度測定のために利用可能な種々の方法の中でも、陽子共鳴周波数（ＰＲＦ）偏移法は、多くの場合、温度変化に対するその優れた線形性、組織タイプからの近独立性、および高空間および時間分解能を伴う温度マップ入手に起因して、選定の方法である。ＰＲＦ偏移法は、水分子中の陽子のＭＲ共鳴周波数が、（有利には、組織タイプ間で比較的に一定である、比例定数で）温度とともに線形に変化する現象に基づく。温度に伴う周波数変化が、小さい、すなわち、バルク水に関しては、－０．０１ｐｐｍ／℃、組織では、およそ－０．００９６～－０．０１３ｐｐｍ／℃のみであるため、ＰＲＦ偏移は、典型的には、撮像が２回、すなわち、最初に、温度変化に先立ってベースラインＰＲＦ位相画像を入手し、次いで、温度変化の後に第２の位相画像、すなわち、治療画像を入手し、それによって、温度の変化に比例するわずかな相変化を捕捉するために実施される、位相感受性の撮像方法を用いて検出される。温度変化のマップが、次いで、ピクセル毎に、ベースライン画像と治療画像との間の位相差を決定することによって（再構築された）画像から算出され、位相差を、静磁場の強度および（例えば、勾配再呼出エコーの）エコー時間（ＴＥ）等の撮像パラメータを考慮しながらＰＲＦ温度依存度に基づいて温度差に転換し得る。図２Ａは、１つのみのエコー２０２が、連続パルスシーケンス（すなわち、繰り返し時間またはＴＲ）間で（ＴＥにおいて）取得される、ＰＲＦ温度測定のための例示的単一エコーＧＲＥパルスシーケンスを描写する。

【0041】

単一エコーＧＲＥパルスシーケンスを使用して入手される、典型的ＭＲ画像では、場の不均質性によって引き起こされる、空間偏移Δｒが、存在し得る。ＭＲ受信機の帯域幅を増大させることは、空間偏移を低減させる傾向にあるが、温度の標準偏差（または測定誤差）を増大させる代償を伴う。結果として、帯域幅の選定は、空間偏移を最小限にすることと測定誤差を低減させることとの間のトレードオフを反映する。本課題に対処するために、いくつかの実施形態では、ＭＲＩシステムは、広い帯域幅を伴う受信機を使用し、ＰＲＦ温度測定のためにマルチエコーＧＲＥパルスシーケンスを印加する。例えば、図２Ｂを参照すると、マルチエコーＧＲＥパルスシーケンスでは、エコートレイン内のＮ個のエコー（描写されるように、Ｎ＝６である）が、各ＴＲ内で入手され得る。各エコーは、広い受信機帯域幅を伴う比較的に短い時間窓内で入手されるため、空間偏移Δｒは、ごくわずかであり得る。加えて、Ｎ個のエコーの全てから受信される信号を結合する（または平均化する）ことによって、受信されるＭＲ信号の信号対雑音比（ＳＮＲ）が、改良され得、これは、ひいては、温度の測定誤差を低減させる。図２Ｃは、マルチエコーＧＲＥパルスシーケンスを使用して入手される撮像スライスの完全に（または「全面的に」）サンプリングされたデータセットを描写する。描写されるように、各行は、１つのＴＲ内で入手される、位相エンコーディングステップｋ_ｙと関連付けられ（および、平面に対して垂直であり、したがって、示されていない、ｋ_ｘが、ＭＲ信号読出の数と関連付けられる）、加えて、複数のエコーＮが、各ＴＲ内で入手され、各エコーＥ_ｉ（ｉ≦Ｎ）が、ＴＥ_ｉの時点において入手される。各撮像スライスが、Ｎ_１個の位相エンコーディングステップと、Ｎ_２回の読出とを含むと仮定して、スライスの完全にサンプリングされるデータセットは、Ｎ_１×Ｎ_２×Ｎ個のデータ点を含む。神経系に関する従来のＭＲ熱走査では、ＴＲ＝２７ミリ秒およびＮ_１＝１２８であり、したがって、これは、ＭＲ画像の単一のスライスを入手するために、最大３．５秒を要し得る。本走査率は、超音波手技を開始することに先立ってベースラインＰＲＦ位相画像を入手するためには許容可能であり得るが、これは、（概して、標的治療領域と一致する）超音波焦点における温度が迅速に上昇するときは、治療の進行を査定し、超音波手技の間にリアルタイムで標的領域を囲繞する組織への損傷を回避するためには、遅すぎる場合がある。本明細書における、用語「完全にサンプリングされるデータセット」は、完全にサンプリングされるデータセットだけではなく、そのために十分なデータセットが発生または推定され得る、アンダーサンプリングされるデータセットも指す。

【0042】

故に、図２Ｃに描写される完全にサンプリングされたｋ空間データセットは、超音波手技に先立ってのみ入手され得、超音波手技の間、部分的な未加工のｋ空間データのみが、入手され、下記にさらに説明されるように、リアルタイムの（または少なくとも加速された）ＭＲ温度測定を提供する。図３は、種々の実施形態による、加速されたＭＲ温度測定への例示的アプローチ３００を図示する。本方法は、着目超音波手技に先立って実施される予備ステップ、および手技の間に実行されるステップを含む。第１の予備ステップ３０２では、１つ以上のＭＲマルチエコーＧＲＥパルスシーケンスが、超音波治療または診断の開始に先立って、着目領域（例えば、治療領域）の入手された複数の完全にサンプリングされたｋ空間画像に対して順に発生され、各画像は、エコー時間ＴＥ_ｉにおいて受信される、エコーＥ_ｉに対応し、異なる画像が、同一または異なるＭＲコイル１０４を使用して入手されてもよい。第２の予備ステップ３０４では、各ｋ空間画像が、逆フーリエ変換され、画像を再構築し、これは、次いで、ベースライン基準画像として利用されることができる。治療／診断の開始後、（ステップ３０６において）１つ以上のＭＲマルチエコーＧＲＥパルスシーケンスが、発生され、着目領域の複数の画像を入手してもよく、再び、各治療画像は、時間ＴＥ_ｉにおいて入手される、エコーＥ_ｉに対応する。１つの実装では、各治療画像は、ベースライン基準画像のものに対応する、エコー時間ＴＥ_ｉを有する。加えて、部分的な未加工のｋ空間データのみが、データ入手時間を短縮させ、それによって、着目領域（例えば、治療領域）内の組織状態（例えば、温度）のリアルタイムの監視を可能にするために、治療画像のエコーのうちの少なくとも１つにおいて入手される。（部分的な未加工のｋ空間データは、ステップ３０２において行われるように、ｋ空間画像データ全体未満を取得するための「アンダーサンプリング」によって入手されてもよい。）例えば、図２Ｃに示される完全にサンプリングされたｋ空間データでは、各ＴＲ内の各エコーが、ｋ_ｙのＮ_１個の値の全てにおいてサンプリングされる。データをアンダーサンプリングするとき、全てのエコーが、各ＴＲの間にサンプリングされるが、異なる位相エンコード値が、欠落している。例えば、図４Ａでは、全６個のエコーが、各ＴＲ内でサンプリングされるが、奇数のｋ_ｙ値のみが、エコー４０２によってサンプリングされ、偶数のｋ_ｙ値のみが、エコー４０８によってサンプリングされる。図４Ａのｋ_ｙエコーマトリクスは、半分空であるため、ＴＲの数（および、したがって、走査時間）は、図２Ｃに描写される、完全にサンプリングされたｋ_ｙエコーマトリクスと比較して、係数Ｒ＝２だけ低減される。

【0043】

図４Ｂは、種々の実施形態による、ｋ空間データをアンダーサンプリングするための別の例示的アプローチを描写する。描写されるように、各ＴＲ４１０内のＮ個のエコーが、それらが全て、異なる位相エンコード値を有し、異なるＴＥ時間に対応するように、交互配置される様式において入手される。本アプローチは、したがって、さらに、データ入手時間を６倍短縮させる（例えば、１８個の位相エンコード値を有する、３つのみのＴＲが、データを入手するために要求される）。図４Ｃを参照すると、種々の実施形態では、Ｇ_ｙ勾配ブリップのセットが、実装され、（図４Ｂに描写されるように）１つのＴＲの中に複数の位相エンコード値を有する、エコーを入手する。示されるように、各勾配ブリップは、２つのエコーの入手時間の間に導入されてもよい。

【0044】

したがって、上記に説明されるアプローチを利用するＭＲ走査は、撮像スライスのデータセット毎にＮ_１／Ｒ（Ｒは、撮像加速係数を表し、それぞれ、図４Ａおよび４Ｂの２および６に等しい）個のみの位相エンコード線の入手を要求し得、これは、Ｒ個のスライスが、３秒以内に入手されることを可能にし、事実上、繰り返し時間をＴＲ×Ｒミリ秒まで増大させる。ＴＲ＝２７ミリ秒およびＮ_１＝１２８である、神経系のための従来のＭＲ熱走査では、図４Ｂに描写されるようなｋ空間データをアンダーサンプリングするステップは、有利なこととして、ＭＲ画像の６個のスライスが、１６２ミリ秒のＴＲを伴って、３秒以内に入手されることを可能にし（各スライスは、２１個の位相エンコード線を有する）、対照的に、従来のアプローチでは、１つのみのスライスが、３秒以内に入手され得る。加えて、上記に説明されるアンダーサンプリングアプローチを使用してＴＲを増大させることは、より低いＴ_１信号の飽和に起因して、身体組織（例えば、脳）からのＭＲ信号の品質を有意に改良し得る。例えば、脳のＭＲ信号のＳＮＲは、（従来のＭＲ走査における）２７ミリ秒から（図４Ｂの）１６２ミリ秒までＴＲを増大させることによって、１．９倍増大され得る。

【0045】

図４Ａおよび４Ｂに図示されるような、各ＴＲ４１０内に複数のエコーを含む、ｋ空間データ軌道が、例示的実施形態を表しているにすぎず、受信されたエコーの位相エンコード値およびＴＥ時間を選択することによって確立された他の軌道が、ある加速率でＭＲデータを入手するために好適であり得、したがって、本発明の範囲内であることに留意されたい。例えば、図４Ｄに描写される軌道が、いくつかの実施形態において使用されてもよい。代替として、各ＴＲ４１０内で入手されるＮ個のエコーは、無作為に選択されてもよい。加えて、各ＴＲ４１０内のデータ軌道は、別のＴＲ４１０内のものと同一または異なってもよい。

【0046】

再び図３を参照すると、種々の実施形態では、反復的推定プロシージャが、下記にさらに説明されるように、１つ以上の制約に基づいて、欠落しているｋ空間データを算出的に決定するために実装される（ステップ３０８）。図５Ａは、本明細書による、反復的推定プロシージャを実施するための例示的アプローチ５００を描写する。例えば、第１のステップ５０２において、欠落しているｋ空間データは、同一のエコートレイン内の１つ以上のエコーから入手されたデータ（例えば、それらの平均）に基づいて、最初に推定されてもよい。加えて、または代替として、欠落しているｋ空間データは、同一のエコー時間ＴＥ_ｉにおけるエコーから入手されるが、異なるエコートレインまたは異なる治療画像内、および／またはベースライン画像内のものであるｋ空間データに基づいて、最初に推定されてもよい。いくつかの実施形態では、欠落しているｋ空間データは、１つ以上のゼロまたは任意の恣意的数で最初に充填されることができる。入手された（またはサンプリングされた）ｋ空間データおよび推定されるデータは、次いで、ｋ空間データＫ_{ｔｒｅａｔｍｅｎｔ}（ｋ_ｙ，ｋ_ｘ，ｅｃｈｏ）の完全（または全面的）なセットを形成し、これは、（ステップ５０４において）逆フーリエ変換され、治療画像Ｘ_{ｔｒｅａｔｍｅｎｔ}（ｙ，ｘ，ｅｃｈｏ）を再構成することができる。

【化1】

【0047】

いくつかの実施形態では、反復的推定プロシージャは、次いで、欠落しているｋ空間データを正確に推定するように、データセット上に１つ以上の制約を課す。一実施形態では、制約は、治療画像についての以前の知識に基づく。大きさρ_ｉ（ｘ，ｙ）と、位相φ_ｉ（ｘ，ｙ）とを有する、ベースライン画像上のピクセル（ｘ，ｙ）におけるエコーＥ_ｉの信号Ｓ_{ｉ，ｂａｓｅｌｉｎｅ}は、以下のように表されることができる。

【化2】

温度が、治療の間に着目領域内で上昇する（ΔＴ）とき、治療画像のピクセル（ｘ，ｙ）におけるエコーＥ_ｉの信号Ｓ_{ｉ，ｔｒｅａｔｍｅｎｔ}は、以下のように表されることができる。

【化3】

式中、Ｋ＝０．０１ｐｐｍ／℃である。したがって、治療画像と対応するベースライン画像との間のエコーＥ_ｉの位相差（すなわち、ΔΦ_ｉ）は、エコーのエコー時間ＴＥ_ｉと正に（例えば、線形に）相関し得る。

【化4】

【0048】

故に、反復的推定プロシージャ５００によって課される第１の制約は、エコー時間ＴＥ_ｉと、治療画像および対応するベースライン画像間の位相差との間の線形関係であり得る。再び図５Ａを参照すると、一実施形態では、反復的推定プロシージャは、治療画像の中のピクセルの新しい位相差が、エコー時間ＴＥに対して線形であり、算出される位相差ΔΦ（ｘ，ｙ，ｅｃｈｏ）に可能な限り近接するように、（ステップ５０６において）最初に、ピクセル毎に、治療画像とその対応するベースライン画像との間の位相差ΔΦ（ｘ，ｙ，ｅｃｈｏ）を算出し、続いて、治療画像のピクセル（ｘ，ｙ）毎に、新しい位相差ΔΦ_ｎｅｗを推定する。

【化5】

式中、Ｆ_０およびＦ_１は、ＴＥから独立した、２つの係数である。例えば、ベースライン画像と関連付けられるエコーＥｊの信号は、以下のように表されることができる。

【化6】

治療画像と関連付けられるエコーＥｊの信号は、以下のように表されることができる。

【化7】

したがって、場所（ｘ，ｙ）毎に、ベクトルｄ_ｊが、エコーＥ_ｊ（ｊ＝１～Ｎ）に関して以下のように定義され得る。

【化8】

方程式（４）に基づいて、

【化9】

は、以下のように表されることができる。

【化10】

続いて、方程式（８）内のベクトルｄ_ｊに基づいて、Ｎ－１個の成分を伴うベクトルｖ_ｋを定義することができる。

【化11】

式中、ΔＴＥは、２つの隣接したエコー間のＴＥの増分を表す。方程式（１０）内の位相は、ｋから独立しているため、Ｆ_１は、以下のように表されることができる。

【化12】

式中、

【化13】

は、４象限アークタンジェント演算子を表し、「ｍｅａｎ」は、ｖの全ての成分にわたる平均値を表す。
加えて、ベクトルｄおよび既知のＦ_１に基づいて、Ｎ個の成分を伴う別のベクトルｕを定義することができる。ｕのｊ番目の成分ｕ_ｊが、以下のように求められる。

【化14】

ｕ_ｊの位相は、ｊから独立しているため、Ｆ_０は、以下のように表されることができる。

【化15】

方程式（１１）および（１２）に基づいて、２つの係数Ｆ_０およびＦ_１が、算出されることができる。下記にさらに説明されるように、Ｆ_０およびＦ_１の正確な値が、推定されるデータの合理的収束が達成された後、取得され得る。種々の実施形態では、対応するベースライン画像内の位相とともに、Ｆ_０およびＦ_１に基づいて推定される、新しい位相差が、次いで、（ステップ５０８において）ピクセル毎に治療画像の中のピクセルの位相を更新するために利用されてもよい。例えば、治療画像の中のエコーＥｊの信号は、以下のように更新され得る。

【化16】

【0049】

加えて、超音波手技から結果として生じる温度上昇が、主に、治療画像の中のピクセルの位相に影響を及ぼすため、反復的推定プロシージャは、（ステップ５１０において）治療画像の中の各ピクセルの大きさを、対応するベースライン画像の中のものと同一になるように設定する、第２の制約を課してもよい。

【化17】

加えて、または代替として、スピン密度ρ_{ｔｒｅａｔｍｅｎｔ}が、（体温を２０℃上回る温度までの加熱から結果として生じる、約２５％のＴ１の延長に起因して）加熱の間にρ_{ｂａｓｅｌｉｎｅ}に対して１０％～１５％減少するという事実が、考慮されることができる。故に、方程式（１３）内の治療の間のエコーｊの信号は、以下のように記述されることができる。

【化18】

式中、Ａ（１より小さい）は、信号振幅の本減少を表す、正の実数である。Ａは、ｘ、ｙに依存するが、エコー時間から独立している。方程式（６）および（１４）、およびＡが、全てのＴＥに関して同一であるという事実から、以下を結論付ける。

【化19】

式中、Ｓ_{ｊ，ｔｒｅａｔｍｅｎｔ}およびＳ_{ｊ，ｂｓｅｌｉｎｅ}は、エコーｊにおける治療信号およびベースライン信号であり、

【化20】

は、括弧の内側における、信号の大きさを表す。方程式（１４ａ）内の「ｍｅａｎ」は、Ａが、全てのエコーおよびコイルに関して同一であるため、エコーおよびコイルの全てにわたる平均値である。各反復の間、方程式（１４ａ）を使用してＡを計算し、次いで、これを方程式（１４ｂ）の中に代入する。Ａは、空間的にゆっくりと変動しつつあるため、標準的な全変動（ＴＶ）演算子が、これを平滑化させ、その正確度を改良するために使用されてもよい。方程式（１４ｂ）内の計算された信号Ｓ_{ｊ，ｔｒｅａｔｍｅｎｔ}が、図５Ａのステップ５１０において課される。反復の間、Ｆ_０、Ｆ_１、およびＡは、補正値に収束する。代替として、または加えて、共鳴周波数を変化させる、加熱によって引き起こされる、化学的偏移の影響が、振幅変化に加えて、またはその代わりに、方程式（１４ｂ）に統合されることができる。加えて、または代替として、モデルは、ベースライン画像および治療画像の両方に関して、「水脂肪分離」の周知の技法を使用することによって、脂肪を考慮するように拡張されることができる。

【0050】

上記の制約が、課された後、反復的推定プロシージャは、（ステップ５１２において）フーリエ変換を更新された治療画像に適用し、更新されたｋ空間データセットを発生させてもよい。いくつかの実施形態では、反復的推定プロシージャは、（ステップ５１４において）最初に入手されたアンダーサンプリングされたｋ空間データによって、更新されたｋ空間データセットを置換する、第３の制約を課す。

【化21】

式中、Ｋ_ｓａｍｐは、種々の場所において最初に入手されたｋ空間データｋ_ｙｓおよびｋ_ｘｓを表す。したがって、ステップ５１４の後、ｋ空間データセットは、エコーから入手されたデータ、および少なくとも部分的に、上記に説明されるようなエコー時間に基づいて推定されるデータを含む。ステップ５１４の後、ｋ空間データは、（ステップ５１６において）逆フーリエ変換され、治療画像を再構築することができる。ステップ５０４－５１６は、欠落しているｋ空間データが適切に推定されていることを示す、終了条件（例えば、２つの連続する反復間の画像差が、所定の閾値を下回る、またはある回数の反復が、実施されている）が満たされるまで、反復的に実施されてもよい。したがって、上記に説明されるような繰り返しの代替様式で制約を課すことによって、欠落しているｋ空間データの適切な推定が、達成されることができる。いくつかの実施形態では、反復的推定プロシージャは、欠落しているｋ空間データを決定するように訓練される、畳み込みまたは再帰型ニューラルネットワークを使用してもよい（例えば、ｈｔｔｐｓ：／／ａｒｘｉｖ．ｏｒｇ／ｐｄｆ／１７１２．０２８６２．ｐｄｆ参照（参照することによって本明細書に組み込まれる））。図５Ｂは、推定アプローチ５００を使用した、反復間での画質の改良を描写する。

【0051】

再び図５Ａを参照すると、最初に入手されたアンダーサンプリングされたｋ空間データによって更新されたｋ空間データセットを置換するために、第３の制約が課された後、随意のステップ５１８が、実施され、好適なノイズフィルタリングアプローチを使用して、再構築された治療画像のＳＮＲを改良してもよい。例えば、画像データセットＸ内の各場所（ｙ，ｘ）が、空間的に非常に相関される（これは、全てのエコーが同一の解剖学的構造に対応することを意味する）、Ｎ個（すなわち、エコーの数）の点を有するため、種々の実施形態では、ノイズフィルタリングアプローチは、空間的相関を活用し、治療画像のノイズレベルを低減させる、局所的に低ランクの（ＬＬＲ）正則化を実装する。画像のＳＮＲを改良するためのＬＬＲ正則化を使用するためのアプローチが、例えば、Ｔ． Ｚｈａｎｇ ｅｔ ａｌ．による「Ａｃｃｅｌｅｒａｔｉｎｇ Ｐａｒａｍｅｔｅｒ Ｍａｐｐｉｎｇ ｗｉｔｈ ａ Ｌｏｃａｌｌｙ Ｌｏｗ Ｒａｎｋ Ｃｏｎｓｔｒａｉｎｔ」Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ ｉｎ Ｍｅｄｉｃｉｎｅ， ｖｏｌｕｍｅ ７３， ６５５－６６１ （２０１５）、およびＪ．Ｆ． Ｃａｉ ｅｔ ａｌによる「Ａ Ｓｉｎｇｕｌａｒ Ｖａｌｕｅ Ｔｈｒｅｓｈｏｌｄｉｎｇ Ａｌｇｏｒｉｔｈｍ ｆｏｒ Ｍａｔｒｉｘ Ｃｏｍｐｌｅｔｉｏｎ」ＳＩＡＭ Ｊｏｕｒｎａｌ ｏｎ Ｏｐｔｉｍｉｚａｔｉｏｎ， ｖｏｌｕｍｅ ２０（４）， １９５６－１９８２ （２０１０）（その開示全体は、参照することによって本明細書に組み込まれる）等の学術刊行物に提供される。一実施形態では、ステップ５１８が、最後の５～１０回の反復においてのみ実施される。典型的には、算出された温度のＳＮＲが、ＬＬＲ正則化を使用して、２倍改良され得る。加えて、または代替として、ウェーブレット変換が、ノイズフィルタリングのために使用されることができる。

【0052】

治療画像の中の各ピクセルの大きさを、対応するベースライン画像の中のものと同一であるように設定する、第２の制約が、いくつかの実施形態では、緩和され得ることに留意されたい。例えば、本制約は、最初の数回の反復（例えば、１～５回の反復）においてのみ課され得、推定されるデータの合理的収束が達成された後、本制約は、除去されてもよい（すなわち、ステップ５１０は、最初の数回の反復の後、スキップされてもよい）。したがって、治療画像の中のピクセルの大きさは、対応するベースライン基準画像の中のものに近接し得るが、必ずしもそれと同じであるわけではない。

【0053】

続いて、治療画像は、対応するベースライン画像に対して比較され、着目領域の温度マップを算出的に発生させ、それによって、超音波手技の間の温度変化のリアルタイムの監視を提供してもよい（ステップ３１０）。一実施形態では、温度マップは、加重関数Ｗ_ｉを使用して、エコートレイン内の各エコーＥ_ｉと関連付けられる、算出された温度変化ΔＴ_ｉに加重することによって発生される。上記に説明されるように、温度変化ΔＴ_ｉは、（エコー時間ＴＥ_ｉにおいて）エコーＥ_ｉから入手された画像とその対応するベースライン画像との間の位相差ΔΦ_ｉに基づいて算出され得る。

【化22】

各エコーと関連付けられる温度変化ΔＴ_ｉの加重された組み合わせが、次いで、治療画像と関連付けられる温度変化を決定するために算出され得る。

【化23】

ステップ３０６－３１０が、手技の間に着目領域を監視するために、超音波手技の全体を通して反復的に実施されてもよい。

【0054】

故に、本発明の種々の実施形態は、（例えば、マルチエコーＧＲＥパルスシーケンスの中の１つ以上のエコー内のｋ空間データをアンダーサンプリングすることによって）ＭＲ撮像率を効率的に向上させながら、（例えば、上記に説明される制約および反復的推定アプローチを使用して、欠落している（または入手されていない）ｋ空間データを適切に決定することによって）超音波手技の間にリアルタイムで領域を監視するために十分な分解能を伴う画像を提供し得る。

【0055】

上記に与えられる実施例では、単一のスライスのみが、所与のＴＲ内で励起される。故に、Ｒの加速度を伴う、３秒の総走査時間が、所望されるとき、ＴＲの数は、Ｒ倍低減され、結果として、ＴＲは、Ｒだけ増大する。例えば、上記に示されるように、ＴＲは、Ｒ＝６であるとき、２７ミリ秒から１４２ミリ秒まで増大し、加えて、ＳＮＲは、約１．９増大する。種々の実施形態による代替アプローチが、各ＴＲの間、全てのスライスを同時に励起する。ＲＦ励起の間、各スライスは、異なる位相によって励起され、ＲＦ波形は、各スライスを別個に励起するために要求される、波形の線形結合である。Ｒ個のスライスが、同時に励起される場合、それらは、相互に重複する。それらを分離するために、Ｒ個の波形が、印加され、Ｒ個のデータセットを入手してもよく、各スライスの位相は、公知の方法において波形間で線形に変化する。最後に、スライスは、Ｒ個のデータセットに対して逆フーリエ変換またはアダマール変換を適用することによって分離され得る。本場合では、２７ミリ秒の短いＴＲが、Ｒ個の波形の全てを励起するために使用されてもよい。しかしながら、全てのスライスが、同時に励起されるため、ＳＮＲは、ＴＲ＝２７ミリ秒を伴う単一のスライスと比較して、

【化24】

だけ増大する。これは、大きいＲに関して非常に有益である。例えば、Ｒ＝８に関して、ＳＮＲ利得は、

【化25】

である。図６は、その間に周波数差Ω_０を有する、２つのＲＦパルスＲＦ_１およびＲＦ_２に適用し、その間にある位相差（例えば、π）が存在するように位相エンコードされる、２つのスライスを同時に入手するための、例示的アプローチを描写する。２つのＭＲ信号Ｓ_ｎ１およびＳ_ｎ２を受信するステップに応じて、２つのスライスが、以下の方程式を使用することによって分離され得る。
スライス_１＝（Ｓ_ｎ１＋Ｓ_ｎ２）／２
スライス_２＝（Ｓ_ｎ１－Ｓ_ｎ２）／２

【0056】

ＭＲ撮像率を加速させるために上記に説明されるアプローチを採用することは、有利なこととして、画像の空間分解を向上させながら、従来のＭＲ走査アプローチと比較して温度測定毎により少ない数のスライスを入手し得る。例えば、Ｒ＝６であるとき、それぞれが、２５６個の位相エンコーディングステップ（すなわち、Ｎ_１＝２５６）を有する、３つのＭＲ撮像スライスを入手することが、可能であり得る。より多くの位相エンコーディングステップを有する画像が、次いで、１２８個の位相エンコーディングステップが使用される、従来のＭＲ画像と比較して、より高い空間分解能を提供し得る。加えて、所望されるとき、上記に説明されるアプローチは、入手された画像の時間分解能を向上（すなわち、温度測定毎に要求される時間を短縮）させ得る。例えば、Ｒ＝６であるとき、３つのスライスを有する撮像データセットを１．６秒で入手することが、可能であり得る。本入手率は、１つのスライスを有する撮像データセットを入手するために３秒が要求される、従来のアプローチと比較して、より迅速である。

【0057】

方程式（１４ａ－１４ｃ）に記載されるアプローチは、拡張されることができる。一実施形態では、ベースライン画像が、場マップの変動または水の移動に起因して安定していない場合、いくつかのベースライン画像が、収集され、治療相が、ベースライン画像の線形結合と適合される。いくつかの実施形態では、本アプローチは、最小化されるべきｃｏｎｓｔ関数として公式化されることができる。

【化26】

ここで、

【化27】

は、コイルｎおよびエコーｊに関してアンダーサンプリングされたデータであり、Ｅ_ｊは、エコーｊに関するアンダーサンプリングされたエンコードパターンであり、Ｔは、実空間からｋ空間への変換演算子であり、ｘは、治療画像に対するベースラインピクセルに関してモデル化された、実空間内のピクセルであり、Ｒｅｇ（ｘ）は、任意の以前の知識（Ｔｉｋｈｏｎｏｖ、全変動、疎性等）に基づく正規化項である。

【0058】

そのようなモデルの一実施例は、ベースライン画像が、（例えば、水移動に起因して）安定しておらず、Ｔ２＊のある弛緩が、加熱に応じて生じる（Ｒ２＊＝１／Ｔ２＊である）場合である。本場合では、ｘは、以下のように算出されることができる。

【化28】

ρは、エコーｊにおけるｌベースラインの大きさであり、Ｎ_ｂは、ベースラインの数である。正則化項が、測定された空間場所（２Ｄまたは３Ｄのいずれか）を横断して、Ｆ_１および／またはＦ_０および／またはＲ_２ ^＊および／または振幅Ａに適用されることができる。ＬＬＲ等の他の正規化項も、エコートレインに沿って実装されることができる。加えて、または代替として、ＬＬＲは、治療データとベースラインとの間の差異に適用されることができる。加えて、または代替として、時空間的相関が、時間的制約付き再構築（ＴＣＲ）を使用して考慮されることができる。

【0059】

欠落しているｋ空間データに適応することに対する別のアプローチは、全結合層を含み得る、畳み込みニューラルネットワーク（ＣＮＮ）または他のニューラルネットワークアーキテクチャ、再帰型ニューラルネットワーク（ＲＮＮ）、「ニューラルインニューラル」（ＮｉＮ）ネットワーク、および他のもの等の訓練されたニューラルネットワークを使用することである。周知であるように、ニューラルネットワークは、ヒトの脳に類似する様式において情報を処理する。ニューラルネットワークは、具体的な問題を解決するように並行して稼働する、多数の非常に相互接続される処理要素（ニューロン）から構成される。ニューラルネットワークは、実施例によって学習し、具体的なタスクを実行するためにプログラムされることはできない。実施例は、慎重に選択されなければならず、そうでなければ、有用な時間が、浪費されるか、さらに悪いことには、ネットワークが、不正確に機能する場合がある。「教師あり学習」の原理に基づいて動作するニューラルネットワークは、標識された入力画像（標識は、例えば、分類等の出力に関連のある特性を規定する）を使用して訓練される。

【0060】

いくつかの実施形態では、ニューラルネットワークは、ｋ空間内のエイリアス処理された画像からｋ空間内のアーチファクトのない画像へのマッピングを学習する。ネットワーク訓練のための１つの画像駆動型の方法では、損失関数Ｌは、正則化項、例えば、

【化29】

を含む。本場合では、ｎは、１または２のいずれかであり、正則化パラメータＲｅ（ｘ，θ）は、確率分布関数についての関連のある以前の知識から導出されることができる。一実施形態では、Ｔｉｋｈｏｎｏｖ正則化が、データの平滑性に関して使用される。

【0061】

他の教師あり学習アプローチでは、エンコード演算子

【化30】

の知識が、ｋ空間内の所与のデータから実空間内の

【化31】

を見出すために使用される。いくつかの実施形態では、アンダーサンプリングされたｋデータのデータ一貫性が、損失関数のための付加的な項として考慮される。

【化32】

式中、

【化33】

は、実空間からｋ空間への変換である。ＤＬ温度測定では、遡及的ベースライン（ｂｌ）データｋ^ｂｌのデータ一貫性を考慮する、付加的な損失関数が、採用され得る。

【化34】

【0062】

ｋ個の軌道を表す、エンコーダ演算子Ｅ_ｊはまた、ニューラルネットワークアーキテクチャ（ＮＮＡ）内に自由パラメータを含むこともでき、学習されることができる。いくつかの実施形態では、異なるエコーのｋ空間データ軌道は、異なり得る。一実施形態では、物理的制約が、損失関数として勾配軌道に追加され、ｋ個のサンプリング軌道が学習されることを可能にする。Ｅ_ｊによってサンプリングされるパターンが、渦巻状、半径方向、または他の幾何学的なパターン軌道から学習されながら、軌道を特徴付けるパラメータが、ニューラルネットワークによって学習されることを可能にすることができる。

【0063】

加えて、または代替として、実際のｋ個の軌道は、渦電流効果および／または付随する勾配効果のため、および／またはＢ_０不均質性のためのパラメータを含む物理モデルに基づいて計算されることができる。ニューラルネットワークは、入力として、（Ｂ_０場マップまたは事前走査勾配遅延測定値等の）事前走査データを受信することができる。完全にサンプリングされたベースライン軌道とともに、ｋ個の軌道のための自動補正が、達成されることができる。一実施形態では、ｋ個の軌道は、治療走査に先立って勾配波形を変動させることによって、補正される。加えて、または代替として、ｋ個の軌道を補正するために使用されるニューラルネットワークは、ＮＮＡ全体内に統合される。

【0064】

加えて、または代替として、別の物理的に駆動される損失関数、すなわち、エコー画像から構築され得る場マップパラメータが、使用されることができる。一実施形態では、方程式（１６）に基づいて計算される、ＴＭＡＰ（ΔＴ）が、ＮＮＡから再構築される、

【化35】

に対して比較される、グランドトゥルースとして使用される。ＴＭＡＰは、各ピクセルにおいて計算されたΔＴを使用して取得される、実空間内の熱マップである。

【化36】

【0065】

一実施形態では、ＵＮｅｔ ＣＮＮアーキテクチャが、採用される。Ｎ_{ｃｏｉｌｓ}によってアンダーサンプリングされたｋデータが、Ｎ_{ｃｏｉｌｓ}ベースラインｋデータとともに、入力データとして使用されることができる。一実施形態では、ゼロ充填が、欠落しているｋ個のデータ点において設定される。非直交ｋサンプリングの場合では、グリッディングが、採用されることができる。次いで、ｋデータから実空間への逆フーリエ変換が、取得される。一実施形態では、ＣＮＮは、損失関数Ｌ_ｉｍａｇが、ＣＮＮステップ内でアクティブであるが、ＣＮＮステップの間において、データの一貫性損失項Ｌ_ｕｓ、Ｌ_ｂｌが、使用される場合、反復的に使用される。別の実施形態では、エンドツーエンドＣＮＮが、直接、ｋ空間からＴＭＡＰに進行することができる。訓練データは、多種多様な現実的特徴、すなわち、ファントム形状、磁場ドリフト、ＲＦ干渉、渦電流、異なる温度パターン、形状、サイズ、スケール、ベースラインと治療との間の移動、および他の特徴を考慮する、現実的かつ詳細なブロッホシミュレーションを使用して発生され得る。

【0066】

別の実施形態は、敵対的生成ネットワーク（ＧＡＮ）を利用する。一実施形態では、生成ネットワーク損失関数は、Ｌ１およびまたはＬ２正則化を伴う、Ｌ_{ｉｍａｇｅ}である。敵対的ネットワークは、データ一貫性損失関数を使用してもよい。さらに別の実施形態は、（ベースライン相を含む）治療相の順次性質を前提として、再帰型ニューラルネットワークを使用する。

【0067】

Ｃａｆｆｅ、Ｋｅｒａｓ、ＰｙＴｏｒｃｈ、Ｔｈｅａｎｏ、およびＴｅｎｓｏｒＦｌｏｗは、好適なニューラルネットワークプラットフォームである（設計選好に従って、実装されたシステムに対してクラウドベースまたはローカルであってもよい）。画像を伴うタスクに関して、ＣＮＮは、典型的には、ＮＮＡ内の主要成分として好ましい。本明細書において提示されるＣＮＮの使用法は、代表的であり、実際の実装は、代替ネットワークおよびＮＮＡとの組み合わせも使用することもできる。ＣＮＮは、画像を入力として処理し、設計に応じて、分類または別の画像を出力してもよい。例えば、Ｕ－ｎｅｔアーキテクチャに基づくＣＮＮネットワークが、入力データを画像に変えるために使用されることができる。ＣＮＮは、典型的には、特に、入力画像に対して動作し、出力を生産する、畳み込み層と、プーリング層とを含む。

【0068】

より一般的には、本発明の実施形態は、ＮＮＡを使用し、ベースラインｋデータまたは複数のベースラインとして収集された、完全データの追加の有無にかかわらず、欠落データを伴う入力ｋ空間画像から熱マップを発生させる、またはそれを識別する。１つのアプローチでは、訓練ライブラリは、ベースラインの欠落データと、完全データとを伴う一連のｋ空間画像を含み、標識は、対応する熱マップである。例えば、標識は、上記に説明されるように、訓練入力としての役割を果たす、欠落データを伴うｋ空間画像から発生されてもよい。

【0069】

本アプローチは、サブサンプルｋ空間から所望の熱マップまでの直接変換を利用する。別の実施形態では、熱マップは、間接的に２つのステップで生成され、すなわち、サブサンプリングされるｋ空間が、最初に、ＮＮＡを使用して複雑な画像（例えば、少なくとも２つのエコーからの画像）に変換され、結果として生じる画像が、次いで、従来のＰＲＦ技法を使用して熱マップを発生させるために使用される。サブサンプリングされたｋ空間からＭＲＩ画像を再構築するためにＮＮＡを活用する第１のステップが、異なる技法を使用して遂行されてもよい。１つのアプローチでは、ｋ空間は、（例えば、ゼロを欠落データに割り当てることによって）最初に完了され、除去するためにＣＮＮが使用されるアーチファクトを呈するであろう、画像に再構築される。別のアプローチでは、ＮＮＡは、サブサンプリングされたｋ空間から直接画像を発生させる。加えて、または代替として、ＮＮＡは、入力として、少なくとも１つの準最適な熱マップ（例えば、アーチファクトを伴うマップ）を受入し、出力として、１つ以上の改良された熱マップ（ＭＬ＿Ｍａｐ）を発生させるように訓練されることができる。いくつかの実施形態では、改良された熱マップは、（例えば、方程式１１を

【化37】

と置換することによって）再構築された画像をさらに補正するために使用されることができる。

【0070】

本明細書に説明されるＭＲ加速アプローチが、任意のマルチエコーＭＲパルスシーケンスに適用され得ること、およびアンダーサンプリングされたｋ空間データが、アンダーサンプリングが、マルチエコートレインのエコー内で行われる限り、任意の「増分された」寸法（例えば、位相エンコード寸法、渦巻状交互配置寸法、半径方向角度寸法等）に沿って収集され得ることに留意されたい。例えば、本明細書に説明されるＭＲ加速アプローチは、２Ｄ ＭＲ撮像スライスを入手するために、マルチスライスのマルチエコーＧＲＥパルスシーケンスおよび／またはマルチ帯域のマルチエコーＧＲＥパルスシーケンス内で適用されてもよく、ｋ空間軌道は、直交性、渦巻状、半径方向、または恣意的軌道であってもよい。例えば、図７Ａ－７Ｃを参照すると、ｋ空間データは、各ＴＲ内の中に３つのエコー７０２－７０６を伴う渦巻状様式において入手され得る。加えて、各エコーは、３つのショットを含むが、それらのうちの１つのみ（例えば、それぞれ、エコー７０２、７０４、７０６内のショット７１２、７１４、７１６）が、入手され、他の２つのショットが、欠落している。上記に説明されるような欠落データの推定も、同様に、本走査アプローチに適用され得る。例えば、欠落ショットが、反復的推定アプローチを使用して、上記に説明される制約に基づいて、数値的に算出されてもよい。同様に、ｋ空間データが、半径方向軌道を含む場合、アンダーサンプリングが、各エコー内の「アーム」の数を低減させるために実施されてもよい。続いて、制約とともに、反復的推定アプローチが、欠落アームを決定するために実装されてもよい。したがって、Ｍ個の走査ショットおよびＮ個のエコー（それによって、Ｍ×Ｎ個の総ショット）が、従来のアプローチにおいて、撮像スライスを入手するために要求されるが、本明細書に説明される加速されたＭＲ撮像アプローチを使用することは、有利なこととして、係数、すなわち、Ｒだけショットの総数（すなわち、Ｍ×Ｎ／Ｒ）を低減させ得る。同一の構想が、同様に、走査ショットの数を低減させるために、マルチショットのマルチエコーのエコー平面撮像（ＥＰＩ）においても使用され得る。

【0071】

いくつかの実施形態では、体積３Ｄ ＭＲ熱マップが、望ましい。本状況では、ｋ空間データの３つの寸法が、入手され得、寸法のうちの２つ（例えば、間接的Ｚエンコード、渦巻状交互配置等）が、増分され得る。故に、いずれか一方または両方の２つの増分された寸法に沿ったデータが、上記に説明される加速されたＭＲ撮像アプローチを使用して、アンダーサンプリングされ、推定されてもよい。例えば、治療に先立って、入手されたｋ空間データセットは、１つ以上の完全な３Ｄベースライン画像を含んでもよい。治療の間、入手されたｋ空間データセットは、各ＴＲ内の増分された寸法に沿った、異なる値に対応する、複数のショットを含んでもよく、加えて、異なるショットが、異なるＴＥにおいてサンプリングされてもよい。結果として、治療画像内の各ショットは、ある欠落データを有してもよい。再び、ベースラインと治療画像との間の主要な変化は、位相内にあるため、治療画像内の欠落データは、制約に基づいて上記に説明される、反復的アプローチを使用して推定されてもよい。加えて、３Ｄデータ入手は、直交入手、渦巻状入手、半径方向入手（例えば、ｈｔｔｐｓ：／／ｏｎｌｉｎｅｌｉｂｒａｒｙ．ｗｉｌｅｙ．ｃｏｍ／ｄｏｉ／ａｂｓ／１０．１００２／ｍｒｍ．２６８６２参照）、２Ｄの渦巻状平面をとの組み合わせにおける、半径方向入手（例えば、ｈｔｔｐｓ：／／ｗｗｗ．ｎｃｂｉ．ｎｌｍ．ｎｉｈ．ｇｏｖ／ｐｕｂｍｅｄ／２８６４３３８３およびｈｔｔｐｓ：／／ｏｎｌｉｎｅｌｉｂｒａｒｙ．ｗｉｌｅｙ．ｃｏｍ／ｄｏｉ／ａｂｓ／１０．１００２／ｍｒｍ．２６８６２参照）、または恣意的軌道であってもよい。一実施形態では、恣意的軌道は、入手されたベースラインに基づいて定義される。

【0072】

さらに、本明細書に説明されるマルチエコーＧＲＥパルスシーケンスを加速させるための方法は、ＭＲ温度測定のみに限定されない。マルチエコーＧＲＥパルスシーケンスを使用する他の用途もまた、加速されたＭＲデータ入手アプローチから恩恵を享受し得る。一実施例は、高速場マップが、場の変化に応じて、またはマルチエコーＧＲＥパルスシーケンスを使用した磁化率強調撮像（ＳＷＩ）の構築のために入手される、場マップの測定である。別の例示的用途は、位相と対照的に、エコー画像（例えば、ｆ－ＭＲＩにおける微小出血検出、ＢＯＬＤＰＳＤ）間のピクセルあたりの大きさのコントラストおよび信号減衰がアンダーサンプリングされ、推定される、Ｔ_２ ^＊の検査である。再び、これらの用途の全てにおいて、１つ以上の完全にサンプリングされたベースライン画像が、手技に先立って要求され得る。ベースライン画像が、治療画像に対して静的ではない場合、物理モデルが、変化を予測し、それによって、それを補償するための、最適化手技に統合されてもよい（例えば、方程式（１８）参照）。

【0073】

いくつかの状況では、入手された信号のＳＮＲが、上記に説明される制約を使用して欠落データの意味のある推定を提供するためには低すぎる場合がある。例えば、脳組織の微小出血を検出するとき、局所的な暗点が、エコートレイン画像に沿って出現し得（高速Ｔ_２＊）、結果として、これらの点におけるＴＥと位相差との間の線形関係が、低ＳＮＲに起因して測定されることができない。一実施形態では、加重マスクが、そこから結果として生じる信頼性を欠く制約が、再構築された画像に影響を及ぼさないように、暗点に対応するピクセルに適用されることができる。

【0074】

加えて、ＭＲＩシステムは、複数のＭＲＩコイルを含んでもよく、各コイルは、標的領域からのＭＲ信号の少なくとも一部を受信し得る。超音波手技に先立って、および／またはその間に入手されるｋ空間画像はそれぞれ、ＭＲＩコイルのうちの１つ以上のものと関連付けられてもよい。結果として、再構築された画像は、上記に説明されるような関連付けられたエコー時間および関連付けられたＭＲＩコイルに基づき得る。一実施形態では、複数のコイルによって受信されたＭＲ信号から抽出される情報が、融合されてもよく、融合された情報は、単一のＭＲコイルを使用して測定される情報と均等であり得る（または少なくともそれに類似し得る）。加えて、または代替として、各コイルによって受信されたＭＲ信号は、上記に説明されるように処理され、画像を再構築する、および／または熱マップを発生させてもよい。複数のコイルと関連付けられる画像および／または熱マップは、次いで、例えば、加重平均を使用して収集され、標的領域の平均化された画像および／または熱マップを生成してもよい。

【0075】

一般に、ｋ空間データをアンダーサンプリングするステップと、ＭＲＩシステムおよび／または超音波システムのコントローラ内に統合されているか、または別個の外部コントローラによって提供されるかどうかにかかわらず、上記に説明されるように、反復的推定アプローチを使用して欠落しているｋ空間データを算出的に決定するステップとを含む、加速されたＭＲ撮像入手を実施するための機能性は、ハードウェア、ソフトウェア、または両方の組み合わせ内に実装される、１つ以上のモジュールに構造化されてもよい。機能が、１つ以上のソフトウェアプログラムとして提供される、実施形態に関して、プログラムは、ＰＹＴＨＯＮ、ＦＯＲＴＲＡＮ、ＰＡＳＣＡＬ、ＪＡＶＡ（登録商標）、Ｃ、Ｃ＋＋、Ｃ＃、ＢＡＳＩＣ、種々のスクリプト言語、および／またはＨＴＭＬ等のいくつかの高水準言語のいずれかで記述されてもよい。加えて、ソフトウェアは、標的コンピュータ（例えば、コントローラ）上に常駐するマイクロプロセッサを対象とするアセンブリ言語で実装されることができ、例えば、ソフトウェアは、これが、ＩＢＭ ＰＣまたはＰＣクローン上で起動するように構成される場合、Ｉｎｔｅｌ ８０ｘ８６アセンブリ言語で実装されてもよい。ソフトウェアは、限定ではないが、フロッピー（登録商標）ディスク、ジャンプドライブ、ハードディスク、光ディスク、磁気テープ、ＰＲＯＭ、ＥＰＲＯＭ、ＥＥＰＲＯＭ、フィールドプログラマブルゲートアレイ、またはＣＤ－ＲＯＭを含む、製造品上で具現化されてもよい。ハードウェア回路網を使用する実施形態は、例えば、１つ以上のＦＰＧＡ、ＣＰＬＤ、またはＡＳＩＣプロセッサを使用して実装されてもよい。

【0076】

加えて、本明細書において広く使用される、用語「コントローラ」は、上記に説明されるような任意の機能性を実施するために利用される、必要なハードウェアコンポーネントおよび／またはソフトウェアモジュールの全てを含み、コントローラは、複数のハードウェアコンポーネントおよび／またはソフトウェアモジュールを含んでもよく、機能性は、異なるコンポーネントおよび／またはモジュールの間に分散されることができる。最後に、（上記に説明されるように）ベースラインの大きさの摂動が存在し、物理モデルが、想定されるとき、上記に説明される最適化アプローチが、同様に、ベースライン画像の調節も含み得る。

【0077】

本発明のある実施形態が、上記に説明される。しかしながら、本発明が、それらの実施形態に限定されず、むしろ、本明細書に明示的に説明されるものへの追加および修正もまた、本発明の範囲内に含まれることに明確に留意されたい。

【図1】

【図2A】

【図2B】

【図2C】

【図3】

【図4A】

【図4B】

【図4C】

【図4D】

【図5A】

【図5B】

【図6】

【図7】