Memanipulasi Proyeksi Azimuthal Equidistant di QGIS

8

Saya memiliki bentuk negara-negara dunia yang diproyeksikan sebagai Azimuthal Equidistant dengan Chicago, AS sebagai pusatnya. Saya ingin mengubah proyeksi sehingga Edinburgh, UK adalah pusatnya. Orang akan berpikir bahwa ini harus sesederhana mengubah 2 angka dalam parameter sistem koordinat ...

Masalah yang sepertinya saya alami adalah bahwa QGIS menyulitkan atau bahkan mustahil untuk melihat kode proj4 lengkap dari parameter sistem koordinat shapefile yang diberikan.

Inilah yang saya coba sejauh ini:

  • Ditambahkan shapefile dari negara-negara dunia yang diproyeksikan sebagai Azimuthal Equidistant dengan Chicago sebagai pusatnya
  • Properti yang dibuka, tab metadata, salin kode parameter proj4
  • Pergi ke pengaturan, proyeksi kustom, tempel dalam kode proj4, mengubah lat / lon dari yang Chicago ke yang Edinburgh
  • Shapefile tersimpan dari negara-negara dunia, menetapkan proyeksi kustom baru saya, menambahkan file baru ke peta
  • File baru saya terlihat seperti peta Pangea, dengan semua benua bersama-sama.

Saya pikir harus ada parameter lain selain yang ditunjukkan pada kode proj4 dari tab metadata. Satu-satunya hal yang saya ubah dalam kode proj4 adalah pusat peta, tetapi perubahan yang saya lihat jauh lebih dari itu. Pengamatan saya menunjukkan bahwa lebih dari satu variabel independen berubah.

Saya akan senang mengirim shapefile saya kepada siapa saja yang ingin bermain dengannya.

Terima kasih,

Daniel Wolf, Penggemar Lingkungan / Geospasial

Dan Wolf
sumber

Jawaban:

8

Menggunakan string proj sederhana:

+proj=aeqd +lat_0=51 +lon_0=7

Saya mendapatkan gambar ini dari batas dunia Bumi Alami, dan grid 10 derajat yang harus membentang dari -150 ° hingga 150 ° E:

aeqd pada ellipsoid WGS84

Sepertinya segala sesuatu di luar lon_0 +/- 90 ° akan terpotong, dan bagian-bagian belahan bumi selatan juga tertekuk.


EDIT

Menurut tiket ini: http://trac.osgeo.org/proj/ticket/211 proj memiliki masalah dengan memproyeksi ulang data aeqd pada ellipsoid.

Solusi sederhana: gunakan bola (dan bukan ellipsoid default) dengan mendefinisikan radius:

+proj=aeqd  +R=6371000 +lat_0=51 +lon_0=7

Ini memberi saya gambar berikut:

aeqd di sebuah bola

Garis kisi dipadatkan untuk menunjukkan kurva yang bagus.

AndreJ
sumber