Programowo identyfikuje połączone pole w ArcMap?

Czy jest możliwe programowe zidentyfikowanie pola łączenia, które jest używane do łączenia tabeli dwóch zestawów danych w ArcMap? Obecnie używam ArcGIS 10.0, SP5 i wolę rozwiązanie ArcPy , jednak nie sprzeciwiałbym się innym rozwiązaniom, jeśli rozwiązanie ArcPy nie jest dostępne.

Jedną z metod, które wypróbowałem, było zapętlenie wszystkich pól i poszukiwanie pasującego „baseName”, ale jest to tylko „wyuczone przypuszczenie”, w którym masz nadzieję, że nazwy pól w obu bazach danych są takie same.

Aby przedstawić graficznie to, czego szukam, zasadniczo chcę zidentyfikować „wejściowe pole łączenia” i „wyjściowe pole łączenia”, jak widać w oknie dialogowym „Dodaj łączenie”, ale oczywiście po fakcie.

To jest pytanie do tagu. Czy program „Dołącz” może zostać wykryty programowo? , ale w tym przypadku chcę rozszerzyć funkcjonalność w celu identyfikacji FIELD (ów) używanych do łączenia dwóch (lub więcej) zestawów danych razem.

Oto podejście ArcObjects, oparte na tym przykładzie , w celu wyliczenia wszystkich złączeń na warstwie i wyświetlenia ich nazw miejsc docelowych i tabel źródłowych oraz kluczy podstawowych i obcych:

1. Uzyskaj odniesienie do elementu, ILayerktóry ma jedno lub więcej złączeń
2. Prześlij ILayerdoIDisplayTable
3. Prześlij IDisplayTable.DisplayTablewłaściwość doIRelQueryTable
4. Podczas gdy bieżąca tabela to IRelQueryTable:
   1. Sprawdzić RelQueryTable„s DestinationTablei SourceTablewłaściwości
   2. Sprawdź OriginPrimaryKeyi OriginForeignKeywłaściwości IRelQueryTable.RelationshipClassnieruchomości.
   3. Ustawić aktualną tabelę do obecnego RelQueryTable„s SourceTablenieruchomości

Ten skrypt w języku Python (wykorzystujący komty i ten moduł pomocniczy ) przejdzie przez wszystkie złączenia, od najnowszego do najwcześniejszego, i wydrukuje nazwy tabeli docelowej i źródłowej, pierwotny klucz źródłowy i początkowy klucz obcy dla każdego złączenia:

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)

Przykładowe dane wyjściowe, biorąc pod uwagę warstwę źródłową z trzema połączeniami:

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

Aby uzyskać więcej informacji, zobacz Jak uzyskać dostęp do ArcObjects z Python?

Umieść wszystkie dane pól w ciągach, (po ich zamówieniu) porównaj je z funkcją fuzzycompare i wybierz te, które dały najlepsze dopasowanie lub dopasowanie, z określoną precyzją.

To rozwiązanie występuje, gdy niektóre dane nie pasują. Jeśli uważasz, że obie kolumny zawsze pasują, po prostu zamów i porównaj, aby uzyskać idealne dopasowanie za pomocą zwykłej funkcji porównania.

Spróbuj tego:

• Użyj narzędzia Przekształcenie XSLT z zestawu narzędzi Metadane, aby zapisać plik metadanych xml / html dla danego zestawu danych.

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

• Użyj parsera HTML, aby odczytać plik metadanych i wyszukać pole złączenia z historii geoprzetwarzania narzędzia Połącz pole

• Przykładowe dane wyjściowe z narzędzia Transformacja XSLT

Połączone nazwy tabel znajdują się w obiekcie IFeatureLayer - IFeatureLayerDefinition jako ciąg znaków, który, jak sądzę, prawdopodobnie zawiera sprzężony SQL, a tym samym nazwy pól.

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

Czy masz na myśli, że nie możesz uzyskać dostępu do tego obiektu?

aby znaleźć pasujące pola niezależnie od nazwy pola, możesz zrobić coś takiego:

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)