蔡曉禹，呂 亮，杜 蕊

(1. 重慶交通大學(xué) 交通運(yùn)輸學(xué)院，重慶 400074；2. 山地城市交通系統(tǒng)與安全重慶市重點(diǎn)實(shí)驗(yàn)室，重慶 400074)

0 引言

隨著區(qū)域一體化進(jìn)程的加快，城市出行需求增長(zhǎng)迅速，私家車作為城市交通的主要組成部分，具有靈活、機(jī)動(dòng)的優(yōu)勢(shì)。在綜合運(yùn)輸競(jìng)爭(zhēng)日益激烈的背景下，研究私家車的出行特征，有助于改善城市交通組織，為城市智能的交通發(fā)展奠定基礎(chǔ)。

近些年，RFID(射頻識(shí)別技術(shù))在交通領(lǐng)域中被廣泛挖掘應(yīng)用，主要體現(xiàn)在公路收費(fèi)、停車場(chǎng)管理、車輛身份識(shí)別等方面[2-3]。根據(jù)相關(guān)外文文獻(xiàn)，美國在RFID的應(yīng)用方面處于相對(duì)領(lǐng)先的地位，英國、德國、瑞士、瑞典、日本等國均具有較為成熟且先進(jìn)的RFID產(chǎn)品，如美國的E-Zpass 收費(fèi)、新家坡的不停車收費(fèi)以及英國倫敦的交通擁堵費(fèi)。此外德國也研究了RFID技術(shù)在交通信息采集方面的重要應(yīng)用，德國柏林交通管理中心在 2000 年建立基于 RFID技術(shù)的交通信息采集系統(tǒng)。21世紀(jì)初，RFID逐漸進(jìn)入中國市場(chǎng)，并得到了政府的大力支持。2019年6月，伴隨著《中國RFID技術(shù)政策白皮書》的發(fā)布，我國的RFID技術(shù)進(jìn)入了國家產(chǎn)業(yè)發(fā)展戰(zhàn)略層面?，F(xiàn)如今，我國的 RIFD技術(shù)發(fā)展相對(duì)成熟的城市有南京、杭州、無錫、深圳和重慶等地。上海、杭州、重慶等多個(gè)城市依靠設(shè)置在道路上的RFID系統(tǒng)對(duì)公交進(jìn)行實(shí)時(shí)監(jiān)控，實(shí)現(xiàn)對(duì)車輛到站、離站時(shí)間的預(yù)測(cè)，以及配合紅綠燈實(shí)現(xiàn)公交優(yōu)先通行等。隨著我國電子車牌國家標(biāo)準(zhǔn)的即將出臺(tái)以及物聯(lián)網(wǎng)技術(shù)的持續(xù)發(fā)展，RFID 技術(shù)將會(huì)在自動(dòng)駕駛、智能交通管理、智慧城市發(fā)展中具有更加廣闊的應(yīng)用前景。

不同時(shí)間段里城市人群移動(dòng)的動(dòng)態(tài)變化規(guī)律，體現(xiàn)了城市交通流體變化規(guī)律，隱含城市地域功能性信息等，其實(shí)質(zhì)是對(duì)城市居民出行特征的研究分析。因此，車輛出行特征研究在城市建設(shè)規(guī)劃等方面具有重要的意義，同時(shí)也是制訂交通政策的科學(xué)依據(jù)。據(jù)相關(guān)國外研究，Sini Guo等[4]人主要研究軌跡數(shù)據(jù)的分割技術(shù)，從而分析車輛行為，提取商務(wù)信息，提出了一種集二元理論和牛頓法相結(jié)合的高效數(shù)據(jù)值算法來獲取最優(yōu)值。Zun Wang等[5]學(xué)者利用GPS數(shù)據(jù)中瞬時(shí)速度分布的變化系數(shù)來測(cè)量行程時(shí)間可靠性，提出旅行時(shí)間可靠性與道路密度的模型。Csáji等[6]基于大約10萬組手機(jī)信令數(shù)據(jù)發(fā)掘了葡萄牙用戶的活動(dòng)模式。MingqiLv等[7]人提出一種序列模式挖掘算法，通過考慮估計(jì)的空間接近度，從單元軌跡數(shù)據(jù)集中挖掘運(yùn)動(dòng)模式。A. Ladino等[8]就提出了一種行程時(shí)間的動(dòng)態(tài)聚類模型，利用聚類的質(zhì)心級(jí)別與趨勢(shì)變化規(guī)律來預(yù)測(cè)出行的旅行時(shí)間。Justin van Dijk[9]以軌跡數(shù)據(jù)提取利用不同的4種機(jī)器學(xué)習(xí)算法來識(shí)別軌跡數(shù)據(jù)中的活動(dòng)點(diǎn)和旅行點(diǎn)，引入局部密度概念并評(píng)估了4種算法的性能。在國內(nèi)的研究中，董宏輝等[10]人依據(jù)軌跡數(shù)據(jù)研究出行者的出行方式。充分挖掘GPS數(shù)據(jù)，提出兩種基于相似性度量和窗口的轉(zhuǎn)換點(diǎn)識(shí)別方法，然后將該方法運(yùn)用到BP神經(jīng)網(wǎng)絡(luò)、決策樹、KNN和支持向量機(jī)(SVM)4種識(shí)別交通出行的算法中，最后比較得出采用SVM識(shí)別交通出行方式能夠取得最優(yōu)結(jié)果。丁玲等[11]研究了在公共交通優(yōu)先政策下對(duì)出行個(gè)體行為的影響。研究根據(jù)分析將出行個(gè)體劃分為了3類，并采用聚類的方法將出行個(gè)體進(jìn)行劃分。再基于非集計(jì)模型對(duì)出行行為進(jìn)行分析，研究綜合考慮了交通狀態(tài)影響以及公共交通優(yōu)先策略影響下，出行個(gè)體的行為變化。

