ver84tover90.m

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

% WGS-84直角坐标转换成PE-90直角坐标
clear;
input('\n\ninput vertical_84:');                                                % 输入WGS-84坐标系下的直角坐标值
X_84=input('\nX_84=');                                                          % X轴分量（米）
Y_84=input('\nY_84=');                                                          % Y轴分量（米）
Z_84=input('\nZ_84=');                                                          % Z轴分量（米）
X_90=1/(1+1.9e-6*1.9e-6)*[X_84+1.9e-6*(Y_84-2.5)];
Y_90=1/(1+1.9e-6*1.9e-6)*[-1.9e-6*X_84+(Y_84-2.5)];
Z_90=1/(1+1.9e-6*1.9e-6)*[(1+1.9e-6*1.9e-6)*Z_84];
fprintf('\nX_90=%.3f \n',X_90);                                          
fprintf('\nY_90=%.3f \n',Y_90);
fprintf('\nZ_90=%.3f \n',Z_90);
fprintf('\n\nDeltaX=%.3f \n',X_84-X_90);
fprintf('\nDeltaY=%.3f \n',Y_84-Y_90);
fprintf('\nDeltaZ=%.3f \n',Z_84-Z_90);