topic Clipping a GRID file (input raster) to a Shapefile (Clip Extent) in Python Questions

Clipping a GRID file (input raster) to a Shapefile (Clip Extent)

AdelaideBaker — Mon, 09 May 2011 21:37:58 GMT

Hi,

I am using ArcGIS 10.0 and am trying to clip GRID files to a shapefile boundary. I can do this easily in ArcMap using Data Management Tools > Raster > Raster Processing > Clip. However when it comes to python I get an error about my inFolder. Because GRID files don't have a .shp, .tif etc. at the end do I have to add something else to my script? Thanks!

import arcpy

inFolder = "C:/Documents and Settings/someone/Desktop/Data/tmin_01"
resultsFolder = "C:/Documents and Settings/someone/Desktop/Data/output"
clipFeature = "C:/Documents and Settings/ajbaker/Desktop/PRISM_Trial/Boundary/Hydrographic_GB_boundary.shp"

arcpy.Clip_management("inFolder", "#", "clipFeature","resultsFolder", "0", "ClippingGeometry")

Re: Clipping a GRID file (input raster) to a Shapefile (Clip Extent)

DarrenWiens2 — Mon, 09 May 2011 21:51:02 GMT

inFolder is a variable that contains a string, so you don't put it in quotes (the way it is now, it's literally looking for a file with the path "inFolder"). Same with the other variables you pass as parameters. If you were directly entering the string in the clip command, you would enclose it in quotes.

Re: Clipping a GRID file (input raster) to a Shapefile (Clip Extent)

AdelaideBaker — Tue, 10 May 2011 14:45:27 GMT

Thanks. Problem solved!

Clip Multiple GRID layers

AdelaideBaker — Thu, 12 May 2011 19:47:31 GMT

I am very new at all of this and have another problem. How do I clip multiple GRID files. Also how might I be able to overwrite the input file? Here is my script so far:

import arcpy

arcpy.env.overwriteOutput = True

inFolder = "C:/Documents and Settings/user/Desktop/Data/tmin_1900_01"
resultsFolder = "C:/Documents and Settings/user/Desktop/Data/tmin_1900_01"
clipFeature = "S:/Projects/Wetlands/Title/Maps/Data layers/GB_boundary.shp"
arcpy.workspace = resultsFolder

#Process: Clip...
for input_feature in inFolder:
for clip_feature in clipFeature:
outfile = arcpy.CreateUniqueName(resultsFolder, inFolder)

#Usage: Clip_management in_raster rectangle out_raster {in_template_dataset} {nodata_value} {NONE | ClippingGeometry}
arcpy.Clip_management(inFolder, "#", outfile, clipFeature, "0", "ClippingGeometry")

print "Clip_management worked!"

print arcpy.GetMessages()