Mais uma vez, encolhendo-me, prendi completamente minha jaqueta e acelerei meu ritmo, correndo para o trabalho. No trabalho, café quente, uma unidade de sistema quente e tarefas interessantes. Mas eles não giraram na minha cabeça, enquanto o vento frio revolveu no meu peito, mas uma leve suspeita vaga sobre as mudanças da estação. Era necessário verificar se o aquecimento e o resfriamento chegaram a tempo, e não antes? Talvez as estações mudem?
Escusado será dizer que, ao chegar ao trabalho, esqueci de pensar na minha vaga suspeita. Mas não se dissipou completamente, mas ocasionalmente surgia na minha cabeça quando, por exemplo, sofria de calor em maio ou via uma árvore amarelada em julho.
Portanto, sem realmente seguir o princípio de que você precisa começar a fazer uma idéia nas primeiras 48 horas, depois de muitos meses, fiz uma pequena pesquisa para testar minha suspeita. Um lembrete foi um artigo sobre Gismeteo .
Dados
O primeiro passo foi obter os dados. Como se viu, encontrar temperaturas médias diárias não é tão fácil quanto parece. Após cinco minutos de pesquisa, me deparei com uma pergunta na Torradeira , onde cheguei ao site www.ncdc.noaa.gov . Mas os dados recebidos não me satisfizeram (muitas omissões nas datas), continuaram pesquisando no Google, consegui o que procurava no site http://aisori-m.meteo.ru/waisori/ .
Dados:
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"
Obtenção de temperaturas anuais com um deslocamento
Antes de tudo, converti o csv com os dados em uma matriz de objetos.
Em seguida, ele escreveu uma função que retornava um ano da matriz com um deslocamento de um determinado número de dias.
Mais detalhadamente, temos uma parte da 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 ...
Se você solicitar todos os dados de 1956 na função, ele simplesmente retornará e, se você solicitar dados com um deslocamento de -1 dia, a temperatura de cada dia será obtida no dia anterior. Ou seja, para 1º de janeiro de 1956, a temperatura será medida em 31 de dezembro 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) ...
Obtendo a temperatura média mensal
Na fatia resultante da matriz para o ano, calculei as temperaturas médias mensais, obtendo assim um vetor de temperatura para o ano.
[ -16.716129032258067, -16.317857142857143, -17.690322580645155, -11.57333333333333, 0.3967741935483872, 9.163333333333332, 17.845161290322586, 17.225806451612904, 15.490000000000004, 9.922580645161291, 1.7233333333333334, -8.26774193548387 ]
Distância entre os anos de compensação
Se eu tiver um vetor de temperatura por um ano, posso comparar a distância entre os vetores.
Ou seja, se o tempo mal tivesse mudado de ano para ano, a distância entre os vetores seria zero. A partir daqui, decidi que compararia, por exemplo, o vetor de 2001 e todos os vetores de 2002 obtidos pelo deslocamento.
... 2001 == 2002 -2 2001 == 2002 -1 2001 == 2002 2001 == 2002 +1 2001 == 2002 +2 ...
Depois de tirar o deslocamento dessa comparação, onde a distância entre os vetores é mínima. Ou seja, se a distância entre 2001 e 2002 +2 mínima, a temporada de 2002 foi deslocada 2 dias à frente em relação a 2001.
Olhando para o futuro, direi que tirei 1956 no ano inicial e observei a mudança por 30 dias.
Programei o algoritmo em js e criei um site pequeno para facilitar a verificação.
Código fonte no github .
Verificando Premissas
Kemerovo
Sim! As estações começaram a mudar mais rapidamente a partir de 2013 por 10 dias e, em média, desde 1956 por 8 dias!
Moscovo
Quase não há deslocamento.
Arkhangelsk
Kaliningrado
São Petersburgo
Tyumen
Vladivostok
Conclusão
A sensação não foi enganada, mas as sensações ainda precisam ser verificadas.
PS Na verdade, eu estava procurando Novosibirsk. No site Novosibirsk não era , mas Kemerovo eu estava completamente satisfeito.
PPS Como se viu um pouco mais tarde, o arquivo estava ausente em 2018, e eu precisava de janeiro deste ano, então baixei o mês que faltava no site pogodaiklimat.ru .