Use os comandos abaixo no raster calculator no menu Spatial analyst no ArcGis
regiongroup(ingrid1, #, EIGHT, WITHIN)
Select([Calculation], 'count > 50')
Nibble( [ingrid1], [Calculation2 ])
onde :
Ingrid1 o nome do grid de entrada, onde deseja-se eliminar as pequenas regiões
count > 50 ---> número de células limite, abaixo do qual as células são eliminadas e trocadas pela classe vizinha, ou seja o menor agrupamento de células que permanecerão no grid
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 5 hectares.
- 5 hectares = 50.000 m2
1 célula = 20 x 20 = 400 m2
numero de células = 50000/400 = 125 células
então substitua o valor 50 por 125 ===> Select([Calculation ], 'count > 125')
Observe que o comando region group só funciona em grids do tipo inteiro. Se o grid que você quer limpar não estiver neste formato converta-o para inteiro usando o comando
InGrid1 = Int(Meu_Grid)
onde :
- Meu_Grid é o nome do seu grid de entrada
- Ingrid é o grid resultante convertido para o formato inteiro
1 comentários:
Não consigo aplicar a formula no ArcGis 10.1, está dando um erro de sintaxe. Consegui transformar os números para inteiros somente. Teve alguma alteração para o ArcGis 10.1 ou tem algum macete para ser feito o procedimento? Obrigado!
Postar um comentário