计算机图形学4.ppt

《计算机图形学4.ppt》由会员分享，可在线阅读，更多相关《计算机图形学4.ppt（86页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、 3.2 Bezier 曲线与曲面曲线与曲面清华大学清华大学 计算机图形学计算机图形学 由于几何外形设计的要求越来越高,传统的曲线曲面表示方法,已不能满足用户的需求。1962年，法国雷诺汽车公司的P.E.Bezier构造了一种以逼近为基础的参数曲线和曲面的设计方法，并用这种方法完成了一种称为UNISURF 的曲线和曲面设计系统，1972年，该系统被投入了应用。清华大学清华大学 计算机图形学计算机图形学Bezier方法将函数逼近同几何表示结合 起来，使得设计师在计算机上就象使用作图工具一样得心应手。典故日本的穗板：天上掉下来 为边向量 清华大学清华大学 计算机图形学计算机图形学剑桥的 Fores

2、t 常庚哲：中国的Bezier，曲面凸性梁友栋：几何连续的浙大学派，梁叶郑马刘鼎元：实用的几何连续条件Hoschek的故事刘汪佳话纪念Bezier的CAGD专辑 清华大学清华大学 计算机图形学计算机图形学 清华大学清华大学 计算机图形学计算机图形学3.2.1 Bezier曲线的定义和性质曲线的定义和性质1定义 给定空间n+1个点的位置矢量Pi（i=0，1，2，n），则Bezier曲线可定义为：清华大学清华大学 计算机图形学计算机图形学 其中，Pi构成该Bezier曲线的特征多边形，Bi,n(t)是n次Bernstein基函数：0=1,0!=1清华大学清华大学 计算机图形学计算机图形学2Betn

3、stein基函数的性质 （1）正性 （2）端点性质 清华大学清华大学 计算机图形学计算机图形学（3）权性 由二项式定理可知：清华大学清华大学 计算机图形学计算机图形学（4）对称性 因为 清华大学清华大学 计算机图形学计算机图形学(5）递推性。即高一次的Bernstein基函数可由两个低一次的Bernstein调和函数线性组合而成。因为，清华大学清华大学 计算机图形学计算机图形学（6）导函数 （7）最大值。在 处达到最大值。清华大学清华大学 计算机图形学计算机图形学（8）升阶公式 清华大学清华大学 计算机图形学计算机图形学（9）积分清华大学清华大学 计算机图形学计算机图形学3Bezier曲线的性

4、质（1）端点性质 a）曲线端点位置矢量 由Bernstein基函数的端点性质可以推得，当t=0时，P(0)=P0；当t=1时，P(1)=Pn。由此可见，Bezier曲线的起点、终点与相应的特征多边形的起点、终点重合。清华大学清华大学 计算机图形学计算机图形学 b)切矢量 因为，所以当t=0时，P(0)=n(P1-P0)，当t=1时，P(1)=n(Pn-Pn-1)，这说明Bezier曲线的起点和终点处的切线方向和特征多边形的第一条边及最后一条边的走向一致。清华大学清华大学 计算机图形学计算机图形学c.）二阶导矢当t=0时，当t=1时，上式表明：2阶导矢只与相邻的3个顶点有关，事实上，r阶导矢只与

5、（r+1）个相邻点有关，与更远点无关。将 、及 、代入曲率公式 ，可以得到Bezier曲线在端点的曲率分别为：清华大学清华大学 计算机图形学计算机图形学d.）k阶导函数的差分表示n次Bezier曲线的k阶导数可用差分公式为：其中高阶向前差分矢量由低阶向前差分矢量递推地定义：例如：清华大学清华大学 计算机图形学计算机图形学（2）对称性。由控制顶点 构造出的新Bezier曲线，与原Bezier曲线形状相同，走向相反。因为：这个性质说明Bezier曲线在起点处有什么几何性质，在终点处也有相同的性质。清华大学清华大学 计算机图形学计算机图形学（3）凸包性由于 ，且 ，这一结果说明当t在0，1区间变化时

