Primeiramente obrigado a todos! Por favor, a função cbind tem um limite para escrita dos dados de um raster, pois meu raster contém 319 linhas x 319 colunas, e a função interrompe no registro 24999, faltando 65001. Esta função é limitada ou gera todos os dados que ela se propõe a ler? A princípio esta função (cbind) atende ao que desejo, pois gera as posições i e j, além de ler o valor do pixel correspondente no raster. Abaixo está uma pequena parte da matriz gerada, após a leitura do meu raster. A coluna com os valores 0 (zero), junto com os valores de row e col são os que desejo num shapefile (tipo polígono ou ponto). matriz row col [1,] 1 1 1 0 [2,] 1 2 2 0 [3,] 1 3 3 0 (....) Por favor, em relação a conversão para shapefile, qual seria a melhor função no R para criar o shapefile com os atributos da matriz acima. Abaixo está o código elaborado por Jônatam, com alterações apenas na quantidade de linhas e colunas, agora resta apenas converter o resultado "matriz" para shapefile. ## criando um raster r <- raster(ncol=319,nrow=319) ## atribuindo valores r[] <- (1:ncell(r)) ## linha, coluna, celula, valor cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r)) # resultado matriz matriz<-cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r)) matriz Mais uma vez agradeço a lista. Att, Alex Santos