Definiowanie Winkel Tripel w wersji 4.8.0?

15

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?

podmrok
źródło
1
Działa dla mnie: echo „34 -34” | proj "+ proj = wintri"
Paul Ramsey
1
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
7
@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
8
Myślę, że powinieneś zapisać to w odpowiedzi, ponieważ jest to ... odpowiedź. Narzędzia wiersza poleceń również odmawiają wykonania odwrotnej projekcji.
Paul Ramsey,
1
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

Odpowiedzi:

5

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
1
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