Definindo Winkel Tripel no proj 4.8.0?

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?

qgis coordinate-system proj

— underdark
fonte

Funciona para mim: eco "34 -34" | proj "+ proj = wintri"

— Paul Ramsey

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

@ 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 #

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

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

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