Na medida em que os dados vai, eu estou trabalhando com NHD arquivos .shp, dems 10m, e alguns dados LIDAR.
Meu objetivo é determinar o gradiente para segmentos de 100m de uma rede de fluxos.
Eu já sou capaz de fazer isso, mas espero que meu fluxo de trabalho não seja ideal, principalmente porque não consigo lidar com redes ramificadas.
Se todos vocês estivessem fazendo isso, que tipo de etapas você usaria?
Além disso, postei sobre o problema aqui , onde acho que fiz um trabalho muito melhor descrevendo quais são meus objetivos.
arcgis-desktop
dem
lidar
hydrography
Jacques Tardie
fonte
fonte
Respostas:
Dado que você possui o LIDAR DEM, você deve usar os fluxos derivados dele. Isso garante um registro perfeito.
O ponto crucial da idéia é estimar as inclinações médias em termos das elevações nas extremidades dos segmentos.
Um dos procedimentos mais fáceis é "explodir" a rede de fluxo em seus arcos não ramificados de componentes. Converta a coleção em uma camada de "rota" com base na distância, tornando-a "mensurável". Agora é fácil gerar uma coleção de "eventos" de rota com base em uma tabela de marcos (a intervalos de 100 m, por exemplo) para cada arco e extrair as elevações de DEM desses pontos de evento. As sucessivas diferenças de elevação ao longo de cada arco, divididas por 100m, estimam as inclinações do segmento médio.
A figura a seguir mapeia os arcos dos fluxos derivados de uma análise de acumulação de fluxo de um DEG de USGS 7,5 minutos (parte de Highland County, VA). É cerca de 10 km de diâmetro (6 milhas).
Como você está procurando uma barragem remanescente, que pode ser indicada por uma mudança de gradiente em apenas algumas dezenas de metros (para uma barragem muito pequena), considere usar segmentos ainda menores . Se o conjunto de dados for muito áspero para fornecer sinais claros, você poderá filtrá-lo facilmente mais tarde (por meio de médias móveis ou de outra forma, como plotagem das elevações e diferenciação da spline). De fato, essa abordagem coloca você no domínio da análise de séries temporais, onde a variável de interesse é a elevação, não o gradiente, e você procura padrões consistindo em seções de nível curto seguidas por mudanças repentinas.
Este é um gráfico de elevações de DEM observadas em intervalos de 100m ao longo da maioria (não todos) dos segmentos de fluxo representados. (O tamanho da célula é de 30m.) Onde necessário, os arcos eram reorientados para fazer a elevação geralmente diminuir da esquerda para a direita. (Se você olhar com atenção, poderá ver onde eu perdi uma: ela sobe da esquerda para a direita.)
Esse detalhe do arco 16 (o segmento longo na parte superior do mapa) mostra o que você pode obter quando os fluxos não estão perfeitamente registrados no DEM: em lugares em que o fluxo parece fluir para cima. No entanto, os segmentos que sugerem características de pool-and-drop são facilmente identificados, especialmente após os marcos 1800 (metros ao longo do segmento), 4000, 4600 e 6500. Essa identificação pode ser automatizada de várias maneiras, principalmente após a limpeza da série de elevação (suavizando isto).
Você pode ver que o intervalo de amostragem de 100 m usado aqui realmente não é bom o suficiente para identificar recursos muito menores, com 400 a 500 metros de comprimento. Portanto, para encontrar uma pequena barragem remanescente, você provavelmente gostaria de experimentar um intervalo de 10 a 25 m no seu LIDAR DEM.
BTW, o que torna um segmento de fluxo "muito pequeno" para esse tipo de trabalho não é de tamanho curto nem de tamanho de célula grande, embora ambos participem da decisão. "Muito pequeno" depende de como você usará as inclinações estimadas e quão incertas essas estimativas podem ser. Para alguns trabalhos, pode até fazer sentido estimar gradientes em intervalos de 10m em uma grade de 10m!
fonte
Estou fazendo uma análise de hidrologia do meu lado e, como eu era para criar minha varredura de Flow Direction, lembrei-me de seu post. Isso é apenas uma facada no escuro, mas, no ArcGIS 10, há uma opção para criar uma varredura de queda de saída. Gostaria de saber se de alguma forma poderia ser usado para resolver o seu problema.
fonte
A resposta de Jakub é boa, porque considera cada célula sem necessidade de divisão de linha adicional. Se você combinasse um raster de fluxo com a acumulação de fluxo ao longo desse raster, você poderia obter a distância ao longo do fluxo e, em seguida, representar graficamente a inclinação no eixo y e a distância do fluxo no eixo x. Você também precisaria considerar a distância diagonal, mas isso poderia ser resolvido usando a Direção Euclidiana.
fonte