有以下两种方法:
1.图像金字塔
利用pyrUp、pyrDown函数,但输入图像必须是2的幂,否则报错。
从下面的博客才知道有“2的幂”这个要求,对博主表示感谢!
https://blog.csdn.net/qq_18343569/article/details/47829863
采用1400*1400尺寸的仿真数据
// 下采样
Mat projectionMatDwn(rows / rowDwn, cols / colDwn, CV_32FC1, projectionDwn);
pyrDown(projectionMat, projectionMatDwn, Size(cols / colDwn, rows / rowDwn));
结果报错。
由于项目中图像尺寸不一定是2的幂,放弃该方法。
注意,pyrDown及下面的resize函数,要求Size变量内先列后行!
2.矩阵插值
利用resize函数
2.1 Opencv resize函数与matlab imresize函数存在区别,见下博客。
https://blog.csdn.net/qq_38676868/article/details/99406533
2.2 resize函数使用注意方法,参考以下博客。
https://blog.csdn.net/u012005313/article/details/51943442
Matlab resize中的插值方法(来自help文档)
为了保持与MATLAB相同的输出,项目采用了双线性插值(默认)
Mat IobjectMD; // 下采样物体透射数据
resize(IobjectM, IobjectMD, Size(IcolsDwn, IrowsDwn), INTER_LINEAR); // 双线性插值(默认)
另外,经验证,opencv将内存指针(如float *)封装为MAT,后续Mat的处理均针对该内存。不必再利用for循环赋值,将Mat数据赋值到内存指针。
MyScatterCorrect(projections, projections0, detU, detV, detUDwn, detVDwn);
cout << endl << "projections" << endl;
cout << projections[0] << " " << projections[1] << " " << projections[2] << endl;
IobjectM = IobjectM.mul(1-SF);
cout << endl << "IobjectM" << endl;
cout << IobjectM.at<float>(0, 0) << " " << IobjectM.at<float>(0, 1) << " " << IobjectM.at<float>(0,2) << endl;
即MyScatterCorrect输出的projections,与函数内变量IobjectM对应相同的内存。