Example - Resampling

This example demonstrates how to reproduce rasterio’s resampling example here.

from rasterio.enums import Resampling

import rioxarray

%matplotlib inline

Load in xarray dataset

See docs for rioxarray.open_rasterio

Notes:

• masked=True will convert from integer to float64 and fill with NaN. If this behavior is not desired, you can skip this.

xds = rioxarray.open_rasterio(
    "../../test/test_data/compare/small_dem_3m_merged.tif",
    masked=True,
)

Resampling

API Reference for rio.reproject:

• DataArray.reproject

• Dataset.reproject

upscale_factor = 2
new_width = xds.rio.width * upscale_factor
new_height = xds.rio.height * upscale_factor

xds_upsampled = xds.rio.reproject(
    xds.rio.crs,
    shape=(new_height, new_width),
    resampling=Resampling.bilinear,
)

xds.shape

(1, 245, 574)

xds_upsampled.shape

(1, 490, 1148)

xds.rio.resolution()

(3.0, -3.0)

xds_upsampled.rio.resolution()

(1.5, -1.5)