ArcPy学习

第1章面向ArcGIS的Python编程语言的基础

第2章管理地图文档和图层

引用当前的地图文档
引用磁盘上的地图文档
获取地图文档的图层列表
限制图层列表
缩放至所选要素
改变地图范围
添加图层到地图文档
插入图层到地图文档
更新图层的符号系统
更新图层属性
操作数据框中的启用时间的图层

引用当前的地图文档2.2-28

import arcpy.mapping as mapping

mxd = mapping.MapDocument("CURRENT") # 引用当前活动的地图文档

print(mxd.title)

mxd.title = "Copy of Crime Project"

mxd.saveACopy("c:/ArcpyBook/Ch2/crime_copy.mxd")

引用磁盘上的地图文档2.3-30

import arcpy.mapping as mapping

mxd = mapping.MapDocument("c:/ArcpyBook/Ch2/crime_copy.mxd")

print(mxd.title)

获取地图文档的图层列表2.4-31

import arcpy.mapping as mapping

mxd = mapping.MapDocument("CURRENT")

layers = mapping.ListLayers(mxd)

for lyr in layers:

    print(lyr.name)

限制图层列表2.5-33

import arcpy.mapping as mapping

mxd = mapping.MapDocument("CURRENT")

for df in mapping.ListDataFrames(mxd):

    if (df.name == 'Crime'):

        layers = mapping.ListLayers(mxd,"Burg*",df)

        for layer in layers:

            print(layer.name)

缩放至所选要素2.6-35

import arcpy.mapping as mapping

mxd = mapping.MapDocument("CURRENT")

df = mapping.ListDataFrames(mxd,"Crime")[0]

layer = mapping.ListLayers(mxd,"Burglaries*",df)[0]

df.extent = layer.getSelectedExtent()

改变地图范围2.7-37

import arcpy.mapping as mapping

mxd = mapping.MapDocument("CURRENT")

for df in mapping.ListDataFrames(mxd):

    if (df.name == 'Crime'):

        layers = mapping.ListLayers(mxd,'Crime Density by School District',df)

        for layer in layers:

            query = '"NAME" = \'Lackland ISD\''

            layer.definitionQuery = query

            df.extent = layer.getExtent()

添加图层到地图文档2.8-39

import arcpy.mapping as mapping

mxd = mapping.MapDocument("CURRENT")

df = mapping.ListDataFrames(mxd)[0]

layer = mapping.Layer(r"C:\ArcpyBook\data\School_Districts.lyr")

mapping.AddLayer(df,layer,"AUTO_ARRANGE")

插入图层到地图文档2.9-42

import arcpy.mapping as mapping

mxd = mapping.MapDocument("CURRENT")

df = mapping.ListDataFrames(mxd, "Crime")[0]

refLayer = mapping.ListLayers(mxd, "Burglaries*", df)[0]

insertLayer = mapping.Layer(r"C:\ArcpyBook\data\CityOfSanAntonio.gdb\Crimes2009")

mapping.InsertLayer(df,refLayer,insertLayer,"BEFORE")

更新图层的符号系统2.10-45

import arcpy.mapping as mapping

mxd = mapping.MapDocument("CURRENT")

df = mapping.ListDataFrames(mxd, "Crime")[0]

updateLayer = mapping.ListLayers(mxd,"Crime Density by School District",df)[0]

sourceLayer = mapping.Layer(r"C:\ArcpyBook\data\CrimeDensityGradSym.lyr")

mapping.UpdateLayer(df,updateLayer,sourceLayer,True)

更新图层属性2.11-48

import arcpy.mapping as mapping

mxd = mapping.MapDocument("CURRENT")

df = mapping.ListDataFrames(mxd, "Crime")[0]

updateLayer = mapping.ListLayers(mxd,"Crimes2009",df)[0]

sourceLayer = mapping.Layer(r"C:\ArcpyBook\data\BurglariesNoForcedEntry.lyr")

mapping.UpdateLayer(df,updateLayer,sourceLayer,False)

操作数据框中的启用时间的图层2.12-53

import arcpy.mapping as mapping, os

mxd = mapping.MapDocument("CURRENT")

df = mapping.ListDataFrames(mxd, "Crime")[0]

dft = df.time

dft.currentTime = dft.startTime

while dft.currentTime <= dft.endTime:

    fileName = str(dft.currentTime).split(" ")[0] + ".pdf"

    mapping.ExportToPDF(mxd,os.path.join(r"C:\ArcpyBook\Ch2", fileName))

    print("Exported " + fileName)

    dft.currentTime = dft.currentTime + dft.timeStepInterval

第3章查找和修复丢失的数据链接

查找地图文档和图层文件中丢失的数据源
使用MapDocument.findAndReplaceWorkspacePaths()方法修复丢失的数据源
使用MapDocument.replaceWorkspaces()方法修复对视的数据源
使用replaceDataSource()方法修复单个图层和表对象
查找文件夹中所有地图文档内丢失的数据源

查找地图文档和图层文件中丢失的数据源3.2-60

import arcpy.mapping as mapping

mxd = mapping.MapDocument(r"c:\ArcpyBook\Ch3\Crime_BrokenDataLinks.mxd")

listBrokenDS = mapping.ListBrokenDataSources(mxd)

for layer in listBrokenDS:

    print(layer.name)

使用MapDocument.findAndReplaceWorkspacePaths()方法修复丢失的数据源3.3-62

import arcpy.mapping as mapping

mxd = mapping.MapDocument(r"c:\ArcpyBook\Ch3\Crime_BrokenDataLinks.mxd")

mxd.findAndReplaceWorkspacePaths(r"C:\ArcpyBook\Ch3\Data\OldData\CityOfSanAntonio.gdb", r"C:\ArcpyBook\Data\CityOfSanAntonio.gdb")

mxd.saveACopy(r"C:\ArcpyBook\Ch3\Crime_DataLinksFixed.mxd")

使用MapDocument.replaceWorkspaces()方法修复对视的数据源3.4-65

import arcpy.mapping as mapping

mxd = mapping.MapDocument(r"c:\ArcpyBook\Ch3\Crime_DataLinksFixed.mxd")

mxd.replaceWorkspaces(r"c:\ArcpyBook\data\CityOfSanAntonio.gdb", "FILEGDB_WORKSPACE",r"c:\ArcpyBook\new_data\CityOfSanAntonio_Personal.mdb","ACCESS_WORKSPACE")

mxd.saveACopy(r"c:\ArcpyBook\Ch3\Crime_DataLinksUpdated.mxd")

使用replaceDataSource()方法修复单个图层和表对象3.5-68

import arcpy.mapping as mapping

mxd = mapping.MapDocument(r"c:\ArcpyBook\Ch3\Crime_DataLinksLayer.mxd")

df = mapping.ListDataFrames(mxd,"Crime")[0]

