# 膨張・収縮

library(imager)

# グレースケール画像
im <- grayscale(load.image("https://www.kitp.org/numericalipb/nameplate.jpg"))

# しきい値
thr <- 0.3

# 0で埋めたcimg変数
im.bin <- imfill(dim = dim(im))

# imの画素値がしきい値以上なら画素値1とする
im.bin[im >= thr] <- 1

# 膨張
im.bin1 <- dilate_square(im.bin, 3)

# 収縮
im.bin2 <- erode_square(im.bin1, 3) 

layout(t(matrix(1 : 4, 2, 2)))
plot(im, interpolate = FALSE, main = "元の画像")
plot(im.bin, interpolate = FALSE, main = "2値化")
plot(im.bin1, interpolate = FALSE, main = "膨張")
plot(im.bin2, interpolate = FALSE, main = "収縮")
layout(1)