Visualize padrões sazonais
Compare precipitação sazonal na região centro oeste e uma sudoeste.
WeatherData fornece precipitação acumulada em cm.
In[1]:=
start = {1999, 1, 1};
end = {2015, 12, 31};In[2]:=
getData[place_] :=
WeatherData[place, "TotalPrecipitation", {start, end, "Month"}]Crie uma TimeSeries de precipitação para Champaign, IL, e Phoenix, AZ. Especifique MissingDataMethod para usar interpolação para completar os dados faltantes.
In[3]:=
makeTS[place_] :=
TimeSeries[getData[place], MissingDataMethod -> "Interpolation"]In[4]:=
tsIL = makeTS[
Entity["City", {"Champaign", "Illinois", "UnitedStates"}]]Out[4]=
In[5]:=
tsAZ = makeTS[Entity["City", {"Phoenix", "Arizona", "UnitedStates"}]]Out[5]=
Crie WeightedData com quantidades de precipitação como ponderações e faça um gráfico DateHistogram calculando a média dos dados sobre o número total de anos.
In[6]:=
dataIL = WeightedData[tsIL["Dates"], tsIL["Values"]];
dataAZ = WeightedData[tsAZ["Dates"], tsAZ["Values"]];In[7]:=
totalyears = (end[[1]] - start[[1]] + 1);In[8]:=
avg := Function[{x, y}, y/totalyears]In[9]:=
makeHistogram[data_, lab_] :=
DateHistogram[data, "Month", avg, DateReduction -> "Year",
PlotRange -> {All, {-1, 9}}, AxesLabel -> "cm", PlotLabel -> lab]In[10]:=
{makeHistogram[dataIL, "IL"], makeHistogram[dataAZ, "AZ"]}Out[10]=
