C++自定义非极大值抑制（Canny边缘检测，亚像素方法）

Mat subPixleRestrainNoMax2(Mat img)//具有更清楚的物理解释，线插分法
{
   
	Mat oriImage = img;
	Mat newImage_x = Mat(img.size(), CV_32SC1, Scalar(0));
	Mat newImage_y = Mat(img.size(), CV_32SC1, Scalar(0));
	Mat newImage_xy = Mat(img.size(), CV_32SC2, Scalar(0));//包含梯度信息
	Mat maxImage = Mat(img.size(), CV_8UC1, Scalar(0));//非极大值抑制后的图像

	//对x方向求Sobel算子梯度
	int width = img.cols;
	int hight = img.rows;

	int max_sobel_x = 0;//记录x方向最大梯度绝对值
	for (int i = 1; i < hight - 1; i++)
		for (int j = 1; j < width - 1; j++)
		{
   
			int nav_value = oriImage.at<uchar>(i + 1, j - 1) + 2 * oriImage.at<uchar>(i, j - 1)
				+ oriImage.at<uchar>(i - 1, j - 1);
			int pos_value = oriImage.at<uchar>(i + 1, j + 1) + 2 * oriImage.at<uchar>(i, j + 1)
				+ oriImage.at<uchar>(i - 1, j + 1);
			int sobel_result = ((pos_value - nav_value) / 6);//相减除以6
			//int sobel_result = (pos_value + 6 * 255 - nav_value) / (6 * 2); //把像素加255 再除以 2
			if (max_sobel_x < abs(sobel_result))
				max_sobel_x = abs(sobel_result);//找到梯度最大值
			newImage_x.at<int>(i, j) = sobel_result;
		}
	float mulit_x = (float)(255 / max_sobel_x);//按倍数扩大像素值，使像素差距放大
	for (int i = 1; i < hight - 1; i++)
		for (int j = 1; j < width - 1; j++)
		{
   
			newImage_x.at<int>(i, j) = (int)(newImage_x.at<int>(i, j) * mulit_x);
		}

	//对y方向求Sobel梯度
	int max_sobel_y = 0;//记录最大y方向梯度绝对值
	for (int i = 1; i < hight - 1; i++)
		for (int j = 1; j < width - 1; j++)
		{
   
			int nav_value = oriImage.at<uchar>(i - 1, j - 1) + 2 * oriImage.at<uchar>(i - 1, j)
				+ oriImage.at