java - 从地理坐标到地磁坐标的转换

Question

我正在尝试在地理坐标和地磁坐标之间进行转换。我找到了以下 Prolog 脚本，但我对它的理解还不够，无法自行转换。目标语言是 Java，但任何可以理解的语言都可以(C、Python、VB 等等)。

http://idlastro.gsfc.nasa.gov/ftp/pro/astro/geo2mag.pro

如果有人可以帮助转换此脚本或解释它到底在做什么(那些数组操作让我感到困惑)，我将不胜感激。

谢谢

Answer 1

根据应用，高度可能是此坐标转换中的重要变量，因为地磁坐标是地球偶极子磁场的映射。

在 Python 中，您可以使用 SpacePy ( http://sourceforge.net/projects/spacepy/ ) 轻松地将地理坐标转换为地磁(反之亦然)。

由于您正在寻找转换为 Java 的源代码，SpacePy 正在实现 Fortran 国际辐射带环境建模 (IRBEM) 库，其源代码可用 (http://irbem.svn.sourceforge.net/viewvc/irbem/web/index.html)

在 Python 中，如果其他人正在寻找快速解决方案:

import spacepy.coordinates as coord
from spacepy.time import Ticktock
import numpy as np
def geotomag(alt,lat,lon):
    #call with altitude in kilometers and lat/lon in degrees 
    Re=6371.0 #mean Earth radius in kilometers
    #setup the geographic coordinate object with altitude in earth radii 
    cvals = coord.Coords([np.float(alt+Re)/Re, np.float(lat), np.float(lon)], 'GEO', 'sph',['Re','deg','deg'])
    #set time epoch for coordinates:
    cvals.ticks=Ticktock(['2012-01-01T12:00:00'], 'ISO')
    #return the magnetic coords in the same units as the geographic:
    return cvals.convert('MAG','sph')