馬金麟，張宗博，謝君平，張其強

(江蘇大學(xué) 汽車與交通工程學(xué)院， 江蘇 鎮(zhèn)江 212013)

?

基于車牌識別數(shù)據(jù)的車輛OD矩陣獲取研究

馬金麟，張宗博，謝君平，張其強

(江蘇大學(xué) 汽車與交通工程學(xué)院， 江蘇 鎮(zhèn)江 212013)

利用高清智能卡口系統(tǒng)采集的車牌數(shù)據(jù)和調(diào)查得到的交通信息建立原始資料數(shù)據(jù)庫，提出車輛OD矩陣的獲取步驟，并利用VBA編程實現(xiàn)對Excel車牌數(shù)據(jù)文件的數(shù)據(jù)挖掘。提出了卡口間的車輛行程時間的獲取方法，運用VBA編程獲取全區(qū)所有的卡口間的行程時間，實現(xiàn)基于行程時間的卡口聚類劃分，從而完成對調(diào)查范圍的交通小區(qū)劃分。以某主城區(qū)為例，通過設(shè)定適當?shù)目陂g行程時間的劃分標準，實現(xiàn)對主城分區(qū)的交通小區(qū)劃分，并根據(jù)交通小區(qū)的車輛OD矩陣獲取流程圖得到車輛OD矩陣。

城市交通；OD矩陣獲取；車牌；行程時間；小區(qū)劃分

傳統(tǒng)的OD調(diào)查方法包括路邊詢問法、發(fā)表調(diào)查法、車輛牌照法、家訪法等。這類方法以抽樣調(diào)查為基礎(chǔ)，通過對調(diào)查后的數(shù)據(jù)進行處理，形成所需的各類OD表。該方法能夠獲取較為精確的OD數(shù)據(jù)，但需要投入大量的人力、物力。近年來許多交通領(lǐng)域的專家、學(xué)者改變了以往的交通調(diào)查理念和方法，而專注于對車牌照進行交通調(diào)查和分析[1]。目前，每天城市路網(wǎng)中的所有卡口設(shè)備記錄的車牌數(shù)據(jù)超過百萬計，但是現(xiàn)今對車牌信息的利用局限于車速測量[2]、出入收費管理[3]、違章處罰[4]等功能。車輛牌照的應(yīng)用局限于維護交通安全和城市治安，防止交通堵塞，沒有從系統(tǒng)全局的角度分析和利用車牌信息。以往的交通調(diào)查和分析的方法和理念具有一定的局限性和低效性，特別是起訖點調(diào)查，已經(jīng)不能滿足現(xiàn)今交通系統(tǒng)分析的需求。因此，針對以上問題，本文利用車牌照信息來獲取公路網(wǎng)的車輛OD矩陣，并用相關(guān)的理論和方法對獲取的OD矩陣進行相關(guān)的分析，這種基于車牌識別數(shù)據(jù)的車輛OD矩陣獲取使得大量車牌數(shù)據(jù)得到充分利用。這種基于車牌照的交通調(diào)查不僅能夠降低調(diào)查的成本，而且對調(diào)查結(jié)果的分析也更加準確和可靠。

1 基于車牌識別數(shù)據(jù)的車輛OD矩陣獲取方法的思路

1) 利用高清智能卡口系統(tǒng)獲取的車牌數(shù)據(jù)和調(diào)查得到的交通信息建立原始資料數(shù)據(jù)庫，提出車牌數(shù)據(jù)處理的步驟，并利用VBA編程實現(xiàn)數(shù)據(jù)挖掘。

2) 提出一種基于車輛行程時間的交通區(qū)劃分方法。通過對卡口之間車輛的行程時間的長短進行卡口的聚類，實現(xiàn)對交通小區(qū)的劃分，并在此基礎(chǔ)上獲取交通小區(qū)的車輛OD矩陣。

3) 對交通矛盾最突出的主城分區(qū)進行交通現(xiàn)狀分析。根據(jù)車輛的行程時間對主城區(qū)進行交通小區(qū)劃分，并在此基礎(chǔ)上獲取車輛OD矩陣表。

2 記錄車牌照法中的數(shù)據(jù)處理

2.1 建立原始資料數(shù)據(jù)庫

