bayes门限去噪主程序及其它相关程序Word格式.docx
- 文档编号:15710564
- 上传时间:2022-11-15
- 格式:DOCX
- 页数:88
- 大小:32.41KB
bayes门限去噪主程序及其它相关程序Word格式.docx
《bayes门限去噪主程序及其它相关程序Word格式.docx》由会员分享,可在线阅读,更多相关《bayes门限去噪主程序及其它相关程序Word格式.docx(88页珍藏版)》请在冰豆网上搜索。
imshow(uint8(t2))
subplot(223)
imshow(uint8(t3))
subplot(224)
imshow(uint8(t4))
t2=bayes_thr(t2);
t3=bayes_thr(t3);
t4=bayes_thr(t4);
%-------------------level2--------------------
[t1_1,t1_2,t1_3,t1_4]=dwt2(t1,wave_name);
t1_2=bayes_thr(t1_2);
t1_3=bayes_thr(t1_3);
t1_4=bayes_thr(t1_4);
%--------------level3---------------------------
[t3_1,t3_2,t3_3,t3_4]=dwt2(t1_1,wave_name);
t3_2=bayes_thr(t3_2);
t3_3=bayes_thr(t3_3);
t3_4=bayes_thr(t3_4);
%--------------level4---------------------------
[t4_1,t4_2,t4_3,t4_4]=dwt2(t3_1,wave_name);
t4_2=bayes_thr(t4_2);
t4_3=bayes_thr(t4_3);
t4_4=bayes_thr(t4_4);
%--------------level5---------------------------
[t5_1,t5_2,t5_3,t5_4]=dwt2(t4_1,wave_name);
%-------------------------------------------------
temp8=idwt2(t5_1,t5_2,t5_3,t5_4,wave_name);
temp9=ones(size(t4_2));
[dep1,wide1]=size(temp8);
[dep2,wide2]=size(t4_2);
for(i=1:
dep2)
for(j=1:
wide2)
temp9(i,j)=temp8(i,j);
end
end
temp6=idwt2(temp9,t4_2,t4_3,t4_4,wave_name);
temp7=ones(size(t3_2));
[dep1,wide1]=size(temp6);
[dep2,wide2]=size(t3_2);
temp7(i,j)=temp6(i,j);
temp4=idwt2(temp7,t3_2,t3_3,t3_4,wave_name);
temp5=ones(size(t1_2));
[dep1,wide1]=size(temp4);
[dep2,wide2]=size(t1_2);
temp5(i,j)=temp4(i,j);
temp2=ones(size(t2));
temp1=idwt2(temp5,t1_2,t1_3,t1_4,wave_name);
[dep1,wide1]=size(temp1);
[dep2,wide2]=size(t2);
temp2(i,j)=temp1(i,j);
temp=idwt2(temp2,t2,t3,t4,wave_name);
%[dep1,wide1]=size(temp1);
%[dep2,wide2]=size(t2);
%if(dep1>
imshow(uint8(temp));
imwrite(uint8(temp),'
f_bayes.bmp'
'
bmp'
)
%------------子程序bayes_thr------------------------
functiont2=bayes_thr(t2);
[dep,wide]=size(t2);
N=dep*wide;
sw=(median(median(abs(t2)))/0.6745).^2
sx=0;
fori=1:
dep
forj=1:
wide
sx=sx+t2(i,j).^2;
sx=sx/N
sx=max(0.0001,(sx-sw))
th_h=1.4142*sw/sqrt(sx)
1:
if(abs(t2(i,j))<
th_h)
t2(i,j)=0;
else
t2(i,j)=sign(t2(i,j))*(abs(t2(i,j))-th_h);
%------------------GCV阈值主程序------------------------------------------
wavename='
n_lena_0.005.bmp'
[dep,wide]=size(a);
r=0.5;
T=40;
%--------dep1--------------------------------------------
[t,t1,t2,t3]=dwt2(a,wavename);
%--------dep2--------------------------------------------
[u,u1,u2,u3]=dwt2(t,wavename);
%--------dep3--------------------------------------------
[v,v1,v2,v3]=dwt2(u,wavename);
%--------dep4--------------------------------------------
[w,w1,w2,w3]=dwt2(v,wavename);
%--------dep5--------------------------------------------
[x,x1,x2,x3]=dwt2(w,wavename);
%--------denoisingprocessing---------------------------
%----------step1-----------------------
t1_1=soft_thr_GCV(t1,r);
t2_1=soft_thr_GCV(t2,r);
t3_1=soft_thr_GCV(t3,r);
%----------step2-----------------------
u1_1=soft_thr_GCV(u1,r);
u2_1=soft_thr_GCV(u2,r);
u3_1=soft_thr_GCV(u3,r);
%----------step3-----------------------
v1_1=soft_thr_GCV(v1,r);
v2_1=soft_thr_GCV(v2,r);
v3_1=soft_thr_GCV(v3,r);
%----------step4-----------------------
w1_1=soft_thr_GCV(w1,r);
w2_1=soft_thr_GCV(w2,r);
w3_1=soft_thr_GCV(w3,r);
%-----------------syn1--------------------------------------
syn1=idwt2(w,w1_1,w2_1,w3_1,wavename);
syn1_1=ones(size(v1));
[dep2,wide2]=size(v1);
syn1_1(i,j)=syn1(i,j);
%-----------------syn2--------------------------------------
syn2=idwt2(syn1_1,v1_1,v2_1,v3_1,wavename);
syn2_1=ones(size(u1));
[dep2,wide2]=size(u1);
syn2_1(i,j)=syn2(i,j);
%-----------------syn3--------------------------------------
syn3=idwt2(syn2_1,u1_1,u2_1,u3_1,wavename);
syn3_1=ones(size(t1));
[dep2,wide2]=size(t1);
syn3_1(i,j)=syn3(i,j);
end
%-----------------syn4--------------------------------------
new_image1=idwt2(syn3_1,t1_1,t2_1,t3_1,wavename);
new_image=ones(size(a));
[dep2,wide2]=size(a);
new_image(i,j)=new_image1(i,j);
figure
imshow(uint8(a))
imshow(uint8(new_image))
imwrite(uint8(new_image),'
GVC.bmp'
)
%-----------------子函数soft_thr_GCV-------------------------------------
functiont_GCV=soft_thr_GCV(t2,r);
th_h=r*(median(median(abs(t2),1))/0.6745)*sqrt(2*
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- bayes 门限 主程序 及其 相关 程序