图像处理去噪课程设计

图像处理去噪课程设计

ID:39826137

大小:1.16 MB

页数:9页

时间:2019-07-12

图像处理去噪课程设计_第1页
图像处理去噪课程设计_第2页
图像处理去噪课程设计_第3页
图像处理去噪课程设计_第4页
图像处理去噪课程设计_第5页
资源描述:

《图像处理去噪课程设计》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、全变分去噪算法摘要图像去噪作为图像的预处理阶段,有着十分重要的意义。其目的就是去除图像的干扰信息,保护图像的局部特征,改善图像的质量,为图像的后续处理(如边缘检测,图像分割,特征提取等)提供可靠的保证,在图像处理领域占有非常重要的位置。而基于偏微分方程的图像去噪方法将数学与工程紧密结合,具有较强的自适应能力和灵活性,并且以此建立的模型容易修改,近年来在去除图像噪声和保护图像的边缘细节等信息等方面得到了快速发展。这种基于PDE的图像去噪方法主要源于约束最优化、能量最小化和变分法,其基本思想是以图像去噪的贝叶斯理论模型为基础,将所研究的问题转化为一个可以运用偏微分方程求解的方法求解

2、图像中相应的能量泛函极小值的过程,该极小值形式为一个偏微分方程,对该方程求解,所得到的解就是所要求的去噪复原图像。PDE图像去噪模型应用较为广泛,该模型将图像的梯度模作为边缘检测算子,通过检测图像的平坦区和边缘区的梯度大小而自适应地选取全变分模型和线性扩散模型,将图像边缘检测和图像去噪两个原本独立地两个过程有机地结合起来。本文将着重介绍基于PDE地图像去噪方法和TV算法。非线性偏微分方程去噪算法假设原图像为:u(x,y)可加性噪声为:n(x,y)(Enx=0,En2x=σ2)则噪声图像可表示为:u0x,y=ux,y+n(x,y)我们希望通过噪声图像得到一个原图像地估计,最小化问

3、题为:minimizeΩ(uxx+uyy)2(1)约束条件为:Ωu=Ωu0噪声Enx=0,En2x=σ2导致的约束:1ΩΩ(u-u0)2=σ2(2)此处Ω为图像区域Ω的面积(1)式和(2)式就定义了约束最优化问题。(由于TV范数的平移不变性:TVu+c=TV[u],c为任意常数,那么第一个约束实际上是已经满足的。因此,我们只需要考虑第二个拟合约束。通过引入拉格朗日乘子λ,可以定义一个新的能量泛函)新的能量泛函:Jux,y=Ω∇u(x,y)dxdy+λ2Ω(ux,y-u0(x,y))2dxdy(3)其中,参数λ对平衡去噪与平滑起重要作用。因此,它依赖于噪声水平。这样就建立了图像去

4、噪的TV模型。这是一个泛函求极值问题,即变分问题。偏微分方程(PDE)的推导:TV复原模型的欧拉-拉格朗日方程为:λu-u0-∇∙∇u∇u=0(5)推到过程如下:根据(1)式可得泛函中的:F=∇ux,y+λ2ux,y-u0x,y2=ux2x,y+uy2(x,y)+λ2(ux,y-u0(x,y))2,(6)此类泛函求极值的必要条件,即欧拉-拉格朗日方程为:FH-Fpx-Fqy=0其中FH=2(u-u0),Fq=λux∇u,Fq=Fp=λuy∇u代入(5)式得:λu-u0-∇∙∇u∇u=0(7)对于光滑函数u,方程(7)中得微分项表示水平线u=u(x0)得曲率,这就揭示了该模型得几

5、何特性。引用虚拟时间变量t,则有:ut=λu-u0-∇∙∇u∇u,(8)在这里我们选取λ(t):λ=-12σ2ux2+uy2-(u0)xuxux2+uy2+(u0)yuyux2+uy2dxdy,(9)由于∇u为分母,为了避免它为零,我们引入一个大于零得小参数,保证它大于零:∇u=ε2+∇u2只要ε保持足够小,就不会影响数字TV复原得性能。数值迭代方法如下:我们让:xi=ih,yi=jh,i,j=0,1,2⋯⋯⋯N,并且Nh=1tn=n∆t,n=0,1,⋯⋯⋯⋯,uijn=uxi,yi,tn,uij0=u0ih,jh+σφih,jh.迭代的初始条件是选定的,φ的均值为零并且二范数

6、为1.迭代公式为:uijn+1=uijn+∆th∆-x∆+xuijn∆+xuijn2+m∆+yuijn,∆-yuijn212+∆-y∆+yuijn∆+yuijn2+m∆+xuijn,∆-xuijn212-∆tλn(uijn-u0(ih,jh))(10)其中i,j=1,2,3⋯⋯⋯N边界条件为:u0jn=u1jn,uNjn=uN-1,jn,ui0n=uiNn=ui,N-1n在这里:∆∓xuij=∓(ui∓1,j-uij),∆∓yuij类似ma,b=(sgna+sgnb2)min⁡(a,b)λn=-h2σ2i,j∆+xuijn2+∆+yuijn2-∆+xuij0∆+xuijn∆+x

7、uijn2+∆+yuijn2-∆+yuij0∆+yuijn∆+xuijn2+∆+yuijn2,图像质量得比较参数:SNRSNR即信噪比,在这里我们定义信噪比得计算方法如下:SNR=-20lg⁡(u-gu)(10)其中u表示噪声图像,g表示复原图像,这里采用上式对SNR进行一个估计。编程实现TV算法,并对复原图像进行性能分析在matlab上编程实现TV复原算法,得到以下结果(γ=20):噪声图像信噪比:SNR=19.9344复原图像信噪比:SNR=20.0875对比发现,复原噪声信噪比增大,表

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。