1 車輛出行軌跡分析

城市車輛出行是出行者利用私家車、公交車以及出租車等方式在不同道路間的出行過程。本研究基于軌跡數(shù)據(jù)分析重點(diǎn)分析私家車及出租車輛的活動(dòng)特點(diǎn)。

對(duì)于私家車軌跡研究，隨機(jī)選取了近千輛私家車個(gè)體進(jìn)行軌跡信息可視化定性分析。通過大量對(duì)私家車軌跡的觀察，發(fā)現(xiàn)私家車出行集中呈現(xiàn)某些統(tǒng)計(jì)特征，即私家車軌跡檢測(cè)數(shù)據(jù)大部分較短，但是大部分完整的RFID點(diǎn)位數(shù)據(jù)存在規(guī)律性。根據(jù)大量私家車出行軌跡分析，定性主要可以分析出存在兩類典型出行特征情況。

根據(jù)出行軌跡可視化分析，可知存在一部分私家車群體在不同工作日內(nèi)的出行會(huì)呈現(xiàn)高度重疊現(xiàn)象，例如周一至周五每天的出行時(shí)間以及軌跡基本一致，僅存在細(xì)微差別，特別是早高峰和晚高峰時(shí)期，例如圖1中給出了車輛編號(hào)為1005406和102563的私家車個(gè)體在2017年19日周四以及2017年20日周五的全天出行軌跡空間變化圖。從圖中可以看出該類型私家車個(gè)體工作日內(nèi)的軌跡是相似的，有極高的重現(xiàn)性。一般這種具有規(guī)律性的車輛出行時(shí)間集中在7:00—8:00以及18:00—19:00之間，即早晚交通高峰時(shí)期，首次出行后會(huì)在固定區(qū)域活動(dòng)或駐留較長(zhǎng)時(shí)間，再進(jìn)行1～2次折返，出行時(shí)間一般為1～2 h，并且這類私家車出行者占總體數(shù)目較多。總體來說，這類私家車用戶出行需求較大，出行時(shí)間較為穩(wěn)定，出行空間分布固定。

圖1 典型私家車Ⅰ軌跡空間Fig.1 Trajectory spaces of typical private car I

還有部分私家車工作日內(nèi)每天的軌跡與之前類型的私家車完全不同，這些私家車用戶的軌跡相比較下十分混亂，且經(jīng)過的卡口點(diǎn)位較高，軌跡路徑隨機(jī)性也較高，但是起訖點(diǎn)基本一致。這種情況的私家車軌跡顯示該車輛出行時(shí)間一般集中在早晨7:00—8:00或者下午13:00—14:00，被檢測(cè)的時(shí)間分布較為均勻，其在道路網(wǎng)絡(luò)的運(yùn)行時(shí)間較大，且逗留位置較多即出行軌跡存在較多的斷點(diǎn)，但駐留時(shí)間不長(zhǎng)，如圖2中車輛編號(hào)為1008533(a)和1021804(b)的私家車個(gè)體在2017年19日周四以及2017年20日周五的全天出行軌跡空間變化圖。從圖上可看出這類私家車用戶出行時(shí)間較為不固定，出行的空間分布較為紊亂，隨機(jī)性高。但從大量類似出行特征的用戶軌跡統(tǒng)計(jì)發(fā)現(xiàn)這類用戶在道路網(wǎng)運(yùn)行時(shí)間較大，雖然出行路徑軌跡隨機(jī)性高重復(fù)率低，但存在區(qū)域集中的現(xiàn)象，說明該類型用戶擁有區(qū)域選擇偏好。

圖2 典型私家車Ⅱ軌跡空間展示圖Fig.2 Trajectory spaces of typical private car II

對(duì)于出租車軌跡研究，本研究隨機(jī)選取了數(shù)百輛出租車個(gè)體進(jìn)行軌跡信息可視化定性分析。發(fā)現(xiàn)出租車工作日內(nèi)出行軌跡空間分布極為廣泛，大多數(shù)出租車出行時(shí)間極早幾乎全天24 h都在路上運(yùn)行，因此在道路網(wǎng)絡(luò)運(yùn)行時(shí)間較長(zhǎng)。且出租車的軌跡空間分布不僅廣泛，一般跨越4～5組團(tuán)區(qū)域，而且軌跡路徑組成較為復(fù)雜、隨機(jī)性高。但根據(jù)大量出租車數(shù)據(jù)的統(tǒng)計(jì)分析可見大量的出租車駕駛員一般會(huì)在熟悉的區(qū)域跑動(dòng)，因此軌跡存在聚集的現(xiàn)象，如圖3(b)172639出租車一周工作日的軌跡多出現(xiàn)在南岸區(qū)、巴南區(qū)。由此可以推測(cè)部分出租車駕駛員出行具有個(gè)體偏好的活動(dòng)區(qū)域。

