R/orp.r

In zarquon42b/Morpho: Calculations and Visualisations Related to Geometric Morphometrics

orp <- function(A, mshape=NULL) { 
  p <- dim(A)[1]
  k <- dim(A)[2]
  n <- dim(A)[3]
  if (is.null(mshape))
      mshape <- arrMean3(A)
  
  m.size <- cSize(mshape)
  Xc <- as.vector(mshape/m.size)
  X <- vecx(A)/m.size
  ##direction along mshape onto plane
  XtoPlane <- t(apply(X,1,function(x){x <- t(c(crossprod(x,Xc))*Xc)}))
  X1 <- X-XtoPlane
  X1 <- t(X1)+Xc
  proj <- array(X1, dim=c(p, k, n))*m.size
  return(proj)
}