Select by location using VBA

JerryBiedenbender — Fri, 01 Oct 2010 15:04:50 GMT

Hello,

   I am trying to write some code that will select a polygon from layer1 "SSESarea" and then select by location all polygons from layer2 "Parcels" with their centroid within the polygon selected from layer1.

I have the first part of the code complete and working but I could use some advice on the select by location part.

The following code selects the polygon that I need from Layer1, then opens the select by location tool in its own window like it does when you open it in arcmap . I would like to have the select by location tool automatically run instead of poping up on the screen.

Private Sub UIButtonControl1_Click()
' SELECT BY ATTRIBUTES
Dim pMxDoc As IMxDocument
Dim pMap As IMap
Dim pActiveView As IActiveView
Dim pFeatureLayer As IFeatureLayer
Dim pFeatureSelection As IFeatureSelection
Dim pQueryFilter As IQueryFilter

Set pMxDoc = Application.Document
Set pMap = pMxDoc.FocusMap
Set pActiveView = pMap

If Not TypeOf pMap.Layer(0) Is IFeatureLayer Then Exit Sub
Set pFeatureLayer = pMap.Layer(0)
Set pFeatureSelection = pFeatureLayer 'QI

'Create the query filter
Set pQueryFilter = New QueryFilter
pQueryFilter.WhereClause = "ProjectName = 'LSP - Arcadia'"

'Invalidate only the selection cache
'Flag the original selection
pActiveView.PartialRefresh esriViewGeoSelection, Nothing, Nothing
'Perform the selection
pFeatureSelection.SelectFeatures pQueryFilter, esriSelectionResultNew, False
'Flag the new selection
pActiveView.PartialRefresh esriViewGeoSelection, Nothing, Nothing

' ZOOM TO SELECTED FEATURES
Dim pCmdItem As ICommandItem

    Set pCmdItem = Application.Document.CommandBars.Find(ArcID.Query_ZoomToSelected)
    pCmdItem.Execute

' MAP SCALE
'Dim pMxDoc As IMxDocument
'Dim pMap As IMap

    Set pMxDoc = ThisDocument
    Set pMap = pMxDoc.FocusMap
    pMap.MapScale = 3000
    pMxDoc.ActiveView.Refresh

' EXPORT TO JPEG
'Dim pMxDoc As IMxDocument (it is necessary while using separately)
'Dim pActiveView As IActiveView (it is necessary while using separately)
Dim pExport As IExport
Dim pPixelBoundsEnv As IEnvelope
Dim exportRECT As tagRECT
Dim iOutputResolution As Integer
Dim iScreenResolution As Integer
Dim hDC As Long

Set pMxDoc = ThisDocument
Set pActiveView = pMxDoc.ActiveView
Set pExport = New ExportJPEG

pExport.ExportFileName = "C:\ExportMap." & Right(pExport.Filter, 3)
iScreenResolution = 150 'default screen resolution is usually 96dpi
iOutputResolution = 150
pExport.Resolution = iOutputResolution

With exportRECT
.Left = 0
.Top = 0
.Right = pActiveView.ExportFrame.Right * (iOutputResolution / iScreenResolution)
.bottom = pActiveView.ExportFrame.bottom * (iOutputResolution / iScreenResolution)
End With

Set pPixelBoundsEnv = New Envelope
pPixelBoundsEnv.PutCoords exportRECT.Left, exportRECT.Top, exportRECT.Right, exportRECT.bottom
pExport.PixelBounds = pPixelBoundsEnv

hDC = pExport.StartExporting
pActiveView.Output hDC, pExport.Resolution, exportRECT, Nothing, Nothing
pExport.FinishExporting
pExport.Cleanup

'*************
'Pops up select by location tool
Set pCmdItem = Application.Document.CommandBars.Find(Query_SelectByLayer)
pCmdItem.Execute

'**********************
End Sub

I found the following code on another forum and think it might be a start but I cannot figure out how to configure it.

************************************************************************
Since "have their center in" is unfortunately not an esri spatial relation that you can use with a spatial filter, I use a work around that will get the same results probably 99.99% of the time in arcobjects. Basically I get a small negative (inside) buffer on the geometry and then create a spatial filter using esriSpatialRelIntersects. The negative buffer distance is very, very small; ~1 meter. You may have to change this value depending on the map units of your data. You can see how this works in the code below, where pGeom is the polygon that I want to search within, and pFeatureLayer is a featurelayer that has the features you want to select

Dim pTopoOp as ITopologicalOperator
Set pTopoOp = pGeom
Dim pSpQuery as ISpatialFilter
Set pSpQuery = New SpatialFilter
With pSpQuery
'the value inside buffer should be very small map unit
Set .Geometry = pTopoOp.Buffer(-1)
.SpatialRel = esriSpatialRelIntersects
.GeometryField = pFeatureLayer.FeatureClass.ShapeFieldName
End With
Dim pFeatSelection as IFeatureSelection
Set pFeatSelection = pFeatureLayer
'Add features to current selection
pFeatSelection.SelectFeatures pSpQuery, esriSelectionResultsAdd, False

**************************************************************************

Re: Select by location using VBA

DuncanHornby — Fri, 01 Oct 2010 20:45:28 GMT

Hi Jerry,

If you want to do what the geo processing tool does why not simply use it? In VBA you access it via the IGeoprocessor interface. The following link takes you to a thread where IGeoProcessor is calling the SelectByLocation tool. You would simply re-jig the parameters array with your layers and selection criteria.

http://forums.esri.com/Thread.asp?c=93&f=983&t=249267&mc=4#msgid762216

Duncan

topic Re: Select by location using VBA in ArcObjects SDK Questions

Select by location using VBA

Re: Select by location using VBA