GEE案例:根据点位提取栅格值(acrgis的多值提取到点)

1. 案例目标

在GEE中实现提取特定点位的栅格值,这一功能有很大的用途,通常为地面站点或者地面采样点的位置,从特定的产品集中提取数据。例如:多光谱反射率,气象参数,植被指数等等。也介绍了像元二分模型计算植被盖度FVC的方法。(max:95%, min:5%)

下面是具体的思路解释,代码也有详细的注释,有点懒没有对代码拆分进行介绍。

  1. 首先导入研究区的ROI和研究区域的点矢量。
    因为是完整的代码,定义参数的时候把需要改的单独列出来了,可以将变量代入具体的数据集进行更改。
  2. 根据自己的需求完善配色盘,也可以自己定义。
  3. 制作数据集,案例提取了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. 提取值到点
    4.1 使用的是sampleRegions,具体含义可以在Docs查询。简要说明: collection是点矢量,properties是需要导出的属性(点矢量可能属性挺多),scale是提取点的时候检索的栅格的范围大小,对应数据集的栅格大小,landast就改为30,sentinel-2就是10
    4.2 reduceRegions貌似也可以,但是我不会
  5. 将结果下载到Google Drive
    Figure 1.
    在这里插入图片描述

在这里插入图片描述

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"</
评论 21
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值