Una vez más, encogiéndome, me abroché completamente la chaqueta y aceleré el paso, apresurándome a trabajar. En el trabajo, café caliente, una unidad de sistema caliente y tareas interesantes. Pero no giraban en mi cabeza, mientras el viento frío revolvía en mi seno, sino una ligera sospecha vaga sobre los cambios de estaciones. ¿Era necesario verificar si el calentamiento y el enfriamiento llegaron a tiempo, y no antes? Tal vez las estaciones cambian?
No hace falta decir que, al llegar al trabajo, olvidé pensar en mi vaga sospecha. Pero no se disipó por completo, sino que ocasionalmente apareció en mi cabeza cuando, por ejemplo, sufrí de calor en mayo o vi un árbol amarillento en julio.
Entonces, sin seguir realmente el principio de que debes comenzar a hacer una idea en las primeras 48 horas, después de muchos meses hice una pequeña investigación para probar mi sospecha. Un recordatorio fue un artículo sobre Gismeteo .
Datos
El primer paso fue obtener los datos. Al final resultó que, encontrar temperaturas promedio diarias no es tan fácil como parece. Después de cinco minutos de buscar en Google, me encontré con una pregunta sobre la tostadora , donde llegué al sitio www.ncdc.noaa.gov . Pero los datos recibidos no me satisficieron (muchas omisiones en las fechas), continué buscando en Google, obtuve lo que estaba buscando en el sitio http://aisori-m.meteo.ru/waisori/ .
Pieza de datos:
29645 1955 1 1 "-31.0" 29645 1955 1 2 "-23.8" 29645 1955 1 3 "-18.7" 29645 1955 1 4 "-13.3" 29645 1955 1 5 "-16.8" 29645 1955 1 6 "-16.8" 29645 1955 1 7 "-18.6" 29645 1955 1 8 "-17.0" 29645 1955 1 9 "-6.7" 29645 1955 1 10 "-10.4"
Obteniendo temperaturas anuales con una compensación
En primer lugar, convertí el csv con los datos en una matriz de objetos.
Luego escribió una función que regresó un año desde la matriz con una compensación de un cierto número de días.
Con más detalle, tenemos una parte de la matriz:
... 1955 12 29 -34.5 1955 12 30 -37.7 1955 12 31 -36.3 1956 1 1 -29.2 1956 1 2 -32.9 1956 1 3 -38.1 1956 1 4 -36.6 ...
Si solicita todos los datos de 1956 en la función, simplemente devuelve, y si solicita datos con un desplazamiento de -1 día, la temperatura de cada día se tomará del día anterior. Es decir, para el 1 de enero de 1956, la temperatura se tomará el 31 de diciembre de 1955.
... 1956 1 1 -36.3 ← ( 1955 12 31) 1956 1 2 -29.2 ← ( 1956 1 1) 1956 1 3 -32.9 ← ( 1956 1 2) 1956 1 4 -38.1 ← ( 1956 1 3) ...
Obtener la temperatura mensual promedio
En el segmento resultante de la matriz para el año, calculé las temperaturas mensuales promedio, obteniendo así un vector de temperatura para el año.
[ -16.716129032258067, -16.317857142857143, -17.690322580645155, -11.57333333333333, 0.3967741935483872, 9.163333333333332, 17.845161290322586, 17.225806451612904, 15.490000000000004, 9.922580645161291, 1.7233333333333334, -8.26774193548387 ]
Distancia entre años compensados
Si tengo un vector de temperatura durante un año, entonces puedo comparar la distancia entre los vectores.
Es decir, si el clima apenas hubiera cambiado de un año a otro, entonces la distancia entre los vectores sería cero. Desde aquí decidí que compararía, por ejemplo, el vector 2001 y todos los vectores 2002 obtenidos por el desplazamiento.
... 2001 == 2002 -2 2001 == 2002 -1 2001 == 2002 2001 == 2002 +1 2001 == 2002 +2 ...
Después de tomar el desplazamiento de esa comparación, donde la distancia entre los vectores es mínima. Es decir, si resulta que la distancia entre 2001 y 2002 +2 mínima, entonces la temporada 2002 se desplazó 2 días más adelante en relación con 2001.
Mirando hacia el futuro, diré que tomé 1956 para el año inicial y miré el cambio durante 30 días de ida y vuelta.
Programé el algoritmo en js e hice un sitio pequeño para una fácil verificación.
Código fuente en github .
Verificación de supuestos
Kemerovo
Si! ¡Las estaciones comenzaron a cambiar más rápido desde 2013 en 10 días, y en promedio desde 1956 en 8 días!
Moscú
Casi no hay desplazamiento.
Arkhangelsk
Kaliningrado
San petersburgo
Tyumen
Vladivostok
Conclusión
La sensación no fue engañada, pero las sensaciones aún necesitan ser verificadas.
PD De hecho, estaba buscando a Novosibirsk. En el sitio de Novosibirsk no estaba , pero Kemerovo estaba completamente satisfecho.
PPS Como resultó un poco más tarde, el archivo faltaba en 2018, y necesitaba este enero, así que descargué el mes faltante del sitio pogodaiklimat.ru .