drdrawfig27.m

来自「beamlet变化的工具箱」· M 代码 · 共 111 行

function dRDrawFig27(OpenNewWindow,HaveTitle,LoadData)

%Test Wavelet Transform on HalfDome and Earlier Ridgelet versus New Ridgelet
% if you don't want to process it, set LoadData=1

if nargin<3
    LoadData=0;
end

if nargin<2
    HaveTitle=1;
end

if nargin<1
    OpenNewWindow=1;
end

if OpenNewWindow
    figure
end


%Test Wavelet Transform on HalfDome and Earlier Ridgelet versus New Ridgelet
%
% Reconstruction from 16 wavelet coefficients.
%

n=128;
if LoadData
    load 'dRDrawFig24.mat' ansel rptheta 
else
ansel = MakeHalfDomeImage(n,pi/8,1.3,3);%n,33, .77,3);
rptheta=FastOrthoRidgeletTransform(ansel);

end
% 2. Take FWT
s8filter = MakeONFilter('Symmlet',8);
wt = FWT2_PO(ansel,3,s8filter);

% 3. Take FRT
%rptheta=FastOrthoRidgeletTransform(ansel);
[theta,anomaly] = FRT_RP(ansel);

%Rearrenged coefficients
wrearr = 128* reverse(sort(abs(wt(:))));
rearr = reverse(sort(abs(theta(:))));
rprearr = reverse(sort(abs(rptheta(:))));  

% Show Partial Reconstruction
nwt = wt .* (128*abs(wt) >= wrearr(16));
anselw = IWT2_PO(nwt,3,s8filter);

ntheta = theta .* (abs(theta) >= rearr(16));
nrptheta = rptheta .* (abs(rptheta) >= rprearr(16));

if LoadData
    load 'dRDrawFig27.mat' anselk anselrp
else
anselk = IFRT_RP(ntheta,0.*anomaly); 
anselrp = abs(Inv_FastOrthoRidgeletTrans(nrptheta));
end

% Show Partial Reconstruction


subplot(2,2,1)
imagesc(abs(ansel))
if HaveTitle
title('HalfDome');
end
axis('image');colormap(1-hot);set (gca,'XTick',[]);
    set (gca,'YTick',[]);
subplot(2,2,2)
imagesc(abs(anselw)); axis('image'); 
if HaveTitle
    title('16 Wavelet Coefficients');
end
set (gca,'XTick',[]);
    set (gca,'YTick',[]);

subplot(2,2,3)
 imagesc(abs(anselk)); axis('image'); 
 if HaveTitle
title('16 Earlier Ridgelet Coefficients');
end
set (gca,'XTick',[]);
    set (gca,'YTick',[]);
subplot(2,2,4)
 imagesc(abs(anselrp)); axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);
if HaveTitle
    title('16 Ridgelet Coefficients');
end
%
%print -depsc RawFig35.eps

% Copyright (c) 2002 Ana Georgina Flesia
%
% Part of BeamLab Version:200
% Built:Friday,23-Aug-2002 00:00:00
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%
%
% Part of BeamLab Version:200
% Built:Saturday,14-Sep-2002 00:00:00
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%