drdrawfig14.m

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

function dRDrawFig14(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

% Normalized Slant Stack of Half Dome
   
%
if LoadData
    load dRDrawFig14
    else
   Img=MakeHalfDomeImage(128,pi/8, 1.3,3);
 	S = FastSlantStack_pre(Img);
end

    


  
  %  figure 
    subplot(1,2,1)
	imagesc(abs(real(S)));
    %
    if HaveTitle
       title('Slant Stack of Half Dome')
    end
	axis('image');%axis off
    set (gca,'XTick',[]);
    set (gca,'YTick',[]);
    colormap(1-hot)
    %colormap(1-gray)
    %print -depsc RawFig16b.eps
    
    s8filter = MakeONFilter('Symmlet',8);
    wt = FWT2_PO(S,3,s8filter);
	zmat = sqrt(abs(wt));
	zmat = 256-3.8*zmat;
    
    n=2*128;
    %figure
    subplot(1,2,2)
    %plot([0 (n) (n) 0 0],[0 0 (n) (n) 0], '-k');
    hold on
	%imagesc(zmat);
    RGBimage(zmat,colormap(1-hot));axis('image');
    if HaveTitle
       title('WT2[Slant Stack of Half Dome]');
    end
	

J = round(log2(n));

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

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

    set (gca,'XTick',[]);
    set (gca,'YTick',[]);%axis off;
    hold off
    %print -depsc RawFig17.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]
%