6、，对某一个t值，P(t)是特征多边形各顶点的加权平均，权因子依次是 。在几何图形上，意味着Bezier曲线P(t)在 中各点是控制点Pi的凸线性组合，即曲线落在Pi构成的凸包之中，如图3.1.9所示。清华大学清华大学 计算机图形学计算机图形学（4）几何不变性。这是指某些几何特性不随坐标变换而变化的特性。Bezier曲线位置与形状与其特征多边形顶点 的位置有关，它不依赖坐标系的选择。清华大学清华大学 计算机图形学计算机图形学（5）变差缩减性。若Bezier曲线的特征多边形 是一个平面图形，则平面内任意直线与C(t)的交点个数不多于该直线与其特征多边形的交点个数，这一性质叫变差缩减性质。此性质反映

7、了Bezier曲线比其特征多边形的波动还小，也就是说Bezier曲线比特征多边形的折线更光顺。清华大学清华大学 计算机图形学计算机图形学（6）仿射不变性对于任意的仿射变换A：即在仿射变换下，的形式不变。清华大学清华大学 计算机图形学计算机图形学3.2.2 Bezier曲线的曲线的递推递推(de Casteljau)算法算法 计算Bezier曲线上的点，可用Bezier曲线方程，但使用de Casteljau提出的递推算法则要简单的多。如下图所示，设 、是一条抛物线上顺序三个不同的点。过 和 点的两切线交于 点,在 点的切线交 和 于 和 ，则如下比例成立：这是所谓抛物线的三切线定理。(示意图见

8、下页)清华大学清华大学 计算机图形学计算机图形学 清华大学清华大学 计算机图形学计算机图形学当P0，P2固定，引入参数t，令上述比值为t:(1-t)，即有：t从0变到1，第一、二式就分别表示控制二边形的第一、二条边，它们是两条一次Bezier曲线。将一、二式代入第三式得：当t从0变到1时，它表示了由三顶点P0、P1、P2三点定义的一条二次Bezier曲线。并且表明：这二次Bezier曲线P20可以定义为分别由前两个顶点（P0，P1）和后两个顶点（P1，P2）决定的一次Bezier曲线的线性组合。依次类推，由四个控制点定清华大学清华大学 计算机图形学计算机图形学义的三次Bezier曲线P30可被

9、定义为分别由(P0，P1，P2)和(P1，P2，P3)确定的二条二次Bezier曲线的线性组合，由(n+1)个控制点Pi(i=0,1,.,n)定义的n次Bezier曲线Pn0可被定义为分别由前、后n个控制点定义的两条(n-1)次Bezier曲线 P0n-1与P1n-1的线性组合：由此得到Bezier曲线的递推计算公式：这便是著名的de Casteljau算法。用这一递推公式，在给定参数下，求Bezier曲线上一点P(t)非常有效。上式中：是定义Bezier清华大学清华大学 计算机图形学计算机图形学曲线的控制点，即为曲线 上具有参数t的点。de Casteljau算法稳定可靠，直观简便，可以编出

10、十分简捷的程序，是计算Bezier曲线的基本算法和标准算法。当n=3时，de casteljau算法递推出的Pki呈直角三角形，对应结果如图3.1.11所示。从左向右递推，最右边点P30即为曲线上的点。清华大学清华大学 计算机图形学计算机图形学 清华大学清华大学 计算机图形学计算机图形学这一算法可用简单的几何作图来实现。给定参数 ，就把定义域分成长度为 的两段。依次对原始控制多边形每一边执行同样的定比分割，所得分点就是由第一级递推生成的中间顶点 ,对这些中间顶点构成的控制多边形再执行同样的定比分割，得第二级中间顶点 。重复进行下去，直到n级递推得到一个中间顶点 即为所求曲线上的点 ，如图3.1

