我有一个相当简单的python循环,它调用一些函数,并将输出写入文件。为此,将创建一个文件夹,并将该文件保存在此文件夹中。
当我第一次使用唯一的文件名运行程序时,它运行正常。但是,如果我尝试再次运行它,它将无法工作,我不明白为什么。我很确定这不是覆盖文件的问题,因为我在重新运行之前删除了该文件夹,这是该文件存储的唯一位置。有一个我误解的概念吗?
有问题的文件是'buff1.shp'。我使用Python 2.5在ArcGIS中运行一些分析
感谢您的任何建议(包括有关如何改进我的编码风格的建议)。另外一个注意事项是我的循环目前只使用一个值,因为我现在正在测试它。
# Import system modules
import sys, string, os, arcgisscripting, shutil
# Create the Geoprocessor object
gp = arcgisscripting.create()
# Load required toolboxes...
gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Spatial Statistics Tools.tbx")
gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Analysis Tools.tbx")
# specify workspace
gp.Workspace = "C:/LEED/Cities_20_Oct/services"
path = "C:\\LEED\\Cities_20_Oct\\services\\"
results = 'results\\'
os.mkdir( path + results )
newpath = path + results
# Loop through each file (0 -> 20)
for j in range(0,1):
in_file = "ser" + str(j) + ".shp"
in_file_2 = "ser" + str(j) + "_c.shp"
print "Analyzing " + str(in_file) + " and " + str(in_file_2)
#Loop through a range of buffers - in this case, 1,2
for i in range(1,2):
print "Buffering....."
# Local variables...
center_services = in_file_2
buffer_shp = newpath + "buff" + str(i) + ".shp"
points = in_file_2
buffered_analysis_count_shp = newpath + "buffered_analysis_count.shp"
count_txt = newpath + "count.txt"
# Buffer size
b_size = 1000 + 1000 * i
b_size_input = str(b_size) + ' METERS'
print "Buffer:" + b_size_input + "\n"
# Process: Buffer...
gp.Buffer_analysis(center_services, buffer_shp, b_size_input, "FULL", "ROUND", "ALL", "")
print "over"
(为了澄清这个问题,我编辑了一些没有其余代码没有意义的部分。错误仍然存在于程序中。)
错误讯息:
ExecuteError: ERROR 000210: Cannot create output C:\LEED\Cities_20_Oct\services\results\buff1.shp Failed to execute (Buffer).
答案 0 :(得分:2)
我无法看到错误消息blahblah \ buff1.shp中的文件名如何来自您的代码。
for i in range(0,1):
buffer_shp = newpath + "buff" + str(i) + ".shp"
gp.Buffer_analysis(center_services, buffer_shp, etc etc)
应该生成blahblah\buff0.shp
而不是blahblah\buff1.shp
...我强烈建议您显示的代码应该是您实际运行的代码。在gp.Buffer_analysis()
调用之前抛出一个print语句来显示i和repr(buffer_shp)的值。显示所有打印结果。
注释#Loop through a range of buffers (1 ->100)
表示你想从1开始,而不是0。如果评论与代码匹配,它会对你有很大帮助。
不要重复自己;而不是
os.mkdir( path + results )
newpath = path + results
这样做:
newpath = path + results # using os.path.join() is even better
os.mkdir(newpath)
您可能希望养成使用os.path.join()构建所有路径的习惯。
您需要在循环外部调用os.mkdir()
,即每次运行脚本时执行一次,而不是每次循环内循环。
不使用这些陈述的结果:
buffered_analysis_count_shp = newpath + "buffered_analysis_count.shp"
count_txt = newpath + "count.txt"
<强>更新强>
使用您的错误消息中的前几个单词(总是一个好主意!)进行Google搜索:troubleshooting geoprocessing errors提供以下信息:
发生时的地理处理错误 读取或写入ArcSDE / DBMS数据 收到一个通用的'全能'错误 消息,如错误00210时 写输出
这继续提出了一些确定问题的方法。如果这对您没有帮助,您可以尝试在相关的ESRI论坛或GIS StackExchange上询问。
答案 1 :(得分:1)
我看到这是一个3岁的帖子,但对于其他人会添加:
当我生成python脚本以使用Arc时,我总是在导入后立即包含:
arcpy.env.overwriteOutput=True # This allows the script to overwrite files.
你还提到你删除了你的“文件夹”吗?这将是您的目录的一部分,我不知道您在脚本中创建目录的位置。你想要清除文件夹,而不是删除它(也许你的意思是删除文件)。
JJH
答案 2 :(得分:0)
我很想再看看
path =“C:\ LEED \ Cities_20_Oct \ services \”
当然你想要双前斜线,而不是双反斜杠?