sexta-feira, 9 de novembro de 2007

ArcGis : Eliminando pequenas regiões em Grids

Apresentaremos a seguir uma rotina que usa a calculadora raster do ArcGis para eliminar regiões no grid cujo tamanho seja menor que a área mínima mapeável. As regiões eliminadas são substituídas pelo vizinho com maior área de contato com o agrupamento de células eliminadas.

Use os comandos abaixo no raster calculator no menu Spatial analyst


regiongroup(ingrid1, #, EIGHT, WITHIN)

Select([Calculation], 'count > XX')

Nibble( [ingrid1], [Calculation2])


onde :

Ingrid1 é o nome do grid de entrada

count > XX ---> limite abaixo do qual as células são eliminadas e trocadas pela classe vizinha.


Exemplo : grid de classificação de uso e cobertura vindo de uma imagem CBERS com célula de 20m.

deixar regiões com área mínima de 1 hectare.

1 hectare = 10.000 m2

1 célula = 20 x 20 = 400 m2

numero de células = 10.000/400 = 25 células

então substitua o XX por 25 ===> Select([Calculation], 'count > 25')


Observe que o comando region group só funciona em grids do tipo inteiro. Se o grid que você quer limpar estiver em outro formato converta-o para inteiro usando o comando

InGrid1 = Int(Meu_Grid)




Nenhum comentário:

Postar um comentário