ArcGIS / ArcToolBoxのポリゴンマップに4色の定理を自動的に適用する方法は？

各色を手動で選択して各領域に配置する必要がないように、4色の定理を多角形に適用する必要があります。ArcGISおよびArcToolBoxで使用して数学またはプログラムで実行できる拡張機能、プラグイン、スクリプト、またはデータベースがあるかどうかを知りたいので、今後作成するすべてのマップで使用できます。

まず第一に、すべての回答とコメントをありがとう。残念ながら、既存のツールはQGISおよびArcGISの最新バージョンと完全に互換性がありませんでした。したがって、@ polygeoで示されるツール、@ AlexandreのQGISプラグイン、および@Jensのアルゴリズム名（4色マップ）を使用して、独自のソリューションを作成しました。

興味のある人向けのコードを次に示します（ArcGIS向けですが、QGISでも2番目の部分を使用できます）。

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)

アルゴリズムは4色のみが使用されることを保証しないことに注意してください。ソリューションが存在することは証明されていますが、それを達成するには「ブルートフォース」が必要です。私の場合、私は十分に小さい7色を得ました。スクリプトには、解決策が見つかるまで追加のループが含まれる場合がありますが、何百ものマップに対してそれを行う必要があり、7色は問題ありません。

ありVB6開発者のサンプルとのArcGIS 9.xのジオプロセシングツールは、この上のコメントからのArcGISアイデア彼らは10.0+で仕事をしません。

おそらく誰かがそれを移植することに興味があるでしょう。

TopoColourと呼ばれるQGISソリューションは、この関連する質問のコメントに記載されています。ポリゴンを色付けして、それぞれが隣接するものと区別されるようにします

QGISを使用している場合、必要なのはColoring a mapプラグインだと思います。

残念ながら、プラグインはQGIS 1.8バージョンでのみ使用できますが、いつでもダウンロードしてコードの動作を確認できます！

これは、@ radouxjuの答えを関数に適合させたものです。入力フィーチャレイヤーにカラーフィールドを追加して計算します。PolygonNeighborsのフィールド名の末尾に関係なく機能するはずです（異なるユーザー/入力/ arcgisバージョン（？）で異なるようです）

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