Skip to content

Instantly share code, notes, and snippets.

@mojodna

mojodna/README.md

Last active Oct 26, 2018
Embed
What would you like to do?
ESACCI-LC Mid Latitude mask

Mid Latitude masking

First, apply the mid latitude mask as mask.

Next, add envelopes for each of the C-CAP sources as masks (otherwise chunky pixels sneak in, per https://github.com/mojodna/marblecutter-land-cover/issues/6):

update land_cover set mask = ST_Union(mask, (select ST_Union(ST_Envelope(geom)) from land_cover where source='C-CAP')) where source='ESACCI-LC Mid Latitudes';

Then, update the footprint to exclude the mask (to ensure that appropriate sources are selected):

update land_cover set geom=ST_Multi(ST_Difference(geom, mask)) where source='ESACCI-LC Mid Latitudes';
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
You can’t perform that action at this time.