<div dir="ltr">According to rasterio's <a href="https://rasterio.readthedocs.io/en/latest/topics/reproject.html">documentation</a> "reproject" is a geospatial-specific analog to SciPy’s <a href="https://docs.scipy.org/doc/scipy/reference/generated/scipy.ndimage.geometric_transform.html#scipy.ndimage.geometric_transform">geometric_transform</a>.<br><div>But the thing is "geometric_transform" doesn't support masked arrays so it can't properly transform images with nodata value.<br></div><div>My question: is it possible to write a custom transformer function and use it in GDAL to transform an image (considering nodata) to overcome the lack of the appropriate functionality in SciPy.</div><div>Thanks!</div></div>