Làm thế nào để tạo một kho lưu trữ GIS?

Văn phòng của tôi sẽ thấy một sự thay đổi lớn trong phần GIS của nó. Phần này đã hoạt động từ những năm 1980 và có một bộ sưu tập dữ liệu GIS khổng lồ (ví dụ: shapefiles, tệp raster, dữ liệu, v.v.) nhưng chưa từng qua bất kỳ kho lưu trữ nào. Bây giờ nó sẽ xảy ra.

Có cách nào tự động để trích xuất tất cả thông tin về dữ liệu GIS (ví dụ: shapefile, bảo hiểm thông tin hồ quang, tệp lớp, * .mxd, gdb, tệp raster, v.v.) từ PC sang tệp Excel không? Thông tin có thể bao gồm ngày tạo, ngày chỉnh sửa lần cuối, tên thư mục hoặc vùng chứa, v.v.

Điều này hoạt động với tôi, sử dụng arcpy.da.Walkchức năng tại ArcGIS 10.1 SP1:

import arcpy, csv, os

workspace = r"c:\GISData"
output = r"C:\temp\test.csv"

with open(output, 'wb') as csvfile:
    csvwriter = csv.writer(csvfile)
    for dirpath, dirnames, filenames in arcpy.da.Walk(workspace):
        for filename in filenames:
            desc = arcpy.Describe(os.path.join(dirpath, filename))
            csvwriter.writerow([desc.catalogPath, desc.name, desc.dataType])

Các csvmô-đun cũng được sử dụng để đơn giản hóa bằng văn bản cho tập tin đầu ra. Excel có thể mở tệp CSV để bạn có thể xem chúng dưới dạng bảng tính.

Xem arcpy.Describechức năng cho các thuộc tính bổ sung mà bạn có thể bao gồm trong đầu ra.

Nếu bạn đặc biệt muốn phân tích thông tin từ siêu dữ liệu thực tế , hãy xem tập lệnh trong câu trả lời này: Tạo bảng chứa tất cả tên tệp (và có thể siêu dữ liệu) trong Cơ sở dữ liệu địa lý tệp

Khi bạn sử dụng Python, bạn phải sử dụng các mô-đun chính xác để làm những gì bạn muốn. Ví dụ, để tìm tất cả các tệp trong một thư mục có phần mở rộng shp, có nhiều giải pháp đơn giản hơn được trình bày mà không bị phá vỡ, điều này thật tệ ... (như giải pháp được đưa ra bởi Nathan W, nhưng có rất nhiều, rất nhiều giải pháp khác, chỉ tìm kiếm trên Internet)

Một số ví dụ với các mô-đun có liên quan:

1) với mô-đun toàn cầu:

shapefiles chỉ:

import glob
import os
os.chdir("mydir")
for files in glob.glob("*.shp"):
    print files

shapefiles và geodatabase:

import glob
types = ('*.shp', '*.gbd') # the tuple of file types
files_grabbed = []
for files in types:
     files_grabbed.extend(glob.glob(files)) #files_grabbed = the list of shp and gbd files

nếu bạn muốn tìm kiếm trong thư mục con:

import glob
for f in glob.iglob("/mydir/*/*.shp"): #search immediate subdirectories 
    print f

2) với os.listdir và hiểu danh sách (trong hai dòng) -> danh sách kết quả

path = 'mydir'
shape_files = [f for f in os.listdir(path) if f.endswith('.shp')]
gdb_files = [f for f in os.listdir(path) if f.endswith('.gdb')]

3) với mô-đun fnmatch:

import fnmatch
for file in os.listdir('path'):
    if fnmatch.fnmatch(file, '*.shp'):
        print file

và nhiều giải pháp khác, đệ quy vv

Cảm ơn tác phẩm nghệ thuật21 và Nathan W đã phản hồi của bạn. Và vâng, mã của Nathen đã làm nên điều kỳ diệu.

import os, arcpy

#create blank text file
with open("C:\\Temp\\GISlayers.txt", "w") as txt:
for root, dirs, files in os.walk("C:\\Temp\\temp"):
    for f in files:
        #look for shapefiles
        if f.endswith('.shp'):
            desc = arcpy.Describe(root + "\\" + f)
            #write info to text file
            txt.write(desc.name + "," + desc.catalogPath + "\n")

        #look for file geodatabases
        if f.endswith('.gdb'):
            desc = arcpy.Describe(root)
            for child in desc.children:
                #write info to text file
                txt.write(child.name + "," + child.path + "\n")

        #look for layer files
        if f.endswith('.lyr'):
            desc = arcpy.Describe(root + "\\" + f)
            #write info to text file
            txt.write(desc.name + "," + desc.catalogPath + "\n")

        #look for img file
        if f.endswith('.img'):
            desc = arcpy.Describe(root + "\\" + f)
            #write info to text file
            txt.write(desc.name + "," + desc.catalogPath + "\n")

