[发明专利]一种FT域的医学图像融合方法

权利要求书

1.一种FT域的医学图像融合方法，其特征在于，包括以下步骤：

步骤1：针对所有待融合医学源图像进行FT变换，在经过L次尺度分解后，每幅源图像均分别得到1幅低频子带图像和8*L幅高频子带图像；

步骤2：计算低频子带图像中各像素点对应的ISML；获得初始融合决策映射图；采用GF模型分别对三种不同尺度下的初始融合决策映射图进行平滑处理，获得平滑处理后的融合决策映射图；针对平滑处理后的融合决策映射图进一步处理，获得最终融合决策映射图；完成低频子带图像的融合；

步骤3：对高频子带图像系数进行处理；计算高频子带图像中每个像素点的SF数值；完成高频子带图像的融合；

步骤4：将最终融合图像的高频子带图像和低频子带图像进行FT逆变换，得到最终融合图像F。

2.根据权利要求1所述的FT域的医学图像融合方法，其特征在于，所述步骤1包括：输入所有待融合医学源图像，并对这些医学源图像分别进行FT变换，尺度分解级数为L，(l,k)为l尺度下的方向分解级数，其中1≤l≤L，1≤k≤8；每幅待融合医学源图像经过FT变换后，分别得到1幅低频子带图像和8幅高频子带图像。

3.根据权利要求1所述的FT域的医学图像融合方法，其特征在于，所述步骤2中计算低频子带图像中各像素点对应的ISML包括：

计算低频子带图像中各像素点对应的ISML，其数学表达式为：

相应地，各像素点对应的ISML数值可以通过式(2)计算获得：

其中，N和T分别表示窗口半径和阈值；N为一个正奇数，T为0；式(2)可以改写为：

将ISML计算过程中的尺度确定为3级，并将3级对应的参数N分别设置为5,9和17；每个像素点对应不同N数值时的ISML数值可进一步改写为：

ISML_nN(i,j)＝ISML(I_n(i,j),N),n＝1,2,...Num,N＝5,9,17 (4)

其中，(i,j)表示图像中像素点的位置，Num为待融合医学源图像的数量。

4.根据权利要求1所述的FT域的医学图像融合方法，其特征在于，所述步骤2中获得初始融合决策映射图包括：

假设有两幅待融合医学源图像，则对应低频子带图像的初始融合决策映射图通过下列数学表达式获得：

选取参数N的数值分别为5,9和17，尺度参数s的取值范围为区间[1,3]内的整数；根据式(5)，若第一幅图像中像素点的ISML数值大于等于第二幅图像中相同位置像素点的ISML数值，则该像素点在初始融合决策映射图map_s中的对应元素数值为1，否则置为0。