Wie wende ich den Vier-Farben-Satz in einer Polygon-Map in ArcGIS / ArcToolBox automatisch an?

Ich muss das Vier-Farben-Theorem in einer polygonalen Form so anwenden , dass ich nicht jede Farbe manuell auswählen muss, um sie in die einzelnen Bereiche einzufügen. Ich möchte wissen, ob es Erweiterungen, Plug-Ins, Skripte oder Datenbanken gibt, die mit ArcGIS und ArcToolBox mathematisch oder programmatisch verwendet werden können, sodass ich sie jetzt für jede Karte verwenden kann, die ich erstellen möchte.

Zunächst einmal vielen Dank für alle Antworten und Kommentare. Leider waren die vorhandenen Tools nicht vollständig mit den neuesten Versionen von QGIS und ArcGIS kompatibel. Aus diesem Grund habe ich meine eigene Lösung mit dem von @polygeo angegebenen Tool, dem QGIS-Plugin von @Alexandre und dem Namen des Algorithmus (vierfarbige Karte) von @Jens erstellt.

Hier ist mein Code für Interessenten (für ArcGIS, aber der zweite Teil könnte auch in QGIS verwendet werden).

arcpy.MakeFeatureLayer_management(fc, fc[:-4]+ "_lyr" )
try:
    arcpy.AddField_management(fc[:-4] + "_lyr", "color", "SHORT")
except:
    print "field alread exists"   
arcpy.CalculateField_management(fc[:-4] + "_lyr", "color",  "10" , "PYTHON")

arcpy.PolygonNeighbors_analysis(fc[:-4] + "_lyr", fc[:-4] + "_tb.dbf" )
graph = []
cursor=arcpy.da.SearchCursor( fc[:-4] + "_tb.dbf" , ("src_FID","nbr_FID") )
for row in cursor:
    graph.append(row)


pols = arcpy.da.UpdateCursor(fc[:-4] + "_lyr", ("OID@","color"))
colored = []
for pol in pols:
    nbrs = [ second for first, second in graph if first == pol[0]]
    usedcolors = []
    for nbr in nbrs:
        usedcolors += [second for first, second in colored if first == nbr]
    pol[1]=[color for color in range(10) if color not in usedcolors][0]
    colored.append(pol)
    pols.updateRow(pol)

Beachten Sie, dass der Algorithmus nicht garantiert, dass nur 4 Farben verwendet werden: obwohl nachgewiesen wurde, dass die Lösung existiert, ist die "rohe Kraft" erforderlich, um dies zu erreichen. In meinem Fall habe ich 7 Farben, was klein genug ist. Das Skript könnte eine zusätzliche Schleife haben, bis die Lösung gefunden ist, aber ich muss es für Hunderte von Karten und 7 Farben tun, ist OK.

Es gibt ein VB6-Entwicklerbeispiel und ein ArcGIS 9.x-Geoverarbeitungswerkzeug, aber aus den Kommentaren zu dieser ArcGIS-Idee geht hervor, dass sie mit 10.0+ nicht funktionieren.

Vielleicht wäre jemand daran interessiert, es zu portieren.

In den Kommentaren zu dieser verwandten Frage wird eine QGIS-Lösung namens TopoColour angegeben : Färben Sie Polygone so, dass sich jedes von seinen Nachbarn unterscheidet

Wenn Sie QGIS verwenden, brauchen Sie meines Erachtens das Coloring a map-Plugin .

Leider ist das Plugin nur für die QGIS 1.8-Version verfügbar, aber Sie können es jederzeit herunterladen und sehen, wie der Code funktioniert!

Dies ist eine Anpassung von @ radouxjus Antwort in eine Funktion. Es wird ein Farbfeld zum Eingabe-Feature-Layer hinzugefügt und berechnet. Es sollte unabhängig von den Feldnamenendungen von PolygonNeighbors funktionieren (sie scheinen für verschiedene Benutzer / Eingaben / arcgis-Versionen unterschiedlich zu sein (?))

def color_me(feature_layer):
    import arcpy
    try:
        arcpy.AddField_management(feature_layer, 'color', 'SHORT')
    except:
        print 'field alread exists'   

    arcpy.CalculateField_management(feature_layer, 'color',  '10' , 'PYTHON')

    arcpy.PolygonNeighbors_analysis(feature_layer, r'in_memory\neighbor_table' )
    graph = []
    neighbor_fields = [f.name for f in arcpy.ListFields(r'in_memory\neighbor_table') if f.name.startswith(('src', 'nbr'))]
    cursor=arcpy.da.SearchCursor(r'in_memory\neighbor_table' , neighbor_fields)
    for row in cursor:
        graph.append(row)

    pols = arcpy.da.UpdateCursor(feature_layer, ('OID@','color'))
    colored = []

    for pol in pols:
        nbrs = [ second for first, second in graph if first == pol[0]]
        usedcolors = []
        for nbr in nbrs:
            usedcolors += [second for first, second in colored if first == nbr]
        pol[1]=[color for color in range(10) if color not in usedcolors][0]
        colored.append(pol)
        pols.updateRow(pol)
    arcpy.Delete_management(r'in_memory\neighbor_table')