基于纹理显著性赋权的遥感影像融合方法及系统与流程

1.本发明涉及遥感影像处理领域，具体地，涉及基于纹理显著性赋权的遥感影像融合方法及系统。


背景技术：

2.随着航空和航天遥感科学技术的发展，携带不同传感器的遥感卫星陆续发射，越来越多的传感器投入使用，多源航空航天遥感影像数据呈现爆发式增长。各类遥感影像数据在资源监测、城市建设规划、自然灾害预警等相关领域的应用越来越多。不同传感器具有其独特的用途，所以每种传感器在获取遥感卫星影像方面有其独特的特征。在通常情况下，不同特点的卫星传感器在经过的相同区域获得大量的遥感数据，而这些影像数据的空间、时间以及光谱分辨率不同，大多数情况下仅代表了设计者所需要获取的遥感数据信息。所以这些传感器也都有其各自的局限性，获取的影像存在较多冗余信息，而如何综合多种传感器获取的影像信息特点、弥补单一传感器在获取影像方面的不足，成为了遥感数据预处理研究的重点内容。
3.在遥感影像融合算法中，基于hsi变换的算法一般应用于异源遥感影像的处理，在对相同传感器获取的不同波段之间的融合也具有较好的效果。通过对同一区域高分全色影像以及低分多光谱影像进行融合，可以获得结合二者优点的影像数据。但基于his变换的影像融合也有着一定的缺点，影像融合结果的提升主要体现在保留更多高空间分辨率影像的纹理信息以及地空间分辨率影像的光谱信息，但hsi影像融合采用的“i分量替换”、“平均加权融合”等融合方法，均从影像全局一致性出发，采用统一融合权重，所以在提高空间分辨率的同时会导致一定程度的光谱畸变，若要保持光谱畸变最小，则需要损失一定的空间分辨率。
4.针对上述问题，目前需要一套能够同步解决以下技术问题的技术方案：（1）如何有效提高遥感影像的空间分辨率；（2）如何在影像融合时确定合适的加权权重；以满足影像融合在提高遥感影像空间分辨率的同时，尽可能的减少光谱畸变。


技术实现要素：

