การระบุฟิลด์การเข้าร่วมในโปรแกรมเป็นแบบ ArcMap?

เป็นไปได้หรือไม่ที่จะระบุเขตข้อมูลการเข้าร่วมทางโปรแกรมที่ใช้กับตารางเข้าร่วมสองชุดข้อมูลเข้าด้วยกันใน ArcMap ขณะนี้ฉันกำลังใช้ ArcGIS 10.0, SP5 และต้องการโซลูชัน ArcPyแต่ฉันจะไม่คัดค้านโซลูชันอื่น ๆ หากไม่มีโซลูชัน ArcPy

วิธีหนึ่งที่ฉันลองใช้คือวนลูปผ่านทุกฟิลด์และมองหา "baseName" ที่ตรงกัน แต่นั่นเป็นเพียง "การศึกษาที่เดาได้" ที่คุณหวังว่าชื่อฟิลด์ในฐานข้อมูลทั้งสองจะเหมือนกัน

สำหรับการแสดงกราฟิกของสิ่งที่ฉันเป็นหลังจากนั้นฉันต้องการระบุ "อินพุตเข้าร่วมฟิลด์" และ "เอาท์พุทเข้าร่วมฟิลด์" ตามที่เห็นในกล่องโต้ตอบ "เพิ่มเข้าร่วม" แต่แน่นอนหลังจากที่จริงแล้ว

นี่เป็นคำถามเกี่ยวกับแท็กที่สามารถตรวจพบ“ เข้าร่วม” โดยทางโปรแกรมหรือไม่ แต่ในกรณีนี้ฉันต้องการขยายฟังก์ชั่นเพื่อระบุ FIELD (s) ที่ใช้ในการเข้าร่วมชุดข้อมูลสองชุด (หรือมากกว่า) เข้าด้วยกัน

นี่คือวิธีการ ArcObjects ตามตัวอย่างนี้เพื่อระบุการรวมทั้งหมดในเลเยอร์และระบุปลายทางและชื่อตารางต้นทางและคีย์หลักและกุญแจต่างประเทศ:

1. รับการอ้างอิงไปยังILayerที่มีหนึ่งหรือมากกว่าเข้าร่วม
2. ส่งILayerให้กับIDisplayTable
3. ร่ายIDisplayTable.DisplayTableทรัพย์สินไปIRelQueryTable
4. ในขณะที่ตารางปัจจุบันคือIRelQueryTable:
   1. ตรวจสอบRelQueryTableของDestinationTableและSourceTableคุณสมบัติ
   2. ตรวจสอบOriginPrimaryKeyและOriginForeignKeyคุณสมบัติของIRelQueryTable.RelationshipClassคุณสมบัติ
   3. ตั้งค่าตารางปัจจุบันRelQueryTableเป็นSourceTableคุณสมบัติปัจจุบัน

สคริปต์ Python นี้ (ใช้comtypesและโมดูลตัวช่วยนี้) จะเข้าร่วมการเชื่อมต่อทั้งหมดตั้งแต่ล่าสุดจนถึงเร็วที่สุดและพิมพ์ชื่อปลายทางและชื่อตารางต้นทางคีย์หลักต้นกำเนิดและ foreign key สำหรับการเข้าร่วมแต่ละครั้ง:

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)

ตัวอย่างเอาต์พุตรับเลเยอร์ซอร์สที่มีสามตัวรวม:

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

สำหรับข้อมูลเพิ่มเติมให้ดูที่ฉันจะเข้าถึง ArcObjects จาก Python ได้อย่างไร

ใส่ข้อมูลทั้งหมดของฟิลด์เป็นสตริง (หลังจากสั่งให้พวกเขา) เปรียบเทียบกับฟังก์ชั่น fuzzycompare และเลือกผู้ที่ให้การจับคู่ที่ดีที่สุดหรือจับคู่เบย์กับความแม่นยำแน่นอน

วิธีนี้คือเมื่อข้อมูลบางอย่างไม่พอดี หากคุณคิดว่าทั้งสองคอลัมน์จะพอดีกันเพียงแค่สั่งซื้อและเปรียบเทียบสำหรับการจับคู่ที่สมบูรณ์แบบกับฟังก์ชันเปรียบเทียบปกติ

ลองสิ่งนี้:

• ใช้เครื่องมือการแปลง XSLTจากชุดเครื่องมือเมทาดาทาเพื่อเขียนไฟล์ข้อมูลเมตา xml / html สำหรับชุดข้อมูลที่เป็นปัญหา

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

• ใช้ตัวแยกวิเคราะห์ HTMLเพื่ออ่านในไฟล์เมตาดาต้าและค้นหาฟิลด์การเข้าร่วมจากประวัติการประมวลผลทางภูมิศาสตร์ของเครื่องมือฟิลด์การเข้าร่วม

• เอาต์พุตตัวอย่างจากเครื่องมือการแปลง XSLT

ชื่อตารางที่เข้าร่วมอยู่ในวัตถุ IFeatureLayer - IFeatureLayerDefinition เป็นสตริง .. ซึ่งฉันคิดว่าอาจมี SQL เข้าร่วมและทำให้ชื่อเขตข้อมูล

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

หรือคุณหมายถึงถ้าคุณไม่สามารถเข้าถึงวัตถุนั้นได้?

เพื่อค้นหาฟิลด์ที่ตรงกันโดยไม่คำนึงถึงชื่อฟิลด์คุณสามารถทำสิ่งนี้:

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)