lyr = mapping.ListLayers(mxd,"Burglary",df)[0]

lyr.replaceDataSource(r"c:\ArcpyBook\data","SHAPEFILE_WORKSPACE","Burglaries_2009")

mxd.saveACopy(r"c:\ArcpyBook\Ch3\Crime_DataLinksNewLayer.mxd")

查找文件夹中所有地图文档内丢失的数据源3.6-72

import arcpy.mapping as mapping, os

f = open('BrokenDataList.txt', 'w')

for root, dirs, files in os.walk("c:\ArcpyBook"):

    for name in files:

        filename = os.path.join(root, name)

        if ".mxd" in filename:

            mxd = mapping.MapDocument(filename)

            f.write("MXD: " + filename + "\n")

            brknList = mapping.ListBrokenDataSources(mxd)

            for brknItem in brknList:

                print("Broken data item: " +  brknItem.name + " in " + filename)

                f.write("\t" + brknItem.name + "\n")

print("All done")

f.close()

第4章自动化地图制图和打印

创建布局元素的python列表
为布局元素指定唯一的名称
使用ListLayoutElements()函数限制返回的布局元素
更新布局元素的属性
获取可用打印机的列表
使用PrintMap()函数打印地图
导出地图为PDF文件
导出地图为图像文件
导出报表
使用数据驱动页面和ArcPy制图模块构建地图册
将地图文档发布为ArcGIS Server服务

创建布局元素的python列表4.2-77

import arcpy.mapping as mapping

mxd = mapping.MapDocument("CURRENT")

for el in mapping.ListLayoutElements(mxd):

    if el.name != "":

        print el.name

为布局元素指定唯一的名称4.3-79
```
None
```

使用ListLayoutElements()函数限制返回的布局元素4.4-83

import arcpy.mapping as mapping

mxd = mapping.MapDocument("CURRENT")

for el in mapping.ListLayoutElements(mxd,"LEGEND_ELEMENT","*Crime*"):

    print el.name

更新布局元素的属性4.5-84

import arcpy.mapping as mapping

mxd = mapping.MapDocument("CURRENT")

elLeg = mapping.ListLayoutElements(mxd, "LEGEND_ELEMENT","*Crime*")[0]

elLeg.title = "Crimes by School District"

for item in elLeg.listLegendItemLayers():

    print item.name

获取可用打印机的列表4.6-87

import arcpy.mapping as mapping

mxd = mapping.MapDocument("CURRENT")

for printerName in mapping.ListPrinterNames():

    print printerName

使用PrintMap()函数打印地图4.7-88

import arcpy.mapping as mapping

mxd = mapping.MapDocument("CURRENT")

for df in mapping.ListDataFrames(mxd):

    if df.name == "Test_Performance":

        mapping.PrintMap(mxd,"",df)

导出地图为PDF文件4.8-90

import arcpy.mapping as mapping

mxd = mapping.MapDocument("CURRENT")

mapping.ExportToPDF(mxd,r"c:\ArcpyBook\Ch4\Map_PageLayout.pdf")

import arcpy.mapping as mapping

mxd = mapping.MapDocument("CURRENT")

for df in mapping.ListDataFrames(mxd):

    if df.name == "Crime":

        df.referenceScale = df.scale

        mapping.ExportToPDF(mxd,r"c:\ArcpyBook\Ch4\DataFrameCrime.pdf",df)

导出地图为图像文件4.9-92

import arcpy.mapping as mapping

mxd = mapping.MapDocument("CURRENT")

for df in mapping.ListDataFrames(mxd):

    if df.name == "Crime":

        mapping.ExportToJPEG(mxd,r"c:\ArcpyBook\Ch4\DataFrameCrime.jpg",df)

导出报表4.10-93

import arcpy

import os 

path = os.getcwd()

#Create PDF and remove if it already exists

pdfPath = path + r"\CrimeReport.pdf"

if os.path.exists(pdfPath):

    os.remove(pdfPath)

pdfDoc = arcpy.mapping.PDFDocumentCreate(pdfPath)

districtList = ["Harlandale", "East Central", "Edgewood", "Alamo Heights", "South San Antonio", "Southside", "Ft Sam Houston","North East", "Northside", "Lackland", "Southwest", "Judson", "San Antonio"]

mxd = arcpy.mapping.MapDocument(path + r"\Crime_Ch4.mxd")

df = arcpy.mapping.ListDataFrames(mxd)[0]

lyr = arcpy.mapping.ListLayers(mxd, "Crime Density by School District")[0]

pageCount = 1

for district in districtList:

    #Generate image for each district

    whereClause = "\"NAME\" = '" + district + " ISD'"

    lyr.definitionQuery = whereClause

    arcpy.SelectLayerByAttribute_management(lyr, "NEW_SELECTION", whereClause)

    df.extent = lyr.getSelectedExtent()

    arcpy.SelectLayerByAttribute_management(lyr, "CLEAR_SELECTION")

    arcpy.mapping.ExportToBMP(mxd, path + "\DistrictPicture.bmp", df) #single file

    #Generate report

    print("Generating report for: " + district + " ISD")

    arcpy.mapping.ExportReport(report_source=lyr,report_layout_file=path + r"\CrimeLayout.rlf",output_file=path + r"\temp" + str(pageCount) + ".pdf", starting_page_number=pageCount)

    #Append pages into final output

    print("Appending page: " + str(pageCount))

    pdfDoc.appendPages(path + r"\temp" + str(pageCount) + ".pdf")

    os.remove(path + r"\temp" + str(pageCount) + ".pdf")

    pageCount = pageCount + 1

pdfDoc.saveAndClose()

del mxd

使用数据驱动页面和ArcPy制图模块构建地图册4.11-98

import arcpy

import os

# Create an output directory variable

outDir = r"C:\ArcpyBook\Ch4"  

# Create a new, empty pdf document in the specified output directory

finalpdf_filename = outDir + r"\MapBook.pdf"

if os.path.exists(finalpdf_filename):

  os.remove(finalpdf_filename)

finalPdf = arcpy.mapping.PDFDocumentCreate(finalpdf_filename)

# Add the title page to the pdf

print("Adding the title page  \n")

finalPdf.appendPages(outDir + r"\TitlePage.pdf")

# Add the index map to the pdf

print "Adding the index page  \n"

finalPdf.appendPages(outDir + r"\MapIndex.pdf")

# Export the Data Driven Pages to a temporary pdf and then add it to the

# final pdf. Alternately, if your Data Driven Pages have already been

# exported, simply append that document to the final pdf.

mxdPath = outDir + r"\Topographic.mxd"

mxd = arcpy.mapping.MapDocument(mxdPath)

print("Creating the data driven pages \n")

ddp = mxd.dataDrivenPages

temp_filename = outDir + r"\tempDDP.pdf"

if os.path.exists(temp_filename):

  os.remove(temp_filename)

