关于arcgis的python脚本编程, shape文件出png图片问题

发布网友 发布时间：2022-04-24 09:24

我来回答

共3个回答

懂视网 时间：2022-04-18 23:11

出图是项目里常见的任务，有的项目甚至会要上百张图片，所以批量出土工具很有必要。arcpy.mapping就是ArcGIS里的出图模块，能快速完成一个出图工具。

arcpy.mapping模块里常用的类有MapDocument、DataFrame、Layer、DataDrivenPages和TextElement。

MapDocument类是地图文档（.mxd文件）对应的类。初始化参数是一个字符串，一般是.mxd文件的路径：

 mxd=arcpy.mapping.MapDocument(r"F:GeoDataChinaAreaChinaVector.mxd")

DataFrame类用于操作地图内的Data Frame（即下图的Layers），能够控制地图的范围、比例尺等。用arcpy.mapping.ListDataFrames(map_document, {wildcard})函数获取。

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

Layer类用于操作具体的图层。能够控制图斑的样式、可见性等。可以用.lyr文件的路径初始化，也可以通过arcpy.mapping.ListLayers(map_document_or_layer, {wildcard}, {data_frame})函数获取。

lyr1=arcpy.mapping.Layer(r" F:GeoDataChinaAreaProvince.lyr")

df.addLayer(lyr1)

lyr2=arcpy.mapping.ListLayer(mxd,"",df)[0]

DataDrivenPages类需要配合ArcMap中的Data Driven Pages工具使用。用于一个矢量文件内的全部或部分图斑每个出一张图的情况。

TextElement类用于操作地图上的文字，比如图名、页数。通过arcpy.mapping.ListLayoutElements (map_document, {element_type}, {wildcard})函数获取。

txtElm=arcpy.mapping.ListLayoutElements(mxd,"TEXT_ELEMENT")[0]

常见的出图模式有两种：一个矢量文件里每个图斑出一张图，一个文件夹下每个矢量文件出一张图。

每个图斑出一张图：

这种情况有Data Driven Pages工具配合最好。打开ArcMap的Customize->Toolbars->Data Driven Pages，设置好图层、名称字段、排序字段、显示范围和比例尺，保存地图。

# coding:utf-8

import arcpy

 

mxd=arcpy.mapping.MapDocument(r"F:GeoDataChinaAreaChinaVector.mxd")

for pageNum in range(1,mxd.dataDrivenPages.pageCount):

 mxd.dataDrivenPages.currentPageID=pageNum

 mapName=mxd.dataDrivenPages.pageRow.getValue(mxd.dataDrivenPages.pageNameField.name)

 print mapName

 arcpy.mapping.ExportToPNG(mxd,r"F:GeoDataChinaAreaProvince"+mapName+".png")

print 'ok'

一个文件夹下的每个矢量文件出一张图：

# coding:utf-8

import arcpy

import os

 

def GetShpfiles(shpdir):

 shpfiles=[]

 allfiles=os.listdir(shpdir)

 for file in allfiles:

  if os.path.isfile(file):

  if file.endswith('.shp'):

   shpfiles.append(file)

  else:

  shpfiles.extend(GetShpfiles(file))

 return shpfiles

 

allshps=GetShpfiles(r"F:GeoDataChinaAreaProvince")

mxd=arcpy.mapping.MapDocument(r"F:GeoDataChinaAreaChinaVector.mxd")

lyr=arcpy.mapping.ListLayer(mxd)[0]

for shp in allshps:

 paths=os.path.split(shp)

 print paths[1]

 lyr.replaceDataSource(paths[0],"SHAPEFILE_WORKSPACE",paths[1])

 arcpy.mapping.ExportToPNG(mxd,r"F:GeoDataChinaAreaProvince"+paths[1]+".png")

print 'ok'

热心网友 时间：2022-04-18 20:19

1、点类型的转栅格：

PointToRaster example 1 (Python 
window)


Converts point features to a raster dataset.

import arcpy
from arcpy import env
env.workspace = "c:/data"
arcpy.PointToRaster_conversion("ca_ozone_pts.shp", "ELEVATION", 
                                "c:/output/ca_elev", "MAXIMUM", "", 2000)

