# 輪郭抽出（ラプラシアン）

library(imager)

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

# ラプラシアンのフィルタ
f.lap <- as.cimg(c(0, 1, 0, 1, -4, 1, 0, 1, 0), x = 3, y = 3)

# 畳み込み演算
im.lap = correlate(im, f.lap)

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

plot(im, main = "元の画像")
plot(im.lap, main = "ラプラシアン")
plot(abs(im.lap), main = "ラプラシアン（絶対値）")

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