通常情況下，在調(diào)查中要求對每輛車記錄下列數(shù)據(jù)：調(diào)查地點、車牌號碼、車型、時間、行駛方向等。

2.1.1 車牌識別數(shù)據(jù)的獲取

通過卡口系統(tǒng)可以實現(xiàn)對路網(wǎng)中的設(shè)置點位全天候?qū)崟r地記錄車輛的出行信息，并將記錄的信息傳至終端的數(shù)據(jù)庫中進行存儲，系統(tǒng)可自動識別車牌號碼且全天車牌的識別準確率都在95%以上[5]。

對于車牌數(shù)據(jù)的輸出，后臺軟件平臺可根據(jù)輸入條件查詢數(shù)據(jù)。出于分析需要，可根據(jù)選擇的導(dǎo)出條件導(dǎo)出車牌數(shù)據(jù)，保存為Excel格式的數(shù)據(jù)文件。

2.1.2 卡口點位的獲取和編碼

調(diào)查地點的獲取對于車輛的OD矩陣的推斷至關(guān)重要。調(diào)查地點從某種意義上代表了機動車的起訖點。為了獲取調(diào)查地點的位置，采用卡口點位代替調(diào)查地點位置。本次實驗調(diào)查以常州市武進區(qū)為例進行說明，通過收集武進區(qū)的各種規(guī)劃資料和網(wǎng)絡(luò)信息繪制武進地圖，將武進區(qū)所有的卡口位置在地圖中標注，得到共計303個卡口點位。卡口點位分布情況見圖1。

圖1 卡口點位分布情況

建立原始數(shù)據(jù)庫后，借助編寫的程序，可以方便地進行數(shù)據(jù)處理。

2.2 車輛OD矩陣的獲取步驟

為了利用調(diào)查的數(shù)據(jù)獲取車輛出行OD矩陣，必須對原始數(shù)據(jù)庫進行數(shù)據(jù)處理。數(shù)據(jù)處理的目的就是找出車輛在直接影響區(qū)的行蹤。

車輛OD矩陣的獲取步驟：

1) 導(dǎo)出車牌數(shù)據(jù)

在數(shù)據(jù)導(dǎo)出界面中，選擇所需的數(shù)據(jù)導(dǎo)出條件，導(dǎo)出某工作日某時段(起始時間為5點，結(jié)束時間為9點)全區(qū)所有卡口的車牌數(shù)據(jù)。此時段的車輛出行行為多是從家出行至單位，除營運車輛外，多次、往返出行的車輛較少。

2) 獲取卡口點之間的車輛OD矩陣

由于導(dǎo)出的車牌數(shù)據(jù)是Excel文件，選擇利用VBA編程對車牌數(shù)據(jù)文件進行數(shù)據(jù)挖掘，獲取卡口與卡口之間的車輛數(shù)據(jù)信息。具體過程如下：

① 刪除無效數(shù)據(jù)

所導(dǎo)出的Excel格式的車牌數(shù)據(jù)文件，數(shù)據(jù)信息是保存在工作表sheet1中。數(shù)據(jù)形式如表1所示。首先，刪除sheet1中未被系統(tǒng)識別的數(shù)據(jù)(如表1中第4行數(shù)據(jù))和車牌號碼列中的唯一項，保留重復(fù)項。其次，在此基礎(chǔ)上刪除由前4位是數(shù)字和最后一位字母是“X”組成的車牌號碼*常州出租車的車牌號是由前四位數(shù)字和末位字母是“X”組成的統(tǒng)一格式，此步驟是剔除營運車輛的短途往返出行。。

表1 數(shù)據(jù)文件

② 獲取每輛車的起點和終點

以該時段某一輛車的出行情況說明車輛出行的起點和終點的含義，表2是車輛(浙K1779S)在路網(wǎng)行駛中，被高清智能卡口記錄到的行車信息。通過表2，可以在路網(wǎng)圖中繪制出該輛車在路網(wǎng)中的出行軌跡。如圖2所示。

表2 出行情況

圖2 出行軌跡

