[論文レビュー] Stability Analysis for Magnetic Resonance Elastography

Habib Ammari, Alden Waters+1 more|arXiv (Cornell University)|Sep 17, 2014

Ultrasound Imaging and Elastography参考文献 22被引用数 3

ひとこと要約

本稿では、磁気共鳴エラストグラフィ（MRE）における内部変位測定値からせん断弾性率分布を再構成するためのリプシッツ安定性推定を、準不可撓性組織をモデル化するストークス系を用いて確立する。高圧縮率モジュラスの極限において弾性方程式をストークス系に還元することで、著者らは逆問題が適切に定義されていることを証明し、2次元および3次元におけるせん断弾性率再構成のランドヴェーバー反復スキームの収束を示す。3次元では安定性を保証するため、追加の測定値が必要となる。

ABSTRACT

We consider the inverse problem of finding unknown elastic parameters from internal measurements of displacement fields for tissues. The measurements are made on the entirety of a smooth domain. Since tissues can be modeled as quasi-incompressible fluids, we examine the Stokes system and consider only the recovery of shear modulus distributions. Our main result is to establish Lipschitz stable estimates on the shear modulus distributions from internal measurements of displacement fields. These estimates imply convergence of a numerical scheme known as the Landweber iteration scheme for reconstructing the shear modulus distributions.

研究の動機と目的

• 磁気共鳴エラストグラフィ（MRE）における内部変位測定値からせん断弾性率分布を回復する逆問題に取り組む。
• 軟組織における測定された変位場からせん断弾性率を再構成するための数学的安定性推定を確立する。
• 適切な正則性および安定性条件の下で、せん断弾性率再構成のためのランドヴェーバー反復スキームの収束を証明する。
• 3次元への分析の拡張を行い、安定性を保証するための追加の内部測定値の必要性を検討する。
• 過決定系の偏微分方程式系を用いたMREにおける数値再構成手法の理論的基盤を提供する。

提案手法

• 無限大への圧縮率モジュラス（λ → ∞）の極限において、時間調和弾性方程式系をストークス系に還元し、準不可撓性組織をモデル化する。
• せん断弾性率μの過決定系の偏微分方程式系として逆問題を定式化する。
• 得られた系が楕円型であり、ロパティンスキー境界条件を満たすことを証明し、適切な定義を保証する。
• 楕円型系の安定性推定を用いて、内部変位データからのせん断弾性率のリプシッツ境界を導出する。
• 計算された変位場と測定された変位場の差のL²ノルムに基づく不一致汎関数J[μ]を導入する。
• 随伴系を用いてフレシェ微分を計算し、勾配降下法を用いてランドヴェーバー反復スキームを適用する。

実験結果

リサーチクエスチョン

• RQ1 MREにおける内部変位測定値からせん断弾性率を再構成する際、リプシッツ安定性推定を確立できるか？
• RQ2 準不可撓性組織において、ストークス系近似のもとでのせん断弾性率の逆問題はどのように振る舞うか？
• RQ3 2次元および3次元におけるランドヴェーバー反復スキームの収束を保証する条件は何か？
• RQ4 なぜ3次元MREでは、安定な再構成を保証するために追加の内部測定値が必要なのか？
• RQ5 多周波数測定値やその他の改良手法を用いて、過決定系の核を除去できるか？

主な発見

• 2次元および3次元の両方において、内部変位データからのせん断弾性率分布の再構成に対してリプシッツ安定性推定が確立された。
• 3次元では、安定性と収束を保証するために、異なる境界条件の下での第2の変位測定値の追加が必要である。
• 初期推定値が十分に近い場合、ランドヴェーバー反復スキームはH⁴(Ω)ノルムにおいて真のせん断弾性率分布μ_trに収束する。
• 収束結果は、過決定系(L, B)の核Kが自明であるという仮定の下で成り立つ。
• 理論的枠組みは、正則性および安定性条件の下で、勾配に基づく数値的再構成スキームの使用を支持する。
• 本分析は、特に高解像度のせん断弾性率画像による早期がん検出を目的としたMREにおける数値的再構成の厳密な数学的基盤を提供する。