广义Bresenham直线绘制算法while循环内变量D的作用是什么

对应实现代码如下：

#define sign(x) ((x > 0)? 1 : ((x < 0)? -1: 0))

x = x1; 
y = y1;
dx = abs(x2 - x1); 
dy = abs(y2 - y1);
s1 = sign(x2 - x1); 
s2 = sign(y2 - y1);
swap = 0;

if (dy > dx) {
    temp = dx;
    dx = dy;
    dy = temp;
    swap = 1;
}

D = 2*dy - dx;
for (i = 0; i < dx; i++) {
    display_pixel (x, y); 
    while (D >= 0) { 
        D = D - 2*dx;
        if (swap)
            x += s1;
        else
            y += s2;
    } 
    D = D + 2*dy;
    if (swap)
        y += s2;
    else
        x += s1; 
} 

变量D的作用 #

D是算法的整数误差累积项，全程不涉及浮点运算，核心功能是跟踪当前已绘制像素和理想直线路径的偏差，判断短轴（步进幅度更小的轴）什么时候需要补一次位移。
算法开头的轴交换逻辑，是先统一把位移差更大的轴设为主步进轴：如果y方向差dy大于x方向差dx，就交换二者的值并把swap设为1，保证后续循环里主长轴每次固定走1步，短轴只在误差达到阈值时才步进，所有斜率的直线都能复用同一套计算逻辑。
初始化时D = 2*dy - dx是对误差做整数缩放的结果：原始浮点误差下，主长轴每走1像素，理想直线在短轴方向偏移dy/dx，给所有误差项同乘2*dx消去分母，就能把所有计算转成整数运算，避免浮点精度损耗和性能开销。

while循环的执行逻辑 #

外层for循环的迭代次数等于长轴的总像素数（也就是交换后的dx值），每轮循环长轴必然前进1像素，对应循环末尾的步进逻辑：swap=1时长轴是y，所以y累加方向步长s2；swap=0时长轴是x，x累加方向步长s1。
内层while循环专门处理短轴的步进修正：

• 当D >= 0时，说明累积的误差已经跨过半个像素的阈值，理想直线已经落在下一个短轴像素的覆盖范围内，必须给短轴补1次位移。
• 每给短轴补1步，就给D减去2*dx，相当于重置误差基准——短轴移动1格后，之前累积的偏差会被抵消一整段。
• 这里用while而非单次判断if，是做鲁棒性兜底：正常斜率下（dy <= dx）每轮循环最多进1次while，极端参数或计算溢出场景下，循环会持续补步进直到误差回到阈值以下，不会出现偏差持续累积导致直线走形的问题。
  while循环结束后执行D = D + 2*dy，是每走1步长轴后给误差累加短轴的偏移量：对应缩放前的逻辑，长轴走1格，理想直线在短轴方向的偏移量按比例累加到误差项里，为下一轮的决策做准备。

举个简单例子验证：绘制从(0,0)到(6,2)的缓坡直线，此时dx=6、dy=2、swap=0、方向步长s1=1、s2=1，初始D=4-6=-2。整个循环走6轮，最终绘制的像素点依次是(0,0)、(1,0)、(2,1)、(3,1)、(4,1)、(5,2)，完全贴合直线轨迹，没有多余像素。

内容的提问来源于stack exchange，提问作者MiguelP