<div dir="ltr"><div class="gmail_default" style="font-family:tahoma,sans-serif">Jônatan, excelente! </div><div class="gmail_default" style="font-family:tahoma,sans-serif">no seu script apenas substitui rasterToPoints para rasterToPolygons no final.</div>
<div class="gmail_default" style="font-family:tahoma,sans-serif">O objetivo do raster é apenas para extrair as informações do índice (i, j). Os valores do pixel funcionam como ID's. </div><div class="gmail_default" style="font-family:tahoma,sans-serif">
Agora, o resultado final já é georreferenciado e as quadrículas estão prontas para uma posterior análise geoestatística.</div><div class="gmail_default" style="font-family:tahoma,sans-serif">Att,</div><div class="gmail_default" style="font-family:tahoma,sans-serif">
Alex Santos</div><div class="gmail_default" style="font-family:tahoma,sans-serif"><br></div><div class="gmail_default"><div class="gmail_default"><div class="gmail_default"><font face="tahoma, sans-serif">### ENTRADA DE DADOS ###</font></div>
<div class="gmail_default"><font face="tahoma, sans-serif">r <- raster(nc=20,nr=10)</font></div><div class="gmail_default"><font face="tahoma, sans-serif">r[] <- rep(1:5,each=40)</font></div><div class="gmail_default">
<font face="tahoma, sans-serif">plot(r)</font></div><div class="gmail_default"><font face="tahoma, sans-serif">res <- data.frame(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r))</font></div><div class="gmail_default">
<font face="tahoma, sans-serif">names(res) <- c("row","col","cell","v")</font></div><div class="gmail_default"><font face="tahoma, sans-serif"><br></font></div><div class="gmail_default">
<font face="tahoma, sans-serif">## criando stack com 4 layers</font></div><div class="gmail_default"><font face="tahoma, sans-serif">s <- stack(r,r,r,r)</font></div><div class="gmail_default"><font face="tahoma, sans-serif">## atribuindo valores para cada layer baseado no res</font></div>
<div class="gmail_default"><font face="tahoma, sans-serif">s[] <- as.matrix(res)</font></div><div class="gmail_default"><font face="tahoma, sans-serif">plot(s)</font></div><div class="gmail_default"><font face="tahoma, sans-serif">#p<-rasterToPoints(s, spatial=T)</font></div>
<div class="gmail_default"><font face="tahoma, sans-serif">pp<- rasterToPolygons(s,dissolve=F)</font></div><div class="gmail_default"><font face="tahoma, sans-serif">#spplot(pp)</font></div><div class="gmail_default"><font face="tahoma, sans-serif"># writeOGR(p,getwd(),"p2",driver="ESRI Shapefile")</font></div>
<div class="gmail_default"><font face="tahoma, sans-serif"> writeOGR(pp,getwd(),"matriz_pl",driver="ESRI Shapefile")</font></div><div class="gmail_default"><font face="tahoma, sans-serif">pescrito <- readOGR(getwd(), "matriz_pl")</font></div>
<div class="gmail_default"><font face="tahoma, sans-serif">pescrito@data</font></div></div></div></div>