Bagaimana cara membuat peta Azimuthal Equidistant untuk radio jarak jauh?

10

Saya mencoba untuk menghasilkan peta skala dunia atau benua yang menggunakan proyeksi Azimuthal Equidistant. Saya seorang amatir radio, dan untuk antena, seluruh dunia tampak seperti menyebar di bidang datar, dengan antena di pusat dan seluruh dunia di (r, Θ) dari titik itu.

Ada perangkat lunak untuk melakukan ini - AZ_PROJ - tetapi menggunakan database sendiri, dan sepenuhnya ditulis dalam PostScript, sehingga tidak bisa dibajak dengan alat geo biasa.

OGR memiliki proyeksi 'aeqd', yang dapat dipanggil seperti ini:

ogr2ogr -t_srs "+proj=aeqd  +lat_0=43.7308 +lon_0=-79.2647" out.shp in.shp

untuk sebuah pusat di 43.7308 N, 79.2647 W. Sayangnya, ini membungkus segalanya di luar belahan bumi saat ini "di belakang" peta. Ini juga memproyeksikan koordinat ke dalam meter kartesius, di mana saya lebih suka kutub.

Ada saran untuk alternatif, tolong?

coordinate-system cartography scruss
sumber
Ini mungkin jawaban parsial, jika ada yang mencari: Proyeksi Ortho menghasilkan artefak
scruss

Jawaban:

2

Anda harus menentukan jari-jari bola, karena proj tampaknya hanya mendukung rumus bola dari proyeksi ini:

+proj=aeqd  +R=6371000 +lat_0=51 +lon_0=7
AndreJ
sumber
1

Saya tidak mengetahui adanya alat proyeksi yang menghasilkan koordinat polar.

Anda dapat mencoba menggunakan Mesin Proyeksi Esri dll (pe.dll) yang dibundel dengan ArcGIS Explorer yang dapat diunduh secara bebas. Ini adalah dll gaya C yang memiliki metode untuk menentukan jarak azimuth dan geodesik antara dua titik. Untuk detail lebih lanjut lihat posting blog Richie Carmichael .

Kirk Kuykendall
sumber
Terima kasih. Saya sudah bisa menghitung jarak azimuth dan geodesik dengan geod , yang merupakan bagian dari PROJ.4 : echo 'lat1 long1 lat2 long2' | geod + ellps = WGS84 -f "% .3f" -p -I + unit = m
scruss
1
Anda dapat palsu koordinat polar: setelah melakukan (oblique) azimut berjarak sama proyeksi, hanya menyatakan bahwa itu benar-benar hasil dari proyeksi azimut berjarak sama polar, maka unproject itu. Peta yang dihasilkan akan terbalik, tetapi refleksi dan terjemahan vertikal (transformasi affine sederhana) akan menyelesaikannya.
Whuber