Determina il numero minimo di poligoni dal file di forma per coprire un'area di interesse

Ho un gran numero di shapefile che rappresentano le aree di interesse per un'analisi che verrà condotta utilizzando varie fonti di immagini satellitari (IKONOS, RapidEye, ecc.). Sfortunatamente, le immagini non usano un sistema di pathrow come Landsat per esempio, quindi le estensioni variano notevolmente.

Ho gli shapefile ritagliati in ogni AOI che rappresentano l'estensione di diverse acquisizioni di immagini, che sono già state ritenute accettabili. Alcuni di questi shapefile hanno 500 o più poligoni.

Devo trovare un approccio, preferibilmente uno che può essere automatizzato (preferibilmente Python e ArcInfo 10, anche FOSS sarebbe accettabile) per determinare il minor numero di poligoni per coprire ciascuna delle mie aree di interesse.

Come notato da un whuber generalizzare questo tipo di problema per trovare una soluzione di alta qualità sarebbe difficile ma questo approccio potrebbe avvicinarti abbastanza senza molto lavoro. Ecco alcuni pseudo codici basati sui seguenti presupposti:

1. Area di interesse A

2. Set di poligoni Y che coprono completamente 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'idea è che stai riducendo in modo iterativo la tua area di interesse con l'estensione satellitare che ha la più grande sovrapposizione con la restante area di interesse. L'AOI si riduce su ogni iterazione fino a quando non viene lasciato nulla. Probabilmente questa non sarà una soluzione ottimale, ma dovrebbe essere ragionevole e dovrebbe funzionare abbastanza rapidamente.

Ok, quindi hai un'area A che rappresenta un'area e un mucchio di estensioni di immagini che possono essere definite come set Y.

Se ho questo corretto, puoi fare un sacco di diverse funzioni:

1. ritagliare le estensioni di immagini dall'area A
2. Esegui una selezione per posizione utilizzando i poligoni Extent e l'opzione contiene completamente

È quindi possibile esaminare le aree di ciascuno e determinare se si dispone di poligoni a selezione automatica eseguendo un ordinamento della geometria spaziale utilizzando ArcPy e i cursori.

Spero che questo aiuto.