Tengo instalado QGIS con proj 4.8.0. Me gustaría añadir la proyección Winkel Tripel que parece estar definida como
+proj=wintriPero esa cadena de definición no funciona. ¿Es un error de mi versión del proyecto o cuál es la cadena de definición correcta?
Resumiendo la discusión anterior:
Aunque la proyección Winkel Tripel está definida en la biblioteca proj y puede ser llamada desde la línea de comandos, no puede ser utilizada como un CRS personalizado en QGIS porque hay ninguna transformación inversa en la biblioteca del proyecto.
La solicitud de mejora para añadir esta funcionalidad se ha cerrado ya que parece que la transformación inversa no es trivial.
¿Podría ser diferente ahora? Lo soy, con pyproj==1.9.5.1 instalado, capaz de hacer: pp = Proj(proj='wintri', lon_0=47.); pp(*pp(0,0), inverse=True) que devuelve (0,0) (a precisión de máquina), es decir, soy capaz de invertir la transformación hacia delante con el inverse=True palabra clave. Dado que este pyproj parece tener la inversa, ¿podría QGIS soportar esto ahora?
I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.
1 votos
A mí me funciona: echo "34 -34" | proj "+proj=wintri"
1 votos
Sí, eso me recuerda algo. Pero si defino un CRS personalizado con esta cadena, no se reproyecta nada en QGIS. ¿Supongo que necesitaría más parámetros?
0 votos
No, el hecho de que se reproyecte en la línea de comandos parece indicar que funciona. ¿Es posible que QGIS tenga una librería proj diferente a la que cree tener? ¿O que QGIS realmente compruebe las cadenas de proyección por sí mismo antes de pasarlas a proj? (¿Y aún no ha oído hablar de Winkel?)
7 votos
@Paul Ramsey: Parece que aunque la proyección Winkel Tripel está definida en el
projy se puede llamar desde la línea de comandos, no se puede utilizar como un CRS personalizado en QGIS porque no hay ninguna transformación inversa en elprojbiblioteca (véase lists.maptools.org/pipermail/proj/2012-May/006302.html ). Y parece que la transformación inversa no es trivial (y siempre será una aproximación...).8 votos
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.
0 votos
Supongo que se puede añadir la inversa a Proj4 o tal vez incluso como un plugin de QGIS (ver github.com/d3/d3-geo-projection/blob/master/src/winkel3.js#L12 para una implementación de la inversa en Javascript)? Tratando de aprender más acerca de estas dos opciones.
1 votos
Como mencionan los compañeros arriba, si tienes también con proj +proj=aitoff tienes los mismos problemas, es cuestión de que falta la proyección inversa ver: (+proj=aitoff; Rel. 4.7.1, 23 septiembre 2009 <invproj>: proyección inversa no disponible) y ( invproj +proj=wintri; Rel. 4.7.1, 23 septiembre 2009 <invproj>: proyección inversa no disponible ). Extraño pero bueno de saber
0 votos
He añadido una solicitud de mejora en el rastreador de problemas de PROJ: trac.osgeo.org/proj/ticket/234
0 votos
... que fue rechazado porque no hay una solución matemática fácil para el proyecto inverso