关于Python示例程序的基础测绘学的介绍

关于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示例程序的部分介绍。