xiaobo_edge.m

关于小波变换的代码

 function xiaobo_edge(infile)

global myfigure

I = imread(infile);
[row,col]=size(I);

mybar=waitbar(0,'请稍等,程序正在运行,去喝杯咖啡再来....');

wname='bior4.4'
%

% subplot('Position',[0.1 0.4 0.3 0.3]);
% imshow(I);
% title('原图');
%
[ca1,ch1,cv1,cd1] = dwt2(I,wname);
waitbar(0.1);
[ca2,ch2,cv2,cd2] = dwt2(ca1,wname);
waitbar(0.2);
[ca3,ch3,cv3,cd3] = dwt2(ca2,wname);
thr=4;a0=1;n=7;
[edge_mf1,grads1] = local_max_mode(cv1,ch1,thr,a0,n);
thr=15;a0=1;n=5;
waitbar(0.3);
[edge_mf2,grads2] = local_max_mode(cv2,ch2,thr,a0,n);
thr=30;a0=1;n=3;
[edge_mf3,grads3] = local_max_mode(cv3,ch3,thr,a0,n);
% 
waitbar(0.4);
edge_mf1=guiyi(edge_mf1);
edge_mf2=guiyi(edge_mf2);
edge_mf3=guiyi(edge_mf3);

waitbar(0.5);
%
ca3=wthresh(ca3,'h',3000000);
ca2=wthresh(ca2,'h',3000000);
ca1=wthresh(ca1,'h',3000000);
ca2_ca3=idwt2(ca3,edge_mf3,edge_mf3,edge_mf3,wname,size(ca2));
ca2_ca3=guiyi(abs(ca2_ca3));
waitbar(0.6);
com_2=mul_c(ca2_ca3,edge_mf2);
ca1_ca3=idwt2(ca2,com_2,com_2,com_2,wname,size(ca1));
ca1_ca3=guiyi(abs(ca1_ca3));
waitbar(0.7);
com_1=mul_c(ca1_ca3,edge_mf1);
n=5;a0=0.01;
com_1=ada_thr(com_1,n,a0);
waitbar(0.8);
com_1=wthresh(com_1,'h',0.00001);
waitbar(0.9);
com_1=dayu_c(com_1,0);
waitbar(1);

close(mybar);

subplot('Position',[0.1 0.4 0.3 0.3]);
imshow(I);
title('原图');
subplot('Position',[0.6 0.4 0.3 0.3])
imshow(com_1*255);
title('多尺度');
subplot('Position',[0.05 0.05 0.3 0.3])
imshow(edge_mf3*255);
title('2层');
subplot('Position',[0.35 0.05 0.3 0.3])
imshow(edge_mf2*255);
title('3层');
subplot('Position',[0.65 0.05 0.3 0.3])
imshow(edge_mf1*255);
title('4层');