[论文解读] Stability Analysis for Magnetic Resonance Elastography

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

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

一句话总结

本文通过将准不可压缩组织的模型建立在Stokes系统之上，从内部位移测量中建立了磁共振弹性成像（MRE）中剪切模量分布重建的Lipschitz稳定性估计。通过在压缩模量趋于无穷大的极限下将弹性方程简化为Stokes系统，作者证明了该反问题的适定性，并在二维和三维情况下展示了Landweber迭代格式在剪切模量重建中的收敛性，其中在三维情况下需要额外的测量以确保稳定性。

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）中内部位移测量数据重建剪切模量分布的反问题。
• 为从软组织中测量的位移场重建剪切模量提供数学稳定性估计。
• 在适当的正则性和稳定性条件下，证明用于重建剪切模量的Landweber迭代格式的收敛性。
• 将分析扩展至三维，其中为确保稳定性需要额外的内部测量。
• 为使用过定偏微分方程组的MRE数值重建方法提供理论基础。

提出的方法

• 在压缩模量趋于无穷大（λ → ∞）的极限下，将时谐弹性系统简化为Stokes系统，以建模准不可压缩组织。
• 将反问题表述为剪切模量μ的过定偏微分方程组。
• 证明所得系统为椭圆型，并满足Lopatinskii边界条件，从而确保适定性。
• 利用椭圆型系统的稳定性估计，从内部位移数据推导出剪切模量的Lipschitz界。
• 基于计算位移场与测量位移场之间L²范数差异，引入一个偏差泛函J[μ]。
• 通过梯度下降法在泛函上应用Landweber迭代格式，并利用伴随系统计算Fréchet导数。

实验结果

研究问题

• RQ1 能否为从MRE中内部位移测量数据重建剪切模量建立Lipschitz稳定性估计？
• RQ2 在准不可压缩组织中，Stokes系统近似下剪切模量的反问题行为如何？
• RQ3 在二维和三维中，何种条件能确保Landweber迭代格式的收敛性？
• RQ4 为何在三维MRE中需要额外的内部测量才能实现稳定重建？
• RQ5 能否通过多频测量或其他增强手段消除过定系统核（kernel）？

主要发现

• 在二维和三维情况下，均建立了从内部位移数据重建剪切模量分布的Lipschitz稳定性估计。
• 在三维情况下，需增加在不同边界条件下测量的第二个位移场，以确保稳定性与收敛性。
• 当初始猜测足够接近时，Landweber迭代格式在H⁴(Ω)范数下收敛于真实的剪切模量分布μ_tr。
• 收敛结果成立的条件是过定系统（L, B）的核K为平凡核。
• 该理论框架支持在满足正则性和稳定性条件下的基于梯度的数值重建方法在MRE中的应用，并保证收敛性。
• 该分析为MRE中的数值重建提供了严格的数学基础，尤其适用于通过高分辨率剪切模量成像实现早期癌症检测。