Raster ExtractByMask with python

[bogdanpalade1]

Hello,

I am trying to clip a raster with a polygon shapefile.The shapefile contains several polygons representing administrative regions. The method is to temporarily export each polygon to a new shapefile and to extract the coresponding piece of raster.

The problem is that the output of

gp.ExtractByMask_sa(raster, shapefile, newRaster) 

is always just one band of the original RGB raster. Do you know how to extract all the 3 bands?

The original raster is .img

Thank you

[ChrisSnyder]

Try naming your shapefile something else than a number (1.shp). Generally it is a bad idea to name a layer or field with a number as the 1st character.

[bogdanpalade1]

Tried changing the shapefile name but still the same error....

[ChrisSnyder]

Not sure what's wrong... Take look at the sample code: http://forums.arcgis.com/threads/9512-Rectangle-for-Clip_management-tool-in-arcpy (my post Aug 2nd @10:25). 

Are you sure that your raster and shapefile are using the same projection?

[JamieWilliams]

I had the same problem. I solved it by adding the full file path to the feature name.
eg. pathFeature = "R:\\folderContainingShpFile\\" + str(feature)

gp.clip_management(raster,"446716,298473 4425996,437405 499153,348855 4492682,68626",outRaster,pathFeature,"255","ClippingGeometry")