基于dsp的图像去噪实现
基于快速曲波变换的图像去噪算法

ti m to ( S F C )nt0l ye sd- o e mae wt i etP a i a-o N i a o v us (S R= 0828 wt hs e d C - D T o ny id e ni d i g i hg s ek Sgl t— os R t a e P N 3 .5 i h l s h h n e i l h
n ie v ra c r 2 ,u s c iv s b s i a u i . os ain e t= 5) b ta o a he e etvs lq a t l u l y
Ke r s a t d s rt u v l t t n fr ; u v lt w a p n g r m ; S F T ag r m ; ice s i n n y wo d :f s i e e c r ee r s m c re e r p i g a o i c a o l h t C — DC o t l i h c r l p n i g
En i e r g a d A piain ,0 7 4 ( ) 3 - 3 gn e i n p l t s 2 0 , 3 6 :1 3 . n c o
A s a t uvls a e e e rsn nst p o bet wt dsot ut s aog eSsb tte uvl 9 t nfr bt c:C re t n bR r r ee ta i r y fr ojcs i i ni ie ln de ,u h cre t9 r s m r e c p oo h c n i e a o
E m i jyn 36 @vp13cm — a : ag 6 8 i.6 . lh o
YANG i-h n XU Ca h i W ANG If - a . a e e n i n l o i m ae n fs u v lt r n frD. mp t r Ja o g, n- u , t D nI g d - os g ag rt o m i h b sd o at c r ee ta so l Co ue
基于小波变换的图像去噪算法研究与应用

基于小波变换的图像去噪算法研究与应用一、引言图像去噪是图像处理领域的重要问题,随着数字图像处理技术的发展与应用,对图像的去噪要求越来越高。
因此,在图像领域中,图像去噪一直是研究的热点之一。
二、小波变换小波变换是一种信号处理方法,可以用于信号的压缩、去噪、特征提取等。
小波变换通过分析信号中的局部细节信息,可以将信号分解为不同频率的子带,从而更好地处理信号中的各个部分。
三、小波变换在图像去噪中的应用1.小波阈值去噪法小波阈值去噪法是一种基于小波分解的图像去噪方法,该方法通过分解图像为不同频率的小波子带,再对各自的子带进行去噪处理,最后将各子带结果合成为一张图像。
该方法的核心在于确定小波子带的阈值,目前常用的方法有软阈值和硬阈值两种。
软阈值和硬阈值的区别在于,软阈值会使小于阈值的子带信号变为0,但不会对大于阈值的信号做限制;硬阈值和软阈值类似,只是会使小于阈值的子带信号全部变为0。
2.双阈值小波去噪法双阈值小波去噪法是一种基于小波变换的两阶段去噪方法,该方法首先通过小波分解将图像分解为不同频率的小波子带,然后采用两个阈值对各子带进行去噪处理,其中一个阈值用于对高频子带进行去噪,另一个阈值用于对低频子带进行去噪。
该方法的主要优点在于,可以有效地去除噪声的同时,尽可能地保留图像中的细节和纹理信息。
四、实验分析与结果本文选择了几组不同的噪声图像进行去噪处理,将分别采用小波阈值去噪法和双阈值小波去噪法进行实验处理。
实验结果表明,采用小波阈值去噪法能够显著地去除高斯噪声和椒盐噪声;双阈值小波去噪法在去除图像噪声的同时,能够有效地保留图像中的细节信息。
五、结论小波变换是一种重要的信号处理方法,在图像去噪方面得到了广泛的应用。
通过实验对比,小波阈值去噪法和双阈值小波去噪法均能达到不错的去噪效果,可根据不同的噪声类型和噪声强度进行选择和应用。
未来,小波变换方法预计将得到更广泛的应用,为图像处理及相关领域的研究提供更有力的工具和技术。
基于DSP机器视觉的道路图像采集和处理系统设计

