Elencare tutte le classi di caratteristiche nel file Geodatabase, comprese le serie di dati delle caratteristiche?

Da Python come posso creare un elenco di tutte le classi di caratteristiche in un file geodatabase (* .gdb), inclusi i set di dati delle caratteristiche interne? L' esempio standard elenca solo le classi di funzionalità al livello superiore del geodatabase:

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

gp.workspace = 'd:\scratch.gdb'
fcs = gp.ListFeatureClasses()

for fc in fcs:
    print fc

Indica a quale versione di ArcGIS Desktop si applica la tua risposta (Sto cercando 9.3 ma potremmo anche raccogliere tutte le versioni in un unico posto).

Questa routine per arcgis10 restituisce tutti i fcs (autonomi OR all'interno di un set di dati di funzionalità) all'interno di un gdb. Basta impostare il tuo arcpy.env.workspace quindi fare un ciclo for

def listFcsInGDB():
    ''' set your arcpy.env.workspace to a gdb before calling '''
    for fds in arcpy.ListDatasets('','feature') + ['']:
        for fc in arcpy.ListFeatureClasses('','',fds):
            yield os.path.join(arcpy.env.workspace, fds, fc)

Ho finito per usare la risposta di Gotchula , ma senza resa perché generalmente riutilizzo le maniglie FC create e le rese sono usate una volta poi scartate , per me è più facile leggere e capire cosa fcs.append()sta facendo fcs = yield(...).

def listFcsInGDB(gdb):
    ''' list all Feature Classes in a geodatabase, including inside Feature Datasets '''
    arcpy.env.workspace = gdb
    print 'Processing ', arcpy.env.workspace

    fcs = []
    for fds in arcpy.ListDatasets('','feature') + ['']:
        for fc in arcpy.ListFeatureClasses('','',fds):
            #yield os.path.join(fds, fc)
            fcs.append(os.path.join(fds, fc))
    return fcs

gdb = sys.argv [1]
fcs = listFcsInGDB(gdb)
for fc in fcs:
    print fc            

risultati:

d:\> python list-all-fc.py r:\v5\YT_Canvec.gdb
Processing  r:\v5\YT_Canvec.gdb
Buildings_and_structures\BS_2530009_0
Buildings_and_structures\BS_2380009_2
Buildings_and_structures\Tower
Buildings_and_structures\Underground_reservoir
...

Questo è ora in un modulo che chiamo arcplus *. Inserire con l'altro codice o PYTHONPATH e quindi:

import arcplus
fcs = arcplus.listAllFeatureClasses('d:\default.gdb')
for fc in fcs:
    print "magic happens with: ", fc

Arcplus aggiunge anche il filtro jolly; per elaborare solo le classi di funzionalità che iniziano con "HD_" all'interno di set di dati di funzionalità contenenti "Hydro"

fcs = arcplus.listAllFeatureClasses(gdb, fd_filter='*Hydro*', fc_filter='HD_*')

. * ora su Github, aggiornato per 10.x. Per arcgis 9.3 vedi qui .

Mi rendo conto che questa domanda è taggata 9.3, ma chiunque cerchi la stessa risposta dalla 10.1 in poi sta meglio usando arcpy.da.Walk . È più veloce e più preciso di ListDatasets / FeatureClasses / Rasters / ecc.

import arcpy
import os

for root, dirs, datasets in arcpy.da.Walk('d:\scratch.gdb'):
    for ds in datasets:
        print os.path.join(root, ds)

La funzione walk funziona allo stesso modo della passeggiata di python . Esegue l'iterazione attraverso le directory nel percorso specificato e ad ogni iterazione, root rappresenta il percorso completo della directory e directory e set di dati sono elenchi delle sottodirectory e dei file contenuti all'interno.

Quando si percorre un geodatabase, i set di dati delle funzioni vengono trattati allo stesso modo delle directory. Se si desidera elencare solo i set di dati e i set di dati delle funzionalità nella cartella principale e non aprire i set di dati delle funzionalità per visualizzare i contenuti, è possibile eseguire:

for root, dirs, datasets in arcpy.da.Walk('d:\scratch.gdb'):
    print 'feature datasets:'
    for fds in dirs:
        print os.path.join(root, fds)
    print 'datasets:'
    for ds in datasets:
        print os.path.join(root, ds)
    break

Il ListDatasets metodo è quello che penso tu stia cercando. Questo FGDB ha un FD chiamato "Wells" e contiene 3 FC. Per 9.3.1:

C:\WINDOWS\system32>python
Python 2.5.1 (r251:54863, Apr 18 2007, 08:51:08) [MSC v.1310 32 bit (Intel)] on
win32
Type "help", "copyright", "credits" or "license" for more information.
>>> import arcgisscripting
>>> gp=arcgisscripting.create(9.3)
>>> d='c:\data\Third_Party_Wells_PRD.gdb'
>>> gp.workspace = d
>>>
>>> fds=gp.ListDatasets('','Feature')
>>> for fd in fds:
...     print fd
...
Wells
>>> for fd in fds:
...     gp.workspace=d + '/' + fd
...     fcs=gp.ListFeatureClasses()
...     for fc in fcs:
...             print fc
...
Third_Party_Wells_BHL
Third_Party_Wells_LAT
Third_Party_Wells_SHL
>>>