2 車輛出行特征指標(biāo)分析

從軌跡數(shù)據(jù)的定性分析可知軌跡數(shù)據(jù)特征中直接或間接得包含了此出行軌跡與其他出行軌跡的差異性和共同處。出行軌跡的特征指標(biāo)指對(duì)原始出行軌跡數(shù)據(jù)進(jìn)行分析研究后，從中得到的可以表征出行特點(diǎn)的一系列具有代表意義的向量值。在對(duì)出行數(shù)據(jù)研究過程中可以提取出很多特征屬性，其中包含一些目前研究不需要的指標(biāo)；特征指標(biāo)相互之間也會(huì)存在各種關(guān)聯(lián)，有的相互補(bǔ)充，有的相互依賴，有的相互重復(fù)。在研究過程中為了減少不必要的重復(fù)，減少建模的復(fù)雜性和可行性，還需要對(duì)數(shù)據(jù)特征指標(biāo)進(jìn)行分析，剔除冗余或與研究不太相關(guān)的特征指標(biāo)，以選取更加適合、完善的特征指標(biāo)序列進(jìn)行研究分析[12-14]。出行特征指標(biāo)提取的詳細(xì)研究分析過程如圖4流程所示。

圖4 特征指標(biāo)提取研究過程Fig.4 Study process of characteristic indicator extraction

出租車類型車輛較其他類型車輛在軌跡重復(fù)率指標(biāo)中總體較低，在網(wǎng)時(shí)間總體分布很高的特點(diǎn)。這是因?yàn)槌鲎廛囋诔鲂熊壽E上具有極大的隨機(jī)性，但基于大數(shù)據(jù)軌跡的統(tǒng)計(jì)分析，發(fā)現(xiàn)出租車出行總體分布會(huì)存在偏好區(qū)域，即出租車會(huì)在偏好的區(qū)域內(nèi)更頻繁地活動(dòng)。實(shí)際情況下，出租車在空載的情況下，一般會(huì)無目的地運(yùn)行在道路上尋找乘客，但在一段時(shí)間內(nèi)出租車需求少，出租車駕駛員會(huì)根據(jù)經(jīng)驗(yàn)前往乘客較多的區(qū)域。另一方面，不同區(qū)域的乘客雖然出發(fā)點(diǎn)和目的地不相同，但根據(jù)概率也會(huì)存在概率高的出發(fā)點(diǎn)以及目的地。因此提出“熱點(diǎn)活動(dòng)區(qū)域”出行特征，分析出租車群體的最大概率活動(dòng)區(qū)域。在研究車輛出行特征指標(biāo)時(shí)提出熱點(diǎn)活動(dòng)區(qū)域指標(biāo)，熱點(diǎn)活動(dòng)區(qū)域指標(biāo)主要針對(duì)軌跡的出現(xiàn)空間選擇偏好進(jìn)行刻畫，該指標(biāo)體現(xiàn)了車輛個(gè)體最大概率出現(xiàn)的區(qū)域。根據(jù)上述對(duì)出行特征指標(biāo)的差異性分析，最終選取出行頻次、在網(wǎng)時(shí)間、軌跡重復(fù)率、以及首次出行時(shí)段、出行結(jié)束時(shí)段等指標(biāo)變量構(gòu)成出行特征的指標(biāo)向量，見表1。

表1 出行特征指標(biāo)Tab.1 Travel characteristic indicators

3 車輛出行特征模式劃分

3.1 改進(jìn)CFSFDP聚類算法原理

CFSFDP聚類算法(clustering by fast search and find of density peaks)的基本步驟是：首先將多維數(shù)據(jù)集進(jìn)行距離和密度計(jì)算，再選取聚類中心，最后對(duì)數(shù)據(jù)集中非聚類中心點(diǎn)進(jìn)行歸類操作[15-16]。

(1)局部密度ρi

局部密度ρi具有兩種計(jì)算方式，即Cut-off kernel (截止內(nèi)核算法)和Gaussian kernel (高斯算法)，其中當(dāng)數(shù)據(jù)規(guī)模較大時(shí)，局部密度應(yīng)該當(dāng)用截止內(nèi)核算法，當(dāng)數(shù)據(jù)規(guī)模較小時(shí)，利用高斯算法計(jì)算局部密度ρi。

Cut-off kernel (截止內(nèi)核算法)：

(1)

其中函數(shù)

(2)

參數(shù)dc>0為截?cái)嗑嚯x(cutoff distance)，需聚類前研究設(shè)定。

由公式(1)可知，實(shí)際上ρi求得S數(shù)據(jù)集中與xi之間的距離小于截?cái)嗑嚯xdc(不包含xi本身)的個(gè)數(shù)。

