dmakessfigure4.m

beamlet变化的工具箱

function MakeSSFigure4(varargin)
% MakeSSFigure4.m -- Figure 4 of Slant Stack Paper
% Usage: MAKESSFIGURE4 will start a new figure window and plot Figure 4 of Slant Stack Paper.
%        MAKESSFIGURE4(n) for any argument type n will plot the figure on the current figure 
%        window and plot Figure 4 of Slant Stack Paper.
if nargin == 0
	 figure;
else
     gcf;
end

% Illustrate Slant Stack --
n1 = 64;
n2 = 128;
ell = 29;
theta = atan(2*ell/n1);
%
Img    = ReadImage('Lincoln');
BigImg = zeros(n2,n1);
BigImg((n1+1 - (n1/2)):(n1+ n1/2),:) = Img;
%
ExactShift  = zeros(n2,n1);

%
% Build matrices of Difference
for i=1:n1,
	vshift = Shift1dSignal(BigImg(:,i),(i-n1/2-1).*tan(theta));
	ExactShift(:,i) = real(vshift(:));
end
%

%figure;
subplot(1,2,1);
imagesc(BigImg);
axis('image')
title('Lines at Slope 58/64, Ordinary Image')
colormap(gray)
subplot(1,2,2);
imagesc(ExactShift);
axis('image')
title(sprintf('Horizontal Lines; Image Sheared atan(58/64)',theta))
colormap(gray)

% superpose families of lines.

subplot(1,2,1)
slope = 2*ell/n1;
hold on;
for k=-n1:(n2-1),
	if rem(k,2)== 0,
		plot([1 n1],[(n1+1+k) (n1+1+k - n1*slope)],'y')
	end
end

subplot(1,2,2)
hold on;
for k=-n1:(n1-1),
	if rem(k,2)==0,
		plot([1 n1],[(n1+1+k) (n1+1+k)],'y')
	end
end
%
% 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]
%