arcpy实现空间查询_空间查询的实现

Dim pMap As IMap

Dim pPoint As IPoint

Set pMap = MapControl1.Map

Dim pFeatureLayer As IFeatureLayer

Set pFeatureLayer = pMap.Layer(1)

Set pPoint = MapControl1.ActiveView.ScreenDisplay.DisplayTransformation.ToMapPoint(x, y)

Dim pEnv As IEnvelope

Set pEnv = New Envelope

Set pEnv = MapControl1.ActiveView.Extent

pEnv.Height = 100

pEnv.Width = 100

pEnv.CenterAt pPoint

MsgBox pPoint.x

MsgBox pEnv.XMax

Dim pGeometry As IGeometry

Set pGeometry = pEnv

Set pGeometry.SpatialReference = pMap.SpatialReference

pFeatureLayer.Selectable = True

Dim pSFilter As ISpatialFilter

Set pSFilter = New SpatialFilter

With pSFilter

Set .Geometry = pGeometry

.GeometryField = pFeatureLayer.FeatureClass.ShapeFieldName

.SpatialRel = esriSpatialRelIntersects

End With

Dim b As Boolean

b = pSFilter.FilterOwnsGeometry

Dim pFeatureCursor As IFeatureCursor

Set pFeatureCursor = pFeatureLayer.Search(pSFilter, False)