基于python的MODIS年内NDVI和NDWI值时序制作可视化

本文介绍了如何使用Python处理从Google Earth Engine(GEE)获取的2021年MODIS农田NDVI数据,通过数据整理和滤波,创建年内时序图,以提升可视化质量。最终展示结果并邀请读者交流。
摘要由CSDN通过智能技术生成

这个帖子是我在做一个水稻用地提取到的NDVI年内时序值做的一个图,思路是现在GEE提取到MODIS2021年某农田NDVI值,MODIS每八天一景,全年共覆盖45景影像,将GEE提取到的数据下载下来,再有python进行数据整理,滤波,已得到更好的可视化效果。具体代码如下:

from scipy.signal import savgol_filter #滤波库-savgol_filter函数是Savitzky-Golay滤波器的一种实现
import matplotlib.pyplot as plt # Python中常用的绘图库
plt.rcParams['font.family'] = ['SimHei', 'DejaVu Sans']# 引用黑体和。。

dates = ['01-01', '01-09', '01-17', '01-25', '02-02', '02-10', '02-18', '02-26', '03-06', '03-14', '03-22', '03-30', '04-07', '04-15', '04-23', '05-01', '05-09', '05-17', '05-25', '06-02', '06-10', '06-18', '06-26', '07-04', '07-12', '07-20', '07-28', '08-05', '08-13', '08-21', '08-29', '09-0
gee python 是 Google Earth Engine 提供的 Python API,可以方便地使用 Google Earth Engine 平台中的遥感数据和功能。下面是一种在 gee python 中按年逐月提取 NDVI 的方法: ```python # 引入必要的库 import ee import datetime # 初始化 ee python API ee.Initialize() # 定义一个区域(例如一个地理范围的边界) region = ee.Geometry.Rectangle([xmin, ymin, xmax, ymax]) # 定义起始年份和结束年份 start_year = 2010 end_year = 2020 # 定义月份列表 months = range(1, 13) # 从1到12,表示12个月 # 遍历年份和月份,提取 NDVI 数据 for year in range(start_year, end_year + 1): for month in months: # 将年份和月份转换为日期对象 date = datetime.date(year, month, 1) # 构建一个 ee.Date 对象 ee_date = ee.Date(date.isoformat()) # 构建一个 ee.Image 对象,表示指定日期的影像 image = ee.ImageCollection("MODIS/006/MOD13Q1") \ .filterBounds(region) \ .filterDate(ee_date, ee_date.advance(1, 'month')) \ .select('NDVI') \ .mean() # 每个月的 NDVI 取平均 # 将影像导出为 GeoTIFF 格式 task = ee.batch.Export.image.toDrive(image=image, description='NDVI_' + str(year) + '_' + str(month), folder='export_folder', fileNamePrefix='NDVI_' + str(year) + '_' + str(month), region=region, scale=1000) task.start() print('任务已创建,请前往 Google Drive 中的指定文件夹查看导出的数据。') ``` 这段代码通过遍历年份和月份,针对每个月提取给定区域内的 MODIS NDVI 数据,并将结果导出为 GeoTIFF 格式。你需要将 `xmin`、`ymin`、`xmax` 和 `ymax` 替换为你感兴趣的地理范围的具体边界,将导出保存的文件夹路径替换为你希望保存数据的目录。数据导出的链接将通过电子邮件发送给你。 以上是一种在 gee python 中按年逐月提取 NDVI 数据的方法,希望对你有帮助!
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值