ddp.exportToPDF(temp_filename, "ALL")

print("Appending the map series  \n")

finalPdf.appendPages(temp_filename)

# Update the properties of the final pdf

finalPdf.updateDocProperties(pdf_open_view="USE_THUMBS",

                             pdf_layout="SINGLE_PAGE")

# Save your result

finalPdf.saveAndClose()

# remove the temporary data driven pages file

if os.path.exists(temp_filename):

    print("Removing the temporary map series file")

    os.remove(temp_filename)

# Delete variables

#del finalPdf, mxd, ddp

将地图文档发布为ArcGIS Server服务4.12-102

import arcpy.mapping as mapping

wrkspc = r'c:\ArcpyBook\ch4'

mxd = mapping.MapDocument(wrkspc + r"\Crime.mxd")

service = 'Crime'

sddraft = wrkspc + service + '.sddraft'

mapping.CreateMapSDDraft(mxd, sddraft, service)

analysis = mapping.AnalyzeForSD(wrkspc + "Crime.sddraft")

for key in ('messages', 'warnings', 'errors'):

    print("----" + key.upper() + "----")

    vars = analysis[key]

    for ((message, code), layerlist) in vars.iteritems():

        print "    ", message, " (CODE %i)" % code

        print("     applies to:")

        for layer in layerlist:

            print(layer.name)

第5章使用脚本执行地理处理工具

查找地理处理工具
查看工具箱别名
使用脚本执行地理处理工具
讲一个工具的输出作为另一个工具的输入

查找地理处理工具5.2-110
```
None
```
查看工具箱别名5.3-114
```
None
```

使用脚本执行地理处理工具5.4-116

import arcpy

in_features = "c:/ArcpyBook/data/CityOfSanAntonio.gdb/Burglary"

clip_features = "c:/ArcpyBook/Ch5/EdgewoodSD.shp"

out_feature_class = "c:/ArcpyBook/Ch5/ClpBurglary.shp"

arcpy.Clip_analysis(in_features,clip_features,out_feature_class)

讲一个工具的输出作为另一个工具的输入5.5-119

import arcpy

arcpy.env.workspace = "c:/ArcpyBook/data/TravisCounty"

try:

    # Buffer areas of impact around major roads

    streams = "Streams.shp"

    streamsBuffer = "StreamsBuffer"

    distance = "2640 Feet"

    schools2mile = "Schools.shp"

    schoolsLyrFile = 'Schools2Mile_lyr'

    arcpy.Buffer_analysis(streams, streamsBuffer, distance,'FULL','ROUND','ALL')

    # Make a layer

    arcpy.MakeFeatureLayer_management(schools2mile, schoolsLyrFile)

    arcpy.SelectLayerByLocation_management(schoolsLyrFile, 'intersect', streamsBuffer)

except Exception as e:

    print e.message

第6章创建自定义地理处理工具

创建自定义地理处理工具
创建Python工具箱

创建自定义地理处理工具6.2-123

#Script to Import data to a feature class within a geodatabase

import arcpy, os

try:

    outputFC = arcpy.GetParameterAsText(0)

    fClassTemplate = arcpy.GetParameterAsText(1)

    f = open(arcpy.GetParameterAsText(2),'r')

    arcpy.CreateFeatureclass_management(os.path.split(outputFC)[0], os.path.split(outputFC)[1],"point",fClassTemplate)

    lstFires = f.readlines()

    with arcpy.da.InsertCursor(outputFC) as cur:

        cntr = 1

        for fire in lstFires:

            if 'Latitude' in fire:

                continue

            vals = fire.split(",")

            latitude = float(vals[0])

            longitude = float(vals[1])

            confid = int(vals[2])

            pnt = arcpy.Point(longitude, latitude)

            feat = cur.newRow()

            feat.shape = pnt

            feat.setValue("CONFIDENCEVALUE", confid)

            cur.insertRow(feat)

            arcpy.AddMessage("Record number" + str(cntr) + "written to feature class")

            cntr = cntr + 1

except:

    print arcpy.GetMessages()

finally:

    f.close()

创建Python工具箱6.3-139

import arcpy

import requests

import json

class Toolbox(object):

    def __init__(self):

        """Define the toolbox (the name of the toolbox is the name of the

        .pyt file)."""

        self.label = "Toolbox"

        self.alias = ""

        # List of tool classes associated with this toolbox

        self.tools = [USGSDownload]

class USGSDownload(object):

    def __init__(self):

        """Define the tool (tool name is the name of the class)."""

        self.label = "USGS Download"

        self.description = "Download from USGS ArcGIS Server instance"

        self.canRunInBackground = False

    def getParameterInfo(self):

        """Define parameter definitions"""

        # First parameter

        param0 = arcpy.Parameter(

            displayName="ArcGIS Server Wildfire URL",

            name="url",

            datatype="GPString",

            parameterType="Required",

            direction="Input")

        param0.value = "http://wildfire.cr.usgs.gov/arcgis/rest/services/geomac_dyn/MapServer/0/query"

        # Second parameter

        param1 = arcpy.Parameter(

            displayName="Output Feature Class",

            name="out_fc",

            datatype="DEFeatureClass",

            parameterType="Required",

            direction="Input")

        params = [param0, param1]

        return params

    def isLicensed(self):

        """Set whether tool is licensed to execute."""

        return True

    def updateParameters(self, parameters):

        """Modify the values and properties of parameters before internal

        validation is performed.  This method is called whenever a parameter

        has been changed."""

        return

    def updateMessages(self, parameters):

        """Modify the messages created by internal validation for each tool

        parameter.  This method is called after internal validation."""

        return

    def execute(self, parameters, messages):

        inFeatures = parameters[0].valueAsText

        outFeatureClass = parameters[1].valueAsText

        agisurl = inFeatures

        payload = { 'where': 'acres > 5','f': 'pjson', 'outFields': 'latitude,longitude,fire_name,acres'}

        r = requests.get(inFeatures, params=payload)

        decoded = json.loads(r.text)

        with arcpy.da.InsertCursor(outFeatureClass, ("SHAPE@XY", "NAME", "ACRES")) as cur:

            cntr = 1

            for rslt in decoded['features']:

                fireName = rslt['attributes']['fire_name']

                latitude = rslt['attributes']['latitude']

                longitude = rslt['attributes']['longitude']

                acres = rslt['attributes']['acres']

                cur.insertRow([(longitude,latitude),fireName, acres])

                arcpy.AddMessage("Record number: " + str(cntr) + " written to feature class")

                cntr = cntr + 1

第7章查询和选择数据

构造正确的属性查询语句
创建要素图层和表现层
使用Select Layer by Attribute 工具选择要素和行
使用Select Layer by Location 工具选择要素
结合空间查询和属性选择要素

构造正确的属性查询语句7.2-149
```
None
```

