小波变换软硬阈值半软阈值图像去噪matlab程序.doc

小波变换软硬阈值半软阈值图像去噪matlab程序
小波变换软硬阈值半软阈值图像去噪matlab程序
小波变换软硬阈值半软阈值图像去噪matlab程序
％软阈值 硬阈值　半软阈值 巴特沃斯滤波
ｃlｃ
cｌose　ａll
ｃlｅar aｌl
　
mａp=gｒay（２56）；
x=ｉmreａd('ｈzw。jpg＇）； 　
ｘ＝rｇb２ｇrａy(x）;
ｓubｐlot（２，3，１）;
iｍａｇe(ｘ）; 　　 　　　
coｌormap（ｍａp）; 　
tiｔle（’原始图片'）;
aｘｉs ｓquａｒｅ； 　
iniｔ=205５615866； 　
rａnｄn(’seｅｄ’，ｉnit）；
x１=5０.＊randn(ｓｉze(x））;　　　　%均值为0 方差50^2

x=ｄouｂｌe(x）
ｎx＝x+x1;
sｕｂｐlot（２,3，2）；
iｍａge(ｎx）；
ｃolｏｒmａp(map）；
ｔitｌｅ(’加噪后的图片＇）;
axis　sｑｕare;
c=nｕｍ2stｒ(ｃ)；
ｔｅxt(1０0，1０0,＇PSNR：’）；
tｅｘｔ（3０0，100，ｃ）;

％硬阈值
［ｔｈｒ，ｓoｒh，ｋeｅｐａpｐ］＝ｄdeｎcmｐ（'den'，'ｗｖ’，ｎx)；
ｎx１=wdｅnｃmp（’ｇｂl’,ｎx，'sym5＇,２，thr，'ｈ’,ｋeｅpaｐｐ)；
suｂｐloｔ（2,3，３);
ｉｍａgｅ（ｎx1）;
ｔitlｅ('Ó²ãÐֵȥÔëºóͼÏñ'）;
ａxis sqｕarｅ；
a1＝ｐsｎｒ（nx1，x）；
ａ１=ｎum２ｓtr(a1）;
teｘｔ（１00,10０，'ＰSNR：’）；
teｘt(300，10０,ａ1）；
％软阈值
小波变换软硬阈值半软阈值图像去噪matlab程序
小波变换软硬阈值半软阈值图像去噪matlab程序
小波变换软硬阈值半软阈值图像去噪matlab程序
ｎx２=wｄｅncmp(’gｂl’，nｘ，’ｓym5’,2，ｔhr，’ｓ'，ｋeepａpｐ);
sｕbｐlot（2，３，4）；
image（nx2）；
titｌe(’ÈíãÐֵȥÔëºóͼÏñ’）;
axiｓ squaｒｅ；
c=psｎr（nx2，x)；　
ｃ=ｎum２ｓtr(c)；
teｘt（10０,１00，＇PＳNR:＇);
ｔext（300,10０，c);

%半软阈值
ｎx3＝hsoｆt（nx,'ｓym５’,2,０．5，thr）；
ｓubplｏt(2，3，５);
iｍａgｅ(nx3);
ｔiｔlｅ（'°ëÈíãÐֵȥÔëºóͼƬ')；
aｘｉｓ sｑuarｅ;
a４＝psnr（nｘ3，x）; 　 　　　
a４=nｕm２ｓｔｒ（ａ４)；
text(10０,10０，’PSNR:');
ｔext（３0０,１00，a４）；

％巴特沃斯
ｇ＝fft２（nx）;
g＝ｆｆtshｉft（g)；　
［ｍ,n］=sｉｚe（g)；
N＝3;ｄ0=６0；
ｎ1＝floor(m/2）；
n2＝ｆｌoｏr(ｎ／２)；
for i=1：m
　 ｆor ｊ=1:n