t.rast.mapcalcgrass-クラウドでのオンライン

プログラム：

NAME

t.rast.mapcalc -時間的にサンプリングされたマップに対して時空間mapcalc式を実行します
時空ラスターデータセットの。

KEYWORDS

時間、代数、ラスター、時間

SYNOPSIS

t.rast.mapcalc
t.rast.mapcalc - 助けて
t.rast.mapcalc [-ns] 入力=名[,名、...] 表現=string [方法=名[,名、...]]
出力=名 ベース名=ベース名 [nproc=整数] [-上書きする] [-助けます] [-詳細]
[-静かな] [-ui]

フラグ：
-n
ヌルマップを登録する

-s
空間的な重なりを確認する

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

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

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

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

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

パラメーター：
入力=名前[、名前、...] [必要]
入力時空ラスターデータセットの名前

表現=string [必要]
時空間mapcalc式

方法=名前[、名前、...]
入力データセットのサンプリングに使用されるメソッド
オプション： 始めて、 その間、 オーバーラップ、 含む、 同等、 続く、 先行
デフォルト： 等しい

出力=名 [必要]
出力時空ラスターデータセットの名前

ベース名=ベース名 [必要]
出力ラスターマップのベース名
アンダースコアで区切られた数字の接尾辞が付加され、一意の
識別子

nproc=整数
並行して実行するr.mapcalcプロセスの数
デフォルト： 1

DESCRIPTION

t.rast.mapcalc 時空間を実行します マップカルク 時間的にサンプリングされたマップ上の式
時空間ラスターデータセット（STRDS）。 空間および時間演算子と内部変数
式の文字列で使用できます。 空間演算子の説明、
関数と内部変数は、r.mapcalcのマニュアルページで入手できます。 側頭葉
関数については、以下で詳しく説明します。

このモジュールにはいくつかのパラメーターが必要です。 参照されるすべての時空間ラスターデータセット
会場は マップカルク 表現 に記載されている必要があります オプション。 ザ· 最初の 時空ラスター
入力としてリストされているデータセットは、他のすべての時空を一時的にサンプリングするために使用されます
ラスターデータセット。 時間的サンプリング方法は、 方法 オプション。 ザ·
mapcalc式でのSTRDSの順序は、
入力オプションのSTRDS。 結果の時空間ラスターデータセットは、で指定する必要があります
　 出力 オプションと一緒に ベース 登録されている生成されたラスターマップの名前
結果のSTRDSで。 マップ計算の結果として生じる空のマップは、によって登録されていません
ディフォルト。 この動作は、 -n 国旗。 旗 -s 保証するために使用することができます
STRDS内の空間関連マップのみが処理されます。 空間関連とは、
時間的に関連するマップは、空間的な範囲で重複しています。

モジュール t.rast.mapcalc 並列処理をサポートします。 オプション nproc を指定します
並行して開始できるプロセスの数。

時間的にサンプリングされたマップを処理するには、mapcalc式を指定する必要があります。 時間的
内部変数に加えて利用可能です r.mapcalc 空間演算子と
機能：

相対時間と絶対時間でサポートされる内部変数：

· td（） -この内部変数は、現在のサンプル時間のサイズを表します
絶対時間の日数と日数の端数、および相対単位での間隔
相対時間の場合。

· 始まる時間（） -この内部変数は、
サンプル空間時間ラスターデータセットの開始時間と
現在のサンプル間隔またはインスタンス。 時間は日数と分数で測定されます
絶対時間の場合は日数、相対時間の場合は相対単位で。

· 終了時間（） -この内部変数は、
サンプル時空間ラスターデータセットの開始時刻と現在の終了時刻
サンプル間隔。 時間は日数で測定され、絶対値の場合は日数の端数で測定されます
時間、および相対時間の場合は相対単位で。 end_time（）は
時間インスタンスの場合はnull（）で表されます。

現在のサンプル間隔またはインスタンスの絶対時間でサポートされている内部変数：

· start_doy（） -開始時刻からの年の日（doy）[1-366]

· start_dow（） -開始時刻[1〜7]からの曜日（ダウ）、週の開始
月曜日です== 1

· start_year（） -開始時刻の年[0-9999]

· start_month（） -開始時刻の月[1-12]

· start_week（） -開始時刻の年の週[1-54]

· start_day（） -開始時刻からの日[1〜31]

· start_hour（） -開始時刻の時間[0-23]

· start_minute（） -開始時刻の分[0-59]

· start_second（） -開始時間の0番目[59-XNUMX]

· end_doy（） -終了時刻からの年の日（doy）[1-366]