创建要素图层和表现层7.3-154

import arcpy

arcpy.env.workspace = "c:/ArcpyBook/data/CityOfSanAntonio.gdb"

try:

    tView = arcpy.MakeTableView_management("Crime2009Table","Crime2009TView")

except Exception as e:

    print e.message

使用Select Layer by Attribute 工具选择要素和行7.4-158

import arcpy

arcpy.env.workspace = "c:/ArcpyBook/data/CityOfSanAntonio.gdb"

try:

    qry = '"SVCAREA" = \'North\''

    flayer = arcpy.MakeFeatureLayer_management("Burglary","Burglary_Layer")

    arcpy.SelectLayerByAttribute_management(flayer, "NEW_SELECTION", qry)

    cnt = arcpy.GetCount_management(flayer)

    print "The number of selected records is: " + str(cnt)

except Exception as e:

    print e.message

使用Select Layer by Location 工具选择要素7.5-161

import arcpy

arcpy.env.workspace = "c:/ArcpyBook/data/CityOfSanAntonio.gdb"

try:

    flayer = arcpy.MakeFeatureLayer_management("Burglary","Burglary_Layer")

    arcpy.SelectLayerByLocation_management(flayer, "COMPLETELY_WITHIN", "c:/ArcpyBook/Ch7/EdgewoodSD.shp")

    cnt = arcpy.GetCount_management(flayer)

    print("The number of selected records is: " + str(cnt))

except Exception as e:

    print(e.message)

import arcpy

arcpy.env.workspace = "c:/ArcpyBook/data/CityOfSanAntonio.gdb"

try:

    flayer = arcpy.MakeFeatureLayer_management("Burglary","Burglary_Layer")

    arcpy.SelectLayerByLocation_management (flayer, "WITHIN_A_DISTANCE", "c:/ArcpyBook/Ch7/EdgewoodSD.shp","1 MILES")

    cnt = arcpy.GetCount_management(flayer)

    print("The number of selected records is: " + str(cnt))

except Exception as e:

    print(e.message)

import arcpy

arcpy.env.workspace = "c:/ArcpyBook/data/CityOfSanAntonio.gdb"

try:

    flayer = arcpy.MakeFeatureLayer_management("Burglary","Burglary_Layer")

    arcpy.SelectLayerByLocation_management (flayer, "WITHIN_A_DISTANCE", "c:/ArcpyBook/Ch7/EdgewoodSD.shp","1 MILES")

    arcpy.CopyFeatures_management(flayer, "c:/ArcpyBook/Ch7/EdgewoodBurglaries.shp")

    ##cnt = arcpy.GetCount_management(flayer)

    ##print "The number of selected records is: " + str(cnt)

except Exception as e:

    print(e.message)

结合空间查询和属性选择要素7.6-165

import arcpy

arcpy.env.workspace = "c:/ArcpyBook/data/CityOfSanAntonio.gdb"

try:

	qry = '"DOW" = \'Mon\''

	flayer = arcpy.MakeFeatureLayer_management("Burglary","Burglary_Layer")

	arcpy.SelectLayerByLocation_management (flayer, "COMPLETELY_WITHIN", "c:/ArcpyBook/Ch7/EdgewoodSD.shp")

	arcpy.SelectLayerByAttribute_management(flayer, "SUBSET_SELECTION", qry)

	cnt = arcpy.GetCount_management(flayer)

	print("The total number of selected records is: " + str(cnt))

except Exception as e:

	print(e.message)

第8章在要素类和表中使用ArcPy数据访问模块

使用SearchCursor检索要素类中的要素
使用where字句筛选记录
使用几何令爱改进游标性能
使用InsertCursor插入行
使用UpdateCursor更新行
使用UpdateCursor删除行
在编辑会话中插入和更新行
读取要素类中的几何信息
使用Walk()遍历目录

使用SearchCursor检索要素类中的要素8.2-171

import arcpy.da

arcpy.env.workspace = "c:/ArcpyBook/Ch8"

with arcpy.da.SearchCursor("Schools.shp",("Facility","Name")) as cursor:

    for row in sorted(cursor):

        print("High school name: " + row[1])

使用where字句筛选记录8.3-173

import arcpy.da

arcpy.env.workspace = "c:/ArcpyBook/Ch8"

with arcpy.da.SearchCursor("Schools.shp",("Facility","Name"), '"FACILITY" = \'HIGH SCHOOL\'') as cursor:

        for row in sorted(cursor):

                print("School name: " + row[1])

使用几何令爱改进游标性能8.4-174

import arcpy.da

import time

arcpy.env.workspace = "c:/ArcpyBook/Ch8"

start = time.clock()

with arcpy.da.SearchCursor("coa_parcels.shp",("PY_FULL_OW","SHAPE@XY")) as cursor:

    for row in cursor:

        print("Parcel owner: {0} has a location of: {1}".format(row[0], row[1]))

elapsed = (time.clock() - start)

print("Execution time: " + str(elapsed))

使用InsertCursor插入行8.5-178

import arcpy

import os

arcpy.env.workspace = "C:/ArcpyBook/Ch8/WildfireData/WildlandFires.mdb"

f = open("C:/ArcpyBook/Ch8/WildfireData/NorthAmericaWildfires_2007275.txt","r")

lstFires = f.readlines()

try:

        with arcpy.da.InsertCursor("FireIncidents",("SHAPE@XY","CONFIDENCEVALUE")) as cur:

                cntr = 1

                for fire in lstFires:

                        if 'Latitude' in fire:

                                continue

                        vals = fire.split(",")

                        latitude = float(vals[0])

                        longitude = float(vals[1])

                        confid = int(vals[2])

                        rowValue = [(latitude,longitude),confid]

                        cur.insertRow(rowValue)

                        print("Record number " + str(cntr) + " written to feature class")

                        cntr = cntr + 1

except Exception as e:

        print(e.message)

finally:

    f.close()

使用UpdateCursor更新行8.6-183

import arcpy

arcpy.env.workspace = "C:/ArcpyBook/Ch8/WildfireData/WildlandFires.mdb"

try:

	#create a new field to hold the values

	arcpy.AddField_management("FireIncidents","CONFID_RATING","TEXT","10")

	print("CONFID_RATING field added to FireIncidents")

	with arcpy.da.UpdateCursor("FireIncidents",("CONFIDENCEVALUE","CONFID_RATING")) as cursor:

		cntr = 1

		for row in cursor:

			# update the confid_rating field

			if row[0] <= 40:

				row[1] = 'POOR'

			elif row[0] > 40 and row[0] <= 60:

				row[1] = 'FAIR'

			elif row[0] > 60 and row[0] <= 85:

				row[1] = 'GOOD'

			else:

				row[1] = 'EXCELLENT'

			cursor.updateRow(row)

			print("Record number " + str(cntr) + " updated")

			cntr = cntr + 1

