rioxarray_geotiff_load_mask_nodata_calculate_mean.py

1import rioxarray
2
3# Load a raster dataset
4rds = rioxarray.open_rasterio("path/to/raster.tif")
5
6# Mask fill values (standard practice to handle NoData)
7rds = rds.where(rds != rds.rio.nodata)
8
9# Calculate the mean across the spatial dimensions
10mean_val = rds.mean().values
11print(f"Mean value of the raster: {mean_val}")