Gaussian kernel (高斯算法)：

(3)

截止內(nèi)核算法為離散值而高斯算法為連續(xù)值，相對(duì)來說高斯算法更小可能產(chǎn)生沖突(即不同數(shù)據(jù)點(diǎn)具有相同的局部密度)。

(2)距離δi

距離δi的計(jì)算公式為：

(4)

ρq1≥ρq2≥...≥ρqN。

(5)

由計(jì)算公式可知，當(dāng)xi是具有局部密度的數(shù)據(jù)點(diǎn)時(shí)，δi表示xi與S數(shù)據(jù)集中距離最大的數(shù)據(jù)點(diǎn)的距離值；否則，δi表示在所有局部密度大于xi的數(shù)據(jù)點(diǎn)中，與xi距離最小的那些數(shù)據(jù)點(diǎn)的距離值。

圖5 決策圖實(shí)例及示意圖Fig.5 Example and schematic diagram of decision-making diagram

3.2 BP神經(jīng)網(wǎng)絡(luò)識(shí)別模型

人工神經(jīng)網(wǎng)絡(luò)(artificial neural network)，是由大量處理單元通過廣泛互聯(lián)而構(gòu)成的網(wǎng)絡(luò)體系。神經(jīng)網(wǎng)絡(luò)的快速發(fā)展使得其成為模式識(shí)別強(qiáng)有力的工具。神經(jīng)網(wǎng)絡(luò)主要分為前饋神經(jīng)網(wǎng)絡(luò)和循環(huán)神經(jīng)網(wǎng)絡(luò)[17]。本研究介紹的BP神經(jīng)網(wǎng)絡(luò)模型屬于多層前饋神經(jīng)網(wǎng)絡(luò)，能夠?qū)崿F(xiàn)從輸入到輸出的映射功能，數(shù)學(xué)理論證明三層的神經(jīng)網(wǎng)絡(luò)以任意精度逼近任何非線性連續(xù)函數(shù)，這使得BP神經(jīng)網(wǎng)絡(luò)具有較強(qiáng)的非線性映射能力。

BP神經(jīng)網(wǎng)絡(luò)可以由BP和前饋神經(jīng)網(wǎng)絡(luò)兩個(gè)部分組成[18]，BP指的是Back Propagation(逆向傳播)，即誤差的逆向傳播，以及信號(hào)的前向傳播。其學(xué)習(xí)規(guī)則是使用最速下降法，通過反向傳播來不斷調(diào)整網(wǎng)絡(luò)的權(quán)值和閾值，使網(wǎng)絡(luò)的誤差平方和最小，BP網(wǎng)絡(luò)得結(jié)構(gòu)示意圖見圖6。

3.3 基于組合模型的出行特征模式劃分

出行特征群體辨識(shí)主要由3部分內(nèi)容組成：(1)RFID軌跡數(shù)據(jù)的預(yù)處理與出行指標(biāo)參數(shù)的選取，主要包括對(duì)出行群體的軌跡分析、預(yù)判斷和出行指標(biāo)的提取；(2)基于出行特征指標(biāo)的聚類分析，主要是針對(duì)選取的出行特征指標(biāo)進(jìn)行基于峰值密度的聚類分析，再根據(jù)分類結(jié)果分析劃分群體的具體性質(zhì)；(3)出行特征群體的識(shí)別，主要介紹BP神經(jīng)網(wǎng)絡(luò)出行特征識(shí)別模型。出行特征群體的辨識(shí)是先分出租車、私家車、貨車3種車輛類型，分別進(jìn)行群體的劃分和識(shí)別訓(xùn)練過程，針對(duì)不同車輛類型的出行特征分布特點(diǎn)進(jìn)行出行特征群體的研究，每一種車輛類型選取的出行特征指標(biāo)因車輛類型的運(yùn)行特點(diǎn)和分析角度有所不同，本研究先將全面考慮所有出行特征可能的不同分布，再利用聚類思想將選取的出行特征指標(biāo)進(jìn)行研究分析，再進(jìn)行群體識(shí)別。具體的車輛出行特征辨識(shí)方法結(jié)構(gòu)如圖7所示。

圖6 構(gòu)建BP神經(jīng)網(wǎng)絡(luò)模型Fig.6 Building a BP neural network model

圖7 出行特征群體辨識(shí)方法結(jié)構(gòu)Fig.7 Structure of travel characteristic group recognition method

4 案例分析

4.1 數(shù)據(jù)來源

本研究試驗(yàn)分析數(shù)據(jù)來源于重慶市城市交通大數(shù)據(jù)工程技術(shù)研究中心。研究數(shù)據(jù)主要來源于道路安裝的RFID設(shè)備檢測(cè)器，數(shù)據(jù)類型包含車輛RFID檢測(cè)數(shù)據(jù)、RFID卡口經(jīng)緯度數(shù)據(jù)，詳見表2。

表2 GIS路網(wǎng)匹配后的車輛RFID數(shù)據(jù)Tab.2 RFID data of vehicle after GIS road network matching

