Preenchendo lacunas NoData em raster usando o ArcGIS Desktop?

9

Eu tenho uma varredura com lacunas. Quero preenchê-los com médias de células circundantes usando o ArcGIS 10.2 for Desktop.

Eu tentei "Estatísticas Focais", mas calcula a média de CADA célula e não apenas as lacunas.

Eu tentei [Mosaic] .IsNull.Con ([Average], [Mosaic]) encontrado aqui Patching, mas não sei como implementá-lo / ele não funcionará para mim

Eu tentei o Map ALgebra, mas acho que eles estavam usando versões anteriores à 10.2.

O que eu quero. Faça uma varredura falsa com as lacunas preenchidas. Use a varredura falsa para preencher as lacunas da varredura real sem substituir os pontos pelos dados reais.

arcgis-desktop arcgis-10.2 spatial-analyst

— Steven Crisp
fonte

Sintaxe incorreta, use ([mosaic] .isnull) .Con ([average], [mosaic])

— FelixIP 19/02/15

1

Apenas observe que o link 'Patching' parece ter 15 anos (de 2000 a 2001).

— Phloem

2

Muita coisa mudou com os rasters recentemente; portanto, muitos exemplos mais antigos têm a intenção correta, mas apresentam uma sintaxe incorreta para ser usada na versão atual. @crmackey tem a tradução correta, embora eu usasse um círculo e não um retângulo; para dimensões menores, há pouca diferença, mas para áreas (relativamente) grandes da vizinhança, há uma diferença significativa entre círculo e retângulo.

— Michael Stimson

1

Mesmo com 15 anos, ainda é muito relevante. Huber conhecia e conhece o assunto extremamente bem.

— FelixIP

1

Exatamente. A operação é a mesma, apenas a sintaxe foi alterada.

— Michael Stimson

9

Isso funciona:

filled = arcpy.sa.Con(arcpy.sa.IsNull(in_raster),arcpy.sa.FocalStatistics(in_raster,
                        arcpy.sa.NbrRectangle(w, h),'MEAN'), in_raster)

Onde "w" e "h" são o raio da pesquisa.

Isso faz apenas estatísticas focais nas áreas NoData. Eu verifiquei apagando dados de um DEM e depois encontrando a diferença. Você só precisa garantir que o raio da pesquisa seja adequado.

— crmackey
fonte

Acabei de adicionar tudo isso à calculadora raster e alterar o "in_raster" ew eh? Por alguma razão, não consigo fazê-lo funcionar, embora o fluxo faça sentido, pois pula a criação manual da varredura focal e o faz em tempo real.

— Steven Crisp

Tenho que trabalhar! Muito obrigado! foi apenas adicionando "é onde eu não era; t suposto Obrigado.!

— Steven Crisp

1

Sem problemas. Um aviso, como apontou @ MichaelMiles-Stimson, usar NbrCircle () seria melhor.

— 21415 Cracker:

O wh está nas unidades de células do mapa?

— Se você não sabe- apenas GIS

Você seria melhor usando o local arcpy.sa.NbrCircle(radius, units)onde as unidades podem estar 'MAP'ou 'CELL'.

— crmackey

3

Use a função de estatísticas focais, conforme descrito para calcular a média de todas as células da varredura. Em seguida, use as funções Con e IsNull na calculadora raster para substituir apenas as células que são nulas.

Con(IsNull(<inputRaster>), <focalRaster>, <inputRaster>)

Uma desvantagem possível com esta solução é que você pode ver alguns valores médios adicionais indesejados de células em torno do perímetro da sua varredura de saída

— Ben Gosack
fonte

Não consigo votar devido à reputação ou não, mas esse método funciona se eu fizer uma varredura de estatísticas focais primeiro e depois aplicar o código acima. Eu gostaria de poder fazer tudo em uma linha, como no comentário acima, para economizar tempo. Eu tenho um par de centenas destas coisas ...

— Steven Crisp

Depende apenas se você deseja manter suas estatísticas focais rasterizadas ou interpretar on-the-fly ... ambos os métodos são essencialmente os mesmos (+1 para ambos).

— Michael Stimson

3

Você deve usar a Calculadora de varredura e há uma linha de código útil e relevante na página de suporte da ESRI : Como: Remover e substituir nenhum valor de dados dentro de uma varredura usando informações estatísticas dos valores de dados circundantes . Da sua pergunta, parece que isso faria exatamente o que você está tentando fazer, mas deixe-me saber se estou errado. Aqui está minha captura de tela da Calculadora Raster usando uma instrução condicional para calcular um valor médio para células sem dados das células circundantes: insira a descrição da imagem aqui

Se você tiver centenas delas, recomendo usar o Raster Iterator no Model Builder, que alimentará todas as suas rasters na Raster Calculator Tool, desta forma:

insira a descrição da imagem aqui

— Leo
fonte

2

Tentei o método da Calculadora de varredura que LeonB usou, mas ocorreu um erro de falta de memória, o que é surpreendente, considerando as especificações do meu sistema. Escreveu um pequeno script python independente - modificado daqui - no IDLE que funcionou para mim ...

import arcpy
from arcpy import env
from arcpy.sa import *

env.workspace="D:/Watershed/"

inRaster="D:/Watershed/tamtr14_1m"

arcpy.CheckOutExtension ("Spatial")

outCon=Con(IsNull(inRaster), FocalStatistics (inRaster, NbrRectangle (5,5, "CELL"), "MEAN"), inRaster)

outCon.save ("D:/Watershed/outCon")

— wiebe712
fonte