如何基于指定天线坐标通过SPICE计算JUICE的方位角与仰角

要实现基于自定义天线坐标的计算，你需要补充地理坐标转惯性系位置、相对视线矢量计算、地平坐标系转换这几个核心步骤，以下是修改后的完整代码及关键说明：

完整代码 #

import os
import spiceypy as spice
from datetime import datetime, timezone
import math

# 加载元内核
meta_kernel_path = os.path.join('SPICE', 'Kernels', 'mk', 'juice_mk.tm')
spice.furnsh(meta_kernel_path) 

# （可选）如果元内核未包含地球PCK内核，手动加载（替换为实际路径）
# spice.furnsh(os.path.join('SPICE', 'Kernels', 'pck', 'pck00010.tpc'))

# 设置观测时间（UTC）
utc_now = datetime.now(timezone.utc).strftime('%Y-%m-%dT%H:%M:%S')
et = spice.str2et(utc_now)  # 转换为星历时间

# 天线地理参数：纬度、经度（度转弧度），海拔高度（米，此处设为0）
antenna_lat = math.radians(47.397490)
antenna_lon = math.radians(8.550440)
antenna_alt = 0.0  # 单位：米

# 1. 地理坐标转地球固定系（ITRF93）位置（单位：公里）
antenna_itrfs_pos, _ = spice.georec(antenna_lat, antenna_lon, antenna_alt/1000, 
                                    'ITRF93', 'EARTH')

# 2. 地球固定系转J2000惯性系
# 获取ITRF93到J2000的转换矩阵
pxform = spice.pxform('ITRF93', 'J2000', et)
# 矩阵乘法得到J2000系下天线相对于地球质心的位置
antenna_j2000_pos = spice.mxv(pxform, antenna_itrfs_pos)

# 3. 获取JUICE相对于地球质心的J2000位置（单位：公里）
juice_pos, light_time = spice.spkpos('JUICE', et, 'J2000', 'NONE', 'EARTH')

# 4. 计算JUICE相对于天线的视线矢量（J2000系）
rel_pos_j2000 = juice_pos - antenna_j2000_pos

# 5. 视线矢量转本地地平坐标系（东北天ENU）
# 将J2000系视线矢量转回ITRF93系
pxform_j2000_to_itrfs = spice.pxform('J2000', 'ITRF93', et)
rel_pos_itrfs = spice.mxv(pxform_j2000_to_itrfs, rel_pos_j2000)

# 构建ITRF93到ENU的转换矩阵
lat = antenna_lat
lon = antenna_lon
east = [-math.sin(lon), math.cos(lon), 0.0]
north = [-math.sin(lat)*math.cos(lon), -math.sin(lat)*math.sin(lon), math.cos(lat)]
zenith = [math.cos(lat)*math.cos(lon), math.cos(lat)*math.sin(lon), math.sin(lat)]
itrfs_to_enu = [east, north, zenith]

# 转换到ENU系
rel_pos_enu = spice.mxv(itrfs_to_enu, rel_pos_itrfs)

# 6. 计算方位角、仰角和距离
az, el, dist = spice.azlcyl(rel_pos_enu[0], rel_pos_enu[1], rel_pos_enu[2])
az_deg = math.degrees(az)
el_deg = math.degrees(el)

# 输出结果
print(f"JUICE相对于天线的方位角: {az_deg:.2f}° (正北顺时针旋转)")
print(f"JUICE相对于天线的仰角: {el_deg:.2f}° (地平线向上)")
print(f"直线距离: {dist:.2f} 公里")

# 卸载内核
spice.kclear()

关键步骤说明 #

1. 地理坐标转惯性系：
   • 用georec将经纬度、海拔转换为地球固定系（ITRF93）的位置矢量，再通过pxform获取的转换矩阵，将其转到J2000惯性系（与JUICE的位置参考系统一）。
2. 相对视线计算：
   • 用JUICE相对于地球质心的位置，减去天线相对于地球质心的位置，得到天线指向JUICE的视线矢量。
3. 地平坐标系转换：
   • 手动构建地球固定系到本地东北天（ENU）坐标系的转换矩阵，将视线矢量转换为本地坐标系下的分量，再用azlcyl直接计算方位角（正北顺时针）和仰角（地平线向上）。

注意事项 #

• 必须确保加载了地球PCK内核（包含地球形状、自转参数），如果元内核未包含，需手动加载pck00010.tpc这类标准内核文件。
• 天线海拔高度需替换为实际值（单位米，代码中已转换为公里）。

内容的提问来源于stack exchange，提问作者Akut Luna