PointToRaster example 2 
(stand-alone script)

Converts point features to a raster dataset.

# Name: PointToRaster_Ex_02.py
# Description: Converts point features to a raster dataset.
# Requirements: ArcInfo

# Import system moles
import arcpy
from arcpy import env

# Set environment settings
env.workspace = "C:/data"

# Set local variables
inFeatures = "ca_ozone_pts.shp"
valField = "ELEVATION"
outRaster = "c:/output/ca_elev02"
assignmentType = "MAXIMUM"
priorityField = ""
cellSize = 2000

# Execute PointToRaster
arcpy.PointToRaster_conversion(inFeatures, valField, outRaster, 
                               assignmentType, priorityField, cellSize)

2、面类型的转栅格：

PolygonToRaster example 1 
(Python window)

Converts polygon features to a raster dataset.

import arcpy
from arcpy import env
env.workspace = "c:/data"
arcpy.PolygonToRaster_conversion("ca_counties.shp", "NAME", 
                                 "c:/output/ca_counties.img", 
                                 "MAXIMUM_AREA", "MALES", 0.25)

PolygonToRaster example 2 
(stand-alone script)

Converts polygon features to a raster dataset.

# Name: PolygonToRaster_Ex_02.py
# Description: Converts polygon features to a raster dataset.
# Requirements: ArcInfo

# Import system moles
import arcpy
from arcpy import env

# Set environment settings
env.workspace = "C:/data"

# Set local variables
inFeatures = "ca_counties.shp"
valField = "NAME"
outRaster = "c:/output/ca_counties"
assignmentType = "MAXIMUM_AREA"
priorityField = "MALES"
cellSize = 0.5

# Execute PolygonToRaster
arcpy.PolygonToRaster_conversion(inFeatures, valField, outRaster, 
                                 assignmentType, priorityField, cellSize)

3、线类型转栅格

PolylineToRaster example 1 
(Python window)

Converts polyline features to a raster dataset.

import arcpy
from arcpy import env
env.workspace = "c:/data"
arcpy.PolylineToRaster_conversion("roads.shp", "CLASS", "c:/output/roads.img",
                                   "MAXIMUM_COMBINED_LENGTH", "LENGTH", 30) 

PolylineToRaster example 2 
(stand-alone script)

Converts polyline features to a raster dataset.

# Name: PolylineToRaster_Ex_02.py
# Description: Converts polyline features to a raster dataset.
# Requirements: ArcInfo

# Import system moles
import arcpy
from arcpy import env

# Set environment settings
env.workspace = "C:/data"

# Set local variables
inFeatures = "roads.shp"
valField = "CLASS"
outRaster = "c:/output/roads.tif"
assignmentType = "MAXIMUM_COMBINED_LENGTH"
priorityField = "LENGTH"
cellSize = 30

# Execute PolylineToRaster
arcpy.PolylineToRaster_conversion(inFeatures, valField, outRaster, 
                                  assignmentType, priorityField, cellSize)

追问就只是使用了工具箱生成的python代码吗？
1、你自己有没有试过生成的python能不能运行？
2、我要生成的是mxd文件，然后导出png文件，并不需要栅格文件

追答咱只是提供点线索。
你可以试试，毕竟每个人都有自己的事，如果恰好做过，而且热心的人肯帮助你那是你运气好。祝你好运。

热心网友 时间：2022-04-18 21:37

上段代码运行出错就是因为程序认为No是个字符串啊，把No直接改成1运行没有问题 Traceback (most recent call last): File "C:\Documents and Settings\Administrator\桌面\a.py", line 18, in select("XYSheet1.shp","XYSheet1_Select1.shp","\"矿区编号\" =No") File "C:\Documents and Settings\Administrator\桌面\a.py", line 16, in select arcpy.Select_analysis(in_features, out_feature_class, where_clause) File "E:\ArcGIS\Desktop10.0\arcpy\arcpy\analysis.py", line 88, in Select raise eExecuteError: ERROR 999999: 执行函数时出错。指定了不存在的列。指定了不存在的列。执行(Select)失败。