r.slope.aspectgrass - クラウドでオンライン

プログラム：

NAME

r.スロープ.アスペクト - 傾斜、アスペクト、曲率、部分的なラスターマップを生成します。
標高ラスター マップから派生したもの。
アスペクトは東から反時計回りに計算されます。

KEYWORDS

ラスター、地形、アスペクト、傾斜、曲率

SYNOPSIS

r.スロープ.アスペクト
r.スロープ.アスペクト - 助けて
r.スロープ.アスペクト [-a] 標高=名 [スロープ=名] [側面=名] [形式でアーカイブしたプロジェクトを保存します．=string]
[精度=string] [曲率=名] [曲率=名] [dx=名] [dy=名]
[dxx=名] [染料=名] [DXY=名] [zスケール=フロート] [min_slope=フロート] [-上書きする]
[-助けます] [-詳細] [-静かな] [-ui]

フラグ：
-a
現在の領域をラスター標高マップに合わせないでください

-上書き
出力ファイルが既存のファイルを上書きできるようにする

- 助けて
使用状況の概要を印刷する

-詳細
冗長モジュール出力

- 静かな
静かなモジュール出力

--ui
GUIダイアログを強制的に起動する

パラメーター：
標高=名 [必要]
入力標高ラスターマップの名前

スロープ=名
出力傾斜ラスター マップの名前

側面=名
出力アスペクト ラスター マップの名前

形式でアーカイブしたプロジェクトを保存します．=string
傾きを報告するためのフォーマット
オプション： 度、 パーセント
デフォルト： 度

精度=string
出力アスペクトのタイプと傾斜マップ
オプション： 細胞、 FCELL、 DCELL
デフォルト： FCELL

曲率=名
出力プロファイル曲率ラスター マップの名前

曲率=名
出力接線曲率ラスター マップの名前

dx=名
出力一次偏導関数 dx (EW 勾配) ラスター マップの名前

dy=名
出力一次偏導関数 dy (NS スロープ) ラスター マップの名前

dxx=名
出力二次偏導関数 dxx ラスター マップの名前

染料=名
出力二次偏導関数 dyy ラスター マップの名前

DXY=名
出力二次偏微分 dxy ラスター マップの名前

zスケール=フロート
標高単位を水平単位に変換するための乗算係数
デフォルト： 1.0

min_slope=フロート
アスペクトを計算する最小傾斜値 (パーセント単位)
デフォルト： 0.0

DESCRIPTION

r.スロープ.アスペクト 傾斜、アスペクト、曲率、および第 XNUMX と第 XNUMX のラスター マップを生成します。
真の標高値のラスター マップから偏導関数を順序付けします。 ユーザーは次のことを行う必要があります
入力を指定する 標高 ラスター マップと少なくとも XNUMX つの出力ラスター マップ。 ユーザーができることは、
も指定します 形式でアーカイブしたプロジェクトを保存します． 勾配 (度、パーセント、デフォルト = 度) の場合、 zスケール:
標高単位を水平単位に変換する乗算係数。 (デフォルトは1.0)。

この 標高 ユーザーが指定した入力ラスター マップには、実際の標高値が含まれている必要があります。
再スケールまたは分類されたデータ。 標高値が他の単位である場合、
水平単位の場合は、パラメータを使用して水平単位に変換する必要があります。 zスケール.
In GRASS GIS 7, 垂直 ユニット 　 想定される 〜へ be メートル どれか もっと。 例、 if
両言語で 垂直 と 水平な ユニット 　 フィート、 パラメーター zスケール 使用しないでください。

この 側面 出力ラスター マップは、斜面が向いている方向を示します。 側面
カテゴリは東の度数を表します。 カテゴリおよびカラーテーブルファイルも
アスペクト ラスター マップ用に生成されます。 アスペクト カテゴリは次数を表します。
東に進むと反時計回りに増加します。90 度が北、180 度が西、270 度が南です
360は東です。
注: これらの値は方位角 (0 は北、90 は東など) の値に変換できます。
r.mapcalc を使用して:
# CCW から北上に角度を変換します
r.mapcalc "方位角アスペクト = (450 - ccw_アスペクト) % 360"

アスペクトは、傾きが XNUMX に等しい場合には定義されていません。 したがって、非常に小さい細胞のほとんどは、
勾配は最終的にカテゴリ 0、45、...、360 になります。 側面 出力。 減らすことは可能です
地形が存在する領域のアスペクトをフィルタリングしてこれらの方向のバイアスを除去します。
ほぼフラット。 オプション min_slope どのアスペクトの最小傾斜を指定するために使用できます。
が計算されます。 傾きのあるすべてのセルのアスペクト min_slope に設定されています ヌル （データなし）。

この スロープ 出力ラスター マップには、からの傾斜度で示された傾斜値が含まれます。
水平方向の場合 形式でアーカイブしたプロジェクトを保存します．=degrees オプション (デフォルト) が選択されており、次の場合はパーセント上昇します。
形式でアーカイブしたプロジェクトを保存します．=パーセントオプションが選択されています。 カテゴリとカラーテーブルファイルが生成されます。

プロファイル曲率と接線曲率は、最も急な斜面の方向の曲率です。
と輪郭接線の方向にそれぞれ。 曲率は次のように表されます。
1/メートル、たとえば曲率 0.05 は曲率半径 20m に相当します。 凸型
形状値は正、凹形状値は負です。