同 时提 出 了道路 图像 处理 算 法流 程 以及 具 体 的程 序框 架 ,最终 在软 件 集 成 开发 环 境 CCS中生 成 应 用程 序 以及 P 1 动程序 ,在 多次试验 中取得 了不错 的效 果 。 C 驱
关 键 词 : 图像 采 集 ; 图 像 处 理 ;DS P;P 总 线 C1
码 芯片 S A 1A 完成 A D 转换 ,S A 1 A允许 四路 模 A 71 / A 71 拟视 频 输入 ,具有 两 个 模 拟处 理 通道 ,可 以对 摄 像 头输
据 通过 C D数字 相机 输入 的模 拟 图像信 号 ,经视频 解码 C 芯 片转 换 成数 字信 号 ,由 C L ( o lx P ormma l P D C mpe rga be L gc D vc ) 控 制 经 过 FF 输 入 D P进 行 图 像 的增 o i eie IO S 强 、去 噪 、边缘 检测 等 ,为实 现 实 时地将 原 始 图像 和 处
S A DA D 和 F T N R W 两 种 工 作 模 式 。 按 照 C I 6 1格 CR 0 式 ,Y UV 图 像 分 辨 率 为 7 0x 7 , 当 按 行 输 出 时 , 2 5 6 S A 1A 输 出 数 据 流 大 小 为 :7 0 l = 4 0 i 因 为 A 71 2 x 6 14 bt 。
I
{
CL 控 器 PD 制
— ]r —1 一
儿
I址 换 功 存I 地 切 码、能 I 储 译等 帧)
i j
} —C — 广 —P — 一 1 — —
l 显不及后 期处理 ) ( 1
用此 开发 过 程检 测新 算法 的可行 性 ,可 以进行 图像 增强 去 噪 、视频 压缩 、边缘检 测 、视 觉定 位等算 法 的研 究 。
基于剪切波变换和拟合优度检验的遥感图像去噪

第61卷第5期吉林大学学报(理学版)V o l.61 N o.5 2023年9月J o u r n a l o f J i l i nU n i v e r s i t y(S c i e n c eE d i t i o n)S e p2023d o i:10.13413/j.c n k i.j d x b l x b.2022197基于剪切波变换和拟合优度检验的遥感图像去噪成丽波,陈鹏宇,李喆,贾小宁(长春理工大学数学与统计学院,长春130022)摘要:针对遥感图像中的高斯白噪声,提出一种基于剪切波变换和拟合优度检验的遥感图像去噪算法.首先,将含噪遥感图像通过剪切波变换多尺度分解得到不同子带,利用剪切波域下高斯白噪声系数的统计关系估计去噪阈值;其次,计算高频子带的拟合优度检验统计量,将统计量与去噪阈值相比较进行去噪;最后,对系数矩阵进行剪切波逆变换重建去噪图像.仿真实验结果表明,该算法能有效去除遥感图像中的高斯噪声,保持图像的边缘纹理信息,并且在不同噪声水平下,均获得了较高的峰值信噪比,其中与剪切波阈值去噪算法相比平均提高0.33d B.关键词:遥感图像;剪切波变换;拟合优度检验;图像去噪中图分类号:T P341.4文献标志码:A 文章编号:1671-5489(2023)05-1187-08 R e m o t e S e n s i n g I m a g eD e n o i s i n g B a s e d o nS h e a r l e tT r a n s f o r ma n dG o o d n e s s o f F i t T e s tC H E N GL i b o,C H E NP e n g y u,L I Z h e,J I A X i a o n i n g(S c h o o l o f M a t h e m a t i c s a n dS t a t i s t i c s,C h a n g c h u nU n i v e r s i t y o f S c i e n c e a n dT e c h n o l o g y,C h a n g c h u n130022,C h i n a)A b s t r a c t:A i m i n g a tw h i t eG a u s s i a nn o i s e i nr e m o t es e n s i n g i m a g e s,w e p r o p o s e dar e m o t es e n s i n g i m a g e s d e n o i s i n g a l g o r i t h m b a s e do ns h e a r l e t t r a n s f o r m a n d g o o d n e s so f f i t t e s t.F i r s t l y,t h en o i s y r e m o t e s e n s i n g i m a g ew a s d e c o m p o s e d i n t o d i f f e r e n t s u b-b a n d s t h r o u g h s h e a r l e t t r a n s f o r ma tm u l t i p l e s c a l e s,a n d t h e d e n o i s i n g t h r e s h o l dw a s e s t i m a t e du s i n g t h e s t a t i s t i c a l r e l a t i o n s h i p o fw h i t eG a u s s i a n n o i s e c o e f f i c i e n t s i n t h e s h e a r l e t d o m a i n.S e c o n d l y,w e c a l c u l a t e d t h e g o o d n e s s o f f i t t e s t s t a t i s t i c s o f h i g h-f r e q u e n c y s u b-b a n d sa n d c o m p a r e di t w i t h t h e d e n o i s i n g t h r e s h o l df o r d e n o i s i n g.F i n a l l y, s h e a r l e t i n v e r s et r a n s f o r m o nt h ec o e f f i c i e n t m a t r i x w a s p e r f o r m e dt o r e c o n s t r u c tt h e d e n o i s e d i m a g e s.T h e s i m u l a t i o ne x p e r i m e n t r e s u l t s s h o wt h a t t h i s a l g o r i t h mc a ne f f e c t i v e l y r e m o v eG a u s s i a n n o i s e i n r e m o t es e n s i n g i m a g e s,m a i n t a i nt h ee d g e t e x t u r e i n f o r m a t i o no f i m a g e s,a n da c h i e v eh i g h p e a ks i g n a l-t o-n o i s e r a t i ou n d e rd i f f e r e n tn o i s e l e v e l s,a m o n g w h i c ht h ea v e r a g e i n c r e a s e i s0.33dB c o m p a r e dw i t h t h e s h e a r l e t t h r e s h o l dd e n o i s i n g a l g o r i t h m.K e y w o r d s:r e m o t e s e n s i n g i m a g e;s h e a r l e t t r a n s f o r m;g o o d n e s s o f f i t t e s t;i m a g e d e n o i s i n g遥感图像目前已成为人们观察㊁分析世界的一种有效工具,在各领域均广泛应用[1-3].遥感图像在收稿日期:2022-05-01.第一作者简介:成丽波(1971 ),女,汉族,博士,教授,从事小波分析及遥感图像处理的研究,E-m a i l:c l b y y@126.c o m.通信作者简介:陈鹏宇(1999 ),男,汉族,硕士研究生,从事机器学习及遥感图像处理的研究,E-m a i l:919089661@q q.c o m.基金项目:国家自然科学基金(批准号:12171054)和吉林省教育厅科学技术研究项目(批准号:J J K H20230788K J).Copyright©博看网. All Rights Reserved.成像过程中会受随机噪声的干扰[4-6].为获取尽可能真实㊁清晰的遥感图像,满足工程实践的需要,如何在尽量不影响遥感图像原始信息的情况下有效对其进行去噪具有重要意义.遥感图像去噪的目的是抑制和消除随机噪声,保持图像的细节信息及良好的视觉效果.遥感图像的稀疏表示和去噪也是对其进一步处理的基础.经典的遥感图像去噪方法如小波变换[7-9]和复小波变换[10-11]是将噪声视为高频信号,过滤较大的变换系数[12-13],进而实现去噪,但高频系数中的图像细节部分通常会与噪声一起被去除[14].为克服上述问题,近年来,剪切波变换(s h e a r l e tt r a n s f o r m ,S T )[15-16]作为一种新的多尺度分析[17-18]技术已成为该领域研究的焦点.在频域中,剪切波变换是逐层细分的.因此,剪切波是性能较优的多维函数稀疏表示方法[19],使它在图像去噪领域有广泛的应用前景.剪切波变换具有良好的各向异性,可以很好地识别和分析图像的边缘和纹理信息[20].遥感图像含有大量的细节信息,因此,在遥感图像去噪中,剪切波是一种理想的工具.遥感图像噪声的主要来源是空气中的微小粒子,粒子对光的透射率及其空间分布均呈高斯分布[21].基于拟合优度(g o o d n e s so f f i t ,G O F )检验的图像去噪方法[22-23]利用高斯白噪声及变换系数的统计相关性估计去噪阈值,该方法能有效去除高斯噪声.相比于其他检验方法,A n d e r s o n -D a r l i n g (A D )检验[24]能在较小的样本条件下,保持稳健的检验性能.因此,本文使用A D 检验进行遥感图像的噪声识别.针对遥感图像去噪问题,本文提出一种基于剪切波和拟合优度检验的遥感图像去噪算法(s h e a r l e tt r a n s f o r ma n d g o o d n e s s o f f i t t e s t ,S T -G O F ),先使用剪切波变换分解含噪遥感图像,将高频系数进行拟合优度检验,再通过计算剪切波系数经验分布函数(e m p i r i c a l d i s t r i b u t i o n f u n c t i o n ,E D F )和模拟高斯噪声的累积分布函数(c u m u l a t i v ed i s t r i b u t i o n f u n c t i o n ,C D F )的统计距离获得A D 检验统计量,将统计量与阈值进行对比实现去噪.与其他相关算法进行仿真对比实验的结果表明了本文算法在视觉和性能上的优越性.1 剪切波变换二维图像函数f ɪL 2(ℝ2)的连续剪切波变换定义为S H ψf (a ,s ,t )=<f ,ψa ,s ,t >,(1)其中ψa ,s ,t (x )=d e t M a ,s1/2ψ(M -1a ,sx -t ),M a ,s =a a 1/2s 0æèçöø÷a ,a >0,s ɪℝ,t ɪℝ2.每个矩阵M 可分解为剪切矩阵S s =1s æèçöø÷01和各向异性矩阵A a =a 00a 1/æèçöø÷2的乘积,其中a >0且s ɪℝ,则连续剪切波可表示为ψa ,s ,t =a -3/4ψ(S s A a (x -t ))=a -3/4ψa -1-s/a 0a -1/æèçöø÷æèçöø÷2(x -t ),(2)其中ψa ,s ,t 是一个尺度㊁方向及位置参数分别为a ,s ,t 的局部性函数集合,此时a 称为尺度参数,s 为剪切参数,t 为平移参数.每个ψa ,s ,t 的频域支撑区间在一个梯形对内,该梯形对在不同的尺度参数a 下关于原点对称,其方向由剪切参数s 确定.通过减小尺度参数a ,剪切波能很好地捕捉图像的轮廓与边缘信息.每个剪切波ψa ,s ,t 在频域的支撑区间为(ξ1,ξ2):ξ1ɪ-2a ,-12éëêùûúa ɣ12a ,2éëêùûúa ,ξ2ξ1-s ɤa 1/{}2.(3)ψa ,s ,t 的频域支撑区间如图1所示.对于给定的一张不含噪声的图像P ,基于剪切波变换对加噪图像P n o i s y 的阈值去噪算法实现过程为P d e n o i s e =S T -1T σS T P n o i s y ,(4)其中S T 为剪切波正变换,T σ为阈值算子,S T -1为剪切波逆变换.剪切波正变换后获得图像的高频和低频系数,而噪声主要分布在高频系数中,因此需对高频系数选取合适的阈值算法进行处理,对处理后的全部系数进行剪切波逆变换得到去噪图像P d e n o i s e .8811 吉林大学学报(理学版) 第61卷Copyright ©博看网. All Rights Reserved.图1 剪切波水平(A )和垂直(B )频域支撑F i g .1 H o r i z o n t a l (A )a n d v e r t i c a l (B )f r e q u e n c y d o m a i n s u p po r t o f s h e a r l e t 2 拟合优度检验检验一个指定的样本与一套给定结果的吻合程度称为拟合优度检验,通常通过计算检验统计量并与临界值进行比较判断样本分布是否与假设分布一致.若样本分布与假设理论分布一致,即将观测值或数据符合指定模型或分布的情况称为零假设,把观测值拒绝指定模型或分布的情况称为替代假设.不同方法定义了不同的检验统计量作为量化观测值与指定分布下预期值之间差异的衡量标准.2.1 A n d e r s o n -D a r l i n g 检验设F (t )=ðt1(z >t )表示具有支持度t 的输入样本z 的经验分布函数,F r (t )=ʏtp (z >t )d z 表示对应于概率密度函数p (z )的假设累积分布函数.A n d e r s o n -D a r l i n g (A D )检验统计量定义为τA D =ʏɕ-ɕ(F r (t )-F (t ))2ψ(F r (t ))d F r (t ),(5)其中ψ(F r (t ))是为分布函数F r (t )的尾部赋予更多权重的加权函数,定义为ψ(F r (t ))=(F r (t )(1-F r (t )))-1,(6)其目的是提高检验时的灵活性.实际计算过程中式(5)统计距离的数值表达式可简写为τA D =-L -q ,(7)其中:L 表示给定观测值x t 的大小或在拟合优度检验时分块(窗口)的大小;q 定义为q =ðLt =1(2t -1)L[l n (F r (z t ))-l n (F r (z L +1-t ))].(8)在拟合优度检验框架内,错误拒绝候选分布的概率称为虚警概率P f a ,定义为P f a =P {τ>T H 0}=ʏ{z s .t .τ>λ}p (z H 0)d z ,(9)其中H 0表示与噪声检测相对应的零假设.文献[25-26]给出了虚警概率P f a 对应阈值的表格.2.2 本文算法本文利用剪切波变换对噪声图像的剪切波子带系数进行拟合优度检验,提出一种剪切波拟合优度检验去噪算法,在多尺度上将信号和噪声分离为单独的系数.令T G O F 表示阈值算子,则尺度j 处的阈值T j 为T j =T G O F (s ηj ,Pf a ),(10)其中s ηj 为输入标准差为σ的高斯白噪声(w h i t eG a u s s i a nn o i s e ,WG N )剪切波系数.通过计算s ηj的累积分布函数F r (t )进行拟合优度检验找到对应给定虚警概率的阈值T j .针对高斯噪声分布,F r (t )为F r (t )=ʏt-ɕ12πσe z2/σ2d z ,(11)取高频系数s ji 的中值估计噪声为9811 第5期成丽波,等:基于剪切波变换和拟合优度检验的遥感图像去噪 Copyright ©博看网. All Rights Reserved.0911吉林大学学报(理学版)第61卷^σ=m e d i a n(s j i)0.6745,(12)并通过s j i=s j i^σ,(13)对剪切波系数进行归一化.对归一化后的剪切波系数s j i执行拟合优度检验,以产生阈值系数^s,其中^s=GO F(s j i,T j).(14)本文算法的去噪二元假设表示如下:H0:τA DɤT j;i.e.s j iɪ噪声,H1:τA D>T j;i.e.s j iɪ信号,(15)其中H0和H1分别表示多尺度下的零假设和替代假设.通过给定阈值T j,首先利用式(7)计算高频系数的经验分布函数和模拟WG N系数的累积分布函数之间的统计距离,以获得检验统计量τA D,然后将其与阈值进行比较.由式(15)可见,噪声(即H0)样本的系数在变换矩阵中被识别出来,而剩余系数被保留为期望信号(即H1).若τA DɤT j,则该系数源于给定的累积分布函数F r(t)为噪声点.此时,将s j i 设置为0;否则,保留原始值,生成^s.最后,将系数矩阵逆归一化并通过剪切波逆变换恢复去噪图像.2.3算法流程本文提出的S T-G O F算法实现过程主要分为7个步骤.图2为S T-G O F算法的实现流程.图2S T-G O F算法流程F i g.2F l o wc h a r t o f S T-G O Fa l g o r i t h m算法1S T-G O F算法.步骤1)输入含噪遥感图像P n;步骤2)对P n进行剪切波变换P nѳS T(X)得到高频和低频系数矩阵;步骤3)模拟输入WG NѳS T(X),利用式(10),(12)估计T j,^σ,本文P f a=0.005;步骤4)利用式(13)对高频系数s j i进行归一化,计算检验统计量τA D:τA D=-L-ðL t=1(2t-1)L[l n(F r(s))-l n(F r(s L+1-n))];步骤5)若τA DɤT j,则s j iѳ0,否则,s j iѳ^s;步骤6)P dѳS T-1(^sˑ^σ);步骤7)输出去噪遥感图像P d.3仿真实验及结果分析为验证本文算法的有效性,选择灰度遥感图像进行主观和客观指标的评价.主观上,通过视觉观察遥感图像的边缘显示能力,并对比细节恢复程度等对本文算法的去噪效果进行评价;客观上,使用峰值信噪比(p e a ks i g n a l t on o i s e r a t i o,P S N R)衡量本文算法的去噪效果.3.1图像质量评价指标P S N R值作为图像去噪问题中的重要衡量指标,数值越大表示去噪程度越好.本文采用P S N RCopyright©博看网. All Rights Reserved.作为评价指标验证算法的优越性:P S N R =10l g25521X Y ðXi =1ðYj =1(p (i ,j )-d (i ,j ))2,(16)其中p (i ,j )表示原始无噪声图像,d (i ,j )表示去噪图像,X ,Y 表示图像尺寸.3.2 实验结果与分析本文从遥感数据集R S S C N 7中选取6张大小为400ˑ400的遥感图像F o r e s t ,P a r k i n g ,R e s i d e n t ,P l a y g r o u n d ,I n d u s t r y,M o u n t a i n s 作为仿真实验图像,如图3所示.将本文算法与离散小波拟合优度算法(DWT -G O F )㊁双树复小波拟合优度算法(D T C WT -G O F )㊁曲波阈值去噪算法(C u r v e l e t )㊁剪切波阈值去噪算法(S h e a r l e t )进行对比,5种算法均使用MA T L A BR 2018b 实现,计算平台为内存16G B的计算机,搭载处理器为I n t e l (R )C o r e (T M )i 7-8750H C P U@2.20G H z .图3 实验采用的无噪声遥感图像F i g .3 N o i s e -f r e e r e m o t e s e n s i n g i m a g e s u s e d i n e x pe r i m e n t 对比发现,添加了标准差约为20的高斯白噪声的遥感图像接近含有真实随机噪声的遥感图像.因此,对实验图像分别添加标准差为10,15,20,25的高斯白噪声.实验所得5种算法的P S N R 值分别列于表1~表4.由表1~表4可见,本文算法在不同噪声强度下的去噪图像P S N R 值均优于其他4种算法,其中与曲波阈值算法相比平均提高2.15d B ,与DWT -G O F 算法相比平均提高1.04d B ,与剪切波阈值算法相比平均提高0.33d B ,与D T C WT -G O F 算法相比平均提高0.23d B .表1 噪声标准差为10时不同算法的P S N R 值T a b l e 1 P S N Rv a l u e s o f d i f f e r e n t a l go r i t h m sw h e n s t a n d a r dd e v i a t i o no f n o i s e i s 10d B算法F o r e s tP a r k i n gR e s i d e n tP l a y g r o u n d I n d u s t r yM o u n t a i n s C u r v e l e t28.810327.428226.575829.509725.971727.3744DWT -G O F 30.229429.513128.603731.124428.185129.2467S h e a r l e t 30.446530.664429.648231.663729.355529.7056D T C WT -G O F 30.945930.761729.537531.868329.336929.9343S T -G O F 31.011730.870429.878731.923729.700730.2571图4为噪声标准差为25时,实验所得5种算法的去噪图像及局部细节放大结果.由图4可见,本文算法的去噪效果明显优于其他4种算法,能更好地保持遥感图像的边缘㊁轮廓和特征信息.在对局部细节进行放大时,本文算法去噪图像的纹理更细腻,视觉效果最优.1911 第5期成丽波,等:基于剪切波变换和拟合优度检验的遥感图像去噪 Copyright ©博看网. All Rights Reserved.2911吉林大学学报(理学版)第61卷表2噪声标准差为15时不同算法的P S N R值T a b l e2P S N Rv a l u e s o f d i f f e r e n t a l g o r i t h m sw h e n s t a n d a r dd e v i a t i o no f n o i s e i s15d B 算法F o r e s t P a r k i n g R e s i d e n t P l a y g r o u n d I n d u s t r y M o u n t a i n s C u r v e l e t27.349825.338224.819227.944724.163625.7178 DWT-G O F28.081127.113826.363829.022025.925626.9795 S h e a r l e t28.329628.178727.122729.549626.761127.2968 D T C WT-G O F28.646228.120127.105429.536326.877527.5937 S T-G O F28.778428.366827.380729.703327.101827.8700图4不同算法对遥感图像的去噪及细节效果F i g.4D e n o i s i n g a n dd e t a i l e f f e c t s o f d i f f e r e n t a l g o r i t h m s o n r e m o t e s e n s i n g i m a g e sCopyright©博看网. All Rights Reserved.表3 噪声标准差为20时不同算法的P S N R 值T a b l e 3 P S N Rv a l u e s o f d i f f e r e n t a l go r i t h m sw h e n s t a n d a r dd e v i a t i o no f n o i s e i s 20d B算法F o r e s tP a r k i n gR e s i d e n tP l a y g r o u n d I n d u s t r yM o u n t a i n s C u r v e l e t26.402524.202223.827527.020323.162724.6409DWT -G O F 26.826025.467124.842827.651824.408625.5641S h e a r l e t 27.053926.417025.447628.076425.025825.8947D T C WT -G O F 27.331126.379725.503428.153924.995326.2388S T -G O F 27.468426.649425.587728.266225.264026.3585表4 噪声标准差为25时不同算法的P S N R 值T a b l e 4 P S N Rv a l u e s o f d i f f e r e n t a l go r i t h m sw h e n s t a n d a r dd e v i a t i o no f n o i s e i s 25d B算法F o r e s tP a r k i n gR e s i d e n tP l a y g r o u n d I n d u s t r yM o u n t a i n s C u r v e l e t25.699123.383622.948926.290422.484523.9223DWT -G O F 25.871323.811823.310126.487222.796424.2511S h e a r l e t 26.228625.376724.272427.338624.034624.9171D T C WT -G O F 26.427625.176224.422727.218024.008225.1921S T -G O F 26.778025.568524.663327.457424.309825.4887综上所述,针对遥感图像中的高斯白噪声,本文提出了一种基于剪切波变换和拟合优度检验的遥感图像去噪算法,将遥感图像的去噪转换成了一个二元假设检验问题.首先,利用剪切波变换在频域分解含噪遥感图像,对分解得到的高频子带进行拟合优度检验;其次,通过对比检验统计量与去噪阈值分类系数,保留期望信号系数,实现去噪;最后,对系数矩阵逆归一化,通过剪切波逆变换进行图像重建.实验结果表明,与剪切波阈值去噪等算法相比,本文算法能显著提升峰值信噪比,有较强的细节留存和边缘保持能力,主观上也有较好的视觉效果.参考文献[1] 刘通,胡亮,王永军,等.基于卷积神经网络的卫星遥感图像拼接[J ].吉林大学学报(理学版),2022,253(1):99-108.(L I U T ,HU L ,WA N G YJ ,e ta l .S a t e l l i t eR e m o t eS e n s i n g I m a g e M o s a i cB a s e do nC o n v o l u t i o n a l N e u r a lN e t w o r k [J ].J o u r n a l o f J i l i nU n i v e r s i t y (S c i e n c eE d i t i o n ),2022,253(1):99-108.)[2] 夏英,黄秉坤.采用改进Y O L O v 3的高分辨率遥感图像目标检测[J ].重庆邮电大学学报(自然科学版),2022,34(3):383-392.(X I A Y ,HU A N GBK.O b j e c tD e t e c t i o no fH i g hR e s o l u t i o nR e m o t eS e n s i n g I m a ge sB a s e do n I m p r o v e dY O L O v 3[J ].J o u r n a lo fC h o n g q i n g U n i v e r s i t y o fP o s t sa n d T e l e c o mm u n i c a t i o n s (N a t u r a lS c i e n c e E d i t i o n ),2022,34(3):383-392.)[3] 宋志娜,李莎,杨建明,等.基于特征与区域定位增强的遥感舰船目标检测[J ].计算机工程,2023,49(8):257-264.(S O N GZ N ,L IS ,Y A N GJ M ,e t a l .R e m o t eS e n s i n g S h i p T a r g e tD e t e c t i o nB a s e do nF e a t u r ea n d R e g i o nL o c a l i z a t i o nE n h a n c e m e n t [J ].C o m p u t e rE n g i n e e r i n g ,2023,49(8):257-264.)[4] 徐华平,贾小宁.基于仿射不变块相似度量的B M 3D 图像去噪算法[J ].吉林大学学报(理学版),2022,253(1):109-118.(X U H P ,J I A X N.B M 3DI m a g eD e n o i s i n g A l g o r i t h m B a s e do n A f f i n eI n v a r i a n tP a t c hS i m i l a r i t y M e a s u r e [J ].J o u r n a l o f J i l i nU n i v e r s i t y (S c i e n c eE d i t i o n ),2022,253(1):109-118.)[5] 张意,阚子文,邵志敏,等.基于注意力机制和感知损失的遥感图像去噪[J ].四川大学学报(自然科学版),2021,58(4):45-55.(Z HA N G Y ,K A N Z W ,S HA O Z M ,e t a l .R e m o t eS e n s i n g I m a g eD e n o i s i n g B a s e do n A t t e n t i o n M e c h a n i s m a n dP e r c e p t u a lL o s s [J ].J o u r n a l o fS i c h u a n U n i v e r s i t y (N a t u r eS c i e n c eE d i t i o n ),2021,58(4):45-55.)[6] 周航,苏延池,李占山,等.基于子空间表示和加权低秩张量正则化的高光谱图像混合噪声去除方法[J ].吉林大学学报(理学版),2023,61(1):118-126.(Z HO U H ,S U Y C ,L IZS ,e t a l .M i x e dN o i s eR e m o v a lM e t h o df o rH y p e r s p e c t r a l I m ag e sB a s e do nS u b s p a c eR e p r e s e n t a t i o na n d W e i gh t e dL o w -R a n kT e n s o rR e g u l a ri z a t i o n [J ].J o u r n a l o f J i l i nU n i v e r s i t y (S c i e n c eE d i t i o n ),2023,61(1):118-126.)[7] F U C L ,L IP F ,S U IR L ,e ta l .H i g h -S p a t i a l -R e s o l u t i o n O F D R D i s t r i b u t e d T e m p e r a t u r eS e n s o rB a s e do n S t e p -b y -S t e p a n d I m a g eW a v e l e tD e n o i s i n g M e t h o d s [J ].S e n s o r s ,2022,22(24):9972-1-9972-8.3911 第5期成丽波,等:基于剪切波变换和拟合优度检验的遥感图像去噪Copyright ©博看网. All Rights Reserved.4911吉林大学学报(理学版)第61卷[8] WO N S E O K K,S O OHWA N Y,D O O C HU NS,e t a l.P u s h-B r o o m-T y p eV e r y H i g h-R e s o l u t i o nS a t e l l i t eS e n s o rD a t aC o r r e c t i o nU s i n g C o m b i n e d W a v e l e t-F o u r i e r a n d M u l t i s c a l eN o n-l o c a lM e a n sF i l t e r i n g[J].S e n s o r s,2015,15(9):22826-22853.[9] P A U L A,K U N D U A,C HA K I N,e ta l.W a v e l e t E n a b l e d C o n v o l u t i o n a l A u t o e n c o d e r B a s e d D e e p N e u r a lN e t w o r k f o rH y p e r s p e c t r a l I m a g eD e n o i s i n g[J].M u l t i m e d i aT o o l s a n dA p p l i c a t i o n s,2021,81(2):2529-2555.[10] L IM H,J I AZ H,Y A N GJ,e t a l.R e m o t eS e n s i n g I m a g eD e n o i s i n g B a s e do n t h eC o m b i n a t i o no f t h e I m p r o v e dB i S h r i n ka n dD TC WT[J].P r o c e d i aE n g i n e e r i n g,2011,24(6):470-474.[11] C H E N G Y,Q I A NSE.D e n o i s i n g o fH y p e r s p e c t r a l I m a g e r y U s i n g P r i n c i p a lC o m p o n e n tA n a l y s i sa n d W a v e l e tS h r i n k a g e[J].I E E ET r a n s a c t i o n s o nG e o s c i e n c e a n dR e m o t eS e n s i n g,2011,49(3):973-980.[12]J UMA H A A.D e n o i s i n g o fa nI m a g e U s i n g D i s c r e t eS t a t i o n a r y W a v e l e tT r a n s f o r m a n d V a r i o u sT h r e s h o l d i n gT e c h n i q u e s[J].J o u r n a l o f S i g n a l a n d I n f o r m a t i o nP r o c e s s i n g,2013,4(1):33-41.[13] A L E S S A N D R O G,F A B R I Z I O G,N I C O L A A,e ta l.A W a v e l e tF i l t e rC o m p a r i s o no n M u l t i p l eD a t a s e t sf o rS i g n a l C o m p r e s s i o na n dD e n o i s i n g[J].M u l t i d i m e n s i o n a l S y s t e m s a n dS i g n a l P r o c e s s i n g,2021,32(2):791-820.[14]郭林,孟旭东.基于偏微分方程与多尺度分析的图像去噪算法[J].吉林大学学报(理学版),2019,57(4):882-888.(G U O L,M E N G X D.I m a g e D e n o i s i n g A l g o r i t h m B a s e d o n P a r t i a l D i f f e r e n t i a l E q u a t i o n a n d M u l t i-s c a l eA n a l y s i s[J].J o u r n a l o f J i l i nU n i v e r s i t y(S c i e n c eE d i t i o n),2019,57(4):882-888.)[15] G U O K,L A B A T ED.O p t i m a l l y S p a r s eM u l t i-d i m e n s i o n a l R e p r e s e n t a t i o nU s i n g S h e a r l e t s[J].S I AMJ o u r n a l o nM a t h e m a t i c a lA n a l y s i s,2007,39(1):298-318.[16] K U T Y N I O K G,L A B A T ED.R e s o l u t i o no f t h eW a v e f r o n t S e tU s i n g C o n t i n u o u sS h e a r l e t s[J].T r a n s a c t i o n so ft h eA m e r i c a n M a t h e m a t i c a l S o c i e t y,2009,361(5):2719-2754.[17] B E N E D E T T OJ J,L I S.T h eT h e o r y o fM u l t i r e s o l u t i o nA n a l y s i sF r a m e sa n dA p p l i c a t i o n s t oF i l t e rB a n k s[J].A p p l i e da n dC o m p u t a t i o n a lH a r m o n i cA n a l y s i s,1998,5(4):389-427.[18]陈毅军,程浩,巩恩普,等.基于S h e a r l e t变换的尺度方向自适应阈值地震数据随机噪声压制方法[J].吉林大学学报(地球科学版),2021,51(4):1231-1242.(C H E N YJ,C H E N G H,G O N G EP,e t a l.R a n d o m N o i s e S u p p r e s s i o no f S e i s m i cD a t aw i t hS c a l e-O r i e n t e dA d a p t i v eT h r e s h o l dB a s e d o nS h e a r l e tT r a n s f o r m[J].J o u r n a l o f J i l i nU n i v e r s i t y(E a r t hS c i e n c eE d i t i o n),2021,51(4):1231-1242.)[19] L I M W Q.N o n s e p a r a b l e S h e a r l e t T r a n s f o r m[J].I E E E T r a n s a c t i o n s o n I m a g e P r o c e s s i n g,2013,22(5):2056-2065.[20] L Y UZY,Z HA N GCK,HA N M.D S T n e t:A N e w D i s c r e t eS h e a r l e tT r a n s f o r m-B a s e dC N N M o d e l f o r I m a g eD e n o i s i n g[J].M u l t i m e d i aS y s t e m s,2021,27(6):1-13.[21]梅树立,李晓飞,赵海英,等.基于剪切波和全变分的农田遥感图像去噪去伪影方法[J].农业工程学报,2017,33(增刊1):274-280.(M E I SL,L IXF,Z HA O H Y,e t a l.M e t h o do fD e n o i s i n g a n dR e m o v i n g A r t i f a c t s f o rF a r m R e m o t eS e n s i n g I m a g eB a s e do nS h e a r l e t a n dT o t a lV a r i a t i o n[J].T r a n s a c t i o n so f t h eC h i n e s eS o c i e t y o fA g r i c u l t u r a l E n g i n e e r i n g,2017,33(S u p p l1):274-280.)[22] N A V E E D K,S HA U K A TB,R E HMA N N U.S i g n a l D e n o i s i n g B a s e d o nD u a l T r e eC o m p l e xW a v e l e tT r a n s f o r ma n dG o o d n e s s o fF i tT e s t[C]//22n dI n t e r n a t i o n a lC o n f e r e n c eo nD i g i t a lS i g n a lP r o c e s s i n g(D S P).P i s c a t a w a y,N J:I E E E,2017:1-5.[23] N A V E E DK,S HA U K A TB,E H S A NS,e t a l.M u l t i-s c a l e I m a g eD e n o i s i n g U s i n g G o o d n e s s-o f-F i t T e s t B a s e d o nE D FS t a t i s t i c s[J].P L o SO n e,2019,14(5):e0216197-1-e0216197-25.[24] A N D E R S O N T W,D A R L I N G D A.A T e s to fG o o d n e s so fF i t[J].P u b l i c a t i o n so f t h eA m e r i c a nS t a t i s t i c a lA s s o c i a t i o n,1954,49:765-769.[25] S T E P H E N S M A.E D F S t a t i s t i c sf o r G o o d n e s so fF i ta n dS o m eC o m p a r i s o n s[J].J o u r n a lo ft h e A m e r i c a nS t a t i s t i c a lA s s o c i a t i o n,1974,69:730-737.[26] P E A R S O N E S,HA R T L E Y H O.B i o m e t r i k a T a b l e sf o rS t a t i s t i c i a n s[J].R o y a lS t a t i s t i c a lS o c i e t y,1967,130(2):257-258.(责任编辑:韩啸)Copyright©博看网. All Rights Reserved.。
电视画质芯片

