Image masking based on a geometry

Add an image operator for applying a mask on a GridCoverage: all pixel values outside a geometry would be replaced by a fill value. The coverage and the geometry are in arbitrary CRS, not necessarily the same.