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