Calculo de áreas

[Guilherme Rosa]

Olá tenho um arquivo Raster contendo as classes de uso e ocupação do solo, originalmente esse arquivo era para o Brasil todo, então recortei para ficar apenas com o estado do Rio Grande do Sul, consigo calcular as áreas desse arquivo tranquilamente usando o r.report. Contudo estou enfrentando um desafio que não consegui encontrar solução ainda.
Se eu recortar consigo por exemplo calcular as áreas de uso e ocupação para cada um dos municípios de forma individual, mas eu quero fazer isso para todos os 497 municípios, e eu quero fazer isso para anos diferentes também logo fica inviável ficar recortando e calculando de forma individual cada um dos municípios, existe uma forma mais fácil de fazer o calculo para múltiplos municípios e obter o resultado das áreas de classes para cada um dos municípios ?   

[Rodrigo Gabriel Bastos]

Não sei se compreendi muito bem, mas parece que você buscar realizar um Histograma Zonal.

Faça assim:

Carregue o raster de uso do solo e o shapefile (vetor) dos municípios.
Abra a Caixa de Ferramentas de Processamento e procure por Histograma Zonal (Zonal Histogram). Não use "Estatísticas Zonais" comum, pois ela calcula média/máximo/mínimo, o que não serve para classes categóricas.
Camada Raster: Seu raster de uso do solo.
Camada Vetorial: Seu shapefile de municípios.

Assim o algoritmo vai criar uma nova camada onde cada classe de uso do solo vira uma coluna.
Se você tem classe 1 e 2, a tabela de atributos dos municípios ganhará colunas chamadas HISTO_1, HISTO_2, etc., contendo o número de pixels ou área daquela classe dentro daquele município.