基于多尺度多特征的影像科用脑部图像融合识别方法技术

本发明专利技术公开了一种基于多尺度多特征的影像科用脑部图像融合识别方法，其包括以下步骤：采集原始脑部医学图像作为样本，使用加权灰度化算法对原始脑部医学图像进行灰度化处理，得到灰度图像；使用直方图均衡化对灰度图像进行处理，得到均衡后灰度图，采使用改进的Isotropic Sobel边缘检测算子对均衡后灰度图进行边缘检测，得到边缘脑部图像，并得到灰度解剖医学图像MRI和伪彩色功能医学图像，初始化卷积神经网络，训练数据，提取出特征数据；对于经过医学图像识别训练模型训练后的图像采用多尺度分解得到平滑图像和细节图像，进行重构得到融合图像，将融合图像输入脑部医将卷积神经网络模型训练提取的特征数据传到支持向量机中进行训练，将提取的测试特征数据输入医学图像识别训练模型中进行判断，最终得到准确的脑部医学图像识别结果。

【技术实现步骤摘要】
基于多尺度多特征的影像科用脑部图像融合识别方法
本专利技术涉及医疗影像科用图像处理领域，具体涉及一种基于多尺度多特征的影像科用脑部图像融合识别方法。
技术介绍
目前，国内外诸多著名学者都致力于图像分割融合算法的研究，图像分割方法主要有以下四种：基于阈值、基于边缘检测、基于区域、基于能量的分割。基于阈值的分割方法主要有直方图凹面分析法、最大类间方差法、阈值插值法等，该类方法直观简单高效，但是由于图像的复杂性，阈值的选取成为该类方法的一大挑战；基于边缘检测的分割方法中比较经典的算法有Sobel、Prewitt、Laplace以及Canny算子等，该类方法由于噪声点和周围像素点间的阶跃性非常明显，所以极易被误判为边缘；基于区域的图像分割方法主要有区域生长法以及分裂合并法，区域生长法分割较大图像容易出现不连续的空洞从而引发过分割，分裂合并法在不断分裂的过程中容易对边界区域产生破坏；基于能量的图像分割方法包括基于水平集的方法、基于图论的方法、基于ICM的方法等等。基于水平集的方法自出现以来，便成为图像分割领域的热点，在国际顶级期刊和国际会议上均有大量的水平集图像分割新方法提出。Caselles和Malladi等人在图像分割的主动轮廓模型的上下文中介绍了水平集方法，但分割结果存在误差，而且分割结果不稳定。为了解决这个问题，Osher等人采用重新初始化水平集函数的方法来确保水平集的稳定性。然而，重新初始化的做法不仅会引起严重的问题，而且还会影响数字的准确性。之后，李纯明提出了一种新的变分水平集。在其数值实现中，可以在有限差分方法中使用相对较大的时间步长以减少迭代次数，同时确保足够的数值精度。王晓峰提出了一种高效的、鲁棒的水平集方法，将多尺度分段思想引入局部区域，其具有很好的效率和鲁棒性。现有医学图像融合技术主要对两种不同模态的医学图像，根据医学图像模态的不同，医学图像融合系统可以分为三种类型：解剖医学图像与解剖医学图像融合、解剖医学图像与功能医学图像融合以及功能医学图像与功能医学图像融合。MRI-PET和MRI-SPECT医学图像融合系统属于解剖医学图像与功能医学图像融合，该系统的输入图像是灰度和伪彩色。飞利浦公司推出的MRI-PET组合一体化套件将商用的MRI成像扫描仪和具有特殊屏蔽的PET结合起来，得到的图像对癌症转移的诊断及术前分期具有重要的临床价值。但是目前针对于脑部医学图像的处理分割因为其复杂性和变化性，导致识别能力较低，对于医生的检查和治疗没有起到很好的识别性能，有时候往往一些微小的变化未能及时观察出来，导致严重的后果。因此有必要提供一种精确的脑部医学图像识别方法。
技术实现思路
有鉴于现有技术的上述缺陷，本专利技术所要解决的技术问题是提供一种融合度高、识别度大大提高、精确度高的基于多尺度多特征的影像科用脑部图像融合识别方法。本专利技术技术方案如下：一种基于多尺度多特征的影像科用脑部图像融合识别方法，包括以下步骤：采集原始脑部医学图像作为样本，使用加权灰度化算法对原始脑部医学图像进行灰度化处理，得到灰度图像；使用直方图均衡化对灰度图像进行处理，得到均衡后灰度图，采用高次样条函数对灰度直方图进行拟合，所述拟合后的曲线有明显的谷值点和峰值点；灰度直方图划分谷值区间，在高次样条函数拟合的基础上得到平滑可导的拟合曲线，并求取该平滑可导拟合曲线的极值点，根据极值点两边的符号将谷值点筛选出来；使用改进的IsotropicSobel边缘检测算子对均衡后灰度图进行边缘检测，得到边缘脑部图像，并得到灰度解剖医学图像MRI和伪彩色功能医学图像，对拟合改进的IsotropicSobel边缘检测算子为所示：使用自适应阈值算法对边缘图像进行二值化，得到二值化后的脑部医学图像；使用形态学运算操作对二值化后的医学图像进行处理，得到脑部医学候选区域图像，并将候选区域图像形成训练数据，初始化卷积神经网络，所述初始化卷积神经网络为：设置卷积神经网络中的参数，其中包括：卷积核的数量、降采样层的数量、卷积核的大小、降采样层的降幅，初始化卷积核的权重和偏置；将训练数据分批均分到输入卷积神经网络中，训练数据分别经过卷积层、降采样层、卷积层、降采样层、多层感知器，完成前向传播；对多层感知器进行误差计算和梯度计算，并判断误差是否收敛；若是，则将得到的误差和梯度用反向传播算法，经过降采样层、卷积层、降采样层、卷积层、输入层逐层传播，并且逐层更新网络的权重，判断是否为输入层，若是则提取出特征数据；将卷积神经网络模型训练提取的特征数据传到支持向量机中进行训练，把经过卷积神经网络的训练特征数据输入支持向量机，同时，用网格搜索的优化方法来优化支持向量机的参数C和δ，确定最优的支持向量机模型，建立医学图像识别训练模型；对于经过医学图像识别训练模型训练后的图像采用多尺度分解得到平滑图像和细节图像，平滑图像采用信息熵进行融合得到平滑图像FD，而对于细节图像D采用多特征进行融合得到细节图像FS，将平滑图像FD和细节图像FS进行重构得到融合图像，将融合图像输入脑部医学图像识别训练模型中进行判断，最终得到准确的医学图像融合识别结果。进一步的，所述将训练数据分批均分到输入卷积神经网络中，训练数据分别经过卷积层、降采样层、卷积层、降采样层、多层感知器，完成前向传播，具体包括：首先从样本集中取一批样本(X，YP)，其中X是样本数字的向量，Y是X对应的期望值，P是0到9的数字，将X输入卷积神经网络，计算出相应的实际输出OP，OP＝Fn(...F2(F1(XPW(1))W(2))W(n))，n是卷积神经网络的第n层，W表示权值，其中卷积运算是用卷积滤波器在上层网络结构中做卷积运算，然后进行非线性变换，而降采样运算只采用最大池化操作，即最大池采样是通过一个滤波器提取上层网络结构的特征数据，不经过非线性运算，每次滤波后的最大值是数据降采样后的一个特征。进一步的，所述反向传播算法具体为：按极小化误差的方法反向传播并且调整卷积神经网络中的权矩阵，首先对样本批量前向传播，计算出卷积神经网络中所有的激活值；然后，针对每层节点，计算其残差，残差是从后向前的求导过程；接着，计算权值的偏导数，并更新权值参数；最后，重复以上方法迭代卷积神经网络参数使代价函数收敛到一个极小值，最终求解得到卷积神经网络模型。有益效果：本专利技术一种融合度高、识别度大大提高、精确度高的基于多尺度多特征的影像科用脑部图像融合识别方法，该方法将卷积神经网络模型和支持向量机模型进行有机结合。结合卷积神经网络和支持向量机的识别模型很深度的描述样本数据和期望数据的相关性，脑部医学图像由于其特殊性、复杂性，需要融合后才能识别，这样实现精确度大大提高，之前并没有任何相关报道将二者融合实现识别，本专利技术对数字模式分类的区分能力很强。而且卷积神经网络模型和支持向量机模型它们的目标是区分性的，使生成的脑部医学图像识别系统的输出更易优化。通过使用加权灰度化算法对原始脑部医学图像进行灰度化处理，得到灰度图像，这样脑部医学图像就被形象的处理成为容易识别的图像；再使用直方图均衡化对灰度图像进行处理，得到均衡后灰度图，使用改进的IsotropicSobel边缘检测算子对均衡后灰度图进行边缘检测，得到边缘图像相比于原来的边缘检测算子，使其边缘检测算子归一化，能有本文档来自技高网...

