[R-br] Fwd: Representação de gráfico tipo torta em um shapefile

[ASANTOS]

Boa tarde Walmes,

         Estou tentando postar minha dúvida abaixo e aparece um erro no 
servidor, saberia me dizer o porque?

Obrigado,

Alexandre


-------- Mensagem original --------
Assunto: 	Representação de gráfico tipo torta em um shapefile
Data: 	Sat, 24 May 2014 11:24:58 -0400
De: 	ASANTOS <alexandresantosbr em yahoo.com.br>
Para: 	r-br em listas.c3sl.ufpr.br



Bom dia Pessoal,

                  Estou querendo inserir um gráfico do tipo torta em um
shapefile de bacias hidrográficas da América do Sul (bacias), fazendo a
inserção de um gráfico torta para cada bacia , onde localizam-se meus
pontos de interesse (pontos$LEVEL2), porém os valores do gráfico
gostaria que fosse a frequência de pontos$CT em cada pontos$LEVEL2, mas
estou tento dificuldades e vou explicar passo a passo:

# Start
------------------------------------------------------------------------------------------------------
#Pacotes
require (maptools)
require(shapefiles)
require(rgdal)
require(mapplots)
  require(RColorBrewer)
#

#Shapefile e dados

links <- c(

"https://www.dropbox.com/s/enrdm739ts4eba0/sa_bas_ll_r500m.shp",
"https://www.dropbox.com/s/fnnumaix0v88jo1/sa_bas_ll_r500m.shx",
"https://www.dropbox.com/s/o152vc0n15w54em/sa_bas_ll_r500m.dbf",
"https://www.dropbox.com/s/2l24sso20645511/indv_atributos.csv")


tokens    <- gsub("^.*/s/","",dirname(links))
fileNames <- basename(links)
newLinks  <- file.path("http://dl.dropbox.com/s", tokens, fileNames);
newLinks

for (a in newLinks) {
           tryCatch(download.file(a, dest=basename(a), mode='wb'),
                               error=function(...) print("Falha no
download!"))}

#



#Abrindo o shapefile
#
CRS.new <- CRS("+proj=longlat +datum=WGS84")
bacias <- readOGR(".", "sa_bas_ll_r500m")
proj4string(bacias) <- CRS.new
plot(bacias)
#

#
#Coordenadas  dos pontos + atributos
pontos<-read.csv("indv_atributos.csv", sep=";",h=T)
coordinates(pontos) <- c("x_long_dec","y_lat_dec")
proj4string(pontos) <- CRS.new
#

# Selecionos os pontos que estão dentro de cada bacia hidrográfica
inside.sa <- !is.na(over(pontos, as(bacias, "SpatialPolygons")))##
inside.sa
#


# Pontos pertencentes a cada bacia nível 2, tenho varios niveis de
classificação e escolhi o LEVEL2
bacias em bbox
pontos em bbox
pontos$LEVEL2 <- over(pontos, bacias)$LEVEL2



## Adiciono o gráfico torta para cada pontos$LEVEL2 a frequencia de
pontos$CT

add.pie(z=pontos$CT, x=pontos$x_long_dec, y=pontos$y_lat_dec, radius=1,
col=c(alpha("orange", 0.6), alpha("blue", 0.6), alpha("black", 0.6),
alpha("red", 0.6), alpha("grey", 0.6)), labels=" ")

#END----------------------------------------------------------------------------------------------------------



Mas não obtive sucesso, alguém poderia me dar um help?

Obrigado,

-- 
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT                      CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM)   (+55) 65 9686-6970 (VIVO)
e-mails:alexandresantosbr em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================



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