R/stackSkeleton.R

In tkatsuki/dipr: Digital Image Processing with R

#' @export

stackSkeleton <- function(img){
  r <- 1
  row <- nrow(img)
  col <- ncol(img)

  img <- rbind(matrix(0, r, col+2*r),
  cbind(matrix(0, row, r), img, matrix(0, row, r)),
  matrix(0, r, col+2*r))

  stack <- array(c(
    img[2:(row+1), 1:col]*2^7,
    img[3:(row+2), 1:col]*2^6,
    img[3:(row+2), 2:(col+1)]*2^5,
    img[3:(row+2), 3:(col+2)]*2^4,
    img[2:(row+1), 3:(col+2)]*2^3,
    img[1:row, 3:(col+2)]*2^2,
    img[1:row, 2:(col+1)]*2^1,
    img[1:row, 1:col]*2^0),
    c(row, col, 8))

  stack
}