xyz地心坐标转经纬度高程

1.我在读取GNSS的o文件时候，发现有的站点没有经纬度和高程。就像下图一样，只给了近似坐标。（下面是例子，bjfs一般是个人都能找到经纬度和高程）

2.没办法，只能自己转换。（注意单位是km）

clc,clear;

% 单位是km，看好了

wgs84 = wgs84Ellipsoid('kilometer');


% 下面是假设的xyz坐标

x = -2148744.6242;

y = 4426641.5679;

z = 4044656.1585;


% 因为单位是km，所以都除1000

x = x/1000;

y = y/1000;

z = z/1000;


% 地心坐标系转成大地坐标系

% 我们o文件，读取没有经纬度，用这个转换进行了

[lat,lon,alt]=ecef2geodetic(wgs84,x,y,z);


alt = alt*1000;% 换算为m