通過圖2和表2可知：在該段時間內(nèi)，這輛車的經(jīng)過所有卡口時刻的最大值和最小值分別為2015-12-20，08:54:22.708和2015-12-20，08:34:40.285，記最大值(2015-12-20，08:54:22.708)所對應(yīng)的卡口(長虹路-西園路)所在地即為該輛車出行的終點，最小值(2015-12-20，08:34:40.285)所對應(yīng)的卡口(S38常合高速常州南區(qū)收費站)所在地為該輛車出行的起點。

在步驟1的基礎(chǔ)上，篩選得到的每輛車經(jīng)過卡口時刻的最大值和最小值，以及時刻的最大值和最小值對應(yīng)的卡口名稱，保存至工作表sheet2中。

③ 計算各卡口之間的車輛OD矩陣

對于工作表sheet3，在A列中，從A2開始存放sheet2中車輛經(jīng)過卡口時刻最小值對應(yīng)的卡口名稱；在第1行中，從B1開始存放sheet2中車輛經(jīng)過卡口時刻最大值對應(yīng)的卡口名稱。列對應(yīng)的卡口名稱與行對應(yīng)的卡口名稱順序要一致。根據(jù)sheet2中的數(shù)據(jù)，在sheet3中行列交叉處填寫從起點至終點的車輛出行數(shù)量，對于不用的車牌號碼但是具有相同的起點和終點的卡口名稱要累加計數(shù)。至此便可得到卡口與卡口之間的OD矩陣。

3) 求解交通小區(qū)的車輛OD矩陣

以所劃分的交通小區(qū)內(nèi)含有的卡口點位為該交通小區(qū)的數(shù)據(jù)信息，即卡口為元素，交通小區(qū)就是由卡口為元素組成的集合。統(tǒng)計以卡口為元素的不同的交通小區(qū)集合，然后計算各交通小區(qū)集合之間的車輛數(shù)據(jù)信息，由此便可得到所需的車輛OD矩陣。

2.3 Excel VBA在數(shù)據(jù)處理中的應(yīng)用

2.3.1 Excel VBA

Visual Basic for Applications(VBA)是Visual Basic的一種宏語言，是微軟開發(fā)出來在其桌面應(yīng)用程序中執(zhí)行通用的自動化(OLE)任務(wù)的編程語言。

2.3.2 Excel VBA實現(xiàn)車輛OD矩陣獲取步驟

Excel VBA在獲取車輛OD矩陣中的應(yīng)用主要體現(xiàn)在卡口點之間的車輛OD矩陣的獲取過程中。在刪除無效數(shù)據(jù)的過程中，工作表sheet1中未被系統(tǒng)識別的數(shù)據(jù)和出租車的車牌號碼數(shù)據(jù)，可以利用Excel本身的數(shù)據(jù)篩選功能進行剔除。其余數(shù)據(jù)處理工作通過VBA編程得以實現(xiàn)。Excel VBA在實現(xiàn)卡口間OD矩陣獲取的程序設(shè)計分為兩步：

1) 原始數(shù)據(jù)的處理。根據(jù)vba基于ado訪問sql的策略，先建立自身的連接，然后構(gòu)造一個sql語句，取得車輛記錄大于1的所有車輛最大經(jīng)過時間對應(yīng)的記錄(車牌號碼Max，經(jīng)過時間Max，卡口名稱Max)，同理，通過構(gòu)造另外一個sql語句，并取得車輛記錄大于1的所有車輛最小經(jīng)過時間對應(yīng)的記錄(車牌號碼Min，經(jīng)過時間Min，卡口名稱Min)，將處理結(jié)果放入sheet2。

2) 卡口間OD矩陣的生成?；谇笆龅腶do處理方式，將sheet2的結(jié)果作為數(shù)據(jù)源，重新構(gòu)造sql，"select distinct 卡口名稱Max from[sheet2$]",提取卡口名稱Max唯一記錄放入首行，通過構(gòu)造sql，"select distinct 卡口名稱Min from[sheet2$] "，提取卡口名稱Min唯一記錄放入首列。根據(jù)構(gòu)造的這個矩陣，循環(huán)遍歷sheet2的所有記錄，根據(jù)卡口名稱Max在sheet3定位到列號，根據(jù)卡口名稱Min在sheet3定位到行號，通過此行號和列號，定位坐標，并記錄數(shù)據(jù)次數(shù)加1，整個sheet2處理完畢，即可得到sheet3所示的二維表單。

