Masalah ini memiliki banyak solusi yang valid. Salah satunya bekerja sedikit seperti deskripsi Anda, tetapi alih-alih mengiris poligon di lokasi "acak" Anda dapat melakukannya dengan sengaja dengan cara yang dirancang untuk meminimalkan jumlah perhitungan.
Berikut adalah algoritma dasar. Inputnya terdiri dari setiap arah sapuan bidang, P poligon dari area bukan nol, area target a antara nol dan area poligon, dan ambang batas nonnegatif t (dalam unit area). Tujuannya adalah untuk perpecahan P dengan tegak lurus garis ke arah menyapu menjadi dua bagian, satu di sebelah kanan garis dan yang lainnya ke kiri dari garis, sehingga perbedaan antara daerah sebelah kanan dan area target yang ada lebih besar dari t .
Biarkan L menjadi garis berorientasi apa pun yang tegak lurus terhadap arah sapuan. Tentukan f (L) sebagai area P yang ditemukan di sebelah kanan L, minus a . Dalam istilah-istilah ini tugasnya adalah menemukan nol dari f . Karena f tidak mungkin dapat dibedakan, tetapi kontinu, gunakan salah satu metode pembagian dua bagian , metode garis potong , atau - favorit saya- -Metode bos . Semua sederhana dan dijamin akan menyatu. Gunakan t untuk toleransi konvergensi untuk argumen.
Itu dia. Mari kita pertimbangkan apa yang masuk ke pengkodean ini. Temuan root bersifat rutin - Anda dapat menggunakan potongan kode umum untuknya - sehingga pekerjaan GIS diturunkan ke pengkodean f . Untuk itu diperlukan
1. Splitting the polygon by a line.
2. Computing the area of the piece(s) to the right of the line.
Kedua operasi diimplementasikan di hampir semua GIS berbasis vektor. Jika tidak, Anda dapat mengganti garis dengan persegi panjang yang sangat besar yang mewakili setengah bidang di sebelah kanan garis. Langkah 1 menjadi
1'. Clip the polygon to the rectangle.
Itu adalah operasi yang sangat mendasar.
Untuk memulai dengan pencarian root, Anda perlu menemukan interval di mana nol f dijamin terletak. Ini mudah: proyeksikan amplop poligon ("kotak pembatas") ke arah sapuan garis. Proyeksi adalah interval yang Anda inginkan.
Pertanyaan ini memiliki sejarah panjang. Saya menerapkan algoritma ini untuk ArcView 3.x sejak lama dan menggambarkannya berkali-kali di forum pengguna ESRI lama. Google
huber split situs poligon: forums.esri.com
untuk diskusi, tautan ke kode, perangkat tambahan dan variasi (seperti membelah poligon menjadi bagian-bagian dengan ukuran yang diinginkan yang sekompak mungkin) dan algoritma untuk data raster.
Berikut adalah seperti apa bentuk negara bagian AS (dalam proyeksi luas yang sama) dengan sepertiga bagian bawah setiap negara yang dinaungi. Jelas arah sapuannya vertikal.
