Re: [R-br] Como extrair índices linha e coluna de um raster

Elias, informo que minhas amostas não são pixels. Minha metodologia está consolidada, preciso apenas de uma diretriz para encontrar está função no R (pois acredito que o mesmo seja capaz para realizar o que desejo. Abaixo, está em linha gerais minha metodologia, mas por favor, não se atente a mesma pois meu foco é encontra/contrair a função que extraia o índice (i,j) de cada pixel. 1 - Eu tenho uma articulação do mapeamento sistemático que recobre minha área de interesse (no formato shapefile - poligono). Sobre estas áreas que houve a amostragem sistemática. (obs: Numa solução paralela utilizando o FME (software comercial) já agreguei as informações de (i,j) em cada quadrícula, podendo utilizar esta informação como X, Y e representa minhas quadrículas de forma contínua. Mas quero substituir o FME (comercial) pelo R (livre), pois acredito no potencial do R e acho ele capaz de executar esta tarefa. 2 - converti a articulação para raster, já com a proposta de encontrar as informações do índice (i,j) de cada pixel. Sendo que cada pixel possui o mesmo tamanho de minhas quadrículas. Mas espero que compreenda, pois o que desejo realmente é uma função que leia um arquivo raster - extraia o índice (i,j) de cada pixel (compondo uma matriz com os índices mais o valor do atributo) e exporte esta matriz para um shapefile. obs.: Adianto que o mesmo não será georreferenciado e nem há necessidade neste momento, pois haverá o georreferenciamento após a krigagem através dos extremos da articulação original. Att, Alex Santos

On 06/20/2013 11:48 AM, Alex wrote:
já agreguei as informações de (i,j) em cada quadrícula, podendo utilizar esta informação como X, Y e representa minhas quadrículas de forma contínua.
De fato, meu interesse é mais metodológico que tecnológico, especificamente: modelo de dados pontuais ou de áreas. O valor agregado e representado geograficamente pelas coordenadas de um centroide (x,y) não é teoricamente a realização do processo espacialmente contínuo (geoestatístico) que você esta tentando estimar. Ao fazer agregação você desconsidera a variabilidade intra pixel. Porém há situações que que não dispomos dos dados pontualmente referenciados, mas agregados por área, mesmo num contexto em que faz mais sentido considerar modelo de superfície contínua (geoestatístico). Por favor veja: Correlating point-referenced radon and areal uranium data arising from a common spatial process. B. J. Smith and M. K. Cowles Appl. Statist., 56 (2007), 313 - 326 Se você deseja trabalhar ao nível de área, porque não usar modelos para dados de áreas? É computacionalmente conveniente. Veja o exemplo com pixels da página 11 em: http://www.maths.lth.se/matstat/staff/finn/lindgren_toledo_part123a.pdf Nessa apresentação acima, apresenta-se uma técnica de como usar modelos de dados de área para processos contínuos. Mas a solução não é agregar os dados em áreas.
participantes (2)
-
Alex
-
Elias Krainski