GEE批量下载逐月MOD13Q1 NDVI数据代码
时间: 2023-10-08 08:04:27 浏览: 248
以下是使用Python和GEE API批量下载逐月MOD13Q1 NDVI数据的示例代码:
```
import ee
import datetime
# 初始化 GEE
ee.Initialize()
# 定义区域范围
region = ee.Geometry.Polygon(
[[[-124.8876953125, 42.014192020743394],
[-120.64697265625, 42.014192020743394],
[-120.64697265625, 45.72017789701005],
[-124.8876953125, 45.72017789701005],
[-124.8876953125, 42.014192020743394]]], None, False)
# 定义时间范围
start_date = datetime.datetime(2000, 1, 1)
end_date = datetime.datetime(2021, 12, 31)
# 定义数据集
dataset = ee.ImageCollection('MODIS/006/MOD13Q1').select('NDVI')
# 遍历每个月份下载数据
for year in range(start_date.year, end_date.year + 1):
for month in range(1, 13):
date_str = f'{year}-{month:02}-01'
current_date = datetime.datetime.strptime(date_str, '%Y-%m-%d')
if current_date < start_date or current_date > end_date:
continue
start_date_str = current_date.strftime('%Y-%m-%d')
end_date_str = (current_date + datetime.timedelta(days=32)).strftime('%Y-%m-%d')
filtered = dataset.filterDate(start_date_str, end_date_str).filterBounds(region)
image = filtered.mean()
task = ee.batch.Export.image.toDrive(image=image,
description=f'MOD13Q1_NDVI_{year}_{month:02}',
folder='GEE MODIS NDVI',
region=region,
scale=500,
maxPixels=1e13)
task.start()
print(f'Starting export task for {year}-{month:02}')
```
此代码将从2000年1月到2021年12月遍历每个月,并使用GEE API从MOD13Q1数据集中筛选出NDVI图像并将其平均化,然后将其导出到Google Drive中的文件夹“GEE MODIS NDVI”中。你可以在代码中修改区域、时间范围、数据集和导出设置,以适应你的需求。
阅读全文