计算机图形学课程论文

计算机图形学期末综合理论

摘要：计算机图形学直线变换算法二维图形算法图形填充算法

引言：计算机图形学是一种使用数学算法将二维或三维图形转化为计算机显示器的栅格形式的科学。简单地说，计算机图形学的主要研究内容就是研究如何在计算机中表示图形、以及利用计算机进行图形的计算、处理和显示的相关原理与算法。虽然通常认为CG是指三维图形的处理，事实上也包括了二维图形及图像的处理。狭义地理解，计算机图形学是数字图象处理或计算机视觉的逆过程：计算机图形学是用计算机来画图像的学科，数字图象处理是把外界获得的图象用计算机进行处理的学科，计算机视觉是根据获取的图像来理解和识别其中的物体的三维信息及其他信息。实际上，计算机图形学、数字图象处理和计算机视觉在很多地方的区别不是非常清晰，很多概念是相通的，而且随着研究的深入，这些学科方向不断的交叉融入，形成一个更大的学科方向，可称之为“可视计算”。计算机图形学主要包含四大部分的内容：建模(Modeling)、渲染(Rendering)、动画(Animation)和人机交互(Human–computer Interaction, HCI)。

目录

1、概论 (4)

1.1什么是计算机图形学 (4)

1.2计算机图形学的研究内容 (4)

1.3什么是计算机图形学 (4)

1.4计算机中表示图形的方法 (5)

2、直线生成算法 (4)

2.1直线数值微分算法 (5)

2.1.1数值微分算法基本原理 (6)

2.1.2数值微分算法基本原理 (6)

2.1.3数值微分算法设计与代码实现 (7)

2.1.3数值微分算法设计与代码实现 (8)

2.1.4小结 (9)

2.2直线中点画线算法 (9)

2.2.1直线中点画线算法基本原理 (9)

2.2.2数值微分算法描述与步骤 (10)

2.2.3中点划线算法设计与代码实现 (12)

2.2.4小结 (10)

2.3直线Breseham画线算法 (12)

2.3.1直线Breseham画线基本原理 (13)

2.3.2直线Breseham画线算法描述与步骤 (13)

2.3.3直线Breseham画线算法设计与代码实现 (17)

2.3.4小结 (18)

3、二维图形变换 (19)

3.1二维图形平移变换 (19)

3.1.1二维图形平移变换基本原理 (19)

3.1.2二维图形平移变换算法描述与步骤 (19)

3.1.3二维图形平移变换算设计与代码实现 (21)

3.2二维图形缩放变换 (21)

3.2.1二维图形缩放变换基本原理 (21)

3.2.2二维图形缩放变换算法描述与步骤 (21)

3.2.3二维图形缩放变换算法设计与代码实现 (22)

3.3二维图形对换变换 (23)

3.3.1二维图形对换变换基本原理 (23)

3.3.2二维图形对换变换算法描述与步骤 (23)

3.3.3二维图形对换变换算设计与代码实现 (25)

3.4二维图形旋转变换 (26)

3.4.1二维图形旋转变换基本原理 (26)

3.4.2二维图形旋转变换算法描述与步骤 (26)

3.4.3二维图形旋转变换算设计与代码实现 (26)

4、图形填充算法 (27)

4.1种子填充算法 (30)

4.1.1种子填充算法基本步骤 (30)

4.1.2种子填充算法设计与代码实现 (30)

4.2边标志填充算法 (30)

4.2.1边标志填充算法基本步骤 (31)

4.2.2种子填充算法设计与代码实现 (31)

4.3小结 (32)

5、总结和展望 (32)

6、参考文献 (33)

1、概论

1.1什么是计算机图形学

计算机图形学是利用计算机研究图形的表示、生成、处理、显示的学科。计算机图形学计算机科学中，最为活跃、得到广泛应用的分支之一。

1.2计算机图形学的研究内容

如何在计算机中表示图形、以及利用计算机进行图形的计算、处理和显示的相关原理与算法，构成了计算机图形学的主要研究内容。

图形硬件、图形标准、图形交互技术、光栅图形生成算法、曲线曲面造型、实体造型、真实感图形计算与显示算法，以及科学计算可视化、计算机动画、自然景物仿真、虚拟现实等。

