drdrawfig23.m

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

function dRDrawFig23(OpenNewWindow,HaveTitle,LoadData)

% if you don't want to process it, set LoadData=1

if nargin<3
    LoadData=0;
end

if nargin<2
    HaveTitle=1;
end

if nargin<1
    OpenNewWindow=1;
end
if OpenNewWindow
    figure
end
LoadData=0;
% DrawFig23: Example of Z_p^2 Ridgelet
%
wcoef = zeros(1,128);
wcoef(45) = 1;
S8Filter = MakeONFilter('Symmlet',8);
%
wavelet = IWT_PO(wcoef,3,S8Filter);
wRad = zeros(127,128);
wRad(:,56) = wavelet(1:127).';
www = InvRadon_Zp2_Fast(wRad);
wRad(:,66)=wavelet(1:127).';
www2 = InvRadon_Zp2_Fast(wRad);


subplot(1,2,1)
imagesc(abs(www)); axis('equal');axis off;
colormap(hsv); 
if HaveTitle
    title('Example of Z_p^2 Ridgelet; direction=56')
end
%print -deps RawFig23a.eps

subplot(1,2,2)
imagesc(abs(www2)); axis('equal');axis off;
colormap(hsv);
if HaveTitle
    title('Example of Z_p^2 Ridgelet; direction=66')
end

%print -deps RawFig23b.eps
%
% Copyright (c) 2002 Ana Georgina Flesia
%
% 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]
%