程序設(shè)計的流程如圖3所示。

圖3 程序設(shè)計的流程

3 基于車輛行程時間的交通小區(qū)劃分

3.1 卡口間行程時間的獲取

首先，利用卡口記錄的車牌數(shù)據(jù)獲取卡口的車輛的時間差。然后，對所獲取的時間差進行區(qū)間分段計數(shù)處理，通過時間差分段區(qū)間的車輛比例刻畫卡口間的車輛行程時間分布情況。

詳細的數(shù)據(jù)處理步驟如下：

1) 在原始車牌數(shù)據(jù)的基礎(chǔ)上，先篩選卡口，然后運行編寫好的宏命令，得到了所有車輛經(jīng)過卡口的時間差值，將結(jié)果輸出到工作表sheet2中。

2) 對時間差值進行分段，求出對應(yīng)時間差分段的車輛數(shù)、車輛比例。取車輛比例最高的時間差分段區(qū)間為卡口車輛行程時間。

3.2 VBA在基于車輛行程時間的交通區(qū)劃分中的應(yīng)用

3.2.1 Excel VBA實現(xiàn)行程時間數(shù)據(jù)處理

基于大數(shù)據(jù)量的考慮，利用excel ado 讀取表外文件的特性，把數(shù)據(jù)源單獨放入一個excel文件，一方面可以減少程序文件的體積，另一方面可以利用ado處理大數(shù)據(jù)的優(yōu)勢，不通過open的方式訪問工作薄，建立SQL查詢語句快速搜索相關(guān)符合要求的記錄。

本程序設(shè)計分為兩步：

第1步，通過構(gòu)造sql語句，"sql="select distinct 卡口名稱 from[源數(shù)據(jù)$]",提取卡口數(shù)據(jù)，采用遞歸算法Sub comb(a,arr,x As Long,y As Long,z As Long,jj As Long)，得到所有卡口組合。

第2步，直接用利用excel ado，通過構(gòu)造sql語句”sql="select 車號，datediff(""s"",min(left(經(jīng)過時間,len(經(jīng)過時間)-4))，max(left(經(jīng)過時間,len(經(jīng)過時間)-4))) as 間隔 from[源數(shù)據(jù)$] where 卡口名稱=′" & Cells(i,"A") & "′ or 卡口名稱=′" & Cells(i,"B") & "′ group by 車號"”，循環(huán)遍歷所有的卡口組合，實現(xiàn)卡口之間的時間差的統(tǒng)計，在data表中預(yù)制好公式(詳見表3)，通過在data工作表預(yù)制的公式，就可以在統(tǒng)計出時間差的同時，得出各區(qū)間時段的車輛數(shù)，繼而求出與各時段車輛數(shù)最大值相對應(yīng)的卡口間的行程時間。最后將每一組卡口間的行程時間保存于另一張工作表中。

用編寫的VBA程序計算車輛通過所有卡口之間的車輛行程時間，程序的流程圖見圖4，運行情況見圖5。

表3 Data表中的公式設(shè)計

圖4 計算車輛通過所有卡口之間的車輛 行程時間的程序流程

圖5 程序運行情況

表3中：第1列和第2列分別存儲經(jīng)過宏運行后獲取的每個卡口組合的車牌號碼和時間差；第3列(時間分段)，起始值為5，步長(區(qū)間長度)設(shè)置為15；第4列預(yù)制公式“=COUNTIFS(B:B,"<="&C2+15,B:B,">= "&C2)”，用于統(tǒng)計每個時間分段的車輛數(shù)；第4列的第一行預(yù)制公式“=OFFSET(D1,MATCH(MAX(E:E), E:E,0)-1,0)”，用于求出與各時段車輛數(shù)最大值相對應(yīng)的卡口間的行程時間。

3.2.2 Excel VBA實現(xiàn)基于行程時間的卡口聚類劃分

交通小區(qū)是具有一定交通關(guān)聯(lián)度和交通相似度的節(jié)點或連線的集合，隨時間、關(guān)聯(lián)度和相似度的變化而變化，反映城市路網(wǎng)交通特征的時空變化特性[6]。車輛的行程時間作為卡口劃分的依據(jù)，卡口間行程時間設(shè)置的大小直接影響到將所有卡口劃分的集合數(shù)量的多與少。根據(jù)不同的研究目的，可以設(shè)置不同的劃分標準，隨之調(diào)查范圍內(nèi)的所有卡口也將劃分為不同的集合，得到不同的劃分結(jié)果。

