梯度图像的互信息测度-Read.PPT

基于互信息的医学图像配准 许向阳 华中科技大学 医学图像信息研究中心 讨论内容 互信息的概念 互信息的改进 搜索策略的优化 1. 互信息的概念 互信息是信息理论中的一个基本概念。 1948年C. E. Shannon提出， A mathematical theory of communication,Bell System Technical Journal, vol. 27, pp. 379–423/623–656,1948. 1995年Colligon和Viola首次用于医学图像配准。 互信息用于描述两个系统间的统计相关性，或者一个系统中包含的另一个系统的信息的多少。 对于来自同一对象的两幅图像，配准后，互信息应为最大。 互信息只是一种相似性测度 基于互信息的配准方法仍然遵循一般的配准框架 1. 互信息的概念 Collignon A, et al, Automated multi-modality registration based on information theory, Information Processing in Medical Imaging, 1995,263-274 Viola P, Wells WM, Alignment by Maximization of mutual information, International Conference on computer Vision, 1995, Los Alamitos,CA: IEEE Computer Science Press,1995,16~23. 1. 互信息的概念 熵 熵（entropy）是信息的测度。熵的概念来源于通讯理论。 Pi = 灰度值为i的像素点数 / 总像素点数 1. 互信息的概念 熵可以有三种解释： 当一个事件发生时，其给出的信息数量； 一个事件的结果的不确定性； 事件发生的概率分布。 Question: 熵大，说明什么问题？ 熵小，又说明什么？ 1. 互信息的概念 熵 Jumarie entropy 不适于多模态图像配准 Renyi entropy 在图像配准中很少使用 1. 互信息的概念 联合熵 联合熵H(A,B)是随机变量A和B相关性的统计量。PAB(a,b)是联合概率分布函数。 1. 互信息的概念 1. 互信息的概念 条件熵 1993年，Woods应用条件熵作为配准的测度，开发了PET-MT图像配准的软件AIR. 1. 互信息的概念 互信息 1. 互信息的概念 互信息 1. 互信息的概念 式(2)含有项H(A,B)，表示最大互信息与最小联合熵有关。 联合熵表明了两幅图的灰度值分布。 在通常情况下，联合熵小，表明对齐的程度高。 但是，采用两幅图像重叠区域的联合熵可能会出现两幅图像完全误配时，联合熵最小的情况。 如一幅图像例的一部分背景区域和另一幅图像的一部分背景区域重叠，其它部分不重叠时，联合熵最小。 采用互信息可以避免这种情况，因为对只有部分背景对齐的情况，H(A)和H(B)很小，因而互信息很小。随着A、B重叠区域的增大，A、B中包含更多的信息，使得H(A)、H(B)增大。尽管联合熵也增大，但是互信息在增加，更好的反映了对齐的程度。 1. 互信息的概念 互信息 1. 互信息的概念 圆代表了图像的边缘熵 两圆合并区域为联合熵 重叠部分为互信息 互信息的优点 分割和特征提取方法中存在的问题 在图像的获取和特征提取过程中会产生不同程度的噪声，使得提取到的特征点位置存在一定的误差，要使一幅图像中的特征点精确匹配另一幅图像中的特征点是很困难的； 从两幅图像中提取到的图像特征点集数目是不等的，确定它们之间的对应关系较难； 一幅图像中的某些特征点在另一幅图像中没有相对应的特征点，即存在着出界点； 特征点集之间的变换可能是刚性的，也可能是非刚性的。 互信息的性质 非负性：I(A,B)=0。B的不确定性不会因为对A的了解而增大。 独立性：当PAB(a,b) = PA(a)PB(b)时，I(A,B)=0。若A和B是独立的，则不能从一幅图像中获得任何关于另一幅图像的信息。 对称性：I(A,B) = I(B,A) 自信息：I(A,A) = H(A) 有界性：I(A,B) =min(H(A),H(B)) 互信息的性质 I(A,B) = min (H(A), H(B)) = (H(A) + H(B)) / 2 = max (H(A), H(B)) = H(A,B) = H(A) + H(B) 2. 互信息测度的改进