geo84tover84.m

本程序设计了常用的几个坐标系的转换

% WGS-84下地理坐标转换成直角坐标
clear;
input('\n\ninput geography_84:');              % 输入WGS-84坐标系下的地理坐标值
phi_84=input('\nphi_84=');                     % WGS-84下的纬度（度）
lamda_84=input('\nlamda_84=');                 % WGS-84下的经度（度）
H_84=input('\nH_84=');                         % WGS-84下的高度（米）
a=6378137;                                     % WGS-84椭球体长半轴
f=1/298.257223563;                             % WGS-84椭球体扁率
e=2*f-f^2;                                     % WGS-84椭球体第一偏心率
N=a/sqrt(1-e*(sin(phi_84/180*pi))^2);          % N为卯酉圈曲率半径
X_84=(N+H_84)*cos(phi_84/180*pi)*cos(lamda_84/180*pi);
Y_84=(N+H_84)*cos(phi_84/180*pi)*sin(lamda_84/180*pi);
Z_84=[N*(1-e)+H_84]*sin(phi_84/180*pi);
fprintf('\nX_84=%.3f \n',X_84);                                                   % 以小数形式输出，保留3位小数，且输出换行
fprintf('\nY_84=%.3f \n',Y_84);
fprintf('\nZ_84=%.3f \n',Z_84);