本发明涉及一种图像处理方法，尤其是一种基于对比度分析的图像融合方法。
背景技术：
图像增强是图像处理必不可少的环节。无论是频域的处理方法还是空域的方法，主要作用都是有目的地强调图像的整体或局部特性，将原来不清晰的图像变得清晰或强调某些感兴趣的特征，扩大图像中不同物体特征之间的差别，抑制不感兴趣的特征，使之改善图像质量、丰富信息量。
但是，通过增强的图像都有一些缺陷，如部分细节信息的丢失，亮度提升过多造成过曝现象，对比度拉伸造成部分区域过暗等等。
技术实现要素：
本发明的目的是提供一种使得图像获得增强处理效果的同时，在细节上也保留了更好的对比度的方法。
本发明的技术方案是：
一种基于对比度分析的图像融合方法，它包括以下步骤：
步骤一：分别构建输入图像及其增强图像的形态学梯度图；
步骤二：将输入图像及其增强图像的形态学梯度图相减，得到输入图像增强前后的对比度差图；
步骤三：根据对比度差图构建融合系数矩阵；
步骤四：根据融合系数矩阵完成输入图像及其增强图像的融合，对融合后的处理图像像素进行范围控制，使像素在图像合理范围内并输出。
进一步地，步骤一具体为：
步骤101：对输入图像及其增强图像进行去噪处理；
步骤102：分别计算输入图像及其增强图像中一定邻域范围内的对比度，得到形态学梯度图一和形态学梯度图二，计算公式为：
ci(x，y)＝max(ni(x，y))-min(ni(x，y))公式1
其中，x，y是输入图像及其增强图像的坐标；i＝1，2分别代表输入图像及其增强图像；n代表邻域范围；c代表形态学梯度图。
进一步地，领域范围选择3×3邻域、5×5邻域或者7×7邻域。
进一步地，步骤二具体为：
步骤201：将步骤一中得到的输入图像及其增强图像的形态学梯度图相减，得到对比度差图d，计算公式为：
d(x,y)＝c2(x,y)–c1(x,y)公式2
步骤202：对对比度差图d进行归一化处理，保证数值的有效性。
进一步地，步骤三具体为：
步骤301：计算输入图像及其增强图像的融合系数矩阵w，计算公式为：
其中：a,b为梯度的陡度调节系数，d表示对比度差图的矩阵；
步骤302：对融合系数矩阵w进行归一化处理，保证数值的有效性。
进一步地，步骤四具体为：
步骤401：根据融合系数矩阵完成输入图像及其增强图像的融合，计算公式如下：
r＝w×i2+(1-w)×i1公式4
其中，r代表融合的结果，i1，i2分别代表输入图像及其增强图像；
步骤402：对步骤401中输出的结果r进行归一化处理，去除溢出数值；
步骤403：输出处理后的图像。
进一步地，步骤301中，a＝5，b＝-min(d)/(max(d)-min(d))；max(d)和min(d)分别表示矩阵中的最大值和最小值。
本发明的有益效果：
本发明的单图融合是针对输入图像与其增强图像的融合，图像信息本质是一样的，并不是不同的图像融合。通过本方法融合输出的图像灰度分布更加均匀，细节完整的保留下来。本方法不仅对灰度图像有效，对彩色图像也同样有效，且对色彩图像在相关色彩通道上的还原更加保真。
附图说明
通过结合附图对本发明示例性实施方式进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施方式中，相同的参考标号通常代表相同部件。
图1示出了本发明的原理框图。
具体实施方式
下面将参照附图更详细地描述本发明的优选实施方式。虽然附图中显示了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。
一种基于对比度分析的图像融合方法，它包括以下步骤：
步骤一：分别构建输入图像及其增强图像的形态学梯度图；(所述的增强图像可以为：针对目标、边缘、分辨率、对比度等的增强图像)
步骤二：将输入图像及其增强图像的形态学梯度图相减，得到输入图像增强前后的对比度差图；
步骤三：根据对比度差图构建融合系数矩阵；
步骤四：根据融合系数矩阵完成输入图像及其增强图像的融合，对融合后的处理图像像素进行范围控制，使像素在图像合理范围内并输出。
本发明的步骤一具体中：
步骤101：对输入图像及其增强图像进行去噪处理；
步骤102：分别计算输入图像及其增强图像中一定邻域范围内的对比度，得到形态学梯度图一和形态学梯度图二，计算公式为：
ci(x，y)＝max(ni(x，y))-min(ni(x，y))公式1
其中，x，y是输入图像及其增强图像的坐标；i＝1，2分别代表输入图像及其增强图像；n代表邻域范围(优选领域范围选择3×3邻域)；c代表形态学梯度图；max表示求邻域的最大值；min表示求邻域的最小值。
本发明的步骤二具体为：
步骤201：将步骤一中得到的输入图像及其增强图像的形态学梯度图相减，得到对比度差图d，计算公式为：
d(x,y)＝c2(x,y)–c1(x,y)公式2
步骤202：对对比度差图d进行归一化处理，保证数值的有效性。
本发明的步骤三具体为：
步骤301：计算输入图像及其增强图像的融合系数矩阵w，计算公式为：
其中：a,b为梯度的陡度调节系数，d表示对比度差图的矩阵；a＝5，b＝-min(d)/(max(d)-min(d))；max(d)和min(d)分别表示矩阵中的最大值和最小值；
步骤302：对融合系数矩阵w进行归一化处理，保证数值的有效性。
本发明的步骤四具体为：
步骤401：根据融合系数矩阵完成输入图像及其增强图像的融合，计算公式如下：
r＝w×i2+(1-w)×i1公式4
其中，r代表融合的结果，i1，i2分别代表输入图像及其增强图像；
步骤402：对步骤401中输出的结果r进行归一化处理，去除溢出数值；输出处理后的图像。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。