(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202210947992.3 (22)申请日 2022.08.09 (71)申请人 中科微影 （浙江） 医疗科技有限公司 地址 317317 浙江省台州市仙居县白塔镇 仙居县经济开发区白塔 工业集聚区 (72)发明人 张源　孙文凤　 (74)专利代理 机构 北京一枝笔知识产权代理事 务所(普通 合伙) 11791 专利代理师 张庆瑞 (51)Int.Cl. G06T 5/50(2006.01) G06V 10/80(2022.01) (54)发明名称 一种核磁共振图像与CT图像融合方法与系 统 (57)摘要 本发明提供了一种核磁共振图像与CT图像 融合方法与系统， 其中该方法包括： 对核磁共振 图像和CT 图像进行小波分解并计算分解后图像 的信息梯度值； 根据信息梯度值构建融合系数； 计算核磁共振小波图像和CT小波图像的信息值； 根据信息值构建融合公式； 根据融合公式和融合 系数构建图像融合模型； 利用图像融合模型对核 磁共振图像和CT 图像进行融合得到融合后的图 像。 本发明通过待融合图像的信息梯度值和 信息 值构建图像融合模型， 并基于此对核磁共振图像 和CT图像进行融合， 不仅可以改善图像的清晰 度， 可以弥补单一的核磁共振图像或者CT图像信 息量不足的缺陷， 大大提高了医生对图像的判读 效果， 为后续病灶定位 提供了有利的影 像依据。 权利要求书3页 说明书7页 附图2页 CN 115471432 A 2022.12.13 CN 115471432 A 1.一种核磁共 振图像与CT图像融合方法， 其特 征在于， 包括： 步骤1： 获取核磁共 振图像和CT图像； 步骤2： 对所述核磁共振图像和所述CT图像进行小波分解得到核磁共振小波图像和CT 小波图像； 步骤3： 根据所述核磁共 振小波图像和所述CT 小波图像得到信息梯度值； 步骤4： 根据信息梯度值构建融合系数； 步骤5： 计算所述核磁共 振小波图像和所述CT 小波图像的信息值； 步骤6： 根据所述信息值构建融合公式； 步骤7： 根据所述融合公式和融合系数构建图像融合模型； 步骤8： 利用所述图像融合模型对所述核磁共振图像和CT图像进行融合得到融合后的 图像。 2.如权利要求1所述的一种核磁共振图像与CT图像 融合方法， 其特征在于， 所述步骤3： 根据所述核磁共 振小波图像和所述CT 小波图像得到信息梯度值， 包括： 步骤3.1： 对所述核磁共振小波图像和所述CT小波图像进行叠加平均得到叠加平均图 像； 步骤3.2： 采用信息梯度值计算模型计算所述核磁共振小波图像、 所述CT小波图像和所 述叠加平均图像的信息梯度值； 其中， 信息梯度值计算模型为： 其中， 表示核磁共振小波图像在以点(i,j)为中心， 邻域窗口大小为N ×M区域 内的信息梯度值， 表示CT小波图像在以点(i,j)为中心， 邻域窗口大小为N ×M区 域内的信息梯度值， 表示叠加平均图像在以点(i,j)为中心， 邻域窗 口大小为N × M区域内的信息梯度值， 表示核磁共振小波图像在邻域窗口大小为N ×M区域内行方向的 一阶差分， 表示核磁共振小波图像在邻域窗口大小为N ×M区域内列方向的一阶差分， 表示CT小波图像在邻域窗口大小为N ×M区域内行方向的一阶差分， 表示CT小波图像 在邻域窗口大小为N ×M区域内列方向的一阶差分， 表示叠加平均图像在邻域窗口大小 为N×M区域内行方向的一阶差分， 表示叠加平均图像在邻域窗口大小为N ×M区域内列 方向的一阶差分。 3.如权利要求2所述的一种核磁共振图像与CT图像 融合方法， 其特征在于， 所述步骤4： 根据信息梯度值构建融合系数， 包括： 步骤4.1： 根据信息梯度值计算加权系数； 其中， 所述加权系数计算公式为：权　利　要　求　书 1/3 页 2 CN 115471432 A 2其中， W1(i,j)表示第一加权系数， W2(i,j)表示第二加权系数， 表示核磁共振小波图 像总的信息梯度值， 表示CT图像总的信息梯度值， 表示叠加平均图像总的信息梯度 值； 步骤4.2： 根据所述加权系数 得到融合系数。 4.如权利要求3所述的一种核磁共振图像与CT图像融合方法， 其特征在于， 所述步骤 4.2： 根据所述加权系数 得到融合系数， 包括： 采用公式： 计算融合系数； 其中， ωB2(i,j)表示融合系数， ωA2(i,j)表示中间系数。 5.如权利要求4所述的一种核磁共振图像与CT图像 融合方法， 其特征在于， 所述步骤5： 计算所述核磁共 振小波图像和所述CT 小波图像的信息值， 包括： 步骤5.1： 根据所述核磁共振小波图像和所述CT小波图像的绝对值计算核磁共振小波 图像和所述CT 小波图像的能量 值； 步骤5.2： 根据所述能量值得到核磁共振小波图像和所述CT小波图像的信息值； 其中， 信息值计算公式为： 式中， ωA2(i,j)表示核磁共振小波图像的信息值， ωB2(i,j)表示CT小波图像的信息值， REA(i,j)表示核磁共 振小波图像的能量 值， REB(i,j)表示CT 小波图像的能量 值。 6.如权利要求5所述的一种核磁共振图像与CT图像 融合方法， 其特征在于， 所述步骤6： 根据所述信息值构建融合公式， 包括： 采用公式： 构建融合公式； 其中， 表示核磁共振小波图像在(i+m,j+n)点的像素值， 表示CT图像在(i+m,j+n)点的像素值。 7.如权利要求6所述的一种核磁共振图像与CT图像 融合方法， 其特征在于， 所述步骤7： 根据所述融合公式和融合系数构建图像融合模型， 包括：权　利　要　求　书 2/3 页 3 CN 115471432 A 3