<html><head></head><body><div class="ydp90cb4a93yahoo-style-wrap" style="font-family:Helvetica Neue, Helvetica, Arial, sans-serif;font-size:13px;"><div><div>Não seria o caso de reduzir a resolução do shapefile? Já teve situaçoes que usei gSimplify() e escolhi parametros que se mostraram adequados ao meu objetivo.<br></div><div><br></div><div class="ydp90cb4a93signature">Elias T. Krainski</div></div>
<div><br></div><div><br></div>
</div><div id="yahoo_quoted_8571321397" class="yahoo_quoted">
<div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;">
<div>
Em quarta-feira, 16 de janeiro de 2019 12:08:04 BRST, Yury Duarte <yurynepomuceno@gmail.com> escreveu:
</div>
<div><br></div>
<div><br></div>
<div><div id="yiv3081619695"><div><div dir="ltr">Bom dia Elias, como vai?<div><br clear="none"></div><div>Interessante sua solução!</div><div>Não estou acostumado a salvar minhas saídas no formato .RData. Irei testar esse formato nos casos onde irei trabalhar os dados/analises apenas em ambiente R.</div><div>Aparentemente, esse formato economiza bastante espaço em disco.<br clear="none">Entretanto, para essa situação em especial, eu realmente tenho a necessidade de salvar as saídas em formato .shp, pois tenho que aplicar o resultado do script em outros ambientes (nos SIGs em geral e no Google Earth). </div><div><br clear="none"></div><div>Caso tenha(m) mais alguma sugestão, ficarei feliz em testar.</div><div><br clear="none"></div><div>Mais uma vez, agradeço pela disponibilidade e apoio!</div><div><br clear="none"></div><div>Att </div><div><br clear="all"><div><div class="yiv3081619695gmail_signature" dir="ltr"><div dir="ltr"><div>Yury Duarte<br clear="none"></div>Engenheiro Agrônomo - ESALQ/USP<br clear="none"></div></div></div><br clear="none"></div></div><br clear="none"><div class="yiv3081619695yqt4394568671" id="yiv3081619695yqt57013"><div class="yiv3081619695gmail_quote"><div class="yiv3081619695gmail_attr" dir="ltr">Em qua, 16 de jan de 2019 às 11:24, Elias T. Krainski <<a rel="nofollow" shape="rect" ymailto="mailto:eliaskrainski@yahoo.com.br" target="_blank" href="mailto:eliaskrainski@yahoo.com.br">eliaskrainski@yahoo.com.br</a>> escreveu:<br clear="none"></div><blockquote class="yiv3081619695gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex;"><div><div class="yiv3081619695gmail-m_5529097057479592537ydpf430d47eyahoo-style-wrap" style=""><div><div>se você for usar apenas em ambiente R, a melhor opção é usar save(...). Frequentemente uso com compress='xz'.<br clear="none"></div><div><br clear="none"></div><div class="yiv3081619695gmail-m_5529097057479592537ydpf430d47esignature">Elias T. Krainski</div></div>
<div><br clear="none"></div><div><br clear="none"></div>
</div><div class="yiv3081619695gmail-m_5529097057479592537yahoo_quoted" id="yiv3081619695gmail-m_5529097057479592537yahoo_quoted_8407030181">
<div style="">
<div>
Em terça-feira, 15 de janeiro de 2019 11:30:21 BRST, Yury Duarte por (R-br) <<a rel="nofollow" shape="rect" ymailto="mailto:r-br@listas.c3sl.ufpr.br" target="_blank" href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>> escreveu:
</div>
<div><br clear="none"></div>
<div><br clear="none"></div>
<div><div id="yiv3081619695gmail-m_5529097057479592537yiv0746425720"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Bom dia colegas listeiros!<div><br clear="none"></div><div>Ultimamente venho tentando manusear arquivos espaciais utilizando o R.</div><div>Em alguns casos, objetos do tipo 'large SpatialPolygons' são gerados dentro do código. Quando isso ocorre (acredito que seja algo que acontece por padrão, devido a grande quantidade de informações armazenadas no objeto), a função writeOGR, que utilizo para salvar os shapefiles, quebra, por só tratar de objetos do tipo 'Spatial' e não os 'large Spatial'.<br clear="none">Sendo assim, gostaria de saber se existe a possibilidade de salvar um shapefile a partir de um objeto do tipo 'large Spatial' ou ainda se existe a possibilidade de transformar esse num objeto do tipo SpatialPolygonsDataFrame, para que a função writeOGR possa funcionar normalmente.</div><div><br clear="none"></div><div>Segue código desenvolvido seguido do erro gerado. <br clear="none">Desde já, agradeço pela ajuda de todos!</div><div><br clear="none"></div><div><div>rodar_bibliotecas = function(necessarias_para_o_projeto){</div><div> </div><div> if(necessarias_para_o_projeto){</div><div> library(rgrass7)</div><div> library(spatstat) </div><div> library(maptools) </div><div> library(shapefiles)</div><div> library(foreign)</div><div> library(magrittr)</div><div> library(formattable)</div><div> library(e1071)</div><div> library(rlang)</div><div> library(rgdal)</div><div> library(rgeos)</div><div> library(raster)</div><div> library(sp)</div><div> library(sf)</div><div> library(RcppCNPy)</div><div> library(deldir)</div><div> library(dismo)</div><div> library(dplyr)</div><div> library(ggplot2)</div><div> library(gstat)</div><div> library(tidyverse)</div><div> library(smoothr)</div><div> library(viridisLite)</div><div> }</div><div> return('Ok')</div><div>}</div></div><div><br clear="none"></div><div><div>rodar_bibliotecas(TRUE)</div><div><br clear="none"></div><div>raiz = '/home/yury/pesquisa/arquivos/areas/'</div><div>clientes = dir(raiz)</div><div><br clear="none"></div><div>for (cliente in 1:length(clientes)) {</div><div> </div><div> analises = dir(paste0(raiz, clientes[cliente]))</div><div> q = list()</div><div> </div><div> for (analise in 1:length(analises)) {</div><div> </div><div> shapefile_contorno = readOGR(paste0(raiz, clientes[cliente], '/', analises[analise], '/', 'vetores.shp'))</div><div> shapefile_contorno$id = 1</div><div> dissolvido = unionSpatialPolygons(shapefile_contorno, IDs = shapefile_contorno$id)</div><div> #plot(dissolvido)</div><div> </div><div> q[analise] = dissolvido</div><div> </div><div> }</div><div> </div><div> fim = q[[1]]</div><div> for (x in 2:length(q)){</div><div> </div><div> fim = bind(fim, q[[x]])</div><div> </div><div> }</div><div>plot(fim)</div><div><br clear="none"></div><div>path_to_save = paste0(raiz, clientes[cliente])</div><div>writeOGR(fim,</div><div> dsn = path_to_save,</div><div> layer = paste0('contorno_', clientes[cliente]),</div><div> driver = 'ESRI Shapefile',</div><div> overwrite_layer = TRUE)</div><div><br clear="none"></div><div>}</div></div><div><span style="color:rgb(197,6,11);font-size:10.4pt;white-space:pre-wrap;"><br clear="none"></span></div><div><span style="color:rgb(197,6,11);font-size:10.4pt;white-space:pre-wrap;">Error in writeOGR(fim, dsn = path_to_save, layer = paste0("contorno_", : obj must be a SpatialPointsDataFrame, SpatialLinesDataFrame or SpatialPolygonsDataFrame</span><br clear="none"></div><div><br clear="all"><div><div class="yiv3081619695gmail-m_5529097057479592537yiv0746425720gmail_signature" dir="ltr"><div dir="ltr"><div>Yury Duarte<br clear="none"></div>Engenheiro Agrônomo - ESALQ/USP<br clear="none"></div></div></div></div></div></div></div></div></div>_______________________________________________<br clear="none">R-br mailing list<br clear="none"><a rel="nofollow" shape="rect" ymailto="mailto:R-br@listas.c3sl.ufpr.br" target="_blank" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br clear="none"><a rel="nofollow" shape="rect" target="_blank" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br clear="none">Leia o guia de postagem (<a rel="nofollow" shape="rect" target="_blank" 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.</div>
</div>
</div></div></blockquote></div></div></div></div></div>
</div>
</div></body></html>