clc;
% !!! 此代码遵循 WGS84 模型 !!!; % 定义纬度为 50 的 X、Y、Z 样本; 经度 100; 高度200M; % http://www.apsalin.com/convert-geodetic-to-cartesian.aspx X = -713345.437320888; Y = +4045583.0097852; Z = +4862942.24652593;
%估计曲率半径; 一 = 6378137; % 仪表; % 半长轴; f = 1/298.257223563; % 椭球展平; b = a *(1-f); % 定义半短轴;
%估算辅助值; P = sqrt(X ^ 2 + Y ^ 2); Theta = atan(Z * a / P * b); e = sqrt(((a^2) - (b^2))/a^2); %地球的第一偏心率; e2 = sq
2021-09-22 21:00:38
2KB
matlab
1