[转载] :经纬度与WGS84坐标转换本代码实现了WGS84系统中大地坐标(BLH )与空间直角坐标(XYZ )的相互转换，符合标准语法，可以直接使用
WGS 84:-2175790.739698914461032.112077343992337.79032463
blh :38.99999998116.00000000000033.000000069718808
模制颜色转移！ //WGS84系统BLH坐标和空间直角坐标转换！ //Fortran Coder http://fcode.cn！//自适应框架' sccodeimplicitnoneinteger， parameter :3360 DP=selected _ real _ kind (12 ) real (kind=DP ) ) ) ) ) 652 parameter 33603360 pi=3.1440 private 3360: a=6378137._ DP real private :3360 f=1.0 _ DP/298.25723563 _ DP real (kind=DP (，parameter )， private 3360: t=a* (1.0 _ DP-f private :3360 e=sqrt (a* a-t * t )/(a* a ) ) containssubroutineXYZ2blh(x，x ) intent ) out ) : B、l、 hreal(kind=DP ) :3360 n integer : il=atan (ABS (y/x ) ) (b=ATAN ) ABS ) z/sqrt (x * xy * y ) ) ) ) ) 1.0_DP-e*e ) b=atan ) ) z ) ) n无言enddoif(x0._DP.and.y0._DP ) L=pi-Lif ) x0._ DP.and.y0 180._ DP/pil=l * 180._ DP/piendsubroutinexyz2blhsubroutineblh2XYZ (b，l，h，x，y，z ) real(kind=DP )，ind zreal ) kind=DP ) : N，Br， lrbr=b * pi/180._ dplr=l * pi/180._ dpn=a/sqrt (1.0 _ DP-e * e * sin(br ) sin(br ) br ) ) n ) x=(NH ) ) ) 652 endsubroutineblh2xyzendmodulecorrtransprogramwww _ fcode _ cnusecode ) h=33._DPcallblH2XYZ(b，l，h，x，y，z ) write(*，* ) ' WGS84: '，x，y，z call XYZ2BLH(x ) x，y，z，b
posted on 2015-07-31 17:13北读(…)评论(…)编辑收藏转载于：3https://www.cn blogs.com/China 3s/p/4692690.html