except Exception as e:

    print(e.message)

使用UpdateCursor删除行8.7-187

import arcpy

import os

arcpy.env.workspace = "C:/ArcpyBook/Ch8/WildfireData/WildlandFires.mdb"

try:

	with arcpy.da.UpdateCursor("FireIncidents",("CONFID_RATING"),'[CONFID_RATING] = \'POOR\'') as cursor:

		cntr = 1

		for row in cursor:

			cursor.deleteRow()

			print("Record number " + str(cntr) + " deleted")

			cntr = cntr + 1

except Exception as e:

    print(e.message)

在编辑会话中插入和更新行8.8-189

import arcpy

import os

arcpy.env.workspace = "C:/ArcpyBook/Ch8/WildfireData/WildlandFires.mdb"

try:

	edit = arcpy.da.Editor('C:/ArcpyBook/Ch8/WildfireData/WildlandFires.mdb')

	edit.startEditing(True)

	with arcpy.da.UpdateCursor("FireIncidents",("CONFIDENCEVALUE","CONFID_RATING")) as cursor:

		cntr = 1

		for row in cursor:

			# update the confid_rating field

			if row[0] > 40 and row[0] <= 60:

				row[1] = 'GOOD'

			elif row[0] > 60 and row[0] <= 85:

				row[1] = 'BETTER'

			else:

				row[1] = 'BEST'

			cursor.updateRow(row)

			print("Record number " + str(cntr) + " updated")

			cntr = cntr + 1

	edit.stopEditing(True)

except Exception as e:

    print(e.message)

读取要素类中的几何信息8.9-193

import arcpy

infc = "c:/ArcpyBook/data/CityOfSanAntonio.gdb/SchoolDistricts"

# Enter for loop for each feature

for row in arcpy.da.SearchCursor(infc, ["OID@", "SHAPE@"]):

    # Print the current multipoint's ID

    print("Feature {0}:".format(row[0]))

    partnum = 0

    # Step through each part of the feature

    #

    for part in row[1]:

        # Print the part number

        #

        print("Part {0}:".format(partnum))

        # Step through each vertex in the feature

        #

        for pnt in part:

            if pnt:

                # Print x,y coordinates of current point

                #

                print("{0}, {1}".format(pnt.X, pnt.Y))

            else:

                # If pnt is None, this represents an interior ring

                #

                print("Interior Ring:")

        partnum += 1

使用Walk()遍历目录8.10-195

import arcpy.da as da

import os

print("os walk")

for dirpath, dirnames, filenames in os.walk(os.getcwd()):

    for filename in filenames:

        print(filename)

print("arcpy da walk")

for dirpath, dirnames, filenames in da.Walk(os.getcwd(),datatype="FeatureClass"):

    for filename in filenames:

        print(os.path.join(dirpath, filename))

第9章获取GIS数据的列表和描述

使用ArcPy列表函数
获取要素类或表中的字段列表
使用Describe()函数返回要素类的描述性信息
使用Describe()函数返回栅格图像的描述性信息

使用ArcPy列表函数9.2-199

import arcpy

arcpy.env.workspace = "C:/ArcpyBook/data/CityOfSanAntonio.gdb"

fcList = arcpy.ListFeatureClasses()

for fc in fcList:

    print(fc)

import arcpy

arcpy.env.workspace = "C:/ArcpyBook/data/CityOfSanAntonio.gdb"

fcList = arcpy.ListFeatureClasses("C*")

for fc in fcList:

    print(fc)

import arcpy

arcpy.env.workspace = "C:/ArcpyBook/data/CityOfSanAntonio.gdb"

fcList = arcpy.ListFeatureClasses("C*","polygon")

for fc in fcList:

    print(fc)

获取要素类或表中的字段列表9.3-202

import arcpy

arcpy.env.workspace = "C:/ArcpyBook/data/CityOfSanAntonio.gdb"

try:

    fieldList = arcpy.ListFields("Burglary")

    for fld in fieldList:

        print("%s is a type of %s with a length of %i" % (fld.name, fld.type, fld.length))

except Exception as e:

    print(e.message)

使用Describe()函数返回要素类的描述性信息9.4-204

import arcpy

arcpy.env.workspace = "c:/ArcpyBook/data/CityOfSanAntonio.gdb"

try:

    descFC = arcpy.Describe("Burglary")

    print("The shape type is: " + descFC.ShapeType)

    flds = descFC.fields

    for fld in flds:

        print("Field: " + fld.name)

        print("Type: " + fld.type)

        print("Length: " + str(fld.length))

    ext = descFC.extent

    print("XMin: %f" % (ext.XMin))

    print("YMin: %f" % (ext.YMin))

    print("XMax: %f" % (ext.XMax))

    print("YMax: %f" % (ext.YMax))

except:

    print(arcpy.GetMessages())

使用Describe()函数返回栅格图像的描述性信息9.5-208

import arcpy

arcpy.env.workspace = "c:/ArcpyBook/data"

try:

    descRaster = arcpy.Describe("AUSTIN_EAST_NW.sid")

    ext = descRaster.extent

    print("XMin: %f" % (ext.XMin))

    print("YMin: %f" % (ext.YMin))

    print("XMax: %f" % (ext.XMax))

    print("YMax: %f" % (ext.YMax))

    sr = descRaster.SpatialReference

    print(sr.name)

    print(sr.type)

except Exception as e:

    print e.message

第10章使用Add-in定制ArcGIS界面

下载并安装Python Add-in Wizard
创建按钮加载项和使用Python加载项模块
安装和测试加载项
创建工具加载项

下载并安装Python Add-in Wizard10.2-212
```
None
```

创建按钮加载项和使用Python加载项模块10.3-214

import arcpy

import pythonaddins

class ButtonClassImportWildfires(object):

    """Implementation for Wildfire_addin.button (Button)"""

    def __init__(self):

        self.enabled = True

        self.checked = False

    def onClick(self):

        layer_files = pythonaddins.OpenDialog('Select Layers to Add', True, r'C:\ArcpyBook\data\Wildfires', 'Add')

        mxd = arcpy.mapping.MapDocument('current')

        df = pythonaddins.GetSelectedTOCLayerOrDataFrame()

        if not isinstance(df, arcpy.mapping.Layer):

            for layer_file in layer_files:

                layer = arcpy.mapping.Layer(layer_file)

                arcpy.mapping.AddLayer(df, layer)

        else:

            pythonaddins.MessageBox('Select a data frame', 'INFO', 0)

安装和测试加载项10.4-223
```
None
```

创建工具加载项10.5-228

import arcpy

import pythonaddins

def __init__(self):

  self.enabled = True

  self.cursor = 3

  self.shape = 'Rectangle'

