# スタック

library(imager)

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

# 9枚のノイズ画像を加算
n <- 9
for (i in 1 : n) {
  # ノイズ画像の読み込み
  im.noise <- load.image(sprintf("https://www.kitp.org/numericalipb/noise%d.png", i))

  # 画像の加算
  iml <- imlist(im.stack, im.noise)
  im.stack <- add(iml)
}
im.stack <- im.stack / n

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

plot(grayscale(boats), interpolate = FALSE, main = "元の画像")
plot(im.noise, interpolate = FALSE, main = "最後のノイズ画像")
plot(im.stack, interpolate = FALSE, main = "スタック")

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