Comparing Raster Extents

JonathanMulder

Greetings,

I'm very new to Python having just taken ESRI's 3-day Python scripting course.

I'm trying to clip a large raster into smaller rasters that are the size of DOQQs.  I found some code that compares raster extents, then if the extents intersect, it clips the raster.

My code has issues with this line: "If DOQQExtent.within(TargetRasterExtent.extent):"  This is my complete code:

import arcpy

# Get directories of DOQQ images
in_raster = arcpy.GetParameterAsText(0)
in_directories = arcpy.GetParameterAsText(1)
out_raster = arcpy.GetParameterAsText(2)
TargetRasterExtent = arcpy.Describe(in_raster).extent
for CurrentDirectory in in_directories:
    DOQQRasters = arcpy.ListFiles("*.tif")
    for CurrentRaster in CurrentDirectory:
        DOQQExtent = arcpy.Describe(CurrentRaster).extent
        If DOQQExtent.within(TargetRasterExtent.extent):
#        if arcpy.Describe(CurrentRaster).extent.overlaps(TargetRasterExtent) == 'True':
            print "Extent Overlaps"

Any idea what's wrong with the line?

Thanks for any help and Happy Thanksgiving!

Jonathan Mulder
California Departement of Water Resources

JonathanMulder

Update: I think I'm getting closer but still not getting the clips.  I've modified my code as follows, but no clipped rasters are generated.  I suspect it's in my line comparing extents and the Clip_management function doesn't fire off.

Any insight?

import arcpy
#from arcpy import env
#from arcpy.sa import*
arcpy.CheckOutExtension("spatial")

# Get name of raster mosaic
in_rastermosaic = arcpy.GetParameterAsText(0)
# Get directories of DOQQ images
in_directories = arcpy.GetParameterAsText(1)
# Get directories of outputted rasters
out_raster = arcpy.GetParameterAsText(2)

TargetRasterExtent = arcpy.Describe(in_rastermosaic).extent
for CurrentDirectory in in_directories:
    DOQQRasters = arcpy.ListFiles("*.tif")
    for CurrentRaster in DOQQRasters:
        DOQQExtent = arcpy.Describe(CurrentRaster).extent
        if DOQQExtent.within(TargetRasterExtent.extent):
            arcpy.Clip_management(in_rastermosaic,"#","1958_" + CurrentRaster,"#",0,"NONE")
    #        if arcpy.Describe(CurrentRaster).extent.overlaps(TargetRasterExtent) == 'True':
            print "Extent Overlaps"
        else:
            print "Extent does not Overlaps"