标签：Mat int C++ ++ points 拟合 乘法 col row

#include <opencv2/opencv.hpp>   
#include <opencv2/highgui/highgui.hpp>

using namespace std;
using namespace cv;
/*
* @brief 得到最小二乘法拟合直线的系数矩阵X    A*X = B
* param[in]  points 单个虚线轮廓内的所有中点
* param[out] 最小二乘法拟合出的直线的系数  直线 y= ax + b，返回结果即为（a,b）
*/
Mat getParameterMatrix(vector<Point2i> points)
{
	// 最小二乘拟合
	// 构建矩阵A
	Mat A = Mat::zeros(2, 2, CV_64FC1);
	for (int row = 0; row < A.rows; row++)
	{
		for (int col = 0; col < A.cols; col++)
		{
			for (int k = 0; k < points.size(); k++)
			{
				A.at<double>(row, col) = A.at<double>(row, col) + pow(points[k].x, row + col);
			}
		}
	}
	//构建矩阵B
	Mat B = Mat::zeros(2, 1, CV_64FC1);
	for (int row = 0; row < B.rows; row++)
	{
		for (int k = 0; k < points.size(); k++)
		{
			B.at<double>(row, 0) = B.at<double>(row, 0) + pow(points[k].x, row) * points[k].y;
		}
	}
	// A*X = B， 其中X为要求解的系数矩阵
	Mat X;
	solve(A, B, X, DECOMP_LU);
	return X;
}

标签：Mat,int,C++,++,points,拟合,乘法,col,row
来源： https://blog.csdn.net/qq_41965957/article/details/121029597

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。