[R-sig-Geo] Principal Directions with surf.ls

[Gabor Grothendieck]

Suppose we plot a quadratic surface as in:

library(spatial)
example(surf.ls)

How can we plot lines on the graph to show the
principal directions/axes of the elliptical contours
and how do we get them?

Thanks.