Definindo Winkel Tripel no proj 4.8.0?

14

Eu tenho o QGIS com o proj 4.8.0 instalado. Eu gostaria de adicionar a projeção de Winkel Tripel, que parece ser definida como

+proj=wintri

Mas essa sequência de definições não funciona. Isso é um bug na minha versão do proj ou qual é a string de definição correta?

underdark
fonte
1
Funciona para mim: eco "34 -34" | proj "+ proj = wintri"
Paul Ramsey
1
Sim, isso ecoa alguma coisa. Mas se eu definir um CRS personalizado com essa sequência, nada será reprojetado no QGIS. Acho que precisaria de mais alguns parâmetros?
Underdark
7
@ Paul Ramsey: Parece que, embora a projeção de Winkel Tripel seja definida na projbiblioteca e possa ser chamada na linha de comando, ela não pode ser usada como um CRS personalizado no QGIS porque não há transformação inversa na projbiblioteca (consulte as listas .maptools.org / pipermail / proj / 2012-May / 006302.html ). E parece que a transformação inversa não é trivial (e sempre será uma aproximação ...).
Jake #
8
Eu acho que você deveria escrever isso em uma resposta, como é ... a resposta. As ferramentas de linha de comando também se recusam a fazer uma projeção inversa.
Paul Ramsey
1
Como os colegas mencionaram acima, se você também tem proj + proj = aitoff com os mesmos problemas, é uma questão de falta de projeção inversa, consulte: (+ proj = aitoff; Rel. 4.7.1, 23 de setembro de 2009 <invproj>: projeção inversa não disponível) e (invproj + proj = wintri; Rel. 4.7.1, 23 de setembro de 2009 <invproj>: projeção inversa não disponível). Estranho, mas bom saber
huckfinn

Respostas:

4

Resumindo a discussão acima:

Embora a projeção Winkel Tripel seja definida na biblioteca do proj e possa ser chamada na linha de comando, ela não pode ser usada como um CRS personalizado no QGIS, porque não há transformação inversa na biblioteca do proj.

A solicitação de aprimoramento para adicionar essa funcionalidade foi encerrada, pois parece que a transformação inversa não é trivial.

underdark
fonte
Isso poderia ser diferente agora? Eu sou, com pyproj==1.9.5.1instalado, capaz de fazer: o pp = Proj(proj='wintri', lon_0=47.); pp(*pp(0,0), inverse=True)que retorna (0,0)(à precisão da máquina), ou seja, sou capaz de reverter a transformação direta com a inverse=Truepalavra - chave. Como esse pyproj parece ter o inverso, o QGIS poderia suportar isso agora?
Ahmed Fasih