如何解决打开具有许多组的远程 zarr 商店并使用 xarray 保持坐标
我想读入 https://hrrrzarr.s3.amazonaws.com/index.html#sfc/20210208/20210208_00z_anl.zarr/ 的远程 zarr 存储。 zarr 商店的信息位于 https://mesowest.utah.edu/html/hrrr/zarr_documentation/zarrFileVariables.html
我能够读入一个变量,但它似乎没有捕获与该变量相关的坐标或属性(我很可能缺少 open_mfdataset
或 open_zarr
的 kwargs)。因为有不同级别的嵌套,我不确定什么是正确的传递路径
import xarray as xr
import s3fs
fs = s3fs.S3FileSystem(anon=True)
uri = "s3://hrrrzarr/sfc/20210208/20210208_00z_anl.zarr/10m_above_ground/UGrd/10m_above_ground"
file = s3fs.S3Map(uri,s3=fs)
ds = xr.open_mfdataset([file],engine="zarr")
>>> ds
<xarray.Dataset>
Dimensions: (projection_x_coordinate: 1799,projection_y_coordinate: 1059)
Dimensions without coordinates: projection_x_coordinate,projection_y_coordinate
Data variables:
UGrd (projection_y_coordinate,projection_x_coordinate) float16 dask.array<chunksize=(150,150),Meta=np.ndarray>
uri = "s3://hrrrzarr/sfc/20210208/20210208_00z_anl.zarr/10m_above_ground/UGrd"
file = s3fs.S3Map(uri,engine="zarr")
>>> ds
<xarray.Dataset>
Dimensions: (projection_x_coordinate: 1799,projection_y_coordinate: 1059)
Coordinates:
* projection_x_coordinate (projection_x_coordinate) float64 -2.698e+06 ......
* projection_y_coordinate (projection_y_coordinate) float64 -1.587e+06 ......
Data variables:
forecast_period timedelta64[ns] ...
forecast_reference_time datetime64[ns] ...
height float64 ...
pressure float64 ...
time datetime64[ns] ...
解决方法
Xarray 无法理解嵌套的 zarr 组。它期望所有变量和坐标都在一个平面组中。我认为您唯一的选择是手动合并数据集。你试过吗
ds = xr.open_mfdataset([file1,file2],engine="zarr")
?
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。