drdrawfig06.m

beamlet变化的工具箱

function dRDrawFig06(OpenNewWindow,HaveTitle,LoadData)
% DrawFig6: Describe layout of ridgelet and ortho-ridgelet coefficients

if nargin<3
    LoadData=0;
end

if nargin<2
    HaveTitle=1;
end

if nargin<1
    OpenNewWindow=1;
end

if OpenNewWindow
    figure
end

n=128; 

% 3. DrawRidgeGrid: Depict subband layout of ridgelet coefficients

%figure;
subplot(1,2,1)
plot([0 (2*n) (2*n) 0 0],[0 0 (2*n) (2*n) 0], '-r'); hold on
plot([n n],[0 (2*n)],'-r');

J = round(log2(n));

%horizontal lines
for js=3:(J-1),
	plot([0 n],[(2^(js+1)) (2^(js+1))],'-r')
	plot([n (2*n)],[(2^(js+1)) (2^(js+1))],'-r')
end
axis ij;  axis tight;

set(gca,'YTick',[8 24 48 96 192]);
set(gca,'YTickLabel','j=3|j=4|j=5|j=6|j=7');
set(gca,'XTick',[]);
set(gca,'XTickLabel','');
text(64,-10,'s=1');
text(192,-10,'s=2');
hold off

% DrawFig6b: Describe layout of ortho-ridgelet coefficients

n=128; 

subplot(1,2,2)

%outside box
plot([0 (2*n) (2*n) 0 0],[0 0 (2*n) (2*n) 0], '-r'); hold on
plot([n n],[0 (2*n)],'-r');

J = round(log2(n));

% vertical lines
for is=3:(J-2),
	plot([(2^is) 2^is],[0 (2^(is+3))],'-r')
	plot([(n + 2^is) (n + 2^is)],[0 (2^(is+3))],'-r')
end

plot([(2^(J-1)) 2^(J-1)],[0 (2^(J+1))],'-r')
plot([(n + 2^(J-1)) (n + 2^(J-1))],[0 (2^(J+1))],'-r')



%horizontal lines
for js=3:(J-1),
	plot([0 n],[(2^(js+1)) (2^(js+1))],'-r')
	plot([n (2*n)],[(2^(js+1)) (2^(js+1))],'-r')
end


axis ij;  axis tight;

set(gca,'YTick',[8 24 48 96 192 ]);
set(gca,'YTickLabel','j=3|j=4|j=5|j=6|j=7');
set(gca,'XTick',[1 12 24 48 96 129 140 152 176 224]);
set(gca,'XTickLabel','i=3|4|5|6|7|3|4|5|6|7');
text(64,-10,'s=1');
text(192,-10,'s=2');
hold off

%print -depsc RawFig6b.eps  

%
% 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]
%