GEE：数据结构——对List的操作(添加、删除、替换、获取长度、遍历、过滤)

### 使用 MCD12Q2 数据集计算植被生长季长度 MCD12Q2 是基于 MODIS 的全球陆地覆盖产品之一，提供了多种分类方案下的年度土地覆盖类型数据。尽管该数据主要用于描述土地覆盖变化，但它可以与其他遥感指标（如 NDVI 或 EVI）结合来估算植被生长季节的长度。 以下是通过 Google Earth Engine (GEE) 实现这一目标的方法： #### 1. 加载必要的库和初始化地图 为了操作地球引擎中的影像集合，首先需要导入 `ee` 库并创建一个基础的地图对象。 ```javascript // 初始化Google Earth Engine Map.setCenter(-98.75, 39.5, 4); // 设置初始视图中心点和缩放级别 var startDate = '2001-01-01'; // 起始日期 var endDate = '2021-12-31'; // 结束日期 ``` #### 2. 导入 MCD12Q2 和辅助数据集 虽然 MCD12Q2 提供的是土地覆盖信息，但通常会结合其他时间序列数据（例如 MODIS NDVI/EVI），以便更精确地定义植被活动周期。 ```javascript // 导入MCD12Q2数据集 var mcd12q2 = ee.ImageCollection("MODIS/006/MCD12Q2"); // 过滤时间和区域范围内的NDVI数据作为补充 var modisNdvi = ee.ImageCollection("MODIS/006/MOD13A2") .filterDate(startDate, endDate); ``` #### 3. 定义植被活跃期检测函数 可以通过分析 NDVI 时间曲线的变化趋势找到每年的生长期起点和终点。这里提供一种简单方法：寻找年内最大值及其对应的上升沿下降沿位置。 ```javascript function calculateGrowingSeason(image){ var ndviBand = image.select('NDVI'); // 找到全年最高NDVI值所在的位置 var maxNdviTime = ndviBand.reduce(ee.Reducer.max()); // 假设春季开始于NDVI快速增加的时间段；秋季结束则对应减少阶段 var springStart = ...; // 自定义逻辑实现此处省略... var autumnEnd = ...; return image.addBands(springStart).addBands(autumnEnd); } ``` #### 4. 映射处理整个时间序列 将上述自定义功能应用于所有年份的数据上，并提取感兴趣的结果字段。 ```javascript var processedImages = modisNdvi.map(calculateGrowingSeason); // 将结果导出或者绘制出来查看效果 print(processedImages.limit(5)); Map.addLayer(processedImages.first(), {min:0,max:9}, "Sample Layer"); ``` 以上代码片段仅展示了一个基本框架[^3]，实际应用中可能还需要考虑更多细节因素比如异常天气影响、不同生态系统特性差异等等。 --- ###