關(guān)于Python示例程序的基礎(chǔ)測(cè)繪學(xué)的介紹
在Python的實(shí)際運(yùn)用中如果你接觸到基礎(chǔ)測(cè)繪學(xué)中的Python示例程序時(shí)，你是否對(duì)其有不解之處？如果你對(duì)Python示例程序想更深的了解，你就可以點(diǎn)擊以下的文章。希望你會(huì)有所收獲。
基礎(chǔ)測(cè)繪學(xué)之 若干Python示例程序
已經(jīng)到IEU學(xué)了半年測(cè)繪，本著短期培訓(xùn)的態(tài)度，以最短的時(shí)間學(xué)了n多的課程，一知半解者甚眾。
好在本科于zju養(yǎng)成了“好讀書，不求甚解”的習(xí)慣，但凡給公式者絕不深究推導(dǎo)，再編成程序，以后使用時(shí)直接調(diào)用，豈不美哉。惜到IEU培訓(xùn)居然不讓用電腦，身邊唯一帶點(diǎn)智能的電子玩意就只有一臺(tái)搭載Symbian S60v3的Nokia N81...
時(shí)恰同班好友Mr. BigHead研究跨平臺(tái)語(yǔ)言python，用于在其WindowsMobile的多普達(dá)上編程，每日在其側(cè)滑全鍵盤上編得不亦樂乎。獨(dú)樂樂不如眾樂樂，遂于我N81上裝py之塞班版，與之同習(xí)。數(shù)月內(nèi)，編“圖幅編號(hào)計(jì)算”(地圖學(xué))、“大地測(cè)量”、“復(fù)合導(dǎo)線近似平差計(jì)算”(現(xiàn)代測(cè)量學(xué))等若干小程，整理于下。
(1)地圖學(xué) 系列比例尺地形圖分幅編號(hào)計(jì)算
[作用] 輸入經(jīng)緯度，計(jì)算該點(diǎn)所在各比例尺地圖的圖幅編號(hào)
- [code]
圖幅編號(hào)，輸入經(jīng)緯度，計(jì)算所在各比例尺圖幅編號(hào)
- 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系列比例尺地形圖分幅編號(hào)計(jì)算
\nPowered by FanZhenzhen IEU Version 1.0\n')- Ld=int(raw_input(u('輸入經(jīng)度L度數(shù):')))
- Lf=int(raw_input(u('輸入經(jīng)度L分?jǐn)?shù):')))
- Lm=int(raw_input(u('輸入經(jīng)度L秒數(shù):')))
- Bd=int(raw_input(u('輸入緯度B度數(shù):')))
- Bf=int(raw_input(u('輸入緯度B分?jǐn)?shù):')))
- Bm=int(raw_input(u('輸入緯度B秒數(shù):')))
- L=dfm2d(Ld,Lf,Lm)
- B=dfm2d(Bd,Bf,Bm)
計(jì)算1:100萬(wàn)圖幅編號(hào)
- 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萬(wàn)圖幅編號(hào) '),
Seq100_hang,'-',Seq100_lie
以上文章就是對(duì)Python示例程序的部分介紹。
【編輯推薦】
