一种基于修正的多尺度retinex算法对医疗图像进行增强的方法

摘要

一种基于修正的多尺度retinex算法对医疗图像进行增强的方法，它包括以下步骤：S1、采用多种医疗图像采集设备获取同一病灶或感兴趣部位的图像序列；S2、采用图像增强方法分别对前述各种医疗图像采集设备获取的图像序列进行处理；S3、对增强后的所有图像序列进行像素级融合获取融合后的图像序列；S4、采用图像增强方法对融合后的图像序列进行处理，得到成品图像序列。本发明通过改进的retinex算法，增加了边缘强度，有利于人体内相近器官的区分，将肿瘤和脏器进一步分别开来，增加了病灶或感兴趣部位的可视性，有助于临床诊断、放射治疗计划的制定和评价。

说明书

技术领域

本发明涉及医疗图像处理领域，尤其是加强多源医学图像的成像结果并适用于在像素级精准融合的图像显示方法，具体地说是一种能够增加病灶或感兴趣部位的可视性，有助于临床诊断、放射治疗计划的制定和评价的基于修正的多尺度retinex算法对MR和CT图像进行增强的方法。

背景技术

随着医学、计算机技术及生物工程技术的发展，医学影像为临床诊断提供了多种模态的医学图像，如CT(计算机X线断层扫描)、MRI(磁共振成像)、SPECT(单光子发射计算机断层成像)、PET(正电子发射计算机断层扫描)、DSA(数字减影血管造影技术)、超声图像、电阻抗图像等。

不同的医学图像提供了相关脏器的不同信息，比如CT图像具有很高的分辨力骨骼成像的能力，非常清晰，对病灶的定位提供了良好的参照，但对病灶本身的显示就较差。MRI虽然它空间分辨力比不上CT图像但是它对软组织成像清晰，有利于病灶范围的确定，可是它又缺乏刚性的骨组织作为定位参照。而PET尽管提供了脏器的新陈代谢功能信息，但对解剖结构的描画却很差。而MRI、CT、X射线成像对人体解剖结构描画得很好，却缺乏人体的功能信息。可见不同模态的医学图像都有各自的优缺点。如果我们能把它们之间的互补信息综合在一起，充分显示形态成像方法的分辨力高、定位准确这一优势，克服功能成像中空间分辨力和组织对比分辨力低的缺点，最大限度地挖掘影像信息把它们作为一个整体来表达，那么在实际临床应用中，就能为医生提供所需要的足够信息，以便了解病变组织或器官的综合信息，从 而作出准确的诊断或制订出合适的治疗方案。

经典的色彩视觉理论认为:人眼的色彩感觉主要取决于光的波长,不同频率的光会给人不同的色彩感觉,物体的颜色是由物体反射光的频率和强度决定的.然而,美国物理学家Edwin Land在20世纪50年代发现有些现象是传统的色彩理论无法解释的,经过近20年的科学实验和分析,Land认为在视觉信息的传导过程中人类的视觉系统对信息进行了某种处理,去除了光源强度和照射不均匀等一系列不确定的因素,而只保留了反映物体本质特征的信息,如反射系数等.当这些描述物体本质特征的相关信息传递到大脑皮层后,经过更为复杂的信息处理,才最终形成人的视觉.基于这样的认识,1977年Edwin Land首次提出了一种被称为Retinex的色彩理论.Retinex这个词本身就是由视网膜Retina和大脑皮层Cortex这两个词组合构成的,于是Retinex理论又被称为视网膜大脑皮层理论。

Land首先提出了Retinex作为人眼感知亮度和色度的视觉模型,其定义理想的图像f(x,y)为：

f(x，y)＝r(x，y)×i(x，y)

即一幅图像f(x,y)可以用环境亮度函数i(x,y)和物体反射函数r(x,y)的乘积来表示。环境亮度函数描述周围环境的亮度，与物体无关；而物体反射函数是指物体反射能力，与照明无关，它包含了景物的细节信息.基于这种模型,得到的环境亮度函数是一种变化缓慢的图像低频信息；而反射函数则包含着图像中的大部分高频细节信息。

