[论文解读] Inverse diffusion problems with redundant internal information

François Monard, Guillaume Bal|arXiv (Cornell University)|Jun 21, 2011

Numerical methods in inverse problems参考文献 10被引用 5

一句话总结

本文提出了两种新颖的重构方法，用于从冗余内部测量 Hᵢⱼ(x) = σ²ᵃ(x)|∇uᵢ|² 中恢复标量扩散系数 σ(x)，其中 uᵢ 是椭圆方程的解。第一种方法求解一个过定的常微分方程组，将二维和三维结果推广至任意维度，并几乎涵盖所有 α ∈ ℝ；第二种方法求解一个线性椭圆方程组，证明在特定条件下具有唯一解，通过基于框架的边界条件和相容性约束实现稳定重构。

ABSTRACT

This paper concerns the reconstruction of a scalar diffusion coefficient $σ(x)$ from redundant functionals of the form $H_i(x)=σ^{2α}(x)| abla u_i|^2(x)$ where $α\in\Rm$ and $u_i$ is a solution of the elliptic problem $ abla\cdot σ abla u_i=0$ for $1\leq i\leq I$. The case $α=\frac12$ is used to model measurements obtained from modulating a domain of interest by ultrasound and finds applications in ultrasound modulated electrical impedance tomography (UMEIT) as well as ultrasound modulated optical tomography (UMOT). The case $α=1$ finds applications in Magnetic Resonance Electrical Impedance Tomography (MREIT). We present two explicit reconstruction procedures of $σ$ for appropriate choices of $I$ and of traces of $u_i$ at the boundary of a domain of interest. The first procedure involves the solution of an over-determined system of ordinary differential equations and generalizes to the multi-dimensional case and to (almost) arbitrary values of $α$ the results obtained in two and three dimensions in \cite{CFGK} and \cite{BBMT}, respectively, in the case $α=\frac12$. The second procedure consists of solving a system of linear elliptic equations, which we can prove admits a unique solution in specific situations.

研究动机与目标

• 将从冗余内部测量中稳定重构扩散系数 σ(x) 的方法推广至任意空间维度及几乎所有 α ∈ ℝ。
• 构建一个重构框架，利用内部泛函 Hᵢⱼ(x) = σ²ᵃ(x)∇uᵢ·∇uⱼ，其建模了如 UMEIT、UMOT、MREIT 和 ImpACT 等混合成像技术。
• 通过确保在定义域内每一点处 n 个向量 Sᵢ = σᵃ∇uᵢ 构成一个框架，建立 σ 重构的唯一性与稳定性。
• 分析并强制执行由冗余测量引发的过定系统中的相容性条件。
• 提出一种关于 uᵢ 的新椭圆系统，其系数由测量值 Hᵢⱼ 导出，与 σ 无关，并在 Fredholm 交替条件下证明其唯一可解性。

提出的方法

• 将原始反问题转化为关于 Sᵢ = σᵃ∇uᵢ 的方程组，从而将 σ 显式地从系统中消除。
• 利用 ∇uᵢ 为无旋场且椭圆方程 ∇·(σ∇uᵢ) = 0 的事实，推导出涉及 Hᵢⱼ = Sᵢ·Sⱼ 的非线性过定系统。
• 构造 uᵢ 的边界条件，使得在定义域内每一点处向量 Sᵢ 构成 ℝⁿ 中的框架，从而保证可逆性与稳定性。
• 在第一种方法中求解过定常微分方程组，将先前在二维和三维中的结果推广至更高维度。
• 构建一个关于 uᵢ 的线性椭圆方程组，其系数仅依赖于测量值 Hᵢⱼ，而不依赖于 σ，并通过 Fredholm 交替条件证明其唯一可解性。
• 通过微分几何约束强制执行相容性条件，如 dF^♭ = 0（其中 F = ∇logσ），从而导出 ∑ₚ,ₚ Mₚₚⁱʲ : ∇uₚ ⊗ ∇uₚ = 0。

实验结果

研究问题

• RQ1 能否将二维和三维中基于冗余内部测量的 σ(x) 稳定重构推广至任意空间维度 n ≥ 2 及几乎所有 α ∈ ℝ？
• RQ2 边界数据需满足何种条件，才能确保在每一点处向量 Sᵢ = σᵃ∇uᵢ 构成 ℝⁿ 中的框架，从而实现 σ 的唯一重构？
• RQ3 如何推导并强制执行过定常微分方程组与椭圆方程组的相容性条件，以确保重构的一致性与路径无关性？
• RQ4 在何种条件下，所提出的关于 uᵢ 的线性椭圆方程组具有唯一解，且与 σ 无关？
• RQ5 能否仅使用测量泛函 Hᵢⱼ = σ²ᵃ|∇uᵢ|² 来重构 σ，而无需依赖已知不稳定的非冗余测量？

主要发现

• 第一种方法在任意维度 n ≥ 2 及几乎所有 α ∈ ℝ 下，提供了 σ 重构的全局唯一性与稳定性结果，前提是 Sᵢ 向量构成框架。
• 由框架条件与内部测量 Hᵢⱼ 导出的过定常微分方程组是可解的，并能实现 σ 的稳定重构。
• 第二种方法构建了一个关于 uᵢ 的线性椭圆方程组，其系数仅依赖于测量值 Hᵢⱼ，且对除离散 (n, α) 值集合外的所有情况证明了唯一可解性。
• 重构的相容性条件编码于 F = ∇logσ 为梯度的条件中，从而导出 ∑ₚ,ₚ Mₚₚⁱʲ : ∇uₚ ⊗ ∇uₚ = 0，该条件必须由数据满足。
• 复几何光学（CGO）解的使用确保了在 n ≥ 3 维下框架无关性，从而可构造合适的边界条件。
• 该框架推广了二维与三维中 α = 1/2 和 α = 1 时的先前结果，覆盖了 UMEIT、UMOT、ImpACT 和 MREIT 的应用。