Dear colleagues,

I am working on fusing MODIS NDVI with Landsat NDVI for the same region and date. I have already reprojected the MODIS rasters using .rio.reproject_match() from rioxarray to perfectly match the Landsat raster in:

• CRS (EPSG code matches exactly),
• Pixel resolution,
• Spatial extent and shape,
• Coordinate grids (x and y arrays).

After reprojecting, the MODIS metadata (transform, bounds, resolution, shape) becomes identical to the Landsat raster. However, when I overlay the two, MODIS NDVI still appears shifted northward, even though everything in the metadata says they are aligned.

Check:

Landsat bounds: (189645.0, 3328725.0, 229275.0, 3368685.0) MODIS (shifted) bounds: (189645.0, 3328725.0, 229275.0, 3368685.0)