[发明专利]一种实时的多模态医学图像融合方法

摘要

本发明公开了一种实时的多模态医学图像融合方法，该方法包括对已配准的若干幅医学源图像进行拉普拉斯金字塔分解得到每一幅源图像对应的低频和高频分量；基于稀疏表示的方法对所有源图像的低频分量进行融合处理，基于系数绝对值的方式对所有源图像的高频分量进行融合处理；将融合后的低频分量与高频分量进行拉普拉斯金字塔重建得到多模态医学融合图像。本发明方法很好地弥补了传统基于拉普拉斯金字塔变换和稀疏表示方法各自的不足，可以获得明显优于传统方法的结果。同时，本发明方法计算效率高，实时性好，在医学临床诊断等场合具有很大的应用价值。

主权项

一种实时的多模态医学图像融合方法，其特征在于，该方法包括：对已配准的若干幅医学源图像进行拉普拉斯金字塔分解得到每一幅源图像对应的低频和高频分量；基于稀疏表示的方法对所有源图像的低频分量进行融合处理，基于系数绝对值的方式对所有源图像的高频分量进行融合处理；将融合后的低频分量与高频分量进行拉普拉斯金字塔重建得到多模态医学融合图像；所述基于稀疏表示的方法对所有源图像的低频分量进行融合处理包括：设对源图像IA和IB的低频分量LA和LB进行融合，则首先采用滑窗技术从低频分量LA和LB中提取重叠的图像块；具体的：对于低频分量LA和LB，均从左上角到右下角移动一个大小为的滑窗，其中，为字典D中原子的尺寸；取出相应位置的图像块保存，再移动滑窗，每次移动的步长为s个像素，直至提取完毕；若源图像IA和IB的大小为H×W，则每一源图像的低频分量中提取的图像块个数T表示为：其中，符号表示向上取整；设和分别表示从低频分量LA和LB中取出的图像块；将低频分量LA和LB中第i个位置的图像块和分别展开成列向量和再归一化使其所有元素的均值为0，处理过程表示为：v^xi=vxi-v‾xi·1;]]>其中，1表示一个所有元素均为1的n×1向量，和分别是和的均值；利用字典D分别对和进行稀疏分解，得到稀疏表示系数和处理过程表示为：αAi=argminα||α||0s.t.||v^Ai-Dα||2<ϵ;]]>αBi=argminα||α||0s.t.||v^Bi-Dα||2<ϵ;]]>其中，ε表示稀疏重构的允许误差，||·||0为向量的零范数；采用比较L1范数的方式对和进行融合，得到融合后的稀疏系数：αFi=αAiif||αAi||1>||αBi||1αBiotherwise;]]>从而获得低频分量LA和LB中第i个位置图像块列向量的融合结果：vFi=DαFi+v‾Fi·1;]]>其中，表示融合向量的均值；对低频分量LA和LB中所有位置的图像块的列向量进行融合，获得融合向量建立一个大小为H×W的空图像LF，对每一个融合结果将其变形为大小的图像块，再放入空图像LF的相应位置，从而完成所有源图像低频分量的融合。