drdrawfig24.m

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

function dRDrawFig24(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
%DrawFig24: Compare reconstruction from 50 and 100 Ortho Ridgelet Z_p and coefficients  on classical HalfDome

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

end

anselp = ansel(1:127,1:127);

% 2. Take FRT
zp2theta = FRT_Zp2(anselp);

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


nzp2theta50 = zp2theta .* (abs(zp2theta)*sqrt(127) >= zp2rearr(50));
anselzp250 = InvFRT_Zp2(nzp2theta50); 

nzp2theta100 = zp2theta .* (abs(zp2theta)*sqrt(127) >= zp2rearr(100));
anselzp2100 = InvFRT_Zp2(nzp2theta100); 

nrptheta100 = rptheta .* (abs(rptheta)./128 >= rprearr(100));
nrptheta50 = rptheta .* (abs(rptheta)./128 >= rprearr(50));

if LoadData
    load 'dRDrawFig24.mat'  anselrp50 anselrp100
else
anselrp50 = abs(Inv_FastOrthoRidgeletTrans(nrptheta50));
anselrp100 = abs(Inv_FastOrthoRidgeletTrans(nrptheta100));
end



subplot(2,2,1);
imagesc(abs(anselzp250)); axis('square'); colormap(1-hot);
    set (gca,'XTick',[]);
    set (gca,'YTick',[]);
if HaveTitle
    title('Reconstructed from 50 Zp2 Coefficients'); 
end
subplot(2,2,2);
imagesc(abs(anselrp50)); axis('square'); 
    set (gca,'XTick',[]);
    set (gca,'YTick',[]);
    if HaveTitle
title('Reconstructed from 50 RP Coefficients'); 
end
subplot(2,2,3);
imagesc(abs(anselzp2100)); axis('square'); 
    set (gca,'XTick',[]);
    set (gca,'YTick',[]);
    if HaveTitle
title('Reconstructed from 100 Zp2 Coefficients'); 
end
subplot(2,2,4);
imagesc(abs(anselrp100)); axis('square');
    set (gca,'XTick',[]);
    set (gca,'YTick',[]);
    if HaveTitle
title('Reconstructed from 100 RP Coefficients'); 
end
 %print -deps RawFig24.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]
%