【GEE学习日记】栅格提取到点（根据经纬度提取栅格值）

一、背景
前面完成了MODIS的LST影像的加载和裁剪（按掩膜提取），接下来需要提取一些点上的LST值。
使用数据：
（1）MOD11A1
（2）研究区范围面矢量
（3）需要提取的点（可以现场写，也可以导入，后面会写出两种情况的代码）
二、参考文献
《Google Earth Engine（像元值提取）》
《GEE案例：根据点位提取栅格值（acrgis的多值提取到点）》
《二. GEE加载矢量数据》
三、代码
先import研究区矢量图

3.1 现场写2个点

//按掩膜提取函数
function ExtractByMask(Image)
{
  return Image.clip(roi)
}
//乘以参数函数
function Multiply(Image)
{
  return Image.multiply(0.02)
}

var roi=table.geometry() //获得研究区
Map.centerObject(roi,5) //设置地图的显示中心
//选择数据集
var LST_dataset=ee.ImageCollection("MODIS/061/MOD11A1")
          .filterDate('2002-01-01','2003-01-01')
          .select('LST_Day_1km')
          .filterBounds(roi);
print(LST_dataset)
var ft=ee.FeatureCollection(ee.List([])) //创建一个要素集合


var p1 = ee.Geometry.Point([93.36, 32.23]) //第1个点的位置
var p2 = ee.Geometry.Point([93.23, 33.14]) //第2个点的位置
//把两个点合在一起组成要素合集
var pts = ee.FeatureCollection(ee.List([ee.Feature(p1).set('name','p1'),
                               ee.Feature(p2).set('name','p2')]))
//设置栅格值提取到点函数
var fill = function(img, ini) {
  var inift = ee.FeatureCollection(ini)
  var ft2 = img.sampleRegions({
  collection:pts,
  properties:ee.List(['name']),
  //scale:30
  });
  var date = img.date().format()
  var ft3 = ft2.map(function(f){return f.set("date", date)})
  return inift.merge(ft3)
}
//先乘以参数再裁剪
var LST_Multiply=LST_dataset.map(Multiply)
var LST_ExtractByMask=LST_Multiply.map(ExtractByMask)
//在ImageCollection上迭代
var newft = ee.FeatureCollection(LST_dataset.iterate(fill, ft))
print(newft)
//导出表格
Export.table.toDrive({
  collection: newft,
  description: 'sample_get3',
  fileFormat: 'CSV'
});
//设置显示属性
var color={
  min:200,
  max:350,
  palette:[
    'FFFFFF','CE7E45','DF923D','F1B555','FCD163','99B718','74A901',
    '66A000','529400','3E8601','207401','056201','004C00','023B01',
    '012E01','011D01','011301'
  ],
};
Map.addLayer(LST_ExtractByMask,color,'LST') //添加LST图层
Map.addLayer(pts) //添加点
print('OK')

结果如下图所示：

打开云盘

3.2 导入点文件
我的研究区和所需要的点位置如下图所示：

先上传点矢量文件：

四、总结

1. 原始数据集是365张影像，但由于只有193张影像在该点处有值，因此得到的表格只有193行，其余天数都是没有nodata