[论文解读] The inverse problem for the local geodesic ray transform

Günther Uhlmann, András Vasy|arXiv (Cornell University)|Oct 7, 2012

Numerical methods in inverse problems参考文献 8被引用 30

一句话总结

该论文在黎曼流形（维度≥3）上，于边界附近满足局部严格凸性条件时，建立了测地线X射线变换的局部可逆性与稳定性。通过微局部分析与Neumann级数重构方法，证明了在指数加权范数下的注入性与稳定性估计，并在由严格凸超曲面全局叶状分解的条件下，将结果推广至全局注入性。

ABSTRACT

Under a convexity assumption on the boundary we solve a local inverse problem, namely we show that the geodesic X-ray transform can be inverted locally in a stable manner; one even has a reconstruction formula. We also show that under an assumption on the existence of a global foliation by strictly convex hypersurfaces the geodesic X-ray transform is globally injective. In addition we prove stability estimates and propose a layer stripping type algorithm for reconstruction.

研究动机与目标

• 解决黎曼流形中维度≥3的严格凸区域上测地线X射线变换的局部反问题。
• 在边界点附近满足局部凸性条件下，建立变换在局部邻域内的稳定注入性。
• 通过Neumann级数提供逆变换的重构公式。
• 在由严格凸超曲面构成的全局叶状分解条件下，将局部结果推广至全局注入性。
• 在小度量扰动及区域紧子集条件下，推导出一致稳定性估计。

提出的方法

• 使用定义函数ρ，并构造光滑函数˜x，使得d˜x(p) = −dρ(p)，以定义局部邻域O_p = {˜x > −c} ∩ X。
• 在散射微局部分析框架下应用加权Sobolev空间H^s_γ，其中权重为e^{γ/(˜x + c)}。
• 通过紧支集光滑函数（接近高斯函数）的傅里叶变换，证明变换算子A = x^{-1}e^{-γ/x}Ae^{γ/x}的边界主符号的椭圆性。
• 通过Neumann级数构造算子A的拟逆，证明其在有限维核的正交补空间上的有界性与可逆性。
• 利用分解A = L ∘ I，其中I为X射线变换，L为拟微分算子，将A的稳定性与I的稳定性关联起来。
• 通过加权Sobolev空间与标准Sobolev空间之间的包含映射，推导出在H^s范数下带指数权重的最终稳定性估计。

实验结果

研究问题

• RQ1 在边界点附近满足局部凸性假设时，测地线X射线变换能否以稳定方式实现局部可逆？
• RQ2 在黎曼度量的小扰动下，变换是否仍保持注入性与稳定性？
• RQ3 能否在由严格凸超曲面构成的全局叶状分解条件下，从局部可逆性推出全局注入性结果？
• RQ4 是否存在一种构造性方法（如Neumann级数）来求解反问题？
• RQ5 在反问题中，稳定性估计的最优正则性与权重结构为何？

主要发现

• 对每个p ∈ ∂X，存在邻域O_p = {˜x > −c} ∩ X，使得在H^s(O_p)（s ≥ 0）上，局部测地线X射线变换是注入的。
• 存在稳定性估计：||f||_{H^{s-1}_γ(O_p)} ≤ C||If|_{M_{O_p}}||_{H^s(M_{O_p})}，其中C对小c及小度量扰动一致。
• 该稳定性估计在指数加权Sobolev空间H^s_γ(O_p) = {f ∈ H^s_loc(O_p) : e^{-γ/(˜x + c)}f ∈ H^s(O_p)}中成立，对任意γ > 0均有效。
• 构造了Neumann级数重构公式，提供了一种稳定且显式的方法，用于从If恢复f。
• 在由严格凸超曲面构成的全局叶状分解条件下，若补集测度为零，则全局X射线变换在L^2(X)上是注入的；若补集为空集，则在H^s(X)（s > n/2）上也是注入的。
• 稳定性常数在c与ρ₀满足小量条件时一致，且结果可推广至小度量扰动。