nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

Count number of features in FeatureClasses

Andreas_ForøTollefsen

Hi all.

What i want to do is to create a list of the number of features in each of my featureclasses.

I have already created a list of all the featureclasses in the workspace.

Now i need to iterate through all of these and count the number of features in each featureclass.

import arcgisscripting
import os
gp = arcgisscripting.create(9.3)

gp.Workspace = "H:\\geodata\\"

print "Generating list of featureclasses in" " " + gp.Workspace
fcList = gp.ListFeatureClasses()

for fc in fcList:
    print fc

Anyone have a suggestion? I assume i need to use the cursor object?

Trouver plus de publications taguées avec

Accepted answers

All comments

JohnRoberts1

Sub main()

Dim strServer As String
Dim strInstance As String
Dim strDB As String
Dim strUser As String
Dim strPass As String
Dim strVersion As String

strServer = "GISSDE"
strInstance = "sde:sqlserver:gissde\production"
strDB = "Utilities"
strUser = "gisreader"
strPass = "gisreader"
strVersion = "sde.DEFAULT"

        Dim pPropset As IPropertySet
        Set pPropset = New PropertySet
            With pPropset
            .SetProperty "Server", strServer
            .SetProperty "Instance", strInstance
            .SetProperty "Database", strDB
            .SetProperty "User", strUser
            .SetProperty "Password", strPass
            .SetProperty "Version", strVersion
            End With

            Dim pWorkspace As IWorkspace
            Dim pFact As IWorkspaceFactory
            Set pFact = New SdeWorkspaceFactory
            'pWOrkspace = pFact.Open pPropset, 0


            Set pWorkspace = pFact.Open(pPropset, 0)

            Dim pFeatws As IFeatureWorkspace
            Set pFeatws = pWorkspace

            Dim pVerWorkspace As IVersionedWorkspace
            Set m_pVersionWork = pWorkspace

            Dim sDatasetName As String
            sDatasetName = "<All>"
            If sDatasetName = "<All>" Then 'Checking all layers in all datasets
                Dim pWork As IWorkspace, pEnumData As IEnumDatasetName
                Dim pFeatDatasetName As IFeatureDatasetName, pDataset As IDatasetName
                Dim pEnumFeatClassName As IEnumDatasetName, pFeatClassName As IDatasetName
                Dim pVersionedObject As IVersionedObject
                Dim pFC As IFeatureClass
                Dim pName As IName
                Dim lFCount As Long
                Set pWork = m_pVersionWork


                'process data outside of feature datasets...
                Set pEnumData = pWork.DatasetNames(esriDTFeatureClass)
                Set pFeatClassName = pEnumData.Next
                Do While Not pFeatClassName Is Nothing 'Loop through the feature classes
                    ' "Checking " & pFeatClassName.Name & " layer for updates . . . "
                    Set pName = pFeatClassName
                    Set pFC = pName.Open
                    lFCount = pFC.FeatureCount(Nothing)

                            Debug.Print (pFeatClassName.Name) & " Has " & lFCount & " Features"

                    Set pFeatClassName = pEnumData.Next
                Loop


                'now process all feature datasets
                Set pEnumData = pWorkspace.DatasetNames(esriDTFeatureDataset)
                Set pDataset = pEnumData.Next
                Do While Not pDataset Is Nothing 'Loop through the dataset names
                    Set pFeatDatasetName = pDataset
                    Set pEnumFeatClassName = pFeatDatasetName.FeatureClassNames
                    Set pFeatClassName = pEnumFeatClassName.Next
                    Do While Not pFeatClassName Is Nothing 'Loop through the feature classes

                        Set pName = pFeatClassName
                        Set pFC = pName.Open
                        lFCount = pFC.FeatureCount(Nothing)
                        Debug.Print (pFeatClassName.Name) & " Has " & lFCount & " Features"
                        Set pFeatClassName = pEnumFeatClassName.Next
                    Loop

                    Set pDataset = pEnumData.Next
                Loop

            End If

End Sub

dgesridgesri

Greetings,

Is there a way to get all the featureclasses in a worksapce through vb.net code?

please advise.....

JustinGinn1

You have no idea how much time I wasted because of this problem... wow. THANK YOU SO MUCH! I thought I was going insane... I had imported the 9.3 geoprocessor on a new version of the script, and I could not figure out why it wasn't working.

BryceStath

