How to write a code that counts the number of points inside each polygon?

Question

Hello,

How to write a code VB ??? Arcobject that counts the number of points inside each polygon?

Please, see the attached image that shows my question

I???ve two layers, point layer and polygon layer, I need to write a code that calculate the number of points inside each polygon

Thank you

kalid

Clip_8.jpg

july2019.PNG

Amarz_0-1761246091827.png

vairaselvan · Answer

i need tat database can u send that

Hello,

How to write a code VB �?? Arcobject that counts the number of points inside each polygon?

Please, see the attached image that shows my question

I�??ve two layers, point layer and polygon layer, I need to write a code that calculate the number of points inside each polygon

Thank you

kalid

Venkata_RaoTammineni · Answer

Dear Venkat,

Thank you so much for the answer. It�??s really helpful

Sure. The code you have developed is working perfectly.

1. how to develop the code so that we don�??t need:
�?� To work in edit session.
�?� To have the field �??name�?� in the district layer.

2. How to develop the code considering the following approach:
�?� to intersect the two layers (districts and hospitals: both are polygons)
�?� And then to make summery in the table of the resultant layer (from the intersection) so that we count all the hospitals in each district.

I�??m supposed to print the result in the form and then to be able to export it to a table to excel

Please, see the attached file

Appreciated

Kalid

Hi,

Whatever the code i have given is the intersect two poly and point layers and count the point within the polygon and add into polygon layer table ...and below is that thread where you can find the code for exporting into excel...Please change the code as per your requirement.

http://forums.esri.com/Thread.asp?c=93&f=992&t=227281&mc=6#msgid692161

JamalNUMAN · Answer

Dear Venkat,

Thank you so much for the answer. It�??s really helpful

Sure. The code you have developed is working perfectly.

1. how to develop the code so that we don�??t need:
�?� To work in edit session.
�?� To have the field �??name�?� in the district layer.

2. How to develop the code considering the following approach:
�?� to intersect the two layers (districts and hospitals: both are polygons)
�?� And then to make summery in the table of the resultant layer (from the intersection) so that we count all the hospitals in each district.

I�??m supposed to print the result in the form and then to be able to export it to a table to excel

Please, see the attached file

Appreciated

Kalid

Hi,

please check below code...

Sub test()
    Dim pMxDoc As IMxDocument
    Set pMxDoc = ThisDocument
    
    Dim pPolyFLayer As IFeatureLayer
    Set pPolyFLayer = pMxDoc.FocusMap.Layer(1)
    
    Dim pPointFlayer As IFeatureLayer
    Set pPointFlayer = pMxDoc.FocusMap.Layer(0)
    
    Dim pEditor As IEditor, pUID As New UID
    pUID.Value = "esricore.Editor"
    
    Set pEditor = Application.FindExtensionByCLSID(pUID)
    If pEditor.EditState = esriStateEditing Then
        pEditor.StartOperation
        CountFeatures pPolyFLayer.FeatureClass, pPointFlayer.FeatureClass, _
                     esriSpatialRelContains, "NAME"
        pEditor.StopOperation "PolygonDensity"
    Else
        MsgBox "no edit session"
    End If
End Sub

Sub CountFeatures(pFC1 As IFeatureClass, _
                 pFC2 As IFeatureClass, _
                 lSpatialRel As esriSpatialRelEnum, _
                 strCountFld As String)
                 
    Dim pFCur As IFeatureCursor
    Set pFCur = pFC1.Search(Nothing, False)
    
    Dim lFld As Long
    lFld = pFC1.FindField(strCountFld)
    If lFld = -1 Then
        MsgBox "field not found: " & strCountFld
        Exit Sub
    End If
    
    Dim pFeat As IFeature
    Set pFeat = pFCur.NextFeature
    Do While Not pFeat Is Nothing
    
        Dim pSF As ISpatialFilter
        Set pSF = New SpatialFilter
        Set pSF.Geometry = pFeat.ShapeCopy
        pSF.SpatialRel = lSpatialRel
        pSF.GeometryField = pFC2.ShapeFieldName
        
        pFeat.Value(lFld) = pFC2.FeatureCount(pSF)
        MsgBox "OID: " & pFeat.OID & ", count: " & pFeat.Value(lFld)
        pFeat.Store
        Set pFeat = pFCur.NextFeature
    Loop

End Sub

is that your looking for ?...

Thanks and Regards,

Venkat

hospita_district.jpg

image 2.png