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