[发明专利]全域尺度的岩溶石漠化信息遥感快速提取方法

权利要求书

1.全域尺度的岩溶石漠化信息遥感快速提取方法，其特征在于：包括如下步骤：

S1、编制南方岩溶区碳酸盐岩分布图：

以南方地区的区域地质资料为依据，根据岩石矿物成分、结构、构造编制南方岩溶区碳酸盐岩分布图；

S2、植被覆盖度FVC定量反演：

选择覆盖南方岩溶区全域的MODIS数据为遥感数据源，采用二分线性像元分解模型定量反演全域的植被覆盖度，得到植被覆盖度表达式为：

FVC＝(NDVI-NDVI_min)/(NDVI_max-NDVI_min)

式中NDVI_min为累计频率为5％时的NDVI值，NDVI_max为累计频率为95％时的NDVI值；

S3、计算植被覆盖演变指数：

以步骤S2是我遥感数据的反演出的FVC值为参数，采用差值法计算不同年份之间的植被覆盖演变指数，具体的表达式为：

FVC_d＝FVC_m-FVC_n

式中FVC_d为n年至m年的植被覆盖度演变指数，其中m的年份在n年份之后，FVC_m为m年的植被覆盖度，FVC_n为n年的植被覆盖度；

S4、剔除特殊地类：

以每年22期植被覆盖度的均值和方差为参数剔除包括水面、城镇在内的非石漠化的特殊地类；

S5、建立岩溶石漠化演变概念模型：

将经过步骤S4除特殊地类后的FVC_d值分级，将FVC_d值在

[-1，-0.6)区间的像元定义为严重恶化

[-0.6，-0.2)区间的像元定义为轻度恶化

[-0.2，0.2)区间的像元定义为基本未变

[0.2，0.6)区间的像元定义为轻度改善

[0.6，1)区间的像元定义为较好改善；

S6、建立岩溶石漠化演变数学模型，实现全域尺度岩溶石漠化快速提取：

将岩溶石漠化演变概念模型转换为方便计算的数学模型，具体是把n年石漠化现状中的无石漠化赋为正整数11、轻度石漠化赋为正整数13、中度石漠化赋为正整数15、重度石漠化赋为正整数17；把n年至m年的植被覆盖度演变指数FVC_d分级中的严重恶化赋为正整数的1、轻度恶化赋为正整数的3、基本未变赋为正整数的5、轻度改善赋为正整数的7、较好改善赋为正整数的9；将两者相乘得到m年的石漠化程度。

2.根据权利要求1所述的全域尺度的岩溶石漠化信息遥感快速提取方法，其特征在于：步骤S2中，根据二分线性像元分解模型计算FVC值，其中小于累计频率5％的NDVI值的像元FVC值赋为0，大于累计频率95％的NDVI值的像元FVC值赋为1，FVC值为正数为0到1之间的正数，计算后的FVC值的小代表植被覆盖的程度。

3.根据权利要求1所述的全域尺度的岩溶石漠化信息遥感快速提取方法，其特征在于：在步骤S4中，根据FVC的取值范围为0到1，相应的FVC_d的取值范围为-1到1，其负值代表植被覆盖度降低，0附近代表植被覆盖度基本未变，正值代表植被覆盖度增加。