一种医学图像融合方法技术

本发明专利技术涉及一种医学图像融合方法，属图像处理领域。本发明专利技术的方法基于NSCT多尺度分解和低秩表示方法，经过图像预处理、图像配准、图像多尺度分解、低秩表示融合等步骤实现。理论与实践证明，本发明专利技术妥善的解决了图像融合中的轮廓不清晰，细节信息丢失等问题。采用本文发明专利技术的方法后，不仅可以高效完成医学图像的融合，而不丢失重要的图像的融合信息，因此可以获得高质量的医学图像融合图像，满足临床医生的诊断阅读需求。

Medical image fusion method

The invention relates to a medical image fusion method, belonging to the image processing field. The method of the invention is based on the NSCT multi-scale decomposition and the low rank representation method, and realizes the steps of image preprocessing, image registration, multi-scale image decomposition and low rank representation fusion. The theory and the practice prove that the invention solves the problems such as the unclear outline and the loss of detail information in the image fusion. Using this method the invention, not only can efficiently complete fusion of medical images, without loss of important image information fusion, so the fusion image is obtained high quality medical images can meet the reading needs of clinicians, diagnosis.

【技术实现步骤摘要】
一种医学图像融合方法
本专利技术涉及医学图像领域，具体涉及基于低秩表示和多尺度分解的医学图像融合。
技术介绍
近年来，随着计算机技术、工业技术以及医学的快速发展，医学影像学为临床医学诊断提供的图像模态越来越多样化。由于成像原理的不同，多模态医学图像使得不同模态的医学影像具备不同的特点和优势，因此将不同特点和优势的图像想融合减少临床诊断的错误率。因此如何提高融合的质量是目前医学融合领域的难点。根据被融合图像成像方式不同，可分为同类方式融合和交互方式融合。图像融合的方式为单模融合、多模融合、模板融合，医学图像分为解剖图像和功能图像。而单模融合是指相同成像方式的图像融合。多模融合是指将形态成像和功能成像技术所获得的图像进行融合，即不同成像方式所获得图像间的融合；模板融合是指利用从健康人的研究中建立一系列模板为标准，将病人的图像与模板图像的融合。如SPECT与MR图像融合，PET与CT图像融合等。现有的基于图像特征的融合方法几乎都是从变换域上的图像编码和压缩技术延伸来的，有Laplacian金字塔法、Gaussian金字塔法、比率低通金字塔法、多分辨率形态滤波法和小波变换法等。现阶段，我国张鑫等提出的Contourlet变换系数加权的医学图像融合，该方法首先是将源图像进过Contourlet变换分解获得不同尺度对个方向下的分解系数。然后再对该系数进行分析来确定融合规则。这个算法可以增加多模态医学图像互补信息，并能够较好地提高医学图像融合的清晰度。Li提出一种单纯形变权的MRI和SPECT的融合方法，是基于GIHS空间，通过非下采样轮廓变换(NSCT)与变权相结合来避免了在MRI和SPECT领域的传统的透明技术的不足。还有一些学者提出了其他的融合方法，例如徐志平于2014提出通过局部极大值的方法来进行融合，该方法保存了图像的细节部分而且极大的提高了融合的质量和效率。SabalanDaneshvar于2014年提出了通过视觉模型和IHS空间的MRI和PET的医学图像融合方法，IHS算法和RIM融合技术能够保存更多的空间特征以及更有效的信息。KavithaC.T等人提出基于混合智能的医学图像融合该方法是将原始图像即PET或者SPECT通过蚁群算法进行边缘检测盒加强，经过这些步骤后，再进过s-PCNN的方法进行处理最后根据融合规则得到最后的同和结果。归纳起来，上述现有的技术只是在考虑在多尺度分解的基础上进行图像融合，但是没有考虑到分频之后系数进行最大限度的进行的利用，融合结果图不能细节信息丢失等问题。
技术实现思路
本专利技术的目的在于如何提高MRI和PET医学图像融合的质量，在低秩表示和多尺度分解的基础上提出了MRI和PET医学图像融合。此方法具有良好的融合效果，对Dicom文件的融合具有更加优异成果，程序具有时间短，能和多种多尺度分解的方法相结合进行融合来客服低秩表示融合方法和基于多尺度分解的图像融合的缺点，具有更好的鲁棒性。本专利技术采用的技术方案如下：该方法分为初始和融合两个阶段，在这两个阶段分别是预处理阶段和融合阶段。预处理阶段是对图像先进行配准然后将配准后的图像进行归一化的处理。融合阶段则是对预处理后的图像进行基于低秩表示的融合算法进行融合得到最终的融合的结果图。具体是依序采用如下步骤：首先对MRI和PET源图像分别记为A，B分别进行归一化处理，得到的结果从RGB空间变换到IHS空间得到强度分量IA，IB和各自的色度和饱和度；得到的IA，IB分量通过多尺度分解，对低频系数使用低秩表示来分成一个稀疏误差项(EA，EB)和一个低秩部分(ZA，ZB)；然后分别对稀疏误差项EA，EB和低秩部分(ZA，ZB)采用1泛数选大法和系数绝对值最大法进行融合得到新的稀疏误差项(EF)和新的低秩部分(ZF)，将这两项重构成新的低频融合系数。高频系数则采用系数绝对值选大法进行融合。对新的分量IF和PET图像的色度和饱和度通过反IHS变换得到最终的融合结果图。本专利技术提出了一种基于低秩表示的医学图像融合方法。该方法是将RGB颜色空间变换到IHS空间是为了更加满足人眼的需求，然后再通过低秩表示的算法来表示特征，Z是反应了图像的主要特征，E是表示了图像的一些辨识特征和显著特征。其有益效果是：1、本专利技术与传统的医学图像融合方法相比，由于低秩模型能够保证了图像的特征信息，因此保证了图像更多的细节信息。2、本专利技术低秩表示方法能和多尺度分解(如：LP，Conturlet，NSCT)相结合来客服多尺度的固有的缺点。3、本专利技术不仅保持了图像更多的细节信息而且运行时间也相比其他方法较好。附图说明下面结合附图和具体实施方式对本专利技术进一步详细说明。图1是本专利技术所述基于低秩表示和多尺度分解MRI和PET的医学图像融合方法流程示意图。图2是本专利技术基于低秩表示的低频融合流程示意图。图3是本专利技术的融合结果，(a)，(b)为MRI和PET的融合结果图，(c)-(h)分别为GIHS，WT，HPF，NSCT，NSCT_PCNN，GAUSSIAN融合方法的融合结果图，(i)为本文提出的算法的融合结果图。具体实施方式以下结合附图对专利技术的实施做出进一步说明。图1为本专利技术提出的基于低秩表示的医学图像融合算法流程。首先将配准后的源图像MRI和PET图像进行归一化后得到初始图像。第一步：根据IHS颜色模型，然后将源图像进行多尺度分解Setp1：将RGB空间转换到IHS空间。Setp2：把强度通过多尺度进行多尺度分解为高频系数和低频系数。Setp3：高频频带每层通过系数绝对值选大法进行融合。Setp4：与此同时低频部分通过低秩表示进行融合得到Lfusion。Step5：将高频融合系数和低频融合系数通过多尺度反变换得到新的融合图像的强度Step6：强IHS空间转换为RGB空间得到最终的融合结果图低频部分的融合算法流程如下：稀疏项(E)和低秩项(Z)计算：Step1：将LA，LB分别通过不重复的原则按4*4窗口进行选取。Step2：选取的区域通过列优先的原则进行排列得到新的向量XA，XB本专利技术的模型算法如下：s.t.X＝Z+EStep1：模型本身是NP难题，所以求解之前需要对其进行松弛变换，变换如下s.t.X＝Z+EStep2：构造拉格朗日函数：Step3：固定其他变量，然后分别对Z，E求其最小化，再更新Y更新算法如下：初始化：Z＝0，E＝0，Y＝0，μ＝10-6，ρ＝1.1，ε＝10-8，ε2＝10-2循环下面步骤直到收敛1.固定其他项更新Z2.固定其他项更新E3.固定其他项更新YYk+1＝Yk+μk(X-Zk+1-Ek+1)4.更新μ5.判断收敛||X-Z-E||∞＜ε2Step4：稀疏项(E)和低秩项(Z)的融合算法：稀疏项(E)采用绝对值1泛数进行融合：低秩项采系数S型权重函数进行融合W(ZA(i，j))＝f(ZA(i，j))W(ZB(i，j))＝f(ZB(i，j))Zfusion＝W(ZA)*ZA+W(ZB)*ZBStep5：通过低秩模型重构得到新的低频系数。Step6：将IHS空间通过逆变换到RGB空间得到融合结果。本专利技术通过结合多尺度和低秩表示的融合方法不仅克服了多尺度分解和低秩表示的固有缺点，同时也提高了图像融合效率以及融合的准确度。最后将本专利技术的结果图和传统的NSCT本文档来自技高网...

