Eu tenho uma varredura de elevação e gostaria de obter algumas variáveis básicas do terreno, como declividade, robustez, aspecto ou curvatura.
O problema é que parece que não tenho a possibilidade de calcular a curvatura com o analista de terreno:
Alguém pode me dizer por que esse indicador não está na lista? Como posso calcular de outra forma?
Respostas:
A curvatura é uma derivada complexa do terreno para calcular, a equação que você usa depende da resolução dos dados de entrada, pois você precisa garantir que os resultados da curvatura que você calcula possam ser diferenciados do ruído nos dados.
Recentemente, muita pesquisa foi feita sobre cálculos de curvatura em dados de alta resolução LiDAR, que mostraram que existe uma quebra de escala em torno de 2 ou 3 metros de resolução e, acima desse ponto, algoritmos diferentes (com os quais não estou familiarizado) precisam ser usados . As melhores informações sobre o cálculo da curvatura topográfica provavelmente vêm de Hurst et al. 2012 e suas referências.
O princípio básico do cálculo da curvatura, como com a inclinação e o aspecto, é passar uma janela em movimento sobre a superfície da elevação e ajustar os valores de elevação a uma função polinomial de 6 termos, cujos coeficientes produzirão a inclinação, o aspecto e a curvatura do centro célula da janela em movimento.
O ArcGIS usa uma janela de pesquisa 3x3 que só produz bons resultados em áreas completamente desprovidas de vegetação, o que torna a ferramenta bastante inútil, a menos que as pessoas estejam cientes dessa limitação, isso pode sugerir por que ela não está presente no QGIS.
A matemática foi derivada originalmente (acho) em Evans (1980) e foi simplificada em algumas páginas em Princípios de sistemas de informação geográfica (link Amazon), que posso recomendar como um bom guia para esse tipo de análise de terreno em um nível básico.
Uma maneira de calcular a curvatura de um DEM é convertê-lo em uma varredura ascii, lê-lo em uma matriz numpy e, em seguida, executar o ajuste polinomial em uma janela móvel que passa pelos dados. Isso é bastante fácil de fazer, mas é muito lento para executar e precisa de uma quantidade razoável de otimização (esse tipo de operação geralmente é portado para c ++ para acelerar).
Para executar a operação no QGIS, você pode usar o plug-in GRASS r.slope.aspect, que também é limitado pela janela fixa 3x3.
Sei que essa não é a resposta simples que você sem dúvida esperava, mas espero que você entenda que a curvatura é complexa para derivar de maneira significativa. Muito bem sucedida.
Evans, I. S. (1980), An integrated system of terrain analysis and slope mapping, Z. Geomorphol., 36, 274–295.
A versão da Raster Analysis da ESRI para calcular a curvatura pode ser útil para desenvolver um plug-in para o QGIS.
Para cada célula, um polinômio de quarta ordem do formulário:
Z = Ax²y² + Bx²y + Cxy² + Dx² + Ey² + Fxy + Gx + Hy + I
é ajustado a uma superfície composta por uma janela 3x3. Os coeficientes a, b, ce etc. são calculados a partir dessa superfície.
As relações entre os coeficientes e os nove valores de elevação para cada célula numerada, como mostrado no diagrama, são as seguintes: Diagrama de valores de curvatura Diagrama de valores de curvatura
A = [(Z1 + Z3 + Z7 + Z9) / 4 - (Z2 + Z4 + Z6 + Z8) / 2 + Z5] / L4
B = [(Z1 + Z3 - Z7 - Z9) / 4 - (Z2 - Z8) / 2] / L3
C = [(-Z1 + Z3 - Z7 + Z9) / 4 + (Z4 - Z6)] / 2] / L3
D = [(Z4 + Z6) / 2 - Z5] / L2
E = [(Z2 + Z8) / 2 - Z5] / L2
F = (-Z1 + Z3 + Z7 - Z9) / 4L2
G = (-Z4 + Z6) / 2L
H = (Z2 - Z8) / 2L
I = Z5
A saída da ferramenta Curvatura é a segunda derivada da superfície - por exemplo, a inclinação da inclinação - de modo que:
Curvatura = -2 (D + E) * 100
Informações completas e fonte:
http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//00q90000000t000000
-2(D+E)*100como (100/L2)*(3*Z5 - [Z2+Z4+Z6+Z8+Z5])mostra como obter este valor como três vezes o valor original, Z5, menos uma soma vizinhança Z2+Z4+Z6+Z8+Z5(usando um círculo de raio 1), tudo por redimensionados 100/L2. São apenas três operações simples da rede.
A curvatura pode ser calculada usando o módulo da SAGA 'Análise de terreno - Morfometria ---> Inclinação, Aspecto, Curvatura'
O cálculo pode ser feito com base em um destes algoritmos:
O LandSerf pode fazer isso. Você pode definir o tamanho das janelas (3 * 3, 5 * 5, 7 * 7, 11 * 11, ...), mas deve ser um número ímpar. análises multiescala. você pode considerar a dependência da escala da inclinação, aspecto e curvatura. http://www.landserf.org/ LandSerf ajustará a equação polinomial co-quadrática a janelas predefinidas específicas, mas se você definir o tamanho de janelas grandes como 50 * 50, leva muito tempo. Depende do tamanho do raster e do tamanho das janelas que você definir. Jo Wood escreveu Landserf para sua tese de doutorado. Está escrito em Java.
Pode ser calculado em diferentes escalas no LandSerf
Você também pode experimentar o SAGA GIS gratuito (http://sourceforge.net/apps/trac/saga-gis/wiki) ou o TAS (http://www.uoguelph.ca/~hydrogeo/TAS/index.html).