隐式曲面形状交互调整及新方法.docVIP

一种隐式曲面交互调整的新方法 周元峰1 伯彭波3 张彩明1，2 1. 山东大学 计算机科学与技术学院，济南，250101 2. 山东经济学院 计算机科学与技术学院，济南，250014 3. 香港大学 计算机科学系，香港 摘 要：提出了一种对隐式曲面形状进行交互调整的新方法, 为隐式曲面的调整两种交互工具, 分别是对曲面上点的位置调整和法向调整. 该方法以调整后的位置和法向为新曲面的插值建立目标函数, 极小化该目标函数求解曲面参数的变化量, 从而确定新的隐式曲面. 从优化角度对极小化问题进行研究, 分别采用牛顿法和SQP（Sequential Quadratic Programming）方法求解该非线性优化问题. 在调整过程中用粒子的方法对隐式曲面进行绘制, 实现了对隐式曲面形状进行实时交互调整. 最后用实例说明了新方法的有效性. 关键词：隐式曲面；交互调整；优化; 粒子系统; 目标函数 该方法约束一组粒子在曲面上这些粒子自动对曲面进行采样粒子之间的排斥力使粒子均匀地分布在曲面上而且该方法能够允许用户交互调整指定粒子的位置从而改变曲面的形状本文提出了交互调整曲面形状的方法，在显示的时候用Witkin[5]的粒子表示方法（如图1所示），这里先简短介绍Witkin的粒子方法。 Witkin的隐式曲面建模方法是基于粒子采样系统来显示和控制一个隐式曲面。在隐式曲面上约束一个圆点使得： 此处可以看作是浮动于隐式曲面上的一个粒子（particle），利用粒子之间的互斥力将粒子均匀的分布于曲面表面，简单的互斥能够使得静态曲面得到很好的粒子分布，但是需要一个较好的采样作为前提条件。Witkin采样方法采用粒子动态分裂/死亡机制将粒子约束在曲面上，初始给定一个随机粒子和一个比较大的互斥半径，用迭代方法对粒子进行分裂或者死亡，若粒子的互斥半径大于设定的半径阈值或者粒子能量大于设定的能量阈值，则将该粒子一分为二，并减少粒子之间的互斥半径；反之，若粒子的互斥半径过小，则该粒子死亡。迭代重复上述步骤，直至所有粒子满足能量和互斥半径约束条件。该方法可以很好的达到收敛效果，这样隐式曲面的调整可以通过调整粒子来实现，达到了交互的目的。 隐式曲面的形变根据隐式曲面系数的变化速度确定： 其中表示隐式曲面系数的变化速度，为约束粒子的个数，为要得到的无约束系数变化速度。 算法基本思想 定义隐式曲面表达式为： （1） 其中定义为曲面上的点坐标，为隐式曲面表达式系数向量。本文提出两种对于任意隐式曲面的形状进行调整的工具，分别是对曲面上指定点的位置调整和法向的调整。在调整过程中，为了实时显示调整结果，我们采用Witkin的粒子方法对隐式曲面进行实时绘制。 隐式曲面表达式（1）可以看作函数的零值面。当用户改变了隐式曲面上的某一个顶点的位置或者法向之后，曲面形状就会随之变化，新的隐式曲面插值由新顶点的位置或法向来指定。本文的方法以新顶点的位置或法向作为曲面形状改变的目标插值，建立新的目标函数，极小化目标函数求解新的隐式曲面系数变化量，从而得到形状改变后的新的隐式曲面。这里，该求解过程看作一个优化过程，以插值条件建立目标函数，进行优化。第3节详细给出了插值目标函数的建立和求解方法。 隐式曲面交互调整新方法 首先介绍曲面交互调整的工具，如图1所示，图中是采用粒子系统绘制的两个blobby球体组成的blobby曲面。图中椭圆区域所标识的蓝色圆盘为曲面上点的位置调整手柄，蓝色圆锥为曲面上点的法向调整手柄。用户选中要调整的圆盘或圆锥，调整到希望的位置或方向。隐式曲面自动调整形状，插值于用户的设定位置或法向。 4.1 基于顶点位置的调整 任意指定隐式曲面上某一顶点，通过交互方式将的位置移动至新位置，这时发生形变后的隐式曲面方程为： 其中为隐式曲面方程系数的变化量。 利用隐式曲面是曲面表达式的零值面的性质，我们定义目标函数如下： （2） 其中为方程未知量。目标函数（2）的几何意义可以看作是调整后的顶点到新曲面的代数距离平方。通过极小化寻找通过新顶点的零值面。即求解: （3） 4.2 基于顶点法向的调整 基于顶点位置的调整对于曲面形状设计来说比较直观，但是调整方法过于单一。本文还给出了基于曲面顶点法向的调整方法。用户指定曲面上的一点，交互地将该点处法向的方向调整到一个新的位置。 （1）无约束法向调整： 假设对的法向调整之后的新法向为，以调整后的曲面在处的法向与用户给定法向同向为目标，建立以下目标函数： （4） 其中表示在处的梯度，为隐式曲面系数变化量。 极小化目标函数（4）