Alan D. Snow

Do you think that their suggestion is worth documenting?

Looks like the data is shifted one cell over.

> However, once a CRSIndex is created, should spatial_ref still be preserved or not? - One of factors in the design of `rioxarray` is to be able to call `to_netcdf`...

Potential alternative: #3239

Yes, `to_epsg` is a very costly call. Removing that or placing it as the fallback instead of happening first would definitely speed things up.

Related: https://github.com/opendatacube/datacube-core/issues/1321

https://github.com/pyproj4/pyproj/discussions/1454#discussioncomment-12160165 After some debugging, it appears that it doesn't matter if the `files.geojson` file is there. You just need to start a new python session for it to work as...

It appears that the context caches the files looked in the context [here](https://github.com/OSGeo/PROJ/blob/77bb0cc0d2d80d14bcfc7921647fa39b144fe0ca/src/proj_internal.h#L813). It gets reset if you set the search paths [code](https://github.com/OSGeo/PROJ/blob/77bb0cc0d2d80d14bcfc7921647fa39b144fe0ca/src/ctx.cpp#L132). I verified that if you call `set_data_dir`...

@rouault, what are your thoughts on this behavior? Do you think that it would be helpful to add a method for the user to invalidate the cache on the context...

Mind debugging internal PROJ? See: https://pyproj4.github.io/pyproj/stable/advanced_examples.html#debugging-internal-proj