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.
Jawaban:
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
sumber
geo2mag.proIDL ke Python . Kode terbalik ini sangat mirip denganmag2geo.pro. Berikut adalah kode port-Python .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."
sumber