Identificando programaticamente campo unido no ArcMap?

É possível identificar programaticamente o campo Join que está sendo usado para unir dois conjuntos de dados juntos no ArcMap? Atualmente, estou usando o ArcGIS 10.0, SP5 e preferiria uma solução ArcPy , no entanto, não me oporia a outras soluções, se uma solução ArcPy não estiver disponível.

Um método que tentei foi percorrer todos os campos e procurar um "baseName" correspondente, mas esse é apenas um "palpite", no qual você espera que os nomes de campo nos dois bancos de dados sejam os mesmos.

Para uma representação gráfica do que estou procurando, basicamente quero identificar o "Campo de junção de entrada" e "Campo de junção de saída", como visto na caixa de diálogo "Adicionar junção", mas depois do fato, é claro.

Esta é uma pergunta de tag para "Uma junção" pode ser detectada programaticamente? , mas neste caso, desejo estender a funcionalidade para identificar os FIELDs que estão sendo usados ​​para unir os dois (ou mais) conjuntos de dados.

Aqui está uma abordagem do ArcObjects, baseada neste exemplo , para enumerar todas as junções em uma camada e listar seus nomes de tabela de destino e origem e chaves primárias e estrangeiras:

1. Obter uma referência a uma ILayerque tenha uma ou mais junções
2. Transmitir ILayerparaIDisplayTable
3. Transmitir a IDisplayTable.DisplayTablepropriedade paraIRelQueryTable
4. Enquanto a tabela atual é uma IRelQueryTable:
   1. Inspeccionar os RelQueryTable's DestinationTablee SourceTablepropriedades
   2. Inspecione as propriedades OriginPrimaryKeye OriginForeignKeyda IRelQueryTable.RelationshipClasspropriedade.
   3. Defina a tabela atual para o atual RelQueryTable's SourceTablepropriedade

Esse script Python (usando comtypes e este módulo auxiliar ) percorrerá todas as junções, do mais recente ao mais antigo, e imprimirá os nomes das tabelas de destino e de origem, chave primária de origem e chave estrangeira de origem para cada associação:

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)

Exemplo de saída, dada uma camada de origem com três junções:

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 mais informações, consulte Como acesso o ArcObjects a partir do Python?

Coloque todos os dados dos campos em seqüências de caracteres (depois de ordená-los) compare-os com uma função de comparação difusa e selecione aqueles que deram a melhor correspondência ou correspondência além de uma certa precisão.

Essa solução é quando alguns dados não cabem. Se você acha que as duas colunas sempre se encaixam, basta pedir e comparar para uma combinação perfeita com uma função de comparação comum.

Tente o seguinte:

• Use a ferramenta Transformação XSLT do conjunto de ferramentas Metadados para gravar um arquivo de metadados xml / html para o conjunto de dados em questão.

  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 um analisador de HTML para ler o arquivo de metadados e procurar o campo de junção no histórico de geoprocessamento da ferramenta Join Field

• Saída de amostra da ferramenta Transformação XSLT

Os nomes das tabelas unidas estão no objeto IFeatureLayer - IFeatureLayerDefinition como uma seqüência de caracteres .. que eu acho que provavelmente contém a junção SQL e, portanto, os nomes dos campos.

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

Ou você quer dizer se você não pode acessar esse objeto?

para encontrar campos correspondentes, independentemente do nome do campo, você pode fazer algo assim:

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)