r.slope.aspectgrass - Online na nuvem

PROGRAMA:

NOME

r.inclinação.aspecto - Gera mapas raster de inclinação, aspecto, curvaturas e parcial
derivados de um mapa raster de elevação.
O aspecto é calculado no sentido anti-horário a partir do leste.

CHAVES

raster, terreno, aspecto, declive, curvatura

SINOPSE

r.inclinação.aspecto
r.inclinação.aspecto --Socorro
r.inclinação.aspecto [-a] elevação=nome [declive=nome] [aspecto=nome] [formato=corda]
[precisão=corda] [pcurvatura=nome] [curvatura=nome] [dx=nome] [dy=nome]
[dxx=nome] [tingido=nome] [dxy=nome] [escala z=flutuar] [inclinação_mínima=flutuar] [-substituir]
[-ajudar] [-detalhado] [-calma] [-ui]

Sinalizadores:
-a
Não alinhe a região atual com o mapa de elevação raster

- sobrescrever
Permitir que os arquivos de saída substituam os arquivos existentes

--Socorro
Imprimir resumo de uso

--verbose
Saída detalhada do módulo

--quieto
Saída silenciosa do módulo

--ui
Forçar o lançamento da caixa de diálogo da GUI

parâmetros:
elevação=nome [obrigatório]
Nome do mapa raster de elevação de entrada

declive=nome
Nome para o mapa raster de declive de saída

aspecto=nome
Nome para mapa raster de aspecto de saída

formato=corda
Formato para relatar a inclinação
opções: graus, por cento
Padrão: graus

precisão=corda
Tipo de aspecto de saída e mapas de declive
opções: CÉLULA, FCELL, Célula DC
Padrão: FCELL

pcurvatura=nome
Nome para o mapa raster de curvatura do perfil de saída

curvatura=nome
Nome para o mapa raster de curvatura tangencial de saída

dx=nome
Nome para o mapa raster de saída de derivada parcial de primeira ordem dx (inclinação EW)

dy=nome
Nome para o mapa raster da derivada parcial de primeira ordem dy (inclinação NS)

dxx=nome
Nome para o mapa raster dxx derivado parcial de segunda ordem de saída

tingido=nome
Nome para o mapa raster dyy derivado parcial de segunda ordem de saída

dxy=nome
Nome para o mapa dxy raster derivado parcial de segunda ordem de saída

escala z=flutuar
Fator multiplicativo para converter unidades de elevação em unidades horizontais
Padrão: 1.0

inclinação_mínima=flutuar
Valor mínimo de inclinação (em porcentagem) para o qual o aspecto é calculado
Padrão: 0.0

DESCRIÇÃO

r.inclinação.aspecto gera mapas raster de declive, aspecto, curvaturas e primeiro e segundo
solicite derivados parciais de um mapa raster de valores reais de elevação. O usuário deve
especifique a entrada elevação mapa raster e pelo menos um mapa raster de saída. O usuário pode
também especifique o formato para inclinação (graus, porcentagem; padrão = graus), e o escala z:
fator multiplicativo para converter unidades de elevação em unidades horizontais; (padrão 1.0).

O elevação o mapa raster de entrada especificado pelo usuário deve conter valores verdadeiros de elevação,
não dados redimensionados ou categorizados. Se os valores de elevação estiverem em outras unidades que não no
unidades horizontais, eles devem ser convertidos em unidades horizontais usando o parâmetro escala z.
In GRAMA GIS 7, vertical unidades e guarante que os mesmos estão não assumiu para be metros qualquer mais. Para a exemplo, if
ambos sua vertical e horizontal unidades e guarante que os mesmos estão pé, parâmetro escala z não deve ser usado.

O aspecto o mapa raster de saída indica a direção para a qual as encostas estão voltadas. O aspecto
categorias representam o número de graus de leste. Arquivos de categoria e tabela de cores também são
gerado para o mapa raster de aspecto. As categorias de aspecto representam os graus de número de
leste e aumentam no sentido anti-horário: 90 graus é o norte, 180 é o oeste, 270 é o sul
360 é o leste.
Observação: esses valores podem ser transformados em valores de azimute (0 é o norte, 90 é o leste, etc.)
usando r.mapcalc:
# converter ângulos de CCW para norte acima
r.mapcalc "azimuth_aspect = (450 - ccw_aspect)% 360"

O aspecto não está definido para declive igual a zero. Assim, a maioria das células com uma pequena
inclinação acabam tendo categoria 0, 45, ..., 360 pol. aspecto saída. É possível reduzir
o viés nessas direções, filtrando o aspecto em áreas onde o terreno é
quase plano. Uma opção inclinação_mínima pode ser usado para especificar a inclinação mínima para qual aspecto
é calculado. O aspecto para todas as células com inclinação inclinação_mínima está definido para nulo (sem dados).

O declive o mapa raster de saída contém valores de inclinação, declarados em graus de inclinação de
a horizontal se formato= opção de graus (o padrão) é escolhida, e em aumento percentual se
formato= opção de porcentagem é escolhida. Arquivos de categoria e tabela de cores são gerados.

As curvaturas de perfil e tangencial são as curvaturas na direção da inclinação mais acentuada
e na direção da tangente do contorno respectivamente. As curvaturas são expressas como
1 / metros, por exemplo, uma curvatura de 0.05 corresponde a um raio de curvatura de 20m. Convexo
os valores da forma são positivos e os valores da forma côncava são negativos.

Exemplo DEM

