Remove rasters from a mosaic dataset using an area of interest polygon in Python

Summary

The Remove rasters from Mosaic Dataset tool does not have the option for spatial query. The code below illustrates the use of ArcPy to remove rasters based on a polygon area of interest.

Procedure

1. 
   

   Code:
   import arcpy
   from arcpy import env
   
   #set variables
   env.workspace = r"D:\Documents\Northwest.gdb"
   clipShape = "parks"
   mosaicDataset = "CityOfPortland"
   
   # Make a mosaic layer to be able to select the footprint
   arcpy.MakeMosaicLayer_management(mosaicDataset, "mdlayer")
   
   # Select the footprint based on polygon clipShape using the intersect method
   arcpy.SelectLayerByLocation_management("mdlayer", "INTERSECT", clipShape)
   
   # Run the remove rasters tool without any query on the selection
   arcpy.RemoveRastersFromMosaicDataset_management("mdlayer", "#", "UPDATE_BOUNDARY", "MARK_OVERVIEW_ITEMS")