费了半天劲,终于把这个WMesh类搞懂了,可惜效果不佳,比Matlab中的mesh差多了。
使用WMesh前,需要有一个Mesh对象,Mesh是三维数据点的基本几何信息、颜色信息、索引信息等集成的对象。Mesh对象无法直接在Viz3d中显示,需要转换为WMesh对象,然后才可以显示。
首先了解Mesh类:
我们先看看Mesh的定义,在Mesh中有四个Mat型的成员变量,分别是cloud、colors和normals。
望名知意,cloud表示点云,它实际上是点云中每个点的坐标组成的二维矩阵,它的通道数为3,所以cloud的数据类型为:CV_32FC3,cloud的三个通道[0]、[1]、[2]分别对应坐标的x、y、z。对cloud还有一个特殊要求,那就是一行、N列,也就是一个一位矢量。需要预先把每个点坐标依次赋值给cloud的一些列元素,比如点P(x,y,z),赋值给cloud的第0行、j列元素:cloud.ptr
在Mesh中的第二个成员colors也是一个一维矩阵,1行、N列,与cloud的尺度一致。colors通常是3通道或四通道,但绝不以是2通道。三通道时,分别对应红、绿、蓝三原色。colors的每个元素与cloud中的元素一一对应,分别给cloud中对应元素赋予颜色。还有一点需要注意,colors的每个通道上的数据类型必须是CV_8U,因此colors的数据类型通常为CV_8UC3。
第三个公有成员是normals,它与colors类似,具体内容还没有研究。
第四个共有成员是polygons,它也是一个一维行矢量,但是列数要比cloud多一列,在polygons的第一个
还有一个公有静态成员函数load,大家自己研究吧!
class CV_EXPORTS Mesh
{
public:
Mat cloud, colors, normals;
//! Raw integer list of the form: (n,id1,id2,...,idn, n,id1,id2,...,idn, ...)
//! where n is the number of points in the poligon, and id is a zero-offset index into an associated cloud.
Mat polygons;
//! Loads mesh from a given ply file
static Mesh load(const String& file);
};
接下来,我们以一个实例来说明Mesh的用法。在本示例中,我读取了一幅图像,将图像的灰度值作为z值,这样就可以以三维的形式观察二维图像数据:
#include <opencv2/viz/vizcore.hpp>
#include <opencv2/opencv.hpp>
#include <iostream>
using namespace cv;
using namespace std;
//#include "myfunction.h"
int main()
{
Mat img=imread("D:/CodeWork/MyImage/opencv_logo.png");
Mat gray;
img.convertTo(gray,CV_32FC1);
int rows=img.rows,cols=img.cols;
//为了赋值方便,这里我用了与img同行、同列的矩阵
Mat cloud(rows,cols,CV_32FC3);
int cx=cols/2,cy=rows/2;
//为cloud矩阵赋值
for(int i=0;i<rows;i++)
{
float y=i-cy;
for(int j=0;j<cols;j++)
{
float x=j-cx;
cloud.ptr<Vec3f>(i)[j][0]=x;
cloud.ptr<Vec3f>(i)[j][1]=-y;
cloud.ptr<Vec3f>(i)[j][2]=(img.ptr<Vec3b>(i)[j][0]
+img.ptr<Vec3b>(i)[j][1]
+img.ptr<Vec3b>(i)[j][2])/3;
}
}
///创建polygons矩阵,并为其赋值,它的第0个元素是三维点的个数
/// 剩下的元素值分别为1、2、3......、N。此次的N表示总的三维点数
Mat polygon(1,cols*rows+1,CV_32SC1,Scalar::all(0));
polygon.ptr<int>(0)[0]=rows*cols;
for(int i=1;i<=polygon.cols;i++)
polygon.ptr<int>(0)[i]=i;
///下面是把二维的cloud转化为一维的cloud
Mat cloud1d=cloud.reshape(0,cloud.rows*cloud.cols);
cloud1d=cloud1d.t();
///创建并赋值Mesh对象
viz::Mesh mesh;
mesh.cloud=cloud1d;
mesh.polygons=polygon;
///创建颜色矩阵
Mat color;
color=img.reshape(0,rows*cols);
color=color.t();
mesh.colors=color;
///mesh本身无法在三维窗口中显示,需要经过WMesh处理,得到WMesh对象wmesh
viz::WMesh wmesh(mesh);
///显示wmesh
imshow("original image",img);
viz::Viz3d my3DWin("My 3D Window");
my3DWin.setBackgroundColor(viz::Color::cyan());
my3DWin.showWidget("wmesh",wmesh);
my3DWin.spin();
return 0;
}