ArcGIS / ArcToolBox içindeki çokgen bir haritada dört renk teoremini otomatik olarak nasıl uygulayabilirim?

Ben uygulamak gerekir dört renkler teoremini elle her bölgede koymak için her renk seçmek gerekmez bir şekilde bir çokgen şeklinde. Matematiksel veya programlı olarak yapmak için ArcGIS ve ArcToolBox ile kullanılabilecek herhangi bir uzantı, eklenti, komut dosyası veya veri tabanı olup olmadığını bilmek istiyorum, bu yüzden oluşturmak için geldim her harita ile şimdi kullanabilirsiniz.

Her şeyden önce, tüm cevaplar ve yorumlar için teşekkürler. Ne yazık ki, mevcut araçlar QGIS ve ArcGIS'in en son sürümleriyle tam uyumlu değildi. Bu nedenle @polygeo tarafından belirtilen aracı, @Alexandre'den QGIS eklentisini ve @Jens'ten algoritmanın adını (dört renk haritası) kullanarak kendi çözümümü yaptım.

İlgilenenler için kodum (ArcGIS için ama ikinci bölüm QGIS'de de kullanılabilir).

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)

Algoritmanın sadece 4 rengin kullanılmasını garanti etmediğini unutmayın: çözümün var olduğu kanıtlanmış olsa da, bunu elde etmek için "kaba kuvvet" gereklidir. Benim durumumda, yeterince küçük 7 renk var. Komut dosyası çözüm bulunana kadar ek bir döngü olabilir, ancak yüzlerce harita ve 7 renk için bunu yapmam gerekiyor.

Bir VB6 geliştirici örneği ve bir ArcGIS 9.x coğrafi işleme aracı var, ancak bu ArcGIS Idea'daki yorumlardan 10.0+ sürümünde çalışmıyorlar.

Belki birisi onu taşımakla ilgilenirdi.

Bu ilgili sorunun yorumlarında TopoColour adlı bir QGIS çözümü verilmiştir: Renk çokgenleri, böylece her biri komşularından farklıdır

QGIS kullanıyorsanız, ihtiyacınız olanın bir harita eklentisi renklendirme olduğuna inanıyorum .

Ne yazık ki, eklenti sadece QGIS 1.8 sürümü için kullanılabilir, ancak her zaman kodun nasıl çalıştığını indirebilir ve görebilirsiniz!

Bu @ radouxju'nun cevabının bir fonksiyona uyarlanmasıdır. Giriş özelliği katmanına bir renk alanı ekler ve hesaplar. PolygonNeighbours'un alan adı sonlarına bakılmaksızın çalışmalıdır (farklı kullanıcılar / girişler / arcgis sürümleri (?) İçin farklı görünüyorlar)

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