drpdrawfig09.m

beamlet变化的工具箱

function dRPDrawFig09(OpenNewWindow,HaveTitle,LoadData)

%DrawFig 9:   basis functions derived from tiles of different heights


if nargin<3
    LoadData=0;
end

if nargin<2
    HaveTitle=1;
end

if nargin<1
    OpenNewWindow=1;
end

if OpenNewWindow
    figure
end


LoadData=0;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% basis functions derived from tiles of different heights
%

% create a btree that has a tile with "width" pi/4 to 3pi/8
% and varying height....

% for the different btrees, plot representative basis functions by
% sampling the dc coefficient of the tile in the DCT domain (coef domain) for
% an image that has energy only in the tile located at rows 34-49 and starting at col 65
% (fill the whole tile with 1 in the RPFT domain)

D = [5,3];
subplot(2,4,1);
for d1 = 5:-1:2
   btree = BuildRPBasisTree2(D,'Congruent',[d1,3]);
   RPFT = zeros(128,128);
   RPFT(65-64/(2^d1)+1:65,65:80)=1;
   im = ifft2_rp(RPFT);
   
   %show the tile
   subplot(2,4,6-d1);
   RPFTGlobal = RPFT;
   RPFTGlobalShow = RPFTGlobal([2:size(RPFTGlobal,1)/2+1,size(RPFTGlobal,1)/2+1:size(RPFTGlobal,1)],:);
   AutoImageInv(abs(RPFTGlobalShow));
   Polrax = axis; hold on;
   PlotRPPartition(btree,'r',Polrax,D); 
   axis square
   set (gca,'XTick',[]);
   set (gca,'YTick',[]);
   if HaveTitle
       title(['D1=' int2str(d1)]);
   end
   
   % show the basis function
   subplot(2,4,5+5-d1);
   coef = FPT2_RPkt(btree,im,D);
   % leave only coef(33,33) == the DC of the tile
   dcc = coef(65,65);
   coef = zeros(size(coef));
   coef(65,65) = dcc;
   im  = IPT2_RPkt(btree,coef,D);
   im_sc = scale_im(im);
   imagesc(abs(1-im_sc))
   if HaveTitle
       title(['D1=' int2str(d1)]);
   end
   axis square
   set(gca,'XTick',[]);
   set(gca,'YTick',[]);
   p = get(gca,'Position');
   set(gca,'Position',p+[0 0.2 0 0])
end;

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