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;
	}