如何用Astropy将ECI坐标转换至ITRF2014等特定ITRS坐标系？

Astropy里的ITRS类默认依赖IERS提供的最新坐标框架，但你可以手动指定具体的ITRF版本（比如ITRF2014），同时有几种方式验证实际使用的框架：

1. 手动指定ITRF2014框架 #

Astropy v4.0及以上版本支持直接通过frame参数指定ITRF版本。转换时明确指定即可：

from astropy.coordinates import GCRS, ITRS, EarthLocation
from astropy.time import Time
import astropy.units as u

# 示例GCRS坐标
gcrs_pos = GCRS(ra=10*u.deg, dec=20*u.deg, distance=1*u.au, obstime=Time('2024-01-01T12:00:00'))

# 转换到ITRF2014，需指定frame参数和观测时间、位置
itrf2014_pos = gcrs_pos.transform_to(
    ITRS(
        frame='itrf2014',
        obstime=gcrs_pos.obstime,
        location=EarthLocation(lon=0*u.deg, lat=0*u.deg, height=0*u.m)
    )
)

2. 确认当前使用的ITRF框架版本 #

方法一：查看ITRS实例的frame属性 #

直接打印转换后坐标的框架标识，确认是否为目标版本：

print(itrf2014_pos.frame.frame)
# 输出应为 'itrf2014'

方法二：验证关联的IERS数据 #

Astropy的ITRF转换依赖IERS提供的极移、岁差章动等参数，不同ITRF版本对应不同的IERS数据基准。可以查看框架关联的IERS表来确认：

print(itrf2014_pos.frame.iers_table)

输出的IERS表版本会对应到具体的ITRF框架（比如ITRF2014对应IERS 2019或更新的表数据）。

3. 注意事项 #

• 若不指定frame参数，Astropy默认使用IERS最新推荐的ITRF框架（如ITRF2020），因此必须手动指定frame='itrf2014'才能确保转换到目标坐标系。
• 老版本Astropy（v4.0以下）不支持直接指定ITRF版本，需升级到最新稳定版。
• 确保IERS数据是最新的：可以通过astropy.utils.iers.IERS_Auto.open()更新本地IERS数据，避免因数据过时导致框架不匹配。

内容的提问来源于stack exchange，提问作者Babul Akhtar