Definire Winkel Tripel in proj 4.8.0?

Ho installato QGIS con proj 4.8.0. Vorrei aggiungere la proiezione di Winkel Tripel che sembra essere definita come

+proj=wintri

Ma quella stringa di definizione non funziona. È un bug nella mia versione proj o qual è la stringa di definizione corretta?

qgis coordinate-system proj

— underdark
fonte

Funziona per me: echo "34 -34" | proj "+ proj = wintri"

— Paul Ramsey,

Sì, echo qualcosa. Ma se definisco un CRS personalizzato con questa stringa, nulla viene riproiettato in QGIS. Immagino che avrei bisogno di qualche altro parametro?

— underdark

@Paul Ramsey: sembra che anche se la proiezione di Winkel Tripel è definita nella projlibreria e può essere chiamata dalla riga di comando, non può essere utilizzata come CRS personalizzato in QGIS perché non c'è trasformazione inversa nella projlibreria (vedi elenchi .maptools.org / pipermail / proj / 2012-mag / 006302.html ). E sembra che la trasformazione inversa non sia banale (e sarà sempre un'approssimazione ...).

— Jake,

Penso che dovresti scriverlo in una risposta, in quanto è ... la risposta. Gli strumenti da riga di comando rifiutano anche di fare una proiezione inversa.

— Paul Ramsey,

Come menzionato sopra dai colleghi, se anche con proj + proj = aitoff hanno gli stessi problemi, è una questione di mancata proiezione inversa vedi: (+ proj = aitoff; Rel. 4.7.1, 23 settembre 2009 <invproj>: proiezione inversa non disponibile) e (invproj + proj = wintri; Rel. 4.7.1, 23 settembre 2009 <invproj>: proiezione inversa non disponibile). Strano ma buono a sapersi

— huckfinn

Riassumendo la discussione sopra:

Mentre la proiezione di Winkel Tripel è definita nella libreria proj e può essere chiamata dalla riga di comando, non può essere utilizzata come CRS personalizzato in QGIS perché non c'è trasformazione inversa nella libreria proj.

La richiesta di miglioramento per aggiungere questa funzionalità è stata chiusa poiché sembra che la trasformazione inversa non sia banale.

— underdark
fonte

Questo potrebbe essere diverso ora? Sono, con pyproj==1.9.5.1installato, in grado di fare: pp = Proj(proj='wintri', lon_0=47.); pp(*pp(0,0), inverse=True)quale ritorna (0,0)(alla precisione della macchina), cioè sono in grado di invertire la trasformazione in avanti con la inverse=Trueparola chiave. Dato che questo pyproj sembra avere l'inverso, QGIS potrebbe supportarlo ora?

— Ahmed Fasih,