加入字段:
>>> arc = ['A','B','C','D','E']
>>> for i in range(5):
... arcpy.AddField_management("idcounty",arc[i],"TEXT")
给idcounty空间数据批量加入五个字段~
Buffer缓冲区
arcpy.Buffer_analysis("thermal","buffer","10 kilometers")
集合面积与几何长度
>>> arcpy.CalculateField_management("idcounty","ID_Area","!shape.area@squarekilometers!","PYTHON_9.3")
>>> arcpy.CalculateField_management("idcounty","ID_Area","!shape.length@kilometers!","PYTHON_9.3")
改变工作空间
>>> arcpy.env.workspace = "F:/Data"
>>> result = arcpy.Buffer_analysis("thermal","t_Buffer","10 kilometers")
>>> print result
F:/Data\t_Buffer.shp
返回要素数目
>>> result = arcpy.GetCount_management("idcounty")
>>> print result.getOutput(0)
44
调用工具的方法,就是工具的英文名称,去掉中间的空格,然后下划线,加入工具集的名称
字段名称:
>>> fieldList = arcpy.ListFields("idcounty")
>>> for field in fieldList:
... print field.aliasname + field.type
获取工作空间内所有要素类的字段名:
>>> fes = arcpy.ListFeatureClasses()
>>> for fe in fes:
... print fe
... nl = ''
... fs = arcpy.ListFields(fe)
... for f in fs:
... nl = nl + ' ' + f.aliasName
... print nl
复制(工作空间中的可以直接写名字)
arcpy.Copy_management("i_Copy.shp","F:/Data/data/qq.shp")
列表的一些函数
>>> a = [1]*10
>>> for i in range(10):
... a[i] = i
...
>>> a
[0, 1, 2, 3, 4, 5, 6, 7, 8, 9]
>>> b = [1]*2
>>> b
[1, 1]
>>> a.extend(b)
>>>
>>> a
[0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 1]
>>> a + b
[0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 1, 1, 1]
>>> a
[0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 1]
>>> a.count(1)
3
>>> a.append(100)
>>> a
[0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 1, 100]
>>> a.index(100)
12
>>> a.insert(0,'I')
>>> a
['I', 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 1, 100]
>>> a.pop(0)
'I'
>>> a
[0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 1, 100]
>>> a.remove(100)
>>> a
[0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 1]
>>> a.reverse()
>>> a
[1, 1, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0]
>>> a.sort()
>>> a
[0, 1, 1, 1, 2, 3, 4, 5, 6, 7, 8, 9]
>>>
Python数据类型
>>> type(3)
>>> type(3.0)
>>> type(1111111111)
>>> type(111111111111111)
>>> type(3.00000000)
>>> type(3.000000000000)
>>> type(1+2j)
>>> type(True)
>>> type('Alex')
>>> type([2,4])
>>> type((3,4))
>>>
详见: http://blog.sina.com.cn/s/blog_4b5039210100e9ya.html
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
------------------------------------------------ ESRI培训 ----------------------------------------------------------------------------------------------------------------------
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
地图文件mxd
>>> mxd = arcpy.mapping.MapDocument("current")
>>> print mxd.filePath
F:\MY_OWN_WORK\Exercise\中国.mxd
数据框架data frame
>>> mxd = arcpy.mapping.MapDocument("current")
>>> dfs = arcpy.mapping.ListDataFrames(mxd)
>>> for df in dfs:
... print df.name
...
图层
Data Frame II
图层layer
>>> mxd = arcpy.mapping.MapDocument("current")
>>> df = arcpy.mapping.ListDataFrames(mxd)[0]
>>> ls = arcpy.mapping.ListLayers(df)
>>> for l in ls:
... print l
...
Cities (population > 5 Million)
Geogrid
Rivers
Lakes
Continents
Ocean
dataframe.extent = layers[0].getSelectedExtent()