但是，上述单尺度Retinex处理方法无法对图像所有灰度范围取得一致的增强效果，在医学图像上，简单应用该方法，无法获取所需的组织信息和病灶特征。

发明内容

本发明的目的是针对人体解剖学成像和功能性成像采用不同工作原理， 对人体骨骼和软组织成像各有优缺点的问题，提出一种基于修正的多尺度retinex算法对医疗图像进行增强的方法。该方法可以加强多源医学图像的成像结果，并适用于在像素级精准融合后再次增强，增加了病灶或感兴趣部位的可视性，有助于临床诊断、放射治疗计划的制定和评价。

本发明的技术方案是：

一种基于修正的多尺度retinex算法对医疗图像进行增强的方法，它包括以下步骤：

S1、采用多种医疗图像采集设备获取同一病灶或感兴趣部位的图像序列；

S2、采用图像增强方法分别对前述各种医疗图像采集设备获取的图像序列进行处理；

S3、对增强后的所有图像序列进行像素级融合获取融合后的图像序列；

S4、采用图像增强方法对融合后的图像序列进行处理，得到成品图像序列。

本发明的步骤S2中，图像增强算法具体为：

S2.1、对步骤S1中获取的各图像序列的灰度图像进行处理，将各像素点灰度值的数据类型转化成float(浮点)型；

S2.2、对各图像序列的灰度图像进行增强处理，具体步骤如下：

S2.2-1、S2.2-1、获取前述灰度图像来源图像的位数，如果是8位，则灰度范围是0～255，如果是16位，则灰度范围是0～4095，在各自对应的整个灰度范围内选取n个尺度因子σ，分别建立对应的高斯环境函数： k表示尺度因子σ即环境函数的编号，k＝1、2、...n；S2.2-2、采用下述公式分别各图像序列中的图像进行处理，得到各自增强后的输出图像R_Mi(x,y)；

其中：i表示任一医疗图像采集设备获取的图像序列中的图像编号；N表示对应医疗图像采集设备获取的图像序列中的图像总数，k表示环境函数的编 号，n表示环境函数的个数；W_k表示与F_k对应的权重系数，I_i(x,y)是第i幅图像的灰度；δ是修正系数(范围是0-4095)，*代表卷积操作；log为对数；F_k(x,y)为步骤S2.2-1所建立的高斯环境函数；(本发明采用高斯函数作为环境函数。其中，卷积操作通过计算机的离散傅里叶变换和反变换实现。实质是将空间域的数学操作变换到频率域简化，然后反变换回来。)

本发明的步骤S2中，图像增强算法还包括：

S2.3、对各种医疗图像采集设备增强后的输出图像R_Mi(x,y)分别进行图像灰度矫正，具体为：

S2.3-1、计算输出图像R_Mi(x,y)的灰度均值avg和标准方差div；

S2.3-2、分别计算出图像R_Mi(x,y)截断拉伸的最大值和最小值；

S2.3-3、计算出图像R_Mi(x,y)的平均强度做为偏移量offset，采用下述公式对图像R_Mi(x,y)的直方图进行拉伸；

S2.3-4、采用步骤(b)的方法进行重心判定，得到拉伸后图像直方图的重心；

用低于该重心的所有灰度即横坐标点数值乘以相应纵坐标值即点数后累计求和，除以图像总点数，得到下限阈值；

用高于该重心的所有灰度即横坐标点数值乘以相应纵坐标值即点数后累计求和，除以图像总点数，得到上限阈值；

将前述上限阈值和下限阈值带入步骤S2.3-3中，作为截断拉伸的最大值和最小值，再次对图像直方图进行线性拉伸。

本发明的步骤S2.3-2和S2.3-3之间，还包括先用gammar拉伸对输出图 像R_Mi(x,y)做处理，使得暗部的对比度得以增强，亮部的对比度受到抑制。

本发明的S2.3-2中，α为比例系数，取值1.5～3。

本发明的步骤S2中，图像增强算法还包括：

S2.4、对灰度矫正后的输出图像进行边缘增强，具体为：

S2.4-1、采用5×5的高斯模板对输出图像以像素为单位进行滑动邻域操作，得到经过高斯模糊的图像矩阵；

