testfastslantstack.m

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

% TestFastSlantStack
%Makes a figure showing Slant stack of  points in image space
	figure;
	X = zeros(32);
	X(6,9) = 1;
	S = FastSlantStack(X);

	subplot(3,2,1)
	imagesc(X)
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);title('X')
	colormap(gray)
	subplot(3,2,2)
	imagesc(real(S))
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);title('S(X)')
	%
	norm(imag(S))
	%
	X = 0 .* X;
	X(19,23) =1;
	S = FastSlantStack(X);
	subplot(3,2,3)
	imagesc(X);title('X')
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);
	colormap(gray)
	subplot(3,2,4)
	imagesc(real(S));title('S(X)')
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);
	%
	norm(imag(S))
	%
	X = 0 .* X;
	X(12,28) =1;
	S = FastSlantStack(X);
	subplot(3,2,5)
	imagesc(X);title('X')
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);
	colormap(gray)
	subplot(3,2,6)
	imagesc(real(S));title('S(X)')
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);
	%
%
%	Copyright (c) 2001 David L. Donoho
%
	
	
%
% Part of BeamLab Version:1.1
% Built:Thursday,14-Feb-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: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]
%