數(shù)據(jù)覆蓋范圍主要為重慶市主城區(qū)，其中覆蓋北碚組團(tuán)、西永組團(tuán)、西彭組團(tuán)、蔡家組團(tuán)、沙坪壩組團(tuán)、空港組團(tuán)、悅來組圖、禮嘉組團(tuán)、人和組團(tuán)、觀音組團(tuán)、渝中組團(tuán)、大楊石組團(tuán)、南坪組團(tuán)、大渡口組團(tuán)、李家沱組團(tuán)、茶園組團(tuán)以及界石組團(tuán)。重點(diǎn)研究區(qū)域選定為機(jī)場(chǎng)路干線影響區(qū)，其實(shí)際地理位置位于人和組團(tuán)內(nèi)，覆蓋范圍實(shí)際是圍繞機(jī)場(chǎng)路以及機(jī)場(chǎng)路所能在交通方面直接影響的區(qū)域，詳細(xì)范圍以及范圍內(nèi)RFID卡口分布如圖8所示。各組團(tuán)區(qū)域劃分以及所在區(qū)域內(nèi)RFID卡口布設(shè)點(diǎn)位分布如圖9所示。

圖8 干線影響區(qū)范圍及RFID卡口Fig.8 Scope of mainline affected area and RFID gates

圖9 組團(tuán)及RFID卡口布設(shè)Fig.9 Grouping and layout of RFID gates

4.2 私家車出行群體辨識(shí)

由于智能交通的發(fā)展，新型智能交通打車APP的出現(xiàn)，使得越來越多的私家車多向性的發(fā)展為商用運(yùn)營(yíng)車輛。合理規(guī)劃管理商用運(yùn)營(yíng)車輛群體，可以更高效率地提高車輛的使用率，提高運(yùn)行效率。根據(jù)相關(guān)文獻(xiàn)可知，商用運(yùn)行車輛的在網(wǎng)時(shí)間一般小于出租車個(gè)體，但比通勤出行需求的車輛運(yùn)行時(shí)間長(zhǎng)，且出行頻次較高，同時(shí)運(yùn)行軌跡與出租車相似隨機(jī)性高，因此軌跡重復(fù)率指標(biāo)值較低。由此可知，選取的出行特征指標(biāo)可以滿足商用運(yùn)行車輛的識(shí)別要求。

私家車群體的試驗(yàn)對(duì)象來源于重慶市主城區(qū)2017年10月17號(hào)—10月22號(hào)的RFID數(shù)據(jù)。由于樣本數(shù)據(jù)過大會(huì)造成程序運(yùn)行需要大量時(shí)間，系統(tǒng)運(yùn)行負(fù)荷大易崩潰現(xiàn)象。因此經(jīng)過多次試驗(yàn)后，隨機(jī)選取26 165組完整有效處理后數(shù)據(jù)樣本，提取出行頻次、在網(wǎng)時(shí)間、軌跡重復(fù)率、首次出行時(shí)段、出行結(jié)束時(shí)段出行特征指標(biāo)進(jìn)行聚類分析。數(shù)據(jù)樣本詳見表3。

表3 私家車指標(biāo)數(shù)據(jù)樣本Tab.3 Samples of private car indicators

(1)私家車群體劃分

由前節(jié)模型說明可知，CFSFDP算法的輸入?yún)?shù)是指標(biāo)相互之間的距離。為便于計(jì)算在計(jì)算指標(biāo)之間的距離前，首先需要對(duì)特征向量組數(shù)據(jù)進(jìn)行編號(hào)，詳見表4。

表4 私家車出行特征指標(biāo)向量Tab.4 Private car travel characteristic indictor vectors

本章節(jié)使用歐式距離公式來計(jì)算指標(biāo)向量之間的距離，詳見式(1)。

(1)

聚類中心A：總包含元素2 669個(gè)，其中核心元素2 669，無噪聲元素。

聚類中心B：總包含元素12 491個(gè)，其中核心元素12 491，無噪聲元素。

聚類中心C：總包含元素11 005個(gè)，其中核心元素11 005，無噪聲元素。

圖10展示了不同類別數(shù)據(jù)在二維圖上的分布情況，需要特別說明的是，圖10中的X軸與Y軸僅代表一種二維空間中的度量，沒有特定的含義。有圖可見，改進(jìn)后的局部密度計(jì)算方法是的每個(gè)數(shù)據(jù)點(diǎn)的值符合實(shí)際，避免了錯(cuò)誤的分配，減少了私家車的噪聲數(shù)據(jù)，且不同類型數(shù)據(jù)存在明顯分離。

圖10 私家車CFSFDP算法分類數(shù)據(jù)二維展示Fig.10 Two-dimensional display of classification data obtained by CFSFDP algorithm

(2)聚類結(jié)果分析及群體定義

由聚類算法結(jié)果可知26 165輛私家車樣本被劃分為3類。為了精準(zhǔn)地認(rèn)識(shí)這3種類型的出行特征，需要對(duì)選取的5個(gè)出行特征指標(biāo)進(jìn)行詳細(xì)的分析。對(duì)聚類結(jié)果的3類私家車樣本分別從5個(gè)指標(biāo)展開對(duì)比分析研究不同的特征規(guī)律，可以準(zhǔn)確定義3類出行特征群體。本研究為比較分析私家車不同群體的特征差異性，繪制參數(shù)比較圖，詳見圖11～圖16。

