これは、Ubuntu Online、Fedora Online、Windows オンライン エミュレーター、MAC OS オンライン エミュレーターなどの複数の無料オンライン ワークステーションの XNUMX つを使用して、OnWorks 無料ホスティング プロバイダーで実行できるコマンド r.terraflowgrass です。
プログラム:
NAME
r.テラフロー - 大規模なグリッドの流れ計算を実行します。
フロートバージョン。
KEYWORDS
ラスター、水文学、流れ、蓄積、シンク
SYNOPSIS
r.テラフロー
r.テラフロー - 助けて
r.テラフロー [-s] 標高=名 埋め=名 方向=名 流域=名
累積=名 tci=名 [d8カット=フロート] [メモリ=整数] [ディレクトリにジョブを開始します。=string]
[統計情報=string] [-上書きする] [-助けます] [-詳細] [-静かな] [-ui]
フラグ:
-s
SFD (D8) フロー (デフォルトは MFD)
SFD: 単一の流れ方向、MFD: 複数の流れ方向
-上書き
出力ファイルが既存のファイルを上書きできるようにする
- 助けて
使用状況の概要を印刷する
-詳細
冗長モジュール出力
- 静かな
静かなモジュール出力
--ui
GUIダイアログを強制的に起動する
パラメーター:
標高=名 [必要]
入力標高ラスターマップの名前
埋め=名 [必要]
出力の塗りつぶし (塗りつぶし) 標高ラスター マップの名前
方向=名 [必要]
出力流れ方向ラスター マップの名前
流域=名 [必要]
出力流域-流域ラスター マップの名前
累積=名 [必要]
出力流量累積ラスター マップの名前
tci=名 [必要]
出力地形収束指数 (tci) ラスター マップの名前
d8カット=フロート
SFD (D8) 方向を使用したルーティング
フローの蓄積がこの値より大きい場合、SFD (D8) 方向を使用してルーティングされます。
(MFD フローのみに意味があります)。 答えが与えられない場合は、デフォルトで無限大が設定されます。
メモリ=整数
使用する最大メモリ(MB単位)
デフォルト: 300
ディレクトリにジョブを開始します。=string
一時ファイルを保持するディレクトリ (サイズが大きくなる可能性があります)
統計情報=string
実行時統計を含むファイルの名前
DESCRIPTION
r.テラフロー ラスター数値標高モデル (DEM) を入力として受け取り、流れを計算します
方向ラスターと流量蓄積ラスター、および浸水標高
ラスター、流域-流域ラスター (流域の周りの流域への分割)、および TCI
(地形収束指数) ラスター マップ。
r.テラフロー はよく知られたアプローチを使用してこれらのラスターを計算しますが、次の点が異なります。
モデリングではなく、アルゴリズムの計算の複雑さに重点が置かれています。
現実的な流れ。 r.テラフロー スケーラブルなソフトウェアの必要性から生まれました。
非常に大きな地形を効率的に処理します。 理論的に最適な値に基づいています
I/O 効率の高いアルゴリズムのフレームワークで開発されたアルゴリズム。 r.テラフロー ました
特に大規模なグリッド向けに設計および最適化されており、地形を処理できます。
他の GIS システムに存在する同様の機能では実用的ではありませんでした。
流れの方向は、MFD (Multiple Flow Direction) モデルまたは
SFD (単一流れ方向、または D8) モデルを以下に示します。 どちらの方法でも計算します
現在のセルの周囲の 3 行 3 列のウィンドウを検査することで、下降勾配の流れの方向を確認します。 SFD
この方法では、最も急な下り勾配の隣接する方向への一意の流れ方向が割り当てられます。 複合機
この方法では、すべての下り勾配の隣接する方向に複数の流れの方向を割り当てます。
最急下降勾配隣接 (SFD) への流れの方向。 すべての下降勾配隣接 (MFD) への流れの方向。
SFD および MFD メソッドは、同じ値を持つセルの流れ方向を計算できません。
すべての隣接セル (平坦なエリア) または下降勾配の隣接セルを持たないセルの高さ
(XNUMXセルピット)。
・台地(はみ出る平らな場所) r.テラフロー ルートがグローバルに流れるようにする
流れはプラトーの流出セルに向かって進みます。
· シンク上 (XNUMX セルのピットを含む、こぼれない平らな場所) r.テラフロー
すべてのシンクが満たされるまで地形を浸水させ、割り当てを行うことで流れを割り当てます。
塗りつぶされた地形上の流れの方向。
地形を水浸しにするために、 r.テラフロー すべてのシンクを特定し、地形を分割します
シンク流域へ (シンク流域には、そのシンクに流入するすべてのセルが含まれます)、
流域の隣接情報を表すグラフを構築し、
この流域-流域グラフは、最小公倍数に沿って流域を互いに結合します。
すべての流域が地形の外側に流路を持つようになるまで、境界を設定します。 洪水により、
すべてのセルが外側に続く下り勾配の流路を持つ、シンクのない地形。
地形、したがって地形内のすべてのセルに SFD/MFD フロー方向を割り当てることができます。
上記。
地形内のすべてのセルに対して流れの方向が計算されると、 r.テラフロー フローを計算します
流れの方向を使用して水を送り、その量を追跡することによって蓄積
水は各セルを通って流れます。
セルの累積流量が、 d8カット オプション、次に
このセルのフローは、SFD (D8) モデルを使用して隣接セルにルーティングされます。 このオプション
フロー累積ラスターにのみ影響し、MFD フローに対してのみ意味があります (つまり、
-s フラグは使用されません); このオプションが SFD フローに使用された場合、無視されます。 デフォルト値
of d8カット is 無限大.
r.テラフロー tci ラスター (地形収束インデックス、
蓄積流量と局所的な傾きの比の対数)。
アルゴリズムの詳細については、以下の [1,2,3、XNUMX、XNUMX] を参照してください。
注意事項
が使用するテクニックの XNUMX つ r.テラフロー は時空のトレードオフです。 特に、
I/O 負荷の高い検索を避けるため、 r.テラフロー を計算して操作します
各セルにその 8 つの関連情報が格納される拡張標高ラスター
隣接するセルごとに合計最大 80B。 結果として r.テラフロー 中間で動作します
最大 80N バイトの一時ファイル。N はセルの数 (行 x
標高ラスター内の列)(より正確には 80K バイト、K は列数)
入力標高ラスター内の有効な (データなしではない) セル)。
これらすべての中間一時ファイルは、 STREAM_DIR
オプション。 ノート: STREAM_DIR 最大 2 つを保存するには、十分な空きディスク容量が必要です
80Nバイト。
この メモリ オプションを使用して、モジュールのメイン メモリ (RAM) の最大量を設定できます。
処理中に使用されます。 実際には、 値 金額を過小評価する必要があります
マシン上の利用可能な (空き) メイン メモリの容量。 r.テラフロー 多くても常に使用します
これだけのメモリがあれば、仮想メモリ システム (スワップ スペース) は決して使用されません。 の
デフォルト値は 300 MB です。
この 統計情報 オプションは、統計情報 (stats) を含むファイルの名前を定義します。
実行されます。
r.テラフロー 行と列の数には制限があります (それぞれ最大 32,767)。
で使用される内部型 r.テラフロー 標高を保存するためのメソッドはコンパイル時に定義できます。
デフォルトでは、 r.テラフロー 標高を浮動小数点数として内部に保存するようにコンパイルされます。 他の
必要に応じてユーザーがバージョンを作成できます。
標高を short として内部に保存するコンパイルに関するヒント:
このようなバージョンでは、使用するスペースが少なくなります (セルあたり最大 60B、中間ファイル最大 60N)。
したがって、スペースと時間の効率が向上します。 r.テラフロー フローティングでの使用を目的としています
ポイント ラスター データ (FCELL)、および r.テラフロー (短い) 整数ラスター データ (CELL) を含む
最大標高が短い SHRT_MAX=32767 の値を超えないこと (これは
標高がメートル単位で保存されている場合、地球の地形データに対する制約にはなりません)。
両方 r.テラフロー と r.テラフロー (短い) 入力標高ラスターを操作します。
整数、浮動小数点、倍精度浮動小数点 (CELL、FCELL、DCELL) のいずれか。 入力ラスターの場合、
許可された内部範囲を超える値が含まれています(の略) r.テラフロー (短い),
フロート用 r.テラフロー)、プログラムは警告メッセージを表示して終了します。 それ以外の場合、すべての値が
入力標高ラスターが範囲内にある場合、それらはに変換 (切り捨て) されます。
内部立面タイプ (の略) r.テラフロー (短い)、浮動小数点 r.テラフロー)。 これで
ケースの精度が失われ、人工的な平らな領域が作成される可能性があります。 たとえば、次の場合
r.テラフロー (短い) 浮動小数点ラスター データ (FCELL または DCELL) で使用され、値
標高の部分は短く切り捨てられます。 これにより、人工的に平らな領域が作成される可能性があり、
の出力 r.テラフロー (短い) のものよりも現実的ではないかもしれません r.テラフロー on
浮動小数点ラスターデータ。 の出力 r.テラフロー (短い) と r.テラフロー
整数ラスター データ (CELL マップ) についても同様です。
例
ノースカロライナ州の狭い地域のサンプル データセットの例:
g.region raster=elev_lid792_1m
r.terraflow 標高=elev_lid792_1m 塗りつぶし=elev_lid792_1m_filled \
方向=elev_lid792_1m_方向 swatershed=elev_lid792_1m_swatershed \
累積=elev_lid792_1m_accumulation tci=elev_lid792_1m_tci
Flow 累積
スピアフィッシュのサンプル データ セット:
g.region raster=elevation.10m -p
r.terraflow elev=elevation.10mfilled=elevation10m.filled \
dir=elevation10m.mfdir swatershed=elevation10m.watershed \
累積=標高10m.accu tci=標高10m.tci
g.region raster=elevation.10m -p
r.terraflow elev=elevation.10mfilled=elevation10m.filled \
dir=elevation10m.mfdir swatershed=elevation10m.watershed \
累積=標高10m.accu tci=標高10m.tci d8cut=500 メモリ=800 \
統計=標高10mstats.txt
onworks.net サービスを使用してオンラインで r.terraflowgrass を使用する