OpenCV idft() 和MATLAB ifft2的结果不匹配

Question

所以我正在 MATLAB 中测试我的算法并且它已经完成了。

然后现在使用 OpenCV 2.4.5 进行 C++ 移植的编码。

问题是两个平台OpenCV和MATLAB的逆傅立叶变换方法。

所以我用简单的矩阵进行了测试。

下面是测试结果。

主题矩阵是... 3 x 3 2-D。

   1 2 3
   4 5 6
   7 8 9

-MATLAB-

test = [ 1, 2, 3;
         4, 5, 6;
         7, 8, 9];

ifft2（测试）；

结果

   5.0000 + 0.0000i  -0.5000 - 0.2887i  -0.5000 + 0.2887i
  -1.5000 - 0.8660i   0.0000 + 0.0000i   0.0000 + 0.0000i
  -1.5000 + 0.8660i   0.0000 + 0.0000i   0.0000 + 0.0000i

-OPENCV-

注意：元素是相同的值。

Mat a = Mat::zeros(3, 3, CV_64FC1);
Mat b = Mat::zeros(3, 3, CV_64FC1);

 a.at<double>(0,0) = 1;
 a.at<double>(0,1) = 2;
 a.at<double>(0,2) = 3;

 a.at<double>(1,0) = 4;
 a.at<double>(1,1) = 5;
 a.at<double>(1,2) = 6;

 a.at<double>(2,0) = 7;
 a.at<double>(2,1) = 8;
 a.at<double>(2,2) = 9;

 idft(a, b, DFT_SCALE, 0);

结果

    4.33333     -4.13077     2.79743
   -2.10313     -0.103134   -2.83518
   -0.563533     2.16852     1.43647

我仍然没有找到解决方案。 即使这也不能给我一个解决方案。

编辑：问题已解决。 我把 CV_64FC1 作为输入和 CV_64FC2 作为输出。 两个矩阵必须具有相同的深度，输入和输出都必须是 64_FC2。 并且标志 DFT+COMPLEX_OUTPUT+DFT_SCALE 与 MATLAB 的 ifft2 相同。

-解决了-

Mat input = Mat::zeros(3, 3, CV_64FC2); 
Mat output = Mat::zeros(3, 3, CV_64FC2);
idft(input, output, DFT_COMPLEX_OUTPUT+DFT_SCALE, 0);

Answer 1

我相信你需要cv::DFT_COMPLEX_OUTPUT+cv::DFT_SCALE因为输入到idft显然是一个复杂的数值矩阵结果。

此外，我认为您将需要一个 2 通道数组用于输出（类型CV_64FC2 ），与输入类似。 与OpenCV 中的任何多通道图像一样，您可以访问具有适当向量类型的元素（例如，对于双精度型， .at<cv::Vec2d>(i,j) ，其中Vec2d将实部和虚部存储在位置i,j )。

Answer 2

我认为如果您使用 2 通道输入矩阵（CV_64FC2），您应该使用

a.at<Vec2d>(0,0)[0] = 1; // Re - part
a.at<Vec2d>(0,0)[1] = 0; // Im - part

代替：

a.at<double>(0,0) = 1;