GEE批量提取站点DN值—以Landsat 8 C2 SR 产品为例

PS：首先声明一下，本人是GEE 0级选手，帖子是参考其它技术贴消化编辑的。主要是为了记录数据处理的过程，欢迎交流。

1. 首先导入站点矢量文件。切记，要注意矢量文件的坐标系，做好与影像的坐标系保持一致！！！
   

2. 加载矢量文件至脚本：Assets 中找到导入的矢量站点文件，点击箭头导入
   

3. 上代码

var table: table  Users/DN/samples
Map.addLayer(table)
var geometry = table.geometry()
var L8dataset = ee.ImageCollection('LANDSAT/LC08/C02/T1_L2')
                  .filterBounds(area)
                  .filterDate('2018-03-01','2018-05-01')
                  .filter(ee.Filter.lt('CLOUD_COVER', 80))
                  .select('SR_B1','SR_B2','SR_B3','SR_B4','SR_B5','SR_B6','SR_B7')//
               //   .map(maskL8sr);
print(L8dataset)
var image = L8dataset.mean().clip(geometry)
var pixels = image.reduceRegion({
  reducer: ee.Reducer.toList(),
  geometry: table,
  scale: 100
})
var numberOfPixels = ee.List(pixels.values().get(0)).size()
var bandNames = image.bandNames()
var features = ee.FeatureCollection(
  ee.List.sequence(0, numberOfPixels.subtract(1))
    .map(function (i) {
      return bandNames.iterate(function (bandName, feature) {
        bandName = ee.String(bandName)
        var pixelValue = ee.List(pixels.get(bandName)).get(i)
        return ee.Feature(feature)
          .set(bandName, pixelValue)
      }, ee.Feature(ee.Geometry.Point([0, 0])))
    })
)

print(features)
Export.table.toDrive({
  collection: features,
  description:'pix_bands3',
  fileFormat: 'CSV'
})

4. 运行完成，Run—下载
   

参考帖子：

1. https://zhuanlan.zhihu.com/p/143963044
2. https://zhuanlan.zhihu.com/p/479033806

致谢：

一路跌跌撞撞，感谢某贾姓小师妹的帮助。