Add bounds fails with coordinates crossing the antimeridian or the Greenwich meridian.
@aulemahal y travaille déjà.
Depuis le 14/12/2022.
Évaluation
Cette issue n'a pas encore été évaluée.
Description
Yet another bug with .cf.add_bounds... This one was already present before the 2D addition, but it would have been quite infrequent with 1D coordinates.
The interpolation / extrapolation estimates are not aware of the cyclic nature of the longitude.
MWE:
import xarray as xr
import cf_xarray
ds = xr.Dataset(
{'data': (('x', 'y'), np.ones((3, 3)))},
coords={
'x': [1, 2, 3],
'y': [1, 2, 3],
'lon': (('x',), [170, -170, -150], {'units': 'degrees_east'}),
'lat': (('y',), [-10, 0, 10], {'units': 'degrees_north'}),
},
)
ds.cf.add_bounds(('lon', 'lat')).lon_bounds
yields:
array([[ 340., 0.],
[ 0., -340.],
[-160., -140.]])
but we would expect:
array([[ 160., 0.],
[ 0., -160.],
[-160., -140.]])
Sadly, a quick internet search hasn't given me any pointers towards tools that already take this into account, except for projection-based tools (like pyproj).
We could have a check for "longitude" that looks for jumps and wraps the coordinates, but that jump threshold would be arbitrary. It's nonetheless the simplest way out I see, before parsing grid_mapping into pyproj objects.
- Langage dominant
- Python
- Étoiles
- 183
- Forks
- 52
- Merge moyen
- 6 h 39 min
- PR mergées (30 j)
- 1
Préparer son environnement
Ce projet ne fournit ni conteneur de développement, ni Dockerfile, ni guide de contribution : l'installation est à votre charge. Commencez par son README, et consultez notre guide de la première contribution pour les étapes générales.
Par où commencer
- Lisez l'issue en entier, puis le guide de contribution du projet.
- Signalez en commentaire que vous la prenez — cela évite que deux personnes fassent le même travail.
- Forkez le dépôt et travaillez sur une branche.
- Ouvrez une pull request qui référence le numéro de l'issue.
Autres issues de xarray-contrib/cf-xarray
-
Difficulté 2/5 1-3 heures Accessibilité débutants 68/100
xarray-contrib/cf-xarray#623 ·
-
Consider adding opt-in support for circular axes when calling `bounds_to_vertices` (related to #594)Ouverte
Difficulté 5/5 Plus d'une semaine Accessibilité débutants 38/100
xarray-contrib/cf-xarray#596 · 1 commentaire ·
-
Difficulté 3/5 1-2 jours Accessibilité débutants 25/100
xarray-contrib/cf-xarray#570 ·
-
Difficulté 4/5 3-5 jours Accessibilité débutants 15/100
xarray-contrib/cf-xarray#568 ·
-
Difficulté 3/5 1-2 jours Accessibilité débutants 35/100
xarray-contrib/cf-xarray#566 · 2 commentaires ·
Toutes les issues de xarray-contrib/cf-xarray
Issues similaires
-
changelog investigate
Difficulté 2/5 1-3 heures Accessibilité débutants 65/100
ramnes/notion-sdk-py#423 ·
-
Difficulté 2/5 1-3 heures Accessibilité débutants 77/100
Les mainteneurs répondent en général sous 1 jour
-
Difficulté 2/5 1-3 heures Accessibilité débutants 78/100
Les mainteneurs répondent en général sous 1 jour
-
Difficulté 1/5 1-3 heures Accessibilité débutants 89/100
duckdb/duckdb-web#7514 ·
Les mainteneurs répondent en général sous 1 jour
-
bug
Difficulté 2/5 1-3 heures Accessibilité débutants 72/100
jbacalso24/mimry#525 ·
Les mainteneurs répondent en général sous 1 jour