计算机图形学（一）DDA画线算法讲解与源代码

很早之前就想写一个计算机图形学系列的讲解，可是只写了2篇，然后就搁置了很长一段时间，现在也算是有时间来继续之前的想法了。

首先介绍一下算法：

源代码的下载：点击打开链接

已知直线过端点P0（x0,y0），P1（x1,y1）的直线段的斜率K=(y1-y0)/(x1-x0)，画线的过程为：从x的左端点x0开始，向x的右端点步进，步长=1（像素），按y=kx+b来计算y坐标，并取像素点（x，round(y)）作为当前坐标点。但这样计算 每一个点需要做一个乘法，一个加法。设步长为∆x，有xi+1=xi+1;于是yi+1=k(xi+1)+b=kxi+1，这样的话就能简化算法，只用加法进行计算。

本文程序画线使用的是橡皮筋，然后可以自定义像素大小，更方便容易看出计算出来的点。而且更能了解这个算法的过程。

其中主要的算法：

 

1. void CDrawlineView::DDAline(int x0, int y0, int x1, int y1) //DDA算法！

2. {

3. int x;

4. float dx,dy,y,k;

5. dx=float(x1-x0);

6. dy=float(y1-y0);

7. k=dy/dx;

8. y=float(y0);

9. x=float(x0);

10. Adjustment(x,x0,x1,y0,y1,y,k);

11. }

按照自己设置的像素，进行相应的调整：

1. void CDrawlineView::Adjustment(int &x,int &x0,int &x1,int y0,int y1,float &y,float &k) //判断后的画线算法DDA的精华

2. {

3. CClientDC dc(this);

4. dc.SetROP2(R2_NOTXORPEN);

5. CBrush bru;

6. bru.CreateSolidBrush(RGB(255,255,0));

7. CBrush *oldbru=dc.SelectObject(&bru);

8. if(Judge(k)==1||Judge(k)==3) //斜率为0~1或-1~0的时候

9. {

10. if(x0>x1) //向左上方画线

11. {

12. for(x=Ajust(x0,gaps);x>=x1;x-=gaps)

13. {

14. dc.Ellipse(x-int(gaps/2.0),Ajust(y,gaps)+int(gaps/2.0),x+int(gaps/2.0),Ajust(y,gaps)-int(gaps/2.0));

15. y=y-k*gaps;

16. }

17. }

18. else //向右下方画线

19. {

20. for(x=Ajust(x0,gaps);x<=x1;x+=gaps)

21. {

22. dc.Ellipse(x-int(gaps/2.0),Ajust(y,gaps)+int(gaps/2.0),x+int(gaps/2.0),Ajust(y,gaps)-int(gaps/2.0));

23. y=y+k*gaps;

24. }

25. }

26. }

27. if(Judge(k)==2||Judge(k)==4) //斜率大于1或小于-1的时候

28. {

29. k=1.0/k;

30. if(y0<y1) //向左上方画线

31. {

32. for(y=Ajust(y0,gaps);y<=y1;y+=gaps)

33. {

34. dc.Ellipse(Ajust(x,gaps)-int(gaps/2.0),y+int(gaps/2.0),Ajust(x,gaps)+int(gaps/2.0),y-int(gaps/2.0));

35. x=x+k*gaps;

36. }

37. }

38. else //向右下方画线

39. {

40. for(y=Ajust(y0,gaps);y>=y1;y-=gaps)

41. {

42. dc.Ellipse(Ajust(x,gaps)-int(gaps/2.0),y+int(gaps/2.0),Ajust(x,gaps)+int(gaps/2.0),y-int(gaps/2.0));

43. x=x-k*gaps;

44. }

45. }

46. }

47. dc.SelectObject(&oldbru);

48. }

 

 

1. int CDrawlineView::Judge(float k) //判断斜率K的大小

2. {

3. if(k>=0&&k<=1)

4. return 1;

5. if(k>1)

6. return 2;

7. if(k<0&&k>=-1)

8. return 3;

9. if(k<-1)

10. return 4;

11. return 4;

12. }

 

1. int CDrawlineView::Ajust(int x,int gaps) //调整x,y到自己的方格点上！

2. {

3. if(x%gaps>gaps/2)

4. x=int(x/gaps+1)*gaps;

5. else

6. x=int(x/gaps)*gaps;

7. return x;

8. }

其中主要关注直线的斜率的4种情况。然后4中不同的情况也可分为2种来画。

 

 

说的也不够详细，大概就是这样子，如若本文中有什么不对的地方，还请大家指出。谢谢。