OpenCV 最小二乘拟合方法求取直线倾角

    工业相机拍摄的图像中，因为摄像质量的限制，图像中的直线通过处理后，会表现出比较严重的锯齿。在这种状况下求取直线的倾角（其实就是直线的斜率），若是是直接选取直线的开始点和结束点来计算，或是用opencv自带的哈夫曼直线方法，都会引发较大的角度误差，通常会达到好几度。偏差这么大，显然达不到工控要求。后来尝试采起直线点集作最小二乘拟合，偏差缩小到0.5如下。如下是算法的代码：

 

	//最小二乘拟合计算直线的倾角
	int pointCount = pointVect.size();
	if (pointCount > 0)
	{
		int xCount = 0;
		int yCount = 0;
		int xyCount = 0;
		int xxCount = 0;
		for (int i = 0; i< pointCount; i++)
		{
			xCount += pointVect.at(i).x;
			yCount += pointVect.at(i).y;
			xyCount += (pointVect.at(i).x * pointVect.at(i).y);
			xxCount += (pointVect.at(i).x * pointVect.at(i).x);
		}
		double k = (double)(pointCount * xyCount - xCount * yCount) / (double)(pointCount * xxCount - xCount * xCount);
		double sinValue = - k / (sqrt(1 + k * k));
		double radian = asin(sinValue);
	        double pi = 3.1415926535;
	        double angle = radian * 180.0 / pi;
	}