【技术保护点】
基于NSCT和低秩表示的医学图像融合方法，其特征在于：所述多尺度分解后的低频信息来获取低秩原矩阵(X)；针对低秩表示的针对图像的细节信息，低秩表示模型得到的稀疏项保证了图像的一些辨识特征和显著特征而低秩部分则反映了分解图像的主要特征，故使用低秩模型对低频系数进行融合。所述低秩表示的具体步骤为：1)将L

【技术特征摘要】
1.基于NSCT和低秩表示的医学图像融合方法，其特征在于：所述多尺度分解后的低频信息来获取低秩原矩阵(X)；针对低秩表示的针对图像的细节信息，低秩表示模型得到的稀疏项保证了图像的一些辨识特征和显著特征而低秩部分则反映了分解图像的主要特征，故使用低秩模型对低频系数进行融合。所述低秩表示的具体步骤为：1)将LA，LB分别通过不重复的原则按4*4窗口进行选取。2)选取的区域通过列优先的原则进行排列得到新的向量XA，XB。2.根据权利2中要求所述的低秩表示模型的图像融合方法步骤为：1)本发明的模型算法如下：s.t.X＝Z+E因为模型本身是NP难题，所...

【专利技术属性】
技术研发人员：刘哲，宋余庆，徐春艳，
申请(专利权)人：江苏大学，
类型：发明
国别省市：江苏,32