[论文解读] On the Efficiency of Strategies for Subdividing Polynomial Triangular Surface Patches

Jean Gallier|ArXiv.org|Jun 13, 2006

Advanced Numerical Analysis Techniques参考文献 9被引用 3

一句话总结

本文提出了基于 de Casteljau 算法的多项式三角形曲面片最优且高效的细分算法。它引入了一种四次调用的规则细分方法（已证明为最优），一种三次调用的菱形细分结构，以及一种四次调用的蜘蛛网状细分，可生成六个子三角形，所有方法均在保持几何精度的同时最小化算法调用次数。

ABSTRACT

In this paper, we investigate the efficiency of various strategies for subdividing polynomial triangular surface patches. We give a simple algorithm performing a regular subdivision in four calls to the standard de Casteljau algorithm (in its subdivision version). A naive version uses twelve calls. We also show that any method for obtaining a regular subdivision using the standard de Casteljau algorithm requires at least 4 calls. Thus, our method is optimal. We give another subdivision algorithm using only three calls to the de Casteljau algorithm. Instead of being regular, the subdivision pattern is diamond-like. Finally, we present a ``spider-like'' subdivision scheme producing six subtriangles in four calls to the de Casteljau algorithm.

研究动机与目标

• 最小化细分多项式三角形曲面片所需的 de Casteljau 算法调用次数。
• 开发比朴素方法计算开销更低的高效细分策略。
• 探索超越标准规则细分的替代细分模式——菱形和蜘蛛网状结构。
• 证明规则细分中 de Casteljau 调用次数的理论下限。
• 通过递归细分示例和视觉模式展示其实际适用性。

提出的方法

• 采用四次标准 de Casteljau 算法调用的规则细分策略，通过证明任何方法都无法在少于四次调用内完成，从而证明其最优性。
• 通过递归分割子三角形并选择性保留控制网，实现仅需三次调用的菱形细分结构，利用重心坐标。
• 通过四次调用实现蜘蛛网细分方案，生成六个子三角形，利用质心并递归分割初始三个子三角形。
• 使用 sdecas3 函数计算子三角形的控制网，并选择性丢弃不需要的子网以优化计算。
• 应用极形式和多重仿射对称映射表示多项式曲面，并推导控制点之间的关系。
• 实现递归细分，当省略中心三角形时，可生成类似 Sierpinski 垫片的视觉模式。

实验结果

研究问题

• RQ1 对多项式三角形曲面片进行规则细分，所需的 de Casteljau 算法调用次数最少是多少？
• RQ2 是否可以构造一种非规则细分模式，在少于四次 de Casteljau 算法调用下仍能生成有效划分？
• RQ3 使用四次 de Casteljau 算法调用可生成哪些几何细分模式？它们在效率和结构上如何比较？
• RQ4 如菱形或蜘蛛网状等替代细分方案如何影响生成子三角形的数量和分布？
• RQ5 当省略中心三角形时，递归细分会涌现出何种视觉或分形效应？

主要发现

• 所提出的规则细分算法恰好使用四次 de Casteljau 算法调用，且已被证明为理论最小值，因此为最优。
• 菱形细分策略仅用三次调用即实现有效划分，优于朴素的十二次调用方法。
• 蜘蛛网细分方案在四次调用内生成六个子三角形，展示了新型高效非规则细分模式。
• 当规则细分方法被递归应用且中心三角形被省略时，可在曲面片上生成类似 Sierpinski 垫片的分形图案。
• 所有提出的算法均可通过控制网技术的直接扩展，适用于有理曲面片。
• 利用重心坐标并选择性保留控制网，可在不重复调用 de Casteljau 算法的前提下实现高效计算。