¿Programar la identificación del campo unido en ArcMap?

¿Es posible identificar mediante programación el campo de unión que se está utilizando para unir dos conjuntos de datos en ArcMap? Actualmente estoy usando ArcGIS 10.0, SP5 y preferiría una solución ArcPy , sin embargo, no me opondría a otras soluciones, si una solución ArcPy no está disponible.

Un método que probé fue recorrer todos los campos y buscar un "baseName" coincidente, pero eso es solo una "suposición educada" en la que esperas que los nombres de campo en ambas bases de datos sean los mismos.

Para una representación gráfica de lo que busco, básicamente quiero identificar el "Campo de unión de entrada" y el "Campo de unión de salida" como se ve en el cuadro de diálogo "Agregar unión", pero después del hecho, por supuesto.

Esta es una pregunta etiquetada para ¿Se puede detectar un "Join" mediante programación? , pero en este caso quiero ampliar la funcionalidad para identificar los CAMPOS que se utilizan para unir los dos (o más) conjuntos de datos.

Aquí hay un enfoque de ArcObjects, basado en este ejemplo , para enumerar todas las uniones en una capa y enumerar sus nombres de tabla de origen y destino y las claves primarias y externas:

1. Obtenga una referencia a una ILayerque tenga una o más uniones
2. Lanzar el ILayeraIDisplayTable
3. Emitir la IDisplayTable.DisplayTablepropiedad aIRelQueryTable
4. Mientras que la tabla actual es un IRelQueryTable:
   1. Inspeccionar los RelQueryTable's DestinationTabley SourceTablepropiedades
   2. Inspeccionar las propiedades OriginPrimaryKeyy OriginForeignKeyde la IRelQueryTable.RelationshipClasspropiedad.
   3. Establecer la tabla actual de la corriente RelQueryTable's SourceTablepropiedad

Este script de Python (usando comtypes y este módulo auxiliar ) pasará por todas las uniones, de la más reciente a la más temprana, e imprimirá los nombres de la tabla de origen y destino, la clave primaria de origen y la clave externa de origen para cada combinación:

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)

Ejemplo de salida, dada una capa fuente con tres combinaciones:

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

Para obtener más información, consulte ¿Cómo accedo a ArcObjects desde Python?

Ponga todos los datos de los campos en cadenas (después de haberlos ordenado) compárelos con una función de comparación difusa y seleccione aquellos que dieron la mejor coincidencia o coincidencia con cierta precisión.

Esta solución es cuando algunos datos no encajan. Si cree que ambas columnas siempre encajarían, simplemente ordene y compare para una combinación perfecta con una función de comparación ordinaria.

Prueba esto:

• Utilice la herramienta de transformación XSLT del conjunto de herramientas de metadatos para escribir un archivo de metadatos xml / html para el conjunto de datos en cuestión.

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

• Use un analizador HTML para leer el archivo de metadatos y busque el campo de unión desde el historial de geoprocesamiento de la herramienta Campo de unión

• Salida de muestra de la herramienta de transformación XSLT

Los nombres de las tablas unidas están en el objeto IFeatureLayer - IFeatureLayerDefinition como una cadena ... que creo que probablemente contiene el SQL de unión y, por lo tanto, los nombres de campo.

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

¿O quieres decir si no puedes acceder a ese objeto?

para encontrar campos coincidentes independientemente del nombre del campo, puede hacer algo como esto:

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)