def onRectangle(self, rectangle_geometry):

    extent = rectangle_geometry

    arcpy.env.workspace = r'c:\ArcpyBook\Ch10'

    if arcpy.Exists('randompts.shp'):

        arcpy.Delete_management('randompts.shp')

        randompts = arcpy.CreateRandomPoints_management(arcpy.env.workspace,'randompts.shp',"",rectangle_geometry)

        arcpy.RefreshActiveView()

    return randompts

第11章异常识别和错误处理

默认的Python错误消息
添加Python异常处理结构(try/except/else)
使用GetMessages()函数获取工具消息
根据严重性级别筛选工具消息
测试和响应特定的错误消息

默认的Python错误消息11.2-235

import arcpy

arcpy.env.workspace = "c:/ArcpyBook/data"

arcpy.Buffer_analysis("Streams.shp","Streams_Buff.shp")

添加Python异常处理结构(try/except/else)11.3-236

import arcpy

try:

    arcpy.env.workspace = "c:/ArcpyBook/data"

    arcpy.Buffer_analysis("Streams.shp","Streams_Buff.shp")

except:

    print("Error")

使用GetMessages()函数获取工具消息11.4-238

import arcpy

try:

    arcpy.env.workspace = "c:/ArcpyBook/data"

    arcpy.Buffer_analysis("Streams.shp","Streams_Buff.shp")

except:

    print(arcpy.GetMessages())

根据严重性级别筛选工具消息11.5-240

import arcpy

try:

    arcpy.env.workspace = "c:/ArcpyBook/data"

    arcpy.Buffer_analysis("Streams.shp","Streams_Buff.shp")

except:

    print(arcpy.GetMessages(2))

测试和响应特定的错误消息116-241

import arcpy

try:

    arcpy.env.workspace = "c:/ArcpyBook/data"

    arcpy.Buffer_analysis("Streams.shp", "Streams_Buff.shp")

except:

    print("Error found in Buffer tool \n")

    errCode = arcpy.GetReturnCode(3)

    if str(errCode) == "735":

        print("Distance value not provided \n")

        print("Running the buffer again with a default value \n")

        defaultDistance = "100 Feet"

        arcpy.Buffer_analysis("Streams.shp", "Streams_Buff", defaultDistance)

        print("Buffer complete")

第12章使用Python实现ArcGIS的高级功能

ArcGIS REST API入门
使用Python构建HTTP请求并解析响应
使用ArcGIS REST API和Python获取图层信息
使用ArcGIS REST API和Python导出地图
使用ArcGIS REST API和Python查询地图服务
使用ESRI World Geocoding Service镜像地理编码
使用FieldMap和FieldMappings
使用ValueTable将多值输入到工具中

ArcGIS REST API入门12.2-245
```
None
```

使用Python构建HTTP请求并解析响应12.3-250

import requests

import json

agisurl = "http://server.arcgisonline.com/arcgis/rest/services?f=pjson"

r = requests.get(agisurl)

decoded = json.loads(r.text)

print(decoded)

#print(r.text)

使用ArcGIS REST API和Python获取图层信息12.4-254

import requests

import json

agisurl = "http://sampleserver1.arcgisonline.com/ArcGIS/rest/services/Demographics/ESRI_Census_USA/MapServer/1"

payload = { 'where': 'STATE_FIPS = \'48\' and CNTY_FIPS = \'021\'','returnCountyOnly': 'false',

            'returnIdsOnly': 'false', 'returnGeometry': 'false',

            'f': 'pjson'}

r = requests.get(agisurl, params=payload)

#r = requests.get(agisurl)

#print(r.text)

decoded = json.loads(r.text)

print("The layer name is: " + decoded['name'])

print("The xmin: " + str(decoded['extent']['xmin']))

print("The xmax: " + str(decoded['extent']['xmax']))

print("The ymin: " + str(decoded['extent']['ymin']))

print("The ymax: " + str(decoded['extent']['xmax']))

print("The fields in this layer: ")

for rslt in decoded['fields']:

    print(rslt['name'])

使用ArcGIS REST API和Python导出地图12.5-257

import requests

import json

agisurl = "http://sampleserver1.arcgisonline.com/ArcGIS/rest/services/Specialty/ESRI_StateCityHighway_USA/MapServer/export"

payload = { 'bbox': '-115.8,30.4,-85.5,50.5','size': '800,600', \

            'imageSR': '102004', 'format': 'gif', 'transparent':'false', \

            'f': 'pjson'}

r = requests.get(agisurl, params=payload)

print(r.text)

使用ArcGIS REST API和Python查询地图服务12.6-260

import requests

import json

agisurl = "http://sampleserver1.arcgisonline.com/ArcGIS/rest/services/Demographics/ESRI_Census_USA/MapServer/1/query"

payload = { 'where': 'STATE_FIPS = \'48\' and CNTY_FIPS = \'021\'','returnCountyOnly': 'false',

            'returnIdsOnly': 'false', 'returnGeometry': 'false', 'outFields':'POP2000,POP2007,BLKGRP',

            'f': 'pjson'}

r = requests.get(agisurl, params=payload)

#print(r.text)

decoded = json.loads(r.text)

#print(decoded)

for rslt in decoded['features']:

    print("Block Group: " + str(rslt['attributes']['BLKGRP']))

    print("Population 2000: " + str(rslt['attributes']['POP2000']))

    print("Population 2007: " + str(rslt['attributes']['POP2007']))

使用ESRI World Geocoding Service镜像地理编码12.7-264

import requests

import json

agisurl = "http://geocode.arcgis.com/arcgis/rest/services/World/GeocodeServer/find"

payload = { 'text': '1202 Sand Wedge, San Antonio, TX, 78258', 'f': 'pjson'}

r = requests.get(agisurl, params=payload)

decoded = json.loads(r.text)

print("The geocoded address: " + decoded['locations'][0]['name'])

print("The longitude: " + str(decoded['locations'][0]['feature']['geometry']['x']))

print("The latitude: " + str(decoded['locations'][0]['feature']['geometry']['y']))

print("The geocode score: " + str(decoded['locations'][0]['feature']['attributes']['Score']))

print("The address type: " + decoded['locations'][0]['feature']['attributes']['Addr_Type'])

使用FieldMap和FieldMappings12.8-266

import arcpy

