`reproject` with `src_geoloc_array` introduces unexpected shift for group of pixels about rasterio HOT 2 CLOSED

Perhaps useful for context, this is the effect on a real-world raster, where the roads become significantly misaligned:

from rasterio.

This bug seems to be caused by rasterio choosing to use an Approximate Transformer - the bUseApproxTransformer flag is always True here unless using RPCs:

With the current main branch, with cdef bint bUseApproxTransformer = True:

And changing it to False

My naive suggestion for a fix would be adding bUseApproxTransformer = False inside the elif src_geoloc_array is not None branch in the _reproject function, although I'm not fully sure of what other side effects using a non-approximate transformer might give (slightly worse performance?):

from rasterio.