添加xy坐標文件生成圖形-度分秒轉換


1.在arcgis中創建的模型如下:
2.說明
(1)分度秒需要轉換成十進制度,需要編寫對經度、緯度的處理。
(2)txt中的坐標無法直接處理度分秒,所以臨時先轉換成點要素,對點要素中的經度、緯度坐標進行處理,然后導出坐標表(dbf),再進行一次展點。
3.說明
對於一些坐標如112.2400,113.0000臨時轉換成點要素時,arcgis自動把這些坐標(文本)當成了浮點型,去掉了所有0,需要編寫額外的腳本處理這種情況。
4.python對經緯度處理的腳本
(1)經度
def calJD(x):
s = str(x) +'0000'
s = s.format("000.0000")
du=s[0:3]
fen=s[4:6]
mia=s[6:8]
degree=float(du)
minix=float(fen)
secdx=float(mia)
res=degree+minix/60+secdx/3600
return '%.6f'%res
(2)緯度
def calWD(y):
s=str(y)+'0000'
s = s.format('00.0000')
du=s[0:2]
fen=s[3:5]
mia=s[5:7]
degree=float(du)
mini=float(fen)
secd=float(mia)
res=degree+mini/60+secd/3600
return '%.6f' % res
 
坐標文件示例如下:
112.4900,37.4345
112.5700,37.4345
112.5700,37.4000
113.0000,37.4000
113.0000,37.3430
112.4900,37.3430


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM