關于Python示例程序的基礎測繪學的介紹
在Python的實際運用中如果你接觸到基礎測繪學中的Python示例程序時,你是否對其有不解之處?如果你對Python示例程序想更深的了解,你就可以點擊以下的文章。希望你會有所收獲。
基礎測繪學之 若干Python示例程序
已經到IEU學了半年測繪,本著短期培訓的態度,以最短的時間學了n多的課程,一知半解者甚眾。
好在本科于zju養成了“好讀書,不求甚解”的習慣,但凡給公式者絕不深究推導,再編成程序,以后使用時直接調用,豈不美哉。惜到IEU培訓居然不讓用電腦,身邊唯一帶點智能的電子玩意就只有一臺搭載Symbian S60v3的Nokia N81...
時恰同班好友Mr. BigHead研究跨平臺語言python,用于在其WindowsMobile的多普達上編程,每日在其側滑全鍵盤上編得不亦樂乎。獨樂樂不如眾樂樂,遂于我N81上裝py之塞班版,與之同習。數月內,編“圖幅編號計算”(地圖學)、“大地測量”、“復合導線近似平差計算”(現代測量學)等若干小程,整理于下。
(1)地圖學 系列比例尺地形圖分幅編號計算
[作用] 輸入經緯度,計算該點所在各比例尺地圖的圖幅編號
- [code]
圖幅編號,輸入經緯度,計算所在各比例尺圖幅編號
- def dfm2d(d,f,m):
- temp=m/60.0
- temptemp=temp+f
- res=d+temp/60.0
- return res
- def u(s):
- return s.decode('utf8')
- def gridnn(n,x,y,xs,xe,ys,ye):
- dx=(xe-xs)/(n+0.0)
- dy=(ye-ys)/(n+0.0)
- for i in range(1,n+1):
- if y>=(ye-i*dy):
- break
- tempy=i
- for i in range(1,n+1):
- if x<=(xs+i*dx):
- break
- tempx=i
- return (tempy-1)*n+tempx
- #main
- print u('歡迎使用\n系列比例尺地形圖分幅編號計算
\nPowered by FanZhenzhen IEU Version 1.0\n')- Ld=int(raw_input(u('輸入經度L度數:')))
- Lf=int(raw_input(u('輸入經度L分數:')))
- Lm=int(raw_input(u('輸入經度L秒數:')))
- Bd=int(raw_input(u('輸入緯度B度數:')))
- Bf=int(raw_input(u('輸入緯度B分數:')))
- Bm=int(raw_input(u('輸入緯度B秒數:')))
- L=dfm2d(Ld,Lf,Lm)
- B=dfm2d(Bd,Bf,Bm)
計算1:100萬圖幅編號
- Seq100_hang=int(B/4)+1
- Seq100_lie=int(L/6)+31
- Le100=(Seq100_lie-30)*6
- Bn100=Seq100_hang*4
- print u('1:100萬圖幅編號 '),
Seq100_hang,'-',Seq100_lie
以上文章就是對Python示例程序的部分介紹。
【編輯推薦】