互信息的数学解释以及matlab编程_matlab求互信息函数

作者：2023面试高手 | 2024-03-15 05:01:30

踩

matlab求互信息函数

1.互信息

互信息是信息论中的一个基本概念，通常用于描述两个系统间的统计相关性，或者是一个系统中所包含另一个系统中信息的多少。

引入互信息对图像进行配准是因为：不需要假定两幅图像灰度的对应关系，不需要对图像进行分割。但是这对Ct-MR配准、或者CT-PET配准是可行的，对于CT-US配准是不可行的。

在概率论中，两个随机变量A和B，他们的边缘概率分布为p_A (a)和p_B (b)，他们的联合概率分布为p_AB (a,b)。当p_AB (a,b)=p_A (a)∙p_B (b)时A与B是相互独立的。

在多模医学图象配准问题中 ,虽然两幅图象来源于不同的成像设备 ,但是它们基于共同的人体解剖信息 ,所以当两幅图象的空间位置完全一致时 ,其中一幅图象表达的关于另一幅图象的信息 ,也就是对应象素灰度的互信息应为最大。互信息I(A,B)通过计算p_AB (a,b)和p_A (a)∙p_B (b)的差距来 得到A和B的依赖程度 。

2.用熵来表示互信息

熵表达的是一个系统的复杂性或者是不确定性。

系统熵的定义：

两个系统的联合熵：

条件熵：

那么，互信息可以表示为：

对于离散的数字图像，联合概率分布P_AB(a,b)可以利用归一化的联合直方图来表示：

那么，边缘概率P_A(i)可以表示为：（即在联合直方图，沿着行或者列求和）

同理，边缘概率P_B(i)可以表示为：（即在联合直方图，沿着行或者列求和）

3.互信息分类

1.互信息： MuInf = H_A+H_B-H_AB;

2.归一化互信息：MuInf = (H_A+H_B)/H_AB;

4.互信息的Matlab编程


<span style="font-size:18px;">function [MuInf] = MutualInfo(R,F,method)
% 互信息函数
%R/F：可以是二值图像，也可以是灰度图像
%method：互信息与归一化互信息选项，多采用归一化互信息
if size(F,3) == 3
    F = rgb2gray(F);
end
if size(R,3) == 3
    R = rgb2gray(R);
end;
% 计算互信息
Hist = JoinHist(R,F);
[rows,columns] = size(R);
N_Hist = Hist./(rows*columns);%联合直方图归一化
Marg_A = sum(N_Hist); %对所有列单独求和
Marg_B = sum(N_Hist,2); %对所有行单独求和
H_A = 0;
H_B = 0;
for i=1:1:size(N_Hist,1)   %计算熵H_A
    if Marg_A(i) ~= 0
        H_A = H_A+(-Marg_A(i)*log2(Marg_A(i)));
    end
end
for i=1:1:size(N_Hist,2)   %计算熵H_B
    if Marg_B(i) ~= 0
        H_B = H_B + (-Marg_B(i)*log2(Marg_B(i)));
    end
end
H_AB = sum(sum( -N_Hist.*log2(N_Hist+(N_Hist == 0)) ));
 if strcmp(method,'MI')
     MuInf = H_A+H_B-H_AB;
 end
 if strcmp(method,'NMI')
     MuInf = (H_A+H_B)/H_AB;     
 end;
end
 
%% 统计F、R两张图像的联合直方图
function Hist = JoinHist(R,F)
[rows,columns] = size(R);
Hist = zeros(256,256);
for i = 1:1:rows
    for j = 1:1:columns
        Hist(R(i,j)+1,F(i,j)+1) = Hist(R(i,j)+1,F(i,j)+1)+1;
    end
end
end</span>

5.实验仿真

声明：本文内容由网友自发贡献，不代表【wpsshop博客】立场，版权归原作者所有，本站不承担相应法律责任。如您发现有侵权的内容，请联系我们。转载请注明出处：https://www.wpsshop.cn/w/2023面试高手/article/detail/239052