《Mathematics and Computers in Simulation》：Detecting and filling holes to reconstructing explicit surfaces using bicubic splines
编辑推荐：
本文研究了一种针对显式曲面在多个未知子域存在信息缺失情况下的空洞检测与填充方法。该方法通过在双三次样条（bicubic splines）这一合适有限维空间中构建最小化问题，实现了对包含缺失信息的显式曲面的近似重建。研究证明了所提问题解的存在性与唯一性，给出了计算方法，并建立了方法的收敛性结果。通过图形与数值算例，验证了该方法在曲面重建领域的有效性，为解决CAD/CAM、地球科学、机器人视觉、三维人体扫描、牙科重建、逆向工程、卫星与雷达图像重建等众多计算图形学领域中的曲面缺失问题提供了新的工具。
在计算图形学、计算机辅助设计与制造（CAD/CAM）、地球科学、机器人视觉、三维人体扫描、牙科重建、逆向工程、以及从卫星和雷达信息中重建图像等诸多领域中，一个常见而棘手的问题是：如何准确、高效地修复三维模型或曲面上的“空洞”？这些空洞，即模型表面一个或多个区域的数据缺失，可能源于扫描过程中的物体遮挡、反射，或原始数据本身的不完整。面对一个已知其部分点云数据、但存在未知缺失区域的显式曲面，我们能否自动检测出空洞的位置，并基于周围已知信息，以高精度、高保真的方式“填补”这些空白，从而重建出一个完整、光滑的表面？这不仅是理论上的挑战，更是驱动众多工程与科学应用发展的核心需求。
以往的研究者们已尝试了多种途径。例如，有些方法专注于在满足特定形状约束条件下填充空洞，通过构造函数来估计缺失区域内部未知的导数；另一些研究则利用参数三次样条函数空间中的二次泛函优化来进行曲线插值，或通过双三次样条函数空间中的二次泛函优化实现曲面或数据集的插值。还有工作致力于寻找能最小化特定泛函的光滑二元样条函数，以满足包括体积限制在内的多种条件。此外，利用双二次样条或变分样条来逼近广义偏移曲面的方法也有所发展。然而，这些现有方法大多侧重于填补过程中的单一环节。A. Kouibia和M. Pasadas在发表于《Mathematics and Computers in Simulation》的这项研究中指出，首次将基于网格划分点分布的孔洞检测、孔洞外的平滑双三次样条拟合以及用于填补孔洞的变分插值这三个核心方面结合起来进行研究，构成了本工作的主要新颖之处。
为了系统解决上述问题，研究人员设计并阐述了一套完整的混合近似方法。整个流程直观地体现在其核心示意图中：首先，给定一个存在孔洞的显式函数的三维点集及其在定义域上的投影；接着，识别出定义域上存在信息缺失的子区域（即孔洞）；然后，基于已知点（函数图像上的给定点），通过平滑方法在孔洞外部重建曲面；最后，根据孔洞外已重建的曲面，通过插值方法填补孔洞内部的函数值，从而完成整个曲面的完全重建。
具体而言，本研究采用了几个关键技术方法：1. 孔洞检测：基于给定的、分布在矩形定义域R = (a, b) × (c, d)上的有限点集A（该点集位于孔洞之外，并满足一定的稠密性条件），通过对定义域进行规则的矩形网格K?划分，将那些与已知点集A无交集的网格单元判定为属于孔洞区域，从而近似检测出未知孔洞H的位置和形状。理论证明，当网格足够细密时，检测出的近似孔洞H?N可以无限逼近真实孔洞H。2. 孔洞外曲面平滑逼近：在检测出孔洞区域后，研究在由两组节点集Tn和Tm生成的双三次样条函数空间S3(Tn×Tm) 中求解一个最小化问题。该问题旨在寻找一个样条函数σ，使其在已知数据点AM上对给定函数值f的拟合误差（通过离散l2范数度量）与函数本身的光滑性（通过涉及二阶偏导的希尔伯特范数|·|2, R\H?N度量）之间取得平衡，其中ρ>0为平滑参数。该问题解的存在唯一性得到了证明，并给出了其计算方式。3. 孔洞内变分插值填充：在获得孔洞外的近似样条σ后，研究转向填补孔洞H?N。通过在索伯列夫空间H3(H?N)中求解一个变分插值问题来实现：寻找函数u，使其在孔洞边界?H?N上与外部近似样条σ的值及其一阶法向导数相匹配（即满足狄利克雷和诺伊曼边界条件），同时使函数在孔洞内部的三阶希尔伯特半范数|u|3, H?N最小化。该内问题同样被证明存在唯一解，并且通过将外部近似σ的边界信息作为输入，可以得到孔洞内的完整填充函数。
研究结果部分通过严谨的数学证明和数值实验展示了方法的有效性：
- •
解的存在唯一性与收敛性：研究首先在严格的数学框架下，分别证明了孔洞外平滑问题和孔洞内插值问题的解的存在性与唯一性。这为整个方法奠定了坚实的理论基础，确保了计算过程的稳定性和结果的确定性。进而，论文建立了一个关键的收敛性结果：当已知数据点集AM变得越来越稠密（M → +∞），并且用于构建样条空间的网格也越来越精细（n, m → +∞）时，由本方法得到的全局近似解（即孔洞内外的组合函数）将在索伯列夫空间H3(R) 的意义下收敛到真实的未知函数f。这从理论上保证了方法的可靠性，即只要有足够多的已知数据，就能获得任意精度的重建曲面。
- •
数值与图形示例验证：为了直观展示方法的实用性和效果，论文提供了具体的数值与图形示例。通过对具有已知解析表达式的测试函数（如f(x, y) = sin(x)cos(y)）人为构造“孔洞”，并应用所述方法进行重建，可以定量比较重建结果与真实函数之间的误差。示例表明，该方法能够有效地检测出不同形状（如矩形、圆形）的孔洞，并在孔洞内外生成光滑的逼近曲面。误差分析显示，重建误差随着数据点数量的增加和网格的细化而减小，这与理论收敛性结论一致。图形结果清晰展示了从带孔洞的离散点云到完整光滑曲面的重建过程，验证了方法在处理实际曲面缺失问题中的潜力。
结论与讨论部分对研究进行了总结并展望其意义。本研究成功提出并分析了一种用于显式曲面空洞检测与填充的混合近似方法。该方法创新性地将基于点分布网格划分的孔洞检测、基于双三次样条空间的孔洞外平滑以及基于变分原理的孔洞内插值这三个关键步骤整合到一个连贯的框架中。研究不仅提供了完整的算法流程，更从数学上严格证明了各步骤解的存在唯一性，并建立了整个方法的收敛性理论，确保了方法的数学严密性和可靠性。
这项工作的意义重大。首先，其方法学上的创新性在于首次实现了检测、外部平滑与内部插值三个环节的协同研究，为解决复杂的曲面缺失问题提供了更系统、更强大的工具。其次，广泛的应用潜力使其能够直接服务于众多前沿科技领域，包括但不限于：提高CAD/CAM系统中模型修复的自动化程度与精度；辅助地球科学研究中从不完整勘测数据重建地形曲面；提升机器人视觉对不完整三维场景的理解能力；改善三维人体与牙科扫描中因遮挡导致的数据缺失问题；优化逆向工程中从破损实物到完整数字模型的重建流程；以及增强从卫星或雷达的不完整探测数据中复原清晰图像的能力。最后，研究所依托的双三次样条工具具有良好的数值计算性质和光滑性（C2连续），保证了重建曲面的高质量，而变分插值方法则能自然地融入边界约束，实现孔洞内外曲面的光滑过渡。
综上所述，这项由A. Kouibia和M. Pasadas完成的研究，为计算几何与图形学中的曲面重建经典难题提供了一个理论坚实、步骤清晰、且极具应用前景的解决方案。其提出的混合框架和严格的数学分析，将为后续研究者开发更高效、更稳健的曲面处理算法提供重要的借鉴和基础。
