Calculate number of hours below a threshold temperature

Question

Paul Fremeau 2018-7-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/411236-calculate-number-of-hours-below-a-threshold-temperature

评论： Paul Fremeau 2018-7-20

I have a big timetable with temperature recorded every 10 minutes. I want to calculate the number of hours per day below a threshold temperature. Any ideas? The timetable is important because my original data is missing some timesteps, and I use retime() to fill in the blanks.

I tried using a for loop along with an if statement, but they don't seem to work with timetables.

I also tried converting the timetable back to an array, but ran into difficulties. My guess is that I'm missing the easy way to do this.

Thanks!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Kelly Kearney 2018-7-19

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/411236-calculate-number-of-hours-below-a-threshold-temperature#answer_329542

I think you should be able to accomplish this by retiming twice; once to fill in any gaps with whatever method you prefer, and the second time to aggregate over each day. If you assume each measurement is representative of a full 10 minutes, then you can simply sum the number of points per day that meet the threshold criteria.

First we build an example timetable with ~10-minute data with some gaps

t = datetime(2018,1,1) + days(0:minutes(10):days(2))';
temp = rand(size(t))*15 - 2;
T = timetable(t, temp);
T = T(rand(size(t))>0.1,:);

Now retime:

dt = minutes(10);
thresh = 0;
T2 = retime(T,'regular', 'linear', 'timestep', dt);
T3 = retime(T2,'daily',@(x) hours(sum(x<thresh)*dt));

4 个评论
显示 2更早的评论隐藏 2更早的评论

Kelly Kearney 2018-7-20

Trying to do run-length calculations in a one-liner gets a bit complicated, so I would write a separate subfunction for that purpose.

T2 = retime(T,'regular', 'linear', 'timestep', dt);
T3 = retime(T2, 'daily', @(x) consectime(x,thresh,dt));
function y = consectime(v, th, dt)
      [b,n] = RunLength(v > th);
      y = hours(max(n(b==1)) * dt);
end

The T3.temp field now holds the length (in hours) of the longest consecutive period over the threshold, per day. I use the RunLength function from the FEX to do a lot of the heavy lifting here.

Paul Fremeau 2018-7-20