수학을 직접하지 않고 어떻게 지자기 좌표를 지리적 좌표로 변환 할 수 있습니까?


10

지자기 좌표에있는 데이터가 있습니다. 즉, 위도 및 경도이지만 지리적 북극이 아닌 캐나다자기 북극 을 기준으로합니다 . 짐작할 수 있듯이 다른 좌표계에있는 다른 데이터와 일치시키고 싶습니다.

(데이터는 오로라 타원의 계산으로 시작되었습니다. 즉, 그것들은 지자기에서 시작되었으며 원칙적으로 상담 할 수있는 지리학 적 소스가 없습니다.)

데이터는 PostGIS에 저장됩니다.

이 작업을 수행하는 쉬운 방법은 지자기 좌표계의 EPSG 코드를 찾고 해당 코드로 데이터에 태그를 지정한 다음 모든 것이 당시에 필요한대로 마술처럼 재 투영됩니다. 그러나 geomagnetic 또는 magnetic을 포함하는 patialreference.org에서 EPSG 코드를 찾을 수 없습니다 . 다양한 Google 검색에서도 아무것도 나타나지 않습니다.

필요한 정밀도 수준을 위해, 다른 북극을 갖도록 수정 된 위도 / 경도 좌표계 (예 : WGS84)로 충분합니다. 따라서 이러한 좌표계를 정의하고 PostGIS에로드하는 쉬운 방법이 있다면 괜찮습니다.

설명

예를 들어, 지자기 좌표에서 지리적으로 변환하는 대신 UTM 15N에서 WGS84로 변환하려고한다고 가정합니다. UTM 객체를 수정되지 않은 PostGIS에 저장하고 EPSG 코드 26915로 태그를 지정할 수 있습니다. 그런 다음 나중에 쿼리 할 때 PostGIS에 "EPSG 코드 4326으로 결과를 반환하십시오"라고 말하면 PostGIS가 자동으로 WGS84로 변환됩니다.

마찬가지로 지자기 데이터에 좌표 시스템으로 태그를 지정하고 PostGIS가 reproject()수동으로 함수 를 호출하지 않고도 장면 뒤에서 재 투영을 수행하도록하고 싶습니다 .


1
지자기 좌표를 어떻게 얻었는지 설명해 주시겠습니까? 나는 그것이 정확한 해결책의 열쇠를 가질 수 있다고 생각한다. 그것이 지리적 좌표로부터 계산 된 것으로 보인다. (그렇지 않으면, 어떻게 지자기 위도를 직접 측정 할 것인가?)
whuber

@ whuber, 질문이 편집되었습니다. 설명 주셔서 감사합니다.
Reid

답변:


7

나는 어떤 오픈 소스가이 지자기에서 지리를 지원하는지 확신하지 못한다.

그러나 좌표가 거의 없다면 pls. 이것을 시도하십시오.
http://wdc.kugi.kyoto-u.ac.jp/igrf/gggm/index.html

NASA는 알고리즘을 발표했습니다. http://idlastro.gsfc.nasa.gov/ftp/pro/astro/mag2geo.pro


1
proj4 팀 티켓을 열었습니다. 이후 버전에서 지원을 추가 할 수 있습니까? trac.osgeo.org/proj/ticket/178
nickves

누군가 geo2mag.pro IDL 코드를 Python으로 이식했습니다 . 이 역 코드는와 매우 유사합니다 mag2geo.pro. 다음은 파이썬으로 작성된 코드 입니다.
Sroo Stroobandt

2

vadivelan 메이크업에 의해 제공되는 링크는 그 취소 지자기 좌표가 영역의 단지 회전이다 는 (a) 북쪽 극은 현재의 자기를 통과하므로 남쪽 본초 자오선은 물리적 자극을 통과 극와 (b).

예를 들어, IGRF-95는 1995 년의 자북 위치를 위도 79.30 및 경도 -71.41 도로 사용합니다.

구면 회전은 일반적으로 3D 지오 센 트릭 직교 좌표로 변환하고, 선형 변환 (3 x 3 행렬 곱셈 또는 쿼터니언 곱셈을 통해)을 적용하고 구면 좌표로 다시 변환하여 여러 가지 방법으로 계산할 수 있습니다 . 그러나 이것을 프로그래밍하고 싶지 않다면 트릭이 있습니다.경사 투사를 지원하는 GIS는 재 투영을 통해 회전 할 수 있습니다. 지구 자기 좌표를 마치 입체 좌표 투영의 극점 (0 도의 중앙 자오선 포함)을 사용하여 지리적 좌표 인 것처럼 투영하면됩니다. 그런 다음 79.30에 원점 위도 및 중앙 자오선 -71.41을 가진 동일한 투영의 경사면을 사용하여 결과를 투영 해제하십시오. 마지막으로 "매직 재 투영"을 얻기 위해 원하는 프로젝션을 사용하여 한 번 더 투사하십시오.


흠. 이러한 GEOS 또는 PROJ 좌표계 또는 PostGIS가 수용 할 다른 것을 정의하고 라이브러리가 명시 적으로 수행하지 않고 라이브러리를 재 투영하도록하는 방법이 있습니까?
Reid

"명시 적으로"의 의미에 따라 다릅니다. 어느 시점에서 (가능하면 가능한 빨리 프로세스 초기에) 좌표를 지자기에서 지리적으로 물리적으로 변환해야합니다. 적절한 GIS (예 : QGIS가 작동하는 것처럼 보임) 또는 재 투영 유틸리티를 사용하여이 작업을 한 번 수행 할 수 있습니다. 내가 설명한 트릭을 사용하는 경우 재 투영 방정식을 알 필요가 없습니다.
whuber

그 대답에 대한 질문을 명확히하겠습니다.
Reid
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.