Saya ingin menghapus piksel terisolasi (atau gumpalan piksel <9) dari gambar raster saya.
library(raster)
# create some raster data
r <- raster(ncols=12, nrows=12)
set.seed(0)
r[] <- round(runif(ncell(r))*0.7 )
rc <- clump(r)
Sebelum menghapus gumpalan piksel <9
Setelah menghapus gumpalan piksel <9 :
Di Erdas, ada alat Saringan untuk melakukan ini, tetapi bagaimana cara menggantinya di R?
Jawaban:
Perhatikan bahwa dari
clump
fungsinya, ID rumpun "4" memiliki 2 sel di sisi kanan yang terhubung dengan sel di sebelah kiri.sumber
Terima kasih kepada @gene dan https://geoscripting-wur.github.io/AdvancedRasterAnalysis/ Sekarang saya dapat menjawab pertanyaan saya (disalin dan dimodifikasi):
perpustakaan (raster)
dan hasil yang diinginkan (dengan satu baris dan satu kolom NA ditambahkan di setiap sisi raster).
sumber
Jawaban tambahan untuk data yang diproyeksikan.
bfastSpatial::areaSieve()
mengambil argumen ambang, sehingga rumpun yang lebih kecil dari ambang dihapus.https://www.rdocumentation.org/packages/bfastSpatial/versions/0.6.2/topics/areaSieve
sumber