我希望ISpatialFilter.SpatialRelDescription与IRelationalOperator.Relation保持一致。但是,当我运行下面的代码时,我得到:
Total features : 44368
Count1: 9 T********
Count2: 9 T********
Count1: 21 F********
Count2: 44359 F********
Count1: 1 ******F**
Count2: 1 ******F**
这是F********
关系中的错误吗?
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;
}
1
抱歉,没有答案,但地理空间幽默+1
—
bFlood 2010年