ArcGIS Desktopを使用して散在するポイントをマップするために動的ページの数を最小化しますか？

時々、興味のあるポイントを表示するためにマップブックを作成する必要があります。通常のメッシュを使用してページを作成する最初のステップ：

a）単一のポイント（25ページなど）が端にあるページがいくつかあり、b）ページが多すぎるため、このソリューションは好きではありません。

最初の問題はコードを使用して簡単に修正できます。ページ範囲の長方形を関連するポイント範囲の中心に移動します。

私はまだそれが好きではありません、ページ数は同じままなので非常に込み合っています。それらはすべて、レポートの複数のコピーで実際のA3用紙ページになることに注意してください。

そこで、ページ数を減らすコードを作成しました。この例では、45から34です。

これが達成できる最良の結果かどうかはわかりませんが、

すべてのポイントをキャプチャするために指定されたサイズの長方形の数を最小限に抑えるためにポイントをシャッフルするための最良の戦略（疑似コード、パブリケーション、Pythonライブラリ）は何ですか？確かに、誰かがゲーム理論、軍事芸術、漁業でそれを発見しました

これは元の質問の更新です：

これは、必要な実際の範囲とページサイズを示しています。

164ページのうち10ページを表示する拡大ズーム：

サンプルポイントフィーチャクラス

長方形のサイズは、制限内に収まるとすぐに変更できます。

これは答えではありません、興味がある人のために私がPythonソリューションを投稿すると思っただけです：

# ---------------------------------------------------------------------------
# PAGE MAKER
# 
# ---------------------------------------------------------------------------
# Import arcpy module
import arcpy, traceback, os, sys
from arcpy import env

width=650
height=500

try:
    def showPyMessage():
            arcpy.AddMessage(str(time.ctime()) + " - " + message)
    mxd = arcpy.mapping.MapDocument("CURRENT")
    points = arcpy.mapping.ListLayers(mxd,"points")[0]
    pgons = arcpy.mapping.ListLayers(mxd,"pages")[0]

    g=arcpy.Geometry()
    geometryList=arcpy.CopyFeatures_management(points,g)
    geometryList=[p.firstPoint for p in geometryList]
    curT = arcpy.da.InsertCursor(pgons,"SHAPE@")
    while True:
        nPoints=len(geometryList)
        small=[geometryList.pop(0)]
        for p in geometryList:
            small.append(p)
            mPoint=arcpy.Multipoint(arcpy.Array(small))
            ext=mPoint.extent
            cHeight=ext.height
            cWidth=ext.width
            if cHeight>height or cWidth>width:
                small.remove(p)
        mPoint=arcpy.Multipoint(arcpy.Array(small))
        ext=mPoint.extent
        xC=(ext.XMin+ext.XMax)/2
        yC=(ext.YMin+ext.YMax)/2
        LL=arcpy.Point (xC-width/2,yC-height/2)
        UL=arcpy.Point (xC-width/2,yC+height/2)
        UR=arcpy.Point (xC+width/2,yC+height/2)
        LR=arcpy.Point (xC+width/2,yC-height/2)
        pgon=arcpy.Polygon(arcpy.Array([LL,UL,UR,LR]))
        curT.insertRow((pgon,))
        short=filter(lambda x: x not in small,geometryList)
        arcpy.AddMessage('Grabbed %i points, %i to go' %(len(small),len(short)))
        if len(short)==0: break
        geometryList=short[:]
    del mxd
except:
    message = "\n*** PYTHON ERRORS *** "; showPyMessage()
    message = "Python Traceback Info: " + traceback.format_tb(sys.exc_info()[2])[0]; showPyMessage()
    message = "Python Error Info: " +  str(sys.exc_type)+ ": " + str(sys.exc_value) + "\n"; showPyMessage()

最近調査計画に適用しました：

更新：

一部のパターンでは、「迷子」ポイントを最初に処理する方法が適しているようです。私は「凸包」の皮を使ってそれらを識別しました、whuberのアイデア、投稿が見つかりません、申し訳ありません。

これは、最大カバレッジ問題の幾何学的バージョンのように見えます。これは、セットカバレッジ問題と密接に関連しており、これら2つはNP完全です。

したがって、それを解決するには、近似を使用できます。私は次のアルゴリズムを試してみると、それは完全に動作するようです。問題の複雑さのために、私たちは最良の答えを見つけることができません。

1. Foreachポイントは、ランダムな距離でN = 10の長方形を生成します。長方形がポイントをカバーしていることを確認するだけです（各長方形には少なくとも1つのポイントが属しており、各ポイントは少なくとも1つの長方形に属しています）。
2. すべてのポイントがカバーされるまで繰り返します。カバーされていないポイントの最大数をカバーする長方形を取得します。ポイントをカバー済みとしてマークします。

このアルゴリズムの実装は、サークル専用です。http：//jsfiddle.net/nwvao72r/3/