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

在Python的实际运用中如果你接触到基础测绘学中的Python示例程序时,你是否对其有不解之处?如果你对Python示例程序想更深的了解,你就可以点击以下的文章。希望你会有所收获。

目前成都创新互联已为近1000家的企业提供了网站建设、域名、网页空间、绵阳服务器托管、企业网站设计、扶绥网站维护等服务,公司将坚持客户导向、应用为本的策略,正道将秉承"和谐、参与、激情"的文化,与客户和合作伙伴齐心协力一起成长,共同发展。

基础测绘学之 若干Python示例程序

已经到IEU学了半年测绘,本着短期培训的态度,以最短的时间学了n多的课程,一知半解者甚众。

 好在本科于zju养成了“好读书,不求甚解”的习惯,但凡给公式者绝不深究推导,再编成程序,以后使用时直接调用,岂不美哉。惜到IEU培训居然不让用电脑,身边唯一带点智能的电子玩意就只有一台搭载Symbian S60v3的Nokia N81...

时恰同班好友Mr. BigHead研究跨平台语言python,用于在其WindowsMobile的多普达上编程,每日在其侧滑全键盘上编得不亦乐乎。独乐乐不如众乐乐,遂于我N81上装py之塞班版,与之同习。数月内,编“图幅编号计算”(地图学)、“大地测量”、“复合导线近似平差计算”(现代测量学)等若干小程,整理于下。

(1)地图学 系列比例尺地形图分幅编号计算

[作用] 输入经纬度,计算该点所在各比例尺地图的图幅编号

 
 
 
  1. [code]

图幅编号,输入经纬度,计算所在各比例尺图幅编号

 
 
 
  1. def dfm2d(d,f,m):
  2. temp=m/60.0
  3. temptemp=temp+f
  4. res=d+temp/60.0
  5. return res 
  6. def u(s):
  7. return s.decode('utf8')
  8. def gridnn(n,x,y,xs,xe,ys,ye):
  9. dx=(xe-xs)/(n+0.0)
  10. dy=(ye-ys)/(n+0.0)
  11. for i in range(1,n+1):
  12. if y>=(ye-i*dy):
  13. break
  14. tempy=i
  15. for i in range(1,n+1):
  16. if x<=(xs+i*dx):
  17. break
  18. tempx=i
  19. return (tempy-1)*n+tempx
  20. #main
  21. print u('欢迎使用\n系列比例尺地形图分幅编号计算
    \nPowered by FanZhenzhen IEU Version 1.0\n')
  22. Ld=int(raw_input(u('输入经度L度数:'))) 
  23. Lf=int(raw_input(u('输入经度L分数:')))
  24. Lm=int(raw_input(u('输入经度L秒数:')))
  25. Bd=int(raw_input(u('输入纬度B度数:')))
  26. Bf=int(raw_input(u('输入纬度B分数:')))
  27. Bm=int(raw_input(u('输入纬度B秒数:')))
  28. L=dfm2d(Ld,Lf,Lm)
  29. B=dfm2d(Bd,Bf,Bm)  

计算1:100万图幅编号

 
 
 
  1. Seq100_hang=int(B/4)+1
  2. Seq100_lie=int(L/6)+31
  3. Le100=(Seq100_lie-30)*6
  4. Bn100=Seq100_hang*4
  5. print u('1:100万图幅编号 '),
    Seq100_hang,'-',Seq100_lie

以上文章就是对Python示例程序的部分介绍。

分享标题:关于Python示例程序的基础测绘学的介绍
本文URL:http://www.hantingmc.com/qtweb/news21/273771.html

网站建设、网络推广公司-创新互联,是专注品牌与效果的网站制作,网络营销seo公司;服务项目有等

广告

声明:本网站发布的内容(图片、视频和文字)以用户投稿、用户转载内容为主,如果涉及侵权请尽快告知,我们将会在第一时间删除。文章观点不代表本网站立场,如需处理请联系客服。电话:028-86922220;邮箱:631063699@qq.com。内容未经允许不得转载,或转载时需注明来源: 创新互联