MATLAB程序北京54转换为WGS84坐标(GPS).doc

《MATLAB程序北京54转换为WGS84坐标(GPS).doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、%北京54转换为WGS84坐标（GPS）%X=.0300Y=.4500X=.0255;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%输入X值Y=.2860;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%输入X值x=X;y=Y-fix(Y/)*-;b=x*3600*180/pi/.4969;T=(cos(b*pi/180/3600))^2;Bf=b+(+(+(2350+22*T)*T)*T)*sin(b*pi/180/3600)*cos(b*pi/180

2、/3600)*3600*180/pi/10^10;Q=(cos(Bf*pi/180/3600))^2;Nf=.902-[21562.267-(108.973-0.612*Q)*Q]*Q;Z=y/(Nf*cos(Bf*pi/180/3600));b2=(0.5+0.*Q)*sin(Bf*pi/180/3600)*cos(Bf*pi/180/3600);b3=0.-(0.-0.*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*Z^2)

3、*Z^2]*Z^2*b2*180*3600/pi;L=[1-(b3-b5*Z^2)*Z^2]*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/)*3+L1;fprintf('')fprintf('p点的WGS84坐标:%ft%ft%ft',B1,B2,B3)fprintf('p点的WGS84坐标:%

4、ft%ft%ft',L4,L2,L3)