2020-06-03 14:20:17 -07:00

30 lines
702 B
Python

#!/usr/bin/env python
# Lint as: python3
"""Reverse mask a region.
Create an image that masks everything except for the specified polygon.
"""
import ee
import ee.mapclient
ee.Initialize()
ee.mapclient.centerMap(-115, 39, 5)
fc = (ee.FeatureCollection('RESOLVE/ECOREGIONS/2017')
.filter(ee.Filter().eq('ECO_NAME', 'Great Basin shrub steppe')))
# Start with a black image.
empty = ee.Image(0).toByte()
# Fill and outline the polygons in two colors
filled = empty.paint(fc, 2)
both = filled.paint(fc, 1, 5)
# Mask off everything that matches the fill color.
result = both.mask(filled.neq(2))
ee.mapclient.addToMap(result, {
'palette': '000000,FF0000',
'max': 1,
'opacity': 0.5
})