drdrawfig12.m

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

function dRDrawFig12(OpenNewWindow,HaveTitle, LoadData)

%Figure 12a: Plot Ridgelet tiling

if nargin<2
    HaveTitle=1;
end
if nargin<1
    OpenNewWindow=1;
end
if OpenNewWindow
    figure
end

LoadData=0;

%figure; 
subplot(1,2,1)
hold on;
J=6;
color='b';
for j=0:J,
	rj = 2^(-j);
	
	ij = (J-j);
	nij = 2^(ij)*2;
	for l=0:(nij-1),
		thetal = 2*pi*l ./ nij;
		poutx = rj * cos(thetal); pouty = rj * sin(thetal);
		pinx  = poutx/2; piny = pouty/2;
		plot([poutx pinx],[pouty piny],color)
	end
    plotcircle(rj);
end
axis('square')
raxp = axis;
axis('off')
if HaveTitle
       title('Ridgelet Tiling')
    end
hold off
%print -deps RawFig12a.eps


% Plot digital ridgelet tiling
	
    deep    = [5 6];
	b_ridge = MakeRPRidgeletBasis(deep);
	
    %figure;
    subplot(1,2,2)
    raxp = axis;
   PlotRectoPolarRPPartition(b_ridge,'k',raxp,deep);axis off;
   if HaveTitle
        title('Digital Ridgelet Tiling')
    end

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