1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
|
\name{Rotation}
\alias{Rotation}
\title{ Rotate a set of point by a certain angle }
\description{
Rotate a set of XY coordinates by an angle (in radians)
}
\usage{
Rotation(xy, angle)
}
\arguments{
\item{xy}{ A 2-columns matrix or data frame containing a set of X and Y coordinates. }
\item{angle}{ Numeric. A scalar giving the angle at which the points should be rotated. The angle is in radians. }
}
\value{
A 2-columns matrix of the same size as \code{xy} giving the rotated coordinates.
}
\author{ F. Guillaume Blanchet }
\examples{
set.seed(1)
### Create a set of coordinates
coords <- cbind(runif(20), runif(20))
### Create a series of angles
rad <- seq(0, pi, l=20)
opar <- par(mfrow=c(5,4), mar=c(3,3,1,1))
for(i in rad){
coords.rot <- Rotation(coords, i)
plot(coords.rot, xlab="", ylab="")
}
par(opar)
### Rotate the coordinates by an angle of 90 degrees
coords.90 <- Rotation(coords, 90*pi/180)
coords.90
plot(coords, xlim=range(rbind(coords.90,coords)[,1]),
ylim=range(rbind(coords.90,coords)[,2]), asp=1)
points(coords.90, pch=19)
}
\keyword{ manip }
|
