R/whiten.R

In stnava/itkImageR: Examples using Rcpp to interface R and ITK

whiten <- function(x) {
  if (nargs() == 0) {
    print("Usage:  x_whitened<-whiten( x ) ")
    return(1)
  }
  svdx <- svd(x %*% t(x))
  dd <- (svdx$d)^(-1/2)
  xw <- ((svdx$u %*% diag(dd)) %*% t(svdx$v)) %*% x
}