O exemplo funcionou perfeitamente! A Rosa dos Ventos ficou demais!
É interessante fazer isto para mapas do Brasil!


De: Pedro Henrique Melo Albuquerque <pedro.melo.albuquerque@gmail.com>
Para: r-br@listas.c3sl.ufpr.br
Enviadas: Sexta-feira, 2 de Março de 2012 20:49
Assunto: [R-br] RES: RES: Rosa dos Ventos

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@listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Mauro Sznelwar
Enviada em: Friday, March 02, 2012 8:44 PM
Para: r-br@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@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.
abraço,
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.

--
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://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
Site Pessoal - http://jodavidferreira.com

 

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.