drdrawfig18.m

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

function dRDrawFig18(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

%Study of the kernel of the ORT, n=128, Subband display
% 

n=128;
par = [ 4  8 2  2 1 ];
theta = zeros(2*n,2*n);
theta_display = zeros(2*n,2*n);
theta(2^par(1)+par(2)+1,(2^par(3)+par(4)+1)+(par(5)-1)*n) = 1;
theta_display(2^par(1)+par(2)+1:2^par(1)+par(2)+3,(2^par(3)+par(4)+1)+(par(5)-1)*n-3:(2^par(3)+par(4)+1)+(par(5)-1)*n) = 1;

if LoadData
    load 'dRDrawFig16.mat' ridgelet theta_new
else
    ridgelet = Inv_FastOrthoRidgeletTrans(theta);
    theta_new = FastOrthoRidgeletTransform(ridgelet);
end

    
J=log2(n);
zmat = sqrt(abs(theta_new));
zmat = 256-3.8*zmat;


if par(5)==1
   SubBand=[zmat(2^par(1)+1:2^(par(1)+1),1:2^3) ...
   zmat(2^par(1)+1:2^(par(1)+1),2^3+1:2^4) ...
   zmat(2^par(1)+1:2^(par(1)+1),2^4+1:2^5)];
   else
   SubBand =[zmat(2^par(1)+1:2^(par(1)+1),n+1:n+2^3) ...
   zmat(2^par(1)+1:2^(par(1)+1),n+2^3+1:n+2^4) ...
   zmat(2^par(1)+1:2^(par(1)+1),n+2^4+1:n+2^5)];
end
%figure
imagesc(abs(SubBand))
colormap(hot)
axis('square');
    set (gca,'XTick',[]);
    set (gca,'YTick',[]);
if HaveTitle
    title('Single SubBand')
end
%print -depsc RawFig18.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]
%