《Optik》：Free-form surface discretization through surface embedded circle involutes
编辑推荐：
自由曲面光学元件在现代光学系统中无处不在。表面离散化，即选取一组能够有效表征曲面的离散点集，是多个表面制造工艺（如加工、计量、检测或表征）中的关键步骤。本文提出了一种方法，该方法通过沿嵌入在曲面几何内的圆渐开线对表面进行离散化来生成点，该渐开线通过求解一组测地
自由曲面光学元件在现代光学系统中无处不在。表面离散化，即选取一组能够有效表征曲面的离散点集，是多个表面制造工艺（如加工、计量、检测或表征）中的关键步骤。本文提出了一种方法，该方法通过沿嵌入在曲面几何内的圆渐开线对表面进行离散化来生成点，该渐开线通过求解一组测地线方程计算得到。为了对比，将该方法与一种更简单的变体进行了比较，其中圆渐开线在一个平面内生成，随后其上的点被正交投影到曲面上。这两种方法的结果通过三个典型的自由曲面实例进行了分析。
论文解读文章
研究背景与动机
自由曲面光学在现代光学系统中应用广泛，但其非旋转对称的特性带来了诸多理论挑战，尤其是在处理大量设计参数方面。表面离散化，即从连续曲面中选取一组有代表性的离散点集，是连接设计与制造、计量、检测等下游工序的关键环节。然而，如何在自由曲面上实现高效、均匀的采样，同时兼顾路径规划的经济性，是一个亟待解决的问题。特别是在计算机数控扫描系统中，螺旋型路径因其效率和良好的表面质量而被广泛采用，而圆渐开线已被证明是在平面上实现均匀采样并最小化路径长度的最优曲线之一。本研究旨在将这一概念推广至任意自由曲面，解决现有方法在处理高度非旋转对称曲面时采样不均匀的局限性。
研究方法与技术路线
研究人员提出了两种表面离散化方法：表面渐开线法（Surface Involute method, SI）和平面渐开线法（Planar Involute method, PI）。SI法是核心技术，其关键在于计算嵌入在曲面几何内的圆渐开线。首先，在参考平面（Monge patch）上选定一个平面圆，并将其投影到曲面上形成一条闭合曲线。随后，通过数值求解二阶常微分方程组——测地线方程，来计算该闭合曲线的渐开线。具体而言，研究人员利用Monge patch参数化下的Christoffel符号推导了测地线方程，并使用MATLAB内置的可变步长Runge-Kutta求解器（ode45）进行求解。最后，对生成的沿渐开线的非均匀点集进行样条插值，以获得弧长均匀的离散点。作为对照，PI法则直接在参考平面上生成圆渐开线并进行均匀采样，然后将这些点沿参考平面的法向（z轴）正交投影到曲面上。两种方法均使用“外壳”（crust）算法将点云转化为三角网格，并通过计算每个三角形面积与平均面积的比值来量化采样的均匀性。
研究结果
测试曲面与采样可视化（图1）：研究人员选用了三种解析函数描述的典型自由曲面进行测试：曲面A为准旋转对称的Alvarez型曲面，曲面B为包含中等非轴对称分量的Alvarez型曲面，曲面C为极端非旋转对称的自由曲面，其曲率符号变化剧烈。图1直观展示了在低采样密度下，SI法（红点）和PI法（黑点）生成的离散点分布。结果表明，对于曲面A，两种方法差异不大；但对于曲面B和C，PI法产生的点集在边缘区域明显偏离均匀分布。
网格均匀性视觉评估（图2）：通过生成的三角网格图可以观察到，对于曲面A，两种方法的网格单元大小都较为均匀，但PI法在边缘处出现了稍大的三角形。而对于曲面B和C，PI法的网格非均匀性显著增强，尤其在曲面C上表现最为明显，而SI法则始终维持了较好的均匀性。
定量统计分析（图3）：研究人员利用箱线图对不同采样密度（由参数h控制）下的网格单元面积比进行了统计分析。结果显示，随着h减小（采样密度增加），SI法的中位数趋近于1，四分位距显著缩小。例如，对于曲面A（h=1mm），第一和第三四分位数分别为0.9965和1.0048；对于曲面C，则为0.8376和1.1612。相比之下，PI法的箱线图并未显示出随h减小而改善的趋势，且其中位数和四分位距均表明其网格均匀性远逊于SI法。
计算成本分析（表1）：研究人员在Intel Core i7-9700处理器、32GB RAM环境下运行MATLAB R2024a代码，记录了不同采样密度下两种方法的执行时间。结果表明，PI法的计算成本极低（毫秒级），而SI法的耗时则显著更长，且强烈依赖于采样密度。例如，当h=0.2mm时，SI法耗时约10-11分钟（约650-680秒）；而当h≥0.7mm时，耗时可降至10秒以内。值得注意的是，不同曲面类型对计算成本的影响不大。
讨论与结论
讨论部分指出，SI法的核心优势在于通过求解局部测地线方程，生成了真正嵌入曲面几何的圆渐开线，从而确保了采样点在曲面上的均匀分布。研究强调，虽然本研究中使用了圆形作为生成渐开线的基曲线，但这并非强制选择。对于极度偏离旋转对称的曲面，选用更能适应曲面自身对称性的生成曲线可能更为合适，尽管这可能会影响中心区域的均匀性。对于偏离旋转对称不极端的情况（如曲面A和B），圆形是一个良好选择。即使在极端情况（曲面C）下，SI法仍能获得合理的离散化结果，但未来可通过选择其他生成曲线来改进算法。
结论部分翻译如下：研究表明，当待离散化的曲面仅轻微偏离旋转对称时，可以采用平面渐开线法（PI），即在参考平面上生成圆渐开线后再投影到曲面上。然而，对于表面渐开线法（SI），均匀性的偏差会沿着曲面均匀分布，尽管由于所选生成圆半径的原因会呈现出螺旋结构。因此，通过减小h（即提高离散化密度）可以显著降低这种效应。相反，在平面渐开线法中，由于忽略了曲面几何，在曲率较大且变化剧烈的区域（如曲面边缘）不可避免地会出现非均匀性问题。未来的工作将致力于深化渐开线方法在高度非旋转对称自由曲面场景中的应用能力。
