Я пытаюсь записать точки в netcdf и поэтому создаю xr.DataArray.
Если бы я хотел написать 2d-массив, я бы просто сделал:
import xarray as xr
da = xr.DataArray(data=tmp, dims=["lon", "lat", "time"], coords={'lon' : lon2d, 'lat' : lat2d, 'time': np.array([np.datetime64(day)])})
Для случая точек я пытаюсь:
da = xr.DataArray(data=tmp, dims=["z", "time"], coords={'time': np.array([np.datetime64(day)])})
Однако я не знаю, как дать координаты широты / долготы 'z', как в случае с массивом с накоплением, как в примере ниже:
<xarray.DataArray 'P' (time: 30, z: 3)>
dask.array<reshape, shape=(30, 3), dtype=float64, chunksize=(30, 3), chunktype=numpy.ndarray>
Coordinates:
* time (time) datetime64[ns] 2017-01-01 2017-01-02 ... 2017-01-30
lat (z) float64 10.12 11.12 12.12
lon (z) float64 47.62 46.62 45.62
Dimensions without coordinates: z
Attributes:
units: mm.day-1
long_name: Precipitation
standard_name: Precipitation