drdrawfig07.m

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

function dRDrawFig07(OpenNewWindow,HaveTitle,LoadData)
% DrawFig7.  Ortho-Ridgelet Transform on HalfDome
% 
% Figure 7. The pseudo ortho-ridgelet coefficients of HalfDome.
% 

if nargin<3
    LoadData=0;
end

if nargin<2
    HaveTitle=1;
end

if nargin<1
    OpenNewWindow=1;
end

if OpenNewWindow
    figure
end

% 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

J=log2(n);
% 2. Take FRT
theta1 = rptheta;

%figure;
subplot(1,2,1)
imagesc(abs(theta1)); colormap(1-gray);
title('Ortho Ridgelet Transform of HalfDome')
hold on

plot([0 (2*n) (2*n) 0 0],[0 0 (2*n) (2*n) 0], '-k');% hold
plot([n n],[0 (2*n)],'-k');

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

%horizontal lines
for js=3:(J-1),
	plot([0 n],[(2^(js+1)) (2^(js+1))],'-k')
	plot([n (2*n)],[(2^(js+1)) (2^(js+1))],'-k')
end
axis ij;  axis square;axis off;hold off
%print -depsc RawFig7a.eps

%figure; 
subplot(1,2,2)
imagesc(sqrt(abs(theta1))); colormap(1-gray);
title('sqrt|Ortho Ridgelet Transform| of HalfDome')
%previous uses log
%imagesc(log(abs(theta1))); colormap(1-gray);colorbar;
%title('log(Ortho Ridgelet Transform) of HalfDome')
hold on

plot([0 (2*n) (2*n) 0 0],[0 0 (2*n) (2*n) 0], '-k');% hold
plot([n n],[0 (2*n)],'-k');

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


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


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