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

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é.

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.

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:

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:

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

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.