11、.12所示。清华大学清华大学 计算机图形学计算机图形学 清华大学清华大学 计算机图形学计算机图形学3.2.3 Bezier曲线的曲线的拼接拼接 几何设计中，一条Bezier曲线往往难以描述复杂的曲线形状。这是由于增加由于特征多边形的顶点数，会引起Bezier曲线次数的提高，而高次多项式又会带来计算上的困难，实际使用中，一般不超过10次。所以有时采用分段设计，然后将各段曲线相互连接起来，并在接合处保持一定的连续条件。下面讨论两段Bezier曲线达到不同阶几何连续的条件。清华大学清华大学 计算机图形学计算机图形学 给定两条Bezier曲线P(t)和Q(t)，相应控制点为Pi(i=0,1,.,n)和

12、Qj(j=0,1,.,m)，且令 ，如图3.1.13所示，我们现在把两条曲线连接起来。图3.1.13 Bezier曲线的拼接清华大学清华大学 计算机图形学计算机图形学（1）要使它们达到G0连续的充要条件是：Pn=Q0；（2）要使它们达到G1连续的充要条件是：Pn-1，Pn=Q，Q1三点共线，即：（3）要使它们达到G2连续的充要条件是：在G1连续的条件下，并满足方程 。我们将 、和 ，、代入，并整理，可以得到：选择 和 的值，可以利用该式确定曲线段 的特征多边形顶点 ，而顶点 、已被 连续条件所确定。要达到 连续的话，只剩下顶点 可以自由选取。清华大学清华大学 计算机图形学计算机图形学如果从上式

13、的两边都减去 ，则等式右边可以表示为 和 的 线性组合：这表明 、和 五点共面，事实上，在接合点两条曲线段的曲率相等，主法线方向一致，我们还可以断定：位于直线 的同一侧。清华大学清华大学 计算机图形学计算机图形学3.2.4 Bezier曲线的升阶与降阶曲线的升阶与降阶1Bezier曲线的升阶 所谓升阶是指保持Bezier曲线的形状与定向不变，增加定义它的控制顶点数，也即是提高该Bezier曲线的次数。增加了控制顶点数，不仅能增加了对曲线进行形状控制的灵活性，还在构造曲面方面有着重要的应用。对于一些由曲线生成曲面的算法，要求那些曲线必须是同次的。应用升阶的方法，我们可以把低于最高次数的的曲线提升

14、到最高次数，而获得同一的次数。曲线升阶后，原控制顶点会发生变化。下面，我们来计算曲线提升一阶后的新的控制顶点。设给定原始控制顶点 ，定义了一条n次Bezier曲线：清华大学清华大学 计算机图形学计算机图形学增加一个顶点后，仍定义同一条曲线的新控制顶点为 ，则有：对上式左边乘以 ，得到：比较等式两边 项的系数，得到：化简即得：其中 。清华大学清华大学 计算机图形学计算机图形学此式说明：新的控制顶点 是以参数值 按分段线性插值从原始特征多边形得出的。升阶后的新的特征多边形在原始特征多边形的凸包内特征多边形更靠近曲线。三次Bezier曲线的升阶实例如图3.1.14所示。清华大学清华大学 计算机图形学

15、计算机图形学 清华大学清华大学 计算机图形学计算机图形学2Bezier曲线的降阶 降阶是升阶的逆过程。给定一条由原始控制顶点 定义的n次Bezier曲线，要求找到一条由新控制顶点 定义的n-1次Bezier曲线来逼近原始曲线。假定 是由 升阶得到，则由升阶公式有：从这个方程可以导出两个递推公式：和清华大学清华大学 计算机图形学计算机图形学两种降阶格式Forrest 格式Farin格式清华大学清华大学 计算机图形学计算机图形学降阶逼近的文献M.A.Watkins and A.J.Worsey,Degree reduction of Bzier curves,Computer Aided Desi

