Tuesday, 6 March 2018

arcobjects - Bug in F******** spatial relation? (I'm not cussing)


I would expect ISpatialFilter.SpatialRelDescription to behave consistently with IRelationalOperator.Relation. However, when I run the code below, I get:


    Total features : 44368
Count1: 9 T********
Count2: 9 T********
Count1: 21 F********
Count2: 44359 F********
Count1: 1 ******F**
Count2: 1 ******F**


Is this a bug in the F******** relation ?


private void TestSpatialFilter()
{
IMxDocument mxdoc = ((IApplication)m_application).Document as IMxDocument;
IFeatureLayer fLayer = mxdoc.FocusMap.get_Layer(0) as IFeatureLayer;
IEnvelope env = ((IActiveView)mxdoc.FocusMap).Extent;
env.Expand(0.5, 0.5, true);
env.Project(((IGeoDataset)fLayer.FeatureClass).SpatialReference);
Debug.Print("Total features : {0}", fLayer.FeatureClass.FeatureCount(null));
IFeatureClass fc = fLayer.FeatureClass;

string[] rels = { "T********", "F********", "******F**" };
foreach (string rel in rels)
{
Debug.Print("Count1: {0} {1}", Count1(fc, env, rel), rel);
Debug.Print("Count2: {0} {1}", Count2(fc, env, rel), rel);
}
}

public static int Count1(IFeatureClass fc, IGeometry geom, string relDescription)
{

ISpatialFilter sf = new SpatialFilterClass();
sf.SpatialRel = esriSpatialRelEnum.esriSpatialRelRelation;
sf.SpatialRelDescription = relDescription;
sf.Geometry = geom;
int i = fc.FeatureCount(sf);
Marshal.FinalReleaseComObject(sf);
return i;
}
public static int Count2(IFeatureClass fc, IGeometry geom, string relDescription)
{

int i = 0;
string relation = String.Format("RELATE(G1,G2,'{0}')", relDescription);
IRelationalOperator relOp = geom as IRelationalOperator;

IFeatureCursor fCur = fc.Search(null, false);
IFeature feat;
while ((feat = fCur.NextFeature()) != null)
{
if(relOp.Relation(feat.Shape,relation))
i++;

}
Marshal.FinalReleaseComObject(fCur);
return i;
}

Answer



That does seem a bit odd. Perhaps it was only tested with few specific operations, Ex: Intersects, Disjoint, etc? What does FF*FF**** (inverse of T********) do for you? I seem to remember Within not working a few versions back so it does not surprise me.


No comments:

Post a Comment

arcpy - Changing output name when exporting data driven pages to JPG?

Is there a way to save the output JPG, changing the output file name to the page name, instead of page number? I mean changing the script fo...