Chỉ tên tệp và vị trí. Máy tính tôi sẽ làm việc với rất nhiều tệp bảo hiểm (tệp thông tin hồ quang), nó cũng sẽ hoạt động trên chúng chứ?

Nếu bạn có ArcGIS Desktop 10.0 (hoặc bất kỳ gói dịch vụ nào của nó), tôi nghĩ rằng cách tốt nhất của bạn là viết một tập lệnh python sử dụng os.walk để xem qua thư mục GIS đã xác định và tìm kiếm các phần mở rộng tệp GIS phổ biến như .shp ,. gdb, .mdb, v.v ... và ghi kết quả vào tệp văn bản được phân tách bằng dấu phẩy. Sau đó, bạn có thể đưa tệp văn bản vào excel, xem ví dụ mã bên dưới:

import os, arcpy

#create blank text file
txt = open("C:\\Temp\\GISlayers.txt", "w")

for root, dirs, files in os.walk("C:\\Temp\\temp"):
    for f in files:

        #look for shapefiles
        foundSHP = f.find(".shp")
        if foundSHP >0:
            checkEXT = f[-3:]
                if checkEXT <> "xml":
                    desc = arcpy.Describe(root + "\\" + f)
                    #write info to text file
                    txt.write(desc.name + "," + desc.catalogPath + "\n")

        #look for file geodatabases
        foundGDB = f.find(".gdb")
        if foundGDB >0:
            desc = arcpy.Describe(root)
            for child in desc.children:
            #write info to text file
            txt.write(child.name + "," + child.path + "\n")
        break
txt.close()

Nếu bạn đang sử dụng ArcGIS 10.1 (hoặc phiên bản mới hơn) cho Máy tính để bàn thì có một Câu trả lời khác ở đây sử dụng arcpy.da.Walk không có sẵn ở 10.0 hoặc trước đó.

Nếu bạn muốn tránh lập trình, đây có thể là phương pháp dễ nhất và nhanh nhất.

Có một tiện ích bổ sung cho Excel có tên là ASAP Utility . Có bản dùng thử miễn phí 90 ngày nhưng sau đó, đó là $ 49 USD cho sử dụng kinh doanh. Nó miễn phí cho sinh viên hoặc sử dụng cá nhân. Phần bổ sung thêm rất nhiều chức năng hữu ích. Một trong số đó là tạo một danh sách các tập tin trong cấu trúc thư mục. Nó cũng cung cấp các thuộc tính tập tin. Bạn có thể giới hạn kết quả theo loại tệp nếu bạn muốn.

Đây là một video về cách làm điều này.

Tôi đã sử dụng tiện ích này trước đây và kết quả rất nhanh.

Lưu ý, tôi không liên kết với công ty phần mềm này.

Tôi không thể có được câu trả lời khác để làm việc đầy đủ.

Trong ví dụ đầu tiên, trong một thư mục có cả cơ sở dữ liệu địa lý và shapefiles, tôi chỉ nhận được một danh sách các lớp tính năng trong cơ sở dữ liệu địa lý, nhưng khi tôi nhận xét phần cơ sở dữ liệu địa lý của tập lệnh, tôi đã nhận được một danh sách các shapefile.

Trong ví dụ thứ hai, phần cơ sở dữ liệu địa lý hoàn toàn không hoạt động, vì vậy tôi đã sao chép trong phần cơ sở dữ liệu địa lý của ví dụ đầu tiên. Một lần nữa, tôi nhận được một danh sách chỉ các cơ sở dữ liệu địa lý.

Sau đó, nó đánh tôi: cơ sở dữ liệu địa lý được đọc trước khi shapefiles và tập lệnh dừng lại ở break phần cơ sở dữ liệu địa lý.

Trở thành một người mới chơi python Tôi không biết tại sao lại breakcần thiết, nhưng không có nó, kịch bản dường như đi vào một vòng lặp vô tận, nhưng vì điều đó breaklà cần thiết đối với tôi khi đặt phần cơ sở dữ liệu địa lý vào vòng lặp của chính nó, sau tập tin khác các loại được liệt kê, sẽ giải quyết vấn đề:

#create blank text file with open("C:\\Temp\\GISlayers.txt", "w") as txt: for root, dirs, files in os.walk("C:\\Temp\\temp"): for f in files: #look for shapefiles, etc. [code...] for f in files: #look for geodatabases [code...]

Khi tôi làm điều đó tôi đã nhận được danh sách đầy đủ của tôi.