· end_dow（） -終了時刻[1-7]からの曜日（ダウ）、週の開始は
月曜日== 1

· end_year（） -終了時刻の年[0-9999]

· end_month（） -終了時刻の月[1-12]

· end_woy（） -終了時間の週（woy）[1-54]

· end_day（） -開始時刻からの日[1〜31]

· end_hour（） -終了時刻の時間[0-23]

· end_minute（） -終了時刻の分[0-59]

· end_second（） -終了時間の0番目[59-XNUMX]
当学校区の 終わり_* 関数は、時間の場合はnull（）内部変数で表されます
インスタンス。

注意事項

の内部作業について説明します t.rast.mapcalc 例を挙げて。 XNUMXつあると想像してください
入力としてのSTRDS。それぞれが月単位の粒度です。 最初のXNUMXつ A 6つのラスターマップがあります（a3..。
a8）XNUMX月からXNUMX月までの時間範囲。 XNUMX番目のSTRDS B 12個のラスターマップがあります（b1
... b12）XNUMX月からXNUMX月まで。 ラスターマップの値は、
インターバル開始時刻からの月。 データセット A データセットのサンプリングに使用されます B 〜へ
データセットを作成する C。 月が
開始時刻はXNUMX月またはXNUMX月です。それ以外の場合は、マップを乗算します。 コマンドは次のようになります
以下：

t.rast.mapcalc input = A、B output = C basename = c method = equal \
expression = "if（start_month（）== 5 || start_month（）== 6、（A + B）、（A * B））"

データセットCで作成されたラスターマップは、次のように一覧表示できます。 t.rast.list:

名前start_timemin max
c_1 2001-03-01 00:00:00 9.0 9.0
c_2 2001-04-01 00:00:00 16.0 16.0
c_3 2001-05-01 00:00:00 10.0 10.0
c_4 2001-06-01 00:00:00 12.0 12.0
c_5 2001-07-01 00:00:00 49.0 49.0
c_6 2001-08-01 00:00:00 64.0 64.0

内部的には、時空間表現は、
サンプルデータセットA、時間関数は数値、名前に置き換えられます
時空ラスターデータセットの一部は、対応するラスターマップに置き換えられます。 The
最終的な式はに渡されます r.mapcalc、結果として6回実行されます。

r.mapcalc expression = "c_1 = if（3 == 5 || 3 == 6、（a3 + b3）、（a3 * b3））"
r.mapcalc expression = "c_2 = if（4 == 5 || 4 == 6、（a4 + b4）、（a4 * b4））"
r.mapcalc expression = "c_3 = if（5 == 5 || 5 == 6、（a5 + b5）、（a5 * b5））"
r.mapcalc expression = "c_4 = if（6 == 5 || 6 == 6、（a6 + b6）、（a6 * b6））"
r.mapcalc expression = "c_5 = if（7 == 5 || 7 == 6、（a7 + b7）、（a7 * b7））"
r.mapcalc expression = "c_6 = if（8 == 5 || 8 == 6、（a8 + b8）、（a8 * b8））"

実施例

次のコマンドは、XNUMX月に新しいラスター時空データセットを作成しています
マップは、温度が0より大きい場合はヌルに設定され、それ以外の場合は設定されます。
元の値。 他の月は元の月としてコピーされます。
t.rast.mapcalc input = tempmean_monthly output = january_under_0 basename = january_under_0 \
expression = "if（start_month（）== 1 && tempmean_monthly> 0、null（）、tempmean_monthly）"
＃XNUMX月のみの最小値または最大値を印刷する
t.rast.list january_under_0 columns = name、start_time、min、max | grep 01-01
名前| start_time | min | max
january_under_0_01|2009-01-01 00:00:00|-3.380823|-7e-06
january_under_0_13|2010-01-01 00:00:00|-5.266929|-0.000154
january_under_0_25|2011-01-01 00:00:00|-4.968747|-6.1e-05
january_under_0_37|2012-01-01 00:00:00|-0.534994|-0.014581
＃これらは元のデータであり、最大値が異なることがわかります
t.rast.list tempmean_monthly columns = name、start_time、min、max | grep 01-01
2009_01_tempmean|2009-01-01 00:00:00|-3.380823|7.426054
2010_01_tempmean|2010-01-01 00:00:00|-5.266929|5.71131
2011_01_tempmean|2011-01-01 00:00:00|-4.968747|4.967295
2012_01_tempmean|2012-01-01 00:00:00|-0.534994|9.69511

onworks.netサービスを使用してオンラインでt.rast.mapcalcgrassを使用する