圖11 私家車出行群體出行頻次比較Fig.11 Comparison of travel frequencies of travel groups

圖12 私家車出行群體軌跡重復(fù)率比較Fig.12 Comparison of travel group trajectory repetition rates

圖13 私家車出行群體在網(wǎng)時(shí)間比較Fig.13 Comparison of online time of travel groups

① 出行頻次

類型A：該分類有較高的出行頻次，基本集中3～5次/d，平均4次/d。

類型B：該分類有較低的出行頻次，基本集中2～3次/d，平均2次/d。

圖14 類型A出行時(shí)段Fig.14 Travel time of Type A

圖15 類型B出行時(shí)段Fig.15 Travel time of Type B

圖16 類型C出行時(shí)段Fig.16 Travel time of Type C

類型C：該分類在出行頻次上集中2～4次/d，平均3次/d，出行頻次分布處于中間水平。

② 軌跡重復(fù)率

類型A：該類型有較低的軌跡重復(fù)率，分布集中在[0.18，0.46]，平均在0.3左右，說明該類私家車駕駛員每天的行程不太一樣，具有多變性和偶然性。

類型B：該分類有較高的軌跡重復(fù)率，基本集中[0.4，0.83]，平均在0.7左右，據(jù)相關(guān)分析軌跡重復(fù)率在0.7左右，即代表該類私家車駕駛員具有極高的軌跡重現(xiàn)性。

類型C：該類型在軌跡重復(fù)率集中[0.27，0.7]，軌跡重復(fù)率處于中間水平。

③ 在網(wǎng)時(shí)間

類型A：該類型有較高的在網(wǎng)時(shí)間，分布集中在[3 609，15 340]，平均在7 503 s，即2 h左右。說明該類私家車駕駛員每天在行程上花費(fèi)大量的時(shí)間。

類型B：該分類有較低的在網(wǎng)時(shí)間，基本集中[1 296，4 642]，平均在2 688 s，即0.7 h左右。代表該類私家車駕駛員在行程上所費(fèi)時(shí)間較少，對(duì)于車輛的依賴性較低。

類型C：該類型在網(wǎng)時(shí)間指標(biāo)集中[2 387,6 773]，平均4 325 s，駕駛員在網(wǎng)時(shí)間處于中間水平。

④ 出行時(shí)段分布

類型A：該類型首次出行時(shí)段較高集中在07:00—09:00時(shí)段，出行結(jié)束時(shí)間基本集中在22:00—24:00 時(shí)段。其中00:00—05:00出行的車輛分布較其他類型多。

類型B：該類型首次出行時(shí)段較高集中在07:00—09:00 時(shí)段，出行結(jié)束時(shí)間基本集中在17:00—19:00 時(shí)段。

類型C：該類型首次出行時(shí)段較高集中在07:00—09:00 時(shí)段，出行結(jié)束時(shí)間基本集中在19:00—22:00 時(shí)段。

表5 私家車出行特征群體指標(biāo)Tab.5 Private car travel characteristic group indicators

類型A出行頻次高，軌跡分布較離散且在網(wǎng)時(shí)間較高的特點(diǎn)看，該類型接近于私家車商用的司機(jī)特性，例如(滴滴車等)，因此類型A在本研究定義為商用私家車出行群體。

類型B出行頻次低，平均2次/d，出行的軌跡重現(xiàn)性高且在網(wǎng)時(shí)間較短的特點(diǎn)看，該類型接近于通勤私家車司機(jī)的特性，例如(即基本上班和下班用車)，因此類型B在本研究定義為通勤私家車出行群體。

類型C各特征指標(biāo)都處于中間水平的特點(diǎn)看，該類型是介于類型A和類型B之間的，因此類型C在本研究定義為其他私家車出行群體。

⑤ 群體識(shí)別結(jié)果

私家車出行群體識(shí)別性能函數(shù)見圖17，網(wǎng)絡(luò)經(jīng)過設(shè)定的296次訓(xùn)練達(dá)到設(shè)定的期望誤差。私家車10 466 組測(cè)試數(shù)據(jù)識(shí)別結(jié)果中準(zhǔn)確識(shí)別出10 173組數(shù)據(jù)，將類型A識(shí)別為類型B的有7組、類型B識(shí)別為類型C的有89組數(shù)據(jù)，將類型A識(shí)別為類型C的有20組，類型C識(shí)別為類型A的有59組，將類型C識(shí)別為類型B的有113組數(shù)據(jù)，將類型B識(shí)別為類型A的有5組，因此測(cè)試集的識(shí)別準(zhǔn)確率高達(dá)97.2%。

圖17 私家車出行群體識(shí)別性能函數(shù)Fig.17 Performance function for private car travel group recognition

4.3 出租車出行群體辨識(shí)

