Regions

[mafla]

I don't request to change this. Just to remember that I find it strange that regions are always "bounding boxes" and not just the regions. Rationale: If you want to compute a statistic over e.g. some SREX region like CEU (which is a triangle) it is indeed actually nice to see the data in the bounding box, but if I want to do some statistics on the whole thing you will need to go via the mapCube-type of aggregation instead of simply mapslicsing a mean.

[meggart]

Can we mark this as solved through the new shapefile interface? It would still be in the responsibility of the user to do the masking then. I could imagine having a list of pre-stored shapefiles and then we could provide a function like getregionmask("CEU") that returns the rasterized shape for CEU over even a function like maskregion("CEU") that turns values outside the region into missings. Any thoughts?

[mafla]

Sounds perfect! 👍 Need an example online ...