drdrawfig25.m

beamlet变化的工具箱

function dRDrawFig01(OpenNewWindow,HaveTitle,LoadData)

% 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


%DrawFig25: Compare reconstruction from 8 ,32,128,512 Ortho Ridgelet Z_p and coefficients  on "object" See Do-Veterli(2001)

% % 1. Make an image
n=256;
ansel = double(imread('object.tif'));
anselp = ansel(1:255,1:255);

% 2. Take FRT
zp2theta = FRT_Zp2(anselp);
if LoadData
    load 'dRDrawFig25.mat' rptheta
else
    rptheta=FastOrthoRidgeletTransform(ansel);

end

%
 %Rearrenged coefficients 
 % Note Normalization;
 
zp2rearr = reverse(sort(abs(zp2theta(:)))) * sqrt(255);
rprearr = reverse(sort(abs(rptheta(:))))./256;  

nzp2theta128 = zp2theta .* (abs(zp2theta)*sqrt(255) >= zp2rearr(128));
anselzp2128 = InvFRT_Zp2(nzp2theta128); 
nrptheta128 = rptheta .* (abs(rptheta)./256 >= rprearr(128));


nzp2theta256 = zp2theta .* (abs(zp2theta)*sqrt(255) >= zp2rearr(256));
anselzp2256 = InvFRT_Zp2(nzp2theta256); 
nrptheta256 = rptheta .* (abs(rptheta)./256 >= rprearr(256));

if LoadData
    load 'dRDrawFig25.mat' anselrp128 anselrp256
else
anselrp128=abs(Inv_FastOrthoRidgeletTrans(nrptheta128));

anselrp256=abs(Inv_FastOrthoRidgeletTrans(nrptheta256));

end


subplot(2,2,1);
imagesc(abs(anselzp2128)); axis('square'); colormap(1-hot);
    set (gca,'XTick',[]);
    set (gca,'YTick',[]);
    if HaveTitle
title('Reconstructed from 128 Zp2 Coefficients'); 
end


subplot(2,2,2);
imagesc(abs(anselzp2256)); axis('square'); set (gca,'XTick',[]);
    set (gca,'YTick',[]);
    if HaveTitle
        title('Reconstructed from 256 Zp2 Coefficients');
    end

subplot(2,2,3);
imagesc(abs(anselrp128)); axis('square');set (gca,'XTick',[]);
    set (gca,'YTick',[]);
    if HaveTitle
title('Reconstructed from 128 RP Coefficients'); 
end

subplot(2,2,4);
imagesc(abs(anselrp256)); axis('square'); set (gca,'XTick',[]);
    set (gca,'YTick',[]);
    if HaveTitle
title('Reconstructed from 256 RP Coefficients'); 
end

% print -depsc RawFig25.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]
%