È possibile utilizzare la proiezione dell'omolosina Goode in QGIS? Sia QGIS che GRASS sembrano avere problemi con le interruzioni, schiacciando le latitudini più elevate insieme. Ho provato + proj = goode + ellps = sphere + datum = wgs84 come CRS personalizzato.
È possibile utilizzare la proiezione dell'omolosina Goode in QGIS?
Risposte:
C'è stata una domanda sulle proiezioni interrotte prima: proiezioni cartografiche discontinue o interrotte?
La conclusione è stata che proj4 (che viene utilizzato per la gestione della proiezione in QGIS e GRASS) non supporta le proiezioni interrotte.
Avrei dovuto menzionare un biglietto sul tracker Proj.4 per il supporto di Goode Homolosine (che è stato creato dopo che quella domanda è stata posta): trac.osgeo.org/proj/ticket/106 . Sembra che non abbia un nuovo proj come pensavo, però; Proverò ad aggiornare il mio proj e vedrò come va. Ci scusiamo per la domanda duplicata.
—
Mattwigway,
Tutte le proiezioni sono interrotte ... vale a dire il contro dominio è uno spazio finito, non infinito;)
—
Luís de Sousa
Con QGIS Lisboa, esiste ora un supporto di base per la proiezione dell'omolosina Goode.
Ho scaricato i dati IGH da http://edc2.usgs.gov/glcc/fao/data_specifications.php
e configurare questo utente 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
per il Nord America.
lon_0 e x_0 provengono dall'Appendice 03 di questo documento: http://proceedings.esri.com/library/userconf/proc98/proceed/TO850/PAP844/P844.HTM
Sovrapponendo i confini del mondo della Terra naturale_m, ottengo questa immagine:
Adatto per il Nord America (circa 20 km fuori luogo) e distorto in altre regioni.
MODIFICARE
Ho appena notato che proj ormai conosce l'omolosina interrotta:
+proj=igh +ellps=sphere +towgs84=0,0,0 +lon_0=100w +x_0=-11119487.43
mi dà questa immagine dei confini del mondo delle terre naturali e una griglia vettoriale di 10 gradi:
Non sono più sicuro di quali parametri a parte +proj=ighsiano più necessari.
Usando la definizione sopra ottengo uno spostamento dei pixel. Per me la definizione ha
—
Thomas,
+proj=igh +towgs84=0,0,0funzionato bene
QGIS 2.18 si sta avvicinando al pieno supporto dell'omolosina per i raster. Devi solo essere consapevole del problema n. 959 : esegui qualsiasi deformazione raster direttamente con GDAL utilizzando i CHECK_WITH_INVERT_PROJparametri e sei a posto.
Per quanto riguarda i vettori, le cose sono ancora complesse. Prevalgono alcuni degli errori che AndreJ mostra nella sua risposta . La causa è semplice: QGis non è consapevole del fatto che si tratta di un programma cartografico e considera il co-dominio di qualsiasi proiezione come uno spazio finito.
Questo ha una soluzione semplice: interseca qualsiasi livello vettoriale importato in QGis con il co-dominio Homolosine. Qualche tempo fa ho reso disponibile un vettore di co-domini discreto definito ad intervalli di un grado di longitudine e latitudine.
