图像维纳滤波实现（3）

这里研究一下psf2otf.m这个函数。主要关注的是函数怎样实现psf的size的变化的。

我们使用的是OTF = PSF2OTF(PSF,OUTSIZE)这个形式。同样的在psf2otf.m中存在parseinputs函数，用于检验参数和赋值变量。

[psf, psfSize, outSize] = ParseInputs(varargin{:});

该函数分别给psf，psfSize，outSize三个变量赋值。如果psf不全为0，则执行下面语句：

padSize = outSize - psfSize;

padSize表示需要填充的大小，这里理解起来还是比较直观，把psf扩充成跟outSize一样大的矩阵就可以了。比如说outSize为[100,100],psfSize[10,10],则padSize就为[90,90].

 psf     = padarray(psf, padSize, 'post');

该句表示填充图像或者填充数组，参数中不包含padval，即填充时默认用0进行填充，‘post’表示direction，表示在每一维后面进行填充。举个例子对[1 ,2;3 ,4]填充到5*5的矩阵就是如下效果：

                                    1   2   0   0   0

                                    3   4   0   0   0

                                    0   0   0   0   0

                                    0   0   0   0   0

                                    0   0   0   0   0

psf    = circshift(psf,-floor(psfSize/2));

接下来对新得到的矩阵进行周期平移circshift，其中-floor(psfSize/2)代表平移的方向，分为纵向和横向。继续上面的例子，平移方向矩阵为[-1,-1],就得到如下结果：

                                   4 0 0 0 3

                                   0 0 0 0 0 

                                   0 0 0 0 0 

                                   0 0 0 0 0

                                   2 0 0 0 1

最后使用otf = fftn(psf);求出我们需要的H(u,v)即可。

 

 

 

 

 

 

 

转载于:https://www.cnblogs.com/ideawu1001/archive/2012/11/15/2772234.html