Fixed it. This is what satisfied the colleague, and was put into production. Lines 30-34 is how I check to see if the current worksheet is a Grid or a workspace containing grids. Had to go through this extra BS because I couldn't use da.walk.
#testclipgridstifsshapes.py
#test trying to walk a directory tree, find grids, and clip to a area of Interest (AOI)
# and make it run in 10.0 so without da.walk.
# Paul Huffman, Yakama Fisheries, 7/21/2014
# Import arcpy module
import arcpy
import os
# Set Geoprocessing environments
#arcpy.env.scratchWorkspace = "c:\\avdata\\PythonTest\\ClipLidar\\Scratch"
arcpy.env.overwriteOutput = True
clipshape = "c:\\avdata\\PythonTest\\ClipLidar\\AOIPoly_SP.shp"
InputFolder = "c:\\avdata\\PythonTest\\ClipLidar\\TestData\\Input"
#find input folder bottom folder
InputString = InputFolder.split("\\")[-1]
#Find AOI extent
AOI = arcpy.mapping.Layer(clipshape)
AOIextent=AOI.getExtent()
#Intialize empty lists
rasterList=[]
#Find the input Grid feature classes
for root, dirs, files in os.walk(InputFolder):
for targetdir in dirs:
targetworkspace = root + "\\" + targetdir
#if the targetworkspace can be initialized as a Raster object, do nothing else
try:
arcpy.Raster(targetworkspace)
#if the targetworkspace fails when an attempt is made to initialize it as a raster object, search for rasters within it
except:
arcpy.env.workspace = targetworkspace
rasterInDir = arcpy.ListRasters("*", "GRID")
if len(rasterInDir) > 0:
for raster in rasterInDir:
rasterList.append(targetworkspace + "\\" + raster)
##for raster in rasterList:
## print raster
#start comparing extents and clipping
for name in rasterList:
FileObj = arcpy.mapping.Layer(name)
FileExtent = FileObj.getExtent()
#print "Processing: " + name
#if image overlap or touches AOI polygon, then clip image, write to output
Extent_Disjoint = str(FileExtent.disjoint(AOIextent))
if Extent_Disjoint == 'False':
print "Clipping: " + name
foldername = os.path.dirname(name)
OutputPath = foldername.replace(InputString,"Output")
if not os.path.exists(OutputPath):
os.makedirs(OutputPath)
if len(os.path.basename(name)) > 10:
Shortname = (os.path.basename(name))[0:11]
print "Grid name shortened to " + Shortname
Clipname = os.path.splitext(Shortname)[0] + "cp"
else:
Clipname = os.path.splitext(name)[0] + "cp"
Clipbasename = os.path.basename(Clipname)
print "Writing " + (os.path.join(OutputPath,Clipbasename))
arcpy.Clip_management(name,"#",(os.path.join(OutputPath,Clipbasename)),clipshape,"#","ClippingGeometry")
#Find the input tifs and shapefiles
for root, dirs, files in os.walk(InputFolder):
for name in files:
if os.path.splitext(name)[1] == ".tif":
FileObj = arcpy.mapping.Layer(os.path.join(root, name))
FileExtent = FileObj.getExtent()
#print "Processing: " + name
#if image overlap or touches AOI polygon, then clip image, write to output
Extent_Disjoint = str(FileExtent.disjoint(AOIextent))
if Extent_Disjoint == 'False':
print "Clipping: " + name
OutputPath = root.replace(InputString,"Output")
if not os.path.exists(OutputPath):
os.makedirs(OutputPath)
Clipname = os.path.splitext(name)[0] + "cp.tif"
print "Writing " + (os.path.join(OutputPath,Clipname))
arcpy.Clip_management((os.path.join(root, name)),"#",(os.path.join(OutputPath,Clipname)),clipshape,"#","ClippingGeometry")
#for shapefiles, use a different clip, write to output.
if os.path.splitext(name)[1] == ".shp":
print(os.path.join(root, name))
FileObj = arcpy.mapping.Layer(os.path.join(root, name))
FileExtent = FileObj.getExtent()
Extent_Disjoint = str(FileExtent.disjoint(AOIextent))
if Extent_Disjoint == 'False':
print "Clipping: " + name
OutputPath = root.replace(InputString,"Output")
if not os.path.exists(OutputPath):
os.makedirs(OutputPath)
Clipname = os.path.splitext(name)[0] + "cp"
print "Writing " + (os.path.join(OutputPath,Clipname))
arcpy.Clip_analysis((os.path.join(root, name)),clipshape,(os.path.join(OutputPath,Clipname)))
