# 平滑化

library(imager)

# グレースケール画像
im <- grayscale(boats)

# ノイズ画像
im.noise <- im + 0.1 * imnoise(width(im), height(im))

# 平滑化（ボックス、メディアン、ガウシアン）を5回繰り返す
im.box <- im.noise
im.med <- im.noise
im.gauss <- im.noise
for (i in 1 : 5) {
  im.box <- boxblur(im.box, 3)
  im.med <- medianblur(im.med, 3)
  im.gauss <- isoblur(im.gauss, sqrt(0.5))
}

# 同時に3枚並べるレイアウト
layout(t(1 : 3))

plot(im.box, interpolate = FALSE, main = "ボックス・フィルタ")
plot(im.med, interpolate = FALSE, main = "メディアン・フィルタ")
plot(im.gauss, interpolate = FALSE, main = "ガウシアン・フィルタ")

# レイアウトを元に戻す
layout(1)