16、gn,20(7),1988,398-405胡事民、孙家广、金通光、汪国昭，Approximate degree reduction of Bezier curves,Tsinghua Science and Technology,No.2,1998,997-1000.雍俊海、胡事民、孙家广、谭新宇，Degree reduction of B-spline curves,Computer Aided Geometric Design,2001,Vol.13,NO.2,2001,117-127.清华大学清华大学 计算机图形学计算机图形学3.2.5 Bezier曲面曲面 基于Bezier曲线的讨论，

17、我们可以方便地可以给出Bezier曲面的定义和性质，Bezier曲线的一些算法也可以很容易扩展到Bezier曲面的情况。清华大学清华大学 计算机图形学计算机图形学 1定义设 为 个空间点列，则 次张量积形式的Bezier曲面定义为：其中 ，是Bernstein基函数。依次用线段连接点列中相邻两点所形成的空间网格，称之为特征网格。清华大学清华大学 计算机图形学计算机图形学Bezier曲面的矩阵表示式是：在一般实际应用中，不大于4。清华大学清华大学 计算机图形学计算机图形学 2性质 除变差减小性质外，Bezier曲线的其它性质可推广到Bezier曲面：（1）Bezier曲面特征网格的四个角点正好是

18、Bezier曲面的四个角点，即 ，。（2）Bezier曲面特征网格最外一圈顶点定义Bezier曲面的四条边界；Bezier曲面边界的跨界切矢只与定义该边界的顶点及相邻一排顶点有关，且 、和 （图3.1.15打上斜线的三角形）；其跨界二阶导矢只与定义该边界的及相邻两排顶点有关；。清华大学清华大学 计算机图形学计算机图形学 （3）几何不变性。（4）对称性。（5）凸包性。清华大学清华大学 计算机图形学计算机图形学 3Bezier曲面片的拼接如图3.1.16所示，设两张mn次Bezier曲面片分别由控制顶点 和 定义。清华大学清华大学 计算机图形学计算机图形学如果要求两曲面片达到 连续，则它们有公共的

19、边界，即：（3.1.10）于是有 如果又要求沿该公共边界达到 连续，则两曲面片在该边界上有公共的切平面，因此曲面的法向应当是跨界连续的，即：（3.1.11）清华大学清华大学 计算机图形学计算机图形学下面来研究满足这个方程的两种方法。（1）鉴于（3.1.10）式，（3.1.11）式最简单的解是：（3.1.12）这相当于要求合成曲面上v为常数的所有曲线，在跨界时有切向的连续性。为了保证等式两边关于v的多项式次数相同，必须取 （一个正常数）。于是有：即 清华大学清华大学 计算机图形学计算机图形学（2）由于（3.1.12）式使得两张曲面片在边界达到 连续时，只涉及曲面 和 的两列控制顶点，比较容易控制

20、。用这种方法匹配合成的曲面的边界，u向和v向是光滑连续的。实际上，该式的限制是苛刻的。是否存在跟简单而合理的方法？清华大学清华大学 计算机图形学计算机图形学 为了构造合成曲面时有更大的灵活性，Bezier在1972年放弃把（3.1.12）式作为 连续的条件，而以 （3.1.13）来满足（3.1.11）式，这仅仅要求 位于 和 所在的同一个平面内，也就是曲面片 边界上相应点处的切平面，这样就有了大得多的余地，但跨界切矢在跨越曲面片的边界时就不再连续了。同样，为了保证等式两边关于v的多项式次数相同，须为任意正常数，是v的任意线性函数。清华大学清华大学 计算机图形学计算机图形学4递推(de Cast

21、eljau)算法 Bezier曲线的递推(de Casteljau)算法，可以推广到Bezier曲面的情形。若给定Bezier曲面特征网格的控制顶点 和一对参数值 ，则：(3.1.14)清华大学清华大学 计算机图形学计算机图形学其中 (3.1.15)或 (3.1.16)清华大学清华大学 计算机图形学计算机图形学 (3.1.15)与(3.1.16)中的下标 的变化范围已在(3.1.14)式中给出。上面给出了确定曲面上一点的两种方案。当按(3.1.15)式方案执行时，先以u参数值对控制网格u向的n+1个多边形执行曲线de Casteljau算法，m级递推后，得到沿v向由n+1个顶点 构成的中间多边

