MIT híbrido v0.9.3 — Piloto ciego de validación meteorológica

Alcance y objetivo

Este módulo no calcula un nuevo MIT, no clasifica automáticamente tormentas y no descarga radar. Genera un dossier de observaciones pluviométricas para una evaluación meteorológica independiente de las predicciones de conectividad (v0.9–v0.9.2).

Selecciona 18 parejas de episodios focal–control = 36 casos, priorizando años recientes, ambos semestres, los tres motivos de revisión y una calidad de emparejamiento estricta: diferencia de duración ≤30 minutos y de año hidrológico ≤2 años. Los controles son episodios diferentes, en fechas distintas, no representaciones alternativas de la misma tormenta. La muestra es dirigida, no aleatoria: sus proporciones no deben extrapolarse a la población de 11.366 huecos.

Los 18 pares se determinan en un orden reproducible y se conserva una tabla confidencial de correspondencias. La documentación entregada al revisor no revela clase de conectividad, condición de focal/control, discrepancia MIT ni la identificación de la pareja. Los gráficos únicamente muestran las observaciones y el hueco de Málaga.

Instalación

Copiar 09c_piloto_radar_22.R a

C:/Users/algar/OneDrive/000-FUNDAMENTOS/DATOS/Hidrosur/Datos_Hidrosur/Calculo_MIT-Hibrido/

Entradas existentes, sin copiar de nuevo:

  1. Resultados/22/barrera_180dias/Propagacion_v092/03_plan_radar_con_controles.csv (v0.9.2; 238 filas).
  2. Resultados/22/barrera_180dias/Espacial_v071/Series_5min/A_regional_aniohid_XXXX.rds, para los años seleccionados.
  3. Resultados/22/barrera_180dias/Espacial_v071/Cache_estaciones/vecina_{19,20,21,101,25}_SENSOR.rds. Estas son las cachés depuradas de las vecinas y permiten extraer lluvia antes y después del hueco sin releer sus CSV originales.
  4. Hidrosur/22.csv y el lector ya usado Calculo_MIT-Hibrido/05_validacion_espacial_22.R. La estación 22 no está en los RDS de vecinas; el código la lee con el control de calidad de v0.6 una vez, guarda únicamente el recorte piloto de las 36 ventanas y lo reutiliza al repetir la ejecución.
  5. Calculo_MIT-Hibrido/Docs/inventario_estaciones.csv, con los metadatos de las seis estaciones.

Requisitos: R >= 4.1 y paquete data.table. Los gráficos usan exclusivamente R base. No se requieren conexiones a internet.

Ejecución

ruta <- "C:/Users/algar/OneDrive/000-FUNDAMENTOS/DATOS/Hidrosur/Datos_Hidrosur/Calculo_MIT-Hibrido"
source(file.path(ruta,"09c_piloto_radar_22.R"))

# Prueba sintética, sin leer grandes archivos
v093_autoprueba()

# Debe anunciar 18 parejas, 36 casos y disponibilidad de entradas
v093_preflight()

# Ejecución (la primera vez procesa solamente 22.csv; no toca las 5 vecinas originales)
res093 <- v093_ejecutar()

Para construir solo los CSV sin dibujar: v093_ejecutar(generar_figuras=FALSE).

El preflight comprueba existencia de entradas y cardinalidad de parejas antes del trabajo costoso. La ejecución compara los RDS originales con las cinco cachés antes de leer 22.csv; cualquier inconsistencia cancela el proceso. Los RDS se consultan año a año. Se verifica que Málaga 22 sea seca dentro de los huecos Q y que la serie y el registro quinqueminutal coincidan exactamente.

Salidas (carpeta nueva, no modifica v0.9.2)

Resultados/22/barrera_180dias/Piloto_radar_v093/

Carpeta PARA_REVISION_CIEGA — única que debe compartirse con revisores externos:

  • 00_metadatos_estaciones.csv: estaciones y coordenadas y altitud originales.
  • 02_ventanas_radar_PROVISIONALES.csv: intervalos para localizar información radar/reanálisis, identificados solo por REV_XXX; no contiene las decisiones del grafo.
  • 03_cronograma_5min_seis_estaciones.csv: matriz larga de seis pluviométricos y huecos (sin columnas de clasificación).
  • Figuras/REV_XXX.pdf: 36 cronogramas individuales con seis paneles de lluvia (mm/5 min), un panel de fracción de vecinas con precipitación y sombreado del periodo seco. Una cruz indica dato ausente, no sequedad.
  • Etiquetas/01_plantilla_revision_CIEGA_NO_SOBRESCRIBIR.csv: hoja de revisión de la identidad meteorológica; las respuestas están vacías. Si la hoja ya existe, no se sobrescribe.

