Skip to content

Reader模块

读取数据。

功能

通过给定的AOI范围： - 读取MinIO服务器上的数据，如era5、gpm、gfs等； - 通过STAC读取landsat、sentinel等遥感数据。

使用

  • 读取era5数据
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
from hydro_opendata.reader.minio import ERA5LReader
import numpy as np

era5 = ERA5LReader()

# 指定开始及结束时间
start_time=np.datetime64("2021-06-01T00:00:00.000000000")
end_time=np.datetime64("2021-06-30T23:00:00.000000000")

# 通过指定四至范围读取
bbox=(121,39,123,40)
ds1 = era5.open_dataset(data_variables=['Total precipitation'], start_time=start_time, end_time=end_time, dataset='wis', bbox=bbox)

# 通过矢量数据文件读取
shp = 'basin.shp'
ds2 = era5.from_shp(data_variables=['Total precipitation'], start_time=start_time, end_time=end_time, dataset='wis', shp=shp)

# 通过已有aoi对象读取
aoi = gpd.read_file(shp)
ds3 = era5.from_aoi(data_variables=['Total precipitation'], start_time=start_time, end_time=end_time, dataset='wis', aoi=aoi)
  • 读取gpm数据
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
from hydro_opendata.reader.minio import GPMReader
import numpy as np

# 指定数据集wis或camels
gpm = GPMReader('wis')

# 指定开始及结束时间
start_time=np.datetime64("2023-06-01T00:00:00.000000000")
end_time=np.datetime64("2023-06-30T23:30:00.000000000")

# 通过指定四至范围读取
bbox=(121,39,123,40)
ds1 = gpm.open_dataset(start_time=start_time, end_time=end_time, dataset='wis', bbox=bbox, time_resolution='1d')

# 通过矢量数据文件读取
shp = 'basin.shp'
ds2 = gpm.from_shp(start_time=start_time, end_time=end_time, dataset='wis', shp=shp, time_resolution='1d')

# 通过已有aoi对象读取
aoi = gpd.read_file(shp)
ds3 = gpm.from_aoi(start_time=start_time, end_time=end_time, dataset='wis', aoi=aoi, time_resolution='1d')
  • 读取gpm数据
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
from hydro_opendata.reader.minio import GFSReader
import numpy as np

gfs = GFSReader()

# 指定日期
creation_date=np.datetime64("2023-06-01")

# 通过指定四至范围读取
bbox=(121,39,123,40)
ds1 = gfs.open_dataset(creation_date=creation_date, creation_time='00', dataset='wis', bbox=bbox)

# 通过矢量数据文件读取
shp = 'basin.shp'
ds2 = gfs.from_shp(creation_date=creation_date, creation_time='00', dataset='wis', shp=shp)

# 通过已有aoi对象读取
aoi = gpd.read_file(shp)
ds3 = gfs.from_aoi(creation_date=creation_date, creation_time='00', dataset='wis', aoi=aoi)
Last update: 2023-11-02