drpdrawfig13.m

beamlet变化的工具箱

function dRPDrawFig13(OpenNewWindow,HaveTitle,LoadData)

%DrawFig 13:   % show the fft domain with the different bases + entropy Measure



if nargin<3
    LoadData=0;
end

if nargin<2
    HaveTitle=1;
end

if nargin<1
    OpenNewWindow=1;
end

if OpenNewWindow
    figure
end


LoadData=0;

% build several types of Bases
D = [3,3];
btree1 = BuildRPBasisTree2(D,'Diadic1',[3]);
btree2 = BuildRPBasisTree2(D,'Vertical',3);
btree3 = BuildRPBasisTree2(D,'FIO',[1]);
btree4 = BuildRPBasisTree2(D,'Congruent',D);
noBases = 4;
plotDim = [2 2];

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%% IMAGE 1: concentrated in angular (oriented texture)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%im1 = readImage
im1 = imread('linesTexture.tif');
im=double(im1(1:64,190:190+63));

im = im ./max(im(:));
img0  = im - mean(mean(im));

Entropy = 'N(eps)';
par = 1.5;

% show the fft domain with the different bases + entropy Measure
FFTGlobal  = fftshift(fft2(im));
FFTGlobalShow = FFTGlobal([2:size(FFTGlobal,1)/2+1,size(FFTGlobal,1)/2+1:size(FFTGlobal,1)],:);
%figure; 
subplot(plotDim(1),plotDim(2),1);
for (i = 1:noBases)
   subplot(plotDim(1),plotDim(2),i);
   AutoImageInv(log(1+abs(FFTGlobalShow)));
   Polrax = axis; hold on;
   eval(['btree = btree' int2str(i) ';']);
   PlotRectoPolarRPPartition(btree,'y',Polrax,D);
   % calculate statistics - already calculated above
   coef  = FPT2_RPkt(btree,img0,D);
   e = CalcEntropy(coef,Entropy,par);
   
   title(['Entropy = ' num2str(e)]);
   set(gca,'XTick',[]);
   set(gca,'YTick',[]);
end;
colormap(1-hot)
%
% 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]
%