1.3基本概念

1.3.1怎样用计算机生成、处理和显示图形的学科

（1）图形的输入

如何开发利用图形输入设备及软件将图形输入到计算机中去，以便作各种处理。

（2）图形的处理

包括对图形进行变换(几何变换，投影变换)和运算（集合运算）,裁减，着色，形变等……

（3）图形的输出

如何将图形特定的表示形式转换成图形输出系统便于接受的表示形式，并将图形在显示屏或打印机等输出设备上输出，或以文件的形式保存在磁盘上。

1.3.2相关的概念

（1）图形：计算机图形学的研究对象，能在人的视觉系统中产生视觉印象的客观对象，包括自然景物、拍摄到的图片、用数学方法描述的图形等等

－构成图形的要素

（2）几何要素：刻画对象的轮廓、形状等

（3）非几何要素：刻画对象的颜色、材质等

1.4计算机中表示图形的方法

1.4.1点阵表示

?枚举出图形中所有的点(强调图形由点构成)

?简称为图像（数字图像）

1.4.2参数表示

?由图形的形状参数(方程或分析表达式的系数，线段的端点坐标等)以及属

性参数(颜色、线型等)来表示图形

?简称为图形

2.直线生成算法

2.1直线数值微分算法

2.1.1数值微分算法基本原理

?设直线的起点坐标为(x0，y0)，终点坐标为(x1, y1)，二者均为整数坐标值，

令dx＝x1- x0, dy＝y1- y0，则直线斜率k为：

k=

dx

dy

目标是能快速地求出能很好地表示直线的像素,涉及到大量的重复性的运算，因此，尽量用加法代替乘法运算，因为计算机中：

加法运算比乘法运算效率高

2.1.2数值微分算法描述与步骤

当|k|≤1时，从x的左端点x0开始，向x右端点步进。步长=1个像素(即Dx =1)，计算相应的y坐标y=kx+B；取像素点(x, round(y))作为当前点的坐标。

注意到：

yi+1 = kxi+1+B

= k(xi+Dx) +B

= kxi+B+kDx

= yi+kDx

= yi+k

即：当x每递增1，y递增k(即直线斜率)；

上述分析的算法仅适用于|k| ≤1的情形。在这种情况下，x每增加1，y最多增加1

当|k| >1时，必须把x，y地位互换

2.1.3数值微分算法设计与代码实现

线段DDA算法伪代码描述下面用伪代码给出DDA算法。Procedure DDA-line(x1, y1, x2, y2)

BEGIN

/求线段在两坐标轴方向改变量的较大者/

IF abs(x2-x1)>=abs(y2-y1) THEN

length=abs(x2-x1)

ELSE

length=abs(y2-y1)

ENDIF

/定义dx或dy中的较大值为1/

dx=(x2-x1)/length dy=(y2-y1)/length x=x1+0.5*Sign(dx) y=y1+0.5*Sign(dy) i=1

