arcgis-desktoparcpy
The example below shows three raster (i.e. in black) and a single county polygon feature class. How can I clip/extract by mask the rasters by their corresponding county polygons (i.e. Edwards, Clark, Barber) in an automated fashion? For this example, the output should include three separate rasters clipped to the county boundaries. A python or ArcGIS solution suitable for batch processing would be ideal. Thanks for the help.
After using Spatial Join or Intersect to get the County name attribute onto each point, try using Dan Patterson's Split Layer by Attributes tool available on the Geoprocessing Model and Script Tool Gallery.
Alternatively you could use ModelBuilder to automate this using a different approach involving Select Layer by Location (click thumbnail for full image).
Lifted from this thread on the ESRI forums: Batch Selecting by Location and Exporting Shapefiles
One thing I notice is that your third parameter is a hard coded output (C:/data/lidar). The way it's written now will loop through each of your features and overwrite the output every time, but since you may not have allowed the automatic overwriting of files, this could potentially be the hang-up. Try creating a unique output name for each iteration:
#creating a workspace and a list of feature classes
arcpy.env.workspace = "C:/data/lidar/lidar.gdb"
fcs = arcpy.ListFeatureClasses()
#looping through each feature class and creating a raster based on the extent of
#feature class
i=0
for fc in fcs:
outputPath = "C:/data/lidar" + str(i)
i+=1
arcpy.Clip_management("perth", "#", outputPath, fc, "", "ClippingGeometry")
Also, I'm not certain that you intended to place the outputs in the C:/data folder named lidar... note that the third parameter in clip is the full path of your output raster, not a folder that it will be placed in. If you don't specify an extension in your output path name and place them in a standard folder, it will be a grid, so right now your program is attempting to create a new grid dataset named 'lidar' in the C:/data folder.
Best Answer
Presumably since you want to automate this, you'll want each raster to figure out which polygon should be used to clip itself? Pseudo code: