After merging this #781 (initial shift for N&K) and moving all operations needed in base.py, a new problem emerges in case of CoregPipeline.
Indeed, we can use reference_elev = reference_elev.translate(-shift_x, -shift_y) as reference_elev was cast from DEM to array in _preprocess_coreg_fit_raster_raster.
@rhugonnet I am pointing you to the new ticket; if you could explain it to me in more detail the trick I need to use, that would be great :)
Thank you