Combinando raspadores de colina de alta resolução e generalizados para produzir alívio sombreado aprimorado?

Costumo adicionar dados LiDAR a um mapa junto com um alívio sombreado. Na maioria das vezes, tenho que generalizar o LiDAR DEM, dependendo da escala do mapa. Ainda assim, muitas vezes há muitos detalhes e os cumes e vales das montanhas não são bem pronunciados.

Encontrei recentemente um artigo que descreve o processo de combinar dois modelos de colina, um detalhado e o outro generalizado. O alívio resultante é muito bom. Eu gostaria de produzir esse efeito em meus mapas: supondo que eu tenha 2 rasters de colina (detalhados e generalizados), qual é a melhor maneira de combiná-los para obter esse efeito?

— Jakub Sisak GeoGraphics
fonte

Respostas:

Parece o trabalho de Tom Patterson sobre o Resolution Bumping GTOPO 30 no Photoshop . A teoria é descrita suficientemente bem para ser adaptável a outros softwares, embora seja necessário trabalhar com as especificações. A idéia básica é generalizar (desfocar) um conjunto de dados, muito, para enfatizar a forma geral e ocultar detalhes específicos e, em seguida, misturar as imagens de alta e baixa resolução usando canais alfa ou porcentagens de opacidade.

Enquanto o processo descrito por Patterson é aplicado aos próprios modelos de elevação, adaptei-o aos relevos da seguinte maneira.

Carregue a imagem de relevo no GIMP:

duplicar Camada
camada superior == a camada de detalhes, deixe em branco, defina Mode como Overlay
bottom == a layer da shape, aplique Gaussian blur com a configuração 20
salve uma cópia (relief_blurred.tif)

Reatribuir a projeção:

gdalcopyproj shaded_relief.tif relief_blurred.tif

Há muito espaço para adaptação na escolha de quais porcentagens usar, se deve usar os filtros Overlay, Screen ou Dark, etc.

— Matt Wilson
fonte

método muito bom

— Mapperz

ohhh eu vou tentar isso, parece legal.

— Nathan W

Eu estava esperando por uma solução ArcGIS, mas vou tentar isso também ...

— Jakub Sisak GeoGraphics

Estou começando a pensar que a única maneira de fazer isso deve ser como você sugeriu em um editor de imagens como GIMP, Photoshop etc. (Como a rampa de cores do ArcGIS não pode ser definida de transparente para preto, mas sempre o padrão é branco para preto)

— Jakub Sisak GeoGraphics

Vejo que "ArcGIS" é uma tag, Jakub. Usando o Spatial Analyst, você simplesmente calcularia uma média ponderada das duas colinas. Por exemplo, a mistura 60-40 pode ser gerada com um cálculo como este:

(60*[Detailed hillshade] + 40*[Generalized hillshade]) / 100

Se você precisar, o desfoque gaussiano pode ser executado executando alguns meios focais circulares sobre uma varredura. O raio efetivo de fazer m meios focais, cada um dos raios r , é r * Sqrt ( m ).

— whuber
fonte

Obrigado Bill. Minha entrada de colina detalhada é um tamanho de célula 5, 5, mas a colina generalizada é 20, 20. O resultado da calculadora raster é, portanto, 20, 20 também. Parece que aprimora as sombras primárias, mas os detalhes são perdidos. O que você recomendaria para manter o tamanho da célula original? Reamostrar a colina generalizada? Qual técnica de reassentamento usar?

— Jakub Sisak GeoGraphics

Eu tentei, tudo funciona, mas o produto final Hillshade não é muito diferente dos outros que um pouco mais suave originais detalhado ... (Em vez de intensificar a Hillshade detalhada original)

— Jakub Sisak GeoGraphics

a idéia central, embora não dita, na resposta de Bill é que praticamente tudo o que você pode fazer no Gimp / Photoshop / etc. no que diz respeito à mistura de camadas e sobreposições, pode ser feito com cálculo de varredura - desde que você saiba a matemática. Não é fácil, e é por isso que não existem muitas soluções prontas. Seguindo apenas a teoria, não sou especialista em matemática, acho que o que está faltando no exemplo acima é uma terceira varredura para atuar como um controle deslizante de opacidade / transparência de um pobre homem.

— mate Wilkie

@ Jakub: Você está certo, você precisa manter o tamanho da célula mais fino, ou então tudo está perdido. Nos produtos ESRI, isso geralmente é feito configurando preliminarmente alguma variável "ambiente" em algum lugar para especificar o tamanho, extensão e projeção da célula de saída. Procure um menu "opções" do Spatial Analyst ou um menu "ambiente" do ArcToolBox no ArcGIS.

— whuber

@ Matt: Eu posso ter perdido algo na descrição de Tom Patterson. (É difícil de seguir, porque ele diz apenas quais botões apertar, não o porquê ou o que eles estão realmente fazendo.) Você se importaria de elaborar um pouco o que o "controle deslizante opacidade / transparência" está fazendo por nós?

— whuber