r.seriesgrass - クラウドでオンライン

プログラム：

NAME

r.シリーズ - 各出力セル値を、セルに割り当てられた値の関数にします。
入力ラスター マップ レイヤー内の対応するセル。

KEYWORDS

ラスター、集約、シリーズ

SYNOPSIS

r.シリーズ
r.シリーズ - 助けて
r.シリーズ [-nz] [ =名[,名、...]] [file=名] 出力=名[,名、...]
方法=string[,string、...] [分位=フロート[,フロート、...]] [重み=フロート[,フロート、...]]
[範囲=こんにちは、こんにちは] [-上書きする] [-助けます] [-詳細] [-静かな] [-ui]

フラグ：
-n
NULL の伝播

-z
ファイルを開いたままにしない

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

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

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

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

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

パラメーター：
=名前[、名前、...]
入力ラスターマップの名前

file=名
XNUMX つのラスター マップ名と、オプションで行ごとに XNUMX つの重み、フィールド セパレーターを含む入力ファイル
名前と体重の間は |

出力=名前[、名前、...] [必要]
出力ラスターマップの名前

方法=文字列[、文字列、...] [必要]
集計操作
オプション： 平均、 カウント、 中央値、 モード、 最小、 min_raster、 最大、 max_raster、
stddev、 範囲、 和、 分散、 多様性、 スロープ、 オフセット、 detcoeff、 値、 quart1、
quart3、 perc90、 分位、 歪度、 尖度

分位=float [、float、...]
Method=quantile で計算する分位数
オプション： 0.0-1.0

重み=float [、float、...]
各入力マップの重み係数。デフォルト値は各入力マップの 1.0 です。

範囲=こんにちは、こんにちは
この範囲外の値を無視する

DESCRIPTION

r.シリーズ 各出力セル値を、に割り当てられた値の関数にします。
入力ラスター マップ レイヤー内の対応するセル。

次の方法を使用できます。

・平均：平均値

· count: 非 NULL セルの数

・median：中央値

· モード: 最も頻繁に発生する値

・最小値：最低値

· maximum: 最高値

· 範囲: 値の範囲 (最大 - 最小)

· stddev: 標準偏差

· sum: 値の合計

· 分散: 統計分散

· 多様性: 異なる値の数

・傾き：線形回帰の傾き

· オフセット: 線形回帰オフセット

· detcoeff: 決定の線形回帰係数

· tvalue: 線形回帰の t 値

· min_raster: 時系列値が最小のラスタ マップ番号

· max_raster: 時系列値が最大のラスタ マップ番号
ほとんどのパラメーターは複数の回答を受け入れるため、複数の集計を許可することに注意してください。
XNUMX 回の実行で計算されます。例:

r.series 入力=map1,...,mapN \
出力=map.mean,map.stddev \
メソッド=平均、標準偏差
または：

r.series 入力=map1,...,mapN \
出力=map.p10,map.p50,map.p90 \
method=分位,分位,分位\
分位数=0.1,0.5,0.9
すべてのオプションに同じ数の値を指定する必要があります。

注意事項

データなし （ヌル） 取り扱い
連絡先 -n フラグ、対応する入力セルのいずれかが NULL であるセルは、
自動的に NULL に設定されます (NULL 伝搬)。 集計関数は呼び出されないため、
すべてのメソッドは、 -n フラグ。

無し -n フラグ、各セル (NULL を含む) の入力の完全なリストが渡されます
集計関数。 個々の集計は、選択したとおりにデータを処理できます。 ほとんどの場合、彼らは
非 NULL 値の集計を計算し、すべての場合にのみ NULL 結果を生成します。
入力は NULL です。

最小 と 分析
この min_raster と 最大ラスター メソッドは、ラスターマップの番号でマップを生成します
時系列の最小値/最大値を保持します。 番号付けは から始まります 0 以下 n
にリストされている最初と最後のラスター 入力=それぞれ。

レンジ 分析
Status 範囲= オプションが指定されている場合、その範囲外の値はすべてとして扱われます
NULL の場合。 の 範囲 パラメータは次のように設定できます 低、高 しきい値: 外側の値
この範囲の
-n が指定された場合、結果は NULL になります)。 の 低、高 しきい値は浮動小数点です。
だから使用する -inf or infファイル 単一のしきい値 (例: 範囲=0,inf 負の値を無視する、または
範囲=-inf、-200.4 -200.4 を超える値を無視するには)。

