¿Es posible usar la proyección Goode Homolosine en QGIS?

11

¿Es posible usar la proyección Goode Homolosine en QGIS? Tanto QGIS como GRASS parecen tener problemas con las interrupciones, aplastando las latitudes más altas. Intenté + proj = goode + ellps = sphere + datum = wgs84 como un CRS personalizado.

Mattwigway
fuente

Respuestas:

0

Ha habido una pregunta sobre proyecciones interrumpidas antes: ¿Proyecciones de mapas discontinuas o interrumpidas?

La conclusión fue que proj4 (que se usa para el manejo de proyecciones en QGIS y GRASS) no admite proyecciones interrumpidas.

bajo oscuro
fuente
Debería haber mencionado un ticket en el rastreador Proj.4 para el soporte de Goode Homolosine (que fue creado después de que se hizo esa pregunta): trac.osgeo.org/proj/ticket/106 . Sin embargo, parece que no tengo un proyecto tan nuevo como pensaba; Intentaré actualizar mi proyecto y veré cómo funciona. Perdón por la pregunta duplicada.
Mattwigway
1
Todas las proyecciones se interrumpen ... es decir, el contradominio es un espacio finito, no infinito;)
Luís de Sousa
6

Con QGIS Lisboa, ahora hay un soporte básico para la proyección Goode Homolosine.

Descargué datos de IGH de http://edc2.usgs.gov/glcc/fao/data_specifications.php

y configure este usuario CRS:

+proj=goode +lon_0=-100 +x_0=-11119487.43 +y_0=0 +a=6370997 +b=6370997 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs

para América del Norte

lon_0 y x_0 son del Apéndice 03 de este documento: http://proceedings.esri.com/library/userconf/proc98/proceed/TO850/PAP844/P844.HTM

Superponiendo los límites naturales del mundo terrestre_m, obtengo esta imagen:

ingrese la descripción de la imagen aquí

Se ajusta bien en América del Norte (unos 20 km fuera de lugar) y distorsionado en otras regiones.


EDITAR

Acabo de darme cuenta de que el proyecto ya conoce la goolesolina interrumpida:

+proj=igh +ellps=sphere +towgs84=0,0,0 +lon_0=100w +x_0=-11119487.43

me da esta imagen de los límites mundiales de las Tierras Naturales y una cuadrícula de vectores de 10 grados:

ingrese la descripción de la imagen aquí

Ya no estoy seguro de qué parámetros +proj=ighson necesarios.

AndreJ
fuente
Usando la definición anterior, obtengo un desplazamiento de píxeles. Para mí, la definición +proj=igh +towgs84=0,0,0funcionó bien
Thomas
3

QGIS 2.18 está a punto de ser totalmente compatible con Homolosine para rásteres. Solo necesita estar al tanto del problema # 959 : realice cualquier deformación de trama directamente con GDAL utilizando los CHECK_WITH_INVERT_PROJparámetros y estará bien.

En cuanto a los vectores, las cosas siguen siendo complejas. Algunos de los errores que AndreJ muestra en su respuesta prevalecen. La causa es simple: QGis no es consciente de que es un programa cartográfico y considera que el codominio de cualquier proyección es un espacio finito.

Esto tiene una solución simple: intersecar cualquier capa vectorial importada en QGis con el codominio Homolosine. Hace algún tiempo, puse a disposición un vector de codominio discreto definido a intervalos de un grado de longitud y latitud.

DOI

ingrese la descripción de la imagen aquí

Luís de Sousa
fuente