寻求将[经度,纬度,海拔]转换为以地面站为原点的[X,Y,Z]坐标系的方案

首先需要注意：你提供的示例数据存在格式错误——经纬度中的点是千分位分隔符（而非小数点），且部分数值超出合理范围（比如飞机经度原数据34.890.851.772明显不符合WGS84经度-180~180°的范围），以下方案会基于修正后的合理数据展开。

Python 实现（基于pyproj库） #

pyproj是Python处理地理坐标转换的常用工具，支持WGS84与ECEF（地心地固坐标系）的直接转换。

步骤与代码 #

1. 安装依赖：

pip install pyproj

2. 坐标转换代码：

from pyproj import Transformer

# 修正后的示例数据（十进制度数+米级海拔）
ground_lat = 41.947694  # 地面站纬度
ground_lon = 3.209083   # 地面站经度
ground_alt = 379.41     # 地面站海拔

plane_lat = 41.9763015750  # 飞机纬度
plane_lon = 3.4890851772   # 修正后的飞机经度
plane_alt = 971.32         # 飞机海拔

# 定义WGS84经纬度 → ECEF坐标系的转换器
transformer = Transformer.from_crs("EPSG:4326", "EPSG:4978", always_xy=True)

# 转换地面站到ECEF
ground_x, ground_y, ground_z = transformer.transform(ground_lon, ground_lat, ground_alt)

# 转换飞机到ECEF
plane_x, plane_y, plane_z = transformer.transform(plane_lon, plane_lat, plane_alt)

# 计算以地面站为原点的相对XYZ坐标
relative_x = plane_x - ground_x
relative_y = plane_y - ground_y
relative_z = plane_z - ground_z

print(f"飞机相对地面站的XYZ坐标：({relative_x:.3f}, {relative_y:.3f}, {relative_z:.3f})")

说明：

• EPSG:4326是WGS84经纬度坐标系，EPSG:4978是WGS84对应的ECEF坐标系
• always_xy=True确保输入顺序为经度、纬度，符合GIS工具的通用规范

Matlab 实现 #

Matlab内置geodetic2ecef函数，无需额外安装工具即可完成转换。

代码示例 #

% 修正后的示例数据
ground_lat = 41.947694;  % 地面站纬度（度）
ground_lon = 3.209083;   % 地面站经度（度）
ground_alt = 379.41;     % 地面站海拔（米）

plane_lat = 41.9763015750;  % 飞机纬度
plane_lon = 3.4890851772;   % 修正后的飞机经度
plane_alt = 971.32;         % 飞机海拔

% 转换为ECEF坐标系（默认使用WGS84椭球参数）
[ground_x, ground_y, ground_z] = geodetic2ecef(ground_lat, ground_lon, ground_alt);
[plane_x, plane_y, plane_z] = geodetic2ecef(plane_lat, plane_lon, plane_alt);

% 计算相对坐标
relative_x = plane_x - ground_x;
relative_y = plane_y - ground_y;
relative_z = plane_z - ground_z;

fprintf('飞机相对地面站的XYZ坐标：(%.3f, %.3f, %.3f)\n', relative_x, relative_y, relative_z);

关键注意事项 #

• 经纬度必须是十进制度数，若你拿到的是度分秒格式，需先转换为十进制
• 转换后的相对XYZ坐标基于ECEF坐标系，完全满足你3D椭球模拟的需求
• 若需要局部坐标系（如东北天ENU），可基于ECEF相对坐标进一步转换，但当前场景无需额外操作

内容的提问来源于stack exchange，提问作者Leon Archinger