电视画质芯片电视画质芯片是指用于提升电视机画面显示质量的特定芯片。
它主要用于增强图像细节、对比度、亮度、色彩等方面,以提供更清晰、更真实的观影体验。
以下将从电视画质芯片的原理、发展历程和优势等方面进行阐述。
一、电视画质芯片的原理电视画质芯片是通过对原始图像进行处理和优化,来改善电视画面质量的。
其原理主要包括以下几个方面:1. 图像去噪。
电视信号传输过程中会受到干扰,从而产生噪点。
电视画质芯片通过去除这些噪点,使画面更加清晰。
2. 锐化处理。
电视信号传输中会丢失一些细节信息,导致画面边缘不够清晰。
电视画质芯片通过增强轮廓和边缘的对比度,让画面更加锐利。
3. 色彩增强。
电视信号传输中,颜色饱和度可能会有所降低,导致画面色彩不真实。
电视画质芯片通过加强色彩饱和度,使画面更加鲜艳。
4. 对比度调整。
电视信号传输中,对比度可能会不够明显,导致画面细节不够丰富。
电视画质芯片通过调整亮度和对比度的比例,让画面更加清晰明亮。
二、电视画质芯片的发展历程随着科技的发展,电视画质芯片也在不断演进。
以下是电视画质芯片发展的几个重要阶段:1. 第一代电视画质芯片:最早的电视画质芯片主要是基于硬件电路实现图像处理,功能较为简单,只能进行一些基本的去噪和锐化处理。
2. 第二代电视画质芯片:随着数字信号处理技术的快速发展,第二代电视画质芯片开始采用数字信号处理器(DSP)实现,能够对图像进行更加精细的处理,同时增加了对色彩和对比度的调整功能。
3. 第三代电视画质芯片:随着高清、超高清电视技术的普及,第三代电视画质芯片开始支持更高的分辨率和更广的色彩空间。
同时,还加入了更高级的图像处理算法,如动态范围增强和局部对比度增强等。
4. 当前电视画质芯片:目前最新的电视画质芯片在提高画面质量的同时,还加入了一些智能功能,如自适应画质调节、智能场景识别等,以提供更加个性化的观影体验。
三、电视画质芯片的优势电视画质芯片的存在,极大地提升了电视机的画面质量,带来了许多优势:1. 提供更好的观影体验。
完整版基于DSP数字滤波器设计及仿真

