bin = [-3:1:25]*3600;
x = histc(time,bin)';

max(time)/3600
max(timestamps)/3600