Python脚本入门学习经典手册.doc
《Python脚本入门学习经典手册.doc》由会员分享,可在线阅读,更多相关《Python脚本入门学习经典手册.doc(67页珍藏版)》请在咨信网上搜索。
1、精品文档Python脚本使用详解目录写在前面的话2前言2一、PYTHON语言基础31数学运算符32字符串操作43模块的使用(Modules)54使用def构建函数65流程控制结构:If,While,For66简单输入和输出9二、ARCGIS&PYTHON101如何创建地理处理对象(geoprocessor object)102获取地理处理帮助102.1举例:如何使用Geoprocessor Programming Model中的Lists113使用地理处理工具Toolboxes和Aliases114在建模中使用脚本(Scripts in ModelBuilder)125 在PythonWin里
2、调试地理处理脚本185.1 调试选择和消息195.2PythonWin的调试工具205.3地理处理工具举例216使用描述(Describe)和存在(Exists)获取数据信息216.1描述226.2存在(Exists)236.3在循环中使用描述和存在237在Python脚本中使用地图代数(Map Algebra)268数据管理和指针(Data Management and Cursors)278.1数据管理(Data Management)278.2指针(Cursors)28附录1:地理处理脚本中输入&输出方法指南31附录2:其他32写在前面的话一直想学习ArcGIS中的Python脚本,大四
3、下半学期终于有了时间,可是想找到这么一本好的教材不容易。茫茫互联网,终于找到了旧金山州立大学Jerry Davis教授的个人主页,对其中Geoprocessing Scripts With Python如获至宝,独乐乐不如众乐乐,现在将其教程翻译并结合自己的学习情况给出总结。希望能够给更多想学习Python的同学一个参考。另外,在我刚开始接触Python时,是看了台湾辅仁大学一位老师的视频课件,在此致谢。我想从两个大部分总结:一、Python语言基础;二、ArcGIS&Python。其中第一部分参考了Python精要参考(第二版)、Python编程金典(读书笔记)等书籍文献。对于多数读者来说,
4、可能或多或少有一些编程基础,所以理解起来应该不成问题。文中多数数据来自Jerry Davis教授的主页,放在“C:prog”目录下,为了直观,我将运算结果一并编辑,方便参考。值得一提的是ArcGIS的在线帮助文档,一个实时更新的GIS宝库,很多专业性知识都可以找到答案,点击链接ArcGIS10中文帮助、ArcGIS9.3.1或9.3英文帮助。获取更过脚本例子来学习 :ESRI的地理处理模型和脚本工具库是个不错的选择。由于我也是初次接触,翻译或者心得难免有纰漏之处,希望同仁们可以多多交流!前言在GIS建模或GIS数据管理中,你可能经常需要处理一系列步骤才可以完成的工作;你可能有一个工作目录下的数
5、据需要重投影、裁剪到研究区域,或者用某种方法组合成期望的结果;我们也经常需要根据不同情形用不同方法处理数据,因此我们需要作出选择,而高质量的决策需要考虑很多低水平的决策,这可以通过脚本程序模型辅助完成。脚本编程的主要目的是使枯燥的处理数据工作自动化,通过逻辑来指挥处理过程。我想自动化和逻辑是关键,它们区别于我们多数使用计算机时的交互活动。我们发E-mail,写文章或者设计地图,都需要和计算机交互,而处理一系列数据,我们需要自动化和利用逻辑来指导自动化。在地理处理脚本逻辑中,我们需要在允许我们做的事情中作出决定,比如,处理栅格数据不同于矢量数据,或为没投影的数据设置投影,或处理仅在特定时间搜集的
6、数据集。对于重要的GIS工作来说,脚本以及其他形式的程序是必需的,而非可有可无。在接下来的联系中,我们会探索Python的使用以及创建脚本来使用ArcGIS里众多的地理处理工具。所有你能在ArcToolbox或Model中使用的工具都能够用在Python脚本中,这些脚本可以生成脚本工具,像其他地理处理工具一样使用。一、Python语言基础安装PythonWin,在ArcGisDesktop9.3.isoDesktopPythonWin目录下可以找到PythonWin的安装程序,默认是不安装的,。同时会安装win32com以及允许任何脚本在基于Dispatch的地理处理过程中工作。ArcGIS1
7、0中引入了全新的Python Window来增强内嵌的Python体验。警告:不要尝试更新随ArcGIS安装的Python到一个新的版本!下面介绍Python的一些简单语法和规则。1数学运算符Python提供了多样化的通用数学运算符多数编程语言的特征,以及许多通过import的modules提供的符号。常用的有+,-,*,/,*(幂),%(取模,即除后的余数)。下面的表格显示了整型(Integer)和浮点型(Float)各种组合运算的结果,记住一条规则,只要参与运算的有浮点型,则结果为浮点型;全为整型时,结果才为整型。输入表达式结果Notes2+35整型结果2.+35.02.是浮点型,结果浮点
8、型2-3-12*36整型结果2.*36.0浮点型5/22整型5./22.55%21取模Az=270Newaz=az+180Print newaz%36090取模的用途之一方位角加180后逆转方向5*22525*0.55.0没有sqrt()功能,除非添加math模块2字符串操作注:使用Python帮助:有超过30种内置方法来处理字符,请到Sequence Types下的String Methods寻找帮助!字符串是一串字母,比如San Francisco,字符串下标从0开始。学习字符串语法的最好方法是自己动手尝试,下标展示之:输入结果Notesprint zhulj.capitalize()Zh
9、uljs.capitalize()即将capitalize()方法用于ss=zhuljprint s.capitalize()print s0zStrings可以像一个字母列表一样处理,第一个字母下标为0,某个字符段可以用1:3来格式化:从第1个的开头到第3个的开头,不包括下标为3的字母;s-1表示倒过来第一个,相当于slen(s)-1s1=s1print s1hprint s-2:ljprint s2:3uprint s2:4ulprint s2:,s:5ulj zhuljs2=s.upper()print s2ZHULJ我们可以将字符串方法的结果赋给新的变量s3=s+s2print s3z
10、huljZHULJ字符串组合用“+”print s*3zhuljzhuljzhulj字符串重复用“*”,后为重复次数selstr=elev1000print selstrelev1000字符串可以使用单引号或双引号,跨行时用双引号。othersel=”elev1000”print otherselelev1000print s.isupper()False一些方法返回值为布尔型(True或False),一些返回索引值(下标值)print s2.isupper()Truep=d:/work/lu.shpprint p.find(.)10print p.find(/)2plist=p.split(
11、/)print plistd:, work, lu.shp你可以用split()方法解析出不同的字符串片段,并创建一个列表(List),我们可以使用其中不同的元素print plist0d:print plist1workp2=d:worksoil.shpprint p2d:worksoil.shp反斜线“”和某些字母一起有特殊用法,如n为换行,“”为转义字符,如“”则表示“”print Jerrys KidsJerrys Kidsprint JerrysnKidsJerrysKidsp3=rd:worksoil.shpprint p3d:worksoil.shp字符串前加“r”则强制“”代表
12、其本身,而非转义字符,这对于文件路径的操作很方便3模块的使用(Modules)Python提供了一系列内置的方法(大量依赖于模块)用于通用编程。Python安装时自带了大量Modules,最常用的有math,sys,random,array以及os.path。当然还有好多Modules可以下载,比如数字处理(Numeric)numpy,可在www.python.org或里搜索。www.python.org/moin/NumericAndScientific页面中列举了一些。使用Module前,必须import之。通常我们会将一行import 放在程序顶部,比如:import arcgisscr
13、ipting当然,这不必成为你程序的第一行,但必须在使用它里面方法之前。当要引用多个模块是,中间用逗号分隔,比如:import arcgisscripting,sys,string,os,math我们也可以自己为频繁使用的方法创建Module,下面,我们开始体验内置的Modules。math和random模块很多常用的数学计算功能都可以通过math找到,比如三角计算或对数计算,如果要使用复杂数字,就使用cmath模块。和之前一样,通过以下表格来体现模块的使用:输入结果Notesimport mathprint math.log10(100)2.0以10为底的对数print math.log(1
14、00)4.60517018599自然对数print math.pi3.14159265359是一个静态常量,所以不需要括弧pi=math.piprint pi3.14159265359如果不想总是输入“math.pi”可以将其赋给一变量pi3.1415926535897931不需要print即可查看变量值print math.sin(radians)print math.cos(radians)print math.tan(radians)三角函数的计算是弧度制degrad=pi/18045*degrad0.78539816339744828度转化为弧度sin=math.sinsin(45*d
15、egrad)sin(90*degrad)0.707106781186547461.0即使功能函数(像sin)都可以赋给一个变量math.e2.7182818284590451math.hypot(3,4)5.0此方法是求三角形的斜边x1=520382;y1=4152373x2=520475;y2=4152963不同赋值语句间用“;”分隔xr=x2-x1yr=y2-y1math.hypot(xr,yr)math.sqrt(xr*2+yr*2)(xr*xr+yr*yr)*0.5597.28468923956189不同的方式,相同的结果import randomrandom.random()0.27
16、281588185756478random()方法,每次结果都不同,值域为0.0,1.0)rnd=random.randomrnd()0.4456393835072503mu=50s=10print random.gauss(mu,s)46.5282021944使用def构建函数有点像Module,但更简单,函数是一个自己定义功能,用在之后的代码中,并且提供任何你想要使用的参数。这个函数从此可像变量那样在程序中使用,结合例子更容易理解。接下来的代码定义了一个将度转换为弧度的简单函数,同时也定义了一个弧度转换为度的函数,它们和Excel内置的函数类似。import mathdef radians
17、(angdeg): return angdeg*math.pi/180def degrees(angrad): return angrad*180/math.piprint math.sin(radians(45)print degrees(math.acos(0.5)运行之,得到结果:0.707106781187 60.0 5流程控制结构:If,While,For任何脚本或编程语言的一个重要特征就是执行一系列不同情形语句的能力。你想要创建一系列山影栅格来代表夏天、冬天和春秋分。山影(hillshade)工具需要有太阳高度角和方位角作为输入参数。重要日期太阳倾角夏至(6月21日)23.44春秋
18、分(3月21日,9月21日)0冬至(12月21日)-23.44接下来是一段相当简单的代码,通过太阳倾角(太阳光线正午垂直照射的纬度)获取太阳角和方位角以及纬度。输入两个参数:lat(研究区域的纬度,南半球为负)和decl(太阳倾角),由此得到sunangle和azimuth:lat=30decl=20sunangle=90-lat+declazimuth=180if sunangle90: sunangle=180-sunangle azimuth=0print sunangle,azimuth上面的例子中lat和decl强制赋了值。有三种流程控制操作:if 仅在一个特定情形下才执行语句;wh
19、ile 当一种情形存在下,持续执行语句for 遍历一系列值这些语法和def有些相似:初始语句后加顿号、需要执行的语句块有缩进。这三个结构的一些重要的公共特征:if、while、for语句均以冒号结尾,接下来是缩进的代码块,用于if、while、for定义的情形。在脚本编写窗口,你会发现,你在一行末尾打上冒号后,下一行自动缩进,在接下来的一行按下退格键取消缩进。如果你只需做一件事情,你可以在冒号后面同一行添加简短的语句,比如:if x0: print x 比0大print 下一行不要缩进了。if(continued)接下来,我们会探索一下另一个方便的模块:os.path:开始之前,在d:/下创建
20、一个“testfolder”文件夹,然后新建一个“test.txt”文件;尝试以下代码段,确保print语句前有缩进。import os.pathif os.path.exists(d:/testfolder/test.txt): print 测试文件夹存在 print txt文件存在elif os.path.exists(d:/testfolder): print 测试文件夹存在 print 测试文件夹存在,但txt文件不存在else: print 两者都不存在可选探索示例接下来的例子做的事情对GIS非常重要,但是实际上不用任何地理处理代码。USGS7.5米分辨率DEM(数字高程模型)是文本
21、文件(USGSDEM文件),投影为UTM,UTM北向和东向单位是米,但是高程单位可能是英尺(feet)或米(meters)。因此在获取垂直或水平距离信息时会有问题,比如坡度可以通过垂直距离/水平距离获得。如果你不在使用Z值之前设置为0.3048,将会出现错误结果。但是不幸的是,你可能不知道DEM文本文件的垂直单位是英尺还是米。这些信息保存在第539个字符里,“1”代表英尺,“2”代表米,所以可以通过读取这个文件判断。下面的脚本演示了上述内容:import fileinputinfile=rc:progpendatawoodside.demfirstline=fileinput.input(in
22、file)0unitchar=firstline539unit=(unknown:not a 7.5 DEM?)if unitchar=1:unit=feetif unitchar=2:unit=metersprint nElevation in+ +unitfileinput.close()输出结果:Elevation in feetwhile(continued) 运行下面的代码,说明了一种while循环:x=1while x10: print x x=x+1屏幕依次输出19 下面说明一下“=”(等于)的概念:x=5z=x=4print z输出Falsex=5z=x=5print z输出T
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Python 脚本 入门 学习 经典 手册
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【w****g】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【w****g】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。