//頭文件
#include?“opencv2/opencv.hpp“
#include?
#include?
#include?
#include?“line_fit.h“

using?namespace?cv;
using?namespace?std;
//宏定義
#define?random（x）?（（rand（））%x）

//全局變量

//函數聲明




int?main（）
{
	/*載入原圖片*/
	//Mat?src?=?imread（“E://pic_F.L//Picture//OK1.jpg“?1）;//原始圖片
	Mat?src?=?imread（“OK1.jpg“）;
	Mat?gray;//原始的灰度圖
	cvtColor（srcgrayCOLOR_RGB2GRAY）;

	/*圖片預處理*/
	//平滑濾波使圖像的背景盡量統一
	Mat?_gray;
	bilateralFilter（gray_gray55*25/2）;?//圖像濾波
	imshow（“平滑處理的圖像“_gray）;
	//圖像增強，突出邊緣
	Mat?ruihua;
	Mat?kernel?=?（Mat_（33）?<???????????????????????????????????-1?5?-1
???????????????????????????????????0?-1?0）;?//定義增強圖像的卷積核
	filter2D（_grayruihuagray.depth（）kernel）;	//先濾波，再增強。高通濾波器銳化邊緣
	imshow（“增強后的圖像“ruihua）;

	/*邊緣檢測*/
	Mat?canny_output;//邊緣檢測的輸出圖像
	Canny（ruihua?canny_output?60?150?3）;
	imshow（“canny邊緣檢測的輸出圖“?canny_output）;

	/*利用標準霍夫變換檢測途中的直線*/
	////檢測線段
	//vector?lines;//存放檢測到的直線
	//Mat?H_result;
	//cvtColor（gray?H_result?CV_GRAY2BGR）;
	//HoughLines（canny_output?lines?1?CV_PI/180?150?0?0）;
	//畫出線段
	//for（size_t?i=0;?i	//{
	//	float?rho?=?lines[i][0]?theta?=?lines[i][i];
	//	Point?pt1pt2;
	//	double?a?=cos（theta）?b=sin（theta）;
	//	double?x0=a*rhoy0?=?b*rho;
	//	pt1.x=cvRound（x0?+?1000*（-b））;
	//	pt1.y=cvRound（x0?+?1000*（a））;
	//	pt2.x=cvRound（x0?-?1000*（-b））;
	//	pt2.y=cvRound（x0?-?1000*（a））;
	//	line（H_result?pt1?pt2?Scalar（255）?1?LINE_AA）;
	//}
	//imshow（“直線檢測的結果“?H_result）;
	
	/*利用累計概率霍夫變換檢測直線*/
	////檢測線段
	//vector?lines;
	//HoughLinesP（canny_output?lines?1?CV_PI/180?80?50?10）;
	//Mat?HP_result;
	//cvtColor（gray?HP_result?CV_GRAY2BGR）;
	////畫出線段
	//for（size_t?i=0;?i	//{
	//	Vec4i?l?=?lines[i];
	//	line（HP_result?Point（l[0]?l[1]）?Point（l[2]l[3]）?Scalar（255）1??LINE_AA）;
	//}
	//imshow（“直線檢測結果“?HP_result）;

	/*fitline（）擬合直線*/
	//獲取輪廓，得到一系列的輸入點集
	vector>?cal_contour;//存放輪廓
	vector?contours;
	vector?hierarchy;
	findContours（canny_output?contours?hierarchy?RETR_TREE?CHAIN_APPROX_NONE）;
	findContours（canny_output?cal_contour?hierarchy?RETR_TREE?CHAIN_APPROX_NONE）;
	cv::Mat?image?=?cv::Mat::zeros（gray.size（）CV_8UC3）;
	vector??lines（contours.size（））;
	double?dis_total=0.0;
	for（size_t?i=0;?i	{
		{
			fitLine（contours[i]?lines[i]?DIST_L2?0?1e-2?1e-2）;//擬合直線
			//畫線
			//獲取點斜式的點和斜率
			cv::Point?point0;
			point0.x?=?lines[i][2];
			point0.y?=?lines[i][3];
?
			double?k?=?lines[i][1]?/?lines[i][0];
?
			//計算直線的端點（y?=?k（x?-?x0）?+?y0）
			cv::Point?point1?point2;
			point1.x?=?0;
			point1.y?=?k?*?（0?-?point0.x）?+?point0.y;
			point2.x?=?640;
			point2.y?=?k?*?（640?-?point0.x）?+?point0.y;

			/*計算線性度指標*/
			double?dis?=?0.0;
			if（contours[i].size?>?0）
			{
				for（int?k=0;?k				{
					for（int?l=0;?l