Inclinação (grau) do exemplo DEM Aspecto (grau) do exemplo DEM

Curvatura tangencial (m-1) do exemplo DEM Curvatura do perfil (m-1) do exemplo DEM

Para algumas aplicações, o usuário desejará usar um mapa raster reclassificado de declive que
agrupa valores de inclinação em intervalos de inclinação. Isso pode ser feito usando r.reclassificar. Um exemplo de
uma reclassificação útil é fornecida abaixo:
categoria faixa categoria rótulos
(em graus) (em porcentagem)
1 0- 1 0- 2%
2 2- 3 3- 5%
3 4- 5 6- 10%
4 6- 8 11- 15%
5 9- 11 16- 20%
6 12- 14 21- 25%
7 15- 90 26% e superior
A seguinte tabela de cores funciona bem com o anterior
reclassificação.
categoria vermelho verde azul
0 179 179 179
1 0 102 0
2 0 153 0
3 128 153 0
4 204 179 0
5 128 51 51
6 255 0 0
7 0 0 0

NOTAS

Para garantir que o mapa de elevação raster não seja reamostrado de maneira inadequada, as configurações para
a região atual é ligeiramente modificada (apenas para a execução do programa): o
resolução é definida para coincidir com a resolução do mapa raster de elevação e as bordas do
região (ou seja, o norte, sul, leste e oeste) são deslocados, se necessário, para alinhar ao longo
bordas das células mais próximas no mapa de elevação. Se o usuário realmente deseja o raster
mapa de elevação reamostrado para a resolução da região atual, o -a sinalizador deve ser especificado.

A máscara atual é ignorada.

O algoritmo usado para determinar a inclinação e o aspecto usa uma vizinhança 3x3 em torno de cada célula
no mapa de elevação raster. Assim, não é possível determinar inclinação e aspecto para
as células adjacentes às bordas na camada do mapa de elevação. Essas células são atribuídas a um
valor de "inclinação zero" (categoria 0) nos mapas de inclinação e varredura de aspecto.

A fórmula de Horn é usada para encontrar as derivadas de primeira ordem nas direções xey.

Somente ao usar modelos de elevação inteiros, o aspecto é tendencioso em 0, 45, 90, 180, 225,
Direções 270, 315 e 360; ou seja, a distribuição das categorias de aspectos é muito desigual,
com picos em 0, 45, ..., 360 categorias. Ao trabalhar com elevação de ponto flutuante
modelos, esse viés de aspecto não ocorre.

EXEMPLOS

Cálculo of declive, aspecto, profile e tangencial curvatura
Neste exemplo, um declive, aspecto, perfil e mapa de curvatura tangencial são calculados a partir de um
mapa raster de elevação (conjunto de dados de amostra da Carolina do Norte):
g.region raster = elevação
r.slope.aspect elevação = elevação declive = declive aspecto = aspecto pcurvatura = pcurv tcurvatura = tcurv
# definir tabelas de cores agradáveis ​​para mapas raster de saída
r.colors -n map = declive cor = sépia
r.colors map = aspect color = aspectcolr
r.colors map = pcurv color = curvature
r.colors map = tcurv color = curvature

Figura: Mapa raster de inclinação, aspecto, perfil e curvatura tangencial (Carolina do Norte
conjunto de dados)

Classificação of principal aspecto direções in bússola orientação
No exemplo a seguir (com base no conjunto de dados de amostra da Carolina do Norte), primeiro geramos
o mapa de aspecto padrão (orientado CCW do leste), em seguida, converta-o para a orientação da bússola,
e, finalmente, classificar quatro direções principais de aspecto (N, E, S, W):
g.region raster = elevação -p
# gerar mapa de aspecto com orientação CCW
r.slope.aspect elevation = elevation aspect = myaspect
# gerar a orientação da bússola e classificar as quatro direções principais (N, E, S, W)
r.mapcalc "aspect_4_directions = eval (\\
bússola = (450 - meuaspect)% 360, \\
if (compass> = 0. && compass <45., 1) \\
+ if (bússola> = 45. && bússola <135., 2) \\
+ if (bússola> = 135. && bússola <225., 3) \\
+ if (bússola> = 225. && bússola <315., 4) \\
+ if (bússola> = 315., 1) \\
)"
# atribuir rótulos de texto
r.category aspect_4_directions separator = regras de vírgula = - << EOF
1, norte
2, leste
3, sul
4, oeste
EOF
# atribuir tabela de cores
r.colors aspect_4_directions rules = - << EOF
1 253,184,99
2 178,171,210
3 230,97,1
4 94,60,153
EOF
Mapa de aspecto classificado em quatro direções principais da bússola (subconjunto ampliado mostrado)

REFERÊNCIAS

· Horn, BKP (1981). Colina Sombreamento e da reflectância Mapa, Procedimentos do
IEEE, 69(1): 14-47.

· Mitasova, H. (1985). Cartográfico aspectos of computador à superfície dos talhos, modelagem. PhD
tese. Universidade Técnica Eslovaca, Bratislava

· Hofierka, J., Mitasova, H., Neteler, M., 2009. Geomorfometria in GRAMA SIG. Dentro:
Hengl, T. e Reuter, HI (Eds), Geomorfometria: Conceitos, Programas,
Aplicações. Desenvolvimentos em Ciência do Solo, vol. 33, Elsevier, 387-410 pp,
http://www.geomorphometry.org

Use r.slope.aspectgrass online usando serviços onworks.net