DEM の例

DEM の例からの傾き (度) DEM の例からのアスペクト (度)

DEM の例からの接線方向の曲率 (m-1) DEM の例からのプロファイル曲率 (m-1)

一部のアプリケーションでは、ユーザーは、再分類された斜面のラスター マップを使用することを希望します。
勾配値を勾配の範囲にグループ化します。 これは次を使用して実行できます r.再クラス化。 の例
有用な再分類を以下に示します。
カテゴリ範囲 カテゴリ ラベル
(度単位) (パーセント単位)
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% 以上
次のカラー テーブルは上記のカラー テーブルとうまく連携します。
再分類。
カテゴリー 赤 緑 青
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

注意事項

ラスター標高マップが不適切にリサンプリングされないようにするには、
現在の領域はわずかに変更されます (プログラムの実行のみ)。
解像度は、標高ラスター マップとそのエッジの解像度と一致するように設定されます。
必要に応じて、地域 (つまり、北、南、東、西) をずらして並べます。
標高マップ内の最も近いセルのエッジ。 ユーザーが本当にラスターを必要とする場合
現在の領域解像度にリサンプリングされた標高マップ、 -a フラグを指定する必要があります。

現在のマスクは無視されます。

傾きとアスペクトを決定するために使用されるアルゴリズムでは、各セルの周囲の 3x3 近傍が使用されます。
ラスター標高マップ内。 したがって、傾斜とアスペクトを決定することはできません。
標高マップ レイヤーのエッジに隣接するセル。 これらのセルには、
勾配ラスター マップとアスペクト ラスター マップの両方の「ゼロ勾配」値 (カテゴリ 0)。

ホーンの公式は、x および y 方向の一次導関数を見つけるために使用されます。

整数標高モデルを使用する場合のみ、アスペクトは 0、45、90、180、225、
270、315、360 方向。 つまり、アスペクト カテゴリの分布は非常に不均一であり、
ピークは 0、45、...、360 カテゴリにあります。 浮動小数点標高を使用する場合
モデルでは、そのようなアスペクトバイアスは発生しません。

例

計算 of スロープ、 側面、 プロフィール と 接線 曲率
この例では、傾斜、アスペクト、プロファイル、および接線曲率マップが次の式から計算されます。
標高ラスター マップ (ノースカロライナ州サンプル データセット):
g.region ラスター=標高
r.slope.aspect 標高=標高 勾配=スロープ アスペクト=アスペクト p曲率=pcurv tcurvature=tcurv
# 出力ラスターマップに適切なカラーテーブルを設定します
r.colors -n マップ=斜面の色=セピア
r.colors マップ=アスペクトカラー=アスペクトカラー
r.colors マップ=pcurv カラー=曲率
r.colors マップ = tcurv カラー = 曲率

図: 傾斜角、傾斜角、断面図、および接線曲率ラスター マップ (ノースカロライナ州)
データセット)

Classification of 主要な 側面 方向 in コンパス オリエンテーション
次の例 (ノースカロライナ州のサンプル データセットに基づく) では、最初に生成します。
標準のアスペクト マップ (東から反時計回り) をコンパスの方向に変換します。
そして最後に XNUMX つの主要なアスペクト方向 (N、E、S、W) を分類します。
g.regionラスター=標高-p
# CCW方向のアスペクトマップを生成
r.slope.aspect 標高=標高アスペクト=myaspect
# コンパスの方向を生成し、XNUMX つの主要な方向 (N、E、S、W) を分類します
r.mapcalc "aspect_4_directions = eval( \\
コンパス=(450 - myaspect ) % 360, \\
if(コンパス >=0. && コンパス < 45., 1) \\
+ if(コンパス >=45. && コンパス < 135., 2) \\
+ if(コンパス >=135. && コンパス < 225., 3) \\
+ if(コンパス >=225. && コンパス < 315., 4) \\
+ if(コンパス >=315., 1) \\
)"
# テキストラベルを割り当てる
r.category アスペクト_4_方向セパレータ=カンマ規則=- << EOF
1、北
2、東
3、南
4、西
EOF
# カラーテーブルを割り当てる
r.colors aspect_4_directions rules=- << EOF
1 253,184,99
2 178,171,210
3 230,97,1
4 94,60,153
EOF
XNUMX つの主要なコンパス方向に分類されたアスペクト マップ (拡大されたサブセットを表示)

参考文献

· ホルン、BKP (1981)。 丘 シェーディング と 　 反射率 地図、議事録
IEEE、 69（ＸＮＵＭＸ）：ＸＮＵＭＸ − ＸＮＵＭＸ。

· Mitasova、H. (1985)。 地図製作 側面 of コンピュータ 表面 モデリング。 博士
論文。 スロバキア工科大学、ブラチスラヴァ

· Hofierka, J.、Mitasova, H.、Neteler, M.、2009 年。 地形計測 in GRASS GIS。 に：
Hengl, T. および Reuter, HI (編著)、 地形計測: コンセプト、 ソフトウェア、
アプリケーション。 土壌科学の発展、vol. 33、エルゼビア、387-410 pp、
http://www.geomorphometry.org

onworks.net サービスを使用してオンラインで r.slope.aspectgrass を使用する