try:

    # Local variables

    arcpy.env.workspace = r"c:\ArcpyBook\data"

    outFeatureClass = r"c:\ArcpyBook\data\AllTracts.shp"

    # Create a fieldmappings adding the three new fields

    fieldmappings = arcpy.FieldMappings()

    fldmap_STFIPS = arcpy.FieldMap()

    fldmap_COFIPS = arcpy.FieldMap()

    fldmap_TRACT = arcpy.FieldMap()

    # List all feature classes that start with 'County' and type Polygon

    fclss = arcpy.ListFeatureClasses("County*", "Polygon")

    # Create a value table with the FC to merge

    vTab = arcpy.ValueTable()

    for fc in fclss:

        fieldmappings.addTable(fc)

        fldmap_STFIPS.addInputField(fc, "STFID")

        fldmap_COFIPS.addInputField(fc, "STFID")

        fldmap_TRACT.addInputField(fc, "STFID")

        vTab.addRow(fc)

    # Set Starting and ending point from the input as well as the name of the output fields 

    # STFIPS field

    for x in range(0, fldmap_STFIPS.inputFieldCount):

        fldmap_STFIPS.setStartTextPosition(x, 0)

        fldmap_STFIPS.setEndTextPosition(x, 1)

    fld_STFIPS = fldmap_STFIPS.outputField

    fld_STFIPS.name = "STFIPS"

    fldmap_STFIPS.outputField = fld_STFIPS

    # COFIPS field

    for x in range(0, fldmap_COFIPS.inputFieldCount):

        fldmap_COFIPS.setStartTextPosition(x, 2)

        fldmap_COFIPS.setEndTextPosition(x, 4)

    fld_COFIPS = fldmap_COFIPS.outputField

    fld_COFIPS.name = "COFIPS"

    fldmap_COFIPS.outputField = fld_COFIPS

    # TRACT field

    for x in range(0, fldmap_TRACT.inputFieldCount):

        fldmap_TRACT.setStartTextPosition(x, 5)

        fldmap_TRACT.setEndTextPosition(x, 12)

    fld_TRACT = fldmap_TRACT.outputField

    fld_TRACT.name = "TRACT"

    fldmap_TRACT.outputField = fld_TRACT

    # Add fieldmaps into the fieldmappings object

    fieldmappings.addFieldMap(fldmap_STFIPS)

    fieldmappings.addFieldMap(fldmap_COFIPS)

    fieldmappings.addFieldMap(fldmap_TRACT)

    # Run the merge tool

    arcpy.Merge_management(vTab, outFeatureClass, fieldmappings)

    print("Merge completed")

except Exception as e:

    print(e.message)

使用ValueTable将多值输入到工具中12.9-273

import arcpy

try:

    arcpy.env.workspace = r'c:\ArcpyBook\data'

    vTab = arcpy.ValueTable()

    vTab.setRow (0, "5")

    vTab.setRow (1, "10")

    vTab.setRow (2, "20")

    inFeature = 'Hospitals.shp'

    outFeature = 'HospitalMBuff.shp'

    dist = vTab

    bufferUnit = "meters"

    arcpy.MultipleRingBuffer_analysis(inFeature,outFeature,dist,bufferUnit, '', 'ALL')

    print("Multi-Ring Buffer Complete")

except Exception as e:

    print(e.message)

附录A 自动化Python脚本

在命令行中运行Python脚本
使用sys.argv[]捕获命令行的输入
添加Python脚本到批处理文件
在规定的时间运行批处理文件

在命令行中运行Python脚本A.2-283
```
None
```

使用sys.argv[]捕获命令行的输入A.3-289

import arcpy

arcpy.env.workspace = "C:/ArcpyBook/data/CityOfSanAntonio.gdb"

try:

    fieldList = arcpy.ListFields("Burglary")

    for fld in fieldList:

        print "%s is a type of %s with a length of %i" % (fld.name, fld.type, fld.length)

except Exception as e:

    print(e.message)

import arcpy

import sys

wkspace = sys.argv[1]

fc = sys.argv[2]

try:

  arcpy.env.workspace = wkspace

  fields = arcpy.ListFields(fc)

  for fld in fields:

    print(fld.name)

except Exception as e:

  print(e.message)

添加Python脚本到批处理文件A.4-290

cd c:\ArcpyBook\Appendix1

python ListFields.py c:\ArcpyBook\data Burglaries_2009.shp

在规定的时间运行批处理文件A.5-292
```
None
```

附录B GIS程序员不可不知的5个Python功能

读取带分隔符的文本文件
发送电子邮件
检索FTP服务中的文件
创建ZIP文件
读取XML文件

读取带分隔符的文本文件B.2-298

f = open('c:/ArcpyBook/data/N_America.A2007275.txt','r')

for fire in f:

    lstValues = fire.split(',')

    latitude = float(lstValues[0])

    longitude = float(lstValues[1])

    confid = int(lstValues[8])

    print("The latitude is: " + str(latitude) + " The longitude is: " + str(longitude) + " The confidence value is: " + str(confid))

f.close()

发送电子邮件B.3-301

import smtplib

from email.MIMEMultipart import MIMEMultipart

from email.MIMEBase import MIMEBase

from email.MIMEText import MIMEText

from email import Encoders

import os

gmail_user = "<username>"

gmail_pwd = "<password>"

def mail(to, subject, text, attach):

   msg = MIMEMultipart()

   msg['From'] = gmail_user

   msg['To'] = to

   msg['Subject'] = subject

   msg.attach(MIMEText(text))

   part = MIMEBase('application', 'octet-stream')

   part.set_payload(open(attach, 'rb').read())

   Encoders.encode_base64(part)

   part.add_header('Content-Disposition',

           'attachment; filename="%s"' % os.path.basename(attach))

   msg.attach(part)

   mailServer = smtplib.SMTP("smtp.gmail.com", 587)

   mailServer.ehlo()

   mailServer.starttls()

   mailServer.ehlo()

   mailServer.login(gmail_user, gmail_pwd)

   mailServer.sendmail(gmail_user, to, msg.as_string())

   mailServer.close()

mail("<email to send to>",

   "Hello from python!",

   "This is an email sent with python",

   "bc_pop1996.csv")

检索FTP服务中的文件B.4-306

import ftplib

import os

import socket

HOST = 'ftp.nifc.gov'

DIRN = '/Incident_Specific_Data/2012 HISTORIC/ROCKY_MTN/Arapaho/GIS/20120629'

FILE = '20120629_0600_Arapaho_PIO_0629_8x11_land.pdf'

try:

    f = ftplib.FTP(HOST)

except (socket.error, socket.gaierror), e:

    print('ERROR: cannot reach "%s"' % HOST)

print('*** Connected to host "%s"' % HOST)

try:

    f.login()

except ftplib.error_perm:

    print('ERROR: cannot login anonymously')

    f.quit()

print('*** Logged in as "anonymous"')

try:

    f.cwd(DIRN)

except ftplib.error_perm:

    print('ERROR: cannot CD to "%s"' % DIRN)

    f.quit()

print('*** Changed to "%s" folder' % DIRN)

try:

    f.retrbinary('RETR %s' % FILE,

          open(FILE, 'wb').write)

except ftplib.error_perm:

    print('ERROR: cannot read file "%s"' % FILE)

    os.unlink(FILE)

else:

    print('*** Downloaded "%s" to CWD' % FILE)

f.quit()

创建ZIP文件B.5-310

import os

import zipfile

#create the zip file

