本发明属于医疗图像处理技术领域，特别涉及一种多尺度多特征的心电图医学图像融合及分类方法。
背景技术：
现有心电图医学图像融合技术主要对两种不同模态的医学图像，根据医学图像模态的不同，医学图像融合系统可以分为三种类型：解剖医学图像与解剖医学图像融合、解剖医学图像与功能医学图像融合以及功能医学图像与功能医学图像融合。MRI-PET和MRI-SPECT医学图像融合系统属于解剖医学图像与功能医学图像融合，该系统的输入图像是灰度和伪彩色。飞利浦公司推出的MRI-PET组合一体化套件将商用的MRI成像扫描仪和具有特殊屏蔽的PET结合起来，得到的图像对癌症转移的诊断及术前分期具有重要的临床价值。
由于像素级多尺度融合方法在不同尺度图像上直接对图像的像素值进行处理，能最大限度地保留输入图像的像素信息和提高融合图像的质量。因此，像素级多尺度方法在医学图像融合领域成为一个研究热点。目前，国内外主要从图像分解与重构方法和图像融合规则方法两个方面出发，提出新的像素级多尺度融合方法。在图像分解与重构方法方面，基于频域的图像融合方法利用傅里叶变换和傅里叶逆变换对图像信号进行不同尺度的分解与重构，但是该类方法具有时间复杂度高和运行时间长的特点，这与高实时性的医疗辅助诊断矛盾，并且对处理的实验平台硬件设施和软件设施要求非常高。后来，研究者们提出利用空域滤波器对图像进行处理来对输入图像进行多尺度分解与重构，该类方法虽然能快速地进行图像分解与重构，但是基于空域的图像融合方法抗噪性差。但是现有方法针对不同模态的医学图像采用同一种特征进行融合。
技术实现要素：
本发明旨在解决以上现有技术的问题。提出了一种降低复杂度、提高融合图像的质量、提高分类精确度的多尺度多特征的心电图医学图像融合及分类方法。
本发明的技术方案如下：
一种多尺度多特征的心电图医学图像融合及分类方法，其包括以下步骤：步骤1：获取数据库中所有的心电图医学图像，进行包括去除所述心电图医学图像中的基线和高频噪声在内的预处理步骤，并将其分为灰度解剖心电图医学图像和伪彩色功能心电图医学图像，分别标记为和
步骤2、在多尺度分解过程中，通过改变传播滤波器中高斯卷积核的大小，对对步骤1的两幅源图像和进行平滑处理，使用平滑处理对图像和进行多尺度分解，获得平滑图像和细节图像；对于平滑图像S采用信息熵进行加权融合得到融合平滑图像，而对于细节图像采用多特征进行融合得到融合细节图像FS；传播滤波器的原型为：Sp＝1/Zp∑t∈N(p)wp,t×It，其中，Sp为平滑后的图像，Zp为归一化因子，N(p)为以p为中心像素点的邻域，i表示像素点序列号，It为在坐标位置t的图像像素值，wp,t为传播滤波器权值；其中，传播滤波器权值的第一项表示利用与像素p和像素t邻接的像素之间的距离和高斯卷积核σa的高斯函数；传播滤波器权值的第二项表示利用与像素p邻接的所有像素距离和高斯卷积核σr的高斯函数；
步骤3：利用像素级上的简单运算对融合后的平滑图像FD和细节图像FS进行基于ADMM重构得到融合图像F；输出最终的融合图像F；
步骤4：并采用Adaboos分类器融合图像F进行分类，分类步骤包括：通过第一训练集训练所述Adaboost分类器训练时需要用到的SVM分类器；使用训练出的SVM分类器提取第二训练集的样本特征；根据提取出的第二训练集的样本特征训练Adaboost分类器。
进一步的，所述步骤3对于平滑图像S采用信息熵进行融合得到平滑图像FD采用的是香农提出的信息熵进行权值计算，其中信息熵函数主要利用像素值的概率分布来进行权重的排序，同一个值的像素总数越多表示该像素越重要即赋予该像素较大的权重，反之则赋予该像素较小的权重。
进一步的，所述通过第一训练集训练所述Adaboost分类器训练时需要用到的SVM分类器，包括：
2-1)选择第一目标窗口的一部分区域，作为小块目标窗口，所述第一目标窗口与所述第一训练集中的正样本或负样本一样大小；
2-2)从所述第一训练集提取所述小块目标窗口对应的样本特征；
2-3)根据所述小块目标窗口对应的样本特征训练所述小块目标窗口的SVM分类器；
2-4)使用训练出的所述小块目标窗口的SVM分类器检测负样本原始图片集，将识别为正样本的目标区域截取出来投入到所述第一训练集以更新所述第一训练集；
2-5)返回步骤2-2)，重复执行预定的训练次数，训练出所述小块目标窗口的SVM分类器；
2-6)遍历不同大小和不同位置的所有小块目标窗口，按步骤2-1)～2-5)训练所述所有小块目标窗口的SVM分类器；
2-7)使用训练出的所述所有小块目标窗口的SVM分类器检测测试库中的测试图片，得到所述所有小块目标窗口的SVM分类器的正检率；
2-8)按照正检率从高到低对所述所有小块目标窗口的SVM分类器进行排序，将窗口重合面积超过各自预定比例的任意两个窗口的SVM分类器进行合并，保留其中正检率大的SVM分类器；
2-9)按照正检率从高到低对合并后的所有的SVM分类器进行排序，选择排在前面的预定个数的SVM分类器。
进一步的，所述步骤1去除所述心电图医学图像中的基线和高频噪声在内的预处理步骤，具体包括：采用200ms和500ms带宽的中值滤波器去除ORS复合波和P&T波，并从原始信号中减去，即得到去除基线的信号图像。
本发明的有益效果
本发明采用滤波器和多特征方法进行MRI和SPECT两种不同模态的医学图像融合，利用不同高斯卷积核的传播滤波器将输入图像进行多尺度分解，该方法对图像进行平滑处理，提高平滑图像对噪声的鲁棒性。相对于传统的融合方法，利用亮度、方向以及相位特征来构造融合图像的权值，能获取输入图像更多的重要信息，进而为医生提供更精确的辅助诊疗信息。本发明提供的Adaboost分类器训练方法，通过第一训练集训练所述Adaboost分类器训练时需要用到的SVM分类器；使用训练出的SVM分类器提取第二训练集的样本特征；根据提取出的第二训练集的样本特征训练Adaboost分类器。与现有技术相比，本发明通过预先训练好的SVM分类器映射获得训练样本的样本特征，避免了大量计算，能够缩短Adaboost分类器的训练时间，同时结合了SVM分类器训练精度高的优点，解决了Adaboost分类器在小样本情况下训练效果不好的问题。进一步地，本发明训练出的Adaboost分类器进行目标检测时，检测速度更快，同时正检率提高，误检率下降，提高了检测效果。
附图说明
图1是本发明提供优选实施例多尺度多特征的心电图医学图像融合及分类方法示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、详细地描述。所描述的实施例仅仅是本发明的一部分实施例。
本发明解决上述技术问题的技术方案是：
如图1所示为一种多尺度多特征的心电图医学图像融合及分类方法，其包括以下步骤：
步骤1：获取数据库中所有的心电图医学图像，进行包括去除所述心电图医学图像中的基线和高频噪声在内的预处理步骤，并将其分为灰度解剖心电图医学图像和伪彩色功能心电图医学图像，分别标记为和
步骤2、在多尺度分解过程中，通过改变传播滤波器中高斯卷积核的大小，对对步骤1的两幅源图像和进行平滑处理，使用平滑处理对图像和进行多尺度分解，获得平滑图像和细节图像；对于平滑图像S采用信息熵进行加权融合得到融合平滑图像，而对于细节图像采用多特征进行融合得到融合细节图像FS；传播滤波器的原型为：Sp＝1/Zp∑t∈N(p)wp,t×It，其中，Sp为平滑后的图像，Zp为归一化因子，N(p)为以p为中心像素点的邻域，i表示像素点序列号，It为在坐标位置t的图像像素值，wp,t为传播滤波器权值；其中，传播滤波器权值的第一项表示利用与像素p和像素t邻接的像素之间的距离和高斯卷积核σa的高斯函数；传播滤波器权值的第二项表示利用与像素p邻接的所有像素距离和高斯卷积核σr的高斯函数；
步骤3：利用像素级上的简单运算对融合后的平滑图像FD和细节图像FS进行基于ADMM重构得到融合图像F；输出最终的融合图像F；
步骤4：并采用Adaboos分类器融合图像F进行分类，分类步骤包括：通过第一训练集训练所述Adaboost分类器训练时需要用到的SVM分类器；使用训练出的SVM分类器提取第二训练集的样本特征；根据提取出的第二训练集的样本特征训练Adaboost分类器。
优选的，所述步骤3对于平滑图像S采用信息熵进行融合得到平滑图像FD采用的是香农提出的信息熵进行权值计算，其中信息熵函数主要利用像素值的概率分布来进行权重的排序，同一个值的像素总数越多表示该像素越重要即赋予该像素较大的权重，反之则赋予该像素较小的权重。
优选的，所述通过第一训练集训练所述Adaboost分类器训练时需要用到的SVM分类器，包括：
2-1)选择第一目标窗口的一部分区域，作为小块目标窗口，所述第一目标窗口与所述第一训练集中的正样本或负样本一样大小；
2-2)从所述第一训练集提取所述小块目标窗口对应的样本特征；
2-3)根据所述小块目标窗口对应的样本特征训练所述小块目标窗口的SVM分类器；
2-4)使用训练出的所述小块目标窗口的SVM分类器检测负样本原始图片集，将识别为正样本的目标区域截取出来投入到所述第一训练集以更新所述第一训练集；
2-5)返回步骤2-2)，重复执行预定的训练次数，训练出所述小块目标窗口的SVM分类器；
2-6)遍历不同大小和不同位置的所有小块目标窗口，按步骤2-1)～2-5)训练所述所有小块目标窗口的SVM分类器；
2-7)使用训练出的所述所有小块目标窗口的SVM分类器检测测试库中的测试图片，得到所述所有小块目标窗口的SVM分类器的正检率；
2-8)按照正检率从高到低对所述所有小块目标窗口的SVM分类器进行排序，将窗口重合面积超过各自预定比例的任意两个窗口的SVM分类器进行合并，保留其中正检率大的SVM分类器；
2-9)按照正检率从高到低对合并后的所有的SVM分类器进行排序，选择排在前面的预定个数的SVM分类器。
优选的，所述步骤1)去除所述心电图医学图像中的基线和高频噪声在内的预处理步骤，具体包括：采用200ms和500ms带宽的中值滤波器去除ORS复合波和P&T波，并从原始信号中减去，即得到去除基线的信号图像。以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。