S2.4-2、然后用输出图像R_Mi(x,y)的图像矩阵减去经过高斯模糊后的图像矩阵，得到图像的边缘矩阵，给该边缘矩阵乘以调整参数，得到增强的边缘矩阵；

S2.4-3、将前述增强的边缘矩阵与输出图像R_Mi(x,y)的矩阵相加，得到经过边缘增强的图像。

本发明的滑动邻域操作具体为：利用滑动邻域对灰度图像做高斯模糊，邻域块空白的地方补0，对输入图形的每一个像素，指定的滑动邻域操作决定输出图像相应的像素值；即每当指定的操作从图像矩阵的一个位置转移到另一个位置时，滑动邻域也以相同的方向运动。

本发明的调整参数的取值大于1。

本发明的有益效果：

本发明通过加强人体器官不同灰阶的细节，将相近器官更好的区别开来，提升了人体骨骼、脏器以及病变结构的对比度。在灰度拉伸阶段运算速度较快。

由于医学图像标准DICOM的灰阶比较多，本发明通过计算直方图面积和重心，更好选择了多尺度因子，还原了人体器官不同灰阶的细节，提升了骨骼、脏器和病变结构对比度，有利于将人体内的病灶和解剖学结构与背景区别开来。

本发明通过改进的retinex算法，增加了边缘强度，有利于人体内相近器官的区分，将肿瘤和脏器进一步分别开来。灰度拉伸使用gamma拉伸和阈值 选定的线性拉伸，运算速度较快，更加有利于医学的灰度图进一步增加对比度。

具体实施方式

下面结合实施例对本发明作进一步的说明。

一种基于修正的多尺度retinex算法对医疗图像进行增强的方法，它包括以下步骤：

S1、采用多种医疗图像采集设备获取同一病灶或感兴趣部位的图像序列；

S2、采用图像增强方法分别对前述各种医疗图像采集设备获取的图像序列进行处理；

S3、对增强后的所有图像序列进行像素级融合获取融合后的图像序列；

S4、采用图像增强方法对融合后的图像序列进行处理，得到成品图像序列。

本发明的步骤S2中，图像增强算法具体为：

S2.1、对步骤S1中获取的各图像序列的灰度图像进行处理，将各像素点灰度值的数据类型转化成float型；

S2.2、对各图像序列的灰度图像进行增强处理，具体步骤如下：

S2.2-1、获取前述灰度图像来源图像的位数，如果是8位，则灰度范围是0～255，如果是16位，则灰度范围是0～4095，在各自对应的整个灰度范围内选取n个尺度因子σ，分别建立对应的高斯环境函数：k表示尺度因子σ即环境函数的编号，k＝1、2、...n；

尺度因子选取方式如下，以n＝5为例；

例如，如果以“面积占比”和“重心判定”方式选择，如下所述：

(a)、获取灰度图像的直方图，沿直方图中横轴方向从左往右获取起点至横坐标对应位置处的灰度面积，除以直方图的总灰度面积得到面积占比，在面积占比达到设定比例A1、A2时，选取对应比例的横坐标点作为第一、 二个尺度因子；(横轴为0～4095灰度，纵轴为点数。我们设每个灰度的底边长为1，乘以纵轴数值为该灰度面积。除以图像点数为面积占比。然后，在面积占比达到一定数值时，以20％，40％为例选取两个横坐标点作为第一、二个尺度因子。)

(b)在灰度图像的直方图中，将所有横坐标点的灰度值乘以对应位置的纵坐标值，然后进行求和，除以图像总点数，得出图像重心作为第三个尺度因子；(第三个尺度因子选择比较特殊，采用的是“重心判定”的方法，先用每个横坐标灰度乘以该坐标下的纵坐标值，求出0～4095范围内的总和(实际计算中要把横坐标值加一位，也就是1～4096)，然后除以图像总点数，得出的数值即为图像重心。)

(c)、在灰度图像的直方图中，沿直方图中横轴方向从右往左获取起点至横坐标对应位置处的灰度面积，除以直方图的总灰度面积得到面积占比，在面积占比达到设定比例A3、A4时，选取对应比例的横坐标点作为第四、五个尺度因子；

