arcgis gp 选择图层_ArcGIS Engine实现图层间空间选择的优化策略

转自：https://github.com/xinying180/QueryByLayer

如果您是ArcGIS Engine开发人员，也许会有这样的困惑：为什么对两个要素图层进行空间选择，ArcMap中瞬间就出结果了，而Engine中则慢很多倍，尤其是当数据量大时，该速率甚至无法忍受。图层间如何实现高效的空间选择呢？相信阅读完下面的文章，答案会迎刃而解。

下面就带着问题来开始今天的讨论吧。

问题：

假如有一个居民点数据和一个建筑物数据，想要知道哪些居民点被建筑物所覆盖，如何实现？

答案：

ArcMap中如何实现？

ArcMap中实现此功能很简单，即使用菜单条上的Select By Location或者Select Layer By Location工具，数据(点要素类中含有600个点，面要素类中含有750个面)和结果如下图所示。

用时：0.06秒

ArcGIS Engine中如何实现呢？

ArcGIS Engine中(相信很多用户是不太喜欢直接调用GP工具的)如何实现该功能呢？一般情况下会考虑使用ISpatialFilter结合遍历面要素来实现：

ISpatialFilter可以进行空间查询和属性查询。其Geometry属性传入要查询的几何，仅可使用high-level geometry(如polygon、polyline、points、multipoints)、envelope和geometry bags，这里使用面要素的geometry；其GeometryField属性用于指定过滤所使用的几何字段，这里使用点要素类的几何字段pointFeatureClass.ShapeFieldName；其SpatialRel属性用于指定空间关系，公式为[query_geometry] [spatial_relationship] [feature]，这里的query_geometry为面，feature为点要素，空间关系就应该是包含，所以应该使用esriSpatialRelContains(注意与ArcMap中有所区别，ArcMap中需设置空间关系为CompletelyWithin)。SubFields可以指定查询后返回的属性字段，以逗号分隔。如果不设置的话会默认为“*”，即返回所有字段。如果需要获取要素属性值得话，最好设置该字段，可以提高性能。当然如果不获取属性的话，那这个字段可设可不设，比如选择要素。此外，如果还想同时进行属性查询，直接设置WhereClause语句即可。主要代码：

//从MapControl中获取的点图层

IFeatureLayer pointFeatureLayer = axMapControl1.get_Layer(0) as IFeatureLayer;

IFeatureSelection pointFeatureSelection = pointFeatureLayer as IFe