Winkel Tripel in Proj 4.8.0 definieren?

14

Ich habe QGIS mit ProJ 4.8.0 installiert. Ich möchte Winkel Tripel Projektion hinzufügen, die als definiert zu sein scheint

+proj=wintri

Diese Definitionszeichenfolge funktioniert jedoch nicht. Ist das ein Fehler in meiner ProJ-Version oder was ist die richtige Definitionszeichenfolge?

qgis coordinate-system proj Underdunkel
quelle
1
Funktioniert bei mir: Echo "34 -34" | proj "+ proj = wintri"
Paul Ramsey
1
Ja, das klingt nach etwas. Wenn ich mit dieser Zeichenfolge ein benutzerdefiniertes CRS definiere, wird in QGIS nichts neu projiziert. Vermutlich brauche ich noch ein paar Parameter?
Underdunkel
7
@Paul Ramsey: Obwohl die Winkel Tripel-Projektion in der projBibliothek definiert ist und über die Befehlszeile aufgerufen werden kann, kann sie in QGIS nicht als benutzerdefiniertes CRS verwendet werden, da es in der projBibliothek keine inverse Transformation gibt (siehe Listen) .maptools.org / pipermail / proj / 2012-May / 006302.html ). Und es scheint, dass die inverse Transformation nicht trivial ist (und immer eine Annäherung sein wird ...).
Jake
8
Ich denke, Sie sollten das in eine Antwort schreiben, da es ... die Antwort ist. Die Befehlszeilen-Tools lehnen auch eine inverse Projektion ab.
Paul Ramsey
1
Wie die Kollegen oben erwähnt haben, wenn Sie auch mit proj + proj = aitoff die gleichen Probleme haben, ist es eine Frage der fehlenden inversen Projektion, siehe: (+ proj = aitoff; Rel. 4.7.1, 23. September 2009 <invproj>: Inverse Projektion nicht verfügbar) und (invproj + proj = wintri; Rel. 4.7.1, 23. September 2009 <invproj>: Inverse Projektion nicht verfügbar). Seltsam, aber gut zu wissen
huckfinn

Antworten:

4

Fassen Sie die obige Diskussion zusammen:

Während Winkel Tripel Projektion in der proj Bibliothek definiert und kann von der Kommandozeile aufgerufen wird, kann sie nicht als eigenes KBS in QGIS verwendet werden , weil es keine inverse Transformation in der proj Bibliothek.

Die Erweiterungsanforderung zum Hinzufügen dieser Funktionalität wurde geschlossen, da die inverse Transformation anscheinend nicht trivial ist.

Underdunkel
quelle
Könnte das jetzt anders sein? Ich bin mit pyproj==1.9.5.1installiertem in der Lage zu tun, pp = Proj(proj='wintri', lon_0=47.); pp(*pp(0,0), inverse=True)was (0,0)(zur Maschinengenauigkeit) zurückgibt , dh ich bin in der Lage, die Vorwärtstransformation mit dem inverse=TrueSchlüsselwort umzukehren . Könnte QGIS dies jetzt unterstützen, da dieses Pyproj das Gegenteil zu haben scheint?
Ahmed Fasih