www.pudn.com > 1111111registrationbasedoncorrelation.rar > HiMI.m

% -- use a 9th vector as a pixel 
% -- cal MI as stastic on pixels 
function HiMI = HiMI(Is0, It0) 
% - read images 
% load ../MItest/t2.mat; 
%Is0 , It0 
[rt ct]=size(It0); 
[rs cs] = size(Is0); 
if rt>rs ||ct>cs 
    x1 = (ct-cs)/2+1; 
    x2 = cs+x1-1;  
    It0=It0([x1:x2],[x1:x2]); 
end 
% - replace images to 9th images, a pixel-> 9th vector 
Is9 = i_9th(Is0); 
It9 = i_9th(It0); 
 
% - cal MI over the new images 
 
% MI0 = calMI_tian_1(double(Is0(2:(r-1),2:(c-1))), double(It0(2:(r-1),2:(c-1)))); 
%MI9 = calMI_tian_1(double(Is9), double(It9)); 
MI9 = MI(double(Is9), double(It9)); 
 
HiMI = MI9; 
% - disp MI 
% figure, imshow(RIs); 
% RIr = get_Rect(Ir, my_Rect); 
% figure, imshow(RIr); 
% figure,plot(MI); 
 
function I9 = i_9th(I0) 
[row col]=size(I0); 
%I9=zeros((row-1)*9,(col-1)); 
%for i = 2:(row-1) 
%    for j=2:(col-1) 
%%         p = I0(i,j); 
 
%        i9 = (i-1)*9; 
%       j9 = j-1; 
%        V = I0((i-1):(i+1),(j-1):(j+1)); 
%        I9((i9-8):i9,j9) = V(:); 
         
      i = [2:(row-1)]; j = [2:(col-1)];     
    Pa1 = I0(i-1, j-1)'; 
    Pa2 = I0(i,   j-1)'; 
    Pa3 = I0(i+1, j-1)'; 
    Pa4 = I0(i-1, j  )'; 
    Pa5 = I0(i,   j  )'; 
    Pa6 = I0(i+1, j  )'; 
    Pa7 = I0(i-1, j+1)'; 
    Pa8 = I0(i,   j+1)'; 
    Pa9 = I0(i+1, j+1)'; 
        I9 = [ Pa1(:)'; Pa2(:)';Pa3(:)'; Pa4(:)';... 
          Pa5(:)'; Pa6(:)';Pa7(:)'; Pa8(:)'; Pa9(:)']; 
    end 
end