Определите минимальное количество полигонов из шейп-файла, чтобы охватить область интереса

У меня есть большое количество шейп-файлов, представляющих интересующие области для анализа, который будет проводиться с использованием различных источников спутниковых изображений (IKONOS, RapidEye и т. Д.). К сожалению, на снимках не используется система траекторий, как, например, Landsat, поэтому экстенты сильно различаются.

У меня есть шейп-файлы, прикрепленные к каждому AOI, представляющие пределы различных приобретений изображений, все из которых уже были сочтены приемлемыми. Некоторые из этих шейп-файлов имеют 500 или более полигонов.

Мне нужно найти подход, предпочтительно тот, который можно автоматизировать (предпочтительно Python и ArcInfo 10, а также FOSS), чтобы определить наименьшее количество полигонов, охватывающих каждую из моих областей интересов.

Как отметил Уубер, обобщение проблемы такого типа для поиска высококачественного решения было бы непросто, но такой подход может приблизить вас к цели без большой работы. Вот некоторый псевдокод, основанный на следующих предположениях:

1. Область интересов A

2. Набор полигонов Y, которые полностью покрывают 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

Идея состоит в том, что вы итеративно уменьшаете свою область интересов с помощью экстента спутника, который имеет наибольшее перекрытие с оставшейся областью интересов. AOI уменьшается на каждой итерации, пока ничего не останется. Это, вероятно, не будет оптимальным решением, но оно должно быть разумным и должно выполняться довольно быстро.

Итак, у вас есть область A, представляющая некоторую область, и набор экстентов изображения, которые можно определить как набор Y.

Если я правильно понял, вы можете выполнять множество различных функций:

1. обрезка экстерьера изображения областью A
2. Выполните выбор по местоположению, используя многоугольники Экстента и опцию полностью содержит

Затем вы можете изучить области каждого из них и определить, есть ли у вас собственные выборки полигонов, выполнив некоторую пространственную сортировку геометрии с использованием ArcPy и курсоров.

Надеюсь, это поможет.