<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=windows-1252">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Caros listeiros,<br>
<br>
Tenho um shapefile com inúmeras áreas e gostaria de selecionar
apenas duas e criar um novo shapefile, no entanto, quanto tento unir
os dois polígonos e visualizar a seleção, recebo o erro:<br>
<br>
Erro em polypath(x = mcrds[, 1], y = mcrds[, 2], border = border,
col = col, : <br>
Caminho gráfico inválido<br>
<br>
<br>
<b>Segue o caminho que estou trilhando:</b><br>
<br>
#Arquivo shapefile original<br>
contorno_line_a <- readShapeLines ("completo.shp",
proj4string=CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m
+no_defs"))<br>
df.data = as.data.frame(contorno_line_a)<br>
<br>
# Seleciono o primeiro polígono<br>
df.serra15 = contorno_line_a[contorno_line_a$Link=='00496015',]##
Seleciono um talhão no mapa inteiro<br>
contorno_serra15<-
PolySet2SpatialPolygons(SpatialLines2PolySet(df.serra15),
close_polys=TRUE) ##<br>
P1 = Polygon (contorno_serra15)<br>
#<br>
<br>
# Seleciono o segundo polígono<br>
df.serra16 = contorno_line_a[contorno_line_a$Link=='00496016',]##
Seleciono um talhão no mapa inteiro<br>
contorno_serra16<-
PolySet2SpatialPolygons(SpatialLines2PolySet(df.serra16),
close_polys=TRUE) ##<br>
P2 = Polygon (contorno_serra16)<br>
#<br>
<br>
#União dos dois polígonos<br>
Ps1 = Polygons (list (P1,P2), ID = "t")<br>
<br>
# Crio uma estrutura em data frame dos polígonos espaciais<br>
SPs = SpatialPolygons (list(Ps1), proj4string=CRS("+proj=utm
+zone=23 +south +datum=WGS84 +units=m +no_defs"))<br>
plot(SPs)<br>
SPDF = SpatialPolygonsDataFrame (SPs, data.frame (N = c ("t"),
row.names = c("t")))<br>
SPDF@data<br>
<br>
#Crio o novo shapefile<br>
writeOGR(SPDF, getwd(), 'parcial', 'ESRI Shapefile')<br>
#<br>
<br>
#Chamo o shapefile novo e visualizo<br>
contorno_line_b <- readShapeLines ("parcial.shp",
proj4string=CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m
+no_defs"))<br>
plot(contorno_line_b)<br>
#<br>
<br>
<br>
<br>
Alguém poderia me dar uma ajuda,<br>
<br>
Obrigado, <br>
<br>
<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>
</body>
</html>