drdrawfig20.m

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

function dRDrawFig20(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, Reconstruction from 10 coefficients Corner Ridgelet
% 

n=128;
par = [ 6  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 dRDrawFig20 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');
    set (gca,'XTick',[]);
    set (gca,'YTick',[]);
colormap(hot) 

if HaveTitle
    (sprintf('Dual Ridgelet s=%i;(j,k)=(%i,%i);(i,l)=(%i,%i)',...
par(5),par(1),par(2),par(3),par(4)));
end
      
      
%Reconstruction form 10 coefficients
    
theta_new10 = theta_new .* (abs(theta_new)./(2*n) >= rprearr(10));

if LoadData
    load dRDrawFig20 ridgelet10 
    else
    ridgelet10 = Inv_FastOrthoRidgeletTrans(theta_new10); 

end
subplot(1,2,2)
imagesc(abs(ridgelet10)); axis('square'); colormap(hsv);
    set (gca,'XTick',[]);
    set (gca,'YTick',[]);
title('Reconstructed from 10 OR Coefficients'); 
%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]
%