数值微分算法(DDA画线算法)

基本思想

    已知过端点 P₀ (x₀   ,y₀ ),P₁ (x₁ ,y₁ ) 的直线段L：y=kx+b
    直线斜率为 k=(y1-y0)/(x1-x0)
    从x的左端点x0开始，向x右端点步进。步长=1(个象素)，计算相应的y坐标 y=kx+b；取象素点(x, round(y))作为当前点的坐标。

 

效果如下：

当|k|<=1时，实现代码：

void DDALine( int x0, int y0, int x1, int y1, int color) 
{ 
     int x; float dx,dy,y,k; 
    dx=x1-x0,dy=y1-y0; 
    k=dy/dx,y=y0; 
     for (x=x0;x<=x1;x++) 
    { 
        Drawpixel(x, int (y+0.5),color); 
        y=y+k; 
    } 
}

当|k|>1时，要将X和Y地位互换，实现代码：

void DDALine( int x0, int y0, int x1, int y1, int color) 
{ 
     int y; float dx,dy,x,k; 
    dx=x1-x0,dy=y1-y0; 
    k=dx/dy,x=x0; 
     for (y=y0;y<=y1;y++) 
    { 
        Drawpixel( int (x+0.5),y,color); 
        x=x+k; 
    } 
}

任意斜率代码实现如下：

[cpp]  view plain copy

1. void DDALine(int x0,int y0,int x1,int y1,int color)  
2. {  
3.     int dx,dy,n,k;float xinc,yinc,x,y;  
4.     dx=x1-x0;dy=y1-y0;  
5.     if (abs(dx)>abs(dy))  
6.         n=abs(dx);  
7.     else  
8.         n=abs(dy);  
9.     xinc=(float)dx/n;   yinc=(float)dy/n;  
10.     x=(float)x0;y=(float)y0;  
11.     for (k=1;k<=n;k++)  
12.     {  
13.         Drawpixel(int (x+0.5),int (y+0.5),color);  
14.         x+=xinc;    y+=yinc;  
15.     }  
16. }  

转自：http://blog.csdn.net/akof1314/article/details/5443257