【GEE平台】FLUXNET批量站点MODIS-LAI下载

德布。

已于 2023-04-19 09:59:35 修改

阅读量938

点赞数 3

文章标签： arcgis 数据分析

于 2023-04-18 16:21:05 首次发布

本文链接：https://blog.csdn.net/weixin_44227199/article/details/130216072

版权

该文介绍了如何存储和处理地理空间数据，包括将站点经纬度存储为CSV，用ArcGIS转换为shp文件，以及在GoogleEarthEngine(GEE)平台上批量下载MODIS数据，特别是LAI和Nadir_Reflectance_Band1参数，并将结果导出为CSV文件。

摘要由CSDN通过智能技术生成

1、存储站点经纬度

将站点的经纬度存储在一个表格或CSV文件中，如图所示：
在这里插入图片描述

2、使用ArcGIS将csv文件转化为shp

添加数据（csv文件格式）—右键—显示XY数据(选择经纬度坐标)—设置坐标系WGS1987—确定—右键—数据—导出
点击图层（勾选，可以在图上展示出散点）—右键—数据—导出数据
注意，x选LONG，y选LAT

3、GEE平台上批量下载站点日尺度的MODIS LAI数据

// Change Date here
var dataset = ee.ImageCollection('MODIS/061/MYD15A2H')
                .select('Lai_500m')                  
                .filter(ee.Filter.date('2007-01-01', '2007-12-31'));
                  
// Upload your point shapefile(Need Change!!)
var sample_point = ee.FeatureCollection("users/luoyi/FIUX");

// Start collect data
var ft = ee.FeatureCollection(ee.List([]))
var fill = function(img, ini) {
  // tran ini 2 FeatureCollection
  var inift = ee.FeatureCollection(ini)
  // Extract pixel value
  var ft2 = img.sampleRegions({
  collection:sample_point,
  // Add "NAME" to feature collection (Your Index)
  properties:ee.List(['NAME'])
  });
  var date = img.date().format()
  // var LAT = 
  var ft3 = ft2.map(function(f){return f.set("date", date)})
  return inift.merge(ft3)
}

// Iterate Image and save result to feature collection
var newft = ee.FeatureCollection(dataset.iterate(fill, ft))
// Export data
Export.table.toDrive({
  collection: newft,
  description: 'LAI_get',
  fileFormat: 'CSV'
});

仅供自己学习使用，参考文章https://blog.csdn.net/sqwev/article/details/127623543
修改后在csv文件中添加站点名称：

// Change Date here
var dataset = ee.ImageCollection('MODIS/006/MCD43A4')
                .select('Nadir_Reflectance_Band1')                  
                .filter(ee.Filter.date('2007-01-01', '2007-12-31'));
                  
// Upload your point shapefile(Need Change!!)
//var sample_point = ee.FeatureCollection("users/luoyi/FIUX");
var sample_point = ee.FeatureCollection("users/luoyi/FIUX").map(function(feature) {
  return feature.set("NAME", feature.get("Site_ID"));
});

// Start collect data
var ft = ee.FeatureCollection(ee.List([]))
var fill = function(img, ini) {
  // tran ini 2 FeatureCollection
  var inift = ee.FeatureCollection(ini)
  // Extract pixel value
  var ft2 = img.sampleRegions({
  collection:sample_point,
  // Add "NAME" to feature collection (Your Index)
  properties:ee.List(['NAME'])
  });
  var date = img.date().format()
  // var LAT = 
  var ft3 = ft2.map(function(f){return f.set("date", date)})
  return inift.merge(ft3)
}

// Iterate Image and save result to feature collection
var newft = ee.FeatureCollection(dataset.iterate(fill, ft))
// Export data
Export.table.toDrive({
  collection: newft,
  description: 'FLUXBand1_get',
  fileFormat: 'CSV'
});