Masking satellite data using shapefile

Question

I am using Landsat data to determine NDVI for a specific region of interest (roi). However, the roi has buildings and other artificial objects that I wish to mask.

I have followed through the example at Masking NetCDF time series data from shapefile using Python.

However, this method masks all data outside of the shapefile, and retains the data from inside of the shapefile.

Has anybody any solutions for doing the reverse of this, and thus masking the data INSIDE the shapefile?

Answer 1

https://corteva.github.io/rioxarray/stable/examples/clip_geom.html

You just need to pass in invert=True:

clipped = xds.rio.clip(......, invert=True)