我有一个使用 OpenCV 的 2D Mat 数据/图像。我想从中绘制一个 3D 图,行作为 x 轴,列作为 y 轴,每个坐标(行,列)上的像素值作为 z 轴。我想绘制类似于 Matlab 中的 surf 函数的图形。我尝试使用 bool loadFromData (double **data, unsigned int columns, unsigned int rows, double minx, double maxx, double miny, double maxy) 函数,但是当我调试到这一行时我的程序崩溃了。我实际上在使用这个 qwtplot3d 和 qt4 方面非常新,所以请帮我看看我哪里出错了。
我所做的是将所有像素值复制到一个数组中,然后尝试将其制成一个由行和列组成的二维数组,并将其传递给 loadFromData,如下所示:
//abs_outcorr is Mat type of opencv
int rows= abs_outcorr.rows;
int columns= abs_outcorr.cols * abs_outcorr.channels();
double fromMat[columns*rows];
for (int j=0; j<rows; j++)
{
double* data= abs_outcorr.ptr<double>(j); //access the pixel value from Mat
for (int i=0; i<columns; i++)
{
fromMat[i] = data[i];
}
}
double** x = new double * [rows];
for(int i=0; i<columns; i++)
x[i] = new double[rows];
for(int i=0; i<rows; i++)
for(int j=0; j<columns; j++)
x[j][i] = fromMat[columns*i + j];
for( int i=0; i<columns; i++)
delete x[i];
delete []x;
/* My main project is based on mainwindow with some buttons.
All these codes I implemented in a button clicked & I would like to pop out the graph,
but not sure the codes below work, please help me */
QWidget window;
window.setWindowTitle("plot 3d");
QHBoxLayout layout;
Qwt3D::SurfacePlot qsp(&window);
qsp.loadFromData(x, columns, rows, 0, columns, 0, rows);
qsp.setRotation(30,0,15);
qsp.setScale(1,1,1);
qsp.setShift(0.15,0,0);
qsp.setZoom(0.9);
for (unsigned i=0; i!=qsp.coordinates()->axes.size(); ++i)
{
qsp.coordinates()->axes[i].setMajors(7);
qsp.coordinates()->axes[i].setMinors(4);
}
qsp.coordinates()->axes[Qwt3D::X1].setLabelString("x");
qsp.coordinates()->axes[Qwt3D::Y1].setLabelString("y");
qsp.coordinates()->axes[Qwt3D::Z1].setLabelString("z");
qsp.setCoordinateStyle(Qwt3D::BOX);
qsp.updateData();
qsp.updateGL();
qsp.show();
layout.addWidget(&qsp);
window.setLayout(&layout);
window.resize(1000,800);
window.show();
//////////////////////////////////////////////////////////////
//These are some of my pixel value of my Mat (32x32)
//////////////////////////////////////////////////////////////
Row 0 Column 0 9.61749e-08
Row 0 Column 1 2.16608e-09
Row 0 Column 2 3.25873e-08
Row 0 Column 3 2.65754e-08
Row 0 Column 4 2.93116e-08
Row 0 Column 5 6.55923e-08
Row 0 Column 6 4.56592e-08
Row 0 Column 7 4.91113e-08
Row 0 Column 8 1.73816e-08
Row 0 Column 9 2.27045e-10
Row 0 Column 10 7.36088e-08
.
.
.
Row 16 Column 16 1
.
.
.
Row 31 Column 23 5.1846e-08
Row 31 Column 24 1.01708e-07
Row 31 Column 25 4.25331e-09
Row 31 Column 26 2.77903e-08
Row 31 Column 27 1.14044e-08
Row 31 Column 28 6.03817e-08
Row 31 Column 29 2.65248e-08
Row 31 Column 30 1.46648e-08
Row 31 Column 31 8.05808e-08
基本上这个 Mat/image 是中心的一个亮点,值为 1,而其余的全是黑色,几乎为 0。预期的结果应该看起来像图中间的一个尖锐的直峰。谢谢。