Definiowanie Winkel Tripel w wersji 4.8.0?

Mam QGIS z zainstalowaną wersją proj 4.8.0. Chciałbym dodać projekcję Winkel Tripel, która wydaje się być zdefiniowana jako

+proj=wintri

Ale ten ciąg definicji nie działa. Czy to błąd w mojej wersji proj, czy jaki jest prawidłowy ciąg definicji?

qgis coordinate-system proj

— podmrok
źródło

Działa dla mnie: echo „34 -34” | proj "+ proj = wintri"

— Paul Ramsey

Tak, to coś echa. Ale jeśli zdefiniuję niestandardowy CRS za pomocą tego ciągu, nic nie zostanie ponownie skierowane w QGIS. Chyba potrzebuję więcej parametrów?

— podmrok

@Paul Ramsey: Wygląda na to, że chociaż projekcja Winkel Tripel jest zdefiniowana w projbibliotece i można ją wywołać z wiersza poleceń, nie można jej użyć jako niestandardowego CRS w QGIS, ponieważ nie ma odwrotnej transformacji w projbibliotece (patrz listy .maptools.org / pipermail / proj / 2012-May / 006302.html ). I wydaje się, że odwrotna transformacja nie jest trywialna (i zawsze będzie przybliżeniem ...).

— Jake

Myślę, że powinieneś zapisać to w odpowiedzi, ponieważ jest to ... odpowiedź. Narzędzia wiersza poleceń również odmawiają wykonania odwrotnej projekcji.

— Paul Ramsey,

Jak wspomniano powyżej, jeśli masz również takie same problemy z proj + proj = aitoff, kwestią jest brak odwrotnej projekcji patrz: (+ proj = aitoff; Rel. 4.7.1, 23 września 2009 r. <invproj>: odwrotna projekcja niedostępna) i (invproj + proj = wintri; Rel. 4.7.1, 23 września 2009 <invproj>: odwrotna projekcja niedostępna). Dziwne, ale warto wiedzieć

— huckfinn

Podsumowując powyższą dyskusję:

Podczas gdy projekcja Winkel Tripel jest zdefiniowana w bibliotece proj i można ją wywoływać z wiersza poleceń, nie można jej używać jako niestandardowego CRS w QGIS, ponieważ nie ma odwrotnej transformacji w bibliotece proj.

Żądanie rozszerzenia o dodanie tej funkcji zostało zamknięte, ponieważ wydaje się, że odwrotna transformacja nie jest trywialna.

— podmrok
źródło

Czy teraz może być inaczej? Jestem, z pyproj==1.9.5.1zainstalowanym, w stanie zrobić: pp = Proj(proj='wintri', lon_0=47.); pp(*pp(0,0), inverse=True)który zwraca (0,0)(do precyzji maszyny), tzn. Jestem w stanie odwrócić transformację do przodu za pomocą inverse=Truesłowa kluczowego. Ponieważ ten pyproj wydaje się mieć odwrotność, czy QGIS może teraz to obsługiwać?

— Ahmed Fasih