testinvfastslantstack.m

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

% Test InvFastSlantStack
%Makes two figures, showing the inverse of the Slant stack of points in image space, and Lenna 

	X = zeros(32);
	X(14,5) = 1;
	S = FastSlantStack(X);
	XX = Inv_FastSlantStack(S);
	norm(X-XX,'fro'),
	

    
    figure;
	X = zeros(32);
	X(6,9) = 1;
	S = FastSlantStack(X);
    XX = Inv_FastSlantStack(S);
    
	subplot(3,3,1)
	imagesc(X)
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);title('X')
	colormap(gray)
	subplot(3,3,2)
	imagesc(real(S))
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);title('S(X)')
    subplot(3,3,3)
	imagesc(real(XX))
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);title('Inv(S(X))')  
	%
	
	%
	X = 0 .* X;
	X(19,23) =1;
	S = FastSlantStack(X);
    XX = Inv_FastSlantStack(S);
	subplot(3,3,4)
	imagesc(X);title('X')
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);
	colormap(gray)
	subplot(3,3,5)
	imagesc(real(S));title('S(X)')
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);
    subplot(3,3,6)
	imagesc(real(XX))
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);title('Inv(S(X))')  
	%
	norm(imag(S))
	%
	X = 0 .* X;
	X(12,28) =1;
	S = FastSlantStack(X);
    XX = Inv_FastSlantStack(S);

	subplot(3,3,7)
	imagesc(X);title('X')
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);
	colormap(gray)
	subplot(3,3,8)
	imagesc(real(S));title('S(X)')
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);
    subplot(3,3,9)
	imagesc(real(XX))
	axis('image');set (gca,'XTick',[]);
    set (gca,'YTick',[]);title('Inv(S(X))')  

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