Déterminer le nombre minimum de polygones du fichier de formes pour couvrir une zone d'intérêt

J'ai un grand nombre de fichiers de formes représentant des domaines d'intérêt pour une analyse qui sera réalisée à l'aide de différentes sources d'imagerie satellite (IKONOS, RapidEye, etc.). Malheureusement, l'imagerie n'utilise pas de système pathrow comme Landsat par exemple, de sorte que les étendues varient considérablement.

J'ai des fichiers de formes coupés sur chaque AOI représentant l'étendue des différentes acquisitions d'images, qui ont toutes déjà été jugées acceptables. Certains de ces fichiers de formes contiennent 500 polygones ou plus.

Je dois trouver une approche, de préférence une qui peut être automatisée (Python et ArcInfo 10 de préférence, FOSS serait également acceptable) pour déterminer le moins de polygones pour couvrir chacun de mes domaines d'intérêt.

Comme indiqué par whuber généraliser ce type de problème pour trouver une solution de haute qualité serait difficile mais cette approche pourrait vous rapprocher suffisamment sans beaucoup de travail. Voici un pseudo-code basé sur les hypothèses suivantes:

1. Zone d'intérêt A

2. Ensemble de polygones Y couvrant complètement A

   Start loop
    Iterate through Y
      Select the polygon x from Y that has greatest area of intersection with A
   
    Clip A with polygon x
    Remove x from Y 
    If A is null then end program

L'idée est que vous réduisez de manière itérative votre zone d'intérêt avec l'étendue du satellite qui a le plus grand chevauchement avec la zone d'intérêt restante. L'AOI diminue à chaque itération jusqu'à ce qu'il ne reste plus rien. Ce ne sera probablement pas une solution optimale mais elle devrait être raisonnable et fonctionner assez rapidement.

Ok, vous avez donc une zone A qui représente une zone et un tas d'étendues d'imagerie qui peuvent être définies comme l'ensemble Y.

Si j'ai ce correct, vous pouvez faire un tas de fonctions différentes:

1. découpage des étendues d'imagerie par la zone A
2. Effectuez une sélection par emplacement à l'aide des polygones d'étendue et l'option contient complètement

Vous pouvez ensuite examiner les zones de chacun et déterminer si vous avez des polygones à sélection automatique en effectuant un tri de géométrie spatiale à l'aide d'ArcPy et de curseurs.

J'espère que cette aide.