Kita harus ingat bahwa data ini adalah sampel dari domain litologis diskrit. Seringkali, batas antara dua domain tersebut tidak dapat diidentifikasi di lapangan dan oleh karena itu tidak sah untuk berharap bahwa banyak lokasi sampel akan terletak tepat di sepanjang batas. Solusi yang benar akan menjadi partisi dari area studi dan setiap poligon di dalam partisi itu dapat (dan seringkali akan) melampaui lokasi sampel yang menentukannya. Kecuali untuk perkiraan kasar, ini mengesampingkan pendekatan apa pun yang menggunakan lokasi sampel sebagai simpul dari poligon yang dihasilkan .
Untuk pekerjaan berkualitas tinggi, metode terbaik adalah menyesuaikan model spasial linier umum untuk proses multinomial. Itu adalah prosedur yang membutuhkan keahlian dan upaya yang cukup besar. Sebagai gantinya, Anda dapat mempertimbangkan untuk memperluas setiap titik sampel ke dalam poligon pengaruhnya (alias poligon Thiessen, poligon Voronoi, atau sel Dirichlet). Membatasi ekspansi ke wilayah darat adalah ide yang bagus; ini bisa dilakukan dengan kotak topeng.
Untuk menggambarkan, pertimbangkan dataset yang jauh lebih kecil ini (dari 14.136 poin) yang mewakili 12 kelas litologi yang dibedakan berdasarkan warna:
Berikut adalah detail dari pusat lobus timur, yang menunjukkan posisi tidak teratur dari titik-titik dan perubahan litologi yang relatif cepat di sana. Menelusuri ini secara manual akan menjadi prosedur yang sulit dan sewenang-wenang:
Saya mencapai perluasan dengan mengubah titik-titik ini menjadi grid (sekitar 800 baris dan 1000 kolom) dan menghitung alokasi Euclidean mereka , menggunakan topeng yang membatasi perhitungan untuk tanah yang tidak gletser. (Skema warna dalam dua gambar berikutnya berbeda dari yang sebelumnya.)
Sebagai perbandingan, berikut adalah peta litologis terperinci dari area yang sama yang digambar dengan skala yang sama dengan simbolisasi yang sama:
Dengan set data yang benar-benar besar atau area studi yang berbelit-belit, mungkin perlu segera membuat ubin dan melakukan prosedur ini secara terpisah pada setiap ubin, merapikan hasilnya menjadi satu raster keluaran jika diinginkan. Agar ini berhasil, ubin harus tumpang tindih sedikit untuk menghindari efek tepi (dan kemudian harus dipangkas secara seragam sebelum mosaik).
Alasan utama untuk pergi ke representasi raster adalah (1) cepat dan mudah untuk dihitung dan (2) solusi berbasis vektor yang akurat akan sulit didapat. Jika Anda mencoba buffer, cembung cangkul, lambung cekung, atau apa pun, Anda akan menemukan bahwa mereka semua saling bersinggungan dan mereka masih meninggalkan celah: dengan kata lain, mereka tidak akan menghasilkan partisi ruang yang konsisten secara topologi ke dalam domain litologi yang berbeda.
Salah satu metode berbasis vektor yang akan bekerja adalah untuk menghitung tessellation Voronoi terbatas dari poin ( metode yang baik mengambil O (n * log (n)) waktu untuk n poin), secara spasial menggabungkan sel-sel Voronoi sesuai dengan atribut litologis terkait mereka poin, dan kemudian pisahkan multi-poligon yang dihasilkan ke dalam komponen yang terhubung (jika Anda mau). Namun, jika semua yang Anda butuhkan adalah keluaran vektor , lebih mudah untuk mengelompokkan hasil raster dan mengubahnya menjadi format vektor.