¿Definiendo Winkel Tripel en el proyecto 4.8.0?

Tengo QGIS con el proyecto 4.8.0 instalado. Me gustaría agregar la proyección Winkel Tripel que parece definirse como

+proj=wintri

Pero esa cadena de definición no funciona. ¿Es esto un error en mi versión de proyecto, o cuál es la cadena de definición correcta?

qgis coordinate-system proj

— bajo oscuro
fuente

Funciona para mí: echo "34 -34" | proj "+ proj = wintri"

— Paul Ramsey

Sí, eso hace eco de algo. Pero si defino un CRS personalizado con esta cadena, nada se reproyecta en QGIS. ¿Supongo que necesitaría algunos parámetros más?

— oscuro

@Paul Ramsey: Parece que aunque la proyección Winkel Tripel se define en la projbiblioteca y se puede llamar desde la línea de comandos, no se puede usar como un CRS personalizado en QGIS porque no hay transformación inversa en la projbiblioteca (ver listas .maptools.org / pipermail / proj / 2012-May / 006302.html ). Y parece que la transformación inversa no es trivial (y siempre será una aproximación ...).

— Jake

Creo que deberías escribir eso en una respuesta, ya que es ... la respuesta. Las herramientas de línea de comandos también se niegan a hacer una proyección inversa.

— Paul Ramsey

Como los colegas mencionaron anteriormente, si también tiene proj + proj = aitoff tiene los mismos problemas, es una cuestión de proyección inversa faltante ver: (+ proj = aitoff; Rel. 4.7.1, 23 de septiembre de 2009 <invproj>: proyección inversa no disponible) y (invproj + proj = wintri; Rel. 4.7.1, 23 de septiembre de 2009 <invproj>: proyección inversa no disponible). Extraño pero bueno saberlo

— huckfinn

Resumiendo la discusión anterior:

Si bien la proyección de Winkel Tripel se define en la biblioteca de proyectos y se puede invocar desde la línea de comandos, no se puede usar como un CRS personalizado en QGIS porque no hay transformación inversa en la biblioteca de proyectos.

La solicitud de mejora para agregar esta funcionalidad se ha cerrado ya que parece que la transformación inversa no es trivial.

— bajo oscuro
fuente

¿Podría ser esto diferente ahora? Estoy, con pyproj==1.9.5.1instalado, capaz de hacer: lo pp = Proj(proj='wintri', lon_0=47.); pp(*pp(0,0), inverse=True)que devuelve (0,0)(a la precisión de la máquina), es decir, puedo revertir la transformación hacia adelante con la inverse=Truepalabra clave. Dado que este pyproj parece tener el inverso, ¿QGIS podría soportar esto ahora?

— Ahmed Fasih