基于行程時間的卡口聚類劃分的VBA程序設(shè)計相對簡單。根據(jù)數(shù)據(jù)字典，通過關(guān)鍵字快速檢索數(shù)據(jù)的特點，引入字典變量，建立檢索庫、總庫字典和分庫字典。本程序設(shè)計分為兩步：

1) 將卡口1和卡口2 的所有組合讀入到檢索庫中，用于后期的數(shù)據(jù)快速讀取。循環(huán)遍歷表中所有卡口(卡口1、卡口2)，判斷此卡口是否存在于總庫中，存在即進入一下記錄的判斷。否則加入總庫，并建立對應(yīng)的分庫，并從當前記錄開始到最后記錄開始循環(huán)，判斷包含此卡口的所有記錄中，另一卡口和查詢卡口對應(yīng)分庫中所有組合，是否滿足檢索條件。若滿足即加入對應(yīng)的分庫，并加入總庫，當所有記錄循環(huán)完畢，即得到所有分庫。

2) 通過遍歷表中所有卡口組合，讀取每個分庫的所有卡口數(shù)據(jù)，建立矩陣，并通過先前的檢索庫，寫入對應(yīng)每個坐標數(shù)據(jù)。

通過循環(huán)輸出的所有分庫,即是輸出對應(yīng)于所有卡口劃分形成的集合，程序的流程見圖6。

圖6 基于行程時間的卡口聚類劃分的 VBA程序流程

最后，需要指出兩點：

1) 在進行實際的應(yīng)用時，若進行卡口劃分的行程時間標準值為m，只需將程序代碼第二行改成“Const minNum& =m”，可實現(xiàn)對城市路網(wǎng)中的卡口的自動聚類劃分，代碼截圖見圖7。

圖7 代碼截圖

2) 運用上述的方法進行交通區(qū)劃分，劃分的結(jié)果中會存在少數(shù)卡口同時滿足多個分庫的檢索條件，也即是說某些交通小區(qū)的集合會存在交集。完整的小區(qū)劃分是將所有的卡口劃分成互不相交的子集之并，所以要結(jié)合專家經(jīng)驗對運用上述方法劃分的交通小區(qū)進行優(yōu)化和修正[7]。

4 基于車牌識別數(shù)據(jù)的車輛OD矩陣獲取

以常州市武進區(qū)的主城區(qū)為例，簡述交通小區(qū)的車輛OD矩陣獲取過程。

4.1 車輛OD矩陣的獲取的思路

采用記錄車輛牌號法進行車輛OD調(diào)查時，機動車O、D點在調(diào)查中并不能直接得到。首先，按照車輛OD矩陣獲取步驟對所建立的原始數(shù)據(jù)庫進行數(shù)據(jù)處理，得到卡口與卡口之間的車輛OD矩陣。再結(jié)合對研究區(qū)域的交通小區(qū)的劃分情況，得到交通小區(qū)的車輛OD矩陣?；谲嚺谱R別數(shù)據(jù)的車輛OD矩陣獲取的過程如圖8所示。

圖8 交通小區(qū)的車輛OD矩陣獲取的流程

4.2 調(diào)查范圍的劃定

調(diào)查范圍的劃定是進行交通小區(qū)的劃分和交通小區(qū)的車輛OD矩陣獲取的前提條件，劃定調(diào)查區(qū)域范圍實際上就是確定境界線，區(qū)域的大小與交通規(guī)劃的目標是密切相關(guān)的。

通過對加權(quán)道路網(wǎng)絡(luò)的多重分形分析可知：主城分區(qū)的路網(wǎng)不僅是網(wǎng)絡(luò)密度最高，而且網(wǎng)絡(luò)覆蓋的深入程度、網(wǎng)絡(luò)分布的均勻程度最大，所以這個區(qū)域是交通矛盾最突出的區(qū)域。因此，將調(diào)查范圍劃定為：最北以312國道為界，南至西湖里；東起夏城路，西至龍江路，見圖9。

