Identifikasi secara terprogram Bidang yang Bergabung dalam ArcMap?

Apakah mungkin untuk secara terprogram mengidentifikasi Field Gabung yang sedang digunakan untuk Tabel Gabung dua dataset bersama di ArcMap? Saat ini saya menggunakan ArcGIS 10.0, SP5 dan lebih suka solusi ArcPy , namun saya tidak akan menentang solusi lain, jika solusi ArcPy tidak tersedia.

Salah satu metode yang saya coba adalah mengulang-ulang semua bidang dan mencari "baseName" yang cocok, tapi itu hanya "tebakan terpelajar" di mana Anda berharap bahwa nama-nama bidang di kedua basis data sama.

Untuk representasi grafis dari apa yang saya cari, pada dasarnya saya ingin mengidentifikasi "Input Join Field" dan "Output Join Field" seperti yang terlihat dalam dialog "Add Join", tapi tentu saja setelahnya, tentu saja.

Ini adalah pertanyaan tag-on untuk Bisakah “Bergabung” dideteksi secara terprogram? , tetapi dalam kasus ini saya ingin memperluas fungsionalitas untuk mengidentifikasi FIELD (s) yang digunakan untuk menggabungkan dua (atau lebih) dataset bersama.

Berikut ini adalah pendekatan ArcObjects, yang didasarkan pada contoh ini , untuk menghitung semua gabungan pada layer dan daftar tujuan dan nama tabel sumber dan kunci utama dan asing:

1. Dapatkan referensi ke ILayeryang memiliki satu atau lebih gabungan
2. Lemparkan ILayerkeIDisplayTable
3. Cast IDisplayTable.DisplayTableproperti keIRelQueryTable
4. Sementara tabel saat ini adalah IRelQueryTable:
   1. Memeriksa RelQueryTable's DestinationTabledan SourceTablesifat
   2. Periksa OriginPrimaryKeydan OriginForeignKeyproperti IRelQueryTable.RelationshipClassproperti.
   3. Mengatur meja saat ini untuk saat ini RelQueryTable's SourceTableproperti

Skrip Python ini (menggunakan comtypes dan modul helper ini ) akan melewati semua gabungan , dari yang terbaru hingga yang paling awal, dan mencetak nama tabel tujuan dan sumber, kunci primer asal dan kunci asing asal untuk setiap bergabung:

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)

Contoh output, diberi layer sumber dengan tiga gabungan:

join_table_3 master_fc_join_table_1_join_table_2 JOIN_ID_3 master_fc.MASTER_ID
join_table_2 master_fc_join_table_1 JOIN_ID_2 master_fc.MASTER_ID
join_table_1 master_fc JOIN_ID_1 MASTER_ID

Untuk info lebih lanjut, lihat Bagaimana cara mengakses ArcObjects dari Python?

Masukkan semua data bidang dalam string, (setelah dipesan) membandingkannya dengan fungsi fuzzcompare dan memilih orang-orang yang memberikan kecocokan terbaik atau kecocokan dengan ketepatan tertentu.

Solusi ini adalah ketika beberapa data tidak cocok. Jika menurut Anda kedua kolom akan selalu pas, maka pesan saja dan bandingkan untuk pasangan yang cocok dengan fungsi perbandingan biasa.

Coba ini:

• Gunakan alat Transformasi XSLT dari perangkat Metadata untuk menulis file metadata xml / html untuk dataset yang dimaksud.

  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')

• Gunakan parser HTML untuk membaca dalam file metadata dan mencari bidang gabungan dari riwayat geoproses alat Gabung Bidang ini

• Sampel keluaran dari alat Transformasi XSLT

Nama tabel yang tergabung ada di objek IFeatureLayer - IFeatureLayerDefinition sebagai string .. yang menurut saya mungkin berisi SQL join dan dengan demikian nama bidang.

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

Atau maksud Anda jika Anda tidak dapat mengakses objek itu?

untuk menemukan bidang yang cocok terlepas dari nama bidang, Anda dapat melakukan sesuatu seperti ini:

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)