DDA算法和Bresenham算法画直线

DDA算法:

步骤:

1、求斜率K

2、x坐标增加

3、y坐标增加K,并进行+0.5取整

4、重复步骤2、3直到绘制完

例如:

Bresenham算法:

步骤:

1、计算Δx Δy

2、e=-Δx x0=x0 y=y0

3、绘制点(x,y)

4、更新e=e+2*Δy,且

1
2
3
4
5
1、e>0 时:
(x,y)变为(x+1,y+1
e变为e-2x
2、e<=0时:
x,y)变为(x+1,y)

5、重复步骤3、4直到绘制完

MATLAB代码:

DDA算法:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
function DDA(x,y)
L=max(abs(x(1)-y(1)),abs(x(2)-y(2)));
k=(x(2)-y(2))/(x(1)-y(1));
X=[];Y=[];
for i=1:L+1
X(i)=x(1);
Y(i)=x(2);
x(1)=x(1)+1;
x(2)=x(2)+k;
end
Y=fix(Y+0.5);
figure,
plot(X,Y);
end

Bresenham算法:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
function Bresenham(x,y)
detax=abs(x(1)-y(1));
detay=abs(x(2)-y(2));
e=-detax;
X=[];Y=[];
L=max(abs(x(1)-y(1)),abs(x(2)-y(2)));
for i=1:L+1
X(i)=x(1);
Y(i)=x(2);
e=e+2*detay;
if(e>0)
x(1)=x(1)+1;
x(2)=x(2)+1;
e=e-2*detax;
else
x(1)=x(1)+1;
x(2)=x(2);
end
end
figure,
plot(X,Y);
end

演示:

1
2
3
4
5
6
7
8
9
10
11
12
clc;
clear;
close all;
x=zeros(1,2);
string=['Enter starting point Coordinates:'];
x=input(string);
y=zeros(1,2);
string=['Enter the coordinates of the end point:'];
y=input(string);
Bresenham(x,y);
DDA(x,y);
Midpoint(x,y);

(0,0)——>(5,3)

感觉和直线差了不少


(0,0)——>(100,40)

有点直线那意思了


(0,0)——>(1000,1000)

已经看不出来是折线

如果有错误欢迎指正。