drdrawfig30.m

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

function dRDrawFig30(OpenNewWindow,HaveTitle,LoadData)

%Test ORT on BandPassed classical Half Dome with n=256; layout of coeff.,plot of rearrenged coeff. and errors, recons. from 100, 50 and 16 coeff.
% 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 ORT on BandPassed classical Half Dome with n=256; layout of coeff.,plot of rearrenged coeff. and errors, recons. from 100, 50 and 16 coeff.
%n=256; Original, takes a lot of time!!

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

end


% 1. Make an image
 %  image = MakeHalfDomeImage(n,pi/8,1.3,3);
   
   
   %1.2Bandpass it
    
   	domain = (-(n/2)):(n/2-1);
	[x,y] = meshgrid(domain,domain);
%
		r = sqrt(x.^2 + y.^2);
        j=5;
		sr = r ./ (2^j);
%

Window=BandWindow(sr,3,'bandpass');
Fimage = fftshift(fft2(image));
NewFimage =  Fimage .* Window;
NewImg= abs(ifft2(fftshift(NewFimage)));



% 2. Take FRT
theta = FastOrthoRidgeletTransform(NewImg);

%3. Rearrange coefficients
rearr = reverse(sort(abs(theta(:))))./n;

% Show coefficients

subplot(1,2,1);
plot(rearr,'k','LineWidth',2); 
if HaveTitle
    title('(a): Rearranged Ridgelet Coefficients');
    ylabel('Rearranged Coefficient');
    xlabel('Rank order, m')
end
xlim([0 50]); 
axis('square');
% Show Approximation Numbers
subplot(1,2,2); plot(reverse(cumsum(reverse(rearr.^2))),'k','LineWidth',2); 
if HaveTitle
    title('(b): m-term Approximation Errors')
	ylabel('Squared Error');
	xlabel('Number of Terms, m')
end
xlim([0 50]); 
axis('square');
 %print -depsc RawFigHD2.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]
%