[R-br] Criar um novo shapefile à partir da seleção de áreas no shapefile original

ASANTOS alexandresantosbr em yahoo.com.br
Quinta Agosto 20 17:56:45 BRT 2015


Caros listeiros,

       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:

Erro em polypath(x = mcrds[, 1], y = mcrds[, 2], border = border, col = 
col,  :
   Caminho gráfico inválido


*Segue o caminho que estou trilhando:*

#Arquivo shapefile original
contorno_line_a <- readShapeLines ("completo.shp", 
proj4string=CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m +no_defs"))
df.data = as.data.frame(contorno_line_a)

# Seleciono o primeiro polígono
df.serra15 = contorno_line_a[contorno_line_a$Link=='00496015',]## 
Seleciono um talhão no mapa inteiro
contorno_serra15<- 
PolySet2SpatialPolygons(SpatialLines2PolySet(df.serra15), 
close_polys=TRUE) ##
P1 = Polygon (contorno_serra15)
#

# Seleciono o segundo polígono
df.serra16 = contorno_line_a[contorno_line_a$Link=='00496016',]## 
Seleciono um talhão no mapa inteiro
contorno_serra16<- 
PolySet2SpatialPolygons(SpatialLines2PolySet(df.serra16), 
close_polys=TRUE) ##
P2 = Polygon (contorno_serra16)
#

#União dos dois polígonos
Ps1 = Polygons (list (P1,P2), ID = "t")

# Crio uma estrutura em data frame dos polígonos espaciais
SPs = SpatialPolygons (list(Ps1), proj4string=CRS("+proj=utm +zone=23 
+south +datum=WGS84 +units=m +no_defs"))
plot(SPs)
SPDF = SpatialPolygonsDataFrame (SPs, data.frame (N = c ("t"), row.names 
= c("t")))
SPDF em data

#Crio o novo shapefile
writeOGR(SPDF, getwd(), 'parcial', 'ESRI Shapefile')
#

#Chamo o shapefile novo e visualizo
contorno_line_b <- readShapeLines ("parcial.shp", 
proj4string=CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m +no_defs"))
plot(contorno_line_b)
#



Alguém poderia me dar uma ajuda,

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

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150820/d98537c0/attachment.html>


Mais detalhes sobre a lista de discussão R-br