create masks of geospatial regions
Determine which geographic region each grid point belongs to.
When working with gridded data such as climate model output or reanalysis data it is often important to create regional averages, e.g., over countries, continents or regions defined in the literature. To do so we need to know for each grid point to which region it belongs to.
regionmask is a python package that:
can be used to create masks of (geo)spatial regions for arbitrary (longitude and latitude) grids. These masks indicate which region a gridpoint belongs to and can then be used to aggregate gridded data over the regions. The masks come in three variants:
takes great care to create consistent masks for
contains a number of defined regions, including:
supports using arbitrary existing or user-defined region definitions:
regions defined as shapefiles can be accessed via geopandas
user-defined regions can be created via numpy or shapely
License
regionmask is published under a MIT license.