Carpeta RESERVADO_ANALISTA — NO entregar a quienes etiqueten:

  • 00_claves_y_seleccion.csv: mapeo de ID ciego a estación, hueco, motivo de selección y configuración, con información sensible para el cegamiento.
  • 01_auditoria_seleccion.csv: control de motivos, estaciones y distancia de emparejamiento.

En raíz e Intermedios: auditoría de RDS frente a cachés; sequedad de Málaga; QC de estación 22; fechas; recorte de lluvia de Málaga en RDS reutilizable; auditoría final.

Advertencia crítica: no confundir la base temporal con UTC

La serie de Hidrosur contiene fechas que en las transformaciones históricas aparecen con un sufijo Z. Ese carácter no es prueba de que el sistema original registrase en UTC. La versión usa segundos en una escala técnica sin cambio horario y no convierte a hora civil por defecto.

  • Los campos ventana_inicio_serie y ventana_fin_serie son válidos como etiquetas del reloj original de análisis, sin asignarles zona horaria.
  • ventana_inicio_UTC y ventana_fin_UTC están vacíos hasta que se confirme documentalmente el huso del proveedor.
  • Solo si se demuestra que todas las observaciones relevantes comparten un desfase UTC constante, establecer v093$ajuste_utc_min <- ... antes de ejecutar. Si fueron guardadas en hora local y experimentaron horarios de verano, no utilizar un desplazamiento constante sin resolver las transiciones: habría que adaptar la conversión fecha a fecha.
  • No comparar directamente estos relojes provisionales con radar de AEMET u otras fuentes hasta resolver esa cuestión.

Protocolo sugerido para revisión externa

  1. Auditar el reloj de Hidrosur 5-min y la resolución temporal/convención del producto radar.
  2. Conseguir información meteorológica independiente para cada ventana: radar con barridos consecutivos, mapas sinópticos y/o reanálisis, dejando constancia de cobertura y resolución. El script no asegura disponibilidad histórica ni inventa enlaces de descarga.
  3. Compartir exclusivamente PARA_REVISION_CIEGA. No proporcionar RESERVADO_ANALISTA ni el análisis v0.9.2 durante la primera lectura.
  4. Anotar misma_tormenta = 1 (mismo sistema), 0 (sistemas distintos), vacío (indeterminado). Documentar confianza, criterios y fuentes. La coincidencia de precipitación en estaciones no basta para etiquetar.
  5. Si es viable, usar dos revisores y resolver discrepancias con adjudicación documentada. No rellenar etiquetas por algoritmos ni por umbrales MIT.
  6. Solo después de cerrar las etiquetas, recuperar el mapeo privado y medir tasas de error de reglas MIT y modelos contextuales mediante evaluación fuera de muestra. El piloto dirigido no sustituye una muestra probabilística para estimar error poblacional.

Protección frente a errores / reinicio

  • Se comprueba que cada pareja tenga exactamente un focal y un control, que el emparejamiento sea estrictamente válido y que las 36 claves sean únicas.
  • Se verifica cada malla temporal contra el RDS anual original y la caché de cada vecina: no se interpola ni se convierte NA en cero.
  • Ante cualquier discrepancia, el programa se detiene antes de emitir una revisión como válida.
  • Si una ejecución posterior encuentra la misma carpeta de resultados y etiquetas existentes, se exige la misma correspondencia de IDs o se detiene para no asociar revisiones a casos distintos.
  • Las etiquetas de revisión no se sobrescriben; las demás auditorías se pueden regenerar.

Comprobaciones externas

Las estructuras de los CSV aportados se han verificado en un entorno de análisis Python; las 238 filas contienen 130 focales, 108 controles y 80 parejas con diferencias <=30 min y <=2 años. El criterio implementado selecciona 18 parejas (36 casos) repartidas por ambos semestres y los tres motivos de investigación, concentrando 16 de 18 focales en 2017 o después.

El script ha pasado inspección estática de delimitadores y se ha contrastado contra la estructura conocida de v0.6 y v0.7.1. No se ha podido ejecutar el programa en R con los cachés/RDS locales, por lo que primero deben ejecutarse la autoprueba y el preflight indicados.