drdrawfig16.m

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

function dRDrawFig16(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 ORT, n=128, Reconstruction form 10 coefficients
% 



    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 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;
    
rprearr = reverse(sort(abs(theta_new(:))))./(2*n);  

%figure
subplot(1,2,1)
imagesc(abs(ridgelet));
axis('square');axis off
colormap(1-gray(256)) 
if HaveTitle
       title(sprintf('Dual Ridgelet s=%i;(j,k)=(%i,%i);(i,l)=(%i,%i)',...
par(5),par(1),par(2),par(3),par(4)));

end
set (gca,'XTick',[]);
set (gca,'YTick',[]);
      
%Reconstruction form 10 coefficients
    
theta_new10 = theta_new .* (abs(theta_new)./(2*n) >= rprearr(10));
if LoadData
    load dRDrawFig16 ridgelet10 
    else
    ridgelet10 = Inv_FastOrthoRidgeletTrans(theta_new10); 

end

subplot(1,2,2)
imagesc(abs(ridgelet10)); axis('square'); colormap(hsv);axis off;
if HaveTitle
       title('Reconstructed from 10 OR Coefficients'); 
end

%print -depsc RawFig16.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]
%