กำหนดจำนวนรูปหลายเหลี่ยมขั้นต่ำจาก shapefile เพื่อครอบคลุมพื้นที่ที่สนใจ

ฉันมีไฟล์รูปร่างเป็นจำนวนมากที่แสดงถึงพื้นที่ที่น่าสนใจสำหรับการวิเคราะห์ที่จะดำเนินการโดยใช้แหล่งที่มาของภาพถ่ายดาวเทียม (IKONOS, RapidEye ฯลฯ ) น่าเสียดายที่ภาพไม่ได้ใช้ระบบ Pathrow เช่น Landsat เป็นต้นขอบเขตจึงแตกต่างกันอย่างมาก

ฉันมีรูปร่างที่ตัดกับ AOI แต่ละตัวซึ่งแสดงถึงขอบเขตของการได้รับภาพที่แตกต่างกันซึ่งทั้งหมดนี้ได้รับการยอมรับแล้ว บางรูปร่างเหล่านี้มีรูปหลายเหลี่ยมมากกว่า 500 รูป

ฉันต้องการค้นหาวิธีที่ควรเป็นแบบอัตโนมัติ (Python และ ArcInfo 10 โดยเฉพาะอย่างยิ่ง FOSS ก็เป็นที่ยอมรับเช่นกัน) เพื่อกำหนดรูปหลายเหลี่ยมจำนวนน้อยที่สุดเพื่อครอบคลุมพื้นที่ที่ฉันสนใจแต่ละแห่ง

ตามที่ระบุไว้โดย whuber generalizing ปัญหาประเภทนี้เพื่อค้นหาวิธีแก้ไขปัญหาคุณภาพสูงจะยุ่งยาก แต่วิธีนี้อาจทำให้คุณปิดพอโดยไม่ต้องทำงานมาก นี่คือรหัสหลอกบางส่วนตามสมมติฐานต่อไปนี้:

1. พื้นที่ที่น่าสนใจ

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. ทำการเลือกตามสถานที่โดยใช้รูปหลายเหลี่ยม Extent และตัวเลือกบรรจุอย่างสมบูรณ์

จากนั้นคุณสามารถตรวจสอบพื้นที่ของแต่ละพื้นที่และพิจารณาว่าคุณมีรูปหลายเหลี่ยมที่เลือกด้วยตนเองหรือไม่โดยทำการเรียงลำดับเชิงพื้นที่โดยใช้ ArcPy และเคอร์เซอร์

หวังว่าความช่วยเหลือนี้