drdrawfig08.m

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

function dRDrawFig08(OpenNewWindow,HaveTitle,LoadData)
% DrawFig8: Inspect some Ridgelets
if nargin<3
    LoadData=0;
end

if nargin<2
    HaveTitle=1;
end

if nargin<1
    OpenNewWindow=1;
end

if OpenNewWindow
    figure
end

n=128; 
% Figure 8. 
% \begin{quote} Figure 9a: $\l =(1,4,8,2,2)$; a narrow
% ridgelet near the center of the picture\\ (as $k=8=2^j/2$
% and $\ell =2$ is small)
% 
% \medskip\ni Figure 8b: $\l=(1,4,10,2,2)$. The parameters
% are similar to
% Figure 9a, with a change in ridge location parameter
% $k$; the object is similar to the one in Figure 9a,
% quasi-translated downwards by
% $(10-8)/2^j$, which is 1/8 of the image extent.
% 
% \medskip\ni Figure 8c: $\l=(2,4,8,2,0)$. The parameters
% $s$ and $\ell$ have changed; the object is similar to
% the one in Figure 8a, only quasi-rotated, both because $s$ has changed,
% and $\ell$ has changed.
% 
% \medskip\ni Figure 8d: $\l=(2,4,10,2,0)$. The parameters
% are similar to those used to produce Figure 9c, only the
% ridge location parameter has changed; we see a
% quasi-translation of the object in Figure 9c rightwards by\\
% $(10-8)/16=1/8$ of the image extent.
% 
% \medskip\ni Figure 8e: $\l=(1,4,8,4,5)$; the angular scale
% parameter has increased; we see a ridgelet
% supported near the edges of the image\\ as $\ell =4$.
% 
% \medskip\ni Figure 8f: $\l=(1,4,8,4,11)$; the angular position
% parameter has increased, and we see a quasi-rotation
% of the ridgelet, by $(11-5)/24\x 1/2$ of the image extent,
% or about $3/32 \x 360^{\circ}$.
% \end{quote}

parlist = [ 4  8 2  2 1 ; ...
			4 10 2  2 1 ; ...
			4  8 2  0 2 ; ...
            4 10 2  0 2 ; ...
            4  8 4  5 1 ; ...
            4  8 4 11 1 ];
fignum = 5;
string = 'abcdefg';
n=64;
npar = 6;


if LoadData
    im=imread('dRDrawFig08.tif');
    imagesc(im);axis image;axis off;
else
    for pp = 1:npar,
	    par = parlist(pp,:);
    
	    theta = zeros(2*n,2*n);
	    theta(dyad2ix(par(1),par(2)),dyad2ix(par(3),par(4))+(par(5)-1)*n) =1;
    
        ridgelet=Inv_FastOrthoRidgeletTrans(theta);
    
        theta_new = FastOrthoRidgeletTransform(ridgelet);

        rprearr = reverse(sort(abs(theta_new(:))))./(2*n);  

	    subplot(2,3,pp)
	    imagesc(abs(ridgelet));
	    axis('square');
	    colormap(1-gray) 
	    title(sprintf('(%s): s=%i;(j,k)=(%i,%i);(i,l)=(%i,%i)',...
		string(pp),par(5),par(1),par(2),par(3),par(4)));
        set (gca,'XTick',[]);
        set (gca,'YTick',[]);
 
    end
end
%print -deps RawFig9.eps
	

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