Regularized Subjective-Surface Flow with Monotone Reaction: Viscosity Well-Posedness and Convergence

技术摘要：具有单调反应项的正则化主观表面流

问题陈述
本研究针对用于显微图像中接触及分裂细胞核分割的正则化主观表面模型的数学分析进行了探讨。该模型定义在光滑有界区域 $\Omega \subset \mathbb{R}^d$ ($d \ge 2$) 上，并结合了依赖于图像的边缘系数和预设的相互作用权重。带有参数 $\varepsilon, \nu > 0$ 的正则化模型的控制方程为：
$$\partial_t u_{\varepsilon,\nu} = \nu \Delta u_{\varepsilon,\nu} + A_\varepsilon(\nabla u_{\varepsilon,\nu}) \text{div} \left( \frac{G(x) \nabla u_{\varepsilon,\nu}}{A_\varepsilon(\nabla u_{\varepsilon,\nu})} \right) - \mu R_\eta(x, u_{\varepsilon,\nu})$$
其中 $A_\varepsilon(p) = p / \sqrt{\varepsilon^2 + |p|^2}$，且 $R_\eta$ 是由阈值化图像数据导出的单调反应项。本研究关注当 $(\varepsilon, \nu) \to (0, 0)$ 时的消失正则化极限（vanishing-regularization limit）。

主要的数学挑战在于空间梯度 $\nabla u$ 消失的点。所代表的极限几何算子（即加权水平集平均曲率流）在零梯度处是未定义的。标准的粘性解理论要求通过上、下半连续扩张（几何方向极限）来定义该算子。然而，正则化算子在 $(\varepsilon, \nu) \to (0, 0)$ 时的直接联合极限通常并不等于这些几何极限。因此，标准的失连续算子稳定性理论无法唯一确定极限方程，这使得必须进行更精细的分析，以确保正则化模型收敛到预期的几何流。

方法论
本文采用了全非线性抛物型方程的粘性解理论。其方法论通过以下几个关键步骤展开：

1. 极限问题的构建： 极限方程被构建为一个 Dirichlet 粘性问题。几何算子 $K(x, p, X)$ 在 $p \neq 0$ 时有定义，其在 $p=0$ 处的行为通过源自水平集平均曲率流理论的特定下、上半连续扩张（$L_{geo,*}$ 和 $L^*_{geo}$）来表征。
2. 正则化极限分析： 作者显式地计算了正则化算子 $L_{\varepsilon,\nu}$ 在 $(\varepsilon, \nu, p) \to (0, 0, 0)$ 时的直接上极限和下极限。他们证明了当 Hessian 矩阵是定定时，这些直接极限（$J_*$ 和 $J^*$）与几何扩张不同，从而证明了标准的算子收敛性是不充分的。
3. 消失导数形式（Vanishing-Derivative Formulation）： 为了弥合正则化形式与几何形式之间的差距，本文引入了一种等价的“消失导数”形式。该形式将 $\nabla u = 0$ 且 $D^2 u = 0$ 的点进行单独处理，规定时间导数仅需满足反应项。命题 3.7 确立了这种形式与标准几何粘性形式之间的等价性。
4. 通过不动点与比较原理证明适定性：
   • 比较原理： 定理 5.2 和 5.3 建立了关于几何方程（带有单调反应项）和正则化方程的比较原理。这些原理依赖于反应项 $R_\eta$ 的单调性和曲率算子的结构性质。
   • Perron 构造： 利用比较原理以及通过边界距离估计（引理 5.9）和初始时刻估计（引理 5.10）构建的严格子/超解，作者应用 Perron 方法（引理 5.8）证明了唯一连续粘性解的存在性。
   • 收缩映射： 通过在加权上确界范数下的收缩映射论证处理反应项（引理 7.1），确保了全局适定性以及对初始数据的非扩张依赖性。
5. 收敛性分析： 利用分离主算子与非齐次反应项的稳定性论证（引理 6.12），证明了正则化解 $u_{\varepsilon,\nu}$ 向唯一粘性解 $u$ 的收敛性。这使得作者能够处理零梯度极限处的差异，通过展示正则化解满足消失导数形式（该形式与几何形式等价）来实现这一点。

主要结果

• 全局适定性： 在假设 3.1（光滑、严格平均凸边界；光滑、正值的边缘系数 $G$；光滑、非负的相互作用权重 $\Lambda$；相容的初始数据）下，极限加权几何方程存在唯一的全局连续粘性解（定理 3.5）。
• 区间保持性： 解保持在 $[0, 1]$ 区间内，这对于将 $u$ 解释为归一化分割函数（其边界为水平集 $\{u=0.5\}$）至关重要。
• 非扩张稳定性： 解在 $L^\infty$ 范数下对初始数据具有非扩张依赖性，确保了初始分割中的扰动不会被放大。
• 零梯度选择： 本文指出，正则化算子在零梯度处的直接极限并不会选择预期的几何算子（命题 3.6）。它确立了只有通过特定的消失导数形式才能恢复正确的几何极限（命题 3.7）。
• 一致收敛性： 对于任何序列 $(\varepsilon_n, \nu_n) \to (0, 0)$，正则化解在 $[0, T] \times \Omega$ 上一致收敛于唯一的粘性解（定理 3.8）。至关重要的是，这种收敛性与 $\varepsilon$ 和 $\nu$ 趋向于零的相对速率无关。
• 边界连续性： 在边界为严格平均凸且初始数据在边界附近消失的假设下，粘性解可以连续扩展至边界并达到齐次 Dirichlet 数据。

意义与主张
本文声称为用于 3D 及 3D+时间显微成像分割的正则化主观表面模型提供了严谨的分析基础。具体而言：

• 它通过证明正则化模型尽管在零梯度处存在直接算子极限不匹配的情况，仍能收敛到预期的几何平均曲率流，从而解决了消失正则化极限中的歧义问题。
• 它确立了该模型在粘性意义下的适定性，保证了存在性、唯一性和稳定性。
• 它证实了 $[0, 1]$ 区间的保持性和非扩张稳定性对于极限几何方程同样成立，验证了该模型在自动及半自动分割（其中物体边界由特定水平集定义）中的适用性。
• 研究结果适用于任意维度 $d \ge 2$，并且独立于正则化参数趋向于零时的具体路径，确保了在 $\varepsilon$ 与 $\nu$ 以不同速率衰减的数值实现中的鲁棒性。

这项工作与一篇分析了固定正参数下经典水平模型的伴随论文 [9] 相辅相成。这两项工作共同为该模型在生物图像分割中的应用提供了完整的理论框架。