OpenCVを使用してHalconアルゴリズムを実装する（9）EdgesColorSubpixのCannyメソッドの実装

免責事項：この記事はインターネット上でオープンソースプロジェクトを共有することのみを目的としており、アルゴリズムはオリジナルではありません。

記事の内容は次の場所に複製されています。

Halcon Mat＆Input、Mat＆Output、doubleSigmaでのEdgesColorSubpixのCannyバージョンのC++実装について）; void NonMaxSuppresion（）; void LinkLine（）; void caclSubPixPosition（）; void HalconCannyEdgesColorAmp（Mat＆Input、Mat＆Out https：// blog .csdn.net / qq_39969166 / article / details / 119277349

void DerivateGauss(Mat &Input,Mat &Output,double Sigma);
void NonMaxSuppresion();
void LinkLine();
void caclSubPixPosition();
void HalconCannyEdgesColorAmp(Mat &Input,Mat &Output,double Sigma) {

	if (Input.channels() != 3) {
		return;
	}

	//Mat B, G, R;
	std::vector<Mat> bgrs;
	split(Input, bgrs);

	Mat Bx, By, Rx, Ry, Gx, Gy;
	//求高斯偏导数X和Y
	DerivateGauss(bgrs[0], Bx, Sigma, "x");
	DerivateGauss(bgrs[0], By, Sigma, "y");
	DerivateGauss(bgrs[1], Gx, Sigma, "x");
	DerivateGauss(bgrs[1], Gy, Sigma, "y");
	DerivateGauss(bgrs[2], Rx, Sigma, "x");
	DerivateGauss(bgrs[2], Ry, Sigma, "y");

	Mat gxx = Bx.mul(Bx) + Rx.mul(Rx) + Gx.mul(Gx);
	Mat gyy = By.mul(By) + Ry.mul(Ry) + Gy.mul(Gy);
	Mat gxy = Bx.mul(By) + Rx.mul(Ry) + Gx.mul(Gy);
	
	Mat temp, eValuesMat, eVectorsMat;
	double myArray[2][2];
	Output.create(Input.size(), CV_32F);
	for (int i = 0; i < Input.rows; i++)
	{
		for (int j = 0; j < Input.cols; j++)
		{
			myArray[0][0] = gxx.ptr<float>(i)[j];
			myArray[0][1] = gxy.ptr<float>(i)[j];
			myArray[1][0] = gxy.ptr<float>(i)[j];
			myArray[1][1] = gyy.ptr<float>(i)[j];
			//获得图像中像素位置的一个2X2梯度矩阵，矩阵中每个数据都表示一种梯度
			temp = cv::Mat(2, 2, CV_32F, myArray);
			//求该梯度矩阵的特征值，使用特征值的最大值代替表示该矩阵
			cv::eigen(temp, eValuesMat, eVectorsMat);
			//xld->attribute "response"
			Output.ptr<float>(i)[j] =(float) eValuesMat.ptr<double>(0)[0]; 
			//特征值最大的特征向量 atan2f(-dy,dx)=->xld attribute "angle"			
		}
	}
}