Bagaimana saya bisa mengubah koordinat geomagnetik menjadi koordinat geografis tanpa melakukan perhitungan sendiri?

10

Saya punya beberapa data yang berada dalam koordinat geomagnetik; yaitu, mereka lintang dan bujur, tetapi mengacu pada kutub utara magnetik di Kanada daripada kutub utara geografis. Seperti yang Anda tebak, saya ingin mencocokkannya dengan data lain yang ada di sistem koordinat lainnya.

(Data berasal sebagai perhitungan oval auroral; yaitu, mereka berasal sebagai geomagnetik dan tidak ada sumber geografis yang dapat saya konsultasikan, bahkan pada prinsipnya.)

Data akan disimpan di PostGIS.

Cara mudah untuk melakukan ini adalah menemukan kode EPSG untuk sistem koordinat geomagnetik, menandai data dengan kode itu, dan kemudian semuanya diproyeksikan secara ajaib ke apa pun yang saya butuhkan pada saat itu. Namun, saya tidak dapat menemukan kode EPSG di spatialreference.org yang mengandung geomagnetik atau magnetik . Berbagai pencarian Google juga tidak menghasilkan apa-apa.

Untuk tingkat presisi yang saya butuhkan, sistem koordinat lintang / bujur (misalnya, WGS84) yang dimodifikasi untuk memiliki kutub utara yang berbeda sudah cukup baik. Jadi jika ada cara mudah untuk mendefinisikan sistem koordinat seperti itu dan memuatnya ke PostGIS, itu juga tidak masalah.

Klarifikasi

Sebagai contoh, misalkan alih-alih mengubah dari koordinat geomagnetik ke geografis, saya hanya ingin mengkonversi dari UTM 15N ke WGS84. Saya bisa menyimpan objek UTM ke PostGIS tanpa modifikasi, dan menandainya dengan kode EPSG 26915 . Kemudian, ketika saya kemudian membuat permintaan, saya bisa mengatakan kepada PostGIS, "tolong kembalikan hasil dalam kode EPSG 4326 ", dan PostGIS akan dikonversi ke WGS84 secara otomatis.

Demikian pula, saya ingin menandai data geomagnetik dengan beberapa sistem koordinat dan membiarkan PostGIS melakukan proyeksi ulang di belakang layar, tanpa harus memanggil reproject()fungsi secara manual.

Reid
sumber
1
Bisakah Anda menjelaskan bagaimana koordinat geomagnetik Anda diperoleh sejak awal? Saya pikir itu mungkin memegang kunci untuk solusi yang akurat: tampaknya mereka dihitung dari koordinat geografis. (Kalau tidak, bagaimana seseorang akan secara langsung mengukur garis lintang geomagnetik ?)
whuber
@whuber, pertanyaan diedit. Terimakasih atas klarifikasinya.
Reid

Jawaban:

7

Saya tidak yakin siapa pun open source mendukung geomagnetik ini menjadi geografis.

Tetapi jika Anda memiliki sedikit koordinat, tolong. coba ini ..
http://wdc.kugi.kyoto-u.ac.jp/igrf/gggm/index.html

NASA telah menerbitkan algoritme, Anda dapat mencoba http://idlastro.gsfc.nasa.gov/ftp/pro/astro/mag2geo.pro

vadivelan
sumber
1
Saya sudah membuka tiket untuk tim proj4. Mereka mungkin menambahkan dukungan untuk itu di versi yang lebih baru? trac.osgeo.org/proj/ticket/178
nickves
Seseorang mengirim kode geo2mag.pro IDL ke Python . Kode terbalik ini sangat mirip dengan mag2geo.pro. Berikut adalah kode port-Python .
Serge Stroobandt
2

Tautan yang disediakan oleh vadivelan memperjelas bahwa koordinat geomagnetik hanyalah rotasi bola sehingga (a) kutub utara melewati kutub selatan magnetik saat ini dan (b) meridian utama melewati kutub fisik dan magnet.

Sebagai contoh, IGRF-95 menggunakan posisi kutub utara magnet pada tahun 1995, di lintang 79,30 dan bujur -71,41 derajat.

Rotasi bola dapat dihitung dengan banyak cara , biasanya dengan mengkonversi ke koordinat Cartesian geosentris 3D, menerapkan transformasi linear (melalui perkalian matriks 3 x 3 atau perkalian angka empat), dan mengonversi kembali ke koordinat bola. Tetapi jika Anda tidak ingin memprogram ini, ada trik:SIG yang mendukung proyeksi miring dapat berputar dengan cara proyeksi ulang. Proyeksikan koordinat geomagnetik seolah-olah itu koordinat geografis, katakanlah dengan cara aspek polar dari proyeksi sterografis (dengan meridian pusat pada 0 derajat). Kemudian jangan memproyeksikan hasilnya menggunakan aspek miring dari proyeksi yang sama dengan garis lintang asal di 79,30 dan meridian tengah -71,41. Akhirnya, proyeksikan sekali lagi menggunakan proyeksi yang diinginkan untuk mendapatkan "proyeksi ulang magis."

whuber
sumber
Hmm. Apakah ada cara bagi saya untuk mendefinisikan sistem koordinat GEOS atau PROJ, atau hal lain yang akan diterima PostGIS, dan biarkan perpustakaan melakukan proyeksi ulang tanpa saya harus melakukannya secara eksplisit?
Reid
Itu tergantung pada apa yang Anda maksud dengan "secara eksplisit." Pada titik tertentu - sebaiknya sedini mungkin dalam proses - Anda harus secara fisik mengubah koordinat dari geomagnetik ke geografis. Anda dapat melakukan ini sekali menggunakan GIS yang sesuai (mis., QGIS kelihatannya akan berfungsi) atau utilitas proyeksi ulang. Anda tidak perlu tahu persamaan proyeksi ulang untuk ini jika Anda menggunakan trik yang saya uraikan.
Whuber
Biarkan saya mengklarifikasi pertanyaan untuk menjawab itu.
Reid