文章分类 - Python
将Geodatabase导出为ShapeFile
摘要:可以将Geodatabase中的Feature Class导出成Shapefile有不同的方式,右键一个Geodatabase可以将其中所有的Feature Class导出到一个文件夹中,右键一个Feature Dataset,可以将其中的所有Feature Class导出到一个文件夹中,右键一个Feature Class可以以Single或Multiple的方式导出到shapefile。 但是,如果想在导出的SHPS中保持Geodatabase中的目录结构,怎么办呢?以下的脚本就是解决的这样问题。—————————————————————————我是一条分割线————————————————
阅读全文
何使用PYTHON 向arcmap中加载shapefile 不是即时python窗口
摘要:如何使用PYTHON 向arcmap中加载shapefile 不是即时python窗口 要求在外部独立运行,加载后保存。即时窗口的加载代码为arcpy.MakeFeatureLayer_management("C:/data/contour.shp","test")通过arcpy,arcobjects,arcengine实现均给分你的这个mxd = arcpy.mapping.MapDocument("CURRENT")是立即窗口 不能在外部运行# 获得当前的documentmxd = arcpy.mapping.MapDocumen
阅读全文
Python与开源GIS:使用OGR中拷贝方法创建新的Shapefile
摘要:摘要: 矢量数据集创建也并没有非常难,就是把读取的时候的顺序重复以便就好了。总体来说这个过程就是构建数据源->构建层->构建要素->构建形状->关闭数据源。我们首先来过看一下如何使用已有的数据来构建我们的数据,也就是复制/粘贴方法。在OGR的矢量数据模型中,矢量数据。我们首先来过看一下如何使用已有的数据来构建我们的数据,也就是复制/粘贴方法。在OGR的矢量数据模型中,矢量数据分为datasource,layer,与feature三个层次,在这三个层次上,OGR有不同的复制数据的方法。在datasource层次创建数据我们使用世界各国的边界数据作为示例,利用这个数据,来创
阅读全文
arcgis中用python代码实现条件赋值
摘要:在字段计算器的代码块对话框写入下列代码:rec=""def autoIncrement(a):global recif (a > 68 and a <=75.65): rec ="4级"elif(a > 75.65 and a <=80.00): rec ="3级"elif(a > 80.00 and a <=84.35): rec ="2级"elif(a > 84.35 and a <=91.50): rec ="1级"else: rec = &
阅读全文
将shapefile文件的属性表用python语言输出成excel形式
摘要:def SaveShpAsCSV(ShpFile,OutDir,OutputName): fields = arcpy.gp.listFields(ShpFile) fieldList2 = [] for field in fields: fieldList2.append(str(field.name)) try: arcpy.ExportXYv_stats(ShpFile,fieldList2,"COMMA",OutDir+"/"+OutputName+".csv","ADD_FIELD_NAMES") exc
阅读全文
使用字段计算器对同一地类图斑自动编号(标记重复记录)
摘要:问题描述:在某个表中把某个字段(如字段一)中具有相同值的记录标出来,并且按照从小到大的排序自动增加一个编号,存储在字段二中,实现如下的效果:FID 字段1 字段21 001 0011 2 001 00123 002 00214 002 0022解决方法:在ArcGIS10中实现方法如下,字段2右键调出字段计算器,选择Python脚本,在代码块中写:UniqueDict = {}def isDuplicateIndex(inValue): UniqueDict.setdefault(inValue,0) UniqueDict[inValue] += 1 return str(inValue)+s
阅读全文
Python提取DEM数据的山脊线代码
摘要:#自动提取山脊线import arcgisscriptingimport osimport sys#创建GP工具gp=arcgisscripting.create(9.3)#DEM参数设置 inputdemlyr=sys.argv[1]#原始DEM最大值demmaxvalue=sys.argv[2]#生成山脊线所在文件夹folderpath=sys.argv[3]gp.CheckOutExtension("Spatial")#生成坡向图层AOutRaster = folderpath+"\sub_demA"try:# Check out ArcGIS S
阅读全文
使用Python开发ArcGIS应用1-获取数据相关信息
摘要:使用Python开发ArcGIS应用1-获取数据相关信息八 15 2010 Published by Jack under Python, arcgis这几天抽空写了一些arcgis python入门类的脚本,详细信息注解以后加入.首先是获取数据信息类的应用. 获取信息函数直接copy的.看起来有点繁琐.先不改了.嗯嗯=.=import arcgisscripting#### get field's info functiondef PrintFiledsInfo(fc): des = gp.Describe(fc) fields = des.Fields fields.r...
阅读全文
ArcGIS遍历区域内要素的属性(python)
摘要:import sysreload(sys)sys.setdefaultencoding(""utf-8"") #应用utf码,默认是ascii码import arcpypath="D:/GISData"arcpy.env.workspace=patharcpy.MakeFeatureLayer_management(path+"/YF/样方分布图.shp",""YF"")#将样方导入图层,定名YFarcpy.MakeFeatureLayer_management(path+&
阅读全文
批量导出mxd至jpg
摘要:批量导出mxd至jpg 原来的时候分别写了两个脚本来用不同的方式来导出图片,后来觉得合成一个为好,于是就对原来的脚本做了一下修改。新功能是可以设定是否遍历子目录以导出其中的mxd成jpg图片(可以设定导出的jpg放到各自mxd相应的文件夹或最顶级文件夹),还可以设定导出图后关机,这样我们就可以下班之后来让电脑来处理这样乏味的工作了。 使用方法是将脚本放到文本文件中,更改后缀为py,用IDLE打开,File - >Open -> run -> run Module。---------------------------------------------------------
阅读全文