¿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.
¿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.
Respuestas:
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.
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:
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:
Ya no estoy seguro de qué parámetros +proj=igh
son necesarios.
+proj=igh +towgs84=0,0,0
funcionó bien
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_PROJ
pará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.