Python的經(jīng)緯度與xy坐標(biāo)系相互轉(zhuǎn)換方式
經(jīng)緯度與xy坐標(biāo)系相互轉(zhuǎn)換
代碼是已經(jīng)運行過的,可以直接輸入?yún)?shù)使用
注意事項
- 經(jīng)緯度參數(shù)不可以寫顛倒,否則會由于math.log參數(shù)不能為負(fù)數(shù)和0的原因報錯
- 由于地球半徑的不確定性,所以(x,y)與(lon, lat)相互轉(zhuǎn)換會有誤差
奉上代碼及運行結(jié)果
- 經(jīng)緯度轉(zhuǎn)換成xy坐標(biāo)
實例經(jīng)維度參數(shù) (114.5456282282352,36.631262731204049)
def millerToXY (lon, lat): xy_coordinate = [] #地球周長 L = 6381372*math.pi*2 #平面展開,將周長視為X軸 W = L #Y軸約等于周長一般 H = L/2 #米勒投影中的一個常數(shù),范圍大約在正負(fù)2.3之間 mill = 2.3 #將經(jīng)度從度數(shù)轉(zhuǎn)換為弧度 x = lon*math.pi/180 # 將緯度從度數(shù)轉(zhuǎn)換為弧度 y = lat*math.pi/180 #這里是米勒投影的轉(zhuǎn)換 y = 1.25*math.log(math.tan(0.25*math.pi+0.4*y)) # 這里將弧度轉(zhuǎn)為實際距離 ,轉(zhuǎn)換結(jié)果的單位是公里 x = (W/2)+(W/(2*math.pi))*x y = (H/2)-(H/(2*mill))*y xy_coordinate.append((int(round(x)),int(round(y)))) return xy_coordinate
- xy坐標(biāo)轉(zhuǎn)換成經(jīng)緯度
def xy_to_coor(x, y): lonlat_coordinate = [] L = 6381372 * math.pi*2 W = L H = L/2 mill = 2.3 lat = ((H/2-y)*2*mill)/(1.25*H) lat = ((math.atan(math.exp(lat))-0.25*math.pi)*180)/(0.4*math.pi) lon = (x-W/2)*360/W # TODO 最終需要確認(rèn)經(jīng)緯度保留小數(shù)點后幾位 lonlat_coordinate.append((round(lon,15),round(lat,15))) return lonlat_coordinate
- 調(diào)用兩個函數(shù)
import math print(millerToXY(114.5456282282352,36.631262731204049)) print(xy_to_coor(32805300, 7107399))
打印結(jié)果
[(32805300, 7107399)]
[(114.54562994611808, 36.63126574363226)]
總結(jié)
以上為個人經(jīng)驗,希望能給大家一個參考,也希望大家多多支持腳本之家。
相關(guān)文章
Python如何通過Flask-Mail發(fā)送電子郵件
這篇文章主要介紹了Python如何通過Flask-Mail發(fā)送電子郵件,文中通過示例代碼介紹的非常詳細(xì),對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值,需要的朋友可以參考下2020-01-01python中的循環(huán)結(jié)構(gòu)問題
這篇文章主要介紹了python中的循環(huán)結(jié)構(gòu)問題,具有很好的參考價值,希望對大家有所幫助。如有錯誤或未考慮完全的地方,望不吝賜教2023-03-03Python實現(xiàn)求取表格文件某個區(qū)域內(nèi)單元格的最大值
這篇文章主要介紹基于Python語言,基于Excel表格文件內(nèi)某一列的數(shù)據(jù),計算這一列數(shù)據(jù)在每一個指定數(shù)量的行的范圍內(nèi)(例如每一個4行的范圍內(nèi))的區(qū)間最大值的方法,需要的朋友可以參考下2023-08-08用python監(jiān)控服務(wù)器的cpu,磁盤空間,內(nèi)存,超過郵件報警
這篇文章主要介紹了如果用python監(jiān)控服務(wù)器的cpu,磁盤空間,內(nèi)存,超過郵件報警,幫助大家更好的理解和使用python,感興趣的朋友可以了解下2021-01-01教你怎么用Python操作MySql數(shù)據(jù)庫
很多小伙伴都在問我能不能出一篇怎么用Python操作MySql數(shù)據(jù)庫的教程,今天特地整理了本篇文章,文中有非常詳細(xì)的圖文示例,需要的朋友可以參考下2021-05-05