Hfrhyu

Using ArcGIS 10.4.1 ArcObjects with geoprocessing. I have an input feature class with 100s of polygons and also one raster. I iterate through each polygon and clip the raster and then do additional work. It has performed well up to today. It seems some our polygons are outside of the raster area.

So my code in VB.NET looks like this:

'create an IVariantArray to hold the parameter values.
Dim parameters As IVariantArray = New VarArray

'populate the variant array with parameter values.
parameters.Add(inSurfaceFullName) 'input surface
parameters.Add("") 'rectangle to be used as the clipper
parameters.Add(outClippedRasterFullName) 'the output raster (ex Raster_3)
parameters.Add(inClippingPolyFullName) 'the input template dataset to be used as a clipper. in this case our sub feature class (NOT the input feature class)
parameters.Add("") 'all the pixels with the specified value will be set to NoData in the output raster dataset.
parameters.Add("ClippingGeometry") 'the raster dataset is clipped based on the perimeter of the polygon shape.

'execute the model tool by name.
Dim res As IGeoProcessorResult = GP.Execute("Clip_management", parameters, Nothing)

When it gets to the polygon that does not intersect the raster I get this error:

System.Runtime.InteropServices.COMException (0x80004005): Error HRESULT E_FAIL has been returned from a call to a COM component.
 at ESRI.ArcGIS.Geoprocessing.GeoProcessorClass.Execute(String Name, IVariantArray ipValues, ITrackCancel pTrackCancel)

I have been looking in the ESRI SDK for a pre-process where I can check if the polygon intersects the raster with a True/False as the result. Is that something that exists? Or another way to check the intersection of a geometry with a raster?

The only thing I have tried is to check if the extents intersect but I can see where this may fail. For example a circle geometry just at the corner of the raster. the envelopes will intersect but the geometry and raster will not.

With help from @MichaelStimson, I ended up testing the raster envelope against the polygon geometry. Seems to work rather well. Here is the function I use to test all features in a feature class against the raster:

''' <summary>
''' this function will determine if the polygon feature class is outside of the raster boundary. it will test all features in the
''' feature class and return one result for the overall disjoint result.
''' <param name="pInSurfaceEnv">the input surface envelope</param>
''' <param name="pPolyDsName">the input polygon Dataset name object</param>
''' </summary>
''' <returns></returns>
''' <remarks></remarks>
Private Function IsFeatureClassDisjointWithRaster(pInSurfaceEnv As IEnvelope, pPolyDsName As IDatasetName) As Boolean
 Try
 Dim pRelOp As IRelationalOperator2 = CType(pInSurfaceEnv, IRelationalOperator2)
 Dim pFeatClass As IFeatureClass = OpenFileFeatClassDataset(pPolyDsName.WorkspaceName.PathName, pPolyDsName.Name)
 Dim pFeatCursor As IFeatureCursor = pFeatClass.Search(Nothing, True)

 Dim IsFeatClassDisjoint As Boolean = True 'we assume the entire feature class is disjoint
 Dim pFeature As IFeature = pFeatCursor.NextFeature
 Do While Not pFeature Is Nothing
 Dim IsFeatureDisjoint As Boolean = pRelOp.Disjoint(pFeature.ShapeCopy)

 'we only need one feature to be non disjoint to make the entire feature class non disjoint
 If Not IsFeatureDisjoint Then IsFeatClassDisjoint = False

 'next feature
 pFeature = pFeatCursor.NextFeature
 Loop

 Return IsFeatClassDisjoint
 Catch ex As Exception
 Throw
 End Try
End Function

First you can start by obtaining the extent of the raster (IRasterDataset) with IGeodataset.Extent which returns an IEnvelope which can be cast as an IRelationalOperator.

For each polygon geometry (IGeometry) then test your relational operator for Disjoint which is quicker than trying to find Overlaps. If IRelationalOperator.Disjoint then the polygon does not intersect the extent of the raster so skip to the next one.