5.针对现有技术中的缺陷，本发明的目的是提供一种基于纹理显著性赋权的遥感影像融合方法及系统。
6.根据本发明提供的一种基于纹理显著性赋权的遥感影像融合方法，包括：步骤s1：对高分辨率影像进行纹理显著性计算得到纹理显著性图像；步骤s2：将纹理显著性图像归一化至0-1范围并进行分段权重运算得到权重矩阵；步骤s3：将高分辨率影像和低分辨率影像转换至hsi颜色空间并提取强度信息i波段；步骤s4：高分辨率影像i波段与低分辨率影像i波段进行灰度直方图匹配；
步骤s5：直方图匹配后的高分辨率影像i波段与低分辨率影像i波段基于权重矩阵进行加权运算得到新强度波段；步骤s6：使用新强度波段替换低分辨率影像i波段，将hsi空间图像逆转换至rgb空间得到融合影像，完成基于纹理显著性赋权的遥感影像融合。
7.优选地，所述步骤s1采用：其中，表示影像坐标处的显著性值，表示影像坐标处的像素值。
8.优选地，所述步骤s2采用：步骤s2.1：根据影像显著性矩阵最大、最小值，将显著性信息矩阵元素归一化至0-1范围；其中，表示影像坐标处的显著性值，表示影像显著性最小值，表示影像显著性最大值；表示归一化处理后的影像坐标处的显著性值；步骤s2.2：对归一化处理后的纹理显著性图像采用分段方式计算融合权重矩阵；其中，表示融合权重矩阵在影像坐标处的值；表示归一化处理后的影像坐标处的显著性值，表示模型系数，表示权重计算的分段阈值。
9.优选地，所述步骤s3采用：将高分辨率影像和低分辨率影像转换至hsi颜色空间；若，则；其中，表示hsi变换后的色度波段；表示影像红色波段；表示影像绿色波段；表示影像蓝色波段；表示hsi变换后的亮度波段；表示hsi变换后的饱和度波段。
10.优选地，所述步骤s4采用：步骤s4.1：由高分辨率影像i波段及低分辨率影像i波段得到灰度级的高分辨率影像i波段的连续概率密度函数以及低分辨率影像i波段的连续概率密度函数，并计算具有如下特征的随机变量；
其中，表示积分假变量；表示影像灰度级总数；表示连续概率密度函数的积分上限；表示高分辨率影像i波段连续概率密度函数在灰度值w处的概率密度值；步骤s4.2：基于预设的概率密度函数求得变换函数；其中，表示积分假变量；表示连续概率密度函数的积分上限，表示低分辨率影像i波段的连续概率密度函数在灰度值处的概率密度值；步骤s4.3：求得反变换函数。
11.优选地，所述步骤s5采用：其中，为加权计算得到的新强度波段波段在影像坐标处的值；为高分辨率影像波段在影像坐标处的值；为低分辨率影像波段在影像坐标处的值；为影像坐标高分辨率影像波段的计算权重。
12.优选地，所述步骤s6采用：hsi颜色空间变换公式为：若且且若若若
。
13.根据本发明提供的一种基于纹理显著性赋权的遥感影像融合系统，包括：模块m1：对高分辨率影像进行纹理显著性计算得到纹理显著性图像；模块m2：将纹理显著性图像归一化至0-1范围并进行分段权重运算得到权重矩阵；模块m3：将高分辨率影像和低分辨率影像转换至hsi颜色空间并提取强度信息i波段；模块m4：高分辨率影像i波段与低分辨率影像i波段进行灰度直方图匹配；模块m5：直方图匹配后的高分辨率影像i波段与低分辨率影像i波段基于权重矩阵进行加权运算得到新强度波段；模块m6：使用新强度波段替换低分辨率影像i波段，将hsi空间图像逆转换至rgb空间得到融合影像，完成基于纹理显著性赋权的遥感影像融合。
14.优选地，所述模块m1采用：其中，表示影像坐标处的显著性值，表示影像坐标处的像素值；所述模块m2采用：模块m2.1：根据影像显著性矩阵最大、最小值，将显著性信息矩阵元素归一化至0-1范围；其中，表示影像坐标处的显著性值，表示影像显著性最小值，表示影像显著性最大值；表示归一化处理后的影像坐标处的显著性值；模块m2.2：对归一化处理后的纹理显著性图像采用分段方式计算融合权重矩阵；其中，表示融合权重矩阵在影像坐标处的值；表示归一化处理后的影像坐标处的显著性值，表示模型系数，表示权重计算的分段阈值；所述模块m3采用：将高分辨率影像和低分辨率影像转换至hsi颜色空间；
若，则；其中，表示hsi变换后的色度波段；表示影像红色波段；表示影像绿色波段；表示影像蓝色波段；表示hsi变换后的亮度波段；表示hsi变换后的饱和度波段。
15.优选地，所述模块m4采用：模块m4.1：由高分辨率影像i波段及低分辨率影像i波段得到灰度级的高分辨率影像i波段的连续概率密度函数以及低分辨率影像i波段的连续概率密度函数，并计算具有如下特征的随机变量；其中，表示积分假变量；表示影像灰度级总数；表示连续概率密度函数的积分上限；表示高分辨率影像i波段连续概率密度函数在灰度值w处的概率密度值；模块m4.2：基于预设的概率密度函数求得变换函数；其中，表示积分假变量；表示连续概率密度函数的积分上限，表示低分辨率影像i波段的连续概率密度函数在灰度值处的概率密度值；模块m4.3：求得反变换函数；所述模块m5采用：其中，为加权计算得到的新强度波段波段在影像坐标处的值；为高分辨率影像波段在影像坐标处的值；为低分辨率影像波段在影像坐标处的值；为影像坐标高分辨率影像波段的计算权重；所述模块m6采用：hsi颜色空间变换公式为：若且
若若若若。
16.与现有技术相比，本发明具有如下的有益效果：本发明通过提取高分辨卫星影像纹理显著性特征，根据影像不同位置纹理显著性进行融合权重值的设置，使得在融合过程中，高纹理显著性区域保留较多的空间纹理信息，低纹理区域保留较多的光谱信息，从而在影像融合过程中降低了传统遥感影像hsi变换融合中非纹理显著区域的光谱畸变量。
17.本发明可实现高分辨率全色波段影像与低分辨率多光谱影像的融合，降低了传统遥感影像hsi变换融合中非纹理显著区域的光谱畸变量。
附图说明
18.通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：图1是本发明的一种基于纹理显著性赋权的遥感影像融合方法的流程示意图；图2a-图2c是本发明中利用该技术进行影像融合与传统融合方法结果对比图。
具体实施方式
19.下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术
人员来说，在不脱离本发明构思的前提下，还可以做出若干变化和改进。这些都属于本发明的保护范围。
20.实施例1根据本发明提供的一种基于纹理显著性赋权的遥感影像融合方法，包括：步骤s1：对高分辨率影像进行纹理显著性计算得到纹理显著性图像；步骤s2：将纹理显著性图像归一化至0-1范围并进行分段权重运算得到权重矩阵；步骤s3：将高分辨率影像和低分辨率影像转换至hsi颜色空间并提取强度信息i波段；若高分辨率影像为单波段影像，直接取该影像为强度信息i波段；步骤s4：高分辨率影像i波段与低分辨率影像i波段进行灰度直方图匹配；步骤s5：直方图匹配后的高分辨率影像i波段与低分辨率影像i波段基于权重矩阵进行加权运算得到新强度波段；步骤s6：使用新强度波段替换低分辨率影像i波段，将hsi空间图像逆转换至rgb空间得到融合影像，完成基于纹理显著性赋权的遥感影像融合。
21.具体地，所述步骤s1采用：通过对高分辨率影像进行显著性计算，得到显著性图像，主要体现在纹理明显区域具有较大数值，而纹理不明显区域即地物单一区域具有较小数值；更为具体地，其中，表示影像坐标处的显著性值，表示影像坐标处的像素值。
22.具体地，基于步骤s1得到的纹理显著性图像，进行分段权重运算，所述步骤s2中分段权重运算的具体描述为，影像融合主要目的为获取同时包含高分辨率影像主要表现为纹理特征的较高空间分辨率，以及低分辨率影像的较高光谱分辨率特征。在不同纹理显著性范围，需要保留的空间分辨率特征及光谱分辨率特征不同。因此，在加权融合时，在纹理显著性强的区域给予高分辨率影像更高权重，反之给予高光谱影像更高权重。
23.具体地，所述步骤s2中分段权重运算的具体描述为，基于步骤s1得到的纹理显著性图像，将其影像值范围归一化至0-1，并根据显著性值大小进行分段权重运算。在加权融合时，在纹理显著性强的区域给予高分辨率影像更高权重，反之给予高光谱影像更高权重。
24.所述步骤s2采用：步骤s2.1：根据影像显著性矩阵最大、最小值，将显著性信息矩阵元素归一化至0-1范围；其中，表示影像坐标处的显著性值，表示影像显著性最小值，表示影像显著性最大值；表示归一化处理后的影像坐标处的显著性值；步骤s2.2：对归一化处理后的纹理显著性图像采用分段方式计算融合权重矩阵；
其中，表示融合权重矩阵在影像坐标处的值；表示归一化处理后的影像坐标处的显著性值，表示模型系数，表示权重计算的分段阈值。
25.具体地，所述步骤s3中将影像转换至hsi空间中进行强度信息波段提取，该过程可以在一定程度上将强度信息集中在i波段，强度波段能够最大程度的表征代表空间分辨率的纹理信息，因此在所述步骤s6影像逆变换过程进行波段替换可以在一定程度上尽可能将空间分辨率信息引入多光谱影像，而减少光谱畸变。
26.将高分辨率影像和低分辨率影像转换至hsi颜色空间；若，则；其中，表示hsi变换后的色度波段；表示影像红色波段；表示影像绿色波段；表示影像蓝色波段；表示hsi变换后的亮度波段；表示hsi变换后的饱和度波段。
27.具体地，所述步骤s4采用：步骤s4.1：由高分辨率影像i波段及低分辨率影像i波段得到灰度级的高分辨率影像i波段的连续概率密度函数以及低分辨率影像i波段的连续概率密度函数，并计算具有如下特征的随机变量；其中，表示积分假变量；表示影像灰度级总数；表示连续概率密度函数的积分上限；表示高分辨率影像i波段连续概率密度函数在灰度值w处的概率密度值；步骤s4.2：基于预设的概率密度函数求得变换函数；其中，表示积分假变量；表示连续概率密度函数的积分上限，表示低分辨率影像i波段的连续概率密度函数在灰度值处的概率密度值；步骤s4.3：求得反变换函数。
28.具体地，所述步骤s5中通过基于纹理显著性分段得到的权重矩阵，对高分辨率影像与低分辨率影像i波段进行加权运算得到新的强度信息波段，该波段在纹理特征明显区域保留了高分辨率影像的强度信息，具有较强的纹理信息特征，在纹理特征不明显的区域保留了低分辨率影像的强度信息，在hsi逆变换后得到的融合影像中具有较低的光谱畸变，从而实现了更好的基于纹理特征显著性赋权的遥感影像融合。
29.所述步骤s5采用：
其中，为加权计算得到的新强度波段波段在影像坐标处的值；为高分辨率影像波段在影像坐标处的值；为低分辨率影像波段在影像坐标处的值；为影像坐标高分辨率影像波段的计算权重。
30.优选地，所述步骤s6采用：hsi颜色空间变换公式为：若且且若若若若。
31.根据本发明提供的一种基于纹理显著性赋权的遥感影像融合系统，包括：模块m1：对高分辨率影像进行纹理显著性计算得到纹理显著性图像；模块m2：将纹理显著性图像归一化至0-1范围并进行分段权重运算得到权重矩阵；模块m3：将高分辨率影像和低分辨率影像转换至hsi颜色空间并提取强度信息i波段；若高分辨率影像为单波段影像，直接取该影像为强度信息i波段；模块m4：高分辨率影像i波段与低分辨率影像i波段进行灰度直方图匹配；模块m5：直方图匹配后的高分辨率影像i波段与低分辨率影像i波段基于权重矩阵
进行加权运算得到新强度波段；模块m6：使用新强度波段替换低分辨率影像i波段，将hsi空间图像逆转换至rgb空间得到融合影像，完成基于纹理显著性赋权的遥感影像融合。
32.具体地，所述模块m1采用：通过对高分辨率影像进行显著性计算，得到显著性图像，主要体现在纹理明显区域具有较大数值，而纹理不明显区域即地物单一区域具有较小数值；更为具体地，其中，表示影像坐标处的显著性值，表示影像坐标处的像素值。
33.具体地，基于模块m1得到的纹理显著性图像，进行分段权重运算，所述模块m2中分段权重运算的具体描述为，影像融合主要目的为获取同时包含高分辨率影像主要表现为纹理特征的较高空间分辨率，以及低分辨率影像的较高光谱分辨率特征。在不同纹理显著性范围，需要保留的空间分辨率特征及光谱分辨率特征不同。因此，在加权融合时，在纹理显著性强的区域给予高分辨率影像更高权重，反之给予高光谱影像更高权重。
34.具体地，所述模块m2中分段权重运算的具体描述为，基于模块m1得到的纹理显著性图像，将其影像值范围归一化至0-1，并根据显著性值大小进行分段权重运算。在加权融合时，在纹理显著性强的区域给予高分辨率影像更高权重，反之给予高光谱影像更高权重。
35.所述模块m2采用：模块m2.1：根据影像显著性矩阵最大、最小值，将显著性信息矩阵元素归一化至0-1范围；其中，表示影像坐标处的显著性值，表示影像显著性最小值，表示影像显著性最大值；表示归一化处理后的影像坐标处的显著性值；模块m2.2：对归一化处理后的纹理显著性图像采用分段方式计算融合权重矩阵；其中，表示融合权重矩阵在影像坐标处的值；表示归一化处理后的影像坐标处的显著性值，表示模型系数，表示权重计算的分段阈值。
36.具体地，所述模块m3中将影像转换至hsi空间中进行强度信息波段提取，该过程可以在一定程度上将强度信息集中在i波段，强度波段能够最大程度的表征代表空间分辨率的纹理信息，因此在所述模块m6影像逆变换过程进行波段替换可以在一定程度上尽可能将空间分辨率信息引入多光谱影像，而减少光谱畸变。
37.将高分辨率影像和低分辨率影像转换至hsi颜色空间；
若，则；其中，表示hsi变换后的色度波段；表示影像红色波段；表示影像绿色波段；表示影像蓝色波段；表示hsi变换后的亮度波段；表示hsi变换后的饱和度波段。
38.具体地，所述模块m4采用：模块m4.1：由高分辨率影像i波段及低分辨率影像i波段得到灰度级的高分辨率影像i波段的连续概率密度函数以及低分辨率影像i波段的连续概率密度函数，并计算具有如下特征的随机变量；其中，表示积分假变量；表示影像灰度级总数；表示连续概率密度函数的积分上限；表示高分辨率影像i波段连续概率密度函数在灰度值w处的概率密度值；模块m4.2：基于预设的概率密度函数求得变换函数；其中，表示积分假变量；表示连续概率密度函数的积分上限，表示低分辨率影像i波段的连续概率密度函数在灰度值处的概率密度值；模块m4.3：求得反变换函数。
39.具体地，所述模块m5中通过基于纹理显著性分段得到的权重矩阵，对高分辨率影像与低分辨率影像i波段进行加权运算得到新的强度信息波段，该波段在纹理特征明显区域保留了高分辨率影像的强度信息，具有较强的纹理信息特征，在纹理特征不明显的区域保留了低分辨率影像的强度信息，在hsi逆变换后得到的融合影像中具有较低的光谱畸变，从而实现了更好的基于纹理特征显著性赋权的遥感影像融合。
40.所述模块m5采用：其中，为加权计算得到的新强度波段波段在影像坐标处的值；为高分辨率影像波段在影像坐标处的值；为低分辨率影像波段在影像坐标处的值；为影像坐标高分辨率影像波段的计算权重。
41.优选地，所述模块m6采用：hsi颜色空间变换公式为：
若且且若若若若。
42.实施例2实施例2是实施例1的优选例为了使本领域的技术人员更好地理解本技术方案，下面将结合本技术实施例中的附图对本技术实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例仅仅是本技术一部分的实施例，而不是全部的实施例；基于本技术中的实施例，本领域普通技术人员在没有做出创造性劳动成果前提下所获得的所有其他实施例，都应当属于本技术保护的范围；需要说明的是，在不冲突的情况下，本技术中的实施例及实施例中的特征可以相互组合，下面将参考附图1-2，并结合实施例来详细说明本技术；其中，图1为基于纹理显著性赋权的遥感影像融合方法流程图；图2为利用该技术进行影像融合与传统融合方法结果对比图，其中，图2（a）为原多光谱影像；图2（b）为传统his融合结果；图2（c）为本发明融合结果；基于此，本发明展现的影像更为清晰，同时降低了传统遥感影像hsi变换融合中非纹理显著区域的光谱畸变量。
43.本发明提出一种基于纹理显著性赋权的遥感影像融合方法，解决当前传统影像融
合方法提高空间分辨率的同时会导致一定程度的光谱畸变的问题，具体涉及影像颜色空间变换、影像显著性提取、影像融合分段权重运算等相关技术，如图1所示，包括以下步骤：步骤1：利用高分辨率卫星影像纹理显著性特征，提出一种通过提取纹理显著性图并进行分段计算融合权重的方法，用于遥感影像融合方法的权重计算。
44.根据高分辨率遥感影像纹理特征，计算各向纹理梯度，得到影像纹理显著性矩阵，表征影像各像素、区域的纹理显著性信息。
45.影像纹理显著性计算公式为：其中，表示影像坐标处的显著性值，影像坐标处的像素值。
46.步骤2：对显著性图像归一化至0-1范围并进行分段权重运算；a)根据影像显著性矩阵最大、最小值，将显著性信息矩阵元素归一化至0-1范围，保证影像融合后像素值不会溢出。
47.b)为保证低纹理信息区域更好的保留源影像的光谱信息，采用分段方式计算融合权重矩阵，使纹理信息丰富区域保留更多纹理特征，纹理信息欠缺区域通常也是光谱特征单一且重要区域更好的与源影像光谱信息一致。
48.步骤3：高、低分辨率影像转换至hsi颜色空间并提取强度信息i波段，若高分辨率影像为单波段影像，直接取该影像为强度信息i波段；hsi颜色空间变换公式为：若，则其中，h表示hsi变换后的色度波段，s表示hsi变换后的饱和度波段，i表示hsi变换后的亮度波段，r表示影像红色波段，g表示影像绿色波段，b表示影像蓝色波段。
49.步骤4：高分辨率影像i波段与低分辨率影像i波段进行灰度直方图匹配；a)由高分辨率影像i波段及低分辨率影像i波段得到灰度级的连续概率密度函数、，并计算具有如下特征的随机变量。
50.其中，为积分假变量，表示影像灰度级总数，表示连续概率密度函数的积分上限，表示高分辨率影像i波段连续概率密度函数在灰度值w处的概率密度值。
51.b) 由下式指定的概率密度函数求得变换函数。
52.其中，为积分假变量。表示连续概率密度函数的积分上限，表示低分辨率影像i波段连续概率密度函数在灰度值处的概率密度值。
53.c)求得反变换函数，由于是由得到的，所以该处理是到的映射，而后者正是我们期望的值。
54.步骤5：高分辨率影像与低分辨率影像i波段进行加权运算得到新的波段i_new；加权计算公式为：其中，为加权计算得到的波段在影像坐标处的值；为高分辨率影像波段在影像坐标处的值；为低分辨率影像波段在影像坐标处的值；为影像坐标高分辨率影像波段的计算权重。
55.步骤6：使用i_new替换i波段，将hsi空间图像逆转换至rgb空间，得到融合影像，完成基于纹理显著性赋权的遥感影像融合；hsi颜色空间变换公式为：若且且若若若
本领域技术人员知道，除了以纯计算机可读程序代码方式实现本发明提供的系统、装置及其各个模块以外，完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统、装置及其各个模块以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同程序。所以，本发明提供的系统、装置及其各个模块可以被认为是一种硬件部件，而对其内包括的用于实现各种程序的模块也可以视为硬件部件内的结构；也可以将用于实现各种功能的模块视为既可以是实现方法的软件程序又可以是硬件部件内的结构。
56.以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本技术的实施例和实施例中的特征可以任意相互组合。

