如何使用Python读取IMERG数据
概述:
此配方显示了如何使用Python从全球降水测量(GPM)任务的IMERG数据集中读取数据。
例:
**示例数据: GPM Level 3 IMERG每月0.1 x 0.1度降水(GPM_3IMERGM)于2015年7月。
预计完成以下程序的时间:20分钟
先决条件:
任务:查看数据
最佳时间:用户想要使用Python读取GPM IMERG数据
要求: Python和免费软件包:numpy,matplotlib,basemap和h5py。Matplotlib和底图仅用于绘图。
程序:
1.下载数据
在GES DISC访问数据之前,用户必须首先注册Earthdata登录,然后通过以下步骤授权访问GES DISC中的数据: data-access。
在“ 搜索”字段中,输入GPM_3IMERGM并按Enter键。
单击最新版本的GPM_3IMERGM数据,当前版本为5。
注意:此配方适用于任何IMERG数据,而不仅仅是月度估算。
图1:GPM_3IMERGM的 GES DISC搜索结果示例。
2.运行以下Python脚本
# This is a sample script that reads and plots the IMERG data.
from mpl_toolkits.basemap import Basemap, cm
import matplotlib.pyplot as plt
import numpy as np
import h5py as h5py
dataset = h5py.File('/path/to/data/3B-MO.MS.MRG.3IMERG.20150801-S000000