線形 回帰
線形回帰 (傾き、オフセット、決定係数、t 値) は等しいと仮定します。
時間間隔。 データの時間間隔が不規則な場合、NULL ラスター マップは
時間間隔を等しくするために時系列に挿入されます (例を参照)。

四分位数
r.シリーズ 任意の分位数を計算できます。

メモリ 消費
メモリ使用量は問題ではありません。 r.シリーズ 一度に各マップから XNUMX つの行のみを保持する必要があります。
時間。

マネジメント of 開いた file 制限
処理されるラスター マップの数は、オペレーティング システムの制限によって決まります。 のために
たとえば、通常、ハード制限とソフト制限は両方とも 1024 です。ソフト制限は変更できます。
たとえば、ulimit -n 1500 (UNIX ベースのオペレーティング システム) を使用しますが、ハードウェアよりも高くはありません。
限界。 値が低すぎる場合は、スーパーユーザーとしてエントリを追加できます。
/etc/security/limits.conf
#
your_username ハード nofile 1500
これにより、ハード制限が 1500 ファイルに引き上げられます。 開くファイルが増えるとさらに多くのファイルが必要になることに注意してください
ラム。 Wiki ページ「大規模ラスター データ処理のヒント」も参照してください。

各マップの重み係数は、 重み オプション。 ウェイトの使用
時間範囲が異なるマップの合計または平均を計算するときに意味があります。
デフォルトの重みは 1.0 です。 重みの数は、の数と同じでなければなりません
入力マップと同じ順序である必要があります。 入力ファイルで重みを指定することもできます。

file 開いているファイルにアクセスせずに大量のラスター マップを分析するオプション
コマンドライン引数の制限とサイズ制限。 計算は
オプション方式。 出力マップの単一の行ごとに、すべての入力マップが開かれます
そして閉じました。 RAM の量は、指定された入力マップの数に比例して増加します。
入力オプションとファイル オプションは相互に排他的です。 入力は、改行を含むテキスト ファイルです。
ラスター マップ名とオプションの重みの分離リスト。 マップ名の区切りとして
および重み文字「|」 使用する必要があります。

例

使い方 r.シリーズ ワイルドカードを使用:
r.series input="`g.list pattern='insitu_data.*' sep=,`" \
出力=insitu_data.stddev メソッド=stddev

注意してください g.リスト スクリプトは、マップ名を選択するための正規表現もサポートしています。

使い方 r.シリーズ NULL ラスタ マップの場合 (「完全な」時系列を考慮するため):
r.mapcalc "ダミー = null()"
r.series in=map2001,map2002,dummy,dummy,map2005,map2006,dummy,map2008 \
out=res_slope,res_offset,res_coeff meth=slope,offset,detcoeff

3 回の実行で複数の集計を計算する例 (XNUMX 回の実行で XNUMX つの集計が生成される)
入力マップ):
r.series in=one,two out=result_avg,res_slope,result_count meth=sum,slope,count

r.series の file オプションの使用例:
猫 > input.txt << EOF
map1
map2
map3
EOF
r.series file=input.txt out=result_sum meth=sum

重みを含む r.series のファイル オプションを使用する例。 重み 0.75 は
map2 に割り当てられます。 他のマップには重みがないため、省略できます。
猫 > input.txt << EOF
map1
マップ2|0.75
map3
EOF
r.series file=input.txt out=result_sum meth=sum

XNUMX 日平均を使用して特定の気温を上回った日数をカウントする例
マップ (DOY ワイルドカードとして「???」):
# シェルベースのシステムへのアプローチ
r.series input=`g.list rast pattern="temp_2003_???_avg" sep=,` \
出力=temp_2003_days_over_25度範囲=25.0,100.0メソッド=カウント
# XNUMX つのステップでアプローチする (例: Windows システムの場合)
g.list rast pattern="temp_2003_???_avg" output=mapnames.txt
r.series ファイル=mapnames.txt \
出力=temp_2003_days_over_25度範囲=25.0,100.0メソッド=カウント

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