]> AND Private Git Repository - these_gilles.git/blob - THESE/codes/wave/allcode/den3.m
Logo AND Algorithmique Numérique Distribuée

Private GIT Repository
v1.2 19 décembre
[these_gilles.git] / THESE / codes / wave / allcode / den3.m
1
2
3 s1 = double(imread('st.tif'));
4 s = s1(:,:,3);
5 x = s + 20*randn(size(s));
6 t = 0:5:50;         % threshold range(0~50),increment by 5
7 e = den2(s,x,t);    % using standard method
8 re = rden2(s,x,t);  % using reduced 2-D dual-tree method 
9 ce = cden2(s,x,t);  % using complex 2-D dual-tree method
10
11 figure(1)
12 plot(t,e)
13 hold on
14 plot(t,re,'r')
15 hold on
16 plot(t,ce,'k')
17 title('RMS error V.S. Threshold Pt.') 
18 xlabel('Threshold pt.');
19 ylabel('RMS error');
20 legend('Standard 2D','Reduced 2D dual','Complx 2D dual',0);
21
22
23 [emin,k] = min(e);
24 T = t(k);
25 y = denS2D(x,T);
26 figure(2)
27 clf
28 image(y/4)
29 colormap(gray)
30 axis image
31 axis off
32 print -deps den3_A
33
34 ax = [280*[1 1] 250*[1 1]]+[0 1 0 1]*120;
35
36 figure(3)
37 clf
38 image(y/4)
39 colormap(gray)
40 axis(ax)
41 axis square
42 axis off
43 print -deps den3_B
44
45 [emin,k] = min(ce);
46 T = t(k);
47 y = denC2D(x,T);
48 figure(4)
49 clf
50 image(y/4)
51 colormap(gray)
52 axis image
53 axis off
54 print -deps cden3_A
55
56 figure(5)
57 clf
58 image(y/4)
59 colormap(gray)
60 axis(ax)
61 axis square
62 axis off
63 print -deps cden3_B
64