GEE 下载MODIS/061/MOD16A2GF

dNBR值（differenced Normalized Burn Ratio）是衡量火灾烈度的指标之一，可以通过MODIS/061/MCD64A1数据计算得到。 以下是基于GEE计算dNBR值的步骤： 1. 打开GEE平台并登录账号。 2. 在左侧菜单栏中选择“地图”选项，进入地图界面。 3. 在地图界面中，点击左上角的“+”按钮，添加MODIS/061/MCD64A1数据。 4. 在左侧图层栏中选择“MODIS/061/MCD64A1_500m_aid0001”图层，并调整时间范围至所需的时间段。 5. 在左侧菜单栏中选择“代码编辑器”，并在编辑器中输入以下代码： ``` // 定义AOI区域 var aoi = ee.Geometry.Polygon( [[[113.0, 33.0], [115.0, 33.0], [115.0, 35.0], [113.0, 35.0], [113.0, 33.0]]], null, false); // 选择计算的时间段 var start_date = ee.Date.fromYMD(2019, 1, 1); var end_date = ee.Date.fromYMD(2019, 12, 31); // 选择计算的MODIS波段 var modis_bands = ['BurnDate','QA']; // 筛选出指定时间段内的MODIS数据 var modis_dataset = ee.ImageCollection('MODIS/006/MCD64A1') .select(modis_bands) .filterDate(start_date, end_date) .filterBounds(aoi); // 定义NDVI计算公式 var nbr_formula = function(image) { var nbr = image.normalizedDifference(['sur_refl_b04', 'sur_refl_b07']); return nbr.rename('NBR'); }; // 计算dNBR值 var pre_fire = ee.Image(modis_dataset.sort('system:time_start', false).first()).select('BurnDate'); var post_fire = ee.Image(modis_dataset.sort('system:time_start', true).first()).select('BurnDate'); var pre_nbr = nbr_formula(pre_fire); var post_nbr = nbr_formula(post_fire); var dNBR = post_nbr.subtract(pre_nbr).rename('dNBR'); // 将结果可视化 Map.centerObject(aoi, 9); Map.addLayer(dNBR, {min:-1,max:1,palette:['FF0000','FFFF00','00FF00']}, 'dNBR'); ``` 上述代码中，首先定义了计算区域（AOI）、计算时间段（start_date和end_date）以及需要计算的MODIS波段（modis_bands）。接着，通过ImageCollection函数筛选出指定时间段内的MODIS数据，并选择需要的波段。定义了计算dNBR值所需的函数（nbr_formula），并使用该函数计算出前火情（pre_fire）和后火情（post_fire）的NBR值，最终得出dNBR值（dNBR）。最后将结果可视化，其中min、max和palette参数可根据需要进行调整。 完成上述代码的输入和运行后，即可在地图上看到计算得到的dNBR值。

<think>好的，我需要回答用户关于如何在Google Earth Engine中使用MODIS/061/MOD09A1数据集计算NDVI的问题。首先，我应该回忆一下NDVI的计算公式，即$(NIR - Red)/(NIR + Red)$。接下来，我需要确认MOD09A1数据集中的波段名称是否与常规的NIR和Red对应。根据MODIS的文档，MOD09A1的波段3是红光（Red，620-670nm），波段2是近红外（NIR，841-876nm），所以正确的波段应该是这两个。 然后，我需要考虑如何在GEE中加载数据集。用户提到的是MODIS/061/MOD09A1，所以代码中应该正确引用这个数据集。接着，要选择时间范围，可能需要设置一个日期过滤器，比如最近一年。然后，应用映射函数对每个影像计算NDVI，并添加到影像集合中。 用户可能还希望可视化结果，所以需要设置可视化参数，比如min和max值，以及调色板。最后，将结果添加到地图上显示，并导出到Drive或者生成图表。另外，根据参考的引用[2]，在循环中添加Map.addLayer可能会有帮助，但这里可能只需要在计算后添加图层即可。 需要注意代码的正确性，比如波段名称是否正确，日期格式是否正确，以及NDVI计算后的重命名是否正确。还要确保用户知道如何调整时间范围和导出结果。如果有相关的应用场景或优化方法，可以作为相关问题提供，比如如何处理云覆盖或时间序列分析。</think>在Google Earth Engine中使用MODIS/061/MOD09A1数据集计算NDVI的步骤如下： ### 1. 加载数据集与波段选择 MOD09A1数据集的地表反射率波段对应关系为： - 近红外（NIR）：波段2（841-876nm） - 红光（Red）：波段3（620-670nm） NDVI计算公式为： $$ \text{NDVI} = \frac{\text{NIR} - \text{Red}}{\text{NIR} + \text{Red}} $$ ### 2. 示例代码 ```javascript // 加载MOD09A1数据集 var modis = ee.ImageCollection('MODIS/061/MOD09A1') .filterDate('2022-01-01', '2022-12-31'); // 定义NDVI计算函数 var addNDVI = function(image) { var ndvi = image.normalizedDifference(['sur_refl_b02', 'sur_refl_b01']) .rename('NDVI'); return image.addBands(ndvi); }; // 应用NDVI计算 var withNDVI = modis.map(addNDVI); // 设置可视化参数 var visParams = { min: -0.2, max: 1, palette: ['FFFFFF', 'CE7E45', 'DF923D', 'F1B555', 'FCD163', '99B718', '74A901', '66A000', '529400', '3E8601', '207401', '056201', '004C00', '023B01', '012E01', '011D01', '011301'] }; // 显示结果 Map.addLayer(withNDVI.mean(), visParams, 'NDVI Composite'); // 导出结果到Google Drive（可选） Export.image.toDrive({ image: withNDVI.mean(), description: 'MOD09A1_NDVI_2022', scale: 500, region: geometry // 需要先定义感兴趣区域 }); ``` ### 3. 关键点说明 - 波段名称注意：原始代码中波段编号可能需要调整，根据GEE最新数据集文档确认[^1] - 时间序列处理：使用`.filterDate()`限定时间范围后，通过`.map()`批量处理影像 - 质量控制：实际应用中需结合QA波段进行数据过滤[^2] - 空间分辨率：MOD09A1原始分辨率为500米 ### 4. 结果验证 可通过控制台打印元数据验证计算： ```javascript print(withNDVI.first()); ```