圖9 調(diào)查范圍

4.3 交通小區(qū)劃分

交通小區(qū)的劃分是建立交通模型的基礎(chǔ)，將調(diào)查范圍分為2個層次進行交通小區(qū)劃分：① 通過對卡口之間車輛的行程時間的長短進行卡口的聚類，從而實現(xiàn)對交通小區(qū)的劃分[8]。這一級的小區(qū)劃分是精細劃分，目的是為進行交通分配等微觀方面研究提供保障。② 在第1層次劃分的基礎(chǔ)上，根據(jù)土地利用、路網(wǎng)結(jié)構(gòu)、地形特征等情況對交通小區(qū)合并，形成交通大區(qū)，這一級別的劃分較粗，主要用于出行分布等較宏觀層次的研究工作。

第1層次的交通區(qū)劃分：基于車輛行程時間的交通區(qū)劃分。

首先，通過高清智能卡口系統(tǒng)的軟件平臺，獲取得到調(diào)查范圍內(nèi)2015-12-20日，5:00—9:00這個時間段內(nèi)所有的卡口車牌數(shù)據(jù)，利用上述車輛行程時間的獲取方法得到調(diào)查范圍內(nèi)所有卡口間的行程時間，結(jié)果如表4所示。

表4 卡口間行程時間

卡口名稱1卡口名稱2行程時間/sS38常合高速常州南(武進高新區(qū))收費站S39江宜高速武進65S38常合高速常州南(武進高新區(qū))收費站東龍路牛塘段305┋┋┋鳴新路—新平路龍江路—人民路275

其次，實現(xiàn)基于車輛行程時間的卡口聚類劃分。將已獲取得到的研究范圍內(nèi)卡口間的行程時間作為源數(shù)據(jù)，利用已經(jīng)編寫好VBA程序進行卡口的聚類劃分。為了能夠更好地對調(diào)查范圍內(nèi)的交通現(xiàn)狀進行分析，將聚類劃分的標準設(shè)置為60s，把程序代碼第二行改成“Const minNum& =60”，運行VBA程序即可實現(xiàn)對城市路網(wǎng)中的卡口的自動聚類劃分。通過劃分將主城分區(qū)的86個卡口按照車輛行程時間劃分為47個交通小區(qū)，編號為1～47。

第2層次的交通區(qū)劃分：結(jié)合主城區(qū)用地布局情況，調(diào)查范圍主要形成牛塘片區(qū)、湖塘片區(qū)、綜合中心區(qū)、智心片區(qū)、高新區(qū)等五大功能片區(qū)[9]。將已經(jīng)劃分好的交通小區(qū)合并到五大功能片區(qū)內(nèi)，形成5個交通大區(qū)，編號為1～5。

研究范圍劃分成直接影響區(qū)。為了體現(xiàn)過境交通的影響，在研究范圍周邊設(shè)虛擬小區(qū)。虛擬小區(qū)是將幾條主要的對外交通道路設(shè)置為界線形成的，對外交通道路之間所夾范圍內(nèi)含有的卡口作為該對外小區(qū)的含有的卡口。據(jù)此，對外小區(qū)劃為10個，編號為48～57。

本次交通區(qū)劃分如圖10、11所示。

圖10 交通小區(qū) 圖11 交通大區(qū)

4.4 車輛OD矩陣的獲取

根據(jù)已劃分好的交通分區(qū)，按照OD矩陣的獲取步驟，首先找到對應(yīng)分區(qū)內(nèi)的卡口編號(表4)，然后通過查表將各交通分區(qū)對應(yīng)的列編號和行編號一一列出來，最后將行編號和列編號分別輸入編好的源VBA程序文件中，便可以計算出對應(yīng)交通分區(qū)的車輛的OD矩陣表，從而計算得到主城分區(qū)交通大區(qū)的車輛OD矩陣表和交通小區(qū)的車輛OD矩陣表。

5 結(jié)束語

