hydro_opendata.stac.mini模块¶
用于从minio中获取数据集范围
数据准备¶
In [1]:
Copied!
from hydro_opendata.catalog.minio import ERA5LCatalog, GPMCatalog, GFSCatalog
import geopandas as gpd
from hydro_opendata.catalog.minio import ERA5LCatalog, GPMCatalog, GFSCatalog
import geopandas as gpd
In [2]:
Copied!
aoi = gpd.read_file('basin.geojson')
aoi
aoi = gpd.read_file('basin.geojson')
aoi
Out[2]:
| id | geometry | |
|---|---|---|
| 0 | 0 | MULTIPOLYGON (((122.44241 39.80139, 122.39342 ... |
获取era5数据清单¶
In [3]:
Copied!
era5 = ERA5LCatalog()
era5.datasets
era5 = ERA5LCatalog()
era5.datasets
Out[3]:
{'wis': {'start_time': numpy.datetime64('2015-07-01T00:00:00'),
'end_time': numpy.datetime64('2021-12-31T23:00:00'),
'bbox': [115, 38, 136, 54]}}
In [4]:
Copied!
e = era5.search(aoi=aoi)
e
e = era5.search(aoi=aoi)
e
Out[4]:
| id | dataset | start_time | end_time | geometry | |
|---|---|---|---|---|---|
| 0 | era5-land | wis | 2015-07-01T00:00:00 | 2021-12-31T23:00:00 | POLYGON ((122.39342 39.81027, 122.31744 39.833... |
获取gpm数据清单¶
In [5]:
Copied!
gpm = GPMCatalog()
gpm.datasets
gpm = GPMCatalog()
gpm.datasets
Out[5]:
{'wis': [{'time_resolution': '30 minutes',
'start_time': numpy.datetime64('2016-01-01T00:00:00.000000000'),
'end_time': numpy.datetime64('2023-10-11T08:30:00.000000000'),
'bbox': [73, 3, 136, 54]},
{'time_resolution': '1 day',
'start_time': numpy.datetime64('2000-06-01T00:00:00.000000000'),
'end_time': numpy.datetime64('2023-09-26T23:59:59.000000000'),
'bbox': [73, 3, 136, 54]}],
'camels': [{'time_resolution': '30 minutes',
'start_time': numpy.datetime64('2022-01-01T00:00:00.000000000'),
'end_time': numpy.datetime64('2023-08-31T23:30:00.000000000'),
'bbox': [-125, 25, -66, 50]},
{'time_resolution': '1 day',
'start_time': numpy.datetime64('2000-06-01T00:00:00.000000000'),
'end_time': numpy.datetime64('2014-12-31T23:59:59.000000000'),
'bbox': [-125, 25, -66, 50]}]}
In [6]:
Copied!
g = gpm.search(aoi=aoi)
g
g = gpm.search(aoi=aoi)
g
Out[6]:
| id | dataset | time_resolution | start_time | end_time | geometry | |
|---|---|---|---|---|---|---|
| 0 | gpm-imerg-early | wis | 30 minutes | 2016-01-01T00:00:00.000000000 | 2023-10-11T08:30:00.000000000 | POLYGON ((122.39342 39.81027, 122.31744 39.833... |
| 0 | gpm-imerg-early | wis | 1 day | 2000-06-01T00:00:00.000000000 | 2023-09-26T23:59:59.000000000 | POLYGON ((122.39342 39.81027, 122.31744 39.833... |
获取gfs数据清单¶
In [7]:
Copied!
gfs = GFSCatalog('tp')
gfs.datasets
gfs = GFSCatalog('tp')
gfs.datasets
Out[7]:
{'wis': [{'start': '2016-07-10T00',
'end': '2022-08-31T18',
'bbox': [115, 38, 136, 54]},
{'start': '2022-09-01T00',
'end': '2023-10-11T06',
'bbox': [73, 3, 136, 54]}]}
In [8]:
Copied!
f = gfs.search(aoi=aoi)
f
f = gfs.search(aoi=aoi)
f
Out[8]:
| id | dataset | start_time | end_time | geometry | |
|---|---|---|---|---|---|
| 0 | gfs_atmos.tp | wis | 2016-07-10T00 | 2022-08-31T18 | POLYGON ((122.39342 39.81027, 122.31744 39.833... |
| 0 | gfs_atmos.tp | wis | 2022-09-01T00 | 2023-10-11T06 | POLYGON ((122.39342 39.81027, 122.31744 39.833... |
Last update:
2023-10-22