1. 案例目标
在GEE中实现提取特定点位的栅格值,这一功能有很大的用途,通常为地面站点或者地面采样点的位置,从特定的产品集中提取数据。例如:多光谱反射率,气象参数,植被指数等等。也介绍了像元二分模型计算植被盖度FVC的方法。(max:95%, min:5%)
下面是具体的思路解释,代码也有详细的注释,有点懒没有对代码拆分进行介绍。
- 首先导入研究区的ROI和研究区域的点矢量。
因为是完整的代码,定义参数的时候把需要改的单独列出来了,可以将变量代入具体的数据集进行更改。 - 根据自己的需求完善配色盘,也可以自己定义。
- 制作数据集,案例提取了MODIS的7光谱,NDVI,EVI,LAI,FVC并将所有数据整合到一块。(Figure 1)
3.1. 使用的是均值mean的植被指数
3.2. 数据均需要乘以0.0001是因为存储数据的时候需要节省空间,使用的都是整数型。我们需要从数据集查询scale factor是多少,然后乘以这个scale
3.3 遥感影像,需要使用质量控制波段去云,代码中有去云函数
3.4 计算FVC的代码也贴了上去,使用的是像元二分模型计算的FVC。马里兰的GLASS数据集应该有近40年的FVC产品,可能准确些吧! - 提取值到点
4.1 使用的是sampleRegions,具体含义可以在Docs查询。简要说明: collection是点矢量,properties是需要导出的属性(点矢量可能属性挺多),scale是提取点的时候检索的栅格的范围大小,对应数据集的栅格大小,landast就改为30,sentinel-2就是10
4.2 reduceRegions貌似也可以,但是我不会 - 将结果下载到Google Drive



2. 完整案例
该案例以欧亚大陆为ROI,进行部分点数据集的提取
//配色盘,可以自行定义,用于图层色阶的展示
var palettes = require('users/gena/packages:palettes');
var palette = palettes.misc.cubehelix[7];
// 定义了函数,这部分是批量化的时候需要改变的参数
var roi = region
var output_name = 'export result'
var point = samples
var start_time = '2018-07-15'
var end_time = '2018-08-15'
var features = ['Sample_ID','Latitude','Longitude','NDVI','EVI','LAI',
'b1','b2','b3','b4','b5','b6','b7', 'FVC']
/***************************************************************
1.定义样本点(经纬度信息)
***************************************************************/
// var roi = region
// var point = sample
Map.centerObject(point,4)
Map.addLayer(point,{
color:'red'},'sample_point')
print(point)
/***************************************************************
2.制作数据集,从GEE中提取需要的产品集。(以MODIS为例进行介绍)
***************************************************************/
// 2.1 选择MODIS的植被指数产品——VIs
// var start_time = '2020-7-01'
// var end_time = '2020-08-31'
var collection_mod1 = ee.ImageCollection("MODIS/061/MOD13A2")
.filterBounds(roi)
.filterDate(start_time, end_time);
//select feature (NDVI,EVI)
var VIs = collection_mod1.select(["NDVI","EVI"</

最低0.47元/天 解锁文章
962

被折叠的 条评论
为什么被折叠?



