【数値情報処理b】第2回の答え

課題

gradient_forward.R を書き換えて、画像（人によって異なる）をダウンロードしてグレースケール画像にした後にSobelフィルタで輪郭抽出を行うプログラムを作って提出しなさい。

解答例

report2.R 
  2. # 輪郭抽出（Sobelフィルタ）
  3.  
  4. library(imager)
  5.  
  6. # グレースケール画像
  7. im <- grayscale(load.image("画像のURL"))
  8.  
  9. # x方向のSobelフィルタ
  10. f.Sx <- as.cimg(c(-1, 0, 1, -2, 0, 2, -1, 0, 1), x = 3, y = 3)
  11. # y方向のSobelフィルタ
  12. f.Sy <- as.cimg(c(1, 2, 1, 0, 0, 0, -1, -2, -1), x = 3, y = 3)
  13.  
  14. # 畳み込み演算
  15. im.Sx = correlate(im, f.Sx)
  16. im.Sy = correlate(im, f.Sy)
  17.  
  18. # 同時に4枚並べるレイアウト
  19. layout(t(1 : 4))
  20.  
  21. plot(im, main = "元の画像")
  22. plot(im.Sx, main = "x方向のSobelフィルタ")
  23. plot(im.Sy, main = "y方向のSobelフィルタ")
  24. plot(sqrt(im.Sx ^ 2 + im.Sy ^ 2), main = "輪郭の強さ")
  25.  
  26. # レイアウトを元に戻す
  27. layout(1)

フィルタを 8 で割ってもよい（結果は変わらない）。 

# x方向のSobelフィルタ
f.Sx <- as.cimg(c(-1, 0, 1, -2, 0, 2, -1, 0, 1), x = 3, y = 3) / 8
# y方向のSobelフィルタ
f.Sy <- as.cimg(c(1, 2, 1, 0, 0, 0, -1, -2, -1), x = 3, y = 3) / 8