Mendefinisikan Winkel Tripel dalam proj 4.8.0?

14

Saya memiliki QGIS dengan proj 4.8.0 diinstal. Saya ingin menambahkan proyeksi Winkel Tripel yang tampaknya didefinisikan sebagai

+proj=wintri

Tapi string definisi itu tidak berfungsi. Apakah ini bug dalam versi proj saya, atau apa string definisi yang benar?

underdark
sumber
1
Bekerja untuk saya: echo "34 -34" | proj "+ proj = wintri"
Paul Ramsey
1
Ya, itu menggemakan sesuatu. Tetapi jika saya mendefinisikan CRS khusus dengan string ini, tidak ada yang diproyeksikan ulang di QGIS. Kira saya perlu beberapa parameter lagi?
underdark
7
@ Paul Ramsey: Tampaknya meskipun proyeksi Winkel Tripel didefinisikan di projperpustakaan dan dapat dipanggil dari baris perintah, itu tidak dapat digunakan sebagai CRS khusus di QGIS karena tidak ada transformasi terbalik di projperpustakaan (lihat daftar .maptools.org / pipermail / proj / 2012-Mei / 006302.html ). Dan tampaknya transformasi terbalik tidak sepele (dan akan selalu menjadi perkiraan ...).
Jake
8
Saya pikir Anda harus menulis itu menjadi jawaban, karena itu ... jawabannya. Alat baris perintah juga menolak untuk melakukan proyeksi terbalik.
Paul Ramsey
1
Seperti yang disebutkan oleh kolega di atas, jika Anda juga memiliki proj + proj = aitoff memiliki masalah yang sama, itu adalah masalah proyeksi invers yang hilang, lihat: (+ proj = aitoff; Rel 4.7.1, 23 September 2009 <invproj>: proyeksi terbalik tidak tersedia) dan (invproj + proj = wintri; Rel 4.7.1, 23 September 2009 <invproj>: proyeksi terbalik tidak tersedia). Aneh tapi bagus untuk diketahui
huckfinn

Jawaban:

4

Ringkas diskusi di atas:

Sementara proyeksi Winkel Tripel didefinisikan dalam pustaka proj dan dapat dipanggil dari baris perintah, itu tidak dapat digunakan sebagai CRS kustom di QGIS karena tidak ada transformasi terbalik di pustaka proj.

Permintaan perangkat tambahan untuk menambahkan fungsionalitas ini telah ditutup karena tampaknya transformasi terbalik tidak sepele.

underdark
sumber
Mungkinkah ini berbeda sekarang? Saya, dengan pyproj==1.9.5.1terinstal, dapat melakukan: pp = Proj(proj='wintri', lon_0=47.); pp(*pp(0,0), inverse=True)yang mengembalikan (0,0)(ke presisi mesin), yaitu, saya dapat membalikkan transformasi maju dengan inverse=Truekata kunci. Karena pyproj ini tampaknya memiliki kebalikannya, dapatkah QGIS mendukung ini sekarang?
Ahmed Fasih