<div dir="ltr"><div style>## Já tinha esquecido das linhas e colunas.</div><div style><br></div><div style>require(raster)</div><div style>r <- raster(nc=20,nr=10)</div><div style>r[] <- rep(1:5,each=40)<br></div><div style>
plot(r)</div><div style>res <- data.frame(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r))<br></div><div style>names(res) <- c("row","col","cell","v")<br></div><div style>
<br></div><div style>## criando stack com 4 layers</div><div>s <- stack(r,r,r,r)<br></div><div style>## atribuindo valores para cada layer baseado no res</div><div>s[] <- as.matrix(res)<br></div><div style>plot(s)</div>
<div style>p<-rasterToPoints(s, spatial=T)<br></div><div style>spplot(p)<br></div><div style> writeOGR(p,getwd(),"p",driver="ESRI Shapefile")<br></div><div style><div style="font-family:arial,sans-serif;font-size:12.800000190734863px">
pescrito <- readOGR(getwd(), "p")</div><div style="font-family:arial,sans-serif;font-size:12.800000190734863px">pescrito@data</div></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
2013/7/6 Alex <span dir="ltr"><<a href="mailto:alxcart@gmail.com" target="_blank">alxcart@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr"><div class="gmail_default" style="font-family:tahoma,sans-serif">Jônatan,</div><div class="gmail_default" style="font-family:tahoma,sans-serif">a função rasterToshapefile é interessante, mas não atribuí o índices i e j ao shapefile criado.</div>
<div class="gmail_default" style="font-family:tahoma,sans-serif"> Encontrei outra solução no script abaixo e espero auxilar outras pessoas.</div><div class="gmail_default" style="font-family:tahoma,sans-serif"> Apenas uma ressalva, ao final gero um shapefile de pontos, com os atributos i, j e o valor da célula, mas depois associo ao shapefile de polígonos no Quantum GIS através da função Join de Tabelas.</div>
<div class="gmail_default" style="font-family:tahoma,sans-serif"> Outra solução possível é juntar a sua solução de rasterToPolygon e depois realizar o Join de tabela com os shapefiles de pontos.</div><div class="gmail_default" style="font-family:tahoma,sans-serif">
Realmente seria interessante que a função rasterToPolygon levasse as informações de linha e coluna para o shapefile criado.</div><div class="gmail_default" style="font-family:tahoma,sans-serif"><br></div><div class="gmail_default" style="font-family:tahoma,sans-serif">
Mais uma vez obrigado Jôntan e Elias, vocês me orientaram bastante na solução deste problema. </div><div class="gmail_default" style="font-family:tahoma,sans-serif"><br></div><div class="gmail_default"><div class="gmail_default">
<span style="font-family:tahoma,sans-serif">require(raster)</span><br></div><div class="gmail_default"><font face="tahoma, sans-serif">require(rgdal)</font></div><div class="gmail_default"><span style="font-family:tahoma,sans-serif">require(shapefiles)</span><br>
</div><div class="gmail_default"><span style="font-family:tahoma,sans-serif">require(maptools)</span><br></div><div class="gmail_default"><span style="font-family:tahoma,sans-serif"><br></span></div><div class="gmail_default">
<span style="font-family:tahoma,sans-serif">### ENTRADA DE DADOS ###</span><br></div><div class="im"><div class="gmail_default"><span style="font-family:tahoma,sans-serif">## criando um raster </span><br></div><div class="gmail_default">
<font face="tahoma, sans-serif">r <- raster(ncol=319,nrow=319)</font></div></div><div class="im"><div class="gmail_default"><font face="tahoma, sans-serif">atribuindo valores</font></div><div class="gmail_default"><font face="tahoma, sans-serif">r[] <- (1:ncell(r))</font></div>
<div class="gmail_default"><br></div></div><div class="im"><div class="gmail_default"><font face="tahoma, sans-serif">## linha, coluna, celula, valor</font></div><div class="gmail_default"><font face="tahoma, sans-serif">cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r))</font></div>
</div><div class="im"><div class="gmail_default"><font face="tahoma, sans-serif">matriz<-cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r))</font></div></div><div class="gmail_default"><font face="tahoma, sans-serif">matriz</font></div>
<div class="gmail_default"><font face="tahoma, sans-serif">dim(matriz)</font></div><div class="gmail_default"><font face="tahoma, sans-serif">tail(matriz)</font></div><div class="gmail_default"><font face="tahoma, sans-serif">summary (matriz)</font></div>
<div class="gmail_default"><font face="tahoma, sans-serif">names(matriz) # a matriz não possui nome de campos</font></div><div class="gmail_default"><font face="tahoma, sans-serif"><br></font></div><div class="gmail_default">
<span style="font-family:tahoma,sans-serif"># Data Frame possui nome dos atributos</span><br></div><div class="gmail_default"><span style="font-family:tahoma,sans-serif">class(matriz)</span><br></div><div class="gmail_default">
<font face="tahoma, sans-serif">require(sp)</font></div><div class="gmail_default"><font face="tahoma, sans-serif"> pontos<-data.frame(matriz)</font></div><div class="gmail_default"><font face="tahoma, sans-serif"> head(pontos)</font></div>
<div class="gmail_default"><font face="tahoma, sans-serif"> class(pontos)</font></div><div class="gmail_default"><font face="tahoma, sans-serif"> dim(pontos)</font></div><div class="gmail_default"><font face="tahoma, sans-serif"> names(pontos)</font></div>
<div class="gmail_default"><font face="tahoma, sans-serif"> coordinates(pontos) <- c("row","col") # c("Norte_Sul","Leste_Oeste")</font></div><div class="gmail_default"><font face="tahoma, sans-serif"> class(pontos)</font></div>
<div class="gmail_default"><font face="tahoma, sans-serif"> dim(pontos)</font></div><div class="gmail_default"><font face="tahoma, sans-serif"> getSlots("SpatialPointsDataFrame")</font></div><div class="gmail_default">
<font face="tahoma, sans-serif"> slot(pontos, "bbox")</font></div><div class="gmail_default"><font face="tahoma, sans-serif"> </font></div><div class="gmail_default"><font face="tahoma, sans-serif">require(maptools)</font></div>
<div class="gmail_default"><font face="tahoma, sans-serif">writePointsShape(pontos,"matriz319pt", factor2char = TRUE, max_nchar=254) </font></div><div class="gmail_default"><font face="tahoma, sans-serif">readShapePoints("matriz319pt", proj4string = CRS(as.character(NA)),verbose =</font></div>
<div class="gmail_default"><font face="tahoma, sans-serif">FALSE,repair=FALSE)</font></div><div class="gmail_default"><font face="tahoma, sans-serif"><br></font></div><div class="gmail_default"><br></div></div></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>
<div dir="ltr"><div><font style="font-family:arial">###############################################################</font><br><span style="font-family:arial;font-size:small">## </span><span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px">Jônatan Dupont Tatsch</span></div>
<div><span style="font-family:arial;font-size:small">## </span><span style="color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;font-size:10.909090995788574px;line-height:12.997159004211426px;text-align:center;background-color:rgb(239,239,239)">Professor do Departamento de Física</span><font face="arial" style="font-size:small"><br>
</font></div><div><span style="font-size:small;font-family:arial">## </span><span style="color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;font-size:10.909090995788574px;line-height:12.997159004211426px;text-align:center;background-color:rgb(239,239,239)">Centro de Ciências Exatas e Naturais (CCNE)</span><br>
</div><div><span style="font-size:small;font-family:arial">## </span><span style="color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;font-size:10.909090995788574px;line-height:12.997159004211426px;text-align:center;background-color:rgb(239,239,239)">Universidade Federal de Santa Maria</span></div>
<div><span style="font-family:arial;font-size:small">## </span><span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px">Faixa de Camobi, Prédio 13 - Campus UFSM - </span><span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px">Santa Maria, RS, Brasil - 97105-900</span><br>
</div><div style="font-family:arial;font-size:small">## <span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px">Telefone: +55(55)33012084</span></div>
<div style="font-family:arial;font-size:small">## <span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px"><a href="http://www.ufsm.br/meteorologia" target="_blank">www.ufsm.br/meteorologia</a></span></div>
<div style="font-family:arial;font-size:small">###############################################################</div></div>
</div>