摘 要 本文针对匀速直线运动模糊图像的参数确定问题，提出在MATLAB7.0下得到模糊图像的频谱图和自相关图，并用Matlab自带工具进行手动的方法确定模糊参数；实验结果证明模糊角度和理论值相差0.9°，模糊长度和理论值相差1个像素，与已有的几种方法相比较精确度更高，更有利于提高图像复原的清晰度。
关键词 傅立叶变换；点扩散函数；频谱图；自相关图
中图分类号 TP3 文献标识码 A 文章编号 1674-6708（2016）160-0153-02
在现实生活中我们遇到的好多图像或多或少都存在质量下降问题，对图像复原技术的研究最终归可结为对退化图像的点扩散函数的研究。目前国内外出现了多种估计运动模糊点扩散函数的算法；文献[1-2]中MasayukiTanaka提出运动模糊图像点扩散函数确定方法是：先构造出一个可以调节的三维锯齿波形作为进行检测的函数，然后把检测后的函数与模糊图像频谱进行运算，在运算过程中进行不断改变波形的周期和退化图像的频谱角度，当进行多次操作后使得函数取最小值时可以认为此时的角度就是暗条纹的角度。文献[3-5]中Mohsen-EbrahimiMoghaddam提出的方法是采用Radon变换方法来确定，其步骤是首先得到模糊图像的频谱图，然后对其进行Radon变换，当Radon变换到使得投影轴与频谱图中暗条纹相互垂直时，然后求出Radon变换取得最大值时所对应的方向，即确定点扩散函数的方向。本文就是以匀速直线运动产生的模糊图像为例，提出采用MATLAB7.0自带的工具进行手动确定模糊参数，此方法具有简单灵活性，而且有很高的准确性和稳定性，从而推广到非匀速直线运动中使用。
1 运动模糊图像的频谱特征分析
在MATLAB7.0环境下对其水平方向匀速直线运动产生模糊频谱图如图1所示。
K1=fft2（double（mf））；
M1=fftshift（K1）；
N1=abs（M1）；
P1=（N1-min（min（N1）））/（max（max（N1））-min（min（N1）））
*225；
subplot（1，2，2）；
imshow（P1）；
由于在实验过程中我们假定运动方向是水平的，运动长度L=32个像素，从图1实验结果的频谱图中可以明显的看出亮条纹，并且，这条亮条纹与水平方向是垂直的。
对于非水平方向运动模糊图像频谱图如图2所示。
图2的实验结果是图像运动方向与x轴之间有一定的角度，所以得到的亮条纹和X轴之间也存在着一定的角度。
2 匀速直线运动模糊方向和模糊长度的确定2.1 模糊方向的确定
在图2的基础上，通过多次实验发现，当模糊图像的角度发生变化时，条纹也随之改变，该亮条纹方向正好反映了运动模糊点扩散函数的方向性。于是我们只要求出这条亮条纹的倾角就可以了；所以将求模糊角度的问题转化为在平面上求一个三角形的角度问题；也就是下图中所示的θ角的度数。
根据频谱图用MATLAB7.0自带的作图工具Show Plot Tools作出三条直线构成一个三角形，然后在MATLAB7.0命令窗口中输入函数ginput（）来完成求解。用[x，y]=ginput（2），获得两点坐标，可以取好多个点，再求平均值然后利用k=（y（2）-y（1））/（x（2）-x（1））得到斜率，按照上面步骤求出度数为θ=69.1°，理论值为θ=70°，基本上吻合。
2.2 模糊长度的确定
本文首先采用Matlab7.0做出运动模糊图像的自相关图，生成的自相关图是一条曲线，具体如下图4所示，在这条曲线上出现了一对比较明显的负峰，并且方向朝下，通过多次实验结果显示这两负峰之间的距离L正好等于运动模糊长度的2倍。因此，只要能够计算出这两个负峰之间的距离就可以得出运动模糊长度L。
在图4的自相关图中利用matlab的Data Cursor工具可以测得两负峰之间的坐标分别为69与131，距离为62个像素点，两相关负峰间的距离等于运动模糊长度的2倍，因运动模糊的长度为31个像素，与理论值32个像素基本上吻合（比较精确）。
3 结 论
上面的实验结果基本上和理论值吻合，采用手动的方法来确定匀速直线运动模糊图像的点扩散函数，不但操作方法简单灵活，而且具有很高的精确度和稳定性。对于非匀速直线运动所造成的图像模糊，在某种程度上可以近似看成是多个小段匀速直线运动合成的。所以可以将此方法推广到非匀速直线运动模糊图像的处理上。
参考文献
[1]陈前荣，陆启生，成礼智.运动模糊图像点扩散函数尺度鉴别[J].计算机工程与应用，200（23）：15-19.
[2]黎玮.运动模糊图像复原的参数估计与算法研究[D].合肥：安徽大学，2011.
[3]贺卫国，黎绍发.匀速直线运动模糊长度的精确估计[J].计算机应用，2005，25（6）：1317-1320.
[4]张明谦.运动模糊图像的运动参数估计及复原研究[D].成都：电子科技大学，2009.
[5]杨文涛.分数傅里叶变换在数字图像处理中的应用研究[D].武汉华中科技大学博士学位论文，2009，12（3）：92-95.
[6]张德丰.MATLAB数字图像处理[M].北京：机械工业出版社，2012，3.
[7]博格斯，马科维奇.小波与傅里叶分析基础[M].2版.芮国胜，康健，译.北京：电子工业出版社，2010-2-1.
- 上合组织(连云港)国际物流园与区域经济发展互动关系研究
- 后疫情背景下服务业发展的现实困境与突破路径
- 跨境电商背景下中欧班列运行现状及优化对策研究
- 中国农产品电商全面转型升级
- 论新个人所得税法下的税务筹划
- 河南省驻马店市人口老龄化问题探析
- “一带一路”国际经济贸易纠纷解决机制探讨
- 基于PEST分析的新疆中小企业跨境电子商务发展研究
- “绿色氧吧”的营销策略分析
- 图书馆营销服务发展态势研究
- “互联网+”时代基于二次传播的直播影响力提升研究
- 非物质文化遗产的活态传承研究
- 古镇旅游开发的营商问题简析
- 承德非物质文化遗产旅游市场的现状研究
- 林业生态旅游发展现状分析及对策探讨
- 人口年龄结构变迁对资本流动的影响机制及对中国的启示
- “KOL化”与“去KOL化”
- 基于电商新零售模式的冷链物流发展分析
- 融媒体时代会展传播理念与方式创新
- 以新发展理念构建广东省会展业新发展格局
- O2O模式下的会展服务分析
- 基于网络文本数据的城市会展主题关注度研究
- “一带一路”倡议下国际铁路联运“一单制”规则构建研究
- 商标在先使用规则中应赋予在先使用人排他权的证伪
- 从税制简化改革到税法简化原则：理论拓补与法治完善
- deˌferred liability
- deˌferred share
- deˌferred taxation
- deˌfined benefit
- deˌfined contribution
- deˌlivered at frontier
- deˌlivered ex quay
- deˌlivered ex ship
- deˌmand-pull inˌflation
- deˌpositary receipt
- deˌpository receipt
- deˌpreciable life
- deˌpreciated cost
- deˌpreciation acˌcount
- deˌpreciation alˌlowance
- deˌpreciation fund
- deˌpreciation method
- deˌpreciation rate
- deˌpressed market
- deˌrivative instrument
- deˌrivative product
- deˌrived demand
- deˌsigner label
- deˌvolve sth into sth
- d-glucose
- 今文经学
- 今日
- 今日三，明日四
- 今日不知明日事，前人田土后人收
- 今日不知来日事
- 今日不谈明日事
- 今日事今日毕，留到明日更着急
- 今日河东，明日河西
- 今日爱东，明日爱西
- 今旦
- 今早
- 今早脱下鞋和袜，未审明朝穿不穿
- 今时
- 今昔
- 今昔之感
- 今昔何昔
- 今是前非
- 今是古，古是今
- 今是昔非
- 今是昨非
- 今晨
- 今有……于此（斯）
- 今朝
- 今朝有酒今朝醉
- 今朝有酒今朝醉，明日愁来明日当