[发明专利]一种植被覆盖管理因子的动态评估方法

摘要

本发明提供一种植被覆盖管理因子的动态评估方法，涉及水土保持领域。一种植被覆盖管理因子的动态评估方法，主要包括以下步骤：S1多源遥感影像数据的收集、S2D2的预处理、S3计算D2的植被覆盖度、S4计算每幅D1的植被覆盖度、S5计算每幅D1的植被覆盖度差异、S6计算高精度月MODND1T数据的植被覆盖度、S7计算月降水量栅格数据、S8计算月降雨侵蚀力和年降雨侵蚀力、S9计算月降雨侵蚀力比、S10计算月土壤流失率、S11植被覆盖管理因子集的获取。本发明的一种植被覆盖管理因子的动态评估方法，实现了MODIS影像的产品数据与高分2号卫星影像数据融合，得到了具有高时空分辨率的植被覆盖管理因子集。

主权项

1.一种植被覆盖管理因子的动态评估方法，其特征在于，主要包括以下步骤：S1、多源遥感影像数据的收集：下载研究区域目标研究年份中每月的MODND1T影像数据D1多幅以及高分2号卫星影像数据D2一幅；S2、D2的预处理：采用ENVI软件对S1中的所述D2进行预处理；S3、计算D2的植被覆盖度：采用ENVI软件中的MNF变换，提高S2中的D2的信噪比，利用PPI法提取D2的纯净像元并获取典型地物光谱，利用LSU计算地物中植被的丰度，得到D2的植被覆盖度；S4、计算每幅D1的植被覆盖度：采用ENVI软件中的像元二分模型分别将S1中的每幅所述D1中的NDVI转化为植被覆盖度，以分别得到每幅D1的植被覆盖度；S5、计算每幅D1的植被覆盖度差异：利用ArcGIS软件中的栅格计算器工具箱将S4的多幅所述D1的植被覆盖度分别与S3中的所述D2的植被覆盖度进行做差运算，得到目标研究年份一年中每幅D1的植被覆盖度差异；S6、计算高精度月MODND1T数据的植被覆盖度：将S5中的多幅D1的植被覆盖度差异分别与S3中的所述D2的植被覆盖度进行相加运算，即得到每月的高精度月MODND1T数据的植被覆盖度。S7、计算月降水量栅格数据：下载研究区域目标研究年份中每月的月降水量数据，利用ArcGIS软件中的克里格插值法对研究区域目标研究年份中每月的月降水量数据进行空间插值，得到研究区域目标研究年份中每月的月降水量栅格数据；S8、计算月降雨侵蚀力和年降雨侵蚀力：根据S7中所述月降水量栅格数据计算研究区域目标研究年份中每月的月降雨侵蚀力和年降雨侵蚀力，其公式为：这式中，R_年‑‑年降雨侵蚀力，MJ·mm/·hm²·h；0.0125P_t^1.6295‑‑某月的月降雨侵蚀力，MJ·mm/·hm²·h；P_t‑‑某月侵蚀性降水总量(≥12mm)，mm；t‑‑月份；S9、计算月降雨侵蚀力比：根据S8中的所述月降雨侵蚀力和所述年降雨侵蚀力，计算研究区域目标研究年份中的每月的月降雨侵蚀力比，其计算方法为：月降雨侵蚀力比＝每月降雨侵蚀力/年降雨侵蚀力；S10、计算月土壤流失率：利用ArcGIS软件中的栅格计算器将S9中的每月的所述月降雨侵蚀力比与S6中得到的对应月份的所述高精度月MODND1T数据的植被覆盖度一一对应的进行加权运算，得到对应月份的月土壤流失率，其运算公式为Y＝log(‑0.043*X),这式中，Y‑‑月土壤流失率，％；X‑‑高精度月MODND1T数据的植被覆盖度，％；S11、植被覆盖管理因子集的获取：利用ArcGIS软件中的栅格计算器将S10中的每月的所述月土壤流失率与S8中得到的对应月份的所述月降雨侵蚀力一一对应的进行加权运算，得到多个月植被覆盖管理因子，并将多个所述月植被覆盖管理因子进行整理，即可得到研究区域目标研究年份的植被覆盖管理因子集，其中，运算公式为C＝SLR₁*Ratio₁+SLR₂*Ratio₂+SLR_n*Ratio_n+…+SLR₁₂*Ratio₁₂，这式中，C‑‑植被覆盖管理因子；n‑‑月份；SLR_n‑‑第n个月的土壤流失率，％；Ratio_n‑‑第n个月的降雨侵蚀力比例，％。