R/deform_elements.R

In jakubnowicki/artificialmorpho: Generator of artificial landmark sets

#' Deform object divided into elements
#'
#' Deforms object with randomized strain matrix, randomization is common for landmarks in the same element.
#'
#' @param data data matrix
#' @param elements.list list of elements (landmark numbers)
#' @param a deformation amount
#' @param theta deformation angle
#' @param random.sd deformation random sd
#' @export
#' @import morphoutils

deform.elements <- function(data,elements.list,a,theta,random.sd) {
  splited <- split.matrix(data,elements.list)
  n.elements <- length(splited)
  tmp.list <- list()
  for (i in 1:n.elements) {
    strain <- strain.matrix(a = a+rnorm(1,0,random.sd),theta = theta + rnorm(1,0,random.sd))
    tmp <- deformacja(data = splited[[i]],strain.matrix = strain)
    tmp.list[[i]]<-tmp
  }
  output <- do.call(rbind,tmp.list)
  output <- output[order(as.numeric(rownames(output))),]
  return(output)
}