我正在尝试在地理坐标和地磁坐标之间进行转换.我发现了以下Prolog脚本,但我不太了解它自己转换它.目标语言是Java,但任何可以理解的都很好(C,Python,VB,等等).

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

如果有人可以帮助转换这个脚本或解释它究竟在做什么(这些数组操作让我感到困惑),我真的很感激.

谢谢

解决方法:

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

在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))/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')