Lập trình xác định Trường tham gia trong ArcMap?

Có thể xác định theo chương trình Trường tham gia đang được sử dụng để Bảng nối hai bộ dữ liệu với nhau trong ArcMap không? Tôi hiện đang sử dụng ArcGIS 10.0, SP5 và muốn sử dụng giải pháp ArcPy , tuy nhiên tôi sẽ không phản đối các giải pháp khác, nếu không có giải pháp ArcPy.

Một phương pháp tôi đã thử là lặp qua tất cả các trường và tìm kiếm một "tên cơ sở" phù hợp, nhưng đó chỉ là một "phỏng đoán có giáo dục" mà bạn đang hy vọng rằng tên trường trong cả hai cơ sở dữ liệu đều giống nhau.

Để thể hiện đồ họa của những gì tôi đang theo đuổi, về cơ bản, tôi muốn xác định "Trường tham gia đầu vào" và "Trường tham gia đầu ra" như trong hộp thoại "Thêm tham gia", nhưng thực tế là vậy.

Đây là một câu hỏi gắn thẻ để Có thể phát hiện ra một người tham gia chương trình có thể lập trình được không? , nhưng trong trường hợp này tôi muốn mở rộng chức năng để xác định (các) FIELD đang được sử dụng để nối hai (hoặc nhiều) bộ dữ liệu lại với nhau.

Dưới đây là cách tiếp cận ArcObjects, dựa trên ví dụ này , để liệt kê tất cả các phép nối trên một lớp và liệt kê tên bảng đích và bảng nguồn của chúng và các khóa chính và khóa ngoài:

1. Nhận một tham chiếu đến một ILayercó một hoặc nhiều tham gia
2. Cast ILayerđểIDisplayTable
3. Đúc IDisplayTable.DisplayTabletài sản đểIRelQueryTable
4. Trong khi bảng hiện tại là IRelQueryTable:
   1. Kiểm tra việc RelQueryTable's DestinationTablevà SourceTablethuộc tính
   2. Kiểm tra OriginPrimaryKeyvà OriginForeignKeycác tài sản của IRelQueryTable.RelationshipClasstài sản.
   3. Đặt bảng hiện nay với dòng RelQueryTable's SourceTabletài sản

Tập lệnh Python này (sử dụng comtypes và mô đun trình trợ giúp này ) sẽ đi qua tất cả các phép nối, từ mới nhất đến sớm nhất và in tên bảng đích và nguồn, khóa chính gốc và khóa ngoại gốc cho mỗi phép nối:

from ESRICOMHelpers import * # helper module from https://gis.stackexchange.com/a/5082/753
esriArcMapUI = GetESRIModule("esriArcMapUI")
esriCarto = GetESRIModule("esriCarto")
esriGeoDatabase = GetESRIModule("esriGeoDatabase")

def listJoins(table):
    while CType(table, esriGeoDatabase.IRelQueryTable):
        relQueryTable = CType(table, esriGeoDatabase.IRelQueryTable)
        destTable = relQueryTable.DestinationTable
        sourceTable = relQueryTable.SourceTable
        destDataset = CType(destTable, esriGeoDatabase.IDataset)
        sourceDataset = CType(sourceTable, esriGeoDatabase.IDataset)
        relClass = relQueryTable.RelationshipClass
        print destDataset.Name, sourceDataset.Name, relClass.OriginPrimaryKey, relClass.OriginForeignKey
        table = sourceTable

if __name__ == "__main__":
    #app = GetCurrentApp() # Use if run in-process
    app = GetApp("ArcMap") # Use if run in a standalone script
    mxd = CType(app.Document, esriArcMapUI.IMxDocument)

    # Gets the first layer in the active data frame
    map = mxd.FocusMap
    lyr = map.Layer[0]

    # Need to get the "display table" to access the joins
    displayTable = CType(lyr, esriCarto.IDisplayTable).DisplayTable

    # List the layer's joined tables
    listJoins(displayTable)

Ví dụ đầu ra, được cung cấp một lớp nguồn với ba phép nối:

tham gia_table_3 master_fc_join_table_1_join_table_2 THAM GIA_ID_3 master_fc.MASTER_ID
tham gia_table_2 master_fc_join_table_1 THAM GIA_ID_2 master_fc.MASTER_ID
tham gia_table_1 master_fc THAM GIA_ID_1 MASTER_ID

Để biết thêm thông tin, hãy xem Làm cách nào để tôi truy cập ArcObjects từ Python?

Đặt tất cả dữ liệu của các trường trong chuỗi, (sau khi đã đặt hàng chúng) so sánh chúng với hàm fuzzycompare và chọn những người đưa ra kết quả khớp tốt nhất hoặc khớp với độ chính xác nhất định.

Giải pháp này là khi một số dữ liệu không phù hợp. Nếu bạn nghĩ rằng cả hai cột sẽ luôn phù hợp, thì chỉ cần đặt hàng và so sánh cho một kết hợp hoàn hảo với một chức năng so sánh thông thường.

Thử đi:

• Sử dụng công cụ Chuyển đổi XSLT từ bộ công cụ Siêu dữ liệu để viết ra tệp siêu dữ liệu xml / html cho tập dữ liệu được đề cập.

  arcpy.XSLTransform_conversion(r'X:\temp\Scratch.gdb\fc_FeatureToPoint',"C:\Program Files\ArcGIS\Desktop10.1\Metadata\Stylesheets\ArcGIS.xsl", r'X:\temp\Metadata.html')

• Sử dụng trình phân tích cú pháp HTML để đọc trong tệp siêu dữ liệu và tìm kiếm trường tham gia từ lịch sử xử lý địa lý của công cụ Tham gia

• Đầu ra mẫu từ công cụ Chuyển đổi XSLT

Tên bảng đã tham gia nằm trong đối tượng IFeatureLayer - IFeatureLayerDefDef dưới dạng chuỗi .. mà tôi nghĩ có thể chứa SQL tham gia và do đó là tên trường.

http://edndoc.esri.com/arcobjects/8.3/diagrams/Map%20Layer%20Object%20Model.pdf

Hay bạn có nghĩa là nếu bạn không thể truy cập vào đối tượng đó?

để tìm các trường khớp bất kể tên trường, bạn có thể làm một cái gì đó như thế này:

import arcpy

fc = r"temp/RiversJoined.shp"

fldList1 = [f.name for f in arcpy.ListFields(fc)]
fldList2 =[g.name for g in arcpy.ListFields(fc)]

for f in fldList1:
    values1 = [f_row[0] for f_row in arcpy.da.SearchCursor(fc, (f))]
    for g in fldList2:
        values2 = [g_row[0] for g_row in arcpy.da.SearchCursor(fc,(g))]
        #compare field values
        #get names of matching fields
        if (fldList2.index(g) != fldList1.index(f) and values1 == values2):
            print "match: " + str(g) + " match: "+ str(f)