03基本图元算法2直线算法.docVIP

第3章 基本图元算法 int x,y; } VPOINT; 已知条件：起点,终点 2 直线算法 2.1 DDA直线算法 一、DDA算法 每个点的计算方法： ①在最大位移方向上，取步长=1； ②计算另一坐标分量； x=x±1, y=y±yincre; 或 y=y±1, x=x±xincre; ③取整作为点坐标。 分析：8种斜线情形、水平线、垂直线 二、算法 void DDALine(VPOINT p1, VPOINT p2, int color) { int delta, i; //delta为x、y方向最大位移的绝对值 float xinc, yinc; float x，y; // 留心①：必须float类型 if(abs(p2.x-p1.x)abs(p2.y-p1.y)) { delta=abs(p2.x-p1.x); if(p2.x-p1.x0) xinc=1; else xinc=-1; // 留心yinc=(float)(p2.y-p1.y)/abs(p2.x-p1.x); // 留心③ } else { delta=abs(p2.y-p1.y); if(p2.y-p1.y0) yinc=1; else yinc=-1; xinc=(float)(p2.x-p1.x)/abs(p2.y-p1.y); } x=p1.x; y=p1.y; for(i=1; i=delta; i++) { putpixel(x+0.5, y+0.5, color); //留心④：四舍五入 x=x+xincre; y=y+yincre; } } 水平线、垂直线的处理？ 缺点：存在实数的除法运算。 测试例程： main() { VPOINT p1={200,200}, p2; int radius=100, angle; for(angle=0; angle360; angle++) { p2.x=p1.x+ radius*cos(angle); p2.y=p1.y+ radius*sin(angle); DDALine(p1, p2, 1); } } 三、增量法 由于图形中的点的相邻性，借助相邻点的坐标计算下一点的坐标，将乘/除法运算转换为加/减法运算。 2.2 Bresenham直线算法(中点画线法) 一、构造判别式 构造判别式,将两个候选像素的中点代入，根据值的正负，确定选择哪个像素点。 取判别式为直线方程： 若，则在直线中； 若，则在直线上方； 若，则在直线下方。 二、计算判别式(考虑情形：且p1.xp2.x) 若，取；若，取。 增量法的体现： 如何取下一个像素点？ 如何计算下一步的值？ 能否用当前的值推导下一步的值？ 的初值: 推导下一步的值： 若，取； 若，取。 三、算法 void MidpointLine1(VPOINT p1, VPOINT p2,int color) { int a，b，d1，d2，d，x，y; a=p1.y-p2.y; b=p2.x-p1.x; d=2*a+b; d1=2*a; d2=2*(a+b); // 技巧：只判断正负 putpixel(p1.x,p1.y,color); for(x=p1.x+1,y=p1.y; x=p2.x; x++) { if(d=0) d=d+d1; else { y++; d=d+d2; } putpixel(x，y，color); } } 优点：全部整型数加减运算。 思考：所有情形的算法 void allMidpointLine(VPOINT p1, VPOINT p2,int color) { int delta_x, delta_y; delta_x=p2.x-p1.x; delta_y=p2.y-p1.y; if(abs(delta_y)=abs(delta_x)) { if(delta_x0 delta_y0) MidpointLine1(p1,p2,color); if(delta_x0 delta_y0) MidpointLine1(p2,p1,color); if(delta_x0 delta_y0) MidpointLine2(p1,p2,color); if(delta_x0 delta_y0) MidpointLine2(p2,p1,color); } else { if(delta_x0 delta_y0) MidpointLine3(p1,p2,colo