Est-il possible d'utiliser la projection Goode Homolosine dans QGIS?


11

Est-il possible d'utiliser la projection Goode Homolosine dans QGIS? QGIS et GRASS semblent avoir des problèmes avec les interruptions, écrasant les latitudes plus élevées ensemble. J'ai essayé + proj = goode + ellps = sphère + datum = wgs84 en tant que CRS personnalisé.

Réponses:


0

Une question a déjà été posée sur les projections interrompues: projections cartographiques discontinues ou interrompues?

La conclusion était que proj4 (qui est utilisé pour la gestion des projections dans QGIS et GRASS) ne prend pas en charge les projections interrompues.


J'aurais dû mentionner un ticket sur le tracker Proj.4 pour le support Goode Homolosine (qui a été créé après que cette question a été posée): trac.osgeo.org/proj/ticket/106 . On dirait que je n'ai pas un projet aussi nouveau que je le pensais, cependant; Je vais essayer de mettre à jour mon projet et voir comment ça se passe. Désolé pour la question en double.
mattwigway

1
Toutes les projections sont interrompues ... c'est-à-dire que le contre-domaine est un espace fini, pas infini;)
Luís de Sousa

6

Avec QGIS Lisboa, il existe désormais un support de base pour la projection Goode Homolosine.

J'ai téléchargé les données IGH depuis http://edc2.usgs.gov/glcc/fao/data_specifications.php

et configurer cet utilisateur 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

pour l'Amérique du Nord.

lon_0 et x_0 proviennent de l'annexe 03 de ce document: http://proceedings.esri.com/library/userconf/proc98/proceed/TO850/PAP844/P844.HTM

Superposant les limites naturelles du monde terrestre_m, j'obtiens cette image:

entrez la description de l'image ici

Convient bien en Amérique du Nord (environ 20 km mal placé) et déformé dans d'autres régions.


ÉDITER

Je viens de remarquer que proj sait maintenant que l'homolosine a été interrompue:

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

me donne cette image des limites du monde des terres naturelles et une grille vectorielle à 10 degrés:

entrez la description de l'image ici

Je ne sais plus quels paramètres à part +proj=ighsont nécessaires.


En utilisant la définition ci-dessus, j'obtiens un décalage de pixel. Pour moi, la définition +proj=igh +towgs84=0,0,0a très bien fonctionné
Thomas

3

QGIS 2.18 est sur le point de prendre pleinement en charge l'Homolosine pour les rasters. Vous avez juste besoin d'être conscient du problème # 959 : effectuez une déformation raster directement avec GDAL en utilisant les CHECK_WITH_INVERT_PROJparamètres et vous êtes bon.

En ce qui concerne les vecteurs, les choses sont encore complexes. Certaines des erreurs qu'AndreJ montre dans sa réponse prévalent. La cause est simple: QGis n'est pas conscient qu'il s'agit d'un programme Cartographique et considère le co-domaine de toute projection comme un espace fini.

Cela a une solution simple: intersecter toute couche vectorielle importée dans QGis avec le co-domaine Homolosine. Il y a quelque temps, j'ai rendu disponible un vecteur de co-domaine discret défini à des intervalles d'un degré de longitude et de latitude.

EST CE QUE JE

entrez la description de l'image ici

En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.