Jodavid,
Segue um
exemplo:
require(maptools)
require(classInt)
columbus
<- readShapePoly(system.file("etc/shapes/columbus.shp",
package="spdep")[1])
mypal <-
c("wheat1", "red3")
f5 <-
classIntervals(columbus$CRIME, n=5, style="fisher")
f5Colours
<- findColours(f5, mypal)
plot(columbus,
col=f5Colours,border ="black", axes = TRUE)
title(main="Criminalidade
em Columbus.",xlab="Longitude",ylab="Latitude")
RosaDosVentos
<- function(loc,size,bearing=0,cols,cex=1,...)
{
if(missing(loc)) stop("loc is missing")
if(missing(size)) stop("size is missing")
if(missing(cols)) cols <- rep(c("white","black"),8)
radii
<- rep(size/c(1,4,2,4),4)
x <-
radii[(0:15)+1]*cos((0:15)*pi/8+bearing)+loc[1]
y <-
radii[(0:15)+1]*sin((0:15)*pi/8+bearing)+loc[2]
for (i in
1:15)
{
x1 <-
c(x[i],x[i+1],loc[1])
y1 <-
c(y[i],y[i+1],loc[2])
polygon(x1,y1,col=cols[i])
}
polygon(c(x[16],x[1],loc[1]),c(y[16],y[1],loc[2]),col=cols[16])
b <-
c("L","N","O","S")
for (i in
0:3) text((size+par("cxy")[1])*cos(bearing+i*pi/2)+loc[1],
(size+par("cxy")[2])*sin(bearing+i*pi/2)+loc[2],b[i+1],cex=cex)
}
RosaDosVentos
(c(5,13),0.4)
Por
favor, informe se funcionou.
Grande
abraço.
Pedro
De: r-br-bounces@listas.c3sl.ufpr.br
[mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Jodavid
Ferreira
Enviada em: Friday, March 02, 2012 3:13 PM
Para: r-br@listas.c3sl.ufpr.br
Assunto: Re: [R-br] Rosa dos Ventos
Olá Pedro, Eu estou até entendendo o seu
código.
Mas ainda não consegui colocar a Rosa
dos Ventos.
O problema está, em como fazer ela
aparecer na plotagem.
se ela vai como legenda...
se você puder me enviar algum exemplo
dela funcionando agradeceria bastante.
Em 27 de fevereiro de 2012 14:01,
Jodavid Ferreira <jodavid.arts@gmail.com>
escreveu:
Olá Pedro,
Rodei seu código,ok. Mas não
estou conseguindo visualizar a Rosa.
Não estou conseguindo coloca-la
no mapa, ou até mesmo, plotá-la.
--