arcgispython脚本
㈠ 关于arcgis的python脚本编程， shape文件出png图片问题
1、点类型的转栅格：
PointToRasterexample1(Python
window)
.
importarcpy
fromarcpyimportenv
env.workspace="c:/data"
arcpy.PointToRaster_conversion("ca_ozone_pts.shp","ELEVATION",
"c:/output/ca_elev","MAXIMUM","",2000)
PointToRasterexample2
(stand-alonescript)
.
#Name:PointToRaster_Ex_02.py
#Description:.
#Requirements:ArcInfo
#Importsystemmoles
importarcpy
fromarcpyimportenv
#Setenvironmentsettings
env.workspace="C:/data"
#Setlocalvariables
inFeatures="ca_ozone_pts.shp"
valField="ELEVATION"
outRaster="c:/output/ca_elev02"
assignmentType="MAXIMUM"
priorityField=""
cellSize=2000
#ExecutePointToRaster
arcpy.PointToRaster_conversion(inFeatures,valField,outRaster,
assignmentType,priorityField,cellSize)
2、面类型的转栅格：
PolygonToRasterexample1
(Pythonwindow)
.
importarcpy
fromarcpyimportenv
env.workspace="c:/data"
arcpy.PolygonToRaster_conversion("ca_counties.shp","NAME",
"c:/output/ca_counties.img",
"MAXIMUM_AREA","MALES",0.25)
PolygonToRasterexample2
(stand-alonescript)
.
#Name:PolygonToRaster_Ex_02.py
#Description:.
#Requirements:ArcInfo
#Importsystemmoles
importarcpy
fromarcpyimportenv
#Setenvironmentsettings
env.workspace="C:/data"
#Setlocalvariables
inFeatures="ca_counties.shp"
valField="NAME"
outRaster="c:/output/ca_counties"
assignmentType="MAXIMUM_AREA"
priorityField="MALES"
cellSize=0.5
#ExecutePolygonToRaster
arcpy.PolygonToRaster_conversion(inFeatures,valField,outRaster,
assignmentType,priorityField,cellSize)
3、线类型转栅格
PolylineToRasterexample1
(Pythonwindow)
.
importarcpy
fromarcpyimportenv
env.workspace="c:/data"
arcpy.PolylineToRaster_conversion("roads.shp","CLASS","c:/output/roads.img",
"MAXIMUM_COMBINED_LENGTH","LENGTH",30)
PolylineToRasterexample2
(stand-alonescript)
.
#Name:PolylineToRaster_Ex_02.py
#Description:.
#Requirements:ArcInfo
#Importsystemmoles
importarcpy
fromarcpyimportenv
#Setenvironmentsettings
env.workspace="C:/data"
#Setlocalvariables
inFeatures="roads.shp"
valField="CLASS"
outRaster="c:/output/roads.tif"
assignmentType="MAXIMUM_COMBINED_LENGTH"
priorityField="LENGTH"
cellSize=30
#ExecutePolylineToRaster
arcpy.PolylineToRaster_conversion(inFeatures,valField,outRaster,
assignmentType,priorityField,cellSize)
㈡ 向精通ArcGIS——python脚本开发的同学请教
上段代码运行出错就是因为程序认为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)失败。
㈢ 怎么用python开发arcgis
在Pyhton写的一些代码，用户交互不方便，用户体验比较差，不方便重用。在ArcGIS中可以将用写的Python代码导入到ToolBox中，这样用起来就比较方便了。这里用按要素裁剪栅格的Python来演示如何导入ArcGIS中。代码如下：
import sys
reload(sys)
sys.setdefaultencoding( "utf-8" )
import arcpy
import string
from arcpy.sa import *
try:
raster = arcpy.GetParameterAsText(0) #要裁剪的栅格
clip_feat = arcpy.GetParameterAsText(1) #裁剪要素类
field = arcpy.GetParameterAsText(2) #命名字段
outworkspace = arcpy.GetParameterAsText(3) #命名字段裁剪后输出目录
for row in arcpy.SearchCursor(clip_feat):
mask=row.getValue("Shape")
outPath=outworkspace+"\\"+str(row.getValue(field))
outExtractByMask = ExtractByMask(raster,mask)
outExtractByMask.save(outPath)
except arcpy.ExecuteError:
print arcpy.GetMessages() 在用户工具箱中新建工具箱（在系统工具箱中不能新建），在工具箱右键，添加脚本；
输入名称、标签、描述等信息。下一步，选择脚本文件。（这里需要注意的是：一定要勾选"存储相对路径名"这个选项）
设置参数，这是最重要的一步。其中参数即代码中GetParameterAsText（n），并选择合适的数据类型，在参数属性中也可以进行相关设置。如果输出栅格名称按裁剪要素中的某个字段，需要设置"获取息"属性，还要进行过滤一下。
设置好一切参数过后，就可以来进行测试了。在本机上测试通！
接下来便是怎样移植到其他电脑了的问题了。前面我说说过要存储为相对路径，否则，将报错00576：脚本工具使用的脚本未处于所需位置。创建脚本工具时，有一个选项用于存储相对路径名（而不是绝对路径名）。设置此选项后，脚本的相对位置和保存脚本工具的工具箱必须保持不变。执行工具时移动两者中的任何一个都会出现此错误。
但改为相对路径后，在其他电脑上也能正常运行，但弹出了警告窗口。原因是客户机上没有勾选扩展模块相应的功能。这个示例需要勾选空间分析模板，因为使用了按掩膜提取这个工具。这个脚本实现的详情参考：使用Python脚本批量裁切栅格
发布给他人，涉及到个人知识产权的问题，怎样让他人使用工具，又不能看到脚本代码？如果编辑调试完成了，在脚本工具右键可导入脚本。
设置密码后，即可。
㈣ arcgis中脚本如何看源代码
aml文件直接用记事本打开。
vba直接在arcmap中，在工具条上右键单击，查看源代码即可。
python文件的源文件是*.py，直接用记事本打开。
toolbox里面也有一些工具是python的脚本，可以看源代码的。
有一些是封装过的无法看源代码。
㈤ arcgis 怎么用python做统计分析
方法/步骤
1
打开arcmap软件，点击工具栏上的ArcToolbox工具箱，在工具箱里面选择要进行数据处理的工具（工具的选择视数据处理的内容而定），本文选择了“Define Projection”定义投影的工具进行说明。
双击此工具，在弹出的工具窗口点击右下角的“Show Help>>”按钮。
2
点击“Tool Help”按钮，此时就打开了软件的帮助文档并定位到了此工具（define projection）对应的页面，往下拉动滚动条，定位到此工具的python脚本的内容：
Syntax 介绍了此工具脚本的语法
Code Sample 介绍了脚本使用的例子
拷贝脚本使用例子里面的代码。
3
在软件的主界面点击“Geoprocessing”菜单-->点击“Python”菜单项，即出现内嵌在arcigis中的python工具窗口。
python工具窗口的左边是代码编辑区域，右边是信息提示区域。
4
把第2步拷贝的代码语句粘贴到python工具窗口的左边区域，修改一下工具脚本的参数（视具体工具和数据而定），连续点击两下回车键。
5
等待脚本程序的执行完成，完成之后在桌面右下角会弹出提示信息，同时pythong工具窗口的右边区域也会有具体的工具脚本执行信息。
6
执行完成之后，在arcmap加载处理过的数据，右键查看属性，可以看到通过此工具脚本（define projection）进行投影定义的数据已经具有了投影参数。
http://jingyan..com/article/eae07827ba96a51fec548513.html