Mendaftar semua kelas fitur di File Geodatabase, termasuk di dalam dataset fitur?

Dari python bagaimana saya bisa membuat daftar semua kelas fitur dalam file geodatabase (* .gdb), termasuk di dalam dataset fitur? Contoh standar hanya mencantumkan kelas fitur di tingkat atas geodatabase:

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

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

for fc in fcs:
    print fc

Tolong tunjukkan versi ArcGIS Desktop mana jawaban Anda berlaku (Saya mencari 9.3 tetapi kami mungkin juga mengumpulkan semua versi di satu tempat).

Rutin ini untuk arcgis10 mengembalikan semua fcs (mandiri atau dalam dataset fitur) di dalam gdb. Cukup atur arcpy.env.workspace Anda lalu lakukan for loop

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)

Saya akhirnya menggunakan jawaban gotchula , tetapi tanpa hasil karena saya biasanya menggunakan kembali pegangan FC yang dibuat dan hasil digunakan sekali kemudian dibuang , lebih mudah bagi saya untuk membaca dan memahami apa fcs.append()yang dilakukan daripada 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            

Hasil:

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

Ini sekarang dalam modul yang saya sebut arcplus *. Tempatkan dengan kode atau PYTHONPATH Anda yang lain dan kemudian:

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

Arcplus juga menambahkan pemfilteran wildcard; untuk hanya memproses kelas fitur yang dimulai dengan "HD_" dalam kumpulan fitur yang berisi "Hydro"

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

. * sekarang di Github, ditingkatkan untuk 10.x. Untuk arcgis 9.3 lihat di sini .

Saya menyadari pertanyaan ini ditandai 9.3, tetapi siapa pun yang mencari jawaban yang sama pada 10.1 dan seterusnya lebih baik menggunakan arcpy.da.Walk . Lebih cepat dan lebih akurat daripada ListDatasets / FeatureClasses / Raster / dll.

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)

Fungsi berjalan bekerja dengan cara yang sama seperti berjalan python . Itu iterates melalui direktori di jalur yang diberikan dan pada setiap iterasi, root mewakili path lengkap dari direktori, dan dirs dan dataset adalah daftar dari subdirektori dan file yang terkandung di dalamnya.

Saat berjalan melalui geodatabase, kumpulan data fitur diperlakukan dengan cara yang sama seperti direktori. Jika Anda hanya ingin mendaftar dataset dan dataset fitur dalam folder root dan tidak membuka dataset fitur untuk melihat konten, Anda dapat melakukan:

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

The ListDatasets Metode adalah apa yang saya pikir Anda cari. FGDB ini memiliki FD di dalamnya yang disebut "Wells" dan memiliki 3 FC di dalamnya. Untuk 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
>>>