[发明专利]一种基于低秩分解和稀疏表示的多模态医学图像融合方法

说明书

本发明公开了一种基于低秩分解和稀疏表示的多模态医学图像融合方法，先将待融合的两幅不同多模态医学图像，进行低秩分解分别得到低秩部分图像和稀疏部分图像；利用KSVD算法对选用的非医学图像集训练低秩字典，利用KSVD算法对选用的非医学图像集进行低秩分解后得到的稀疏部分图像集训练稀疏字典；利用稀疏表示方法对低秩部分图像和稀疏部分图像进行稀疏重构，分别得到低秩重构图像和稀疏重构图像；利用稀疏表示方法对低秩重构图像和稀疏重构图像进行稀疏融合，得到融合图像；计算两幅不同多模态医学图像与稀疏重构图像、低秩重构图像之间的差值；将差值加入到融合图像中，得到最终稀疏融合图像。本发明在主观和客观评价指标上都优于传统的融合方法。

技术领域

本发明涉及一种基于低秩分解和稀疏表示的多模态医学图像融合方法，属于图像处理领域。

背景技术

在图像处理领域，图像融合是一项很有发展前景的研究。图像融合技术通过综合同一场景的多传感器图像或者不同场景的同一传感器图像的有用信息来合成一幅融合图像。而合成的融合图像具有前者所有的特征信息，更适合用于后期的处理和研究。一个高效的融合方法可以根据实际需求处理多通道信息。这些优势使得图像融合在很多领域里都特别受欢迎。尤其，在医学影像学为临床提供了超声图像、X射线、电子计算机体层扫描(CT)、磁共振成像(MRI)、数字减影成像(DSA)、正电子发射体层扫描(PET)、单光子发射断层成像(SPECT)等多种模态影像信息。不同的医学影像可以提供人体相关脏器和组织的不同信息，如CT和MRI提供解剖结构信息，而PET和SPECT提供功能信息在实际临床应用中，单一模态图像往往不能提供医生所需要的足够信息，通常需要将不同模态图像融合在一起，得到更丰富的信息以便了解病变组织或器官的综合信息，从而做出准确的诊断或制订出合适的治疗方案。例如，CT利用各种组织器官对X射线吸收系数的不同和计算机断层技术对人体进行成像，它对于骨、软组织和血管的组合成像效果很好，而对软组织则近乎无能为力。MRI利用水质子信息成像，对软组织和血管的显像灵敏度比CT高得多，但对骨组织则几乎不显像。由此可见不同成像技术对人体同一解剖结构所得到的形态和功能信息是互为差异、互为补充的。因此对不同影像信息进行适当的集成便成为临床医生诊断和治疗疾病的尤其重要。

随着图像融合研究的不断深入，很多的研究者提出了很多的图像融合方法。传统的图像融合方法有：线性加权法，它是一种最简单的图像融合方法，它直接对多幅原图像的对应像素点进行加权叠加。IIPF法（高通滤波法），它用高通滤波器算子提取出高分辨率图像的细节信息，然后采用像数相加的方法，将提取出的细节信息叠加到低分辨率图像上，实现低分辨率图像和高分辨率图像之间的数据融合。PCA（主成分分析）法，首先对图像进行主分量变换，由相关矩阵求特征值和特征向量求得各主分量，通过这种融合，可以尽可能地保留全色图像的细节信息，然后对融合后的图像进行反变换，得到包含丰富细节信息的融合图像。IIIS变换法，它是遥感图像融合中比较常用到的一种方法，它将多光谱的TM图像数据变换到IIIS空间，得到3个独立的分量。它首先将SPOT图像数据经过适当的伸缩，以保证它的均值和方差变换到和I分量一致，然后用伸缩后的数据替换I分量，得到II分量，再分别对II和S进行反差扩展，得到III和SI，最后将II，III和SI变换回RGB空间，得到融合图像。和以上四类算法相比，其余四种融合方法功能比较单一。

发明内容

本发明提供了一种基于低秩分解和稀疏表示的多模态医学图像融合方法，以用于对医学图像进行融合，以获得更好融合质量图像。

本发明的技术方案是：一种基于低秩分解和稀疏表示的多模态医学图像融合方法，所述方法的具体步骤如下：

Step1、将待融合的两幅不同多模态医学图像A、B，使用矩阵低秩分解理论进行低秩分解分别得到低秩部分图像A₁、B₁和稀疏部分图像A₂、B₂；