1. //Ax+by+cz=D  
  2. void cvFitPlane(const CvMat* points, float* plane){  
  3.     // Estimate geometric centroid.  
  4.     int nrows = points->rows;  
  5.     int ncols = points->cols;  
  6.     int type = points->type;  
  7.     CvMat* centroid = cvCreateMat(1, ncols, type);  
  8.     cvSet(centroid, cvScalar(0));  
  9.     for (int c = 0; c<ncols; c++){         //用地址的形式才能实现先列后行的操作，如果用Mat型只能先行后列
  10.         for (int r = 0; r < nrows; r++)  
  11.         {  
  12.             centroid->data.fl[c] += points->data.fl[ncols*r + c];  
  13.         }  
  14.         centroid->data.fl[c] /= nrows;  
  15.     }  
  16.     // Subtract geometric centroid from each point.  
  17.     CvMat* points2 = cvCreateMat(nrows, ncols, type);  
  18.     for (int r = 0; r<nrows; r++)  
  19.         for (int c = 0; c<ncols; c++)  
  20.             points2->data.fl[ncols*r + c] = points->data.fl[ncols*r + c] - centroid->data.fl[c];  
  21.     // Evaluate SVD of covariance matrix.  
  22.     CvMat* A = cvCreateMat(ncols, ncols, type);  
  23.     CvMat* W = cvCreateMat(ncols, ncols, type);  
  24.     CvMat* V = cvCreateMat(ncols, ncols, type);  
  25.     cvGEMM(points2, points, 1, NULL, 0, A, CV_GEMM_A_T);  //points也要为CVMat型
  26.     cvSVD(A, W, NULL, V, CV_SVD_V_T);  
  27.     // Assign plane coefficients by singular vector corresponding to smallest singular value.  
  28.     plane[ncols] = 0;  
  29.     for (int c = 0; c<ncols; c++){  
  30.         plane[c] = V->data.fl[ncols*(ncols - 1) + c];  
  31.         plane[ncols] += plane[c] * centroid->data.fl[c];  
  32.     }  
  33.     // Release allocated resources.  
  34.     cvReleaseMat(¢roid);  
  35.     cvReleaseMat(&points2);  
  36.     cvReleaseMat(&A);  
  37.     cvReleaseMat(&W);  
  38.     cvReleaseMat(&V);  
  39. }  

调用的方式：

[cpp]  view plain  copy
 
  1. CvMat*points_mat = cvCreateMat(X_vector.size(), 3, CV_32FC1);//定义用来存储需要拟合点的矩阵   
  2.         for (int i=0;i < X_vector.size(); ++i)  
  3.         {  
  4.             points_mat->data.fl[i*3+0] = X_vector[i];//矩阵的值进行初始化   X的坐标值  
  5.             points_mat->data.fl[i * 3 + 1] = Y_vector[i];//  Y的坐标值  
  6.             points_mat->data.fl[i * 3 + 2] = Z_vector[i];<span style="font-family: Arial, Helvetica, sans-serif;">//  Z的坐标值</span>  
  7.   
  8.         }  
  9.         float plane12[4] = { 0 };//定义用来储存平面参数的数组   
  10.         cvFitPlane(points_mat, plane12);//调用方程   
我们拟合出来的方程：Ax+By+Cz=D

其中 A=plane12[0],    B=plane12[1],   C=plane12[2],   D=plane12[3],

这是要注意的方程的表示

GitHub 加速计划 / opencv31 / opencv
77.38 K
55.71 K
下载
OpenCV: 开源计算机视觉库
最近提交(Master分支：2 个月前 )
c3747a68 Added Universal Windows Package build to CI. 10 天前
9b635da5 - 10 天前
Logo
GitCode 开源社区

旨在为数千万中国开发者提供一个无缝且高效的云端环境，以支持学习、使用和贡献开源项目。

更多推荐

[转载]在Windows环境下安装GNU Radio

转自：在Windows环境下安装GNURadio_恐弱智_新浪博客GNU Radio是用Python开发的，大部分开源的工程能够在Linux环境下运行良好，而Windows下却运行的很勉强，而且安装配置都很复杂。GNU Radio算是个例外了，不光提供了Windows的二进制安装，还有比较详细的说明。我是Python小白，所以折腾了好久才弄好，特意记录下来，免得以后再装还折腾。GNU Radio的

avatar GitCode 开源社区

centOS 8 使用dnf安装Docker

DNF是什么？CentOS 8使用YUM软件包管理器版本v4.0.4。现在，该版本使用DNF(已删除YUM)。DNF是软件包管理器。它会在Linux发行版上安装，执行更新并删除软件包。使用DNF安装Docker跳过具有损坏依赖性的程序包一个有效的解决方案是使您的CentOS 8系统使用以下--nobest命令安装最符合条件的版本：sudo dnf install docker...

avatar GitCode 开源社区

定时同步数据库表(mysql+linux+crontab)

sync.sh里面的参数需要改变，ip/username/password/database/tablesync.sh#!/bin/sh# Please change the IP and password of the data source db.# Then change the table name.filename=/home/nington/db/$(date +%Y-%m

avatar GitCode 开源社区