You can see the counts of each featureclass in ArcCatalog without any Python scripting.
Add features count as a column in the options for display and then touch each featureclass to update the metadata.
then when you select the geodatabase or folder the counts will display in the listing.

thanks
to elaborate a bit..

In ArcCatalog (9.3.1), Tools -> Options -> Contents(tab) -> Which Metadata... (2nd pane) -> # of Features

KevinHibma

The help topic is here:
http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/Result/000v000000n7000000/

Sample code:

import arcpy
inTable = arcpy.GetParameterAsText(0)
result = arcpy.GetCount_management(inTable)
print result.getOutput(0)

BrianWilson

Thanks Joel for pointing that out. My python scripts seem to break with each new ArcGIS release! It is very time consuming updating even a small number of mission critical python scripts with testing, etc. But with each year I have more and more scripts to maintain.

ESRI doesn't think code we wrote last year is important apparently. It also makes exchanging scripts with our clients impossible. Oh well.

Does anyone know where the "gcResult" object returned by GetCount is documented? The page for "Get Count" has no indication of what is being returned.

http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//0017000000n7000000.htm

Thanks Brian

RandyKreuziger

I found this out when my code failed when we switched to 9.3. I don't remember ESRI announcing this change like they did the change to python lists so it was a little bit of a suprise.

Joel

Thanks Joel for pointing that out. My python scripts seem to break with each new ArcGIS release! It is very time consuming updating even a small number of mission critical python scripts with testing, etc. But with each year I have more and more scripts to maintain.

KimOllivier

You can see the counts of each featureclass in ArcCatalog without any Python scripting.
Add features count as a column in the options for display and then touch each featureclass to update the metadata.
then when you select the geodatabase or folder the counts will display in the listing.

featurecount.png

JoelCalhoun

Per the help example, the getCount method changed from 9.2 to 9.3.
So if you are using the old 9.2 geoprocessor you can pass the resulting count directly to a variable.
If you are using the 9.3 geoprocessor you have to extract the value from a result object as the previous poster mentioned.

ArcGIS Script example:

# Using a 9.2-version geoprocessor
#
import arcgisscripting
gp = arcgisscripting.create()

gp.workspace = "D:/Workspace"
count = gp.GetCount_management("roads")

# Using a 9.3-version geoprocessor
#
import arcgisscripting
gp = arcgisscripting.create(9.3)

gp.workspace = "D:/Workspace"
result = gp.GetCount_management("roads")
count = int(result.GetOutput(0))

I found this out when my code failed when we switched to 9.3. I don't remember ESRI announcing this change like they did the change to python lists so it was a little bit of a suprise.

Joel

ChrisMathers

I asked the same recently and was told I would need to do count=gp.getcount(file) then count=int(Count.getOutput(0)) because what you are given as a result of the getcount tool is a geoprocessing result object not a number. I didnt actually need to though. For some reason doing count=gp.getcount(file) then featurecount=count worked fine. I really dont know why as that should make featurecount just another link to the result object.

Andreas_ForøTollefsen

Thanks Dan.

That should be an easier way of doing it.
I wrote the script below, and it did work. However, I think getCount would be better.
Thanks.

import arcgisscripting
import os
gp = arcgisscripting.create(9.3)

gp.Workspace = "H:\\geodata\\"

print "Generating list of featureclasses in " + gp.Workspace
fcList = gp.ListFeatureClasses()

print "Listing all featureclasses in " + gp.Workspace
for fc in fcList:
    print fc

print "listing fields in featureclasses"
for fc2 in fcList:
    fn = gp.ListFields(fc2)
    print str(fc2) + " " + str(len(fn))

print "Counting number of records in each featureclass"
for fc3 in fcList:
    cur = gp.SearchCursor(fc3)
    row = cur.Next()
    alist = []
    while row:
        row = cur.Next()
        alist.append(row)
    alist.sort()
    print "number of features in " + fc3 + " are: " + str(len(alist))
                

del fcList, fc, fc2, fn, cur, row, alist

DanPatterson_Retired

perhaps, getcount
http://webhelp.esri.com/arcgisdesktop/9.3/index.cfm?id=1967&pid=1961&topicname=Get_Count_%28Data_Management http://webhelp.esri.com/arcgisdesktop/9.3/index.cfm?TopicName=Get_Count_(Data_Management)

Rejoindre la conversation

Les membres connectés peuvent publier, suivre les mises à jour, et plus encore. Nouveau ici ? Inscrivez-vous gratuitement.

Liens rapides

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.