Bug when using `bounds` in `merge_arrays`

[OnnoEbbens]

I use the rioxarray.merge.merge_arrays function to merge 4 tiles into a single DataArray. When I use the bounds argument I get an incorrectly merged DataArray. Please sea code and figures below.

tiffs.zip

Code Sample

import rioxarray as rxr
from rioxarray.merge import merge_arrays
import matplotlib.pyplot as plt
bounds = [117607.88 , 498001.379, 123698.88 , 502618.78]
das = []
for i in range(4):
    das.append(rxr.open_rasterio(f'da{i}.tiff'))
out = merge_arrays(das)
plt.figure()
out.plot()
plt.title('merged array without bounds')

out_with_bounds = merge_arrays(das, bounds=bounds)
plt.figure()
out_with_bounds.plot()
plt.title('merged array using bounds')

Environment Information

rioxarray (0.19.0) deps: rasterio: 1.4.3 xarray: 2025.3.1 GDAL: 3.9.3 GEOS: 3.13.0 PROJ: 9.5.0 PROJ DATA: c:\Users\oebbe\anaconda3\envs\nlmod2025\Lib\site-packages\rasterio\proj_data GDAL DATA: C:\Anaconda3\Library\share\gdal

Other python deps: scipy: 1.15.2 pyproj: 3.7.1

System: python: 3.12.9 | packaged by Anaconda, Inc. | (main, Feb 6 2025, 18:49:16) [MSC v.1929 64 bit (AMD64)] executable: c:\Users\oebbe\anaconda3\envs\nlmod2025\python.exe machine: Windows-11-10.0.22631-SP0