riemdist function

Riemannian shape distance

Calculates the Riemannian shape distance rho between two configurations

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.

Returns

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)

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")

Author(s)

Ian Dryden