在实际工程环境下，可以采用控件的方法目视不同尺度因子选择产生的增强效果，以最优化输出结果；A1的范围是15％-25％，；A2的范围是35-45％；A3的范围是35％-45％；A4的范围是15-25％。

图像F(x,y)与高斯函数卷积的本质在于特殊灰阶段通透滤波。高斯滤波器的特性。

S2.2-2、采用下述公式分别各图像序列中的图像进行处理，得到各自增强后的输出图像R_Mi(x,y)；(用来控制环境函数范围的尺度；R>

其中：i表示任一医疗图像采集设备获取的图像序列中的图像编号；N表示对应医疗图像采集设备获取的图像序列中的图像总数，k表示环境函数的编号；W_k表示与F_k对应的权重系数，I_i(x,y)是第i幅图像的灰度；δ是修正系数(范围是0-4095)，*代表卷积操作；log为对数；F_k(x,y)为步骤S2.2-1所建立的高斯环境函数；(本发明采用高斯函数作为环境函数。其中，卷积操作通过计算机的离散傅里叶变换和反变换实现。实质是将空间域的数学操作变换到频率域简化，然后反变换回来。)

进一步的，所述的步骤S2中，图像增强算法还包括：

S2.3、对各种医疗图像采集设备增强后的输出图像R_Mi(x,y)分别进行图像灰度矫正，具体为：

S2.3-1、计算输出图像R_Mi(x,y)的灰度均值avg和标准方差div；

S2.3-2、分别计算出图像R_Mi(x,y)截断拉伸的最大值和最小值；

α为比例系数，取值1.5～3；

S2.3-3、计算出图像R_Mi(x,y)的平均强度做为偏移量offset，采用下述公式对图像R_Mi(x,y)的直方图进行拉伸；

S2.3-4、采用步骤(b)的方法进行重心判定，得到拉伸后图像直方图的重心；

用低于该重心的所有灰度即横坐标点数值乘以相应纵坐标值即点数后累 计求和，除以图像总点数，得到下限阈值；

用高于该重心的所有灰度即横坐标点数值乘以相应纵坐标值即点数后累计求和，除以图像总点数，得到上限阈值；

将前述上限阈值和下限阈值带入步骤S2.3-3中，作为截断拉伸的最大值和最小值，再次对图像直方图进行线性拉伸。

所述的步骤S2.3-2和S2.3-3之间，还包括先用gamma拉伸对输出图像R_Mi(x,y)做处理，使得暗部的对比度得以增强，亮部的对比度受到抑制。

进一步的，所述的步骤S2中，图像增强算法还包括：

S2.4、对灰度矫正后的输出图像进行边缘增强，具体为：

S2.4-1、采用5×5的高斯模板对输出图像以像素为单位进行滑动邻域操作，得到经过高斯模糊的图像矩阵；

S2.4-2、然后用输出图像R_Mi(x,y)的图像矩阵减去经过高斯模糊后的图像矩阵，得到图像的边缘矩阵，给该边缘矩阵乘以调整参数，得到增强的边缘矩阵；

S2.4-3、将前述增强的边缘矩阵与输出图像R_Mi(x,y)的矩阵相加，得到经过边缘增强的图像。

滑动邻域操作具体为：利用滑动邻域对灰度图像做高斯模糊，邻域块空白的地方补0，对输入图形的每一个像素，指定的滑动邻域操作决定输出图像相应的像素值；即每当指定的操作从图像矩阵的一个位置转移到另一个位置时，滑动邻域也以相同的方向运动。

普通图像是自然光反射图像，灰度变化范围从0-255，一般采用RGB三尺度处理。医疗图像主要是照射图像，工业上主要采用了dicom、HL7标准，多为12位图像，其灰度达到了0---4095，我们在多尺度选择上根据运算负荷，自适应的调整尺度因子和各尺度权重。

我们采用不同于分段灰度拉伸的方法，取得了更好的对比度和细节。

在灰度矫正之后，我们引入基于高斯模糊模板的边缘检测，并将边缘增 强效果叠加至原图，进一步优化人体解破学结构。

本发明未涉及部分均与现有技术相同或可采用现有技术加以实现。