数值微分算法(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;
    }
}

任意斜率代码实现如下：

void DDALine(int x0,int y0,int x1,int y1,int color) { int dx,dy,n,k;float xinc,yinc,x,y; dx=x1-x0;dy=y1-y0; if (abs(dx)>abs(dy)) n=abs(dx); else n=abs(dy); xinc=(float)dx/n; yinc=(float)dy/n; x=(float)x0;y=(float)y0; for (k=1;k<=n;k++) { Drawpixel(int (x+0.5),int (y+0.5),color); x+=xinc; y+=yinc; } }