[R-br] RES: RES: Rosa dos Ventos

[Pedro Henrique Melo Albuquerque]

Mauro,

Instale o pacote spdep e tente novamente.

Caso queira fazer de maneira direta você pode baixar o shapefile nesse link:
http://geodacenter.org/downloads/data-files/columbus.zip

 

 

De: r-br-bounces em listas.c3sl.ufpr.br
[mailto:r-br-bounces em listas.c3sl.ufpr.br] Em nome de Mauro Sznelwar
Enviada em: Friday, March 02, 2012 8:44 PM
Para: r-br em listas.c3sl.ufpr.br
Assunto: Re: [R-br] RES: Rosa dos Ventos

 

Eu tentei rodar este exemplo e não consegui! Já deu erro no terceiro
comando:

columbus <- readShapePoly(system.file("etc/shapes/columbus.shp",
package="spdep")[1])
Erro em getinfo.shape(filen) : Error opening SHP file

 

 

Jodavid,
Rodei o exemplo do Pedro no meu PC e funcionou. Observe que a função
"RosaDosVentos (c(5,13),0.4)" substitua o valor 5 por 6.5 que aparecera a
Rosa dos ventos

Humberto


Em 3/2/2012 6:30 PM, Pedro Henrique Melo Albuquerque escreveu: 

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 em listas.c3sl.ufpr.br
[mailto:r-br-bounces em listas.c3sl.ufpr.br] Em nome de Jodavid Ferreira
Enviada em: Friday, March 02, 2012 3:13 PM
Para: r-br em 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.

abraço,

Em 27 de fevereiro de 2012 14:01, Jodavid Ferreira <jodavid.arts em 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.




-- 
Jodavid Ferreira
Responsável pela comissão de Comunicação do MJ da RCCPB -
http://jovem.rccpb.com.br
Responsável pela comissão de Comunicação e do Ministério Jovem do Grupo de
Oração Boa Esperança - http://grupoboaesperanca.com
Bacharelando em Estatística - UFPB
Estagiário no Laboratório de Estatística Aplicada ao Processamento de
Imagens e Geoprocessamento - http://de.ufpb.br/~leapig
<http://de.ufpb.br/%7Eleapig> 
http://jodavidferreira.com





-- 

Jodavid Ferreira
Responsável pela comissão de Comunicação do MJ da RCCPB -
http://jovem.rccpb.com.br
Responsável pelo Ministério Jovem e de Comunicação do Grupo de Oração Boa
Esperança - http://grupoboaesperanca.com
Bacharelando em Estatística - UFPB
Estagiário no Laboratório de Estatística Aplicada ao Processamento de
Imagens e Geoprocessamento - http://de.ufpb.br/~leapig
<http://de.ufpb.br/%7Eleapig> 
Site Pessoal - http://jodavidferreira.com


 

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120302/3bb35595/attachment.html>