Kendimi matematik yapmadan jeomanyetik koordinatları coğrafi koordinatlara nasıl dönüştürebilirim?


10

Jeomanyetik koordinatlarda verilerim var; yani, enlem ve boylamdır, ancak coğrafi kuzey kutbundan ziyade Kanada'daki manyetik kuzey kutbuna referansla . Tahmin edebileceğiniz gibi, bunları diğer koordinat sistemlerindeki diğer verilerle eşleştirmek istiyorum.

(Veriler auroral ovalin bir hesaplaması olarak ortaya çıkmıştır; yani, jeomanyetik olarak ortaya çıkmıştır ve prensipte bile başvurabileceğim bir coğrafi kaynak yoktur.)

Veriler PostGIS'de saklanacaktır.

Bunu yapmanın kolay yolu, jeomanyetik koordinat sistemi için bir EPSG kodu bulmak, verileri bu kodla etiketlemektir ve o zaman her şey sihirli bir şekilde o zaman ihtiyacım olan her şeye yeniden yönlendirilir. Ancak, spatialreference.org üzerinde jeomanyetik veya manyetik içeren herhangi bir EPSG kodu bulamıyorum . Çeşitli Google aramaları da hiçbir şey yapmaz.

İhtiyacım olan hassasiyet seviyesi için, farklı bir kuzey kutbuna sahip olacak şekilde değiştirilmiş herhangi bir enlem / boylam koordinat sistemi (örn. WGS84) yeterince iyidir. Dolayısıyla, böyle bir koordinat sistemi tanımlamanın ve PostGIS'e yüklemenin kolay bir yolu varsa, bu da iyi olur.

açıklama

Örneğin, jeomanyetik koordinatlardan coğrafi konuma dönüştürmek yerine, UTM 15N'den WGS84'e dönüştürmek istediğimi varsayalım. UTM nesnesini değiştirmeden PostGIS'e kaydedebilir ve EPSG kodu 26915 ile etiketleyebilirim . Daha sonra bir sorgu yapılmış Sonra, ben "EPSG kodunda sonuçlarını iade edin, PostGIS için söyleyebiliriz 4326 " ve PostGIS otomatik WGS84 dönüştürmek istiyorum.

Benzer şekilde, bazı koordinat sistemi ile jeomanyetik verileri etiketlemek ve PostGIS'in bir reproject()işlevi manuel olarak çağırmak zorunda kalmadan sahnelerin arkasındaki yeniden projeksiyonu yapmasına izin vermek istiyorum .


1
İlk etapta jeomanyetik koordinatlarınızın nasıl elde edildiğini açıklayabilir misiniz? Bunun doğru bir çözümün anahtarını tutabileceğini düşünüyorum: Görünüşe göre coğrafi koordinatlardan hesaplanıyorlar. (Aksi takdirde, doğrudan bir jeomanyetik enlem nasıl ölçülür ?)
whuber

@whuber, soru düzenlendi. Açıklama için teşekkürler.
Reid

Yanıtlar:


7

Herhangi bir açık kaynağın bu jeomanyetiği coğrafi olarak desteklediğinden emin değilim.

Ancak birkaç koordinatınız varsa, lütfen. bunu dene ..
http://wdc.kugi.kyoto-u.ac.jp/igrf/gggm/index.html

NASA algoritmayı yayınladı, http://idlastro.gsfc.nasa.gov/ftp/pro/astro/mag2geo.pro deneyebilirsiniz


1
Proj4 ekibi için bir bilet açtım. Daha sonraki bir sürümde destek ekleyebilirler mi? trac.osgeo.org/proj/ticket/178
nickves

Birisi geo2mag.pro IDL kodunu Python'a taşıdı . Bu ters kod çok benzer mag2geo.pro. İşte Python tarafından taşınan kod .
Serge Stroobandt

2

Vadivelan marka tarafından sağlanan bağlantılar o temizlemek jeomanyetik koordinatlar kürenin sadece bir rotasyon vardır ki (a) kuzey kutbu akım manyetik geçer böylece güney meridyeni fiziksel ve manyetik kutuplar geçer kutbu ve (b).

Örneğin IGRF-95, 1995'te manyetik kuzey kutbunun pozisyonunu 79.30 enleminde ve -71.41 derece boylamında kullanır.

Küresel dönüşler , tipik olarak 3B jeosantrik Kartezyen koordinatlarına dönüştürülerek, doğrusal bir dönüşüm uygulanarak (3 x 3 matris çarpımı veya bir kuaterniyon çarpımı yoluyla) ve küresel koordinatlara geri dönüştürülerek hesaplanabilir . Ama bunu programlamak istemiyorsanız , bir hile var:eğik çıkıntıları destekleyen bir CBS, bir yeniden projeksiyon yoluyla dönebilir. Jeomanyetik koordinatları, örneğin sterografik izdüşümün (0 derecede merkezi bir meridyen ile) kutupsal bir yönü aracılığıyla coğrafi koordinatlar gibi yansıtın. Ardından 79.30'da orijin enlemi ve -71.41'in merkezi meridyeni ile aynı projeksiyonun eğik yönünü kullanarak sonucu çıkarın. Son olarak, "büyülü yeniden projeksiyon" elde etmek için istenen herhangi bir projeksiyon kullanarak bir kez daha yansıtın.


Hmm. Böyle bir GEOS veya PROJ koordinat sistemi veya PostGIS'in kabul edeceği başka bir şey tanımlamamın bir yolu var mı?
Reid

"Açıkça" ile ne demek istediğinize bağlıdır. Bir noktada - tercihen mümkün olduğunca erken - koordinatları fiziksel olarak jeomanyetikten coğrafi konuma dönüştürmelisiniz. Bunu bir kez uygun bir CBS (örneğin, QGIS işe yarayacak gibi) veya yeniden projelendirme yardımcı programını kullanarak yapabilirsiniz. Belirttiğim hileyi kullanırsanız, bunun için yeniden tahmin denklemlerini bilmenize gerek yoktur.
whuber

Bunu cevaplamak için soruyu netleştireyim.
Reid
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.