Skip to contents

Las Funciones locales de LCZ4r son herramientas especializadas diseñadas para analizar grandes conjuntos de datos, como lecturas de temperatura del aire por hora. Este tutorial se centra en el uso de datos de Berlín, recopilados en 23 estaciones meteorológicas como parte del Observatorio del Clima Urbano (UCO).

Descripción general

Las funciones locales en LCZ4r, ofrecen poderosas herramientas para manejar y analizar datos climáticos urbanos. Permiten análisis de series temporales, mapeo de anomalías térmicas, interpolación espacial y más. A continuación se muestra un resumen de estas funciones:

Función Descripción Datos requeridos Se requiere acceso a Internet
lcz_ts() Analizar series temporales LCZ No es necesario
lcz_anomaly() Calcular anomalías térmicas LCZ No es necesario
lcz_anomaly_map() Mapa Anomalías Térmicas LCZ No es necesario
lcz_interp_map() Realizar interpolación LCZ No es necesario
lcz_plot_interp() Visualice la interpolación LCZ No es necesario
lcz_interp_eval() Evaluar la interpolación LCZ No es necesario
lcz_uhi_intensity() Evaluar LCZ para determinar la intensidad de la isla de calor urbana No es necesario

💡 Consejos:

  1. Utilice el help(lcz_*) función para acceder a documentación detallada para cada función LCZ4r. Por ejemplo, para aprender sobre el lcz_ts función, tipo help("lcz_ts") en la consola.

  2. Cada función LCZ4r admite la imputación para manejar valores faltantes dentro de los marcos de datos. Para obtener más información sobre los métodos de imputación, consulte el argumento “imputar” en la documentación de cada función.

Requisitos de entrada de datos

Para garantizar el buen funcionamiento de las funciones locales, los datos de entrada deben estructurarse como un marco de datos con las siguientes columnas:

  1. fecha: esta columna debe contener información de fecha y hora. Esta columna debe estar etiquetada como “fecha”, “hora”, “marca de tiempo” o “fecha y hora”. Asegúrese de que el formato de fecha y hora se alinee con las convenciones de R (por ejemplo, AAAA-MM-DD HH:MM:SS). Los formatos aceptables incluyen “1999-02-01”. Para obtener más detalles, consulte formato de fechas y horas en la paquete al aire libreR.

  2. Identificador de estación: una columna para identificar estaciones meteorológicas.

  3. Temperatura del aire u otra variable: Al menos una columna que represente la temperatura del aire o cualquier otra variable objetivo.

  4. Latitud y Longitud: dos columnas para las coordenadas geográficas de cada estación. Asegúrese de que la columna tenga el nombre “lat” o “latitude” y “lon”, “long” o “longitude”.

Para simplificar esta configuración, LCZ4r proporciona un marco de datos de muestra, que puede cargar con el siguiente comando:

library(LCZ4r)

# Cargar datos del paquete LCZ4r
data("lcz_data")

# Verificar la estructura de datos
str(lcz_data)
Structure of the LCZ sample dataset

Estructura del conjunto de datos de muestra LCZ que muestra columnas de fecha, estación, temperatura y coordenadas

Importa tus datos

if(!require(data.table)) install.packages("data.table")

# Replace the path with the actual path to your CSV file
my_data <- data.table::fread("PC/path/file_name.csv")

head(my_data)

Personalización de funciones locales en LCZ4r

Puede personalizar aún más las funciones locales para satisfacer necesidades de análisis específicas:

1. Flexibilidad en la elección del horario

Las Funciones Locales tienen un argumento ... que proporciona opciones para filtrar datos por años, meses, días y horas específicos para limitar el período de análisis. Ejemplos de cómo utilizar estos argumentos incluyen:

  • Año(s): Valor(es) numérico(s) que especifican el(los) año(s) a seleccionar. Por ejemplo, year = 1998:2004 selecciona todos los años entre 1998 y 2004 (inclusive), mientras que year = c(1998, 2004) selecciona sólo los años 1998 y 2004.

  • Mes(es): Valor(es) numérico(s) o de caracteres que especifican los meses a seleccionar. Ejemplos numéricos: month = 1:6 (enero a junio), o ejemplos de personajes: month = c("January", "December").

  • Día(s): Valor(es) numérico(s) que especifican los días a seleccionar. Por ejemplo, day = 1:30 selecciona los días del 1 al 30, o day = 15 selecciona sólo el día 15 del mes.

  • Hora(s): Valor(es) numérico(s) que especifican las horas a seleccionar. Por ejemplo, hour = 0:23 selecciona todas las horas del día, mientras hour = 9 selecciona sólo la novena hora.

  • Fecha de inicio: una cadena que especifica la fecha de inicio en cualquiera de los dos start = "DD/MM/YYYY" (p. ej., “2/1/1999”) o formato “AAAA-mm-dd” (p. ej., “1999-02-01”).

  • Fecha de finalización: una cadena que especifica la fecha de finalización en cualquiera de los dos casos. end = "DD/MM/YYYY" (p. ej., “2/1/1999”) o formato “AAAA-mm-dd” (p. ej., “1999-02-01”).