出租車群體的試驗(yàn)對(duì)象來源于重慶市主城區(qū)2017年10月17日至10月22日的RFID數(shù)據(jù)，經(jīng)數(shù)據(jù)處理后，共獲取14 866組完整有效處理后數(shù)據(jù)樣本。根據(jù)前文對(duì)出行特征的差異分析，本研究最終選取軌跡重復(fù)率、在網(wǎng)時(shí)間、熱點(diǎn)活動(dòng)區(qū)域及干線影響區(qū)活動(dòng)偏好4個(gè)指標(biāo)構(gòu)成模型的特征向量。

(1)出租車群體劃分

將出租車出行特征群體分類提取的4種特征指標(biāo)(軌跡重復(fù)率、在網(wǎng)時(shí)間、熱點(diǎn)活動(dòng)區(qū)域及干線影響區(qū)活動(dòng)偏好)，進(jìn)行聚類分析，得到結(jié)果如下：

① 聚類中心F：總包含元素5 562個(gè)，其中核心元素5 199，噪聲元素有363。

② 聚類中心G：總包含元素9 304個(gè)，其中核心元素8 669，噪聲元素有635。

由圖18可知聚類的兩種群體具有部分噪聲點(diǎn)，這可能由于部分出行車輛出行獨(dú)特性造成，但是這種獨(dú)特體較少，且可有不同，不能代表為一類群體，只能證明存在特殊情況。

圖18 出租車分類數(shù)據(jù)二維展示Fig.18 Two-dimensional display of taxi classification data

(2)聚類結(jié)果分析及群體定義

① 在網(wǎng)時(shí)間

類型F：相比類型G，該分類有存在較多樣本的在網(wǎng)時(shí)間較短，從箱圖比較上看略微低于類型G，但總體分布相差并不大。

類型G：同理對(duì)比類型F，該分類總體上略微比類型F的在網(wǎng)時(shí)間長(zhǎng)。

出租車的兩種分類群體對(duì)于在網(wǎng)時(shí)間指標(biāo)上，相差不大都集中位于[14,20]區(qū)間。這與前章分析也較為符合，由于出租車時(shí)商業(yè)運(yùn)營(yíng)車輛，基本全天都處于路網(wǎng)運(yùn)營(yíng)，極少由于特殊原因在網(wǎng)時(shí)間短的，見圖19。

圖19 出租車出行群體在網(wǎng)時(shí)間比較Fig.19 Comparison of online time of taxi travel groups

② 軌跡重復(fù)率

類型F：對(duì)比類型G，該類型的軌跡重復(fù)率略低，且部分重復(fù)率<0.05的樣本比類型G多，分布集中在[0.203，0.278]區(qū)間，平均處于0.242左右，見圖20。

類型G：對(duì)比類型F，該分類軌跡重復(fù)率略高，基本集中[0.238，0.308]區(qū)間，平均在0.272左右，整體分布略高于類型F。類型F相對(duì)比類型G在軌跡重復(fù)率的特征上較低，說明類型G更偏向于在固定區(qū)域活動(dòng)，但總體都處于很低的狀況，這是由于出租車的運(yùn)營(yíng)性質(zhì)決定出租車軌跡的隨機(jī)性以及廣泛性，極少存在重復(fù)率很高或者重復(fù)率很低的情況，因此單獨(dú)來看也不足以劃分為一種常規(guī)群體。

圖20 出租車出行群體軌跡重復(fù)率比較Fig.20 Comparison of trajectory repetition rates of taxi travel groups

③ 熱點(diǎn)活動(dòng)區(qū)域

類型F：該類型出租車熱點(diǎn)活動(dòng)區(qū)域?yàn)榇髼钍M團(tuán)、沙坪壩組團(tuán)、李家沱組團(tuán)、大渡口組團(tuán)、空港組團(tuán)、西永組團(tuán)以及北碚組團(tuán)。其中該群體出租車的熱點(diǎn)活動(dòng)區(qū)域集中分布在大楊石組團(tuán)，占全部的61%，其次是李家沱組團(tuán)，占據(jù)了18%的比例；以及沙坪壩組團(tuán)，占比15%，見圖21。

類型G：該類型出租車熱點(diǎn)活動(dòng)區(qū)域?yàn)橛逯薪M團(tuán)、南坪組團(tuán)、觀音橋組團(tuán)、人和組團(tuán)、禮嘉組團(tuán)以及空港組團(tuán)。其中該群體出租車的熱點(diǎn)活動(dòng)區(qū)域集中分布在觀音橋組團(tuán)，占全部的68%；以及南坪組團(tuán)，占據(jù)了21%的比例，見圖22。

圖21 類型F熱點(diǎn)活動(dòng)區(qū)域分布Fig.21 Distribution of hotspot activity area type F

圖22 類型G熱點(diǎn)活動(dòng)區(qū)域分布Fig.22 Distribution of hotspot activity area Type G

類型G與類型F兩個(gè)分類在熱點(diǎn)活動(dòng)區(qū)域分布上有相交區(qū)域(觀音橋組團(tuán)、禮嘉組團(tuán))，但是在大部分區(qū)域上不是相同的。類型F的熱點(diǎn)活動(dòng)區(qū)域集中在重慶城區(qū)中部和北面，類型G則是集中在西面。

④ 干線影響區(qū)選擇偏好

