《MATLAB程序北京54转换为WGS84坐标(GPS)》由会员分享，可在线阅读，更多相关《MATLAB程序北京54转换为WGS84坐标(GPS)（1页珍藏版）》请在金锄头文库上搜索。
1、%北京54转换为WGS84坐标（GPS）%X=3459174.0300 Y=36503163.4500X=3459181.0255; %输入X值Y=36503206.2860; %输入X值x=X;y=Y-fix(Y/1000000)*1000000-500000;b=x*3600*180/pi/6367558.4969;T=(cos(b*pi/180/3600)2;Bf=b+(50221746+(293622+(2350+22*T)*T)*T)*sin(b*pi/180/3600)*cos(b*pi/180/3600)*3600*180/pi/1010;Q=(cos(Bf*pi/180/360
2、0)2;Nf=6399698.902-21562.267-(108.973-0.612*Q)*Q*Q;Z=y/(Nf*cos(Bf*pi/180/3600);b2=(0.5+0.003369*Q)*sin(Bf*pi/180/3600)*cos(Bf*pi/180/3600);b3=0.333333-(0.166667-0.001123*Q)*Q;b4=0.25+(0.16161+0.00562*Q)*Q;b5=0.2-(0.1667-0.0088*Q)*Q;B=Bf-1-(b4-0.12*Z2)*Z2*Z2*b2*180*3600/pi;L=1-(b3-b5*Z2)*Z2*Z*180*3600/pi;B1=fix(B/3600);B2=fix(B-B1*3600)/60);B3=B-B1*3600-B2*60;L1=fix(L/3600);L2=fix(L-L1*3600)/60);L3=L-L1*3600-L2*60;L4=fix(Y/1000000)*3+L1;fprintf(nn)fprintf(p点的WGS84坐标:%ft%ft%ftn,B1,B2,B3) fprintf(p点的WGS84坐标:%ft%ft%ftn,L4,L2,L3)