Ejemplos

# Seleccione un rango de años (e.g., 1998:2004) o años específicos (e.g., c(1998, 2004))
lcz_ts(
  lcz_map,
  year = 1998:2004
)

# Filtrar por mes, ya sea por valores numéricos. (e.g., 1:6 for January to June) o por nombres
lcz_anomaly(
  lcz_map,
  year = 2012, month = 9
)

# Indique una fecha específica, como Septiembre 1, 2019
lcz_interp_map(
  lcz_map,
  year = 2012, month = 9, day = 1 
)

# Durante una hora específica
lcz_anomaly_map(
  lcz_map,
  year = 2012, month = 9, day = 1, hour = 5
)

# Delimite un período específico utilizando las fechas de inicio y fin.
lcz_uhi_intensity(
  lcz_map,
  var = "airT", 
  station_id = "station",
  start = "1/9/2019", end = "30/10/2019"
)

Para obtener más información, consulte la funciones de utilidad del paquete openair.

2. División de series temporales LCZ por ventana temporal o sitio

Puede segmentar los datos de la serie temporal LCZ utilizando el by argumento, permitiendo el análisis a través de diferentes ventanas temporales o sitios. Las opciones de división disponibles incluyen:

  1. Segmentos temporales: “año”, “temporada”, “temporadaaño”, “mes”, “mesaño”, “día laborable”, “fin de semana”.

  2. División de luz diurna: “luz diurna”, que divide los datos en períodos diurnos y nocturnos. Tenga en cuenta que la opción Luz diurna puede dar como resultado que las horas diurnas y nocturnas se representen en UTC. Ver NOAA y argumento type en paquete al aire libre.

  3. Combinaciones: También puedes usar combinaciones como c("daylight", "month") o c("daylight", "season").

# Dividido por mes
lcz_ts(lcz_map, by = "month")

# Dividido por temporada
lcz_ts(lcz_map, by = "season")

# Dividido por las horas de luz diurna
lcz_anomaly(lcz_map, by = "daylight")

# Combinar múltiples divisiones
lcz_anomaly(lcz_map, by = c("daylight", "season"))

Para obtener más detalles, consulte la escriba el paquete al aire libre.

Flujos de trabajo de ejemplo

A continuación se muestran algunos ejemplos prácticos de cómo utilizar las funciones locales:

Análisis de series de tiempo

# Serie temporal básica para un año específico
lcz_ts(
  lcz_map,
  var = "airT",
  station_id = "station",
  year = 2020
)
LCZ time series analysis showing temperature variations across different LCZ classes

Análisis de series temporales de LCZ que muestran variaciones de temperatura en diferentes clases de LCZ

Mapeo de anomalías térmicas

# Mapeo de anomalías térmicas para los meses de verano
lcz_anomaly_map(
  lcz_map,
  var = "airT",
  station_id = "station",
  month = 6:8,
  year = 2020
)
Thermal anomaly map showing temperature differences across urban zones

Mapa de anomalías térmicas que muestra las diferencias de temperatura entre zonas urbanas

Interpolación espacial

# Interpolar la temperatura para una fecha específica
lcz_interp_map(
  lcz_map,
  var = "airT",
  station_id = "station",
  year = 2020, month = 7, day = 15, hour = 12
)
Spatial interpolation of air temperature across the study area

Interpolación espacial de la temperatura del aire en el área de estudio.


¿Tiene comentarios o sugerencias?

¿Tiene una idea para mejorar o detectó un error? ¡Nos encantaría saber de usted! Haga clic en el botón a continuación para crear una nueva edición (GitHub) y compartir sus comentarios o sugerencias directamente con nosotros.

Abrir incidencia en GitHub