一种基于腹部CT医学图像融合分类的检测装置
技术领域
本发明属于医疗图像处理技术领域,特别涉及一种基于腹部CT医学图像融合分类的检测装置。
背景技术
现有腹部医学图像融合技术主要对两种不同模态的医学图像,根据医学图像模态的不同,医学图像融合系统可以分为三种类型:解剖医学图像与解剖医学图像融合、解剖医学图像与功能医学图像融合以及功能医学图像与功能医学图像融合。MRI-PET和MRI-SPECT医学图像融合系统属于解剖医学图像与功能医学图像融合,该系统的输入图像是灰度和伪彩色。飞利浦公司推出的MRI-PET组合一体化套件将商用的MRI成像扫描仪和具有特殊屏蔽的PET结合起来,得到的图像对癌症转移的诊断及术前分期具有重要的临床价值。
由于像素级多尺度融合方法在不同尺度图像上直接对图像的像素值进行处理,能最大限度地保留输入图像的像素信息和提高融合图像的质量。因此,像素级多尺度方法在医学图像融合领域成为一个研究热点。目前,国内外主要从图像分解与重构方法和图像融合规则方法两个方面出发,提出新的像素级多尺度融合方法。在图像分解与重构方法方面,基于频域的图像融合方法利用傅里叶变换和傅里叶逆变换对图像信号进行不同尺度的分解与重构,但是该类方法具有时间复杂度高和运行时间长的特点,这与高实时性的医疗辅助诊断矛盾,并且对处理的实验平台硬件设施和软件设施要求非常高。后来,研究者们提出利用空域滤波器对图像进行处理来对输入图像进行多尺度分解与重构,该类方法虽然能快速地进行图像分解与重构,但是基于空域的图像融合方法抗噪性差。但是现有方法针对不同模态的医学图像采用同一种特征进行融合。
发明内容
本发明旨在解决以上现有技术的问题。提出了一种降低复杂度、提高融合图像的质量、提高检测精确度的基于腹部CT医学图像融合分类的检测装置。
本发明的技术方案如下:
一种基于腹部CT医学图像融合分类的检测装置,用于处理患者的腹部CT医学图像而得到所述患者的疑似病灶区域,其包括:获取模块,获取所述患者的腹部CT扫描图像,启动CT扫描器对人体的腹部进行电影模式CT扫描,并将所获得的每一张CT图片按人体进入扫描区的先后顺序给CT图像的每一图片赋予床位号和层位号,按扫描成像的时间顺序给CT图像的每一图片赋予相位号;将CT图像的每一图片进行阈值分割并获取人体体表轮廓曲线;以通过人体体表轮廓曲线在垂直方向的最低点的水平线为基准,对每一图片自左向右依次竖向逐点扫描,当图片的一列像素点的像素值全为零时记为零,当图片的一列由上而下有一个像素点的像素值不为零时,便记下第一个像素值不为零的像素点的高度值;然后将所记录的数据按自左向右的顺序排列,再以所得到的排列为行,按床位层位号的顺序由上向下排列,得到每一个床位的轮廓矩阵;将每一轮廓矩阵中所有元素的高度值累加求和,从一个床位中挑出所求和的最大值所对应的一个相位的一组图片,然后以该组图片为起点将相位号大于该组图片相位号的每一组图片按相位号由小到大排列,得到每一个床位内一个完整呼吸周期的图片阵列;预处理模块,对获取每一个床位内一个完整呼吸周期的图片阵列进行预处理得到二值化灰度图像;形态学腐蚀与膨胀模块,对所述二值化灰度图像进行形态学腐蚀运算得到腐蚀后的图像,对所述二值化灰度图像进行膨胀运算得到膨胀后的图像;运算模块,对所述腐蚀后的图像进行开运算得到开运算图,对所述膨胀后的图像进行闭运算得到闭运算图;傅里叶变换模块,对所述开运算图和所述闭运算图进行傅里叶变换得到傅里叶变换值;腹部医学图像融合分类模块,用于对腹部CT扫描图像进行融合分类区分,并将区分结果传输给病灶检测部进行判断;病灶检测部,用于对所述傅里叶变换值结合腹部融合分类区分结果进行判断得到所述患者的疑似病灶区域;所述腹部医学图像融合分类模块包括以下步骤:
数据预处理步骤:获取所述患者的腹部CT扫描图像,进行包括去除所述腹部CT扫描图像中的基线和高频噪声在内的预处理步骤,并将其分为腹部CT扫描灰度解剖图像和伪彩色功能腹部CT扫描图像,分别标记为
和
多尺度分解模块:在多尺度分解过程中,通过改变传播滤波器中高斯卷积核的大小,对两幅源图像
和
进行平滑处理,使用平滑处理对信号
和
进行多尺度分解,获得平滑图像和细节图像;对于平滑图像S采用信息熵进行加权融合得到融合平滑图像,而对于细节图像采用多特征进行融合得到融合细节图像F
S;传播滤波器的原型为:S
p=1/Z
p∑
t∈N(p)w
p,t×I
t,其中,S
p为平滑后的图像,Z
p为归一化因子,N(p)为以p为中心像素点的邻域,i表示像素点序列号,I
t为在坐标位置t的图像像素值,w
p,t为传播滤波器权值;
其中,传播滤波器权值的第一项
表示利用与像素p和像素t邻接的像素之间的距离
和高斯卷积核σ
a的高斯函数;传播滤波器权值的第二项
表示利用与像素p邻接的所有像素距离和高斯卷积核σ
r的高斯函数;
重构融合步骤:利用像素级上的简单运算对融合后的平滑图像F D和细节图像F S进行基于ADMM重构得到融合图像F;输出最终的融合图像F;
分类步骤:并采用改进的SVM算法对融合图像F进行分类,SVM模型如下:
其中β∈R
P表示β是P维的模型向量,x
n表示样例,y
n表示类标签,N表示样本数,n∈{1,…,N}表示样本的索引,P表示特征数,p∈{1,…,P}表示特征的索引,β
p表示第p个特征对应的模型向量分量,f(x
n)表示分类面方程,
表示MCP罚函数,φ(β)表示基于相关系数的罚函数,1(y
n,f(x
n))表示铰链损失函数,λ
1、λ
2和γ表示可调参数。
进一步的,所述多尺度分解模块对于平滑图像S采用信息熵进行融合得到平滑图像F D采用的是香农提出的信息熵进行权值计算,其中信息熵函数主要利用像素值的概率分布来进行权重的排序,同一个值的像素总数越多表示该像素越重要即赋予该像素较大的权重,反之则赋予该像素较小的权重。
进一步的,所述预处理模块包括图像平滑去噪、图像锐化、图像背景分割以及图像二值化步骤。
本发明的有益效果
本发明采用滤波器和多特征方法进行MRI和SPECT两种不同模态的医学图像融合,利用不同高斯卷积核的传播滤波器将输入图像进行多尺度分解,该方法对图像进行平滑处理,提高平滑图像对噪声的鲁棒性。相对于传统的融合方法,利用亮度、方向以及相位特征来构造融合图像的权值,能获取输入图像更多的重要信息,进而为医生提供更精确的辅助诊疗信息。本发明最终构造的SVM分类器可以在保证组特征选择准确性的同时,具备oracle性质和组效应性质,可以实现更为有效的特征选择和类标签分类。可以排除一些不必要的训练数据特征,关注于关键的训练数据,并且按照不同的主题过滤策略对数据进行过滤,把所有分类错误的可能性都考虑进去,很大程度上降低误判的风险。
附图说明
图1是本发明提供优选实施例多尺度多特征的腹部医学图像融合及分类方法示意图;
图2是基于腹部医学图像融合分类的检测装置示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、详细地描述。所描述的实施例仅仅是本发明的一部分实施例。
本发明解决上述技术问题的技术方案是:如图2所示为一种基于腹部CT医学图像融合分类的检测装置,用于处理患者的腹部CT医学图像而得到所述患者的疑似病灶区域,其包括:获取模块,获取所述患者的腹部CT扫描图像,启动CT扫描器对人体的腹部进行电影模式CT扫描,并将所获得的每一张CT图片按人体进入扫描区的先后顺序给CT图像的每一图片赋予床位号和层位号,按扫描成像的时间顺序给CT图像的每一图片赋予相位号;将CT图像的每一图片进行阈值分割并获取人体体表轮廓曲线;以通过人体体表轮廓曲线在垂直方向的最低点的水平线为基准,对每一图片自左向右依次竖向逐点扫描,当图片的一列像素点的像素值全为零时记为零,当图片的一列由上而下有一个像素点的像素值不为零时,便记下第一个像素值不为零的像素点的高度值;然后将所记录的数据按自左向右的顺序排列,再以所得到的排列为行,按床位层位号的顺序由上向下排列,得到每一个床位的轮廓矩阵;将每一轮廓矩阵中所有元素的高度值累加求和,从一个床位中挑出所求和的最大值所对应的一个相位的一组图片,然后以该组图片为起点将相位号大于该组图片相位号的每一组图片按相位号由小到大排列,得到每一个床位内一个完整呼吸周期的图片阵列;预处理模块,对获取每一个床位内一个完整呼吸周期的图片阵列进行预处理得到二值化灰度图像;形态学腐蚀与膨胀模块,对所述二值化灰度图像进行形态学腐蚀运算得到腐蚀后的图像,对所述二值化灰度图像进行膨胀运算得到膨胀后的图像;运算模块,对所述腐蚀后的图像进行开运算得到开运算图,对所述膨胀后的图像进行闭运算得到闭运算图;傅里叶变换模块,对所述开运算图和所述闭运算图进行傅里叶变换得到傅里叶变换值;腹部医学图像融合分类模块,用于对腹部CT扫描图像进行融合分类区分,并将区分结果传输给病灶检测部进行判断;病灶检测部,用于对所述傅里叶变换值结合腹部融合分类区分结果进行判断得到所述患者的疑似病灶区域;如图1所示,所述腹部医学图像融合分类模块包括以下步骤:
数据预处理步骤:获取所述患者的腹部CT扫描图像,进行包括去除所述腹部CT扫描图像中的基线和高频噪声在内的预处理步骤,并将其分为腹部CT扫描灰度解剖图像和伪彩色功能腹部CT扫描图像,分别标记为
和
多尺度分解模块:在多尺度分解过程中,通过改变传播滤波器中高斯卷积核的大小,对两幅源图像
和
进行平滑处理,使用平滑处理对信号
和
进行多尺度分解,获得平滑图像和细节图像;对于平滑图像S采用信息熵进行加权融合得到融合平滑图像,而对于细节图像采用多特征进行融合得到融合细节图像F
S;传播滤波器的原型为:S
p=1/Z
p∑
t∈N(p)w
p,t×I
t,其中,S
p为平滑后的图像,Z
p为归一化因子,N(p)为以p为中心像素点的邻域,i表示像素点序列号,I
t为在坐标位置t的图像像素值,w
p,t为传播滤波器权值;
其中,传播滤波器权值的第一项
表示利用与像素p和像素t邻接的像素之间的距离
和高斯卷积核σ
a的高斯函数;传播滤波器权值的第二项
表示利用与像素p邻接的所有像素距离和高斯卷积核σ
r的高斯函数;
重构融合步骤:利用像素级上的简单运算对融合后的平滑图像F D和细节图像F S进行基于ADMM重构得到融合图像F;输出最终的融合图像F;
分类步骤:并采用改进的SVM算法对融合图像F进行分类,SVM模型如下:
其中β∈R
P表示β是P维的模型向量,x
n表示样例,y
n表示类标签,N表示样本数,n∈{1,…,N}表示样本的索引,P表示特征数,p∈{1,…,P}表示特征的索引,β
p表示第p个特征对应的模型向量分量,f(x
n)表示分类面方程,
表示MCP罚函数,φ(β)表示基于相关系数的罚函数,1(y
n,f(x
n))表示铰链损失函数,λ
1、λ
2和γ表示可调参数。
优选的,所述多尺度分解模块对于平滑图像S采用信息熵进行融合得到平滑图像FD采用的是香农提出的信息熵进行权值计算,其中信息熵函数主要利用像素值的概率分布来进行权重的排序,同一个值的像素总数越多表示该像素越重要即赋予该像素较大的权重,反之则赋予该像素较小的权重。
优选的,所述预处理模块包括图像平滑去噪、图像锐化、图像背景分割以及图像二值化步骤。
以上详细描述了本发明的较佳具体实施例。应当理解,本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此,凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案,皆应在由权利要求书所确定的保护范围内。