WHILE (i

PLOT(Integer(x), Integer(y))

x=x+dx

y=y+dy

i=i+l

END WHILE END

2.1.4小结

这种方法计算方法的缺点是计算量大。考虑到在计算机上实时编辑修改图形时, 连续不断的需要在显示器上显示出大量的直线段, 画线的速度就会非常地慢。因此需要对上述结果详细分析已给出较快的算法。

2.2直线中点画线算法

2.2.1数值微分算法基本原理

?原理

●假设直线的斜率0

●当前像素点为P(x p, y p) ，下一个像素为P1或P2

●设M=(x p+1, y p+0.5)，为p1与p2之中点，Q为理想直线与x=x p+1垂

线的交点。将Q与M的y坐标进行比较：

●M在Q的下方，则P2为下一点；

●M在Q的上方，则P1为下一点。

图2.1a 图2.1b

2.2.2数值微分算法描述与步骤

?构造直线方程

F(x,y) = ax + by +c = 0

?调整方程参数，使得：

●点(x,y)位于直线上方：F(x,y)>0;

●点(x,y)位于直线下方：F(x,y)<0;

?把中点M的坐标带入方程，构造判别式：

d = F(M)

= F(x

p +1,y

p

+0.5)

= a(x

p +1) + b(y

p

+0.5) + c

?因此：

●当d>0时，点M位于Q点上方，取P1为下一点;

●当d<0时，点M位于Q点下方，取P2为下一点；

●d=0时，Q和M重叠，任取P1、P2。如果取P1，则可合并为d≥0

时取P

1

。

2.2a

2.2b

?判别式d的增量算法：

●设确定P i(x i,y i)位置时判别式的值为d i,则：

●如果d i≥0，则x i = x i-1+1，y i = y i-1

d i = F(x i-1+1,y i-1+0.5) = F(x i,y i+0.5)

= ax i + b(y i+0.5) + c

d i+1 = F(x i+1,y i+0.5)

= a(x i+1) + b(y i+0.5) + c

= ax i + b(y i+0.5) + c + a

=d i + a

di≥0→di+1= di+a

◆如果d i<0，则x i = x i-1+1，y i = y i-1+1

d i = F(x i-1+1,y i-1+0.5) = F(x i,y i-0.5)

= ax i + b(y i-0.5) + c

d i+1 = F(x i+1,y i+0.5)

= a(x i+1) + b(y i+0.5) + c

= ax i + b(y i-0.5) + c + a + b

=d i + (a + b)

d i<0→d i+1= d i+a+b

判别式d的初值：

直线的起始端点坐标为(x0, y0)，则：

d1 = F(x0+1,y0+0.5)

= ax0 + by0 + c + a + 0.5b

=F(x0,y0) + a + 0.5b (F(x0,y0) = 0)

= a + 0.5b （存在0.5b，为浮点运算）

整数算法：

算法中仅用到判别式d的符号，所以可用2d代替d进行递增运算：

2.2.3中点划线算法设计与代码实现

void Midpoint Line (int x0,int y0,int x1, int y1,int color)

{

int a, b, d1, d2, d, x, y;

a=y0-y1, b=x1-x0, d=2*a+b;

d1=2*a, d2=2* (a+b);

x=x0, y=y0;

drawpixel(x, y, color);

while (x

{

if (d<0)

{x++, y++, d+=d2; }

else

{x++, d+=d1;}

drawpixel (x, y, color);

}

}

2.2.4小结

如上面代码所示，DDA算法中的y和k都必须用浮点数表示，并且每一步运算都要对y进行舍入取整，这不利于硬件实现。中点画线法只包含整数变量，并且不含乘除法，因而解决了上述问题

2.3直线Breseham画线算法

2.3.1直线Breseham画线基本原理

Bresenham直线算法最初是为数字绘图仪而设计的。它的目标是选择表示直线的最佳像素点位臵. 为此,算法根据直线的斜率确定在x或y方向上每次递增一个单位,而另一方向上根据理论直线段与最近像素点的距离每次的增量为0或1。我们首先讨论直线斜率m ∈ [0, 1]且x2 > x1时的整数Bresenham算法, 然后再推广到画任意线段的算法。

2.3.2直线Breseham画线算法描述与步骤

当直线斜率m ∈ [0, 1],且x2 > x1时,根据式(2.2.9)有:

又由于显示直线的像素点只能取整数值坐标,可以假设直线上第i个像素点坐标为(x i, y i),它是直线上点(x i, y i)的最近似并且x i = x i。于是,要表达的直线上下一个(x i+1, y i+1) = (x i + 1, y i + m)的最近似的像素点的可能位置是(x i +1,y i)或(x i +1,y i +1)。

在x = x i +1处,直线上点的坐标y = m(x i +1)+b. 该点与点(x i +1,y i) 和(x i + 1, y i + 1)的距离分别是d1和d2:

d1 =y?y i =m(x i +1)+b?y i?

d2 =(y i +1)?y=(y i +1)?m(x i +1)?b

这两个距离的差为:

d1 ?d2 =2m(x i +1)?2y i +2b?1

这个差有如下几何意义:?

(1)当此值为正时,真正的直线上点离像素点(x i + 1, y i + 1)较近,说明下一个直线上的像素点应取(x i + 1, y i + 1)。

(2)当此值为负时,真正的直线上离像素点(x i + 1, y i)较近,说明下一个直线

像素点应取(x i + 1, y i)。

(3)当此值为零时,真正的直线上离上、下两个像素点的距离相等,我们规

定取(x i + 1, y i)作为下一个直线像素点。?因此,只要利用(d1 ? d2)的符号就可以决定下一个像素点的选择。如果我们

定义一个新的判别式:?p i =?x?(d1 ?d2)=2?y·x i ?2?x·y i +c (2.2.18)

因此式中的?x = (x2 ? x1) > 0,p i与(d1 ? d2)有相同符号; ?y = y2 ? y1;常数c =

2?y + ?x(2b ? 1)。p i的一个优点是省去了(d1 ? d2)中为了计算m所需要的除法运算。我们知道除法运算用硬件实现是比较复杂的。

现在我们要进一步化简上述误差判别式以得到递推公式,消除常数c. 以i + 1代换此式中的i,得到:

p i+1 =2?y·x i+1 ?2?x·y i+1 +c 与前式相减, 并利用x i+1 = x i + 1,可得,

p i+1 = p i + 2?y ? 2?x · (y i+1 ? y i) 再假设直线的初始端点恰好是其像素点的坐标,即满足:

于是可得p i的初值

y1 = mx1 + b p1 = 2?y ? ?x

这样,利用误差判别变量,并注意到每一步x的增量为x i+1 ? x i = 1就可得到如下算法表示:

(1).

(2).如果p i ≥ 0,

p1 =2?y??x x i+1 = x i + 1 (2.2.23)

y i+1 = y i + 1?p i+1 =p i +2(?y??x) (2.2.24)

2.2. 直线点阵转换算法(3).如果p i < 0,

y i+1 = y i?p i+1 = p i + 2?y

从式(2.2.20)可以看出,第i + 1步的判别变量p i+1仅与第i步的判别变量p i、直线的两个端点坐标分增量?x和?y有关,计算也很简单,只用整数相加和乘2运算,没有四舍五入处理,而乘2可利用左移一位来实现,因此这个算法速度快并易

于硬件实现。

该算法的主要步骤如下:?(1)输入线段的两个端点分别存于(x s, y s)和(x e, y e)中, (2)将第一点作为起始点, 即有(x1, y1) = (x s, y s), (3)分别计算?x、?y及p1,若p1 < 0,下一点为(x1 + 1, y1),否则,取(x1 +

1,y1 +1),?(4)以单位步长增加x坐标,按式((2.2.24))或((2.2.25))计算p i。若p i < 0,下

一点的y坐标不变,否则y坐标加1。?(5)重复步骤(4)直到x逐步增加到x e为此。?在前面的处理中, 我们假设m ∈ [0, 1]。这一假设的几何意义是要绘的直线段与x?轴的夹角不超过与y?轴的夹角。或者更直观地, 直线段的方向更靠近x?轴方向,如果有m > 1, 则说明直线段的方向更靠近y?轴方向。如果我们互换两个坐标轴, 则直线段就满足前面的假设条件. 因此, 只要将算法中的x和y 对换,则上述两个公式依然有效。如果m < 0,则相应的只是x或y的

方向的改变后, 上述两个公式依然有效。这时实际上只是改变?x或?y的符号即可.

2.3.3直线Breseham画线算法设计与代码实现

Procedure Bresenham-line(xs, ys, xe, ye)

BEGIN

dx=ABS(xe-xs);dy=ABS(ye-ys);x=xs;y=ys

s1=SIGN(xe-xs), s2=SIGN(ye-ys)

If dy>dx THEN

temp=dx;dx=dy:dy=temp

interchange=1

ELSE

interchange=O

ENDIF

p=2*dy-dx

FOR i=1 TO dx

PLOT(x, y)

IF p>=O THEN

IF interchange=1

x=x+s1

ELSE y=y+s2

ENDIF

p=p-2*dx

ENDIF

IF interchange=1 THEN

y=y+s2 ELSE

x=x+s1

ENDIF

p=p+2*dy

NEXT i

END

2.3.4小结

以下是近似图形和真实直线之间的差距，可以看得出来，其实实际上图形差距已经很小了。