Latest Contributions by JamalNUMAN

thank you for the answer. sure, i did...please, see the attached images Jamal The clipped raster has the same cell numbers as the original one!!! Why? Please, see the attached image that shows the case (particularly the tables!!!) The clipped raster has the same cell numbers as the original one (exactly the same) despite the fact that the clipped one is just a part of the original raster!!! Of course�?� really weird!!!!! What might be the problem? Thank you Regards Jamal ... View more

Hello Venkat, Thank you very much for the code. It would be very much appreciated if you develop the code so that: �?� First to intersect the layer of hospital with the layer of district (both are polygon) �?� Then to have the frequency in the table of the intersected layer so that we know the number of hospitals in each district In this case: We don�??t need to start editing and to have the field �??name�?� in the district layer in order to get the result (as it�??s in the provided code) Thank you Regards Kalid hi, How do you intersect point and polygon...You are looking for Points within the polygon ?if that is true i have already given you the code that perfectly working...please check previous forms Thanks ... View more

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 ... View more