<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="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 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 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 class="gmail_default"><font face="tahoma, sans-serif">matriz<-cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r))</font></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>