[发明专利]基于Sentinel-2光学卫星影像的蓝色参棚遥感识别方法

权利要求书

1.一种基于Sentinel-2光学卫星影像的蓝色参棚遥感识别方法，其特征在于，其步骤如下：

S1、收集蓝色参棚生命周期内的任意一期Sentinel-2光学卫星影像，并对Sentinel-2光学卫星影像进行预处理，得到Sentinel-2地表反射率图像；

S2、利用GNSS定位工具采集蓝色参棚和其他地物的地理坐标，然后利用光谱采集工具在步骤S1中得到的Sentinel-2地表反射率图像中采集各个地理坐标位置上的蓝色参棚和其他地物的光谱样本；

S3、分析步骤S2中得到的光谱样本，确定蓝色参棚与其他地物不同的有效光谱波段，即蓝色波段、红色波段、第二红边波段、近红外波段、第一短波红外波段、第二短波红外波段；

S4、基于步骤S3中得到的有效光谱波段，构造蓝色参棚第一遥感指数和蓝色参棚第二遥感指数；

S5、基于步骤S1中得到的Sentinel-2地表反射率图像和S4中得到的蓝色参棚第一遥感指数计算得到蓝色参棚第一遥感指数图像；

S6、基于步骤S1中得到的Sentinel-2地表反射率图像和S4中得到的蓝色参棚第二遥感指数计算得到蓝色参棚第二遥感指数图像；

S7、对步骤S1中得到的Sentinel-2地表反射率图像中的待识别区域进行蓝色参棚识别，如果待识别区域中像元位置i上的红色波段上的像元值小于第一阈值γ，且近红外波段上的像元值大于第二阈值δ，那么执行步骤S8，否则，该像元位置的属性赋值为其他地物；

S8、如果待识别区域中像元位置i上的蓝色参棚第一遥感指数图像上的像元值大于第三阈值α，那么执行步骤S9，否则，该像元位置的属性赋值为其他地物；

S9、如果待识别区域中像元位置i上的蓝色参棚第二遥感指数图像上的像元值大于第四阈值β，那么该像元位置的属性为蓝色参棚，否则，该像元位置的属性赋值为蓝色屋顶建筑；

S10、循环执行步骤S7至步骤S9，直至遍历完待识别区域中的所有像元位置，即完成蓝色参棚的遥感自动识别。

2.根据权利要求1所述的基于Sentinel-2光学卫星影像的蓝色参棚遥感识别方法，其特征在于，所述第一阈值γ的确定方法为：对步骤S2中得到的蓝色参棚光谱样本进行统计分析，得到蓝色参棚在红色波段上的像元值小于γ；第二阈值δ的确定方法为：对步骤S2中得到的蓝色参棚光谱样本进行统计分析，得到蓝色参棚在近红外波段上的像元值的大于δ；第三阈值α的确定方法为：基于步骤S2中得到的蓝色参棚光谱样本和步骤S4中得到的蓝色参棚第一遥感指数，计算蓝色参棚在蓝色参棚第一遥感指数中的取值范围，即蓝色参棚的蓝色参棚第一遥感指数大于α；第四阈值β的确定方法为：基于步骤S2中得到的蓝色参棚光谱样本和步骤S4中得到的蓝色参棚第二遥感指数，计算蓝色参棚在蓝色参棚第二遥感指数中的取值范围，即蓝色参棚的蓝色参棚第二遥感指数大于β。

3.根据权利要求1所述的基于Sentinel-2光学卫星影像的蓝色参棚遥感识别方法，其特征在于，所述蓝色参棚第一遥感指数的表达式为：

式中，BGS₁表示蓝色参棚第一遥感指数值，ζ_blue表示蓝色波段的地表反射率，ζ_red表示红色波段的地表反射率，ζ_redEdge2表示第二红边波段的地表反射率。

4.根据权利要求1所述的基于Sentinel-2光学卫星影像的蓝色参棚遥感识别方法，其特征在于，所述蓝色参棚第二遥感指数的表达式为：

式中，BGS₂表示蓝色参棚第二遥感指数值，ζ_nir表示近红外波段的地表反射率，ζ_swir1表示第一短波红外波段的地表反射率，ζ_swir2表示第二短波红外波段的地表反射率。