Làm cách nào để áp dụng định lý bốn màu trong bản đồ đa giác trong ArcGIS / ArcToolBox tự động?

Tôi cần áp dụng định lý bốn màu theo hình đa giác theo cách mà tôi không cần phải chọn thủ công từng màu để đặt vào từng vùng. Tôi muốn biết liệu có bất kỳ tiện ích mở rộng, trình cắm, tập lệnh hoặc cơ sở dữ liệu nào có thể được sử dụng với ArcGIS và ArcToolBox để thực hiện theo cách toán học hoặc lập trình hay không, vì vậy tôi có thể sử dụng nó ngay bây giờ với mọi bản đồ tôi tạo.

Trước hết, cảm ơn tất cả các câu trả lời và ý kiến. Thật không may, các công cụ hiện có không tương thích hoàn toàn với các phiên bản mới nhất của QGIS và ArcGIS. Do đó, tôi đã thực hiện giải pháp của riêng mình bằng công cụ được chỉ định bởi @polygeo, plugin QGIS từ @Alexandre và tên của thuật toán (bản đồ bốn màu) từ @Jens.

Đây là mã của tôi cho những người quan tâm (đối với ArcGIS nhưng phần thứ hai cũng có thể được sử dụng trong QGIS).

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)

Lưu ý rằng thuật toán không đảm bảo rằng chỉ có 4 màu được sử dụng: mặc dù đã được chứng minh rằng giải pháp tồn tại, "lực lượng vũ phu" là cần thiết để đạt được nó. Trong trường hợp của tôi, tôi có 7 màu đủ nhỏ. Kịch bản có thể có một vòng lặp bổ sung cho đến khi tìm thấy giải pháp, nhưng tôi cần thực hiện nó cho hàng trăm bản đồ và 7 màu là OK.

Có một mẫu nhà phát triển VB6 và một công cụ xử lý địa lý ArcGIS 9.x nhưng từ các nhận xét về Ý tưởng ArcGIS này, họ không hoạt động ở 10.0+.

Có lẽ ai đó sẽ quan tâm đến việc chuyển nó.

Một giải pháp của QGIS có tên TopoColour được đưa ra trong các nhận xét của câu hỏi liên quan này: Đa giác màu sao cho mỗi đa dạng khác biệt với các nước láng giềng

Nếu bạn đang sử dụng QGIS, tôi tin rằng thứ bạn cần là plugin Tô màu bản đồ .

Thật không may, plugin chỉ có sẵn cho phiên bản QGIS 1.8, nhưng bạn luôn có thể tải xuống và xem cách mã hoạt động!

Đây là bản chuyển thể của câu trả lời của @ radouxju thành một hàm. Nó sẽ thêm một trường màu vào lớp tính năng đầu vào và tính toán. Nó nên hoạt động bất kể kết thúc tên trường của PolygonNeighbor (chúng dường như khác nhau đối với các phiên bản người dùng / đầu vào / arcgis khác nhau (?))

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