<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Boa noite Éder,<br>
<br>
Tentei seguir seu roteiro usando
maptools::unionSpatialPolygons() e não deu certo a dissolução de
bacias, fiz:<br>
<br>
require(shapefiles)<br>
library(maptools) <br>
library(gpclib) <br>
library(rgdal) <br>
library(PBSmapping)<br>
<br>
<br>
<br>
<br>
#Shapefile e dados<br>
<br>
links <- c(<br>
<br>
<a class="moz-txt-link-rfc2396E" href="https://www.dropbox.com/s/3ph2630xicpexo0/dem.crop.final.tif">"https://www.dropbox.com/s/3ph2630xicpexo0/dem.crop.final.tif"</a>,<br>
<a class="moz-txt-link-rfc2396E" href="https://www.dropbox.com/s/3ph2630xicpexo0/dem.crop.final.tif">"https://www.dropbox.com/s/3ph2630xicpexo0/dem.crop.final.tif"</a>,<br>
<a class="moz-txt-link-rfc2396E" href="https://www.dropbox.com/s/3ph2630xicpexo0/dem.crop.final.tif">"https://www.dropbox.com/s/3ph2630xicpexo0/dem.crop.final.tif"</a>,<br>
<a class="moz-txt-link-rfc2396E" href="https://www.dropbox.com/s/3ph2630xicpexo0/dem.crop.final.tif">"https://www.dropbox.com/s/3ph2630xicpexo0/dem.crop.final.tif"</a>)<br>
<br>
<br>
<br>
if (sum(sapply(basename(links), file.exists))<4) {<br>
<br>
lapply(links, function(a) tryCatch(download.file(a,
dest=basename(a), mode='wb'),<br>
<br>
error=function(...)
message("ERRO!")))} else message('OK!')<br>
#<br>
<br>
#abrido shapefile<br>
#<br>
bacias <- readOGR(".", "sa_bas_ll_r500m")<br>
CRS.new <- CRS("+proj=longlat +datum=WGS84")<br>
proj4string(bacias) <- CRS.new<br>
plot(bacias)<br>
str(bacias)<br>
<br>
## Separando os ID dos atributos do LEVEL2, pois tenho 6 níveis de
bacias<br>
<br>
IDlevel2<- cut(bacias@data$LEVEL2, range(bacias@data$LEVEL2),
include.lowest=TRUE) <br>
<br>
## Bacias "dissolvidas" para o nível 2<br>
<br>
Dissolve_bacia <- unionSpatialPolygons(bacias ,Dissolve_bacia)<br>
<br>
Dissolve_bacia_2 <- SpatialPolygons2PolySet(DissolveResult)<br>
<br>
## Representação gráfica<br>
<br>
plotPolys(Dissolve_bacia_2 , proj =
TRUE,col="wheat1",xlab="longitude",ylab="latitude")<br>
<br>
<br>
<br>
E a figura resultante em plotPolys não são minhas bacias no nível 2.
Poderia me dizer o que estou fazendo de errado,<br>
<br>
Obrigado,<br>
<pre class="moz-signature" cols="72">--
======================================================================
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>
<br>
<br>
<br>
<div class="moz-cite-prefix">Em 26/05/2014 15:12, Éder Comunello
escreveu:<br>
</div>
Alexandre, boa tarde!
<div><br>
</div>
<div>Sugiro que você trabalhe com o mapa de bacias "dissolvidas"
para o nível 2. Para isso você pode fazer uso de
maptools::unionSpatialPolygons(). Para fins de visualização você
pode ainda utilizar outras funções de simplificação.</div>
<div><br>
</div>
<div>Feito isso, além de visualizar mais rápido, você poderá obter o
centróide dos polígonos dissolvidos pra locar os gráficos de
torta. Depois de rodar o sp::over() você precisa tabular as
frequências por bacia e adicionar a coordenada de cada uma.</div>
<div><br>
</div>
<div>A entrada do add.pie() é individual, então você poderá pensar
num laço ou uma função da família apply pra entradas múltiplas.</div>
<div><br>
</div>
<div>Atte.,</div>
<div><br>
</div>
<br clear="all">
<div>
<div dir="ltr">Éder Comunello <<a
href="mailto:comunello.eder@gmail.com" target="_blank">c</a><a
href="mailto:omunello.eder@gmail.com" target="_blank">omunello.eder@gmail.com</a>>
<br>
Dourados, MS - [22 16.5'S, 54 49'W]<br>
</div>
</div>
<pre class="moz-signature" cols="72">
</pre>
</body>
</html>