An error was encountered during an IFeatureSet.Intersection operation. The code to duplicate the error is as follows:
```
//Open the three data sets
IFeatureSet siteDesignArea = FeatureSet.OpenFile(pathToSiteDesignArea);
IFeatureSet bufRecons= FeatureSet.OpenFile(pathToBufRecons);
IFeatureSet fullHucCatchments = FeatureSet.OpenFile(pathToFullCatchmentLayer);
//Project the HUC Catchments to the same projection as siteDesignArea file
fullHucCatchments.Reproject(siteDesignArea.Projection);
//Clip the HUC Catchments to the BufRecons shapefile
IFeatureSet ClippedFeatures = fullHucCatchments.Intersection(bufRecons, FieldJoinType.LocalOnly, null);
```
The error code message is:
no outgoing dirEdge found [ DotSpatial.Topology.Coordinate ]
I've attached the three required shapefiles (shapefiles.zip) and some images showing how the result is supposed to look (images.zip). Note that the result should be two large polygons with several smaller sliver polygons. .
Note that the "FullHUCCatchment" is in the wrong projection, so we require to be able to reproject it first then intersect it.
Comments: Seems that the initial issue was addressed through a workaround. There is a merge problem that could be looked at on a separate issue.
```
//Open the three data sets
IFeatureSet siteDesignArea = FeatureSet.OpenFile(pathToSiteDesignArea);
IFeatureSet bufRecons= FeatureSet.OpenFile(pathToBufRecons);
IFeatureSet fullHucCatchments = FeatureSet.OpenFile(pathToFullCatchmentLayer);
//Project the HUC Catchments to the same projection as siteDesignArea file
fullHucCatchments.Reproject(siteDesignArea.Projection);
//Clip the HUC Catchments to the BufRecons shapefile
IFeatureSet ClippedFeatures = fullHucCatchments.Intersection(bufRecons, FieldJoinType.LocalOnly, null);
```
The error code message is:
no outgoing dirEdge found [ DotSpatial.Topology.Coordinate ]
I've attached the three required shapefiles (shapefiles.zip) and some images showing how the result is supposed to look (images.zip). Note that the result should be two large polygons with several smaller sliver polygons. .
Note that the "FullHUCCatchment" is in the wrong projection, so we require to be able to reproject it first then intersect it.
Comments: Seems that the initial issue was addressed through a workaround. There is a merge problem that could be looked at on a separate issue.