图形学笔记1——直线段扫描转换算法

原创 2007年10月11日 02:50:00
 

DDA

已知过端点P0x0, y0, P1 (x1, y1)的直线段; 直线斜率k = (y1-y0) / (x1-x0) .

画线过程为: x的左端点x0开始,向x右端点 步进,步长=1(像素),按y=k x + b 计算相应的y坐标,并去像素点(x, round (y) )作为当前点的坐标。

可以提高效率:

设步长为Δx x i+1 = xi +Δx, 于是: yi+1=kxi+1 +b = kxi + kΔx +b= yi +kΔx

当Δx=1时, yi+1 = yi +k 这样计算每一点,只需做一次加法运算。

note:

round (y) = int (y+0.5) K>1时,必须把x, y的地位互换:y增加1x相应增加1/k

缺点:由于yk须用浮点数表示,每一步都要对y进行四舍五入后取整,使得该算法不利于硬件实现。

 

改进 

目标:进一步将一个加法改为一个整数法加。

 

中点画线法:

通过观察发现,在画线段的过程中,当前像素点为(xp, yp),下一个像素点有两种选择P1xp+1, yp)和P2xp+1, yp+1)。若M=xp+1, yp+0.5)为P1P2之间中点。Q为理想直线与x = xp +1垂线的交点。 则当MQ的下方时,P2为下一个像素点;当MQ的上方时,应取P1为下一点。(如下图)

http://p.blog.csdn.net/images/p_blog_csdn_net/touzani/303255/o_middle.JPG

对直线段Lp0x0, y0, p1(x1, y1) ),采用方程

F (x, y ) = 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

要判断MQ点上方还是下方,只要把M点坐标代入F(x, y),并判断符号即可。

即判断 d = F(M) = Fxp+1, yp+0.5= a (xp+1)+ b (yp+0.5) + c 的符号

dxp, yp的线性函数,为提高计算效率,可以采用增量计算。

l       d>=0,则取正右方象素P1 (xp+1, yp), 要判断下一个象素,应计算

     d1= F(xp +2, yp +0.5)=a (xp+2)+ b (yp+0.5) + c = d + a 增量为a

l       d<0,则取右上方象素P2 (xp+1, yp +1)。要判断再下一象素,应计算

     d2= F(xp +2, yp +1.5)=a (xp+2)+b (yp+1.5) + c= d + a + b ;增量为ab

设从点(x0, y0)开始画线,d的初值d0 = F (x0+1, y+0.5) = F ( x0, y0 ) + a + 0.5b.

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

进一步,由于d的增量都是整数,只有初始值包含小数,可以用2d代替d来摆脱浮点运算。 这是不影响符号判断的。

 

 

Bresenham算法

Bresenham算法类似于中点法,由误差项符号决定下一个像素取右边点还是右上点。

算法原理:

设直线方程为 y = kx +b, 则有 yi+1= yi + k (xi+1 – xi) = yi + k .

x列的坐标为(xi, yi, 那么下一个像素的列坐标为xi+1,

而行坐标要么仍为yi,要么递增1yi+1

http://p.blog.csdn.net/images/p_blog_csdn_net/touzani/303255/o_bresenham.JPG

是否增1取决于误差项d的值(如上图)

d的初值为0.

x下标每增加1d的值相应增加k.

d0.5时,直线与xi+1列的交点更接近于

像素(xi +1, yi +1),即该像素在y方向递增1.

同时作为下一次计算的基点,因此d值减1.

d<0.5时,交点更接近于像素(xi +1, yi)。

为方便计算,令e = d - 0.5, e的初值为-0.5. 增量为k.

void Bresenhamline (int x0,int y0, int x1, int y1, int color)

/**//* 起点(x0, y0) 终点(x1, y1)*/

...{   int x, y, dx, dy;

    
float k, e;

    dx 
= x1-x0, dy = y1- y0, k=dy/dx;

    e
=-0.5, x=x0, y=y0;

    
for (i=0; i <= dx; i++)

    
...{    drawpixel (x, y, color);

          x
=x+1,e=e+k;

          
if (e³0)

          
...{ y++, e=e-1;}

    }


}

 

为避免小数与除法,可以改用整数: 做如下替换 e1 = 2 * e * dx

Bresenham画线法程序

void IntegerBresenhamline (int x0,int y0, int x1, int y1, int color)

/**//* 起点(x0, y0) 终点(x1, y1)*/

...{   int x, y, dx, dy, e;

dx 
= x1-x0,  dy = y1- y0,  e = -dx;

    x
=x0, y=y0;

    
for (i=0; i <= dx; i++)  ...{  

 drawpixel (x, y, color);

          x
=x++,e=e+2*dy;

          
if (e >= 0 )  ...{ y++, e = e - 2*dx; }

    }


}

 

直线的扫描转换算法

这是第一篇博客,顺便深入学习一下markdown语法。 之后的图形学文章都使用C++实现。 目录 目录 数值微分法 算法理论 算法实现 中点画线法 算法理论 算法实现 Br...
  • chenge0925
  • chenge0925
  • 2018-01-31 19:37:58
  • 52

直线段扫描转换算法-数值微分法(DDA)

在数学上,直线上的点有无穷多个。担当在计算机光栅显示器屏幕上表示这条直线时需要做一些处理。 为了在光栅显示器上用这些离散的像素点逼近这条直线,需要知道这些像素点的x,y坐标。 求出过p0,p1的直...
  • lyp_1020k
  • lyp_1020k
  • 2017-12-02 23:05:21
  • 295

c++ 的直线扫描转换算法源程序

  • 2009年10月28日 13:04
  • 7.37MB
  • 下载

计算机图形学基础 : 基本图形生成算法之直线的扫描转换

数值微分法(DDA)、中点画线法和Bresenham画线算法的学习.
  • JY_95
  • JY_95
  • 2016-08-14 01:35:07
  • 2078

扫描转换直线段--DDA算法

Void LineDDA(int x0,int y0,int x1,int y1,int color) //假定x0 { int x;   float dx,dy,a,y;   dx=x1-x...
  • u012358638
  • u012358638
  • 2013-10-11 00:33:55
  • 776

直线段的扫描转换算法

  • 2011年03月18日 22:31
  • 271KB
  • 下载

图形学扫描转换直线段

  • 2009年05月06日 10:25
  • 1.86MB
  • 下载

动态扫描直线vc++

  • 2011年03月18日 22:43
  • 426B
  • 下载

计算机图形学——直线的三种扫描转换算法

计算机图形学是算比较抽象的一门课程吧,而且内容也比较枯燥,如果没有比较好的耐心,一时半会儿是看不见计算机图形学究竟有什么作用的,但是其中有些内容呢跟C语言有关联,比如直线的扫描转换算法。大部分时候C语...
  • philip2345
  • philip2345
  • 2016-11-30 23:27:03
  • 3086

直线生成算法的实现:分别利用DDA算法、中点Bresenham算法和改进的Bresenham算法扫描转换直线段P1P2

直线生成算法的实现:分别利用DDA算法、中点Bresenham算法和改进的Bresenham算法扫描转换直线段P1P2,其中P1为(0, 0), P2为(8, 6)。   // fhk.cpp :...
  • tingary
  • tingary
  • 2017-05-18 16:20:50
  • 1292
收藏助手
不良信息举报
您举报文章:图形学笔记1——直线段扫描转换算法
举报原因:
原因补充:

(最多只允许输入30个字)