riemdist {shapes}

Riemannian shape distance

Description

Calculates the Riemannian shape distance rho between two configurations

Usage

riemdist(x, y, reflect=FALSE)

Arguments

x

k x m matrix (or complex k-vector for 2D data) where k = number of landmarks and m = no of dimensions

y

k x m matrix (or complex k-vector for 2D data)

reflect

Logical. If reflect = TRUE then reflection invariance is included.

Value

The Riemannian shape distance rho between the two configurations. Note 0 <= rho <= pi/2 if no reflection invariance. (for the Riemannian size-and-shape distance use ssriemdist)

Author(s)

Ian Dryden

References

Kendall, D. G. (1984). Shape manifolds, Procrustean metrics and complex projective spaces, Bulletin of the London Mathematical Society, 16, 81-121.

See Also

procOPA,procGPA

Examples

data(gorf.dat)
data(gorm.dat)
gorf<-procGPA(gorf.dat)
gorm<-procGPA(gorm.dat)
rho<-riemdist(gorf$mshape,gorm$mshape)
cat("Riemannian distance between mean shapes is ",rho," \n")