共生互信息量的医学图像配准
1图像配准
医学图像可以提供有病变组织或器官的大小、形状、空间关系等详细信息,比如 CT图像可以显示骨骼结
构和组织密度分布情况;MR图像和超声(US)图像提供的则是软组织的信息;PET,SPECT能反映人体的功能
和代谢信息.在实际临床应用中,单一模态的图像往往不能提供医生所需要的足够的信息,通常需要将不
同模态的图像融合在一起,得到更丰富的信息以便了解病变组织或器官的综合情况,从而做出准确的诊断
或制订出合适的治疗方案,而配准则是进行融合的前提.图像配准[123]是指对一幅图像进行一定的几何
变换而映射到另一幅图像中,使得两幅图像中的相关点达到空间上的一致.多模态医学图像配准是目前生
物医学工程中的一个热点问题,也是一个难点问题,对于临床诊断和治疗具有重要意义.互信息量[426]技
术是最近几年提出的进行多模态图像配准一种有效方法.互信息量源于信息论,用于度量两个随机变量之
间的相似性.它是一种自动的,基于像素灰度的方法,它不需要选择标志点或提取图像特征,不需要假设图
像中的灰度值存在某种线性关系,因而在世界范围内得到广泛应用,特别是医学图像处理领域[6].但是,
以互信息量作为相似性测度的方法也存在一定的缺点,因为互信息量对噪声、采样点个数比较敏感,当图
像空间分辨率比较低,有噪声影响和图像部分缺损时容易出现误配.为此很多学者进行了改
进,Studholme[7]提出了归一化的互信息量方法,可以消除对图像间覆盖程度较敏感的问题,但该方法仍
然忽略了图像的空间与方向信息.荷兰学者 Pluim[8]提出了将互信息量与图像梯度相结合的方法,在一
般情况下都能够达到较好的效果,但由于图像梯度本身就对噪声十分敏感,当图像中存在一定的噪声时,
该方法的配准成功率就变得比较低,而且求梯度的过程增加了计算时间,因而不能满足临床需要.本文考
虑到像素点及其邻域内不同方向上的点的关系,将空间与方向信息引入到配准的过程中,并利用多元统计
学中的方法进行快速计算,从而克服了前面提到的问题,达到了较好的配准效果.
2基于互信息量的图像配准方法
图像 X的熵[6]可以度量图像的不确定性,其定义为 H(X)=-∑ip(i)logp(i)(1)其中 p(i)为灰度值 i出现
的概率.当某一个灰度出现的概率为 0时,有 0log0=0.由熵的定义可知熵的大小并不依赖于灰度本身,而
只依赖于这些灰度出现的概率.(a)是一幅 MR图像,(b),(c)分别是对其旋转 90°和置乱变换后的图像.显
然平移、旋转和置乱变换只是改变了像素点的空间位置,并没有改变像素点灰度值的大小,也就没有改变
灰度值出现的概率.所以以上三幅图像具有相同的熵,H=5•0875.由此可见图像灰度值空间位置的改变对
熵并没有影响.在多模态医学图像配准问题中,虽然两幅图像来源于不同的成像设备,但是它们基于人体
共同的解剖信息,所以当两幅图像的空间位置完全一致时,它们的对应像素的灰度互信息达到最大值,即
一幅图像表达的关于另外一幅图像的信息最多.这个信息就是互信息量.标准互信息量和归一化的互信息
量的定义分别为 MI=H(X)+H(Y)-H(X,Y)(2)和 NMI=H(X)+H(Y)H(X,Y)(3)其中 H(X,Y)为图像 X,Y的联合
熵.基于互信息量的刚性配准可以表示为 T0=argmaxTMI(X,T(Y))(4)即找到一个最优变换 T使互信息量
达到最大.在配准过程中 H(X),H(Y)并没有发生太大的变化,变化较大的只是 H(X,Y),也就是说配准过程
中真正起作用的只是联合熵,并没有体现出图像 X,Y的熵.为了充分利用它们的信息本文提出了共生熵和
共生互信息量的概念.
3基于共生互信息量的图像配准方法
对于一阶邻域系统(又称 4邻域系统),每一个像素点都有 4个相邻元素.如(a)所示,x表示像素点,y表示
不同方向的邻域.而对于二阶邻域系统(又称 8邻域系统),每一个像素点都有 8个相邻元素,如(b)所示.
(c)~(f)为 4种不同的邻域结构,即像素点 x在不同方向的邻域关系.灰度共生矩阵
(Co2occurrenceMatrix)正是考虑了以上 4种邻域结构.如果一个像素的灰度值为 i,其相邻像素灰度值
为 j,则这两个灰度像素对同时出现的概率记为 p(i,j).灰度共生矩阵不仅包含图像的灰度统计信息,而
且反映灰度分布的空间与方向信息,即不同方向的一对相邻像素出现的概率.由灰度共生矩阵可以得到共
生熵的概念:HCo(X)=-∑i∑jp(i,j)logp(i,j)(5)它反映了区域内的随机程度,包含了图像灰度关于方
向、相邻元素及幅度变化的综合信息.当所有 p(i,j)都相等时,熵达到最大值.为所示 MR图像在不同方
向的联合直方图两幅图像 X,Y的联合共生熵可以定义为 HCo(X,Y)=-
∑i∑j∑k∑lp(i,j,k,l)logp(i,j,k,l)(6)其中 p(i,j,k,l)表示图像 X中一点灰度值为 i、邻域灰度值
为 j,同时对应图像 Y中的一点灰度值为 k、邻域灰度值为 l同时出现的概率.p(i,j,k,l)可以通过四维
联合直方图计算得到,但是对于一幅灰度值范围从 0~255,256个灰度级的图像来说,四维联合直方图的大
小为 256×256×256×256,在实际操作中是很难实现的.下面从多元分布的角度考虑如何快速得到图像 X,Y
的联合共生熵.n维正态分布随机变量 x=(x1,x2,…,xn)T的联合概率密度为 p(x)=p(x1,x2,…,xn)
=1(2π)n2C12exp-12(x-u)TC-1(x-u)(7)其中 u=(u1,u2,…,un)T为变量 x的均值,C为协方差矩阵.它的
熵[9]为 H(x)=log[(2πe)n2C12](8)由文献[10]可知,对于图像 X,Y,将(i,j,k,l)看作变量,可以认为它
近似服从正态分布,从而可以求出 X,Y的联合共生熵:HCo(X,Y)=log[(2πe)2C12](9)这样就得到了一种
新的相似性准则―――共生互信息量:Co2MI=HCo(X)+HCo(Y)-HCo(X,Y)(10)显然新的准则中不仅包含了
图像灰度的统计信息,而且包含了图像的空间与方向信息.
4实验分析
为了验证本算法速度快、精度高、鲁棒性强的特点,我们考虑二元邻域关系,采用 Powell优化算法.用
Matlab7•0在 PC机(P43•0GHzCPU,512MB内存,WindowsXP操作系统)上进行了实验,并与传统的互信息
量方法以及改进的与梯度信息相结合的 GMI[8]方法进行了比较.最后以 B样条[11]为变形函数,以共生
互信息量为相似性准则,将该方法应用到医学图像的弹性配准中.实验 1.验证本文方法对图像大小的鲁
棒性.由于互信息量是对图像概率分布的计算,在数据量减少时就容易出现误配.下面是对大小为 64×64
的 MR与 PET图像进行配准的实验.由于图像尺寸较小,图像联合直方图的统计总数会大量减少,这就像统
计样本变小一样.这样联合直方图就会成为一个稀疏矩阵,对配准参数的变化变得极为敏感.同时由联合
直方图近似得到的配准图像灰度间的联合分布和边缘分布的精度就会变差,从而严重影响到目标函数的
平滑程度,导致产生明显的局部极值,使配准的优化搜索失败.如所示.将配准后 MR图像的边缘加到 PET
图像上,以检测在边缘处的配准效果.(c)为利用标准互信息量进行配准的结果,很明显该方法在边缘处有
较大的误差.(d),(e)分别为 GMI和利用共生互信息量进行配准的结果,由于这两种方法都考虑了图像的
空间与方向信息,相当于将图像进行了扩展,从而增加了有效信息.可以看到两种方法配准后的图像在边
缘处已经基本对齐.但是 GMI方法所需时间约为 Co2MI方法的两倍.(a)~(d)是完全对准的磁共振 PD加权
和 T1加权图像,及分别加入方差为 0•01和 0•02的高斯噪声后的图像.以含噪声的 T1加权图像为浮动
图像进行旋转,计算旋转后图像与含噪声的 PD加权图像的互信息量.(e)为标准互信息量的变化曲线,横
坐标为旋转角度(±30°),纵坐标为互信息量大小.由于图像噪声较大,削弱了图像间的相关性,得到的曲线
波动较大.(f)是 GMI的变化曲线,正如前面分析到的,由于图像梯度本身就对噪声十分敏感,当图像中存
在一定的噪声时,梯度信息反而恶化了互信息的光滑性,使目标函数波动变大,出现了更多的局部极值点.
(g)是利用共生互信息量得到的曲线,由于共生互信息量考虑的是相应灰度点及其邻域出现的概率,从而
抑制了噪声的影响,增强了配准的准确性和鲁棒性.可以看出利用本文方法求出的曲线比较平滑,并保持
了良好的凸性.由于互信息量只依赖于图像灰度信息,当待配准图像本身包含的有效信息较少时,配准效
果较差.(a)所示的 CT图像反映的主要是骨骼结构等硬组织的信息,灰度变化范围较小,与 MR图像配准时
主要是利用的边缘处信息.如果图像有较大缺损时,就很容易出现误配.(c)为利用标准互信息量配准后的
结果.可以看出两幅图像的边缘没有重合.求出的平移,旋转量为Δx=9,Δy=7,Δr=11,与正确的
Δx=10,Δy=7,Δr=10 存在一定偏差.GMI和 Co2MI不仅包含图像的灰度统计信息,而且反映灰度分布的
空间信息,得到的参数与正确值相吻合,但 GMI方法所需时间大约为 Co2MI方法的两倍.本实验参考
Kybic提出的基于 B样条的弹性配准算法[11],以 B样条为变形函数,以共生互信息量为相似性准则,寻
找最优变换 T.(a)为一幅心脏灌注的 MR图像,(b)为变形后图像.(c)为未配准前两图直接相减的结果,
(d)为利用 Kybic算法配准后图像相减的结果,明显地可以看出在头颈及心脏周围仍然没有配准,而利用
Co2MI算法基本上可以配准,能够满足临床医生的需要.
5结论
本文详细分析了基于互信息量的图像配准算法中出现的鲁棒性问题,提出了共生互信息量的概念,并引入
多元统计的方法来计算图像的联合共生熵.实验结果表明,本方法充分利用了图像的灰度信息和空间信
息,很好地解决了基于互信息量的图像配准中出现的鲁棒性问题.更为重要的是,作为一种一般性的配准
算法,Co2MI算法如同 MI一样,可以应用到图像配准以外的更广阔的领域,如经济学、运筹学、模式识别
等.
