使用MathNet傅里叶变换做卷积时图像象限互换问题

你的问题核心是FFT的频谱布局和核的位置不匹配导致的循环移位：

• 标准FFT的输出中，直流分量（频率为0的分量）默认位于矩阵的左上角(0,0)，而你把核的唯一非零值放在了中心(256,256)。
• 根据FFT的移位定理：空间域的循环移位对应频域的相位调制。当你对中心有冲激的核做FFT时，相当于把原本在左上角的直流分量“移”到了中心，和图像FFT相乘后，逆FFT就会给原图做一次反向的循环移位，最终表现为四个象限完全互换。

解决方法 #

方法1：调整核的位置 #

直接把核的非零值放在左上角(0,0)，而非中心：

Matrix<Complex> kernel = new Matrix<Complex>(512, 512);
kernel[0, 0] = 1.0; // 把冲激放在FFT默认的直流分量位置

这样执行你原本的FFT-相乘-逆FFT流程，结果就不会出现象限互换。

方法2：添加频谱移位操作 #

如果必须保留核的中心位置，需要在FFT后对频谱做移位，对齐直流分量：

Matrix<Complex> image = GetImage();
Matrix<Complex> kernel = GetKernel();

// 正向FFT
Fourier.Forward2D(image, FourierOptions.NoScaling);
Fourier.Forward2D(kernel, FourierOptions.NoScaling);

// 对核的FFT结果做频谱移位，将直流分量移回中心
Fourier.Shift2D(kernel);

// 频域相乘
Matrix<Complex> convolved = image.PointwiseMultiply(kernel);

// 逆移位恢复，再做逆FFT
Fourier.Shift2D(convolved);
Fourier.Inverse2D(convolved, FourierOptions.NoScaling);

方法3：逆FFT后手动修正象限 #

如果上述方法不可行，也可以在逆FFT完成后，直接对结果图像做象限互换：

// 执行完原流程后，对convolved做象限交换
SwapQuadrants(convolved); // 自定义函数，交换四个象限：1↔4，2↔3

自定义SwapQuadrants的逻辑：将四个256x256的子矩阵互相交换位置即可。

内容的提问来源于stack exchange，提问作者redcurry