Questions tagged «arcpy»

由Esri创建的网站包,用于将Python桌面(ArcMap和ArcGIS Pro)的地理处理,地图绘制等内容暴露给Python。

3
python.exe已停止工作
大约18个月前,一个已经离开的人写了一个python脚本。然后产生了所需的输出。我被要求再次运行它,但是使用不同的(更高分辨率)数据输入。输入数据集已分为20个子集,每个子​​集约有2700个数据点。但是,在处理了大约300个数据点(范围从295到306,并且并不总是在同一条记录上失败)之后,脚本崩溃(“ python.exe已停止工作”)。 由于脚本年代久远,它是使用arcgisscripting而不是arcpy编写的。广泛地,它使用游标执行以下操作: 对于给定的点,以60分钟的旅行时间作为起点计算成本距离(使用gp.CostDistance_sa)。 调用gp.ExtractValuesToPoints_sa提取每个数据点上的所有单个值,并将要素类输出到文件地理数据库。 读取上面b)中创建的要素类,然后将值写入CSV文件(忽略带有“无数据”的任何点(值-9999))。 对输入文件中的所有剩余数据点重复1、2和3。 处理时间约为 每个数据点平均1分钟。以下是一些相关的技术规格: 该PC具有运行于2.20GHz的四核Intel i7-2720QM CPU和运行Windows 7(64位)的8GB RAM。 Python版本为2.6.6(shell在win32上还声明“ [MSC v,1500 32 bit(Intel)]”)。 还安装了ArcMap 10.0(SP4)。 我试过在另一台PC上运行它(到目前为止没有崩溃)。当前,该作业在较旧的PC上运行成功(但速度较慢),并且已达到419条记录而没有崩溃。该机器的相关规格为: 英特尔酷睿2 DUO E7500处理器,运行频率为2.93 GHz,具有4 GB RAM和64位Windows 7。 Python版本2.5.1(shell在win32上还声明了“ [[MSC v,1310 32位(Intel)]]”)。 已安装ArcMap 9.3(未提及任何Service Pack)。 有人可以提供一些有关脚本为何会工作一段时间然后崩溃的建议,以及如何解决该问题的建议吗? (到目前为止)出现了另一台PC处理该脚本的事实,这表明有些“环境”问题。 作为更新,运行ARCGIS 9.3的PC仍在成功处理数据,并且已经处理了1300个数据点(并且仍在计数)。一位同事还在运行ARCGIS 10.1的PC上运行数据-在两次分别记录267条记录后,它崩溃了。尽管不是结论性的,但通用线程似乎是Arc 9.3将处理数据,而Arc 10.x将不处理数据。

2
如何实现二元Ripley的K函数?
随附的图像显示了一个森林间隙,赤松表示为圆圈,白松表示为十字架。我有兴趣确定两种松树之间是否存在正关联或负关联(即它们是否在同一地区生长)。我知道R spatstat套件中的Kcross和Kmulti。但是,由于我需要分析50个空白并且比R更熟悉python编程,因此我想找到一种使用ArcGIS和python的迭代方法。我也愿意接受R解决方案。 如何实现二元Ripley的K函数?

1
'#'在python'if'中的作用:
我最近将ArcGIS模型导出到python,并正在查看脚本,并被以下行感到困惑: if TRU_W_DatableFeatures == '#' or not TRU_W_DatableFeatures: “#”的作用是什么?
9 arcpy 


2
如何计算表格中的零值并将其写入新字段?
这是一个非常简单的任务,但是我无法理解正确的语法。 我有一个shapefile,其属性类似于以下内容: FID Shape FIELD1 FIELD2 FIELD3 ... 0 Polygon 0 1 0 1 Polygon 3 0 7 2 Polygon 3 4 7 ... 字段的数量及其名称总是不同的。 我需要创建一个新字段(将其命名为NUM),并计算每行中的零个数。 输出示例: FID Shape FIELD1 FIELD2 FIELD3 NUM 0 Polygon 0 1 0 2 1 Polygon 3 0 7 1 2 Polygon 3 4 7 …


3
使用ArcPy检索起点和终点坐标?[关闭]
关闭。这个问题是题外话。它当前不接受答案。 想改善这个问题吗? 更新问题,使其成为地理信息系统堆栈交换的主题。 10个月前关闭。 如何使用ArcPy检索折线要素类的起点和终点坐标? 我希望将段标识符传递给子例程,并使其传回开始和结束坐标。字段计算器方法对我不起作用,因为我需要在其中无法执行的其他计算的值。(我也不想更改数据以将这些坐标存储为属性。)我正在尝试计算“中心出”寻址方案的地址范围。地址值取决于到“县中心”的距离。
9 arcpy  point 

1
使用Python从ArcMap删除图层
问题: 我正在尝试遍历所有图层,以查找名为的图层"CADAnnotation"。 如果该层存在,则从mxd中删除该层 笔记: 从独立脚本运行(即不在Arcmap中) "CADAnnotation" 数据类型是CAD注释要素类 "CADAnnotation" 不是在地理数据库中,而是从AutoCAD .dwg创建的 "CADAnnotation" 在名为 "ACAD" 如果"ACAD"可以删除组图层,那也可以删除它"CADAnnotation"。 到目前为止的代码: for item in mxds: print (item) mxd = arcpy.mapping.MapDocument(item) df=arcpy.mapping.ListDataFrames(mxd,"Project Area")[0] for lyr in arcpy.mapping.ListLayers(mxd, "*",df): if lyr.name == "CADAnnotation": print(lyr.dataSource) arcpy.Delete_management("CADAnnotation") print("Layer Deleted") else: pass 代码注释: 我可以找到图层没问题 该行arcpy.Delete_management("CADAnnotation")不起作用会引发错误。 题: 如何删除"CADAnnotation"和/或"ACAD"分组图层?
9 arcpy  arcmap  layers 

