本文基于吉林一号遥感卫星高空间分辨率数据，选择适用于荒漠草原的NDVI、MSAVI、RVI和PVI四种植被指数，结合地面调查数据，对牧户尺度植被覆盖度进行反演模型的构建和验证。通过对比分析得出，基于MSAVI反演模型优于NDVI、RVI和PVI；其中，4 × 4 m样方拟合精度最高，EA为48.49%，优于1 × 1 m和2 × 2 m样方。 This paper is based on the high spatial resolution data of Jl1 remote sensing satellite. Four vegetation indices, NDVI, MSAVI, RVI and PVI, suitable for desert grasslands, were selected and combined with the ground survey data to construct and verify vegetation coverage inversion model in herds-scale. Through comparative analysis, the inversion model based on MSAVI is superior to NDVI, RVI and PVI. Among them, 4 × 4 m quadrat had the highest fitting accuracy, with EA of 48.49%, which was better than 1 × 1 m and 2 × 2 m quadrat.
本文基于吉林一号遥感卫星高空间分辨率数据，选择适用于荒漠草原的NDVI、MSAVI、RVI和PVI四种植被指数，结合地面调查数据，对牧户尺度植被覆盖度进行反演模型的构建和验证。通过对比分析得出，基于MSAVI反演模型优于NDVI、RVI和PVI；其中，4 × 4 m样方拟合精度最高，EA为48.49%，优于1 × 1 m和2 × 2 m样方。
植被指数，遥感反演，植被覆盖度，吉林一号卫星数据
Youzhi Wang, E’erdemutu Jin*, Sarenna Wang
College of Geographical Science, Inner Mongolia Normal University, Hohhot Inner Mongolia
Received: Nov. 29th, 2023; accepted: Feb. 22nd, 2024; published: Feb. 28th, 2024
This paper is based on the high spatial resolution data of Jl1 remote sensing satellite. Four vegetation indices, NDVI, MSAVI, RVI and PVI, suitable for desert grasslands, were selected and combined with the ground survey data to construct and verify vegetation coverage inversion model in herds-scale. Through comparative analysis, the inversion model based on MSAVI is superior to NDVI, RVI and PVI. Among them, 4 × 4 m quadrat had the highest fitting accuracy, with EA of 48.49%, which was better than 1 × 1 m and 2 × 2 m quadrat.
Keywords:Vegetation Index, Remote Sensing Inversion, Vegetation Coverage, JL1 Satellite Data
Copyright © 2024 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
荒漠草原是草原向荒漠过渡的脆弱生态系统，对气候变化和人类活动的干扰具有敏感性 [
研究区位于内蒙古自治区乌兰察布市四子王旗查干补力格苏木格日乐图雅嘎查，地理位置41˚41′14.97″N ，111˚46′43.77″E，牧户面积230 hm2。草场平均海拔1040 m，为波状起伏地貌，气候属于中温带半干旱大陆性气候，草地类型为短花针茅 + 冷蒿 + 无芒隐子草荒漠草原，植被覆盖度较低。
遥感数据为吉林一号遥感卫星2022年8月16日无云数据(Level 1T)。多光谱数据空间分辨率为2 m，全色波段为0.5 m。地面调查数据为2022年7月28、29日采集的植被盖度数据，共采集70个样点、210个样方数据，每个样点分别采集1 m × 1 m、2 m × 2 m、4 m × 4 m样方3组数据，每个样点使用照相机法拍摄样方照片，手持GPS记录样点地理坐标。研究区样点分布如图1所示。
1) 数据预处理
在ENVI软件中对遥感数据进行辐射校正、大气校正、多光谱和全色波段的融合处理，得到亚米级高空间分辨率多光谱数据；在ArcGIS软件中对融合数据进行投影和几何纠正。
EXCEL中录入地面调查数据的地理坐标、海拔高度、植被覆盖度等属性数据；在ARCGIS软件中将属性数据转成空间数据。
图1. 研究区样点分布图
2) 植被指数运算
结合研究区草地植被覆盖的特征，选取适用于监测中低覆盖度植被的归一化植被指数(NDVI)、修改型调节植被指数(MSAVI)、比值植被指数(RVI)和垂直植被指数(PVI)四种植被指数模型，基于四种植被指数构建植被覆盖度反演模型。植被指数及其计算公式如下，其中NIR为近红外波段，RED为红外波段。
a) NDVI (Normalized Difference Vegetation Index)公式如下：
NDVI = ( NIR − RED ) / ( NIR + RED ) (1)
b) MSAVI (Modified Soil Adjusted Vegetation Index)公式如下：
MSAVI = ( 2 NIR + 1 − ( 2 NIR + 1 ) 2 − 8 ( NIR − RED ) ) / 2 (2)
c) 比值植被指数RVI (Relative Vegetation Index)公式如下：
RVI = NIR / RED (3)
d) PVI (Perpendicular Vegetation Index)公式如下：
PVI = ( NIR − a RED − b ) / a 2 + 1 (4)
3) 精度验证
从样方数据中随机抽样70%的地面数据与植被指数进行回归模型的建立，并利用未参与模型建立的30%样方作为验证数据进行精度计算，计算公式如下：
R 2 = ∑ i = 1 n ( X ⌢ i − X ¯ ) 2 ∑ i = 1 n ( X i − X ¯ ) 2 (5)
RMSE = ∑ i = 1 n ( X ⌢ i − X i ) 2 N (6)
EA = ( 1 − ∑ i = 1 n ( X ⌢ i − X i ) 2 N X ¯ ) × 100 % (7)
回归系数(R2)、均方根误差(RMSE)和精度(EA)公式中，N表示检验样本的数量， X i 表示实测值， X ⌢ i 表示预测值， X ¯ 表示实测值的均值。
选取的样点分别与像元NDVI、MSAVI、RVI和PVI构建样方1 m、2 m和4 m植被覆盖度的线性反演模型。
表1为四种植被指数与1 m样方植被盖度所建立的回归方程。可以看出，植被指数与1 m样方植被覆盖度回归模型中基于MSAVI建立的线性回归模型 y = 33.66 x + 6.0177 拟合度最好，R2值为0.1141。其次为PVI建立的线性回归模型 y = 0.0305 x + 8.8659 ，R2值为0.0849。
| 植被指数 | 曲线形式 | 回归方程 | R2 |
|---|---|---|---|
| NDVI | 线性 | y = 34.493x + 10.768 | 0.0828 |
| MSAVI | y = 33.66x + 6.0177 | 0.1141 | |
| RVI | y = 5.6304x + 10.026 | 0.0417 | |
| PVI | y = 0.0305x + 8.8659 | 0.0849 |
表1. 基于四种植被指数的植被覆盖度回归模型及拟合精度(1 m样方)
表2为四种植被指数与2 m样方植被盖度所建立的回归方程。可以看出，植被指数与2 m样方植被覆盖度回归模型中基于MSAVI建立的线性回归模型 y = 41.073 x + 0.2271 拟合度最高，R2值为0.2915。其次为NDVI建立的线性回归模型 y = 45.301 x + 5.1597 ，R2值为0.2449。
| 植被指数 | 曲线形式 | 回归方程 | R2 |
|---|---|---|---|
| NDVI | 线性 | y = 45.301x + 5.1597 | 0.2449 |
| MSAVI | y = 41.073x + 0.2271 | 0.2915 | |
| RVI | y = 8.6709x + 1.9127 | 0.1695 | |
| PVI | y = 0.0334x + 5.0981 | 0.1749 |
表2. 基于四种植被指数的植被覆盖度回归模型及拟合精度(2 m样方)
表3为四种植被指数与4 m样方植被盖度所建立的回归方程。从表中可以看出，植被指数与4 m样方植被覆盖度回归模型中基于MSAVI建立的线性回归模型y = 41.015x + 5.1516拟合精度最高，R2值为0.1766。其次为NDVI建立的线性回归模型y = 42.864x + 10.716，R2值为0.1332。
| 植被指数 | 曲线形式 | 回归方程 | R2 |
|---|---|---|---|
| NDVI | 线性 | y = 42.864x + 10.716 | 0.1332 |
| MSAVI | y = 41.015x + 5.1516 | 0.1766 | |
| RVI | y = 7.272x + 9.3033 | 0.0724 | |
| PVI | y = 0.0314x + 10.759 | 0.0934 |
表3. 基于四种植被指数的植被覆盖度回归模型及拟合精度(4 m样方)
将预留的样点数据分别代入各样方植被覆盖度拟合精度最高的回归模型中，1 m样方植被覆盖度拟合精度最高的回归模型为FVC = 33.66MSAVI + 6.0177，2 m样方植被覆盖度拟合精度最高的回归模型为FVC = 41.073MSAVI + 0.2271，4 m样方植被覆盖度拟合精度最高的回归模型为FVC = 41.015MSAVI + 5.1516。结果表明(表4)，1 m样方植被覆盖度估算模型预测R2为0.0713，RMSE为10.32，EA为42.69%；2 m样方植被覆盖度估算模型预测R2为0.072，RMSE为10.44，EA为35.56%；4 m样方植被覆盖度估算模型预测R2为0.0851，RMSE为11.13，EA为48.49%；三种样方中4 m样方的植被覆盖度估算模型拟合精度最高。
| 样方植被覆盖度类型 | 验证R2 | 验证RMSE | 精度(EA) |
|---|---|---|---|
| 1m植被覆盖度 | 0.0713 | 10.32 | 42.69% |
| 2m植被覆盖度 | 0.072 | 10.44 | 35.56% |
| 4m植被覆盖度 | 0.0851 | 11.13 | 48.49% |
表4. 植被覆盖度估算模型验证表
对牧户草场植被盖度划分五个等级，分别是0%~20%为低覆盖度，20%~40%为中低覆盖度，40%~60%为中覆盖度，60%~80%为高覆盖度，>80%为高覆盖度。基于模型估算结果分析可得(图2)，研究区中低值被覆盖度分布面积最大，约占总面积的61%；其次是中植被覆盖度分布面积，约占总面积的37%；而中低植被覆盖度、中高覆盖度和高覆盖面积较小，占比均不到2%。因此，草场植被覆盖较差，主要以低植被覆盖为主，是典型的荒漠草原类型的草场。
吉林一号卫星数据具有高空间分辨率的优势，适于反演小尺度植被覆盖情况。在荒漠草原牧户尺度草场研究区，构建NDVI、MSAVI、RVI和PVI四种植被指数回归模型，对比得出基于吉林一号卫星数据的MSAVI反演模型适于研究区植被覆盖度的估算，精度验证得出4 m样方估算模型精度最高。
图2. 研究区植被覆盖度空间分布图
在反演植被覆盖度研究中，应用广泛的是归一化差值植被指数(NDVI)，其中利用基于NDVI的像元二分模型的应用研究较多 [
本文章由“内蒙古师范大学引进高层次人才科研启动经费项目(2016YJRC016)”和“内蒙古师范大学校基金项目(2015YBXM014)”基金资助。
王由志,金额尔德木吐,王萨仁娜. 基于吉林一号卫星数据的荒漠草原植被覆盖度遥感反演研究Remote Sensing Inversion Study on Fractional Vegetation Coverage of Desert Grassland Based on JL1 Satellite Data[J]. 地理科学研究, 2024, 13(01): 34-40. https://doi.org/10.12677/GSER.2024.131004
https://doi.org/10.2111/REM-D-09-00033.1
https://link.cnki.net/urlid/11.3362.s.20231025.1059.002