Reproject indexing errors #762
Replies: 3 comments 2 replies
-
If I had a file, this is what I would try: dataset[["rotated_latitude_longitude", "grid_latitude", "grid_longitude", "tasmax"]].isel(
ensemble_member=0,
time=0,
).rio.set_spatial_dims(
x_dim="grid_longitude",
y_dim="grid_latitude",
inplace=True,
).rio.reproject("'EPSG:27700") Steps:
|
Beta Was this translation helpful? Give feedback.
-
@snowman2 looks like this data is using location arrays for geo referencing. Is this something Would be nice to agree on a format for dealing with location arrays. |
Beta Was this translation helpful? Give feedback.
-
Thanks again @Kirill888 and @snowman2. FYI: looking another route: https://gis.stackexchange.com/questions/312084/transform-coordinate-from-rotated-lat-lon-to-normal-lat-lon-wgs84 |
Beta Was this translation helpful? Give feedback.
-
Thanks for the library. I'm attempting to process https://catalogue.ceda.ac.uk/uuid/d5822183143c4011a2bb304ee7c0baf7, reprojected to
EPSG:27700
. The (I think) related examples I've found thus far include:My guess is I need to figure the exact coordinate structure to solve the
undefined
components, and I wonder if it relates toEURO-CORDEX
and/orcartopy
RotatedGeodetic
. While I try to figure that (at which point I assume I could then usedataset.rio.write_crs()
, this is how far I've got via other routes:I've also tried using
set_spatial_dims
but that eventually leads to the same1-dimensional
error.Update
I think this relates to #772. I've ended up going between
rioxarray
andosgeo.gdal
across various steps, hope to post below if I reach a solution.Beta Was this translation helpful? Give feedback.
All reactions