testadjslowslantstack.m

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

% Test Adj_SlowSlantStack
%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 = SlowSlantStack(X);

	subplot(2,3,1)
	imagesc(X);title('X')
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);
	colormap(gray)
	subplot(2,3,2)
	imagesc(real(S));title('S(X)')
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);
	
	XX = Adj_SlowSlantStack(S);
	subplot(2,3,3)
	imagesc(real(XX))
	axis('image');title('Adj(S(X))')
    set (gca,'XTick',[]);
    set (gca,'YTick',[]);

    X = zeros(32);
	X(14,28) = 1;
	S = SlowSlantStack(X);

	subplot(2,3,4)
	imagesc(X);title('X')
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);
	colormap(gray)
	subplot(2,3,5)
	imagesc(real(S));title('S(X)')
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);
	
	XX = Adj_SlowSlantStack(S);
	subplot(2,3,6)
	imagesc(real(XX))
	axis('image');title('Adj(S(X))')
    set (gca,'XTick',[]);
    set (gca,'YTick',[]);
	%
	norm(imag(S))
	%
%
    figure
	S = zeros(64);
	S(44,28) = 1;
	subplot(2,3,1)
	imagesc(S)
	axis('image');title('S1')
	colormap(gray);set (gca,'XTick',[]);
    set (gca,'YTick',[]);
	subplot(2,3,4)
	XX = Adj_SlowSlantStack(S);
	imagesc(real(XX));title('Adj(S1)')
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);
%	
	S = zeros(64);
	S(44,14) = 1;
	subplot(2,3,2)
	imagesc(S)
	axis('image');title('S2')
	colormap(gray);set (gca,'XTick',[]);
    set (gca,'YTick',[]);
	subplot(2,3,5)
	XX = Adj_SlowSlantStack(S);
	imagesc(real(XX));title('Adj(S2)')
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);
%
	S = zeros(64);
	S(28,55) = 1;
	subplot(2,3,3)
	imagesc(S);title('S3')
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);
	colormap(gray)
	subplot(2,3,6)
	XX = Adj_SlowSlantStack(S);
	imagesc(real(XX));title('Adj(S3)')
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);
	
%
% 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]
%