Netcdf转raster出现的错误

tianlangxing

我想把netcdf文件转为raster，结果出现了如下错误：

Traceback (most recent call last):
  File "1.netcdf_to_raster.py", line 86, in <module>

  File "C:\Program Files (x86)\ArcGIS\Desktop10.5\ArcPy\arcpy\management.py", line 14213, in CopyRaster
    raise e
arcgisscripting.ExecuteError: ERROR 999999: Error executing function.
Failed to create raster dataset
Failed to execute (CopyRaster).

Google了好久也没找到答案，请论坛上高手帮忙看看错误在哪里。我的源代码如下：

1. #Make rasters *.tif file for all time steps (time) in a netcdf file  
   
2. 
   
3. ##
   
4. 
   
5. #Import modules
   
6. 
   
7. import arcpy, os, time, datetime, calendar
   
8. 
   
9. import arcpy
   
10. 
    
11. from arcpy import env  
    
12. 
    
13. from arcpy.sa import *  
    
14. 
    
15.   
    
16. #########Path to folder that contains .nc data
    
17. 
    
18. #########This may need to be changed depending on how data are stored
    
19. 
    
20. arcpy.env.workspace = "D:/Example/GSMD" # Change folder++  
    
21. 
    
22. #########  
    
23. 
    
24.   
    
25. arcpy.env.overwriteOutput = True  
    
26. 
    
27. arcpy.env.scratchWorkspace = "D:/Example/GSMD"  
    
28. 
    
29.   
    
30. 
    
31. # Determines the number of netcdf files in the workspace
    
32. 
    
33. 
    
34. #Data processing starts
    
35. 
    
36. filename = "pr_2017.nc"
    
37. 
    
38. 
    
39. 
    
40. print("Processing: " + filename)  
    
41. 
    
42.     # Set local variable (grabs netcdf file(i))
    
43. 
    
44. inNetCDF = arcpy.env.workspace + "/" + filename  
    
45. 
    
46.   
    
47. 
    
48. ######### This changes depending on variable of interest (tmmx, tmmn, pr, pdsi)  
    
49. 
    
50. variable =  "precipitation_amount"
    
51. 
    
52.     #"surface_downwelling_shortwave_flux_in_air"
    
53. 
    
54.     #"specific_humidity"# #"precipitation_amount"# "relative_humidity" ###
    
55. 
    
56.     #"palmer_drought_severity_index"##"air_temperature"   
    
57. 
    
58.   
    
59. 
    
60. #########Change the variable ###  
    
61. 
    
62. InMemory_netcdf_raster = "precipitation_amount"
    
63. 
    
64.     #"surface_downwelling_shortwave_flux_in_air"#"relative_humidity"
    
65. 
    
66.     #"palmer_drought_severity_index"#"air_temperature"   
    
67. 
    
68. #########  
    
69. 
    
70.       
    
71. 
    
72. #########Where the output tiffs will be stored   
    
73. 
    
74. outRasterFolder = "D:/Example/GSMD/PR_rasters"  
    
75. 
    
76. #########
    
77. 
    
78. #######change the lengths depending on variable (pr is shorter than tmmx) -
    
79. 
    
80. ########This variable assists with naming convention of out file  
    
81. 
    
82. fileroot = filename[3:(len(filename)-3)] #change the name lengths too  
    
83. 
    
84. #######
    
85. 
    
86.       
    
87. 
    
88. Yr = int(fileroot)   
    
89. 
    
90. print "Processing --> " + str(Yr)  
    
91. 
    
92.       
    
93. 
    
94. #Month info  
    
95. 
    
96. allmnths = range(1,13)  
    
97. 
    
98. for mnths in allmnths:  
    
99. 
    
100.     Lastday = calendar.monthrange(Yr, mnths)[1]  
     
101. 
     
102.     MRange = range(1,Lastday+1)  
     
103. 
     
104.     for dyy in MRange:  
     
105. 
     
106.         dyys =int(dyy)  
     
107. 
     
108.         mnthsss = int(mnths)  
     
109. 
     
110.               
     
111. 
     
112.         a = str(mnthsss)+ "/" + str(dyys) +"/" + str(Yr)  
     
113. 
     
114.   
     
115. 
     
116.         #IF Days or Months less than 10, we want to put zero at the front so it is easier to sort.  
     
117. 
     
118.   
     
119. 
     
120.         if (mnthsss < 10) and (dyys < 10):   
     
121. 
     
122.             b = str(Yr)+ "_" + "0" + str(mnthsss)+ "_" + "0" + str(dyys)  
     
123. 
     
124.         elif (mnthsss < 10) and (dyys >=10):  
     
125. 
     
126.             b = str(Yr)+ "_" + "0" + str(mnthsss)+ "_" + str(dyys)  
     
127. 
     
128.   
     
129. 
     
130.         elif (mnthsss >= 10)and (dyys < 10):  
     
131. 
     
132.             b = str(Yr)+ "_" + str(mnthsss)+ "_" + "0" + str(dyys)                     
     
133. 
     
134.         else:  
     
135. 
     
136.             b = str(Yr)+ "_" + str(mnthsss)+ "_" + str(dyys)  
     
137. 
     
138.                   
     
139. 
     
140.         dimensionValues = "day " + a
     
141. 
     
142.               
     
143. 
     
144. ##########CHANGE the word too##              
     
145. 
     
146.                                         #################  
     
147. 
     
148.         outTIFF = outRasterFolder+"/"+"pr_"+b+".tif"  
     
149. 
     
150.                                         #################  
     
151. 
     
152.   
     
153. 
     
154.               
     
155. 
     
156.             # Execute MakeNetCDFRasterLayer         ##############
     
157. 
     
158.         arcpy.MakeNetCDFRasterLayer_md(inNetCDF,"precipitation_amount" ,"lon","lat",InMemory_netcdf_raster,"",dimensionValues,"BY_VALUE")  
     
159. 
     
160.             #print "Created NetCDF Layer for " + a   ############
     
161. 
     
162.               
     
163. 
     
164.             # Process: Copy Raster  
     
165. 
     
166.         arcpy.CopyRaster_management(InMemory_netcdf_raster,outTIFF,"","","","NONE","NONE")  
     
167. 
     
168.         print "Created Raster for pr_" + b  
     
169. 
     
170.                                 #################  
     
171. 
     
172.         b=""
     
173. 
     
174.         dimensionValues=""
     
175. 
     
176.             #InMemory_netcdf_raster=""
     
177. 
     
178.             #inNetCDF=""
     
179. 
     
180. print "\nEnd of processing!!!"  
     
181. 
     

复制代码

在D盘建了一个文件夹，名为"Example"。里面包含一个子文件夹GSMD和上述Python程序文件。GSMD中包含输入文件pr_2017.nc和子文件夹PR_rasters.
如果需要测试程序的话，pr_2017.nc可以在这里下载：http://northwestknowledge.net/metdata/data/pr_2017.nc. 大概有500多兆.

先谢谢大家。

tianlangxing

有朋友说在他的电脑上测试没问题，现在我也不知道问题在哪里。我重新安装过ArcGIS，还是不管用。

tianlangxing

还没有解决，继续请教