全国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值。

相关推荐
Cx330❀10 分钟前
【LangChain】LangChain 核心技术全景指南:从基础入门到 LCEL 链式编程
大数据·elasticsearch·搜索引擎·性能优化·langchain·全文检索
xiaoye-duck13 分钟前
《Linux网络编程》TCP 网络编程(下):基于 TCP 多线程通用字典服务 + 远程命令执行实战
linux·服务器·tcp/ip
陈让然1 小时前
VS Code C/C++ 插件 cpptools 内存占用异常
linux·运维
字节数据平台3 小时前
iDA:从 ChatBI 到专业数据分析助手的演进之路
大数据·人工智能·机器学习·数据分析
2401_894915533 小时前
GEO 优化源码全解析:从搜索引擎到 AI 引擎的底层改写逻辑
java·服务器·前端·数据库·人工智能·分布式·搜索引擎
qq_316411038 小时前
商用冷柜售后报修运维管理系统开发案例
运维
huainingning8 小时前
R4930-G7服务器(RAID-LSI-9560-LP-8i-4GB)网页方式做Raid方法
服务器
智码看视界8 小时前
Day59-阿里云ACK实战:生产级K8s集群部署与运维
运维·阿里云·kubernetes·k8s·日志监控·容器运维·阿里云ack
正儿八经的少年8 小时前
布隆过滤器(解决redis缓存穿透步骤之一)
数据库·redis·缓存
万岳软件开发小城8 小时前
同城O2O系统源码核心架构分析:多商户、订单、支付及运营体系设计
大数据·同城外卖系统源码·同城o2o系统源码·外卖app开发·同城o2o小程序开发·同城o2o软件开发