Bresenham算法的实现思路

提示：扫一扫查出行
【扫一扫了解最新限行尾号】
复制提示

条件已知两个点的坐标p1(x0,y0)，p2(x1,y1)要求画出这条直线

之后的e代表每次的误差积累，初始值为0，可以计算出斜率为k=dy/dx=(y0-y1)/(x0-x1)

1、x为阶跃步长（直线光栅化）适用于0<k<1的情况

即x每次增加1，但是y的坐标根据其是靠近该点所处的单元格的距离来决定，如果离上边近则y加1，如果离下边近则还是y

可以知道机器在画每一个点的时候都会有误差，则画出的第一个点的坐标（x0,y0）也相当于是（x0,y0+e），那么可以知道下个点的坐标为红色的p(x0+1,y0+e+k),但是机器在画p点时，如果p点的纵坐标不是整数时，应该进行判断这个点是接近黄点(x0+1,y0+1)还是接近蓝点(x0+1,y0)，最后如果接近黄点则e累积k-1,如果接近蓝点则e累积k

在判断条件的两边同乘2得：

那么程序为：

void LineB2(int x0,int y0,int x1,int y1,int color,CDC *p){
int i,x,y,dx,dy;
float k,e;
dx=x1-x0;
dy=y1-y0;
k=(float)dy/dx;
e=-0.5;x=x0;y=y0;
for(i=;i<=dx;i++){
p->SetPixel(x,y,color);
x=x+;
e=e+k;
if(e>=){
y+=;
e=e-;
}
}
}

上面的程序中e使用了多次除法运算，算法效率低。

优化

将上面的k使用dy/dx替换，并且同乘dx后，令ξ=dx*e则：

则程序为：

//在程序中e为上面的ξ
void LineB3(int x0,int y0,int x1,int y1,int color,CDC *p){
int i,x,y,dx,dy;
float e=;
dx=x1-x0;
dy=y1-y0;
e=-0.5;x=x0;y=y0;
for(i=;i<=dx;i++){
p->SetPixel(x,y,color);
x=x+;
if(*(e+dy)>=dx){
y+=;
e=e+dy-dx;
}else{
e=e+dy;
}
}
}

2、y为阶跃步长（直线光栅化）适用于k>1的情况

同上可以推出：

可以知道机器在画每一个点的时候都会有误差，则画出的第一个点的坐标（x0,y0）也相当于是（x0+e,y0），那么可以知道下个点的坐标为红色的p(x0+e+1/k,y0+1),但是机器在画p点时，如果p点的横坐标不是整数时，应该进行判断这个点是接近右边(x0+1,y0+1)还是接近蓝左边(x0,y0+1)，最后如果接近右边则e累积1、k-1,如果接近左边则e累积1/k