drpdrawfig05.m

beamlet变化的工具箱

function dRPDrawFig05(OpenNewWindow,HaveTitle,LoadData)

% TestRPDisc, find best partition


if nargin<3
    LoadData=0;
end

if nargin<2
    HaveTitle=1;
end

if nargin<1
    OpenNewWindow=1;
end

if OpenNewWindow
    figure
end


LoadData=0;

% TestRPDisc

	img = MakeDiscImage(128,30,15,15,3);
%
	
    subplot(1,2,1)
	imagesc(img)
	axis image;set (gca,'XTick',[]);
   set (gca,'YTick',[]);
   
   if HaveTitle
       title('Disc');
   end

    MaxDeep = [4 4];
    %Entropy = 'N(eps)';
    %par = 1;%2%1.5;

    img0  = img - mean(mean(img));
    RPFT  = fft2_rp(img0);
	CFT = abs(fftshift(fft2(img0)));
    RPPkt = CalcRPPktTable(RPFT,MaxDeep,'Sine');
	RPtree  = CalcRPStatTree(RPPkt,MaxDeep,'l^p',1)%Entropy,par);%'l^p',1);
	bb    = BestRPBasis(RPtree,MaxDeep); 
	subplot(1,2,2)
    imagesc(log(1 + CFT));
 	Cartax = axis; 
 	hold on;
	PlotRectoPolarRPPartition(bb,'k',Cartax,MaxDeep); 
    if HaveTitle
        title('Best RectoPolar Partition of FFT of Disc');
    end
    axis off
    colormap(1-hot)
	hold off;

   
% Copyright (c) 2002 David Donoho and 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]
%