Mais

Suavizar ou reamostrar o DEM para encontrar estradas íngremes?

Suavizar ou reamostrar o DEM para encontrar estradas íngremes?


Estou tentando encontrar trechos de estrada com declives acentuados.

Meu processo é:

  1. puxar uma cadeia de linha dos dados do Openstreetmap
  2. reamostrar para criar um ponto a cada 10 metros
  3. procure os dados de altura para cada ponto do raster ASTER Global Digital Elevation
  4. calcule o gradiente em intervalos de 10m e procure qualquer valor superior a 14%

Isso funciona, mas gera MUITAS anomalias que devem ser devido aos dados de altura conforme eu andei de bicicleta ao longo dessas seções. Uma seção de 10 m em um segmento de estrada local para mim mostra uma diferença de altura de 14 m entre dois pontos separados por 10 m na linha.

Agora, claramente, isso se deve aos dados de altura. É declarado ter uma precisão de:

20 metros com 95% de confiança para dados verticais e 30 metros com 95% de confiança para dados horizontais.

Como devo proceder para eliminar essas anomalias e ter uma melhor estimativa de onde realmente estão os gradientes mais íngremes da estrada?

Tudo foi escrito em PGPLSQL / POSTGIS até agora. Eu tenho um conjunto de cordas de linha com pontos a cada 10m e uma altura para cada ponto.


Como eu baguncei isso nos comentários e não posso editar agora, vou tentar novamente aqui.

Você tem duas opções:

  1. Divida suas estradas da mesma forma que a resolução de seu DEM. Como está agora, você tem várias amostras (10m de distância) para um único pixel (30m x 30m). Quanto mais amostras por pixel, maior será o salto relativo entre os pixels. Pegue 2 pixels, 30m x 30m, diferença de altura de 1m. Com uma amostra no centro de cada pixel, a inclinação é de 1 m de subida / 30 m de execução (0,033 = raso). Agora, divida uma linha imaginária cruzando o limite dos mesmos dois pixels, espaçados de 10m: 1m de subida / 10m de corrida (0,1 = mais íngreme). Diminua o espaçamento para 1m: subida de 1m / corrida de 1m (1,0: íngreme). E assim por diante.
  2. Reamostrar seu DEM para uma resolução artificialmente fina. Isso não vai economizar muito trabalho, pois agora você terá vários exemplos semelhantes. Por exemplo, divida os pixels em # 1 acima em pixels de 10m x 10m e espace as amostras 10m ao longo da linha. Agora você tem três amostras, todas 0,333 subida / 10m corrida (0,033 = raso, igual ao original).

Não posso comentar, então responderei como abordaria isso no ArcMap. Talvez você possa fazer um fluxo de trabalho semelhante.

Use a ferramenta de extração por máscara em seu DEM com suas linhas de estrada como máscara. Agora você tem uma versão raster das estradas com suas elevações.

Use a ferramenta de inclinação para calcular as inclinações ao longo de seu raster de estrada.

Fiz uma rápida tentativa em alguns dados de fluxo e parece funcionar.


Assista o vídeo: Praca z traumą w 4 prostych krokach