ENH: Storing RPCs in xarray

[remi-braun]

Hi there,

Currently, the RPCs are not stored into xarray (whereas the GCPs since #376). Moreover, rasterio offers the possibility to retrieve them.

So even if I have opened the xarray, I need to add this annoying bit of code to retrieve them:

with rasterio.open("my/raster.tif") as ds:
    rpcs = ds.rpcs

Moreover, if I'm right, you are handling them in reprojection, so they should - in my opinion - be retrievable without calling rasterio.

Is there any specific reason to not storing them? Or maybe I've overlooked something?

[snowman2]

Just hasn't been implemented. A MR is welcome.

[remi-braun]

Thanks for the quick answer. I don't really know what it implies in rioxarray, as I dont fully understand all the whereabouts of the current GCP loading. Should it be done with the same pattern?

[snowman2]

It likely should be done following a similar pattern.