Apakah mungkin untuk menggunakan proyeksi Goode Homolosine di QGIS?

11

Apakah mungkin untuk menggunakan proyeksi Goode Homolosine di QGIS? Baik QGIS dan GRASS tampaknya memiliki masalah dengan interupsi, menekan garis lintang yang lebih tinggi bersama. Saya mencoba + proj = selamat tinggal + ellps = bola + datum = wgs84 sebagai CRS kustom.

mattwigway
sumber

Jawaban:

0

Ada pertanyaan tentang proyeksi yang terputus sebelumnya: Proyeksi peta yang terputus atau terputus?

Kesimpulannya adalah bahwa proj4 (yang digunakan untuk penanganan proyeksi di QGIS dan GRASS) tidak mendukung proyeksi yang terganggu.

underdark
sumber
Saya seharusnya menyebutkan tiket pada pelacak Proj.4 untuk dukungan Goode Homolosine (yang dibuat setelah pertanyaan itu diajukan): trac.osgeo.org/proj/ticket/106 . Sepertinya saya tidak memiliki proj yang baru seperti yang saya pikirkan; Saya akan mencoba memperbarui proj saya dan melihat bagaimana hasilnya. Maaf untuk pertanyaan rangkap.
mattwigway
1
Semua proyeksi terputus ... yaitu domain penghitung adalah ruang terbatas, bukan tak terbatas;)
Luís de Sousa
6

Dengan QGIS Lisboa, sekarang ada dukungan dasar untuk proyeksi Goode Homolosine.

Saya mengunduh data IGH dari http://edc2.usgs.gov/glcc/fao/data_specifications.php

dan atur CRS pengguna ini:

+proj=goode +lon_0=-100 +x_0=-11119487.43 +y_0=0 +a=6370997 +b=6370997 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs

untuk Amerika Utara.

lon_0 dan x_0 berasal dari Lampiran 03 dokumen ini: http://proceedings.esri.com/library/userconf/proc98/proceed/TO850/PAP844/P844.HTM

Melapisi batas dunia Bumi alami, saya mendapatkan gambar ini:

masukkan deskripsi gambar di sini

Cocok baik di Amerika Utara (sekitar 20 km salah tempat), dan terdistorsi di daerah lain.


EDIT

Hanya memperhatikan bahwa proj yang tahu menginterupsi goode homolosine sekarang:

+proj=igh +ellps=sphere +towgs84=0,0,0 +lon_0=100w +x_0=-11119487.43

memberi saya gambar batas dunia Bumi Alam dan kisi-kisi vektor 10 derajat:

masukkan deskripsi gambar di sini

Saya tidak yakin lagi parameter mana +proj=ighyang diperlukan.

AndreJ
sumber
Menggunakan definisi di atas saya mendapatkan pergeseran piksel. Bagi saya definisi itu +proj=igh +towgs84=0,0,0berfungsi dengan baik
Thomas
3

QGIS 2.18 hampir mendukung sepenuhnya Homolosine untuk raster. Anda hanya perlu mewaspadai masalah # 959 : melakukan raster warpping langsung dengan GDAL menggunakan CHECK_WITH_INVERT_PROJparameter dan Anda baik.

Mengenai vektor, semuanya masih rumit. Beberapa kesalahan yang ditunjukkan AndreJ dalam jawabannya menang. Penyebabnya sederhana: QGis tidak menyadari bahwa ini adalah program Kartografi dan menganggap co-domain dari setiap proyeksi sebagai ruang yang terbatas.

Ini memiliki perbaikan sederhana: memotong setiap layer vektor yang diimpor ke QGis dengan co-domain Homolosine. Beberapa waktu yang lalu saya membuat tersedia vektor co-domain diskrit yang didefinisikan pada interval satu derajat bujur dan lintang.

DOI

masukkan deskripsi gambar di sini

Luís de Sousa
sumber