t.rast.mapcalcgrass: en línea en la nube

PROGRAMA:

NOMBRE

t.rast.mapcalc - Realiza expresiones mapcalc espacio-temporales en mapas muestreados temporalmente
de datasets ráster de espacio-tiempo.

PALABRAS CLAVE

temporal, álgebra, raster, tiempo

SINOPSIS

t.rast.mapcalc
t.rast.mapcalc --ayuda
t.rast.mapcalc [-ns] entradas=nombre [,nombre , ...] expresión=cadena [Método=nombre [,nombre , ...]]
salida=nombre nombre base=nombre base [nproces=entero] [-exagerar] [-ayuda] [-verboso]
[-tranquilo] [-ui]

Banderas
-n
Registrar mapas nulos

-s
Verificar la superposición espacial

--Sobrescribir
Permitir que los archivos de salida sobrescriban los archivos existentes

--ayuda
Resumen de uso de impresión

--verboso
Salida del módulo detallado

--tranquilo
Salida de módulo silencioso

--ui
Forzar el inicio del cuadro de diálogo GUI

parámetros:
entradas=nombre nombre, ...] [requerido]
Nombre de los datasets ráster de espacio-tiempo de entrada

expresión=cadena [requerido]
Expresión mapcalc espacio-temporal

Método=nombre nombre, ...]
El método que se utilizará para muestrear el conjunto de datos de entrada.
Opciones: comienzo, durante, superposición, Contiene, igual, sigue, precede
Por defecto: igual

salida=nombre [requerido]
Nombre del dataset ráster de espacio-tiempo de salida

nombre base=nombre base [requerido]
Nombre base para mapas ráster de salida
Se adjuntará un sufijo numérico separado por un guión bajo para crear un
identificador

nproces=entero
Número de procesos r.mapcalc para ejecutar en paralelo
Por defecto: 1

DESCRIPCIÓN

t.rast.mapcalc realiza espacio-temporal mapacalc Expresiones en mapas de muestreos temporales
datasets ráster de espacio-tiempo (STRDS). Operadores espaciales y temporales y variables internas
están disponibles en la cadena de expresión. La descripción de los operadores espaciales,
funciones y variables internas está disponible en la página de manual de r.mapcalc. El temporal
Las funciones se describen en detalle a continuación.

Este módulo espera varios parámetros. Todos los datasets ráster de espacio-tiempo a los que se hace referencia
en la categoría Industrial. mapacalc expresión debe figurar en el Las opciones de entrada opción. La first ráster de espacio-tiempo
El conjunto de datos que se enumera como entrada se utilizará para muestrear temporalmente el resto del espacio-tiempo.
datasets ráster. El método de muestreo temporal se puede elegir utilizando el Método opción. La
El orden de los STRDS en la expresión mapcalc puede ser diferente al orden de los
STRDS está en la opción de entrada. El dataset ráster de espacio-tiempo resultante debe especificarse en
los salida opción junto con la bases nombre de los mapas ráster generados que están registrados
en el STRDS resultante. Los mapas vacíos resultantes del cálculo de mapas no son registrados por
defecto. Este comportamiento se puede cambiar con el -n bandera. La bandera -s se puede utilizar para asegurar
que solo se procesan los mapas relacionados con el espacio en los STRDS. Relacionado espacialmente significa que
Los mapas relacionados temporalmente se superponen en su extensión espacial.

El módulo t.rast.mapcalc admite procesamiento paralelo. La opción nproces especifica el
número de procesos que se pueden iniciar en paralelo.

Se debe proporcionar una expresión mapcalc para procesar los mapas muestreados temporales. Temporal
Las variables internas están disponibles además de las r.mapcalc operadores espaciales y
funciones:

Variables internas admitidas para tiempo relativo y absoluto:

· td () - Esta variable interna representa el tamaño del tiempo de muestra actual
intervalo en días y fracción de días para el tiempo absoluto, y en unidades relativas en
caso de tiempo relativo.

· hora de inicio() - Esta variable interna representa la diferencia de tiempo entre el
la hora de inicio del dataset ráster de espacio-tiempo muestral y la hora de inicio del
intervalo de muestra actual o instancia. El tiempo se mide en días y fracción de
días para el tiempo absoluto y en unidades relativas en el caso del tiempo relativo.

· hora de finalización() - Esta variable interna representa la diferencia de tiempo entre el
la hora de inicio del dataset ráster de espacio-tiempo muestral y la hora de finalización del
intervalo de muestra. El tiempo se mide en días y fracción de días para valores absolutos.
tiempo, y en unidades relativas en caso de tiempo relativo. El end_time () será
representado por null () en el caso de una instancia de tiempo.

