<html>
  <head>
    <meta content="text/html; charset=windows-1252"
      http-equiv="Content-Type">
  </head>
  <body bgcolor="#FFFFFF" text="#000000">
    Caros listeiros,<br>
    <br>
          Ainda estou enrolado com o meu problema e gostaria de ver se
    alguém conseguiria me ajudar, quando seleciono 1 polígono
    (df.serra15) e faço data.frame(ID=1:length(df.serra15)), seguido da
    conversão para SpatialPolygonsDataFrame funciona a transformação
    para shapefile, mas se seleciono dois polígonos
    df.serra15e16<-c(list(df.serra15),list(df.serra16)), recebo o
    erro:   tentativa de obter um slot "polygons" de um objeto de uma
    classe básica ("list") sem slots,<br>
    <br>
         Segue o exemplo que fiz abaixo, sem CRM pois não queria que
    ninguém tivesse que fazer download de nenhum shapefile, obrigado<br>
    <br>
    <br>
    #Tentativa 2
    -------------------------------------------------------------------<br>
    #Pacotes utilizados<br>
    require(sp); require(rgdal); require(shapefiles)<br>
    <br>
    #Arquivo shapefile original<br>
    contorno_line_a <- readShapeLines ("Mosaico_Geral.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 1 polígono<br>
    df.serra15 = contorno_line_a[contorno_line_a$Link=='00496015',]##<br>
    p.df <- data.frame(ID=1:length(df.serra15)) <br>
    p <- SpatialPolygonsDataFrame(contorno_serra15,p.df)<br>
    writeOGR(srdf, "serra15", layer="serra15", driver="ESRI Shapefile")<br>
    <br>
    # Seleciono 2 polígono<br>
    df.serra15 = contorno_line_a[contorno_line_a$Link=='00496015',]##<br>
    df.serra16 = contorno_line_a[contorno_line_a$Link=='00496016',]<br>
    df.serra15e16<-c(list(df.serra15),list(df.serra16))<br>
    p.df2 <- data.frame(ID=1:length(df.serra15e16)) <br>
    p2 <- SpatialPolygonsDataFrame(df.serra15e16,p.df2)<br>
    writeOGR(srdf, "serra15e16", layer="serra15e16", driver="ESRI
    Shapefile")<br>
    #<br>
    Erro em SpatialPolygonsDataFrame(df.serra15e16, p.df2) : <br>
      tentativa de obter um slot "polygons" de um objeto de uma classe
    básica ("list") sem slots<br>
    <br>
    <br>
       <br>
    <br>
           Ainda<br>
    <br>
    <div class="moz-cite-prefix">Em 20/08/2015 16:56, ASANTOS escreveu:<br>
    </div>
    <blockquote cite="mid:55D63F0D.4000509@yahoo.com.br" type="cite">
      <meta http-equiv="content-type" content="text/html;
        charset=windows-1252">
      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 moz-do-not-send="true" class="moz-txt-link-abbreviated" href="mailto:e-mails:alexandresantosbr@yahoo.com.br">e-mails:alexandresantosbr@yahoo.com.br</a> 
        <a moz-do-not-send="true" class="moz-txt-link-abbreviated" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a> 
Lattes: <a moz-do-not-send="true" class="moz-txt-link-freetext" href="http://lattes.cnpq.br/1360403201088680">http://lattes.cnpq.br/1360403201088680</a>                   
======================================================================
</pre>
      <br>
      <fieldset class="mimeAttachmentHeader"></fieldset>
      <br>
      <pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
    </blockquote>
    <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>