Thiago e Elias,

primeiramente obrigado!

Os polígonos são uma matriz com quadriculas que cobrem uma região e sendo que algumas foram selecionadas de forma sistemática. Em cada quadrícula selecionada houve uma análise espacial formando um índice (por exemplo: um percentual de erro).

Como vou aplicar a geoestatística, para estimar quadrículas não amostradas (constatei pelo variograma que há correlação espacial) preciso que as quadrículas sejam contínuas, logo o índice das células é crucial.

Ainda não testei a sugestão do Elias e em relação a função rasterToPolygons não encontrei um parâmetro para extrair a informação do índice das células.

Mas considero ser uma demanda relevante, extrair todas as posições i e j de um raster (junto com o valor do respectivo pixel), como já mencionei inexistente em muitos softwares.

Att,

Alex Santos