Python繪制交通流折線圖詳情

更新時間：2022年06月15日 10:28:15 作者：居安_

這篇文章主要介紹了Python繪制交通流折線圖詳情，文章基于python的相關資料展開折線圖繪制的實現流程，感興趣的小伙伴可以參考一下

一、數據集下載

加州高速公路PEMS數據集

這里繪制PEMS04中的交通流量數據。該數據集中包含舊金山2018年1月1日至2月28日的29條道路上307個探測器每五分鐘收集的數據。

二、折線圖繪制

1、解壓npz文件

npz是一種numpy文件存儲的壓縮格式，可使用numpy進行讀取。

allow_pickle=True用于防止numpy版本過高帶來的錯誤。

data.files查看壓縮文件下的所有文件。

import numpy as np
data = np.load(數據集存放地址, allow_pickle=True)
print(data.files)

可以看到壓縮文件下只有data一個文件：

通過 data['data'] 即可對該數據集進行讀取。我們也可以查看一下該數據的維度。

print(data['data'].shape)

其維度如下：

16992 = 59天×24小時×12（每五分鐘統(tǒng)計一次流量數據），307為探測器數量，3為特征數。

2、折線圖繪制

數據中包含的三個特征為（交通流量，平均速度，平均占用率），取出繪圖需要的第一個特征（這里只繪制其中一個探測器的）。

flow = data['data'][:, 0, 0]

導入matplotlib包，通過pyplot來繪制最簡單的折線圖。

import matplotlib.pyplot as plt

由于希望觀測到工作日的流量變化，我們將橫坐標替換為周一至周日。舊金山2018年1月1日為周一。

繪制的全部代碼如下：

import numpy as np
import matplotlib.pyplot as plt
data = np.load(文件路徑, allow_pickle=True)
flow = data['data'][:, 0, 0]
print(len(flow))
week = ['Sun', 'Mon', 'Tue', 'Wed', 'Thu', 'Fri', 'Sat']
x = []
for i in range(59):
    x.append(week[(i + 1) % 7])
y = []
sum = 0
for i in range(len(flow)):
    if i == 0 or (i + 1) % 288 != 0:
        sum += flow[i]
    else:
        y.append(sum)
        sum = 0
 
fig = plt.figure(figsize=(15, 5))   # 圖片寬度設置的大一些
plt.title('traffic flow in San Francisco')
plt.xlabel('day')
plt.ylabel('flow')
plt.xticks(np.arange(59), x)
plt.plot(np.arange(59), y, linestyle='-')
fig.autofmt_xdate(rotation=45)  # x軸的刻度標簽逆時針旋轉45度
plt.show()

繪制結果如下：