2009-10-29 39 views
1

我有一个相当简单的调用几个函数的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

“它不会工作”是什么意思? –

+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.shpblahblah\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" 

更新

与您的错误消息的前几个字谷歌搜索(总是一个好主意!)带来了这样的:troubleshooting geoprocessing errors它提供如下信息:所发生时 读取或写入的ArcSDE/DBMS数据 接收一个通用的“全捕获”错误 消息

地理处理的错误,如错误00210时 写输出

这继续建议一些方法来确定您确切的问题是什么。如果这对您没有帮助,您可以尝试在相关的ESRI论坛或GIS StackExchange上询问。

+0

我注意到的一个问题是,即使指定了不同的源,ArcMap有时也会更改脚本位置。我没有意识到这一点,我想如果你打开并关闭程序,验证它使用你指定的脚本,它应该没问题。 – womble

0

我会忍不住在

路径再看看= “C:\ LEED \ Cities_20_Oct \服务\”

当然,你想要双前斜线,不是双反斜线?

1

我看到这是一个3岁的帖子,但别人会增加:

arcpy.env.overwriteOutput=True # This allows the script to overwrite files. 

另外:

正如我产生python脚本弧形工作,我总是有我的导入之后你提到你删除了你的“文件夹”?这将是你的目录的一部分,我不知道你在脚本中创建目录的位置。你会想清除文件夹,而不是删除它(也许你的意思是你删除了文件)。

JJH

相关问题