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

相关推荐
渣渣盟16 小时前
Flink 流式写入文件终极指南:从 Row Format 到 Bucket 分区的深度剖析
大数据·flink
Wang's Blog16 小时前
PostgreSQL笔记33: 索引底层原理、扫描类型与属性体系深度解析
数据库·笔记·postgresql
地衣君16 小时前
从 profile 到 kanban:Hermes 多 Agent 协作解析及示例
网络·数据库·tcp/ip
秋田君16 小时前
【无标题】
数据库·oracle
Logintern0916 小时前
数据库分区表的索引实际上包含本地索引和全局索引两种类型
数据库·分区索引
Elastic 中国社区官方博客16 小时前
跳过 mapping 爆炸:ES|QL 无需动态 mapping 即可查询无 schema JSON key
大数据·人工智能·sql·elasticsearch·搜索引擎·json·全文检索
阿沐沐,16 小时前
PowerShell 里改了 config.toml,Codex CLI 仍用旧值:先分清该改哪一层
java·服务器·数据库·人工智能·ai·ai编程
秣宇17 小时前
银河麒麟服务器操作系统关闭 Swap 分区
linux·运维·服务器·github·kylin
RisunJan17 小时前
Linux命令-usernetctl(已废弃 - 通过 usermode-helper 控制网络接口的包装器)
linux·运维·服务器