arcpy.env.overwriteOutput = True

最新推荐文章于 2023-11-08 15:41:30 发布
最新推荐文章于 2023-11-08 15:41:30 发布
阅读量2.9k 收藏 8
点赞数 1
CC 4.0 BY-SA版权
分类专栏： Python for arcgis
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.youkuaiyun.com/qq_36196621/article/details/90237104
探讨如何使用 ArcPy 的环境设置中的 overwriteOutput 参数来覆盖已存在的输出文件，简化地理信息系统 (GIS) 数据处理流程。

摘要生成于 C知道 ，由 DeepSeek-R1 满血版支持， 前往体验 >


arcpy.env.overwriteOutput = True：输出文件夹里面已经有内容的，就覆盖掉。

 

浅谈使用ArcPy执行大数据量处理任务
慕名ArcGIS
04-28 1万+
Python功能强大而易于学习。对于ArcGIS for Desktop用户来讲，Python是提高工作效率的不二选择。 Arcpy是esri提供的用于高效数据处理分析、制图等的Python站点包。 利用ArcPy，我们可以在ArcMap的Python窗口交互执行脚本，还可以创建自定义脚本工具或脚本工具箱，还可以在ArcGIS之外运行独立脚本，享受更纯正的python体验。这一篇说说如何利用Pyth
Arcpy操作
遇事不决，问春风
10-04 2261
Arcgis-Arcpy 持续更新 安装 安装步骤 基础介绍 地理信息系统（geography information system，GIS）是在计算机硬件支持下，对整个或者部分地球表层空间中的有关地理分布数据进行采集，存储，管路，运算，分析，显示和描述的技术系统。地理信息系统处理和管理的对象是多种地理空间实体数据及其关系，包括空间定位数据，图形数据，遥感图像数据，属性数据，主要用于分析和处理一定地理区域内分布的各种现象和过程，解决复杂的规划，决策的管理问题。（有些官方，看看就行） 其他：略（嘿嘿嘿，可以去
arcgis arcpy发布地图服务，并且覆盖已有的服务
huangguo1688的博客
12-07 660
# -*- coding: utf-8 -*- #以防中文报错 import sys reload(sys) sys.setdefaultencoding('utf-8') import arcpy,os import arcpy.mapping as mapping import xml.dom.minidom as DOM arcpy.env.overwriteOutput = True #发布地图服务 folder = r"E:\python\数据" mxdpath=r"E:\python\数..
[ArcPy] 代码汇总
geodoer
09-22 6394
初识ArcPy 使用ArcPy的方法 arcpy批处理示例 地理处理工具 简单示例 工具输入框内自动选择图层和数据 读取XML文件 要素属性操作 复制数据库 处理Excel数据生成多边形
ArcPy常用基础功能
漓涂
06-14 3771
先来一段简单的ArcPy脚本： # -*-coding:utf-8-*- import arcpy try: source_path_name = r"F:/GIS测试数据/测试.gdb" source_name = "DLTB" target_name = "Target_copy" arcpy.env.workspace = source_path_...
ERA5数据处理，简单记录。Arcpy
新手，不断自学
11-17 6202
第一次写这个，还有点小激动，新手，简单记录一下，地学同一问题方法上可能有不同，代码和方法如果有改进的，还请多多指出！ hhh 其实代码还在跑着，不过目前没问题，先写上，写的是ERA5数据下载完经过mosaic后的处理过程，下次有时间再写下载和拼接的过程吧。 以提取中国2015年小时级的lai_lv（Leaf area index, low vegetation）变量为例（ERA5的数据提取的时候是按变量的简写来提取的，详见ERA5-Land: data documentation - Cop...
import arcpy import os def batch_kml_to_shp(input_folder, output_folder, keep_intermediate=False): """ 批量将 KML 文件转换为 SHP 格式 :param input_folder: 包含 KML 文件的输入文件夹 :param output_folder: 输出 SHP 文件的目录 :param keep_intermediate: 是否保留中间文件 (默认 False) """ try: # 设置工作环境 arcpy.env.overwriteOutput = True # 创建临时地理数据库存放中间结果 temp_gdb = arcpy.CreateFileGDB_management( out_folder_path=output_folder, out_name="TempKMLConversion.gdb" )[0] # 遍历所有 KML 文件 kml_files = [f for f in os.listdir(input_folder) if f.lower().endswith('.kml')] for kml_file in kml_files: try: kml_path = os.path.join(input_folder, kml_file) base_name = os.path.splitext(kml_file)[0] # 步骤 1: 将 KML 转换为图层文件 (临时地理数据库) arcpy.AddMessage(f"正在处理: {kml_file}") output_layer = os.path.join(temp_gdb, base_name) arcpy.KMLToLayer_conversion( kml_path, output_layer, "NO_GROUNDOVERLAY" ) # 步骤 2: 从地理数据库提取要素类 feature_classes = arcpy.ListFeatureClasses(f"{base_name}*", "ALL", output_layer) for fc in feature_classes: # 跳过注记和栅格图层 if "Placemarks" not in fc: continue # 步骤 3: 转换为 Shapefile shp_name = f"{base_name}_{fc.split('_')[-1]}.shp" shp_path = os.path.join(output_folder, shp_name) arcpy.FeatureClassToFeatureClass_conversion( os.path.join(output_layer, fc), output_folder, shp_name ) arcpy.AddMessage(f"已创建: {shp_name}") except Exception as e: arcpy.AddWarning(f"处理 {kml_file} 时出错: {str(e)}") # 清理中间文件 if not keep_intermediate: arcpy.Delete_management(temp_gdb) arcpy.AddMessage("已清理临时文件") except arcpy.ExecuteError: arcpy.AddError(arcpy.GetMessages(2)) except Exception as e: arcpy.AddError(f"脚本执行错误: {str(e)}") if __name__ == "__main__": # 脚本工具参数 input_folder = arcpy.GetParameterAsText(0) output_folder = arcpy.GetParameterAsText(1) keep_intermediate = arcpy.GetParameter(2) # 布尔型参数 batch_kml_to_shp(input_folder, output_folder, keep_intermediate)修复运行成果没有成果
06-15
folder:输出SHP文件的目录:paramkeep_intermediate:是否保留中间文件(默认False)"""try:#设置工作环境arcpy.env.overwriteOutput=True#验证输入输出路径ifnotarcpy.Exists(input_folder):arcpy.AddError("输入文件夹...
import arcpy import os def batch_kml_to_shp(input_folder, output_folder, keep_intermediate=False): # 环境检查 arcpy.AddMessage("=== 环境验证 ===") arcpy.AddMessage(f"ArcGIS Pro版本: {arcpy.GetInstallInfo()['Version']}") # 路径验证 if not arcpy.Exists(input_folder): arcpy.AddError("错误: 输入文件夹不存在") return if not os.access(output_folder, os.W_OK): arcpy.AddError("错误: 输出文件夹不可写") return # 许可检查 if arcpy.CheckExtension("kml") != "Available": arcpy.AddError("错误: 缺少KML转换扩展许可") return else: arcpy.CheckOutExtension("kml") try: arcpy.env.overwriteOutput = True arcpy.AddMessage("环境设置完成") # 创建临时GDB temp_gdb_path = os.path.join(output_folder, "TempKMLConversion.gdb") if arcpy.Exists(temp_gdb_path): arcpy.Delete_management(temp_gdb_path) temp_gdb = arcpy.CreateFileGDB_management(output_folder, "TempKMLConversion.gdb")[0] arcpy.AddMessage(f"临时GDB创建于: {temp_gdb}") # 获取KML文件列表 kml_files = [f for f in os.listdir(input_folder) if f.lower().endswith('.kml')] if not kml_files: arcpy.AddError("错误: 未找到KML文件") return arcpy.AddMessage(f"找到{len(kml_files)}个KML文件") for kml_file in kml_files: try: kml_path = os.path.join(input_folder, kml_file) base_name = os.path.splitext(kml_file)[0] # KML转换 output_layer = os.path.join(temp_gdb, base_name) arcpy.AddMessage(f"正在转换: {kml_file} → {output_layer}") arcpy.KMLToLayer_conversion(kml_path, output_layer, "NO_GROUNDOVERLAY") # 定位Placemarks要素类 arcpy.AddMessage("定位转换结果...") feature_classes = arcpy.ListFeatureClasses("*", "ALL", output_layer) placemarks = [fc for fc in feature_classes if "Placemarks" in fc] if not placemarks: arcpy.AddWarning(f"警告: {output_layer}中未找到Placemarks要素类") continue # 输出SHP for fc in placemarks: shp_name = f"{base_name}_{fc.split('_')[-1]}.shp" shp_path = os.path.join(output_folder, shp_name) arcpy.AddMessage(f"生成: {shp_name}") arcpy.FeatureClassToFeatureClass_conversion( os.path.join(output_layer, fc), output_folder, shp_name ) except arcpy.ExecuteError: arcpy.AddWarning(f"处理失败: {kml_file}\n{arcpy.GetMessages(2)}") except Exception as e: arcpy.AddWarning(f"未知错误: {kml_file}\n{str(e)}") finally: # 清理和许可归还 if not keep_intermediate and arcpy.Exists(temp_gdb_path): arcpy.Delete_management(temp_gdb_path) arcpy.AddMessage("临时文件已清理") arcpy.CheckInExtension("kml") arcpy.AddMessage("处理完成") if __name__ == "__main__": input_folder = arcpy.GetParameterAsText(0) output_folder = arcpy.GetParameterAsText(1) keep_intermediate = arcpy.GetParameter(2) batch_kml_to_shp(input_folder, output_folder, keep_intermediate)运行显示无kml许可但是用软件自带的kml转shp可以运行成功
06-15
arcpy.env.overwriteOutput = True try: # ===== 3. 创建临时地理数据库 ===== temp_gdb = os.path.join(output_folder, "TempKML.gdb") if arcpy.Exists(temp_gdb): arcpy.management.Delete(temp_gdb) ...
# -*- coding: utf-8 -*- import arcpy import numpy as np from arcpy.sa import * # 启用Spatial Analyst许可 arcpy.CheckOutExtension("Spatial") # 设置环境 arcpy.env.workspace = r"D:\YourFolder" # 替换为你的实际路径 arcpy.env.overwriteOutput = True # 输入数据（检查是否存在） if not arcpy.Exists("light_pro.tif") or not arcpy.Exists("veg_stack.tif"): arcpy.AddError("输入栅格不存在！请检查文件名") else: light_raster = Raster("light_pro.tif") veg_raster = Raster("veg_stack.tif") # ===== 1. 计算相关系数栅格 ===== # 方法：使用局部统计（3x3窗口） try: arcpy.AddMessage("正在计算相关系数栅格...") out_corr = Correlate(light_raster, veg_raster, "Rectangle 3 3 CELL") # 修改为矩形窗口 out_corr.save("Correlation_1_2.tif") # ===== 2. 分类协同/权衡关系 ===== arcpy.AddMessage("正在进行分类...") corr_raster = Raster("Correlation_1_2.tif") out_class = Con( Abs(corr_raster) > 0.5, Con(corr_raster > 0, 1, 2), # 1=协同, 2=权衡 3 # 3=无显著相关 ) out_class.save("Final_Classification.tif") # ===== 3. 自动符号化 ===== mxd = arcpy.mapping.MapDocument("CURRENT") df = arcpy.mapping.ListDataFrames(mxd)[0] # 添加结果图层 result_layer = arcpy.mapping.Layer("Final_Classification.tif") arcpy.mapping.AddLayer(df, result_layer, "TOP") # 设置分类颜色 sym = result_layer.symbology sym.renderer = "Classified" sym.classificationField = "Value" sym.classBreakValues = [1, 2, 3] sym.classBreakLabels = ["协同", "权衡", "无显著相关"] # 使用预设色带（ArcMap 10.8自带） color_ramp = arcpy.mapping.ListColorRamps(mxd, "Yellow-Green-Light to Dark")[0] # 修改为存在的色带名 sym.colorRamp = color_ramp result_layer.symbology = sym # 刷新视图 arcpy.RefreshActiveView() arcpy.RefreshTOC() mxd.save() arcpy.AddMessage("处理成功！结果已保存为 Final_Classification.tif") except Exception as e: arcpy.AddError(f"处理失败: {str(e)}") # 释放许可 arcpy.CheckInExtension("Spatial")第65行有问题
06-28
arcpy.env.overwriteOutput = True # 允许覆盖输出 ``` #### ⚠️ 注意事项 1. 拓扑验证前需确保： ```python # 检查拓扑是否存在 if arcpy.Exists(topology_path): # 执行验证 else: print(f"错误：拓扑 {...
# -*- coding: utf-8 -*- import arcpy import sys reload(sys) # Python 2.7 需要重载sys sys.setdefaultencoding("utf-8") # 解决中文路径问题[^3][^4] def intersect_layers(input_layers, output_feature_class): """ 执行多个图层的相交操作（适配ArcGIS 10.8） :param input_layers: 输入图层路径列表 :param output_feature_class: 输出要素类完整路径 """ # 1. 坐标系检查（关键步骤） first_ref = arcpy.Describe(input_layers[0]).spatialReference for i, layer in enumerate(input_layers[1:]): layer_ref = arcpy.Describe(layer).spatialReference if first_ref.name != layer_ref.name: arcpy.AddWarning("坐标系不匹配: {} ({}) vs {} ({})".format( input_layers[0], first_ref.name, layer, layer_ref.name )) # 自动投影到第一个图层的坐标系 projected = "in_memory/projected_{}".format(i) arcpy.Project_management(layer, projected, first_ref) input_layers[i+1] = projected # 替换为投影后的临时图层 # 2. 构建相交输入格式 in_features = [[layer, ""] for layer in input_layers] # 3. 执行相交分析 arcpy.Intersect_analysis( in_features=in_features, out_feature_class=output_feature_class, join_attributes="ALL", cluster_tolerance="", output_type="INPUT" ) return output_feature_class def add_layer_to_map(output_path): """将结果图层添加到当前地图""" aprx = arcpy.mp.ArcGISProject("CURRENT") active_map = aprx.activeMap active_map.addDataFromPath(output_path) arcpy.AddMessage(u"结果已添加到当前地图: {}".format(output_path)) if __name__ == "__main__": try: # ========== 环境初始化 ========== arcpy.env.overwriteOutput = True arcpy.AddMessage(u"工具初始化完成") # ========== 许可检查 ========== if arcpy.CheckExtension("3D") == "Available": arcpy.CheckOutExtension("3D") arcpy.AddMessage(u"3D Analyst 许可已启用") else: raise Exception(u"3D Analyst 许可不可用") # ========== 参数获取 ========== input_layers_str = arcpy.GetParameterAsText(0) # 多值输入参数 output_feature_class = arcpy.GetParameterAsText(1) # 输出路径 # 拆分多值参数（分号分隔） input_layers = [lyr.strip() for lyr in input_layers_str.split(";")] arcpy.AddMessage(u"输入图层数: {}".format(len(input_layers))) # ========== 输入验证 ========== if len(input_layers) < 2: raise ValueError(u"需要至少两个输入图层") for layer in input_layers: if not arcpy.Exists(layer): raise IOError(u"图层不存在: {}".format(layer)) # ========== 工作空间设置 ========== # 从输出路径提取工作空间[^3] workspace = os.path.dirname(output_feature_class) if not workspace: workspace = arcpy.env.workspace or "in_memory" arcpy.env.workspace = workspace.replace("\\", "/") # 路径标准化[^3] arcpy.AddMessage(u"工作空间: {}".format(arcpy.env.workspace)) # ========== 执行相交操作 ========== result = intersect_layers(input_layers, output_feature_class) # ========== 结果验证 ========== count = int(arcpy.GetCount_management(result).getOutput(0)) arcpy.AddMessage(u"输出要素数量: {}".format(count)) # ========== 添加到地图 ========== add_layer_to_map(result) arcpy.SetParameterAsText(2, result) # 设置输出参数 arcpy.AddMessage(u"相交操作成功完成！") except arcpy.ExecuteError: err = arcpy.GetMessages(2) arcpy.AddError(u"地理处理错误: {}".format(err.decode('utf-8'))) except Exception as e: arcpy.AddError(u"错误: {}".format(str(e).decode('utf-8'))) # 获取错误行号 exc_type, exc_obj, exc_tb = sys.exc_info() arcpy.AddError(u"错误行号: {}".format(exc_tb.tb_lineno)) finally: # ========== 清理资源 ========== arcpy.CheckInExtension("3D") arcpy.AddMessage(u"3D Analyst 许可已释放") arcpy.Delete_management("in_memory") # 清理临时数据 这个代码里的参数怎么设置
最新发布
07-09
arcpy.env.overwriteOutput = True # 允许覆盖现有数据 arcpy.env.outputCoordinateSystem = input_layers[0].spatialReference ``` #### 三、完整函数调用示例 ```python # 定义输入输出 input_data = ["roads...
import arcpy import os def intersect_layers(input_layers, output_feature_class): """ 执行多个图层的相交操作 :param input_layers: 输入图层列表 :param output_feature_class: 输出要素类路径 """ # 构建相交工具所需的输入格式 in_features = [[layer, ""] for layer in input_layers] # 执行相交分析 arcpy.Intersect_analysis( in_features=in_features, out_feature_class=output_feature_class, join_attributes="ALL", cluster_tolerance="", output_type="INPUT" ) return output_feature_class def add_layer_to_map(output_path): """将结果图层添加到当前地图""" aprx = arcpy.mp.ArcGISProject("CURRENT") active_map = aprx.activeMap active_map.addDataFromPath(output_path) arcpy.AddMessage("结果已添加到当前地图: {}".format(output_path)) if __name__ == "__main__": # 设置环境 arcpy.env.overwriteOutput = True try: # 获取输入参数 input_layers_str = arcpy.GetParameterAsText(0) # 多个输入图层（分号分隔） output_feature_class = arcpy.GetParameterAsText(1) # 输出要素类路径 # 拆分多值参数为图层列表（修正缩进） input_layers = input_layers_str.split(";") if input_layers_str else [] # 验证至少有两个输入图层 if len(input_layers) < 2: arcpy.AddError("需要至少两个输入图层进行相交操作！") raise ValueError("需要至少两个输入图层") # 验证所有输入图层存在 for layer in input_layers: if not arcpy.Exists(layer): arcpy.AddError(f"输入图层不存在: {layer}") raise FileNotFoundError(f"图层 {layer} 不存在") # 调用相交函数 result = intersect_layers(input_layers, output_feature_class) # 将结果添加到地图 add_layer_to_map(result) arcpy.AddMessage("相交操作成功完成！") arcpy.AddMessage(f"输入图层数: {len(input_layers)}") arcpy.AddMessage(f"输出结果: {result}") except arcpy.ExecuteError: arcpy.AddError("工具执行错误: " + arcpy.GetMessages(2)) except Exception as e: arcpy.AddError("发生错误: " + str(e)) 在此基础帮我看看我的代码=都缺少什么，我是用的是10.8版本
07-09
- 建议开启覆盖输出：`arcpy.env.overwriteOutput = True`[^3][^4] #### 二、原代码（引用[4]）的改进点 ```python # 改进后的相交分析脚本（兼容 ArcGIS 10.8） # -*- coding: utf-8 -*- import arcpy import sys...
ArcGisPro脚本工具【3】——高德POI爬取
热门推荐
规划GIS会
12-29 1万+
本次Arcpy脚本工具用的是ArcGisPro2.8版本，尽量使用同一版本，其它版本可以参照。
【ArcPy】创建栅格金字塔
hanv79的博客
11-08 589
在加载数据量比较大的栅格影像时，通常构建栅格金字塔可以实现快速加载；下面就是ArcPy中构建金字塔的核心代码。
【Arcpy】实验一：创建缓冲区
小熊的博客
07-25 912
【Arcpy】Calling ArcGIS Pro geoprocessing tools from scripts
（十八）arcpy开发&利用arcpy实现在arcgis中对字段的删除
yGIS
02-20 2011
本节，我们将带领大家一起来学习一下关于字段删除方面的函数。我们将会删除字段中除了FID shape外的其他字段。其中利用到函数字段列举函数。arcpy.ListField,下面定义的这个函数将返回所有的字段。 def getFieldNames( fc): fieldList = [] for f in arcpy.ListFields(fc): ...
使用ArcPy在独立脚本中使用克里金插值
老靳的WebGIS
09-09 2270
使用ArcPy在独立脚本中使用克里金插值 import arcpy from arcpy.sa import * arcpy.env.overwriteOutput = True # 设置输入点要素 inFeatures = '你的输入要素' # 设置输出结果 outRaster = "你的输出栅格" # 设置输出字段 field = "DRP" # 设置输出栅格像元大小 cellSize = 0.0230704 # 设置半变异函数属性 kModelOrdinary = KrigingModelOrdin
arcpy批量化出图浅试
Hzz的地球调查组
09-08 1003
arcpy批量化出图浅试
利用arcpy按字段值筛选要素工具
github_35700281的博客
10-24 6751
ArcGIS按字段值筛选要素做缓冲区分析时要先进行字段筛选然后再进行缓冲区分析。而不能直接在缓冲区工具中直接筛选，于是乎就有了本工具的编写。 工具中一共有7个参数： 一.输入要素 （shapeFilel类型） 二.SQL筛选语句（SQL语句） 三.输出要素（shapeFile类型input） 四.缓冲半径（字符串） 五,侧类型（字符串） 六.末端类型（字符串） 七.融合类型（字符串
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值