<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Boa tarde Walmes,<br>
<br>
Estou tentando postar minha dúvida abaixo e aparece um erro
no servidor, saberia me dizer o porque?<br>
<br>
Obrigado,<br>
<br>
Alexandre<br>
<div class="moz-forward-container"><br>
<br>
-------- Mensagem original --------
<table class="moz-email-headers-table" cellpadding="0"
cellspacing="0" border="0">
<tbody>
<tr>
<th valign="BASELINE" align="RIGHT" nowrap="nowrap">Assunto:
</th>
<td>Representação de gráfico tipo torta em um shapefile</td>
</tr>
<tr>
<th valign="BASELINE" align="RIGHT" nowrap="nowrap">Data: </th>
<td>Sat, 24 May 2014 11:24:58 -0400</td>
</tr>
<tr>
<th valign="BASELINE" align="RIGHT" nowrap="nowrap">De: </th>
<td>ASANTOS <a class="moz-txt-link-rfc2396E" href="mailto:alexandresantosbr@yahoo.com.br"><alexandresantosbr@yahoo.com.br></a></td>
</tr>
<tr>
<th valign="BASELINE" align="RIGHT" nowrap="nowrap">Para: </th>
<td><a class="moz-txt-link-abbreviated" href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a></td>
</tr>
</tbody>
</table>
<br>
<br>
<pre>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(
<a class="moz-txt-link-rfc2396E" href="https://www.dropbox.com/s/enrdm739ts4eba0/sa_bas_ll_r500m.shp">"https://www.dropbox.com/s/enrdm739ts4eba0/sa_bas_ll_r500m.shp"</a>,
<a class="moz-txt-link-rfc2396E" href="https://www.dropbox.com/s/fnnumaix0v88jo1/sa_bas_ll_r500m.shx">"https://www.dropbox.com/s/fnnumaix0v88jo1/sa_bas_ll_r500m.shx"</a>,
<a class="moz-txt-link-rfc2396E" href="https://www.dropbox.com/s/o152vc0n15w54em/sa_bas_ll_r500m.dbf">"https://www.dropbox.com/s/o152vc0n15w54em/sa_bas_ll_r500m.dbf"</a>,
<a class="moz-txt-link-rfc2396E" href="https://www.dropbox.com/s/2l24sso20645511/indv_atributos.csv">"https://www.dropbox.com/s/2l24sso20645511/indv_atributos.csv"</a>)
tokens <- gsub("^.*/s/","",dirname(links))
fileNames <- basename(links)
newLinks <- file.path(<a class="moz-txt-link-rfc2396E" href="http://dl.dropbox.com/s">"http://dl.dropbox.com/s"</a>, 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)
<a class="moz-txt-link-abbreviated" href="mailto:e-mails:alexandresantosbr@yahoo.com.br">e-mails:alexandresantosbr@yahoo.com.br</a>
<a class="moz-txt-link-abbreviated" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a>
Lattes: <a class="moz-txt-link-freetext" href="http://lattes.cnpq.br/1360403201088680">http://lattes.cnpq.br/1360403201088680</a>
======================================================================
</pre>
<br>
</div>
<br>
</body>
</html>