Variables internas admitidas para el tiempo absoluto del intervalo o instancia de muestra actual:

· start_doy () - Día del año (doy) desde la hora de inicio [1 - 366]

· start_dow () - Día de la semana (dow) desde la hora de inicio [1 - 7], el inicio de la semana
es lunes == 1

· año de inicio() - El año de la hora de inicio [0 - 9999]

· mes_inicio () - El mes de la hora de inicio [1 - 12]

· semana_inicial () - Semana del año de la hora de inicio [1 - 54]

· día_inicio () - Día del mes a partir de la hora de inicio [1 - 31]

· hora_inicial () - La hora de la hora de inicio [0 - 23]

· minuto_inicio () - El minuto de la hora de inicio [0 - 59]

· start_second () - El segundo de la hora de inicio [0 - 59]

· end_doy () - Día del año (doy) desde la hora de finalización [1 - 366]

· end_dow () - Día de la semana (dow) desde la hora de finalización [1 - 7], el inicio de la semana es
lunes == 1

· año final() - El año de la hora de finalización [0 - 9999]

· end_month () - El mes de la hora de finalización [1 - 12]

· end_woy () - Semana del año (woy) del tiempo de finalización [1 - 54]

· día final() - Día del mes a partir de la hora de inicio [1 - 31]

· end_hour () - La hora de finalización [0 - 23]

· end_minute () - El minuto de la hora de finalización [0 - 59]

· end_second () - El segundo de la hora de finalización [0 - 59]
El fin_* las funciones están representadas por las variables internas null () en caso de tiempo
instancias.

NOTAS

Discutiremos el trabajo interno de t.rast.mapcalc con un ejemplo. Imagina que tenemos dos
STRDS como entrada, cada uno con granularidad mensual. El primero A tiene 6 mapas ráster (a3 ...
a8) con un rango temporal de marzo a agosto. El segundo STRDS B tiene 12 mapas ráster (b1
... b12) desde enero a diciembre. El valor de los mapas ráster es el número de
el mes desde la hora de inicio del intervalo. Conjunto de datos A se utilizará para muestrear el conjunto de datos B a
crear un conjunto de datos C. Queremos agregar todos los mapas con las mismas marcas de tiempo si el mes del
la hora de inicio es mayo o junio, de lo contrario multiplicamos los mapas. El comando se verá como
manera:

t.rast.mapcalc entrada=A,B salida=C nombre base=c método=igual
expresión = "si (mes_inicio () == 5 || mes_inicio () == 6, (A + B), (A * B))"

Los mapas ráster resultantes en el dataset C se pueden enumerar con t.rast.lista:

nombre start_time min 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

Internamente se analizará la expresión espacio-temporal para cada intervalo de tiempo del
conjunto de datos de muestra A, las funciones temporales serán reemplazadas por valores numéricos, los nombres
de los datasets ráster de espacio-tiempo serán reemplazados por los mapas ráster correspondientes. los
la expresión final se pasará a r.mapcalc, lo que da como resultado 6 ejecuciones:

r.mapcalc expresión = "c_1 = if (3 == 5 || 3 == 6, (a3 + b3), (a3 * b3))"
r.mapcalc expresión = "c_2 = if (4 == 5 || 4 == 6, (a4 + b4), (a4 * b4))"
r.mapcalc expresión = "c_3 = if (5 == 5 || 5 == 6, (a5 + b5), (a5 * b5))"
r.mapcalc expresión = "c_4 = if (6 == 5 || 6 == 6, (a6 + b6), (a6 * b6))"
r.mapcalc expresión = "c_5 = if (7 == 5 || 7 == 6, (a7 + b7), (a7 * b7))"
r.mapcalc expresión = "c_6 = if (8 == 5 || 8 == 6, (a8 + b8), (a8 * b8))"

EJEMPLO

El siguiente comando está creando un nuevo dataset ráster de espacio-tiempo donde en enero
Los mapas son si la temperatura es más de 0, se establece en nulo; de lo contrario, se establece el
Valor original. Los otros meses se copian como el original.
t.rast.mapcalc entrada=tempmean_monthly salida=enero_bajo_0 nombre base=enero_bajo_0
expresión = "if (start_month () == 1 && tempmean_monthly> 0, null (), tempmean_monthly)"
# imprimiendo los valores mínimos o máximos solo para enero
t.rast.list january_under_0 columnas = nombre, hora_inicio, mínimo, máximo | grep 01-01
nombre | hora_inicio | 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
# estos son los datos originales, puedes ver que el máximo es diferente
t.rast.list tempmean_monthly columnas = nombre, hora_inicio, mínimo, máximo | 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

Utilice t.rast.mapcalcgrass en línea utilizando los servicios de onworks.net