類型F：干線影響區(qū)偏好小于類型G,其集中在[0.012,0.067]之間，平均值為0.037，見圖23。

圖23 出租車出行群體干線影響區(qū)選擇偏好比較Fig.23 Comparison of selection preferences of taxi travel groups affected by main line

類型G：該分類較多的經(jīng)過干線影響區(qū)，偏好指標(biāo)在[0.053,0.146]，分布高于類型F。

綜上所述，類型F與類型G在干線影響區(qū)選擇偏好特征上存在較大的差別，干線影響區(qū)選擇偏好雖然也存在加大的隨機(jī)性，與出租車的熱點(diǎn)活動(dòng)區(qū)域以及出租車駕駛員住址有著一定的關(guān)系。相臨近的區(qū)域偏向選擇干線影響區(qū)。

類型F出行區(qū)域集中在重慶中部以及南部，且干線影響區(qū)選擇偏好較小的特點(diǎn)看，該類型出租車較少經(jīng)過干線影響區(qū)，因此類型F在本研究定義為其他區(qū)域偏好出租車群體。

類型G出行區(qū)域集中在重慶西部，且干線影響區(qū)選擇偏好較大的特點(diǎn)看，該類型干線影響區(qū)的影響想對(duì)較大，因此類型G在本研究定義為干線影響區(qū)偏好出租車群體。出租車的出行特征群體分類主要存在的差異性主要體現(xiàn)在運(yùn)行活動(dòng)區(qū)域上，由表6分析可知出租車中類型G(干線影響區(qū)偏好群體)的活動(dòng)熱點(diǎn)區(qū)域集中在渝中組團(tuán)、南坪組團(tuán)、觀音橋組團(tuán)、人和組團(tuán)以及空港組團(tuán)。

表6 出租車出行特征群體指標(biāo)Tab.6 Taxi travel characteristic group indicators

(三)群體識(shí)別結(jié)果

出租車特征群體得識(shí)的性能函數(shù)見圖24，從圖中可以看出實(shí)際網(wǎng)絡(luò)經(jīng)1 000次可以較精準(zhǔn)達(dá)到期望目標(biāo)。出租車5 946組測(cè)試數(shù)據(jù)識(shí)別結(jié)果見圖25，其中準(zhǔn)確識(shí)別出5 897組數(shù)據(jù)，將類型F識(shí)別為類型G的有32組數(shù)據(jù)，將類型G識(shí)別為類型F的有17組，因此測(cè)試集的識(shí)別準(zhǔn)確率高達(dá)99.18%。

圖24 出租車出行群體識(shí)別性能函數(shù)Fig.24 Performance function for taxi travel group recognition

圖25 出租車特征群體識(shí)別結(jié)果Fig.25 Taxi characteristic group recognition result

5 結(jié)論

本研究主要開展軌跡數(shù)據(jù)驅(qū)動(dòng)的出行特征指標(biāo)分析以及出行特征群體辨識(shí)建模兩大部分研究工作。通過RFID軌跡數(shù)據(jù)，基于數(shù)據(jù)統(tǒng)計(jì)分析研究了出行特征指標(biāo)提取方法，由此建立車輛出行特征辨識(shí)指標(biāo)體系。針對(duì)私家車、出租車的出行特點(diǎn)進(jìn)行出行特征指標(biāo)的定制提取，采用基于密度峰值的聚類算法(CFSFDP)以及BP組合模型，對(duì)特征指標(biāo)進(jìn)行聚類分析，再集合BP神經(jīng)網(wǎng)絡(luò)算法，建立出行特征群體識(shí)別模型。以重慶主城區(qū)內(nèi)私家車、出租車為試驗(yàn)對(duì)象，基于RFID數(shù)據(jù)分別提取出行頻次，軌跡重復(fù)率、首次出行時(shí)間、出行結(jié)束時(shí)間等指標(biāo)進(jìn)行聚類分析和識(shí)別建模，最終實(shí)現(xiàn)車輛的不同出行特征模式的辨識(shí)，具有較好的應(yīng)用前景，具體體現(xiàn)在以下幾個(gè)方面：

(1)本研究對(duì)車輛出行特征進(jìn)行了詳細(xì)的研究及分類，能夠有效辨識(shí)不同的車輛群體，為下一步預(yù)測(cè)不同群體的出行時(shí)長(zhǎng)、路徑選擇，以及掌握車輛的OD出行等方面提供了數(shù)據(jù)支撐；

(2)城市交通流的變化規(guī)律，其實(shí)質(zhì)就是對(duì)城市居民出行特征的研究分析。不同的群體有著不同的出行特性，只有追蹤車輛出行軌跡的溯源，才能從根本上對(duì)交通流進(jìn)行管控，為解決干線擁堵問題奠定了重要基礎(chǔ)；

(3)探究車輛運(yùn)動(dòng)軌跡的溯源，挖掘車輛的出行特征，不僅是精準(zhǔn)掌握車輛出行OD和進(jìn)行合理的路網(wǎng)資源配置的前提，同時(shí)也能夠更加緊密地結(jié)合城市道路交通管控，為城市道路交通個(gè)性化智能管控提供理論依據(jù)。