Monday, 12 August 2019

coordinate system - LSIB feature does not align with island boundary in Google Earth Engine?


It looks like there is an alignment issue that arises in EarthEngine with the LSIB polygons. It exists for Greenland in both the Detailed LSIB 2013 and Simple LSIB 2017 datasets, and can be verified with both Map and Satellite layers in the map pane. At low latitudes the problem does not exist but at high latitudes the polygon seems to be poorly drawn, and fails to cover regions that are clearly terrestrial but does encapsulate marine regions.


Is this a projection issue, a polygon issue, or both?


It doesn't seem to exist in other high-Arctic polygons including Russia, Canada, and Svalbard.



// Greenland
var GL = ee.FeatureCollection('USDOS/LSIB/2013')
.filter(ee.Filter.eq('cc','GL'));
Map.setCenter(-25, 80, 4); // Center on NE GL.
Map.addLayer(GL, {color: 'FF3333', opacity:0.99}, 'Greenland_Detailed');

var GL2 = ee.FeatureCollection('USDOS/LSIB_SIMPLE/2017')
.filter(ee.Filter.eq('country_co','GL'));
Map.addLayer(GL2, {color: 'FF3333', opacity:0.99}, 'Greenland_Simple');


// Russia
var RU = ee.FeatureCollection('USDOS/LSIB/2013')
.filter(ee.Filter.eq('cc','RS'));
print(RU);
Map.addLayer(RU, {color: '00FF00', opacity:0.99}, 'Russia');

// Svalbard
var SV = ee.FeatureCollection('USDOS/LSIB/2013')
.filter(ee.Filter.eq('cc','SV'));
Map.addLayer(SV, {color: 'FFFF00', opacity:0.99}, 'Svalbard');


// Canada
var CA = ee.FeatureCollection('USDOS/LSIB/2013')
.filter(ee.Filter.eq('cc','CA'));
Map.addLayer(CA, {color: '0000FF', opacity:0.99}, 'Canada');


No comments:

Post a Comment

arcpy - Changing output name when exporting data driven pages to JPG?

Is there a way to save the output JPG, changing the output file name to the page name, instead of page number? I mean changing the script fo...