Reproyección del motor de Google Earth a CRS no definidos por la EPSG

Question

¿Es posible reproyectar una imagen en el motor terrestre a coordenadas no definidas por EPSG? En particular, Goode Homosline (SR-ORG:7619) o Mollweide (ESRI:54009). Pasando cualquiera de las dos a .reproject() devuelve un error que indica que no se puede analizar el CRS.

Answer 1

Sí, el Motor Tierra puede reproyectar imágenes a sistemas de referencia de coordenadas no EPSG.

Earth Engine utiliza el lenguaje Java GeoTools para transformaciones del sistema de coordenadas que admite la Sistema de coordenadas Texto conocido (WKT) y un conjunto de proyecciones de base.

Mollweide

En ESRI:54009 Mundo Mollweide sistema de referencia de coordenadas (SRC) tiene un Representación WKT que incluye una entrada de proyección de base PROJECTION["Mollweide"] que es compatible con GeoTools. El siguiente código de Earth Engine muestra cómo reproyectar a World Mollweide:

// Get a sample image (Landsat 8) which uses a UTM projection.
var image = ee.Image('LANDSAT/LC08/C01/T1_TOA/LC08_044034_20170614');

// Get the bounding coordinates of the Landsat image.
var bounds = image.geometry().bounds().getInfo()['coordinates'];

print('Original image (UTM) URL:',
  image.getThumbURL({
    bands: ['B4', 'B3', 'B2'],
    min: 0,
    max: 0.3,
    dimensions: '400'
  })
);

// Project the image to Mollweide.
var wkt = ' \
  PROJCS["World_Mollweide", \
    GEOGCS["GCS_WGS_1984", \
      DATUM["WGS_1984", \
        SPHEROID["WGS_1984",6378137,298.257223563]], \
      PRIMEM["Greenwich",0], \
      UNIT["Degree",0.017453292519943295]], \
    PROJECTION["Mollweide"], \
    PARAMETER["False_Easting",0], \
    PARAMETER["False_Northing",0], \
    PARAMETER["Central_Meridian",0], \
    UNIT["Meter",1], \
    AUTHORITY["EPSG","54009"]]';

var proj_mollweide = ee.Projection(wkt);
var image_mollweide = image.reproject({
  crs: proj_mollweide,
  scale: 300
});

print('Projected image (World Mollweide) URL:',
  image_mollweide.getThumbURL({
    bands: ['B4', 'B3', 'B2'],
    min: 0,
    max: 0.3,
    region: bounds,
    dimensions: '400'
  })
);

(Enlace al código: https://code.earthengine.google.com/5b50ad250e83ec994335615c4db72a00 )

La imagen original en una proyección UTM tiene este aspecto:

Mientras que las imágenes proyectadas World Mollweide parece:

Goode Homosline

En cambio, el SR-ORG:7619 (Interrupted Goode Homosline) tiene un CRS WKT que incluye una entrada de proyección de base PROJECTION["Interrupted_Goode_Homolosine"] que no es compatible con GeoTools, por lo que no es compatible con Earth Engine.