topic Re: getting started with Python in Python Questions

getting started with Python

WadeGivens — Tue, 09 Aug 2011 03:16:11 GMT

Well, finally breaking down and trying to get going with Python. My first attempt is to create a batch raster clip tool that loops through a set of rasters, clips them, and saves then to another directory with "_clip" appended to the file name. The script below runs with no errors. Problem is no results either. Any ideas, or thoughts to get me going in the right direction?

import arcpy
import os

inDIR = "C:/Temp/imagery"
outDIR = "C:/Temp/imagery/clipped"
clipFC = "C:/Temp/imagery/clip_feature.shp"

arcpy.env.workspace = inDIR

rasters = arcpy.ListRasters("*", "TIFF")

for raster in rasters:
 print raster
 outname = raster[:-4] + "_clip.tif"
 arcpy.Compression= "LZW"
 arcpy.Clip_Management(raster, "#", os.path.join(outDIR, raster), clipFC, "#", "ClippingGeometry")

Thanks!
Wade

Re: getting started with Python

JakeSkinner — Sun, 12 Dec 2021 07:13:03 GMT

I believe the problem is with the 'arcpy.ListRasters' function. You will want to specify "TIF" instead of "TIFF". Also, I believe you will want to specify 'os.path.join(OutDir, outname)' within the 'arcpy.Clip_management' function. Ex:

rasters = arcpy.ListRasters("*", "TIF")

for raster in rasters:
 print raster
 outname = raster[:-4] + "_clip.tif"
 arcpy.Compression= "LZW"
 arcpy.Clip_management(raster, "", os.path.join(outDIR, outname), clipFC, "", "ClippingGeometry")

Re: getting started with Python

WadeGivens — Tue, 09 Aug 2011 12:35:00 GMT

Jake:

Thanks! That did the trick. Been kinda hard to get going with some of this. Lots of posts out there that give examples using gp. calls. Not so many using arcpy. Trying to get going with arcpy since that is preferred syntax now.

Thanks again!

Wade

Re: getting started with Python

WadeGivens — Sun, 12 Dec 2021 07:13:05 GMT

Jake:

Now that I have this going, I've added the script to a toolbox in ArcMap and I'm trying to pass the workspace directories and clip shape as parameters (see below). I get the following error when trying to run the tool from the toolbox:

Executing: Clip C:\Temp\imagery C:\Temp\imagery\clipped\clipped C:\Temp\imagery\all_corn_IMPACT_2011.shp
Start Time: Tue Aug 09 09:50:33 2011
Running script Clip...
ERROR 000714: Error in script Clip.
Error in executing: cmd.exe /C C:\GIS_PR~1\PYTHON~1\BATCH_~1 "C:\Temp\imagery" "C:\Temp\imagery\clipped\clipped" "C:\Temp\imagery\all_corn_IMPACT_2011.shp"

Failed to execute (Clip).
Failed at Tue Aug 09 09:50:33 2011 (Elapsed Time: 0.00 seconds)

import arcpy
import os

inDIR = arcpy.GetParameterAsText(0)
outDIR = arcpy.GetParameterAsText(1)
clipFC = arcpy.GetParameterAsText(2)

# inDIR = "C:/Temp/imagery"
# outDIR = "C:/Temp/imagery/clipped"
# clipFC = "C:/Temp/imagery/all_corn_IMPACT_2011.shp"

arcpy.env.workspace = inDIR

rasters = arcpy.ListRasters("*", "TIF")

for raster in rasters:
 print raster
 outname = raster[:-4] + "_clip.tif"
 arcpy.Compression= "LZW"
 arcpy.Clip_management(raster, "", os.path.join(outDIR, outname), clipFC, "", "ClippingGeometry")

Any ideas on what may be causing this error?

Thanks,
Wade

Re: getting started with Python

JakeSkinner — Tue, 09 Aug 2011 14:14:56 GMT

Make sure you have the correct data types setup for the parameters.

Re: getting started with Python

WadeGivens — Sun, 12 Dec 2021 07:13:08 GMT

Still getting the same error:

Executing: Script C:\GIS_Projects\2011\CE_BU_RS\RapidEye\Aug_1st_week\images\5_band C:\GIS_Projects\2011\CE_BU_RS\RapidEye\Aug_1st_week\images\5_band\clipped C:\GIS_Projects\2011\CE_BU_RS\RapidEye\shapefile_bnds\all_corn_IMPACT_2011.shp
Start Time: Tue Aug 09 10:30:53 2011
Running script Script...
ERROR 000714: Error in script Script.
Error in executing: cmd.exe /C C:\GIS_PR~1\PYTHON~1\BATCH_~1 "C:\GIS_Projects\2011\CE_BU_RS\RapidEye\Aug_1st_week\images\5_band" "C:\GIS_Projects\2011\CE_BU_RS\RapidEye\Aug_1st_week\images\5_band\clipped" "C:\GIS_Projects\2011\CE_BU_RS\RapidEye\shapefile_bnds\all_corn_IMPACT_2011.shp"

Failed to execute (Script).
Failed at Tue Aug 09 10:30:53 2011 (Elapsed Time: 0.00 seconds)

Screenshot of my script parameters attached. I can copy and paste the paths into the script and it works fine


import arcpy
import os

inDIR = arcpy.GetParameterAsText(0)
outDIR = arcpy.GetParameterAsText(1)
clipFC = arcpy.GetParameterAsText(2)

arcpy.env.workspace = inDIR

rasters = arcpy.ListRasters("*", "TIF")

for raster in rasters:
 print raster
 outname = raster[:-4] + "_clip.tif"
 arcpy.Compression= "LZW"
 arcpy.Clip_management(raster, "", os.path.join(outDIR, outname), clipFC, "", "ClippingGeometry")

Re: getting started with Python

WadeGivens — Tue, 09 Aug 2011 14:36:53 GMT

Well I feel stupid. My script wasn't saved with a .py extension. That solved the issue

Wade