技术特征：
1.一种基于纹理显著性赋权的遥感影像融合方法，其特征在于，包括：步骤s1：对高分辨率影像进行纹理显著性计算得到纹理显著性图像；步骤s2：将纹理显著性图像归一化至0-1范围并进行分段权重运算得到权重矩阵；步骤s3：将高分辨率影像和低分辨率影像转换至hsi颜色空间并提取强度信息i波段；步骤s4：高分辨率影像i波段与低分辨率影像i波段进行灰度直方图匹配；步骤s5：直方图匹配后的高分辨率影像i波段与低分辨率影像i波段基于权重矩阵进行加权运算得到新强度波段；步骤s6：使用新强度波段替换低分辨率影像i波段，将hsi空间图像逆转换至rgb空间得到融合影像，完成基于纹理显著性赋权的遥感影像融合。2.根据权利要求1所述的基于纹理显著性赋权的遥感影像融合方法，其特征在于，所述步骤s1采用：其中，表示影像坐标处的显著性值，表示影像坐标处的像素值。3.根据权利要求1所述的基于纹理显著性赋权的遥感影像融合方法，其特征在于，所述步骤s2采用：步骤s2.1：根据影像显著性矩阵最大、最小值，将显著性信息矩阵元素归一化至0-1范围；其中，表示影像坐标处的显著性值，表示影像显著性最小值，表示影像显著性最大值；表示归一化处理后的影像坐标处的显著性值；步骤s2.2：对归一化处理后的纹理显著性图像采用分段方式计算融合权重矩阵；其中，表示融合权重矩阵在影像坐标处的值；表示归一化处理后的影像坐标处的显著性值，表示模型系数，表示权重计算的分段阈值。4.根据权利要求1所述的基于纹理显著性赋权的遥感影像融合方法，其特征在于，所述步骤s3采用：将高分辨率影像和低分辨率影像转换至hsi颜色空间；若，则；其中，表示hsi变换后的色度波段；表示影像红色波段；表示影像绿色波段；表
示影像蓝色波段；表示hsi变换后的亮度波段；表示hsi变换后的饱和度波段。5.根据权利要求1所述的基于纹理显著性赋权的遥感影像融合方法，其特征在于，所述步骤s4采用：步骤s4.1：由高分辨率影像i波段及低分辨率影像i波段得到灰度级的高分辨率影像i波段的连续概率密度函数以及低分辨率影像i波段的连续概率密度函数，并计算具有如下特征的随机变量；其中，表示积分假变量；表示影像灰度级总数；表示连续概率密度函数的积分上限；表示高分辨率影像i波段连续概率密度函数在灰度值w处的概率密度值；步骤s4.2：基于预设的概率密度函数求得变换函数；其中，表示积分假变量；表示连续概率密度函数的积分上限，表示低分辨率影像i波段的连续概率密度函数在灰度值处的概率密度值；步骤s4.3：求得反变换函数。6.根据权利要求1所述的基于纹理显著性赋权的遥感影像融合方法，其特征在于，所述步骤s5采用：其中，为加权计算得到的新强度波段波段在影像坐标处的值；为高分辨率影像波段在影像坐标处的值；为低分辨率影像波段在影像坐标处的值；为影像坐标高分辨率影像波段的计算权重。7.根据权利要求1所述的基于纹理显著性赋权的遥感影像融合方法，其特征在于，所述步骤s6采用：hsi颜色空间变换公式为：若且且若
若若。8.一种基于纹理显著性赋权的遥感影像融合系统，其特征在于，包括：模块m1：对高分辨率影像进行纹理显著性计算得到纹理显著性图像；模块m2：将纹理显著性图像归一化至0-1范围并进行分段权重运算得到权重矩阵；模块m3：将高分辨率影像和低分辨率影像转换至hsi颜色空间并提取强度信息i波段；模块m4：高分辨率影像i波段与低分辨率影像i波段进行灰度直方图匹配；模块m5：直方图匹配后的高分辨率影像i波段与低分辨率影像i波段基于权重矩阵进行加权运算得到新强度波段；模块m6：使用新强度波段替换低分辨率影像i波段，将hsi空间图像逆转换至rgb空间得到融合影像，完成基于纹理显著性赋权的遥感影像融合。9.根据权利要求8所述的基于纹理显著性赋权的遥感影像融合系统，其特征在于，所述模块m1采用：其中，表示影像坐标处的显著性值，表示影像坐标处的像素值；所述模块m2采用：模块m2.1：根据影像显著性矩阵最大、最小值，将显著性信息矩阵元素归一化至0-1范围；其中，表示影像坐标处的显著性值，表示影像显著性最小值，表示影像显著性最大值；表示归一化处理后的影像坐标处的显著性值；模块m2.2：对归一化处理后的纹理显著性图像采用分段方式计算融合权重矩阵；
其中，表示融合权重矩阵在影像坐标处的值；表示归一化处理后的影像坐标处的显著性值，表示模型系数，表示权重计算的分段阈值；所述模块m3采用：将高分辨率影像和低分辨率影像转换至hsi颜色空间；若，则；其中，表示hsi变换后的色度波段；表示影像红色波段；表示影像绿色波段；表示影像蓝色波段；表示hsi变换后的亮度波段；表示hsi变换后的饱和度波段。10.根据权利要求8所述的基于纹理显著性赋权的遥感影像融合系统，其特征在于，所述模块m4采用：模块m4.1：由高分辨率影像i波段及低分辨率影像i波段得到灰度级的高分辨率影像i波段的连续概率密度函数以及低分辨率影像i波段的连续概率密度函数，并计算具有如下特征的随机变量；其中，表示积分假变量；表示影像灰度级总数；表示连续概率密度函数的积分上限；表示高分辨率影像i波段连续概率密度函数在灰度值w处的概率密度值；模块m4.2：基于预设的概率密度函数求得变换函数；其中，表示积分假变量；表示连续概率密度函数的积分上限，表示低分辨率影像i波段的连续概率密度函数在灰度值处的概率密度值；模块m4.3：求得反变换函数；所述模块m5采用：其中，为加权计算得到的新强度波段波段在影像坐标处的值；为高分辨率影像波段在影像坐标处的值；为低分辨率影像波段在影像坐标处的值；为影像坐标高分辨率影像波段的计算权重；所述模块m6采用：hsi颜色空间变换公式为：若且
若若若若。

技术总结
本发明提供了一种基于纹理显著性赋权的遥感影像融合方法及系统，包括：步骤S1：对高分辨率影像进行纹理显著性计算得到纹理显著性图像；步骤S2：将纹理显著性图像归一化至0-1范围并进行分段权重运算得到权重矩阵；步骤S3：将高分辨率影像和低分辨率影像转换至HSI颜色空间并提取强度信息I波段；步骤S4：高分辨率影像I波段与低分辨率影像I波段进行灰度直方图匹配；步骤S5：直方图匹配后的高分辨率影像I波段与低分辨率影像I波段基于权重矩阵进行加权运算得到新强度波段；步骤S6：使用新强度波段替换I波段，将HSI空间图像逆转换至RGB空间得到融合影像，完成基于纹理显著性赋权的遥感影像融合。像融合。像融合。


技术研发人员：李现虎 杭君 王婧 郑顺安 周京一
受保护的技术使用者：上海航天空间技术有限公司
技术研发日：2023.08.18
技术公布日：2023/9/20