高斯克呂格投影

在網(wǎng)上搜集整理有關(guān)高斯克呂格投影的代碼并用 python 改寫。

from math import * p = 180.0 / pi * 3600def deg2rad(deg):dms = deg*3600/preturn dms def gauss_zhengsuan(B,L,L0): # 緯度、經(jīng)度、中央經(jīng)線經(jīng)度l = L-L0 l = deg2rad(l) #經(jīng)度B = deg2rad(B) #緯度a = 6378245.0 # Beijing1954 Krasovsky_1940',6378245.0,298.3f = 1.0/298.3b = a - a*fc = a**2/be = sqrt(a**2 - b**2)/ae1 = sqrt(a**2 - b**2)/bp = 3600*180/piBeta0 = 1 - (3.0/4)*e1**2 + (45.0/64)*e1**4 - (175.0/256)*e1**6 + (11025.0/16384)*e1**8Beta2 = Beta0-1Beta4 = (15.0/32)*e1**4 - (175.0/384)*e1**6 + (3675.0/8192)*e1**8Beta6 = -(35.0/96)*e1**6 + (735.0/2048)*e1**8Beta8 = (315.0/1024)*e1**8Xb0 = c*(Beta0*(B) + sin(B)*(Beta2*cos(B) + Beta4*cos(B)**3 + Beta6*cos(B)**5 + Beta8*cos(B)**7))Np = a/sqrt(1 - (e*sin(B))**2)m0 = l*cos(B)t = tan(B)in2 = (e1)**2*cos(B)**2X = Xb0 + (1.0/2)*Np*t*m0**2 + (1.0/24)*(5-t**2 + 9.0*in2 + 4.0*in2**2)*Np*t*m0**4 + (1.0/720)*(61 - 58.0*t**2)*Np*t*m0**6.0Y = Np*m0 + (1.0/6)*(1 - t**2 + in2)*Np*m0**3 + (1.0/120)*(5 - 18.0*t**2 + t**4 + 14.0*in2 - 58.0*t**2*in2)*Np*m0**5 + 500000return Y,X # X表示縱坐標(biāo)，Y表示橫坐標(biāo)result = gauss_zhengsuan(32.34, 120.123, 120)

總結(jié)

以上是生活随笔為你收集整理的高斯克吕格投影 python2.x 版本的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。