Bresenham算法C++实现 任意斜率

算法原理参考参考链接

#include 
#include 
#include 
#include 
#include 
#include 
#include 
using namespace std;

void DrawBresenhamline(int x0, int y0, int x1, int y1)
{
  cv::Mat im(100,100,CV_8UC3, cv::Scalar(255,255,255));
  int dx = x1 - x0;//x偏移量
  int dy = y1 - y0;//y偏移量
  int ux = dx >0 ?1:-1;//x伸展方向
  int uy = dy >0 ?1:-1;//y伸展方向
  dx = abs(dx);
  dy = abs(dy);
  int dx2 = dx <<1;//x偏移量乘2
  int dy2 = dy <<1;//y偏移量乘2
  if(dx > dy)
    {//以x为增量方向计算
      int e = -dx; //e = -0.5 * 2 * dx,把e 用2 * dx* e替换
      int x = x0;//起点x坐标
      int y = y0;//起点y坐标

      for (x = x0; x != x1;x+=ux)
        {
          printf ("%d,%d\n",x, y);
          im.at(x,y) = cv::Vec3b(255,0,0);
          e= e + dy2;//来自 2*e*dx= 2*e*dx + 2dy  （原来是 e = e + k）
          if (e > 0)//e是整数且大于0时表示要取右上的点（否则是右下的点）
            {
              y += uy;
              e= e - dx2;//2*e*dx = 2*e*dx - 2*dx  (原来是 e = e -1)
            }
        }


    }
  else
    {//以y为增量方向计算
      int e = -dy; //e = -0.5 * 2 * dy,把e 用2 * dy* e替换
      int x = x0;//起点x坐标
      int y = y0;//起点y坐标
      for (y = y0; y != y1; y += uy)
        {
          printf ("%d,%d\n",x, y);
          im.at(x,y) = cv::Vec3b(255,0,0);
          e=e + dx2;//来自 2*e*dy= 2*e*dy + 2dy  （原来是 e = e + k）
          if (e > 0)//e是整数且大于0时表示要取右上的点（否则是右下的点）
            {
              x += ux;
              e= e - dy2;//2*e*dy = 2*e*dy - 2*dy  (原来是 e = e -1)
            }
        }

    }
  printf("%d, %d\n",x1,y1);
  im.at(x1,y1) = cv::Vec3b(255,0,0);
  cv::imshow("dd",im);

}

int main()
{
  cout << "Hello World!" << endl;

  DrawBresenhamline(80,0,0,20);

  cv::waitKey(0);
  return 0;
}

注意：行和列注意一下。