22、形。再以v参数值对它执行曲线的de Casteljau算法，n级递推以后，得到一个 ，即所求曲面上的点清华大学清华大学 计算机图形学计算机图形学 也可以按(3.1.16)式方案执行，先以v参数值对控制网格沿v向的m+1个多边形执行n级递推，得沿u向由m+1个顶点 构成的中间多边形。再以u参数值对它执行n级递推，得所求点 。清华大学清华大学 计算机图形学计算机图形学 之所以有这样的算法，是因为清华大学清华大学 计算机图形学计算机图形学3.1.2.6 三边三边Bezier曲面片曲面片 与上一节定义在矩形域上的Bezier曲面片不同，本节介绍的三边Bezier曲面片是定义在三边形域上的，如图3.1.

23、17所示，为了便于区分，我们把上一节介绍的Bezier曲面片称为四边Bezier曲面片。三边曲面片能够较好地适应不规则与散乱数据的几何造型及适合有限元分析中的三边元素的需要。清华大学清华大学 计算机图形学计算机图形学1三角域内点的表示 三角域内一点可以用面积坐标（或重心坐标）来表示，如图3.1.18所示。G是三角形ABC内的任意一点，其面积坐标为 。令三角形ABC面积为s，三角形GBC面积为 ，三角形GCA面积为 ，三角形GAB面积为 ，则：，这里所指的三角形面积是有向面积，按顶点字母顺序，顺时针旋转为正，逆时针旋转为负。清华大学清华大学 计算机图形学计算机图形学 若A=，B=，C=，则：三个

24、坐标分量u、v和w只有两个是独立的，因为u+v+w=1.三角形ABC称为域三角形，或称为三角域。清华大学清华大学 计算机图形学计算机图形学 2三角域上的Bernstein基 单变量的n次的Bernstein基 由 的二项式展开各项组成。双变量张量积的Bernstein基由两个单变量的Bernstein基各取其一的乘积组成。而定义在三角域上的双变量n次的Bernstein基由 的展开式各项组成。清华大学清华大学 计算机图形学计算机图形学Bernstein基函数：其中i+j+k=n，且i，j，k0。可见，三角域上n次Bernstein基共包含了 个基函数，可以用一个三角阵来排列这些基函数。例如，n

25、=2时如图3.1.19所示。清华大学清华大学 计算机图形学计算机图形学 清华大学清华大学 计算机图形学计算机图形学 三角域按Bernstein基的三角阵列相应划分成子三角域，其中诸直线交点同样地称为节点。节点与基函数一一对应。每个结点也由三个指标确定，如图3.1.20所示，它们分别与三参数u，v，w相联系。三角域上Bernstein基同样具有规范性、非负性与递推性。其递推关系为：清华大学清华大学 计算机图形学计算机图形学3三边Bzier曲面片的定义 使一个基函数联系一个控制顶点，一张n次三边Bezier曲面片必须由构成三角阵列的 个控制顶点 定义。因此，n次Bezier曲面可以定义为：按下标顺

26、序用直线连接控制顶点，就形成了曲面的控制网格，它由三角形组成，网格顶点与三角域的节点一一对应。图3.1.21给出了三次三边Bezier曲面片的一个例子。清华大学清华大学 计算机图形学计算机图形学 清华大学清华大学 计算机图形学计算机图形学 当固定三参数之一时，将得到曲面片上一条等参数线。例如，当w固定，让u独立地变化，则得到一条u线；若让v独立地变化，则得到v线，两者实际是同一条曲线。因此，曲面片上有三族等参数线。当三参数之一为零时，则得曲面片的一条边界线，它由相应那排边界顶点定义，就是一般所指的一条非有理n次Bezier曲线。当三参数之一为1时，则得三边曲面片的一个角点，就是控制网格三角顶点

