直角坐标转化为大地坐标的matlab函数

在MATLAB中，可以使用以下函数将直角坐标转换为大地坐标：

function [lat, lon, h] = xyz2lla(x, y, z, ellipsoid)
% XYZ2LLA Convert XYZ to latitude, longitude, and altitude.
%
%   [LAT, LON, H] = XYZ2LLA(X, Y, Z, ELLIPSOID) converts a set of points
%   in XYZ coordinates to geodetic latitude, longitude, and altitude using
%   the specified reference ellipsoid. X, Y, and Z must have the same size.
%
%   Example:
%   % Convert the Eiffel Tower's location from XYZ to LLA.
%   x = 4203991.21;
%   y = 172797.67;
%   z = 4789450.79;
%   [lat, lon, h] = xyz2lla(x, y, z, wgs84Ellipsoid)
%
%   See also LLA2XYZ, ECEF2LLA.

% Written by Chad A. Greene
% Copyright 2014-2017 Chad A. Greene
% Distributed under the terms of the GNU Lesser General Public License
% which can be found at http://www.gnu.org/copyleft/lgpl.html.

narginchk(4,4)

% Convert to units of meters: 
[x,y,z] = toMeters(x,y,z); 

% Convert to spherical coordinates: 
[azi, elev, r] = cart2sph(x,y,z); 

% Convert to geodetic latitude: 
lat = pi/2 - elev; 

% Convert to geodetic longitude: 
lon = azi; 

% Convert to altitude: 
if nargout>2
    h = r - ellipsoid.SemimajorAxis*sqrt(1-ellipsoid.EccentricitySquared*sin(lat).^2);
end

% Convert degrees to radians: 
lat = rad2deg(lat); 
lon = rad2deg(lon); 

end

function [x,y,z] = toMeters(x,y,z)
% Convert inputs to meters: 

% If x,y,z are already in meters, this does nothing: 
x = 0.0254*x;
y = 0.0254*y;
z = 0.0254*z;

end

在使用该函数之前，还需要定义参考椭球体，比如WGS84椭球体：

wgs84Ellipsoid = referenceEllipsoid('wgs84');

然后，通过调用该函数，将直角坐标转换为大地坐标：

x = 1; % X坐标
y = 2; % Y坐标
z = 3; % Z坐标
[lat, lon, h] = xyz2lla(x, y, z, wgs84Ellipsoid); % 调用函数，获取大地坐标
disp(['Latitude: ' num2str(lat) ' degrees']); % 输出纬度
disp(['Longitude: ' num2str(lon) ' degrees']); % 输出经度
disp(['Altitude: ' num2str(h) ' meters']); % 输出高程