Saya memiliki sejumlah besar shapefile yang mewakili bidang yang diminati untuk analisis yang akan dilakukan menggunakan berbagai sumber citra satelit (IKONOS, RapidEye, dll.). Sayangnya, citra tidak menggunakan sistem jalur seperti Landsat misalnya, sehingga luasannya sangat bervariasi.
Saya memiliki shapefile yang terpotong untuk setiap AOI yang mewakili luasan akuisisi citra yang berbeda, yang semuanya sudah dianggap dapat diterima. Beberapa shapefile ini memiliki 500 poligon atau lebih.
Saya perlu menemukan pendekatan, lebih disukai yang dapat diotomatiskan (Python dan ArcInfo 10 lebih disukai, FOSS juga dapat diterima) untuk menentukan jumlah poligon paling sedikit untuk mencakup setiap bidang yang saya minati.
arcgis-10.0
python
shapefile
geoprocessing
Chad Hawkins
sumber
sumber
Jawaban:
Seperti dicatat oleh whuber generalisasi jenis masalah ini untuk menemukan solusi berkualitas tinggi akan sulit tetapi pendekatan ini mungkin membuat Anda cukup dekat tanpa banyak kerja. Berikut adalah beberapa kode pseudo berdasarkan asumsi berikut:
Set poligon Y yang sepenuhnya menutupi A
Idenya adalah Anda secara iteratif mengurangi Area Minat Anda dengan tingkat satelit yang memiliki tumpang tindih terbesar dengan Area Minat yang tersisa. AOI semakin kecil pada setiap iterasi sampai tidak ada yang tersisa. Ini mungkin bukan solusi yang optimal tetapi harus masuk akal dan harus berjalan cukup cepat.
sumber
Ok, jadi Anda memiliki area A yang mewakili beberapa area dan sekelompok luasan pencitraan yang dapat didefinisikan sebagai set Y.
Jika ini benar, Anda dapat melakukan banyak fungsi berbeda:
Anda kemudian dapat memeriksa area masing-masing dan menentukan apakah Anda memiliki poligon memilih sendiri dengan melakukan beberapa penyortiran geometri spasial menggunakan ArcPy dan kursor.
Semoga bantuan ini.
sumber