【技术保护点】
1.一种基于多尺度多特征的影像科用脑部医学图像融合识别方法，其特征在于，包括以下步骤：采集原始脑部医学图像作为样本，使用加权灰度化算法对原始脑部医学图像进行灰度化处理，得到灰度图像；使用直方图均衡化对灰度图像进行处理，得到均衡后灰度图，采用高次样条函数对灰度直方图进行拟合，所述拟合后的曲线有明显的谷值点和峰值点；灰度直方图划分谷值区间，在高次样条函数拟合的基础上得到平滑可导的拟合曲线，并求取该平滑可导拟合曲线的极值点，根据极值点两边的符号将谷值点筛选出来；使用改进的Isotropic Sobel边缘检测算子对均衡后灰度图进行边缘检测，得到边缘脑部图像，并得到灰度解剖医学图像MRI和伪彩色功能医学图像，对拟合改进的Isotropic Sobel边缘检测算子为所示：

【技术特征摘要】
1.一种基于多尺度多特征的影像科用脑部医学图像融合识别方法，其特征在于，包括以下步骤：采集原始脑部医学图像作为样本，使用加权灰度化算法对原始脑部医学图像进行灰度化处理，得到灰度图像；使用直方图均衡化对灰度图像进行处理，得到均衡后灰度图，采用高次样条函数对灰度直方图进行拟合，所述拟合后的曲线有明显的谷值点和峰值点；灰度直方图划分谷值区间，在高次样条函数拟合的基础上得到平滑可导的拟合曲线，并求取该平滑可导拟合曲线的极值点，根据极值点两边的符号将谷值点筛选出来；使用改进的IsotropicSobel边缘检测算子对均衡后灰度图进行边缘检测，得到边缘脑部图像，并得到灰度解剖医学图像MRI和伪彩色功能医学图像，对拟合改进的IsotropicSobel边缘检测算子为所示：使用自适应阈值算法对边缘图像进行二值化，得到二值化后的脑部医学图像；使用形态学运算操作对二值化后的医学图像进行处理，得到脑部医学候选区域图像，并将候选区域图像形成训练数据，初始化卷积神经网络，所述初始化卷积神经网络为：设置卷积神经网络中的参数，其中包括：卷积核的数量、降采样层的数量、卷积核的大小、降采样层的降幅，初始化卷积核的权重和偏置；将训练数据分批均分到输入卷积神经网络中，训练数据分别经过卷积层、降采样层、卷积层、降采样层、多层感知器，完成前向传播；对多层感知器进行误差计算和梯度计算，并判断误差是否收敛；若是，则将得到的误差和梯度用反向传播算法，经过降采样层、卷积层、降采样层、卷积层、输入层逐层传播，并且逐层更新网络的权重，判断是否为输入层，若是则提取出特征数据；将卷积神经网络模型训练提取的特征数据传到支持向量机中进行训练，把经过卷积神经网络的训练特征数据输入支持向量机，同时，用网格搜索的优化方法来优化支持向量机...

【专利技术属性】
技术研发人员：唐天才，
申请(专利权)人：李家菊，
类型：发明
国别省市：贵州,52