2
使用arcpy为x个邻居生成近表
我正在尝试为以下距离捕获3种场景: 换乘地铁站,具有2个或更多的相邻站。即,所讨论的站点连接2条或更多的主要路线,并具有2条或更多的相邻站点。 终点站地铁站,只有1个相邻站。这是行尾的车站。 嵌入式地铁站,其中有2个相邻站,两种方法之一。 我正在尝试计算一个可能称为“相邻站点之间的平均距离”的值 该arcpy.GenerateNearTable_analysis()可处理两种选择:距离最接近的特征,并且距离所有功能之间。 有没有人有解决这些情况的聪明方法?注意,在属性表的“ StationType”字段下,每个站都被指定为“交换”,“端子”或“内联”。 添加: 这是一些基于@whuber在评论中建议的伪代码。我还没有时间弄清楚这个问题,因此,如果有人想刺中它,您将获得复选标记!;) 我看了一下NetworkX库,它似乎可以按我的意愿工作。 给定图: A —― B ―― C ―― D | E 以及节点和链接: Nodes = ["A", "B", "C", "D", "E"] Links = [("A", "B"), ("B", "C"), ("C", "D"), ("B", "E")] def myFunction(node): identify the links that node belongs to count the …
9 python  arcpy  network 

2
创建多个MXD中所有数据集的符号系统摘要表?
这与Mapinfo TAB中的以下“ 保留ESRI符号”有关。作为第一步(在arcpy中),我想访问所有数据集的符号系统,并创建一个文档,其中包含如下图所示的符号系统详细信息。 这是否有可能(我应该花点时间来解决),或者如果不可能,我将不得不在数百个MXD上手动进行破解(:- 要分解这些步骤-希望我可以使过程的至少一部分自动化。 滚动浏览mxd并获取打开的数据集没有问题 lyrFile = arcpy.mapping.Layer,如果lyrFile.visible == True: 找到数据集后,应创建一个.lyr文件(因为大多数只是shp / gdb文件),然后创建AttribXXX = lyrFile.xxx XXX和XXX在哪里 BRIGHTNESS = lyrFile.brightness (格式如下) CONTRAST —contrast DATASETNAME —datasetName DATASOURCE —dataSource DESCRIPTION —description LABELCLASSES —labelClasses LONGNAME —longName NAME —name SHOWLABELS —showLabels TRANSPARENCY —transparency WORKSPACEPATH —A layer's workspace or connection file path 问题是我找不到能获得RGB,线/点/多边形类型,线宽,字体等的任何东西。如果可以找到,则下一步 将上述每个值写入CSV(简单)...例如(必须扩展) f …

1
在Python IDLE中导入arcpy时出现问题
我想在我的python代码中使用arcpy。我可以通过ARCGIS桌面python控制台导入它。但是我无法将其导入Python IDLE。我收到以下错误 Traceback (most recent call last): File "C:\Python26\lib\site-packages\win32com\client\dynamic.py", line 246, in __getitem__raise TypeError("This object does not support enumeration") 我在我的python路径中提到了以下文件夹。 C:\Program Files\ArcGIS\Desktop10.0\arcpy; C:\Program Files\ArcGIS\Desktop10.0\arcpy\arcpy; C:\Python26\ArcGIS10.0\Tools\Scripts 我发现有人提到同样的问题。 导入arcpy会产生“ TypeError:此对象不支持枚举” 用户说这是由于以下原因: 我在同一目录中为Arcgis 9和10编写了代码。关于导入arcpy中的初始化逻辑的一些知识,找到了用于将arcgis 9.3和python 2.6放在一起时使用的自定义arcgisscripting.py。 我只有ARCGIS 10,我不知道arcpy如何查找其模块。

2
如何使用arcpy将向量要素裁剪到现有栅格数据集的范围
我有一个裁剪的栅格数据集,上面覆盖了折线(来自轮廓工具)(请参见附件)。我正在尝试将折线裁剪到栅格数据集的范围,这似乎很简单,但是我尝试了很多事情,但似乎无法使其工作。我最终将需要使用arcpy对此进行编码,但是我认为现在,只是朝着正确的方向轻推一下哪种工具可以完成此操作,这对我现在有帮助。 谢谢。 PS我正在使用ArcInfo 10.0


5
python.multiprocessing和“致命错误(INFADI)缺失目录”
尝试使用arcpy进行多处理时,偶尔会遇到此错误: FATAL ERROR (INFADI) MISSING DIRECTORY 我不知道是什么触发了此错误,它使python进程崩溃,无法对其进行追溯。它是在写入冗长的声音模型的最终栅格输出时发生的。 有时会伴有错误 Unable to write BND file for %TEMP%\ras#### 其中%Temp是正确解析的,而####是一些随机的4位数字。这是不寻常的,因为每个进程都有自己的工作区,应在其中写入大多数文件。 问题不在于输入数据...我可以在失败的输入上重新运​​行该程序,它将正确运行。

1
从经纬度计算UTM区域?
我正在尝试将latlong点转换为UTM。 要定义UTM投影,我需要计算该点的区域。 我正在尝试找出最佳方法。 一种方法是使用经度值找到合适的区域。 这将需要大量的编码。 我很好奇是否有更好的方法吗?

By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.