本文以常州武進區(qū)的主城區(qū)為例，在交通小區(qū)的劃分中，按照研究目的的不同，將調(diào)查區(qū)域分為2個層次進行交通小區(qū)的劃分：① 采用了基于車輛的行程時間對調(diào)查范圍進行了微觀層次的小區(qū)劃分；② 在第1層次劃分的基礎(chǔ)上，根據(jù)土地利用、路網(wǎng)結(jié)構(gòu)、地形特征等情況對交通小區(qū)合并成，實現(xiàn)調(diào)查范圍交通大區(qū)的劃分。根據(jù)交通小區(qū)的車輛OD矩陣獲取的流程，實現(xiàn)了對交通分區(qū)的車輛OD矩陣的獲取。

[1] 王龍飛.基于車牌照的車輛出行軌跡分析方法與實踐研究[D].西安：長安大學(xué)，2011.

[2] 王煒，過秀成.交通工程學(xué)[M].南京：東南大學(xué)出版社，2000.

[3] 趙丹，孫云蓮.基于ICA的車牌牌照字符識別系統(tǒng)[J].武漢大學(xué)學(xué)報(理學(xué)版)，2005,51(S2)：208-210.

[4] Nanne van der Zijpp.Dynamic OD-Matrix Estimation on Motorway Networks[D].Holand:Department of Civil Engineering Delft University of Technology,1996.

[5] 杭州?？低晹?shù)字技術(shù)股份有限公司.高清智能卡口系統(tǒng)方案[EB/OL].[2015-05-20].http://wenku.baidu.com/link?url=mFTDhZBbshS-YY1kILDtSdOc

Fu_HwNOMPQA2hxX3RVxcjB8tKddd4LRAT-XeCli8MJ

zBRayH6wSQP8s7UVSFG3q5Vz8WQGsgKOL5grQClzC.

[6] 蔡毅.基于出租車GPS數(shù)據(jù)的南京市交通小區(qū)的劃分[J].科技創(chuàng)新導(dǎo)報，2015(23)：218.

[7] 李曉丹，楊曉光，陳華杰.城市道路網(wǎng)絡(luò)交通小區(qū)劃分方法研究[J].計算機工程與應(yīng)用，2009，45(5)：19-22.

[8] 蔡毅.基于出租車GPS數(shù)據(jù)的南京市交通小區(qū)的劃分[J].科技創(chuàng)新導(dǎo)報，2015(23)：218.

[9] 南京大學(xué)城市與區(qū)域規(guī)劃系.武進中心城區(qū)空間發(fā)展戰(zhàn)略與優(yōu)化提升規(guī)劃[EB/OL].[2015-06-28].http://wj.czghj.gov.cn/info/WJ2335.html.

(責任編輯 楊黎麗)

Research on Obtaining the Vehicle OD Matrix Based on License Plate Number

MA Jinlin, ZHANG Zongbo, XIE Junping, ZHANG Qiqiang

(School of Automotive and Traffic Engineering, Jiangsu University, Zhenjiang 212013, China)

The original database based on theobtained license plate number and investigation on the traffic information was established. And then, data processing steps for the license plate were put forward with the use of VBA Excel programming license plate data files for data mining. The method of acquiring the travel time between bayonets is proposed by using the VBA programming to obtain all of the travel time between bayonets. The division of the traffic area through clustering analysis the length of the vehicle’s travel time between bayonets was achieved. Taking the main district of Wujin as an example, the main city zone is divided into pieces of traffic areas by setting the appropriate standard of travel time between bayonets. Obtaining OD matrix of vehicle is demonstrated by flowchart of achieving OD matrix of traffic zone.

urban traffic; obtain OD matrix; license plate; travel time; zone division

2016-09-11

國家自然科學(xué)基金資助項目(71101072)

馬金麟(1973—)，男，山東臨清人，碩士，副教授，主要從事交通運輸規(guī)劃與管理研究，E-mail:mjl@ujs.edu.cn。

馬金麟，張宗博，謝君平，等.基于車牌識別數(shù)據(jù)的車輛OD矩陣獲取研究[J].重慶理工大學(xué)學(xué)報(自然科學(xué))，2017(7):48-55.

format：MA Jinlin,ZHANG Zongbo,XIE Junping, et al.Research on Obtaining the Vehicle OD Matrix Based on License Plate Number[J].Journal of Chongqing University of Technology(Natural Science)，2017(7):48-55.

10.3969/j.issn.1674-8425(z).2017.07.007

U491

A

1674-8425(2017)07-0048-08