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@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 ======================================================================