2.1 系统功能介绍一个本质的应用系统中,总存在各种搅乱。
数字滤波器在语音信号办理、信号频谱估计、信号去噪、无线通信中的数字变频以及图像信号等各种信号办理中都有广泛的应用,数字滤波器也是使用最为广泛的信号办理算法之一。
在本设计中,使用MATLAB模拟产生合成信号,尔后利用CCS 进行滤波。
设定模拟信号的采样频率为48000Hz,。
设计一个 FIR 低通滤波器,其参数为:滤波器名称:FIR 低通滤波器采样频率:Fs=48000Hz通带截止频率:15000Hz阻带截止频率:16000Hz通带最大衰减:阻带最少衰减:80dB滤波器系数:由 MATLAB依照前述参数求得。
2.2 整体设计方案流程图使用 MATLAB 编写程序产使用MATLAB滤波器的设生待滤波数据计/解析工具 FDATOOL设计 FIR 滤波器使用 CCS编写程序实现滤波功能查察程序运行后的输入输出波形并比较滤波收效图 1 整体设计方案主要内容和步骤3.1 滤波器原理对于一个 FIR 滤波器系统,它的冲击响应总是又限长的,其系统函数可记为:N 1H zh n z n其中 Nn 01是 FIR 的滤波器的阶数, z n 为延时结, h n 为端口信号函数。
最基本的 FIR 滤波器可用下式表示:N 1y nh k x nk其中 x nk 0k 输入采样序列, h k 是滤波器系数, N 是滤波器的阶数 Y n 表示滤波器的输出序列,也可以用卷积来表示输出序列 y n 与 x n 、 h n 的关系,如下:y n x n h n3.2 操作步骤( 1)打开 FDATOOL , 依照滤波要求设置滤波器种类、通带截止频率、指定阶数、采样频率等。
指定完设计参数后单击按钮Design Filter ,生成滤波器系数。
( 2)把生成的滤波器系数传到目标 DSP 。
选择菜单 Targets->Export to CodeComposer Studio(tm)IDE ,打开 Export to C Header File 对话框,选择 C headerfile ,指定变量名 ( 滤波器阶数和系数向量 ) ,输出数据种类可选浮点型或 32 b ,16 b整型等,依照自己安装选择目标板板号和办理器号,单击OK ,保存该头文件,需指定文件名 (filtercoeff .h)和路径 ( 保存在 c : \ti\myprojects\fir 工程中 ) 。
数字图像去噪典型算法及matlab实现

