drdrawfig15.m

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

function dRDrawFig15(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
if LoadData
   load HalfDome2 
end
% Draw HalfDome & Pseudopolar FFT of HalfDome

if LoadData
    img = HalfDome2;
    else
    img = MakeHalfDomeImage(128,30,pi/3,3);
end
%

  %  figure;
    subplot(1,2,1)
	imagesc(img); %
    if HaveTitle
       title('HalfDome ');
    end
    axis('image');
    set (gca,'XTick',[]);
    set (gca,'YTick',[]);colormap(1-hot)
	
    %print -depscRawFig15a.eps
    
   
    subplot(1,2,2)
    Y = PseudopolarFFT(img);
    imagesc(sqrt(abs(Y)));colormap(1-hot);axis xy;axis square
    xlabel('theta')
    ylabel('sqrt(lambda)')
    if HaveTitle
       title('PseudoPolar Fourier transform')
    end
    %
    set (gca,'XTick',[]);
    set (gca,'YTick',[]);
    %print -depscRawFig15b.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]
%