drdrawfig21.m

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

function dRDrawFig21(OpenNewWindow,HaveTitle,LoadData)

% 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

%Study of the kernel of the ORT, n=128, Synthesis ans Analysis
%             
n=128;
par = [6  8 2  2 1  ];%
theta = zeros(2*n,2*n);
theta_display = zeros(2*n,2*n);
theta(2^par(1)+par(2)+1,(2^par(3)+par(4)+1)+(par(5)-1)*n) = 1;
theta_display(2^par(1)+par(2)+1:2^par(1)+par(2)+3,(2^par(3)+par(4)+1)+(par(5)-1)*n-3:(2^par(3)+par(4)+1)+(par(5)-1)*n) = 1;

if LoadData
    load dRDrawFig20 ridgelet theta_new
    else
    ridgelet = Inv_FastOrthoRidgeletTrans(theta);
    theta_new = FastOrthoRidgeletTransform(ridgelet);
end

J=log2(n);
subplot(1,2,1)
RGBimage(abs(theta_display),colormap(1-gray(256)));
%axis off;
if HaveTitle
title('Synthesis of Dual Ridgelet ');
end
    set (gca,'XTick',[]);
    set (gca,'YTick',[]);
hold on;

%Box
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 image;axis off;hold off

%print -depsc RawFig17a.eps
zmat = sqrt(abs(theta_new));
zmat = 256-3.8*zmat;
    
subplot(1,2,2)
RGBimage(zmat, colormap(hot));
    set (gca,'XTick',[]);
    set (gca,'YTick',[]);
    if HaveTitle
title('Analysis of Dual Ridgelet ');
end
hold on;
%Box
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 image;axis off;hold off

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