图像去噪是数字图像处理中的重要环节与步骤。
去噪效果的好坏直接影响到后续的图像处理工作如图像分割、边缘检测等。
图像信号在产生、传输过程中都可能会受到噪声的污染,一般数字图像系统中的常见噪声主要有:高斯噪声(主要由阻性元器件内部产生)、椒盐噪声(主要是图像切割引起的黑图像上的白点噪声或光电转换过程中产生的泊松噪声)等;目前比较经典的图像去噪算法主要有以下三种:均值滤波算法:也称线性滤波,主要思想为邻域平均法,即用几个像素灰度的平均值来代替每个像素的灰度。
有效抑制加性噪声,但容易引起图像模糊,可以对其进行改进,主要避开对景物边缘的平滑处理。
中值滤波:基于排序统计理论的一种能有效抑制噪声的非线性平滑滤波信号处理技术。
中值滤波的特点即是首先确定一个以某个像素为中心点的邻域,一般为方形邻域,也可以为圆形、十字形等等,然后将邻域中各像素的灰度值排序,取其中间值作为中心像素灰度的新值,这里领域被称为窗口,当窗口移动时,利用中值滤波可以对图像进行平滑处理。
其算法简单,时间复杂度低,但其对点、线与尖顶多的图像不宜采用中值滤波。
很容易自适应化。
Wiener维纳滤波:使原始图像与其恢复图像之间的均方误差最小的复原方法,是一种自适应滤波器,根据局部方差来调整滤波器效果。
对于去除高斯噪声效果明显。
实验一:均值滤波对高斯噪声的效果代码I=imread('C:\Documents and Settings\Administrator\桌面\1.gif');%读取图像J=imnoise(I,'gaussian',0,0.005);%加入均值为0,方差为0.005的高斯噪声subplot(2,3,1);imshow(I);title('原始图像');subplot(2,3,2); imshow(J);title('加入高斯噪声之后的图像');%采用MATLAB中的函数filter2对受噪声干扰的图像进行均值滤波K1=filter2(fspecial('average',3),J)/255; %模板尺寸为3K2=filter2(fspecial('average',5),J)/255;% 模板尺寸为5K3=filter2(fspecial('average',7),J)/255; %模板尺寸为7K4= filter2(fspecial('average',9),J)/255; %模板尺寸为9 subplot(2,3,3);imshow(K1);title('改进后的图像1');subplot(2,3,4); imshow(K2);title('改进后的图像2');subplot(2,3,5);imshow(K3);title('改进后的图像3');subplot(2,3,6);imshow(K4);title('改进后的图像4');PS:filter2用法fspecial函数用于创建预定义的滤波算子,其语法格式为:h = fspecial(type)h = fspecial(type,parameters)参数type制定算子类型,parameters指定相应的参数,具体格式为:type='average',为均值滤波,参数为n,代表模版尺寸,用向量表示,默认值为[3,3]。
图像处理中的图像去噪方法与效果评估

