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
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
|
\name{vectorField}
\alias{vectorField}
\title{Display magnitude/direction vectors}
\description{Display magnitude/direction vectors as arrows on an existing plot.}
\usage{
vectorField(u,v,xpos=NA,ypos=NA,scale=1,headspan=0.1,
vecspec=c("lonlat","rad","deg"))
}
\arguments{
\item{u,v}{x (longitude) and y (latitude) offsets OR
orientation and magnitude in either radians or degrees. See details.}
\item{xpos,ypos}{The centers of the vectors in user units.}
\item{scale}{The proportion of each cell that the maximal vector will fill.
See details.}
\item{headspan}{The extent of the heads of the arrows as a proportion of
cell size.}
\item{vecspec}{How the vectors are described. See details}
}
\value{ nil }
\details{
\samp{vectorField} displays arrows on an existing plot. Each arrow is specified
by a position on the plot \samp{xpos,ypos} and either x/y offsets or
orientation and magnitude. The default is x/y offsets, and the user must
specify whether radians or degrees are used if the orientation/magnitude option
is used.
If the first four arguments are matrices, there must be no missing values. If
these arguments are vectors, the calculation of the scaling of the magnitudes
and length of the arrowheads may be slightly different.
}
\author{Jim Lemon (original code by Robin Hankin and Brian Ripley)}
\seealso{\link{arrows}}
\examples{
\dontrun{
# this requires the maps package, and just wouldn't pass check
require(maps)
map("world",xlim=c(110,155),ylim=c(-40,-10))
par(xpd=TRUE)
text(132,-5,"Approximate magnetic deviation - Australia",cex=1.5)
par(xpd=FALSE)
long<-rep(seq(117.5,152.5,by=5),6)
lat<-rep(c(-12.5,-17.5,-22.5,-27.5,-32.5,-37.5),each=8)
# just show the direction, don't have a magnitude difference
mag<-rep(1,48)
devdeg<-c(110,98,85,65,65,65,65,65,
115,100,90,80,72,66,63,55,
130,100,90,82,72,67,62,54,
122,111,95,86,70,67,56,48,
118,116,110,87,74,68,62,45,
128,115,107,90,78,66,53,45)
vectorField(devdeg,mag,long,lat,scale=0.7,vecspec="deg")
}
# do a magnitude/direction plot with radians
plot(1:10,type="n",main="Random vectors")
mag<-runif(100)+1
dir<-runif(100)*2*pi
xpos<-rep(1:10,10)
ypos<-rep(1:10,each=10)
vectorField(dir,mag,xpos,ypos,scale=0.8,vecspec="rad")
}
\keyword{misc}
|
