我正在用Python修改纽约市地铁十字转门数据,以变成每个车站入口/出口的可视化。
到目前为止,我有一个基于开始(03-24-15
)和结束(03-27-15
)日期的入口/出口计数列表:
{
'endTime': '03-25-14T21:40:30',
'entriesDuringPeriod': 158,
'exitsDuringPeriod': 597,
'startTime': '03-25-14T17:03:23'
},
{
'endTime': '03-26-14T01:00:00',
'entriesDuringPeriod': 29,
'exitsDuringPeriod': 235,
'startTime': '03-25-14T21:00:00'
},
我遇到的问题是不同的时间段没有标准化,有时重叠。我希望能够浏览并创建另一个列表,将这些数字规范化为一小时增量。
我对 Python 时间处理不是很熟悉,我想知道是否有人可以提供一些有关如何开始获取字符串、将它们转换为日期对象以及根据时间划分值的信息。
最终可视化将使用 d3 进行可视化.js如果这很重要的话。
将字符串放入日期时间对象还不错:
from datetime import datetime
from time import time, mktime, strptime
def get_datetime( instr ):
return datetime.fromtimestamp(mktime(strptime(instr, '%m-%d-%yT%H:%M:%S')))
# eg: get_datetime( '03-25-14T21:20:30' ) => datetime.datetime(2014, 3, 25, 21, 20, 30)
分箱/规范化数据在很大程度上取决于您希望如何处理重叠的持续时间......例如。您是否要假设人们以线性方式到达和退出,因此如果时间戳是一个半小时,66%将进入整个小时,33%进入另一个部分小时?
编辑:根据OP的评论,这是完全功能的代码:
from datetime import timedelta
from collections import defaultdict
def add_datum( dd, v ):
end_dt = get_datetime(v['endTime'])
start_dt = get_datetime(v['startTime'])
total_duration = end_dt - start_dt
hour_start = datetime( year = start_dt.year,
month = start_dt.month,
day = start_dt.day,
hour = start_dt.hour )
hour_end = hour_start + timedelta( hours = 1 )
while hour_start < end_dt:
dt = min([hour_end, end_dt]) - max([ hour_start, start_dt ])
fraction = 1.0 * dt.total_seconds() / total_duration.total_seconds()
dd[ hour_start ]['hour'] = hour_start
dd[ hour_start ]['entries'] += v['entriesDuringPeriod'] * fraction
dd[ hour_start ]['exits'] += v['exitsDuringPeriod'] * fraction # exits
hour_start = hour_end
hour_end = hour_end + timedelta( hours = 1 )
return dd
dd = defaultdict(lambda: {'entries':0,'exits':0})
all_data = [{ 'endTime': '03-25-14T21:40:30',
'entriesDuringPeriod': 158,
'exitsDuringPeriod': 597,
'startTime': '03-25-14T17:03:23' },
{ 'endTime': '03-26-14T01:00:00',
'entriesDuringPeriod': 29,
'exitsDuringPeriod': 235,
'startTime': '03-25-14T21:00:00' }]
[ add_datum( dd, i ) for i in all_data ]
res = dd.values()
res.sort( key = lambda i: i['hour'] )
print res
# [{'entries': 32.28038732182594,
# 'exits': 121.97083057677271,
# 'hour': datetime.datetime(2014, 3, 25, 17, 0)},
# {'entries': 34.209418415829674,
# 'exits': 129.25963793829314,
# 'hour': datetime.datetime(2014, 3, 25, 18, 0)},
# {'entries': 34.209418415829674,
# 'exits': 129.25963793829314,
# 'hour': datetime.datetime(2014, 3, 25, 19, 0)},
# {'entries': 34.209418415829674,
# 'exits': 129.25963793829314,
# 'hour': datetime.datetime(2014, 3, 25, 20, 0)},
# {'entries': 30.34135743068503,
# 'exits': 146.00025560834786,
# 'hour': datetime.datetime(2014, 3, 25, 21, 0)},
# {'entries': 7.25,
# 'exits': 58.75,
# 'hour': datetime.datetime(2014, 3, 25, 22, 0)},
# {'entries': 7.25,
# 'exits': 58.75,
# 'hour': datetime.datetime(2014, 3, 25, 23, 0)},
# {'entries': 7.25,
# 'exits': 58.75,
# 'hour': datetime.datetime(2014, 3, 26, 0, 0)}]
解析给定时间戳的方法是:
import datetime
datetime.datetime.strptime("03-25-14T21:08:12", "%m-%d-%yT%H:%M:%S")
这将为您提供一个datetime
对象,该对象可与其他datetime
对象相媲美。datetime
库本身包含用于创建任意datetime
对象的所有函数,然后可用于排序/装箱目的。
鉴于数据中的时间段不统一且包含重叠,因此有几种可能的方法。如果你对进入和退出进行线性平均,你可以利用每个时间段并计算平均每小时发生多少次进入和退出,然后,给定一个小时,你可以遍历所有数据点,找出数据点与该小时的重叠程度(即 15 分钟或整小时),并将数据点每小时的平均进入/退出数由重叠百分比修改应用于累加器。
在伪代码方面:
assume we have an hour to data accumulation bag called htdab
for given hour in time series
for data point in all data points
if data point's start <= given hour <= data point's end
calculate percentage of overlap on top of given hour
add data point's entries per hour * percentage of overlap to our htdab's entry for given hour
add data point's exits per hour * percentage of overlap to our htdab's entry for given hour