Définir Winkel Tripel dans le projet 4.8.0?

J'ai QGIS avec proj 4.8.0 installé. Je voudrais ajouter la projection Winkel Tripel qui semble être définie comme

+proj=wintri

Mais cette chaîne de définition ne fonctionne pas. Est-ce un bug dans ma version proj, ou quelle est la chaîne de définition correcte?

qgis coordinate-system proj

— obscur
source

Fonctionne pour moi: écho "34 -34" | proj "+ proj = wintri"

— Paul Ramsey

Oui, cela fait écho à quelque chose. Mais si je définis un CRS personnalisé avec cette chaîne, rien n'est reprojeté dans QGIS. Je suppose que j'aurais besoin de plus de paramètres?

— underdark

@Paul Ramsey: Il semble que même si la projection Winkel Tripel est définie dans la projbibliothèque et peut être appelée à partir de la ligne de commande, elle ne peut pas être utilisée comme CRS personnalisé dans QGIS car il n'y a pas de transformation inverse dans la projbibliothèque (voir listes .maptools.org / pipermail / proj / 2012-May / 006302.html ). Et il semble que la transformation inverse n'est pas anodine (et sera toujours une approximation ...).

— Jake

Je pense que vous devriez écrire cela dans une réponse, car c'est ... la réponse. Les outils en ligne de commande refusent également de faire une projection inverse.

— Paul Ramsey

Comme les collègues l'ont mentionné ci-dessus, si vous avez également avec proj + proj = aitoff les mêmes problèmes, il s'agit de manquer la projection inverse voir: (+ proj = aitoff; Rel. 4.7.1, 23 septembre 2009 <invproj>: projection inverse non disponible) et (invproj + proj = wintri; Rel. 4.7.1, 23 septembre 2009 <invproj>: projection inverse non disponible). Étrange mais bon à savoir

— huckfinn

Résumant la discussion ci-dessus:

Bien que la projection Winkel Tripel soit définie dans la bibliothèque proj et puisse être appelée à partir de la ligne de commande, elle ne peut pas être utilisée comme CRS personnalisé dans QGIS car il n'y a pas de transformation inverse dans la bibliothèque proj.

La demande d'amélioration pour ajouter cette fonctionnalité a été fermée car il semble que la transformation inverse n'est pas anodine.

— obscur
source

Cela pourrait-il être différent maintenant? Je suis, avec pyproj==1.9.5.1installé, capable de faire: pp = Proj(proj='wintri', lon_0=47.); pp(*pp(0,0), inverse=True)qui retourne (0,0)(à la précision de la machine), c'est à dire, je suis capable d'inverser la transformation directe avec le inverse=Truemot - clé. Puisque ce pyproj semble avoir l'inverse, QGIS pourrait-il le supporter maintenant?

— Ahmed Fasih