《计算机图形学资料》由会员分享,可在线阅读,更多相关《计算机图形学资料(38页珍藏版)》请在金锄头文库上搜索。
1、文档供参考,可复制、编制,期待您的好评与关注! 基本图形的生成1.1 直 线数学上,理想的直线是由无数个点构成的集合,没有宽度。计算机绘制直线是在显示器所给定的有限个像素组成的矩阵中,确定最佳逼近该直线的一组像素,并且按扫描线顺序,对这些像素进行写操作,实现显示器绘制直线,即通常所说的直线的扫描转换,或称直线光栅化。由于一图形中可能包含成千上万条直线,所以要求绘制直线的算法应尽可能地快。本节介绍一个像素宽直线的常用算法:数值微分法(DDA)、中点画线法、Bresenham 算法。1.1.1 DDA(数值微分)算法DDA算法原理:图1-1 DDA方法扫描转换连接两点图1-2 中点画线法每步迭代涉
2、及的像素和中点示意图如图1-1所示,已知过端点的直线段;直线斜率为,从的左端点开始,向右端点步进画线,步长=1(个像素),计算相应的坐标;取像素点, round(y)作为当前点的坐标。计算,当,即当x每递增1,y递增k(即直线斜率)。注意:上述分析的算法仅适用于k1的情形。在这种情况下,x每增加1, y最多增加1。当时,必须把x,y地位互换,y每增加1,x相应增加1/k(请参阅后面的Visual C+程序)。1.1.2 生成直线的中点画线法中点画线法的基本原理如图1-2所示。在画直线段的过程中,当前像素点为P,下一个像素点有两种选择,点P1或P2。M为P1与P2中点,Q为理想直线与X=Xp+1
3、垂线的交点。当M在Q的下方时,则P2应为下一个像素点;当M在Q的上方时,应取P1为下一点。中点画线法的实现:令直线段为Lp0(x0,y0), p1(x1, y1),其方程式F(x, y)=ax+by+c=0。其中,a=y0y1, b=x1x0, c=x0y1x1y0;点与L的关系如下。在直线上,F(x, y)=0;在直线上方,F(x, y)0;在直线下方,F(x, y)0。把M代入F(x, y),判断F的符号,可知Q点在中点M的上方还是下方。为此构造判别式d=F(M)=F(xp+1, yp+0.5)=a(xp+1)+b(yp+0.5)+c。当d0,L(Q点)在M下方,取P1为下一个像素。当d=
4、0,选P1或P2均可,取P1为下一个像素。其中d是xp, yp的线性函数。1.1.3 Bresenham算法Bresenham算法是计算机图形学领域使用最广泛的直线扫描转换算法。由误差项符号决定下一个像素取右边点还是右上方点。图1-3 第一象限直线光栅化Bresenham算法设直线从起点(x1, y1)到终点(x2, y2)。直线可表示为方程y=mx+b,其中b=y1mx1, m=(y2y1)/(x2x1)=dy/dx;此处的讨论直线方向限于第一象限,如图1-3所示,当直线光栅化时,x每次都增加1个单元,设x像素为(xi,yi)。下一个像素的列坐标为xi+1,行坐标为yi或者递增1为yi+1,
5、由y与yi及yi+1的距离d1及d2的大小而定。计算公式为 y = m(xi + 1) + b (1.1) d1 = y yi (1.2) d2=yi+1y (1.3)如果d1d20,则yi+1=yi+1,否则yi+1=yi。式(1.1)、(1.2)、(1.3)代入d1d2,再用dx乘等式两边,并以Pi=(d1d2),dx代入上述等式,得Pi = 2xidy2yidx+2dy+(2b1)dx(1.4)d1d2是用以判断符号的误差。由于在第一象限,dx总大于0,所以Pi仍旧可以用做判断符号的误差。Pi+1为Pi+1 = Pi+2dy2(yi+1yi)dx(1.5)求误差的初值P1,可将x1、y1
6、和b代入式(1.4)中的xi、yi,而得到P1 = 2dydx综述上面的推导,第一象限内的直线Bresenham算法思想如下:(1)画点(x1, y1),dx=x2x1,dy=y2y1,计算误差初值P1=2dydx,i=1。(2)求直线的下一点位置xi+1 = xi + 1,如果Pi0,则yi+1=yi+1,否则yi+1=yi。(3)画点(xi+1, yi+1)。(4)求下一个误差Pi+1,如果Pi0,则Pi+1=Pi+2dy2dx,否则Pi+1=Pi+2dy。(5)i=i+1;如果idx+1则转步骤(2);否则结束操作。1.1.4 程序设计1程序设计功能说明为编程实现上述算法,本程序利用最基
7、本的绘制元素(如点、直线等),绘制图形。如图1-4所示,为程序运行主界面,通过选择菜单及下拉菜单的各功能项分别完成各种对应算法的图形绘制。图1-4 基本图形生成的程序运行界面2创建工程名称为“基本图形的生成”单文档应用程序框架(1)启动VC,选择“文件”|“新建”菜单命令,并在弹出的新建对话框中单击“工程”标签。(2)选择MFC AppWizard(exe),在“工程名称”编辑框中输入 “基本图形的生成”作为工程名称,单击“确定”按钮,出现Step 1对话框。(3)选择“单个文档”选项,单击“下一个”按钮,出现Step 2对话框。(4)接受默认选项,单击“下一个”按钮,在出现的Step 3St
8、ep 5对话框中,接受默认选项,单击“下一个”按钮。(5)在Step 6对话框中单击“完成”按钮,即完成“基本图形的生成”应用程序的所有选项,随后出现工程信息对话框(记录以上步骤各选项选择情况),如图1-5所示,单击“确定”按钮,完成应用程序框架的创建。图1-5 信息程序基本3编辑菜单资源设计如图1-4所示的菜单项。在工作区的ResourceView标签中,单击Menu项左边“+”,然后双击其子项IDR_MAINFRAME,并根据表1-1中的定义编辑菜单资源。此时VC已自动建好程序框架,如图1-5所示。表1-1 菜单资源表菜单标题菜单项标题标示符ID直线DDA算法生成直线ID_DDALINEB
9、resenham算法生成直线ID_BRESENHAMLINE中点算法生成直线ID_MIDPOINTLINE4添加消息处理函数利用ClassWizard(建立类向导)为应用程序添加与菜单项相关的消息处理函数,ClassName栏中选择CMyView,根据表1-2建立如下的消息映射函数,ClassWizard会自动完成有关的函数声明。表1-2 菜单项的消息处理函数菜单项ID消 息消息处理函数ID_DDALINECONMMANOnDdalineID_MIDPOINTLINECONMMANOnMidpointlineID_BRESENHAMLINECONMMANOnBresenhamline5程序结构
10、代码,在CMyView.cpp文件中相应位置添加如下代码:/ DDA算法生成直线void CMyView: OnDdaline()CDC* pDC=GetDC();/获得设备指针int xa=100,ya=300,xb=300,yb=200,c=RGB(255,0,0);/定义直线的两端点,直线颜色 int x,y; float dx, dy, k; dx=(float)(xb-xa), dy=(float)(yb-ya); k=dy/dx, y=ya; if(abs(k)1) for (x=xa;xSetPixel (x,int(y+0.5),c); y=y+k; if(abs(k)=1)
11、for (y=ya;ySetPixel (int(x+0.5),y,c); x=x+1/k; ReleaseDC(pDC);说明:(1)以上代码理论上通过定义直线的两端点,可得到任意端点之间的一直线,但由于一般屏幕坐标采用右手系坐标,屏幕上只有正的x, y值,屏幕坐标与窗口坐标之间转换知识请参考第3章。(2)注意上述程序考虑到当k1的情形x每增加1,y最多增加1;当k1时,y每增加1,x相应增加1/k。在这个算法中,y与k用浮点数表示,而且每一步都要对y进行四舍五入后取整。/中点算法生成直线void CMyView:OnMidpointline() CDC* pDC=GetDC();int x
12、a=300, ya=200, xb=450, yb=300,c=RGB(0,255,0); float a, b, d1, d2, d, x, y; a=ya-yb, b=xb-xa, d=2*a+b; d1=2*a, d2=2* (a+b); x=xa, y=ya; pDC-SetPixel(x, y, c); while (xxb) if (dSetPixel(x, y, c); ReleaseDC(pDC);说明:(1)其中d是xp, yp的线性函数。为了提高运算效率,程序中采用增量计算。具体算法如下:若当前像素处于d0情况,则取正右方像素P1(xp+1, yp),判断下一个像素点的位置,应计算d1=F(xp+2, yp+0.5)=a(xp+2)+b(yp+0.5)=d+a;,其中增量为a。若d=0) s1=1; else s1=-1;if(y2-y1=0) s2=1; else s2=-1;if(deltaydeltax)temp=deltax;deltax=deltay;deltay=temp;interchange=1; else interchange=0; f=2*deltay-deltax; pDC-SetPix