Identification systématique du champ joint dans ArcMap?

Est-il possible d'identifier par programme le champ de jointure utilisé pour regrouper en table deux ensembles de données dans ArcMap? J'utilise actuellement ArcGIS 10.0, SP5 et préférerais une solution ArcPy , mais je ne serais pas opposé à d'autres solutions, si une solution ArcPy n'est pas disponible.

Une méthode que j'ai essayée consistait à parcourir tous les champs et à rechercher un "baseName" correspondant, mais ce n'est qu'une "supposition éclairée" où vous espérez que les noms de champs dans les deux bases de données sont les mêmes.

Pour une représentation graphique de ce que je recherche, je veux essentiellement identifier le "Champ de jointure d'entrée" et le "Champ de jointure de sortie" comme on le voit dans la boîte de dialogue "Ajouter une jointure", mais après coup, bien sûr.

Il s'agit d'une question de balise à Est-ce qu'une «jointure» peut être détectée par programme? , mais dans ce cas, je souhaite étendre la fonctionnalité pour identifier le (s) CHAMP (s) utilisé (s) pour joindre les deux (ou plus) ensembles de données.

Voici une approche ArcObjects, basée sur cet exemple , pour énumérer toutes les jointures sur une couche et répertorier leurs noms de table de destination et source et leurs clés primaire et étrangère:

1. Obtenir une référence à un ILayerqui a une ou plusieurs jointures
2. Cast the ILayertoIDisplayTable
3. Cast la IDisplayTable.DisplayTablepropriété àIRelQueryTable
4. Alors que le tableau actuel est un IRelQueryTable:
   1. Vérifier que les RelQueryTable« s DestinationTableet SourceTablepropriétés
   2. Inspectez les propriétés OriginPrimaryKeyet OriginForeignKeyde la IRelQueryTable.RelationshipClasspropriété.
   3. Mettre la table en cours au courant RelQueryTablede » SourceTablepropriété

Ce script Python (à l'aide de comtypes et de ce module d'assistance ) passera par toutes les jointures, de la plus récente à la plus ancienne, et imprimera les noms des tables de destination et source, la clé primaire d'origine et la clé étrangère d'origine pour chaque jointure:

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)

Exemple de sortie, étant donné une couche source avec trois jointures:

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

Pour plus d'informations, voir Comment accéder à ArcObjects à partir de Python?

Mettez toutes les données des champs dans des chaînes (après les avoir commandées), comparez-les avec une fonction fuzzycompare et sélectionnez ceux qui ont donné la meilleure correspondance ou la correspondance au-delà d'une certaine précision.

Cette solution est lorsque certaines données ne correspondent pas. Si vous pensez que les deux colonnes conviendraient toujours, alors commandez et comparez pour une correspondance parfaite avec une fonction de comparaison ordinaire.

Essaye ça:

• Utilisez l' outil de transformation XSLT du jeu d' outils Métadonnées pour écrire un fichier de métadonnées xml / html pour le jeu de données en question.

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

• Utilisez un analyseur HTML pour lire le fichier de métadonnées et recherchez le champ de jointure à partir de l'historique de géotraitement de l'outil Join Field

• Exemple de sortie de l'outil de transformation XSLT

Les noms des tables jointes sont dans l'objet IFeatureLayer - IFeatureLayerDefinition sous la forme d'une chaîne .. qui, je pense, contient probablement la jointure SQL et donc les noms de champ.

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

Ou voulez-vous dire si vous ne pouvez pas accéder à cet objet?

pour trouver des champs correspondants quel que soit le nom du champ, vous pouvez faire quelque chose comme ceci:

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)