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.