caise.rar_matlab可视化_可视化_图像可视化_图像彩色化_彩色化

%Test for dual tree 1D %Test of perfect reconstruction %***************************************************** %* author: Caroline Chaux * %* institution: Institut Gaspard Monge - CNRS * %* Universit? de Marne la Vall�e * %* date: Thursday, November, 2nd 2006 * %***************************************************** clear all clc %description disp('This program decomposes an image via a 1D dual-tree M-band transform and reconstructs the obtained coeffients.'); %number of bands of the filter bank M=2; %decomposition levels Jmin=4; rm=4; %signal to decompose x=load_image('lena',256); x=x(1,:); % % %Size adaptation % sigma = 0.03 * (max(x)-min(x)); % noise level % xn = x + sigma*randn(1,256); % noisy signal y = perform_wavelet_transform(x, Jmin,+1); % T = 80; % % perform the thresholding % y =y.* (abs(y)>T); z1 = perform_wavelet_transform(y, Jmin,-1); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% x=x(1:M^rm*floor(length(x)/M^rm)); %decomposition [wp,wm,H,HH,G,Gh]=dualtreeM1D_freq(x,M,rm); % wp=[wp{1}{2} wp{2}{2} wp{3}{2} wp{4}{2} wp{4}{1}]; % wp=[wp{4}{1} wp{4}{2} wp{3}{2} wp{2}{2} wp{1}{2}]; % wm=[wm{4}{1} wm{4}{2} wm{3}{2} wm{2}{2} wm{1}{2}]; % wp =wp.* (abs(wp)>T); % wm =wm.* (abs(wm)>T); % wp1{4}{1}=wp(1:16);wp1{4}{2}=wp(17:32);wp1{3}{2}=wp(33:64);wp1{2}{2}=wp(65:128);wp1{1}{2}=wp(129:256); % wm1{4}{1}=wm(1:16);wm1{4}{2}=wm(17:32);wm1{3}{2}=wm(33:64);wm1{2}{2}=wm(65:128);wm1{1}{2}=wm(129:256); z2= idualtreeM1D_freq(wp,wm,rm,H,HH,G,Gh); M=4; [wp,wm,H,HH,G,Gh]=dualtreeM1D_freq(x,M,rm); % % wp1=[wp{1}{4} wp{1}{3} wp{1}{2} wp{2}{4} wp{2}{3} wp{2}{2} wp{3}{4} wp{3}{2} wp{3}{4} wp{4}{4} wp{4}{3} wp{4}{2} wp{4}{1}]; % wp=[wp{4}{1} wp{4}{2} wp{4}{3} wp{4}{4} wp{3}{2} wp{3}{3} wp{3}{4} wp{2}{2} wp{2}{3} wp{2}{4} wp{1}{2} wp{1}{3} wp{1}{4}]; % wm=[wm{4}{1} wm{4}{2} wm{4}{3} wm{4}{4} wm{3}{2} wm{3}{3} wm{3}{4} wm{2}{2} wm{2}{3} wm{2}{4} wm{1}{2} wm{1}{3} wm{1}{4}]; % wp =wp.* (abs(wp)>T); % wm =wm.* (abs(wm)>T); % wp2{4}{1}=wp(1);wp2{4}{2}=wp(2);wp2{4}{3}=wp(3);wp2{4}{4}=wp(4); % wp2{3}{2}=wp(5:8);wp2{3}{3}=wp(9:12);wp2{3}{4}=wp(13:16); % wp2{2}{2}=wp(17:32);wp2{2}{3}=wp(33:48);wp2{2}{4}=wp(49:64); % wp2{1}{2}=wp(65:128);wp2{1}{3}=wp(129:192);wp2{1}{4}=wp(193:256); % wm2{4}{1}=wm(1);wm2{4}{2}=wm(2);wm2{4}{3}=wm(3);wm2{4}{4}=wm(4); % wm2{3}{2}=wm(5:8);wm2{3}{3}=wm(9:12);wm2{3}{4}=wm(13:16); % wm2{2}{2}=wm(17:32);wm2{2}{3}=wm(33:48);wm2{2}{4}=wm(49:64); % wm2{1}{2}=wm(65:128);wm2{1}{3}=wm(129:192);wm2{1}{4}=wm(193:256); z3= idualtreeM1D_freq(wp,wm,rm,H,HH,G,Gh); subplot(141) plot(x) title('zaoshengtu') xlabel('(a)') subplot(142) plot(z1) title('wav') xlabel('(b)') subplot(143) plot(z2) title('2 bands') xlabel('(c)') subplot(144) plot(z3) title('4 bands') xlabel('(d)') p1=psnr(x,z1) p2=psnr(x,z2) p3=psnr(x,z3)