-------- 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@yahoo.com.br>
Para:	r-br@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@bbox
pontos@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@yahoo.com.br
        alexandre.santos@cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================