《计算机图形学》课程实验.ppt

《计算机图形学》课程实验 张玉娟 实验一 直线段生成算法实现 1、数值微分法（DDA） 基本思想： 从直线段的起点像素出发, 依据直线的微分方程依次确定描述直线的各个像素点。 微分方程: dy/dx=m, yi+1=yi+m(xi+1-xi) 若 xi+1=xi+1 yi+1=yi+m (0m1) 算法原理： 设一直线段的起点和终点坐标分别为(xs , ys)和(xe , ye)。则直线段在X和Y方向的增量分别为： △x=xe- xs ，△y= ye- ys 设△t=max(|△x|,|△y|) 取时间步长为1/△t,若当前像素点坐标为(xi, yi),则下一个像素点的坐标可由以下两式确定： 　　　　　　　　　　　　　 　 xi+1=xi+dx=xi+Dx/Dt 　 yi+1=yi+dy=yi+Dy/Dt 2、中点直线生成算法: 假定直线斜率0m1，且已确定点亮象素点P（Xi ,Yi ）,则下一个与直线最接近的像素只能是P1点或P2点。设M为中点，Q为交点，现需确定下一个点亮的象素。 当M在Q的下方- P2离直线更近更近-取P2 。 M在Q的上方- P1离直线更近更近-取P1 M与Q重合， P1、P2任取一点。 设直线方程为：ax+by+c=0 其中a=y0-y1, b=x1-x0, c=x0y1-x1y0 F(x,y)=0 点在直线上面 F(x,y)0 点在直线上方 F(x,y)0 点在直线下方 欲判断M点是在Q点上方还是在Q点下方，只需把M代入F(x,y),并检查它的符号。 令d=F(M)=F(xp+1,yp+0.5) =a(xp+1)+b(yp+0.5)+c 当d0，M在直线下方，取右上方P2； 当d0，M在直线上方，取右方P1； 当d=0，选P1或P2均可，约定取P1； 若d0-M在直线上方-取P1; 此时再下一个象素的判别式为 d1=F(xp+2, yp+0.5)=a(xp+2)+b(yp+0.5)+c = a(xp +1)+b(yp +0.5)+c +a =d+a； 增量为a 若d0-M在直线下方-取P2; 此时再下一个象素的判别式为 d2= F(xp+2, yp+1.5)=a(xp+2)+b(yp+1.5)+c = a(xp +1)+b(yp +0.5)+c +a +b =d+a+b ； 增量为a＋b 画线从(x0, y0)开始，d的初值 d0=F(x0+1, y0+0.5)= a(x0 +1)+b(y0 +0.5)+c = F(x0, y0)+a+0.5b = a+0.5b 由于只用d 的符号作判断，为了只包含整数运算, 可以用2d代替d来摆脱小数，提高效率。 同时d每次的增量也已相应乘2操作。 3、Bresenham画线算法 基本思想： 根据直线的斜率确定选择X或者Y方向作为计长方向, 在此方向上每次递增一个单位步长(或者一个像素单位), 另一个方向上是否同时产生一个单位增量由一个计算量很小的判别式来判断。 算法原理： 以斜率在0~1之间的直线段为例，在第一象限内 这种情况下, 选择X方向为计长方向，即增量dx=1。，设Pi (xi,yi)是已选定的离直线最近的像素, 现在要决定Pi+1是D还是C 显然, 若d0.5, 则C靠近直线, 应选C； 若d=0.5, 则应选D。(m= △y/△x) 令e=d-0.5(初值为m-0.5), 即有： e 0时, 选Pi+1(xi+1, yi), 更新e=e+m; e =0时, 选Pi