27、之一。可见，三边Bezier曲面片与四边Bezier曲面片具有类似的性质。清华大学清华大学 计算机图形学计算机图形学思考题：如何计算Bezier曲面上一点的值？清华大学清华大学 计算机图形学计算机图形学 de Casteljau递推。清华大学清华大学 计算机图形学计算机图形学与定义在矩形域上的四边Bezier曲面片的差别在于：1）定义域不同；2）控制网格不同，后者由呈矩形阵列的控制顶点构成；3）同样是两个独立参数，但最高次数不同，后者两个参数的最高次数是互相独立的，可以不同。而三边Bezier曲面片的三个参数的最高次数都是相同的；4）四边Bezier曲面片是张量积曲面，三边Bezier曲面片是

28、非张量积曲面，这是本质差别。清华大学清华大学 计算机图形学计算机图形学4三边Bzier曲面片与四边Bzier曲面片的转化 由于三边Bzier曲面与四边Bzier曲面有不同的基函数和定义方法，当在同一个CAD系统中使用这两种类型的曲面片时，会带来不相容的困难。下面我们给出了两种曲面片的转化方法。清华大学清华大学 计算机图形学计算机图形学 一种方法是将一张三边Bzier曲面片转化为三张相同次数的四边Bezier曲面片，且各曲面片之间能够很好地匹配。如图3.1.22（a），假定三边Bzier曲面片定义域为D，在D的三条边上各取一点（不包括三个顶点）、和 ，再在三角形 内取一点 ，则线段 、和 将D分

29、成三个四边形 、和 ，与三线段对应的曲线将该三边Bzier曲面片分成三张四边Bzier曲面片。下面，我们给出 上的四边Bzier曲面片的表示，清华大学清华大学 计算机图形学计算机图形学 和 上曲面片的表示则与在上类似。我们选取 如图3.1.22（b）所示。为了方便起见，并将 写为 写为 ，则三边Bzier曲面片方程可以写为：清华大学清华大学 计算机图形学计算机图形学清华大学清华大学 计算机图形学计算机图形学如果我们还引入下列运算符：不变运算符 移位运算符 差分运算符则三边Bezier曲面片方程可进一步表示为：清华大学清华大学 计算机图形学计算机图形学 定义在 上的 的裁剪曲面可由nn次的四边B

30、zier曲面片表示，其控制顶点为：其中 分别是 、和 的坐标。于是，定义在 上的四边曲面片可表示为：清华大学清华大学 计算机图形学计算机图形学 第二种方法是将三角Bzier曲面转化为一张退化的张量积Bzier曲面。清华大学清华大学 计算机图形学计算机图形学 退化的张量积曲面的控制顶点可以由下式计算清华大学清华大学 计算机图形学计算机图形学 其中 是升阶算子清华大学清华大学 计算机图形学计算机图形学 证明:对三角曲面T(u,v,w)应用重新参数化三角域变为矩形域 ,于是有清华大学清华大学 计算机图形学计算机图形学 于是有：清华大学清华大学 计算机图形学计算机图形学 应用：在一些 CAD 应用中，

31、我们经常需设计一些三边曲面，但一般的系统只支持矩形曲面片，很自然的，我们会用一条边退化成一点的矩形片表示。而退化曲面是不受欢迎的.因此我们需解决将该退化Bzier曲面非退化矩形曲面的问题。清华大学清华大学 计算机图形学计算机图形学 而利用前面的讨论，用 次三角Bzier曲面 逼近该退化的 次Bzier矩形片 ，非常简单，只需将如下 次Bzier曲线降阶为 次曲线。(8)清华大学清华大学 计算机图形学计算机图形学 因此我们首先将退化曲面逼近为一张三角Bzier 曲面,在利用我们前面提到的方法将该三角Bzier 曲面转化为三张非退化矩形域Bzier曲面。清华大学清华大学 计算机图形学计算机图形学 清华大学清华大学 计算机图形学计算机图形学 如何提高精度？清华大学清华大学 计算机图形学计算机图形学