pso图像配准 经典算法

pso图像配准 经典算法 最优值查找 论文参考算法

%Caculate MI of a and b in the region of the size（a）

function [mi] = EMI（ab）
%
[MN] = size（a）;
% change a b to range 1:n
if min（a（:））<1
    a = a-min（a（:））+1;
end
if min（b（:））<1
    b = b-min（b（:））+1;
end
    

% ncolora = max（a（:））-min（a（:））+1;
% ncolorb = max（b（:））-min（b（:））+1;
ncolora = max（a（:））;
ncolorb = max（b（:））;

hab = zeros（ncolorancolorb）;
ha = zeros（1ncolora）;
hb = zeros（1ncolorb）;

%count begins
[MN]=size（a）;
for i=1:M
    for j=1:N
       indexx =  a（ij）;
       indexy =  b（ij）;

       hab（indexxindexy） = hab（indexxindexy）+1;
       ha（indexx） = ha（indexx）+1;
       hb（indexy） = hb（indexy）+1;

   end
end

%cal hahbhab
hsum = sum（sum（hab））;
index = find（hab~=0）;
p = hab/hsum;
Hab = sum（sum（-p（index）.*log（p（index））））;

hsum = sum（sum（ha））;
index = find（ha~=0）;
p = ha/hsum;
Ha = sum（sum（-p（index）.*log（p（index））））;

hsum = sum（sum（hb））;
index = find（hb~=0）;
p = hb/hsum;
Hb = sum（sum（-p（index）.*log（p（index））））;

mi = Ha+Hb-Hab;
% cal EMI
% feature efficiency coefficient
n =2/5;
if Hab ~= 0
    e = （mi^n） / （ Hab^（1-n））;
else% not usual happen just for test
    e = -1;%（mi^n） / （ Hab^（1-n）+1e-10）;
end
mi = e;
% % 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2005-11-27 20:08  pso图像配准\
     文件        1218  2005-10-19 20:09  pso图像配准\EMI.m
     文件         697  2005-10-19 20:09  pso图像配准\FPMI.m
     文件        3338  2005-10-19 20:10  pso图像配准\GMI.m
     文件         228  2005-10-19 20:10  pso图像配准\gobad.m
     文件        1404  2005-10-19 20:10  pso图像配准\HiMI.m
     文件        1034  2005-10-19 20:10  pso图像配准\MI.m
     文件         653  2005-10-19 20:10  pso图像配准\myMI.m
     文件         439  2005-10-19 20:10  pso图像配准\oneDimSearch.m
     文件        2528  2005-10-19 20:10  pso图像配准\POWELL.m
     文件       35888  2005-10-19 20:10  pso图像配准\proj.fig
     文件       12544  2005-10-19 20:10  pso图像配准\proj.m
     文件        5902  2005-10-19 20:10  pso图像配准\PSO.m
     文件         390  2005-10-19 20:10  pso图像配准\restore.m
     文件        3719  2005-11-27 20:04  pso图像配准\RMI.m
     文件          67  2005-10-19 20:10  pso图像配准\test.m
     文件         238  2005-10-19 20:10  pso图像配准\test2.m
     文件        7648  2005-10-19 20:10  pso图像配准\testListBoxlj.fig
     文件        7936  2005-11-27 20:02  pso图像配准\testListBoxlj.m
     文件         326  2005-10-19 20:10  pso图像配准\说明.txt