Question

我今天才意识到，MATLAB的imfilter默认卷积方法根本不是卷积，而是相关性。

我想知道，这会如何改变deconvblind的结果？它使用最大似然算法，但对于与相关卷积的图像或通过卷积卷积的图像是这样做的吗？

请帮帮我！我在MathWorks网站上找不到任何相关信息。

Answer 1

它或多或少使用卷积。但是，这在空间域中不，但在频域中。通过在MATLAB命令提示符中使用deconvblind查看open deconvblind的实际源代码，可以瞥见方法的作用。具体来说，你想知道的相关部分是在源头的中间：

% 2.c Determine next iteration image & apply positivity constraint
J{3} = J{2};
H = psf2otf(P{2},sizeI);
scale = real(ifftn(conj(H).*fw)) + sqrt(eps);
J{2} = max(Y.*real(ifftn(conj(H).*CC))./scale,0);
clear scale;
J{4} = [J{2}(:)-Y(:) J{4}(:,1)];
clear Y H;

单元阵列J包含去卷积图像，您可以看到图像是使用频域中两个图像的乘积的逆FFT获得的，这基本上是卷积。

如果您回想起傅立叶变换的属性，则空间域中的卷积是频率的倍增。我不会通过实际的算法，因为你只是想知道它是否是卷积或相关，但是MathWorks确实引用了deconvblind算法使用的一些相关论文。

请参阅：

＆＃34;由D.S.C.加速迭代图像恢复算法。比格斯和M. Andrews，Applied Optics，Vol。 36，第8期，1997年。
＆＃34;哈勃太空望远镜图像和光谱的解卷积＆＃34;，R.J。 Hanisch，R.L.White和R.L. Gilliland。在＆＃34;反卷积图像和Spectra＆＃34;，Ed。 P.A. Jansson，第二版，Academic Press，CA，1997。
＆＃34;通过最大似然重建的光学显微图像 Deconvolution＆＃34;，Timothy J. Holmes等。在＆＃34;手册中生物共聚焦显微镜＆＃34;，Ed。詹姆斯·B·鲍利，全会 Press，New York，1995

FWIW，过滤图像时的相关性和卷积，如果卷积模糊是对称的，则。卷积采用内核并在元素乘积和求和之前将其旋转180度。相关不会执行此旋转，但元素方式的乘积和求和是相同的。如果内核是对称的，这意味着将蒙版旋转180度会给你相同的蒙版。如果掩码不是这样，您可以在'conv'中指定imfilter标志以指定您要执行卷积。

imfilter的默认值是相关性。 deconvblind如何解卷积？

1 个答案: