drdrawfig28.m

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

function dRDrawFig28(OpenNewWindow,HaveTitle,LoadData)

%Test ORT on BandPassed classical Half Dome with n=256;
% 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

%Test ORT on BandPassed classical Half Dome with n=256;
% Original with n=256, takes a lot of time!!;

% 1. Make an image
n=128;
if LoadData
    load 'dRDrawFig24.mat' ansel 
    image = ansel;
else
image = MakeHalfDomeImage(n,pi/8,1.3,3);%n,33, .77,3);

end
%image = MakeHalfDomeImage(n,pi/8,1.3,3);


   
%1.2Bandpass it
    
domain = (-(n/2)):(n/2-1);
[x,y] = meshgrid(domain,domain);
%
r = sqrt(x.^2 + y.^2);
j=5;
sr = r ./ (2^j);
%

Window=BandWindow(sr,3,'bandpass');
Fimage = fftshift(fft2(image));
NewFimage =  Fimage .* Window;
NewImg= abs(ifft2(fftshift(NewFimage)));


subplot(2,2,1)
imagesc(image);colormap(1-hot);axis image
set (gca,'XTick',[]);
set (gca,'YTick',[]);
if HaveTitle
    title('HalfDome')
end
   
subplot(2,2,2);
imagesc(NewImg);axis image;
set (gca,'XTick',[]);
set (gca,'YTick',[]);
if HaveTitle
    title('BandPassed HalfDome')
end


subplot(2,2,3)
imagesc(sqrt(abs(Fimage)));axis image;
set (gca,'XTick',[]);
set (gca,'YTick',[]);
if HaveTitle
    title('FT of HalfDome on sqrt scale')
end

subplot(2,2,4)
imagesc(sqrt(abs(NewFimage)));axis image;
set (gca,'XTick',[]);
set (gca,'YTick',[]);
if HaveTitle
     title('Windowed FT of HalfDome')
end

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