基于 GEE 提取指定点的土壤湿度数据

最新推荐文章于 2025-05-30 20:09:15 发布

@HNUSTer

最新推荐文章于 2025-05-30 20:09:15 发布

阅读量355

点赞数 9

分类专栏： Google Earth Engine（GEE）文章标签： GEE 云计算云平台遥感大数据数据集

本文链接：https://blog.csdn.net/ZK180531/article/details/145578358

版权

Google Earth Engine（GEE）专栏收录该内容

87 篇文章

订阅专栏

1 数据介绍

2 完整代码

3 运行结果

1 数据介绍

数据集ee.ImageCollection("NASA/GLDAS/V021/NOAH/G025/T3H")

数据集NASA/GLDAS/V021/NOAH/G025/T3H，包含很多数据集，包括气象、土壤、气温等数据，例如：CanopInt_inst植物冠层地表水含量、ECanop_tavg树冠水分蒸发、ESoil_tavg裸土蒸发量、Snowf_tavg雪速率、SoilMoi0_10cm_inst土壤湿度、SoilTMP0_10cm_inst土壤温度、Tair_f_inst气温、Wind_f_inst风速等等。

2 完整代码

var roi = /* color: #d63000 */ee.Geometry.Point([121.3664979892787, 37.45796001127996]);
// Create a chart for SoilMoi100_200cm_inst trend in point 1.
Map.centerObject(roi,7)
var dataset = ee.ImageCollection('NASA/GLDAS/V021/NOAH/G025/T3H')
          .filterBounds(roi)
          .filterDate('2006-01-01', '2016-12-31')
          .select('SoilMoi100_200cm_inst');

var radius_lst = ee.List([250, 500, 1000, 1500, 2000]);

var colByRadius = ee.List(radius_lst).map(function (radius) {

  var months = ee.List.sequence(1, 12);
  var start_year = 2006;
  var start_date = '2006-01-01';
  var end_year = 2016;
  var end_date = '2016-12-31';

  var years = ee.List.sequence( start_year, end_year);

  var ptsbuff = ee.FeatureCollection(roi).map(function (radius, bounds) {
    return function(pt) {
      pt = ee.Feature(pt);
    return bounds ? pt.buffer(radius).bounds() : pt.buffer(radius);
  };
  }(radius, false));

  var byMonthYear =  years.map(function (y) {
      return months.map(function(m) {
        var w = dataset.filter(ee.Filter.calendarRange(y, y, 'year'))
                    .filter(ee.Filter.calendarRange(m, m, 'month'))
                    .mean();
        var pointMean = w.reduceRegion(
          {reducer:ee.Reducer.first(), 
          geometry:ptsbuff, scale:500});  
     
      // set the dictionary as property and cast to an ee.Image, as setMulti returns an ee.Element.
      return ee.Feature(null).set("value",
                                  pointMean
                                    .get("SoilMoi100_200cm_inst"))
                                    .set("year",y)
                                    .set("month",m)
                                    .set("radius", radius);
      

      });
    }).flatten();

  return byMonthYear;

}).flatten();

Export.table.toDrive({collection:ee.FeatureCollection(colByRadius), description:"csvExport"})