• 大小: 19KB
    文件类型: .rar
    金币: 1
    下载: 0 次
    发布日期: 2021-06-13
  • 语言: Matlab
  • 标签: 图像配准  

资源简介

这是一个关于基于互信息的图像配准程序,是以matlab编写

资源截图

代码片段和文件信息

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

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        653  2005-10-19 20:10  peizhun\myMI.m

     文件        697  2005-10-19 20:09  peizhun\FPMI.m

     文件       3338  2005-10-19 20:10  peizhun\GMI.m

     文件        228  2005-10-19 20:10  peizhun\gobad.m

     文件       1404  2005-10-19 20:10  peizhun\HiMI.m

     文件       1034  2005-10-19 20:10  peizhun\MI.m

     文件       1218  2005-10-19 20:09  peizhun\EMI.m

     文件        439  2005-10-19 20:10  peizhun\oneDimSearch.m

     文件       2528  2005-10-19 20:10  peizhun\POWELL.m

     文件      35888  2005-10-19 20:10  peizhun\proj.fig

     文件      12544  2005-10-19 20:10  peizhun\proj.m

     文件       5902  2005-10-19 20:10  peizhun\PSO.m

     文件        390  2005-10-19 20:10  peizhun\restore.m

     文件       3719  2005-11-27 20:04  peizhun\RMI.m

     文件         67  2005-10-19 20:10  peizhun\test.m

     文件        238  2005-10-19 20:10  peizhun\test2.m

     文件       7648  2005-10-19 20:10  peizhun\testListBoxlj.fig

     文件       7936  2005-11-27 20:02  peizhun\testListBoxlj.m

     文件        326  2005-10-19 20:10  peizhun\说明.txt

     目录          0  2005-11-27 20:08  peizhun

----------- ---------  ---------- -----  ----

                86197                    20


评论

共有 条评论