MATLAB实现大地坐标系统和空间直角坐标系统互转写代码

在MATLAB中，实现大地坐标系统（如WGS84经纬度）到空间直角坐标系（如XYZ笛卡尔坐标）的转换，以及反之，通常涉及到地理坐标到地球球面坐标系的转换，然后进一步转换到三维空间坐标。以下是基本的步骤和示例代码：

1. **从大地坐标（经度、纬度、海拔高度）到XYZ**:
使用`geodetic2cartesian`函数，它将地理坐标转换为大地坐标系下XYZ坐标。

function [X,Y,Z] = latlonToXYZ(lat, lon, altitude)
    ellipsoid = 'WGS84'; % 假设使用的是WGS84椭球体
    XYZ = geodetic2cartesian(lat, lon, altitude, ellipsoid);
    X = XYZ(1);
    Y = XYZ(2);
    Z = XYZ(3) + altitude; % 高度通常包括地形高度
end

2. **从XYZ到大地坐标**:
可以使用`cartesian2geodetic`函数反向操作。

function [lat, lon] = XYZtoLatLon(X, Y, Z, altitude, ellipsoid)
    latlon = cartesian2geodetic(X, Y, Z - altitude, ellipsoid); % 减去海拔高度
    lat = latlon(1);
    lon = latlon(2);
end

请注意，以上代码假设了地球是一个完美的球体（即使用了`ellipsoid`选项），但在实际应用中，需要考虑地球的扁平化影响。