Example - Merge

import rioxarray # for the extension to load
import xarray
from rioxarray.merge import merge_arrays
# Note: You can merge datasets with the merge_datasets method

%matplotlib inline

Load in xarray dataset

API reference:

• rioxarray.open_rasterio

• rioxarray.merge.merge_arrays

• rioxarray.merge.merge_datasets

dem_test = "../../test/test_data/input/MODIS_ARRAY.nc"
rds = rioxarray.open_rasterio(dem_test)
arrays = [
    rds.isel(x=slice(100), y=slice(100)),
    rds.isel(x=slice(100, 200), y=slice(100, 200)),
    rds.isel(x=slice(100), y=slice(100, 200)),
    rds.isel(x=slice(100, 200), y=slice(100)),
]
merged = merge_arrays(arrays)

rds.where(rds!=rds.rio.nodata).plot();

merged.where(merged!=merged.rio.nodata).plot()

<matplotlib.collections.QuadMesh at 0x7f84f86f13c8>