博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
图像维纳滤波实现(3)
阅读量:6353 次
发布时间:2019-06-22

本文共 852 字,大约阅读时间需要 2 分钟。

这里研究一下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

你可能感兴趣的文章
ts 格式化日期输出
查看>>
Optional
查看>>
sed 命令编辑文本
查看>>
LRUCache 具体解释
查看>>
Activity调用isDestroyed()方法报出,java.lang.NoSuchMethodError
查看>>
使用AFNetworking第三方下载类
查看>>
fhq-treap小结
查看>>
about porting
查看>>
MySQL事务及ACID特性
查看>>
Hadoop_31_MapReduce参数优化
查看>>
linux运维常见英文报错中文翻译(菜鸟必知)
查看>>
[原][osgEarth]添加自由飞行漫游器
查看>>
代码审查 Code Review
查看>>
fastjson如何指定字段不序列化
查看>>
[日常] Go语言圣经--示例: 并发的Echo服务
查看>>
BZOJ1969: [Ahoi2005]LANE 航线规划(LCT)
查看>>
linux内存管理之malloc、vmalloc、kmalloc的区别
查看>>
GreenDao 数据库升级 连接多个DB文件 或者指定不同的model&dao目录
查看>>
M1卡破解(自从学校升级系统之后,还准备在研究下)【转】
查看>>
vue 访问子组件示例 或者子元素
查看>>