# リサイズ（線形補間法）

library(imager)

im <- grayscale(boats)

# 倍率
ax <- 2
ay <- 2

# 元画像の次元
dim0 <- dim(im)

# 新画像の次元
dim1 <- dim0
# 横幅のリサイズ
dim1[1] <- round(ax * dim0[1])
# 縦幅のリサイズ
dim1[2] <- round(ay * dim0[2])

# cimg変数で処理すると遅いので配列で処理する
g0 <- array(im, dim0)
g1 <- array(0, dim1)

# 新画像の座標でループ
for (x1 in 1 : dim1[1]) {
  # 逆変換 x1→x0
  x0 <- x1 / ax
  xi <- as.integer(x0)
  if (xi < 1 || xi + 1 > dim0[1]) next
  p <- x0 - xi
  for (y1 in 1 : dim1[2]) {
    # 逆変換 y1→y0
    y0 <- y1 / ay
    yi <- as.integer(y0)
    if (yi < 1 || yi + 1 > dim0[2]) next
    q = y0 - yi
    g1[x1, y1, ,] <-
      (1 - p) * (1 - q) * g0[xi, yi, ,] +
      p * (1 - q) * g0[xi + 1, yi, ,] +
      (1 - p) * q * g0[xi, yi + 1, ,] +
      p * q * g0[xi + 1, yi + 1, ,]
  }
}

# 配列→cimg変数
im.resize <- as.cimg(g1)

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

xl <- c(1, width(im.resize))
yl <- c(height(im.resize), 1)
plot(im, xlim = xl, ylim = yl, interpolate = FALSE, main = "元の画像")
plot(im.resize, xlim = xl, ylim = yl, interpolate = FALSE, main = "線形補間法")

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