drpdrawfig03.m

beamlet变化的工具箱

function dRPDrawFig03(OpenNewWindow,HaveTitle,LoadData)

%% Plot digital Ridgelet,FIO tiling and wavelet like in Fourier domain

% 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;
n=64;
wplist = [2 3 2; 2 2 1 ; 3 2 1];% 4 3 2 ];
sz = size(wplist); nr = sz(1); 
d = wplist(1,1);
b = wplist(1,2);
k = wplist(1,3);	
wp = MakeWaveletPacket(d,b,k,'Coiflet',3,n);
wp = circshift(wp,[1 20]);

t = (0:(n-1)) ./n;
%wpstr = sprintf('WaveletPacket(%i,%i,%i)',d,b,k);
wFather = zeros(1,64);
wFather([1 9]) =[1 0];
QMFSymm8 = MakeONFilter('Symmlet',8);
x = IWT_PO(wFather,4,QMFSymm8);
B = repmat(wp',1,64);
A = repmat(x,64,1);
C=A.*B;
    
S=zeros(128,128);
S(32:95,32:95) =C; 
    
X =  Inv_FastSlantStack(S);


subplot(2,3,1)
plot(wp,t,'LineWidth',2);axis square;set(gca,'YDir','reverse');axis off
if HaveTitle
    title('Mother Wavelet')
end
subplot(2,3,2);
imagesc(log(1+abs(C)));axis image;
set (gca,'XTick',[]);
set (gca,'YTick',[]);
if HaveTitle
    title('Radon domain')
end
    
subplot(2,3,3)
imagesc(log(1+abs(X)));axis image;
set (gca,'XTick',[]);
set (gca,'YTick',[]);
if HaveTitle
    title('Spatial Domain')
end
subplot(2,3,5)
plot(t,x,'r','LineWidth',2); set(gca,'YDir','reverse');axis off
if HaveTitle
    title('Father Wavelet')
end



%
% 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]
%