I am am trying to merge raster files from the Sentinel 2 Constellation (after some processing) with rioxarray for the UTM Zones 32 and 33. It's looking good, execpt when an larger area for example the City of Berlin(Germany) in included in Scences in UTM32 (EPSG:32632) and UTM33(EPSG:32633). Before merging the images I project them to UTM33. After the mergeing Everything is getting to look unsharp, like having the street twice, but with a quite a distance away.
I am merging those with the given code:
def merge_min_intensity(folder: Path, save_file_name: Path, epsg_code: int) -> xarray.DataArray: rasters: List[xarray.DataArray] =  for file in folder.iterdir(): if str(file).endswith(".tif"): raster = rxa.open_rasterio(file) if int(raster.rio.crs.data['init'].split(":")) != epsg_code: raster = raster.rio.reproject(epsg_code) rasters.append(raster) ret_raster = merge_arrays(rasters, method="min") ret_raster.rio.to_raster(save_file_name) return ret_raster
What might be the reason for this?