Применение фильтров к временным рядам с датами
Получим температуру для заданных дат с помощью WeatherData. Создадим новый временной ряд и воспользуемся MissingDataMethod для восполнения недостающих данных.
In[1]:=
![Click for copyable input](assets.ru/filter-time-series-with-dates/In_93.png)
data = TimeSeries[
WeatherData["Chicago",
"Temperature", {{2015, 1, 1}, {2015, 12, 31}}],
MissingDataMethod -> "Interpolation"]
Out[1]=
![](assets.ru/filter-time-series-with-dates/O_49.png)
In[2]:=
![Click for copyable input](assets.ru/filter-time-series-with-dates/In_94.png)
DateListPlot[data, FrameLabel -> Automatic]
Out[2]=
![](assets.ru/filter-time-series-with-dates/O_50.png)
Bоспользуемся фильтром MinFilter с месячным диапазоном, чтобы определить даты первых и последних заморозков в конкретном году.
In[3]:=
![Click for copyable input](assets.ru/filter-time-series-with-dates/In_95.png)
minF = MinFilter[data, Quantity[1, "Month"]];
код на языке Wolfram Language целиком
Out[5]=
![](assets.ru/filter-time-series-with-dates/O_51.png)
Первый день после весенних заморозков.
In[6]:=
![Click for copyable input](assets.ru/filter-time-series-with-dates/In_98.png)
spring = TimeSeriesWindow[minF, {"15 April 2015", "15 May 2015"}];
In[7]:=
![Click for copyable input](assets.ru/filter-time-series-with-dates/In_99.png)
DateObject[
First@FirstCase[
spring["Path"], _?(#[[2]] > Quantity[0, "DegreesCelsius"] &)]]
Out[7]=
![](assets.ru/filter-time-series-with-dates/O_52.png)
Первый день осенних заморозков.
In[8]:=
![Click for copyable input](assets.ru/filter-time-series-with-dates/In_100.png)
fall = TimeSeriesWindow[minF, {"1 Oct 2015", "31 Oct 2015"}];
In[9]:=
![Click for copyable input](assets.ru/filter-time-series-with-dates/In_101.png)
DateObject[
First@FirstCase[
fall["Path"], _?(#[[2]] <= Quantity[0, "DegreesCelsius"] &)]]
Out[9]=
![](assets.ru/filter-time-series-with-dates/O_53.png)