zfile = zipfile.ZipFile("shapefiles2.zip", "w", zipfile.ZIP_STORED)

files = os.listdir("c:/ArcpyBook/data")

for f in files:

    if f.endswith("shp") or f.endswith("dbf") or f.endswith(".shx"):

        zfile.write("C:/ArcpyBook/data/" + f)

#list files in the archive

for f in zfile.namelist():

        print("Added %s" % f)

zfile.close()

import os

import zipfile

#create the zip file

zfile = zipfile.ZipFile("shapefiles2.zip", "w", zipfile.ZIP_DEFLATED)

files = os.listdir("c:/ArcpyBook/data")

for f in files:

    if f.endswith("shp") or f.endswith("dbf") or f.endswith(".shx"):

        zfile.write("C:/ArcpyBook/data/" + f)

#list files in the archive

for f in zfile.namelist():

        print("Added %s" % f)

zfile.close()

读取XML文件B.6-313

from xml.dom import minidom

xmldoc = minidom.parse("WitchFireResidenceDestroyed.xml")

childNodes = xmldoc.childNodes

eList = childNodes[0].getElementsByTagName("fire")

for e in eList:

    if e.hasAttribute("address"):

        print e.getAttribute("address")

『基于ArcGIS的Python编程秘籍（第2版）』书本源码的更多相关文章

Java编程思想第四版随书源码官方下载方法
见不少人在找net.mindview.util.Print,CSDN上有下载,收积分,以下是官网的下载方法,免费: 官网链接:http://mindview.net/ 电子书下载地址:http://w ...
Python编程导论第2版|百度网盘免费下载|新手学习
点击下方即可免费下载百度网盘免费下载:Python编程导论第2版提取码:18g5 豆瓣评论: 介绍: 本书基于MIT 编程思维培训讲义写成,主要目标在于帮助读者掌握并熟练使用各种计算技术,具备用计 ...
Python编程入门(第3版) PDF|百度网盘下载内附提取码
Python编程入门(第3版)是图文并茂的Python学习参考书,书中并不包含深奥的理论或者高级应用,而是以大量来自实战的例子.屏幕图和详细的解释,用通俗易懂的语言结合常见任务,对Python的各项基 ...
❤️❤️新生代农民工爆肝8万字，整理Python编程从入门到实践(建议收藏)已码：8万字❤️❤️
@ 目录开发环境搭建安装 Python 验证是否安装成功安装Pycharm 配置pycharm 编码规范基本语法规则保留字单行注释多行注释行与缩进多行语句数据类型空行等待用户输 ...
Python编程第四版中文上下册完整版pdf|网盘下载附提取码
点击此处下载提取码:drjh 作者简介 Mark Lutz是Python培训的世界的领先者,他是最初和最畅销的Python著作的作者,从1992年起就是Python社区的先锋人物.Mark有25年的 ...
ArcGIS Python编程案例-电子资料链接
ArcGIS Python编程案例(1)-Python语言基础 https://www.jianshu.com/p/dd90816d019b ArcGIS Python编程案例(2)-使用ArcPy编 ...
基于海龟编辑器python少儿编程
Python 少儿教程为什么要学习编程扫地机器人.物流机器人.自动泊车系统.无人超市.3D打印.微信.支付宝等等,随着人工智能时代的到来,越来越多的岗位将被机器人所替代. 所以,学习编程的最终目的 ...
Python编程规范(PEP8)
Python编程规范(PEP8) 代码布局缩进对于每一次缩进使用4个空格.使用括号.中括号.大括号进行垂直对齐,或者缩进对齐. 制表符还是空格? 永远不要将制表符与空格混合使用.Python最常用 ...
Python 编程规范-----转载
Python编程规范及性能优化 Ptyhon编程规范编码所有的 Python 脚本文件都应在文件头标上 # -*- coding:utf-8 -*- .设置编辑器,默认保存为 utf-8 格式. ...

随机推荐

The Second Week lucklyzpp
The Second Week 文件通配符模式在Linux系统中预定义的字符类 1.显示/etc目录下,以非字母开头,后面跟了一个字母以及其它任意长度任意字符的文件或目录 2.复制/etc目录下 ...
01_Keil与Proteus联合仿真的注意事项
01. 关于keil5和Proteus8的联合仿真的操作步骤,这里就不细说,给个链接,步骤差不多是最齐全的 CSDN博客链接:https://blog.csdn.net/wzk456/article/ ...
Android App性能测试之adb命令
本篇文章总结了Android App性能测试过程中常用的adb命令.通过这些adb命令,可以查看App的性能数据,为评判性能好坏作参考. CPU相关显示占用CPU最大的5个应用 adb shell ...
java基础之反射类型Type
Java在加入泛型之后,仅仅Class已经不足以描述数据的类型了,比如List<String>类型的数据,其Class的类型为List.class,但是其类型包含了泛型参数,所以java引 ...
【MyBatis】几种批量插入效率的比较
批处理数据主要有三种方式: 反复执行单条插入语句 foreach 拼接 sql 批处理一.前期准备基于Spring Boot + Mysql,同时为了省略get/set,使用了lombok,详见p ...
用Java写了一个程序，将一个Mysql库中的表，迁移到另外一个server上的Mysql库中
用Navicat做数据迁移,因为数据量比较大,迁移过过程中一个是进展不直观,另外就是cpu占用率高的时候,屏幕跟死机了一样点不动按钮,不好中断. 想了想,干脆自己写一个. 在网上找了一个sqllite ...
WPF listbox中Checkbox横向排列
<ListBox Height="220" Margin="0" ItemsSource="{Binding RightCollection}& ...
CentOS linux系统将UTC时间修改为CST时间
1.编辑时间配置文件 1 2 3 4 # vi /etc/sysconfig/clock ZONE="Asia/Shanghai" UTC=false ...
数学相关函数在PHP中的应用简介
对于数学计算来说,最常见的其实还是我们使用各种操作符的操作,比如说 +加.-减之类的.当然,PHP 中也为我们提供了一些可以方便地进行其他数学运算的操作函数.这些函数都属于 Math 扩展.这个扩展 ...
ubuntu系统安装docker
系统版本:Ubuntu 18.04 # 更新apt update # 安装依赖apt install apt-transport-https ca-certificates curl software ...

『基于ArcGIS的Python编程秘籍（第2版）』书本源码

ArcPy学习

第1章 面向ArcGIS的Python编程语言的基础

第2章 管理地图文档和图层

第3章 查找和修复丢失的数据链接

第4章 自动化地图制图和打印

第5章 使用脚本执行地理处理工具

第6章 创建自定义地理处理工具

第7章 查询和选择数据

第8章 在要素类和表中使用ArcPy数据访问模块

第9章 获取GIS数据的列表和描述

第10章 使用Add-in定制ArcGIS界面

第11章 异常识别和错误处理

第12章 使用Python实现ArcGIS的高级功能