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

权利要求书

1.一种多模态医学图像的融合方法，其特征在于，包括：

针对多幅待融合的医学源图像分别进行非下采样剪切波变换，每幅所述医学源图像分别得到对应的低频子带图像和多幅高频子带图像；

采用引导核加权引导滤波对所述低频子带图像进行融合处理，得到低频融合子图像；

采用侧窗滤波对多幅所述高频子带图像进行融合处理，得到高频融合子图像；

利用非下采样剪切波变换的逆变换将所述低频融合子图像和所述高频融合子图像进行融合，得到最终融合图像。

2.根据权利要求1所述多模态医学图像的融合方法，其特征在于，所述针对多幅待融合的医学源图像分别进行非下采样剪切波变换，每幅所述医学源图像分别得到对应的低频子带图像和多幅高频子带图像的步骤包括：针对多幅待融合的医学源图像分别进行非下采样剪切波变换，在经过s级尺度分解和ds方向分解后，每幅所述医学源图像分别得到对应的1幅低频子带图像和多幅高频子带图像，其中，S≥1，s＝1,2,…,S；ds表示第s级尺度下的方向分解数。

3.根据权利要求1所述多模态医学图像的融合方法，其特征在于，所述采用引导核加权引导滤波对所述低频子带图像进行融合处理，得到低频融合子图像的步骤包括：

分别计算每幅所述低频子带图像中每个像素点的拉普拉斯能量和，所述每个像素点的拉普拉斯能量和的计算公式包括：

其中，f(x,y)表示位于医学源图像对应的低频子带图像中空间位置(x,y)处的像素点；N表示以像素点f(x,y)为中心像素点的局部区域的半径；

step表示步长；

分别从每幅所述低频子带图像中选取多个以不同像素点为中心半径的局部区域，得到每幅所述低频子带图像对应的所述拉普拉斯能量和；并通过对比多幅所述低频子带图像中相同局部区域半径的所述拉普拉斯能量和，得到多幅低频初始融合决策映射图；

采用引导核加权分别对每幅所述初始低频融合决策映射图进行引导滤波处理，得到多幅低频融合决策引导映射图，并从多幅所述低频融合决策引导映射图中选取相同像素点位置处的最大值作为元素值，从而确定低频最终融合决策映射图；

依据所述低频最终融合决策映射图，对多幅所述低频子带图像进行融合，得到所述低频融合子图像。

4.根据权利要求3所述多模态医学图像的融合方法，其特征在于，分别从每幅所述低频子带图像中选取多个以不同像素点为中心半径的局部区域，得到每幅所述低频子带图像对应的所述拉普拉斯能量和；并通过对比多幅所述低频子带图像中相同局部区域半径的所述拉普拉斯能量和，得到多幅低频初始融合决策映射图的步骤包括：

分别从每幅所述低频子带图像中选取多个以不同像素点为中心半径N的局部区域，其中，N∈(N₁,N₂,N₃,...,N_n)，得到每幅所述低频子带图像对应的所述拉普拉斯能量和SML的计算公式转变成：

其中，A和B分别表示源图像A和源图像B；

所述低频初始融合决策映射图的公式包括：

其中，f_A(x,y)表示位于医学源图像A对应的低频子带图像中空间位置(x,y)处的像素点；f_B(x,y)表示位于医学源图像B对应的低频子带图像中空间位置(x,y)处的像素点。