testadjfastslantstack.m

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

% Test Adj_FastSlantStack
%Makes two figures, showing the adjoint of the Slant stack of points in image space, and Lenna 
	figure;
	X = zeros(32);
	X(6,9) = 1;
	S = FastSlantStack(X);

	subplot(1,3,1)
	imagesc(X)
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);title('X')
	colormap(gray)
	subplot(1,3,2)
	imagesc(real(S));set (gca,'XTick',[]);
    set (gca,'YTick',[]);title('S(X)')
	axis('image')
	
	XX = Adj_FastSlantStack(S);
	subplot(1,3,3)
	imagesc(real(XX))
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);title('Adj(S(X))')
	%
	norm(imag(S))
	%
%
	Lenna = ReadImage('Lenna');
	SLenna = FastSlantStack(Lenna);
	figure;
	subplot(1,3,1)
	imagesc(Lenna)
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);title('Lenna')
	colormap(gray)
	subplot(1,3,2)
	imagesc(real(SLenna))
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);title('S(Lenna)')
	
	XLenna = Adj_FastSlantStack(SLenna);
	subplot(1,3,3)
	imagesc(real(XLenna))
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);title('Adj(S(Lenna))')
	%
	norm(imag(SLenna))
	norm(imag(XLenna))
	%
	
%
% 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]
%