[发明专利]获取耕地土壤盐渍化程度等级分布的方法和装置

权利要求书

1.一种获取耕地土壤盐渍化程度等级分布的方法，其特征在于，所述方法包括：

获取目标区域的卫星影像中组成耕地区域的耕地像元；

根据每个所述耕地像元的不同波段光线的反射率获取每个所述耕地像元的植被指数和盐分指数；

构建所述耕地区域关于所述植被指数和所述盐分指数的二维特征空间；

根据所述二维特征空间内的所述植被指数和所述盐分指数的分布获取每个所述耕地像元的土壤电导率；

根据每个所述耕地像元的土壤电导率获取每个所述耕地像元对应的土壤含盐量；

根据每个所述耕地像元的土壤含盐量确定每个所述耕地像元的盐渍化程度等级。

2.根据权利要求1所述的获取耕地土壤盐渍化程度等级分布的方法，其特征在于，所述获取目标区域的卫星影像中组成耕地区域的耕地像元包括：

获取所述目标区域的卫星影像；

对所述目标区域的卫星影像进行预处理，得到预处理后的目标区域的卫星影像，所述预处理至少包括几何校正或大气校正；

通过最大似然监督分类法确定所述预处理后的目标区域的卫星影像中每个像元的土地利用类型，得到所述目标区域的卫星影像中组成耕地区域的各个耕地像元。

3.根据权利要求1所述的获取耕地土壤盐渍化程度等级分布的方法，其特征在于，所述根据每个所述耕地像元的不同波段光线的反射率获取每个所述耕地像元的植被指数和盐分指数包括：

根据下式获取每个所述耕地像元的植被指数和盐分指数：

和

式中，MSAVI和SI分别为所述耕地像元的植被指数和盐分指数；G、R和NIR分别为所述耕地像元对应的绿光波段、红光波段和近红外波段的反射率值。

4.根据权利要求1所述的获取耕地土壤盐渍化程度等级分布的方法，其特征在于，在所述根据每个所述耕地像元的不同波段光线的反射率获取每个所述耕地像元的植被指数和盐分指数后，根据下式对所述植被指数和所述盐分指数进行正规化处理：

和

式中，MSAVI和SI分别为所述耕地像元的植被指数和盐分指数；MSAVI′和SI′分别为经过正规化处理后的所述植被指数和所述盐分指数；MSAVI_max和MSAVI_min分别为所述耕地像元的植被指数的最大值和最小值；SI_max和SI_min分别为所述耕地像元的盐分指数的最大值和最小值。

5.根据权利要求1所述的获取耕地土壤盐渍化程度等级分布的方法，其特征在于，所述构建所述耕地区域关于所述植被指数和所述盐分指数的二维特征空间包括：

根据正规化处理后的植被指数和盐分指数构建所述耕地区域关于所述植被指数和所述盐分指数的二维特征空间，所述二维特征空间的横轴为所述正规化处理后的植被指数，所述二维特征空间的纵轴为所述正规化处理后的盐分指数。

6.根据权利要求1所述的获取耕地土壤盐渍化程度等级分布的方法，其特征在于，所述根据所述二维特征空间内的所述植被指数和所述盐分指数的分布获取每个所述耕地像元的土壤电导率包括：

根据所述耕地像元在所述二维特征空间内的分布拟合出所述耕地区域的土壤基线：

SI'＝kMSAVI'+b

式中，MSAVI′和SI′分别为经过正规化处理后的所述植被指数和所述盐分指数，k为所述土壤基线的斜率，b为常量；

过所述二维特征空间横轴上的固定坐标点作所述土壤基线的垂线，所述固定坐标点为(1,0)；

根据下式获取所述耕地像元到所述耕地区域的土壤基线的垂线的距离：

式中，E_MSAVI-SI为所述耕地像元到所述耕地区域的土壤基线的垂线的距离；

根据下式获取每个所述耕地像元的土壤电导率：

式中，EC为所述耕地像元的土壤电导率；m和n为常量。