全国30米分辨率逐年植被覆盖度(FVC)数据集

本数据集包括1986至2023年间,全国植被覆盖度数据,FVC范围值为0-1,数据为浮点型,GeoTIFF格式。GeoTIFF文件均可用ArcGIS软件和GDAL读取和打开。数据源来自LANDSAT系列遥感卫星,通过Google Earth Engine云计算平台进行处理(包括去云、云影去除、以及条带错位的修复)获得。本数据集以30米的空间分辨率捕捉植被的生长状态,为从长时间尺度观测植被变化提供了宝贵的数据资源。

植被覆盖度是指植被(包括叶、茎、枝)在地面的垂直投影面积占统计区总面积的百分比。是刻画地表植被覆盖的一个重要参数, 也是指示生态环境变化的重要指标之一。植被覆盖度是描述地表植被的重要参数与生态环境的基本指标,在大气圈、土壤圈、水圈和生物圈中都占据着重要的地位,区域及全球范围的植被覆盖度估算对植被及相关领域的研究具有重要意义。

目前在多个领域内,植被覆盖度进行分级处理后使用。往往参照《土壤侵蚀分类分级标准》,将植被覆盖度划分为4个等级:<30%为低覆盖度、30%---45%为中低覆盖度、45%---60%为中覆盖度、>60%为高覆盖度。 该数据利用Landsat数据,根据植被盖度与NDVI之间的经验关系计算得到。

该产品时间分辨率为1年,空间分辨率30米。算法从当年所有观测数据中,以低云、低探测角度和最高NDVI值为标准,选择最佳的可用像元值,并进行换算。

数据采集和处理

数据集整合了LANDSAT系列多型号卫星(如LANDSAT 5 TM, LANDSAT 7 ETM+, LANDSAT 8 OLI/TIRS,LANDSAT9等)的表面反射率产品。采用先进算法对原始遥感影像进行云屏蔽,并基于已处理影像生成FVC。

(1) 计算NDVI

FVC的计算基于NDVI,NDVI最大值,NDVI最小值,计算NDVI代码模块如下:

var get_NDVI = function(image)

{ var NDVI=image.normalizedDifference('B4','B3').rename('NDVI');

image=image.addBands(NDVI) return image.select("NDVI") };

(2) 选取NDVI最小值和最大值

取两个百分位5%和95%分别设置为NDVI的最小值和最大值来模拟NDVI的土壤值和纯植被覆盖值,代码如下:

function calFVC(BestVI,region,scale)

{ var num = BestVI.reduceRegion({ reducer:ee.Reducer.percentile(5,95),

geometry:region, scale:scale, maxPixels:1e13 });

var min = ee.Number(num.get("NDVI_p5"));

var max = ee.Number(num.get("NDVI_p95"));

//print(top_min); //print(top_max);

//quantile and combine var greaterPart = BestVI.gt(max);

var lessPart = BestVI.lt(min);

var middlePart =ee.Image(1).subtract(greaterPart).subtract(lessPart);

//calculate FVC var tempf1=BestVI.subtract(min).divide(max.subtract(min));

(3) FVC的计算

分别对LANDSAT 5/7/8系列数据计算FVC,其中Landsat5 和Landsat 7NDVI计算利用 B4和B3波段,Landsat8 和Landsat9 NDVI计算利用 B5和B4波段。

FVC计算公式为:

FVC=( NDVI-NDVImin)/(NDVImax-NDVImin)

NDVImax 和NDVImin分别为区域内最大和最小的NDVI值。

相关推荐
XiaoMaqqqq6 小时前
市面上正规的IP驱动产业新场景新工具哪家强
运维·python·网络协议·tcp/ip
阿部多瑞 ABU6 小时前
复杂社会关系、三层结构与自反性
大数据·人工智能·ai写作
许彰午6 小时前
01-Oracle数据库简介与版本选择
数据库·oracle
菜鸟~noob2336 小时前
【电子战】第三季预告—波武器与电子系统强力毁伤
开发语言·数据库·matlab·电子战
huainingning6 小时前
RJ SW Console口忘记密码处理方法
linux·运维·服务器
布裘7 小时前
【银河麒麟】桌面系统循环登录,无法进入桌面?
linux·运维·服务器
慧都小项7 小时前
程序到了Linux才出错?用CLion把调试接到目标环境
linux·运维·服务器
鲲穹AI种草7 小时前
小红书 AI 创作工具怎么选?鲲穹 RedNote 功能实测与横向对比
大数据·人工智能
Lsetea7 小时前
OpenSSL verify报unhandled critical extension:error 34与关键扩展排查
运维·https·ssl证书·openssl·证书链
guo_wen_qiang8 小时前
docker的备份与迁移
运维·docker·容器