Bug dans la relation spatiale F ********? (Je ne plaisante pas)

13

Je m'attendrais à ce que ISpatialFilter.SpatialRelDescription se comporte de manière cohérente avec IRelationalOperator.Relation . Cependant, lorsque j'exécute le code ci-dessous, j'obtiens:

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

Est-ce un bug dans la 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;
}
Kirk Kuykendall
la source
1
désolé, pas de réponse mais +1 pour l'humour géospatial
bFlood

Réponses:

3

Cela semble un peu étrange. Peut-être qu'il n'a été testé qu'avec quelques opérations spécifiques, Ex: Intersects, Disjoint, etc? Que fait FF*FF****(inverse de T********) pour vous? Il me semble que Within ne fonctionnait pas avec quelques versions en arrière donc ça ne me surprend pas.

Dandy
la source
1
Cela semble fonctionner. Si l'espace est si spécial, il semble bien qu'il devrait recevoir une meilleure assurance qualité.
Kirk Kuykendall