//Ax+by+cz=D
void cvFitPlane(const CvMat* points, float* plane){
// Estimate geometric centroid.
int nrows = points->rows;
int ncols = points->cols;
int type = points->type;
CvMat* centroid = cvCreateMat(1, ncols, type);
cvSet(centroid, cvScalar(0));
for (int c = 0; c<ncols; c++){
for (int r = 0; r < nrows; r++)
{
centroid->data.fl[c] += points->data.fl[ncols*r + c];
}
centroid->data.fl[c] /= nrows;
}
// Subtract geometric centroid from each point.
CvMat* points2 = cvCreateMat(nrows, ncols, type);
for (int r = 0; r<nrows; r++)
for (int c = 0; c<ncols; c++)
points2->data.fl[ncols*r + c] = points->data.fl[ncols*r + c] - centroid->data.fl[c];
// Evaluate SVD of covariance matrix.
CvMat* A = cvCreateMat(ncols, ncols, type);
CvMat* W = cvCreateMat(ncols, ncols, type);
CvMat* V = cvCreateMat(ncols, ncols, type);
cvGEMM(points2, points, 1, NULL, 0, A, CV_GEMM_A_T);
cvSVD(A, W, NULL, V, CV_SVD_V_T);
// Assign plane coefficients by singular vector corresponding to smallest singular value.
plane[ncols] = 0;
for (int c = 0; c<ncols; c++){
plane[c] = V->data.fl[ncols*(ncols - 1) + c];
plane[ncols] += plane[c] * centroid->data.fl[c];
}
// Release allocated resources.
cvReleaseMat(￠roid);
cvReleaseMat(&points2);
cvReleaseMat(&A);
cvReleaseMat(&W);
cvReleaseMat(&V);
}
呼叫的方式：
[cpp] view plain copy
CvMat*points_mat = cvCreateMat(X_vector.size(), 3, CV_32FC1);//定義用來存盤需要擬合點的矩陣
for (int i=0;i < X_vector.size(); ++i)
{
points_mat->data.fl[i*3+0] = X_vector[i];//矩陣的值進行初始化 X的坐標值
points_mat->data.fl[i * 3 + 1] = Y_vector[i];// Y的坐標值
points_mat->data.fl[i * 3 + 2] = Z_vector[i];<span style="font-family: Arial, Helvetica, sans-serif;">// Z的坐標值</span>
}
float plane12[4] = { 0 };//定義用來儲存平面引數的陣列
cvFitPlane(points_mat, plane12);//呼叫方程
我們擬合出來的方程：Ax+By+Cz=D
其中 A=plane12[0], B=plane12[1], C=plane12[2], D=plane12[3],
在網上看到這個擬合平面的程式，不是很懂X_vector.size()，X_vector[i]，Y_vector[i]，Z_vector[i]這三個引數的含義，是要把x,y,z的坐標分別放進那三個矩陣中嗎？那data.fl[]又是什么意思？求大神解答
要是有別的實作演算法也萬分感謝
uj5u.com熱心網友回復：
應該是要放進去，data是總體的uj5u.com熱心網友回復：
data是一個X_vector.size行3列的矩陣是嗎？uj5u.com熱心網友回復：
data是一個X_vector.size行3列的矩陣是嗎？轉載請註明出處，本文鏈接：https://www.uj5u.com/gongcheng/102535.html
標籤：圖形處理/算法