图像处理中的图像去噪方法与效果评估图像去噪是数字图像处理中的一项关键任务,它旨在从图像中去除噪声,使其更清晰、更易于分析和理解。
在图像处理的众多应用中,图像去噪是一个必备的步骤,它可以用于医学图像、卫星图像、摄影图像等领域。
目前,有许多图像去噪方法可供选择,这些方法可以根据去噪原理、去噪效果和计算效率等方面进行分类。
下面将介绍几种常用的图像去噪方法,并对它们的效果进行评估。
1. 统计滤波方法统计滤波是一种基于统计原理的去噪方法,它通过对图像的像素值进行统计分析来判断噪声像素和信号像素,并通过滤波操作来抑制噪声。
常用的统计滤波方法包括中值滤波、高斯滤波和均值滤波。
中值滤波是一种简单有效的统计滤波方法,它通过对图像中的每个像素周围的邻域进行排序,然后取中间值作为该像素的新值。
中值滤波对于椒盐噪声和斑点噪声有较好的去除效果,但对于高斯噪声和高频噪声效果较差。
高斯滤波是一种基于高斯函数的滤波方法,它将像素的值与其周围像素的值进行加权平均,权值由高斯函数确定。
高斯滤波可以有效地平滑图像,并且保持边缘信息,但对于噪声的去除效果较差。
均值滤波是一种简单的滤波方法,它将像素的值与其邻域像素的平均值进行替换,可以有效地降低噪声的影响,但会导致图像模糊。
2. 小波变换方法小波变换是一种多尺度分析方法,可以将图像分解为不同频率的子带,然后根据子带的特征对噪声进行去除。
小波变换方法具有良好的去噪效果和较高的计算效率,在图像压缩、细节增强等应用中得到了广泛的应用。
小波去噪方法通常包括两个步骤:小波分解和阈值处理。
在小波分解阶段,图像被分解为不同频率的子带;在阈值处理阶段,对每个子带的系数进行阈值处理,然后通过逆小波变换将图像重建。
常用的小波去噪方法包括基于软阈值和硬阈值的去噪方法。
软阈值方法将小于某个阈值的系数置零,大于阈值的系数乘以一个缩放因子;硬阈值方法将小于阈值的系数置零,大于等于阈值的系数保持不变。
这两种方法在去除噪声的同时也会对图像细节造成一定的损失。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
I 摘 要 数字图像在产生、传输、处理、储存的过程中,不可避免地受到各类噪声的干扰导致信息难以获取,这就直接影响后期处理的效果。因此在对图像进行后续操作前必须进行提前加工处理,而图像去噪就是一种重要的方法之一。 图像噪声有很多种类,本文主要研究椒盐噪声和高斯噪声。中值滤波是一种非线性数字滤波器技术,中值滤波的基本原理是把数字图像或数字序列中一点的值用该点的一个邻域中各点值的中值代替,从而消除孤立的噪声点。 本文以灰度图像去噪为研究课题,通过使用DSP平台对图像分别加入椒盐噪声和高斯噪声,再进行中值滤波处理对比仿真结果,得出中值滤波能有效去除椒盐噪声的结论。也提出了一些中值滤波的优化改进算法思想。
关键词: 图像去噪;DSP;椒盐噪声;高斯噪声;中值滤波 目录 摘要......................................................错误!未定义书签。 1 设计方案................................................1 1.1 设计目的..............................................1 1.2 设计内容及要求........................................1 2 数字图像基础 .......................................1
2.1 数字图像 ............................................1 2.2 数字图像灰度化.......................................2 2.3 噪声的分类与特点.....................................3 2.4 灰度图像噪声的清除...................................7 2.5 图像去噪效果的评价方法..............................10 3 中值滤波...............................................11 3.1 标准中值滤波........................................11 3.2 中值滤波的改进算法..................................14 4 算法及DSP仿真.........................................18 4.1 算法................................................18 4.2 仿真................................................22 5 设计总结...............................................30 参考文献.................................................31 致谢.....................................................32 1
1 设计方案 1.1 设计目的 通过对课程设计任务的完成,使学生理解课程教学的理论内容,并且能够掌握和熟悉DSP的开发流程和基本的编程方法。同时,由于设计中涉及到各种器件的使用,可以提高学生综合运用各种技术和知识的能力。 1.2 设计内容及要求 基于DSP技术完成图像取反设计,具体要求如下: (1)总体方案设计。 (2)设计出软件编程方法,并写出源代码。 (3)仿真与结果分析。 (4)论文格式要符合学院的统一规定,结构要合符逻辑,表达要得体。 2 数字图像基础 2.1 数字图像 在人们的视界里,图像分为彩色和黑白。对于彩色图像,按照色度学理论:任何颜色都可由红、绿、蓝三种基本颜色按不同的比例混合得到。这样,自然界的图像可以用基于位置坐标的三维函数来表示,即:利用 f 表示空间坐标错误!未找到引用源。点的颜色函数,fred、fgreen、fblue分别表示该点颜色的红、绿、蓝三种原色的分量值。由于平面上每一点仅包括两个坐标(x, y),所以我们可以将空间三维函数转化为二维函数。对于黑白图像,就比彩色图像简单多了,只需要用错误!未找到引用源。表示该点图像的灰度(强度)即可。
我们所说的数字图像是相对于模拟图像而言的。是将图像按空间坐标和明暗程度的连续性分类得到的:
(1) 模拟图像指空间坐标和明暗程度都是连续变化的图像,计算机无法对其直接处理。即错误!未找到引用源。是空间的连续函数,错误!未找到引用源。为连续的空间,即在连续空间内,每一个点都有一个精确的值与之相对应。 2
(2) 数字图像是一种空间坐标和灰度均不连续的、用离散的数字表示的图像,这样的图像才能被计算机处理。
由于计算机仅能处理二进制数据,如果要用计算机来处理图像的话,必须把连续的图像函数转化为离散的数据集。经过离散化处理,原来连续的模拟图像就变成了由许多像素依据一定规则形成的数字图像来近似表示,一般用一个矩阵表示那些离散的数据,每一个离散的矩阵元素表示一个像素的颜色值。把像素依据不同的方式组织和存储,就形成了不同的图像格式,把图像数据存成文件就会得到相应的图像文件。
由于是数字图像,那么对于一幅黑白图像来说,只要把各个像素赋值为0或1即可,我们用1 表示白色,用 0 表示黑色,于是我们把一幅黑白图像称为二值图像,彩色图像或其它图像转化为黑白图像的过程叫做二值化。对于一幅彩色图像,每个像素我们都需要用3个取值范围为[0,255]之间的整数值来分别表示红、绿、蓝三原色分量,且这些分量都是用整型数据表示,称之为像素颜色的R, G, B值。表示一个取值范围为错误!未找到引用源。0,255]的整型数据,需要占用 8bit 空间,三个 R, G, B这样的整型数据就需要用24bit 来存储,所以,我们常把一幅真彩色位图称为 24 位位图。
在对数字图像进行处理的过程中,一般先把 24 位的位图转化为灰度图像,然后再进行相关的处理。所谓灰度图像,就是把图像矩阵中每个像素的值用1个错误!未找到引用源。0,255] 内的整型数据表示。如果彩色图像的R, G, B三个分量相等就成为灰度图像。灰度图中的每个像素只需要占有 8bit 空间,所以,灰度图像通常也被称为 8 位位图。
2.2 数字图像灰度化
对数字图像进行灰度化处理,主要方法有以下三种: 1.最大值法:使 R、G、B 的值等于三个值中最大的一个,即: R=G=B=max(R,G,B)
最大值法会使灰度图像亮度变高。 2.平均值法:使 R、G、B 的值等于三个的平均值,即: R=G=B=(R,G,B)/3 平均值法会形成比较柔和的灰度图像。 3
3.加权平均值法:根据重要性或敏感度等其他指标的要求给R,G,B 设定不同的权值,并使R、G、B值加权,即: 其中,错误!未找到引用源。, 错误!未找到引用源。, 错误!未找到引用源。分别为R, G, B的权值。错误!未找到引用源。, 错误!未找到引用源。, 错误!未找到引用源。取不同的值,加权平均值法将形成不同的灰度图像。人眼对绿色的敏感度最高,对红色的敏感度次之,对蓝色的敏感度最低,因此让错误!未找到引用源。可以得到合理的灰度图像。实验和理论推导都表明,当错误!未找到引用源。, 错误!未找到引用源。, 错误!未找到引用源。时,即:
此时,R, G, B, 的取值就是该像素的亮度值,此时得到的灰度图像最合理。 2.3 噪声的分类与特点
噪声可以理解为“妨碍人们感觉器官对所接收的信源信息理解的因素”。噪声通常是一种随机信号,不可能预测,因此它只能用概率统计的方法来分析。使用某种方法从被摄体或信息源把信息传递给受看者时,如果通过视觉接收平面二维亮度分布,那么对这种接收起干扰作用的亮度分布就叫图像噪声。噪声存在于图像的获取、传输和处理的各个环节。因此,去噪己经成为图像处理中极其必要的步骤之一,同时也是图像处理领域中一个重点研究的课题。本节先对噪声的来源和性质进行简要的介绍。
黑白图满足二维亮度分布错误!未找到引用源。函数,则噪声就是对亮度的干扰,用错误!未找到引用源。来表示。因为噪声具有随机性,所以我们需要用随机过程来描述,也就是说要求知道噪声的分布函数和密度函数。然而,在许多情况下,这些函数难以确定和描述,甚至根本无法得到,所以我们只能用统计方法来描述噪声,例如均值、方差、自相关函数等。
描述噪声的总功率:错误!未找到引用源。
噪声的交流功率利用方差描述:错误!未找到引用源。 噪声的直流功率利用均值的平方表示: 错误!未找到引用源。 4
一般在图像处理技术中常见的噪声有如下几种:白噪声、椒盐噪声、脉冲噪声、量化噪声等,我们重点介绍高斯噪声和椒盐噪声。
2.3.1 椒盐噪声 椒盐噪声(Pepper Noise) :是一种脉冲噪声,它在图像中产生黑色点和白色点。这种噪声在图像中表现非常明显,对图像具有严重的破坏作用,极大地妨碍了图像分割、边缘检测、特征提取等后续处理。它的噪声概率密度函数可表示为:
因为脉冲信号的强度通常比图像信号强度要大的多,因此脉冲噪声可以用灰度极限值(黑或白)来量化。我们一般假设a 和b 作为所取图像允许的最大灰度和最小灰度值,即“饱和”值。如果错误!未找到引用源。, 灰度 b在图像中表示为白点,而灰度 a 在图像中表示为黑点。如果错误!未找到引用源。或错误!未找到引用源。等于0,这时脉冲噪声被称之为单极性噪声。如果错误!未找到引用源。和错误!未找到引用源。都不等于0,但两者大小接近,这时脉冲噪声就像椒盐随机洒在图像上一样,因此,我们把双极性的脉冲噪声也叫做椒盐噪声。
2.3.2 高斯噪声 所谓高斯噪声是指它的概率密度函数服从高斯分布(即正态分布)的一类噪声。
高斯分布,也称正态分布,又称常态分布,记为错误!未找到引用源。,其中μ, 错误!未找到引用源。为分布的参数,分别为高斯分布的期望和方差。当有确定值时,p(x)也就确定了,特别当μ=0,错误!未找到引用源。=1时,X的分布为标准正态分布。
高斯随机变量z的PDF如下式: 错误!未找到引用源。