Mais

Soma os valores em um raster (25m) usando outro raster (1km) ArcGIS Desktop 10.2

Soma os valores em um raster (25m) usando outro raster (1km) ArcGIS Desktop 10.2


Estou resumindo alguns dados populacionais usando uma sobreposição raster de 1km. Converti dados de população, camada poligonal, em um raster de 25 metros que contém o valor de população / m ^ 2. Esse valor é então multiplicado por 625 m ^ 2 (a área de cada célula), então eu descobri que poderia fazer um resumo zonal usando minha camada raster de 1km (que se ajusta perfeitamente à camada raster de 25 metros), mas agora eu estou preso. Se eu pudesse apenas somar os valores de cada célula de 25 metros que ocorre em uma célula de 1 km, então eu só precisaria multiplicar por 625 para obter a população total aproximada nessa célula ... Eu poderia usar uma reamostragem majoritária, mas isso deveria dê uma estimativa melhor.

Recebi ótimos conselhos sobre como codificar uma solução, mas gostaria de descobrir seguindo essa linha de pensamento, se possível.


Se você quiser ficar completamente em raster, seria bom se Resample tivesse um método de Soma, mas contanto que você tenha Spatial Analyst disponível, Zonal Statistics deve ser capaz de resolver seu problema.

Se estiver usando outro raster para suas zonas, cada zona (que será cada célula neste caso) precisa de um valor único. Todos os valores das células de 1km são iguais? Nesse caso, você precisará alterar o raster e / ou adicionar uma tabela de atributos de raster. Alternativamente, em vez de usar um raster para entrada de zona, você pode apenas usar sua rede arrastão (assumindo que sejam polígonos). Veja o link de ajuda de Estatísticas Zonais acima - há várias considerações para o raster ou vetor usado para definição de zona (tamanho da célula, snap, número de bandas, presença de uma tabela de Atributos Raster, etc.).

A saída resultante deve ser um raster de célula de 1km com cada célula a soma dos valores de células raster de 25m. Você pode então executar este raster de saída através da Calculadora de Raster para sua operação [625 x cada célula] para um raster de resumo final.

Se você quiser trabalhar em tabelas, há também a ferramenta Estatísticas zonais como tabela, que cria uma tabela em vez de uma saída raster. Observe que tanto ele quanto a criação de uma tabela de atributos raster já resumem células idênticas como registros únicos na tabela com uma contagem de frequência, não registros de células individuais. Se você quer uma mesa com cada célula como um registro (assumindo que eles ainda não sejam todos valores únicos) você pode consultar Como obter os valores de cada célula na tabela de atributos raster?


Você já tem seus polígonos de sobreposição, agora transforme seu raster de dados em polígonos usando raster para polígono, selecione 'não simplificar' e seus polígonos corresponderão exatamente às suas células. Existem limites para o tamanho de um raster que pode ser processado, pois o arquivo de forma intermediário deve ser inferior a 2 GB para qualquer um dos arquivos; Devido a essas limitações de arquivo, você pode agrupar a poligonização usando extrair por máscara e selecionar metades ou quartos de seus polígonos sobrepostos como a máscara. Essa abordagem também irá acelerar as coisas, pois cada lote pode ser feito simultaneamente com várias instâncias do ArcCatalog.

Antes de realizar a união, seria melhor usar Integrar para minimizar / aliviar ligeiras incompatibilidades entre as camadas do polígono que resultariam em polígonos lascados. Certifique-se de ler os avisos sobre esta ferramenta e mantenha um backup; com parâmetros excessivamente entusiasmados, pode destruir dados confiáveis.

Una cada parte por vez com os recursos de sobreposição; A união é limitada a duas classes de recursos com níveis de licença básicos ou padrão que devem ser adequados para esta operação. Após a união, todos os valores dos dados de sobreposição e raster estarão presentes, para resumir esta estatística de resumo de uso com um campo de caso do (s) atributo (s) importante (s) da sobreposição e soma dos atributos de raster.

Nota: este processo funcionará igualmente bem se você executá-lo em pontos gerados a partir do raster, a menos que você precise de uma área poligonal de sobreposição, caso em que a integração seria mais prejudicial do que benéfica.


Assista o vídeo: Resampling Raster ArcGis. changing the cell size of Raster dataset in ArcGis