# Discovering Public Transit Riders’ Travel Pattern from GPS Data: A Case Study in Harbin

Authors: Shi An; Lei Wang; Haiqiang Yang; Jian Wang
Journal: Journal of Sensors
 (2017)

---

## Abstract

This paper proposes a public transit riders’ travel pattern measuring method based on divided cells and public transit vehicle’s GPS data. The method consists of two parts: detecting urban origin and destination areas and measuring the public transit riders’ travel pattern. Moreover, a series of indicators are proposed to reflect the public transit riders’ travel pattern. A case study is carried out to evaluate the methods, which use the GPS data collected from taxis and buses in Harbin, China. The study is expected to provide a better understanding of public transit riders’ travel patterns.


---

## Body

## 1. Introduction

The dramatic increase of urban vehicles leads to many serious problems including traffic congestions, road accidents, and air pollutions, which become general conundrums in many metropolitans all over the world, or even in smaller cities. Public transit is considered to be one of the most effective solutions for these general conundrums. Public transit includes various services that provide mobility to the general public, including buses, trains, ferries, shared taxi, and their variations [1, 2]. Public transit has many obvious advantages, such as less expenses, more effective mobility, and saving travel time, which makes more and more urban residents turn to public transit service. The city travelers and commuters proportion of Beijing who takes public transit system (including bus, taxi, and rail transit) has continuously increased in last decade. And this number came up to 54.2% in 2015, which represents more than 15.5 million trips per day [3]. This trend also occurred in other metropolitans and smaller cities all over the world.Public transit riders usually exhibit a fixed travel pattern. That is to say, at a macrolevel, usually a fixed number of Origin-Destination pairs locate in same place of urban area, and the number of trips between these OD pairs stays steady every day [4–6]. The usage ratios of each kind of public transit stay steady every day. And the trips at morning or evening rush hours occupy a large proportion of the total trips every day [3]. On the other hand, at a microlevel, a single traveler moves from resident area to work place in the morning and moves back in the evening for each day. If the travel pattern of public transit riders could be identified, the urban public transit manager can benefit from it. For example, analyzing those riders who prefer to choose public transit vehicles rather than private cars helps transit authorities to improve the strategies and even make new policy to attract new riders [7]. With a better understanding of the transfer behavior of public riders, transit agencies can adjust the bus route to make transfer easier, which can enhance the riders’ satisfaction [8]. By calculating the shortest path lengths between all station pairs, the original-destination matrix, and trip lengths, transit agencies can develop fare change plans to manage demand or raise revenue [9, 10].The Original-Destination (OD) matrix is a typical representation of residents’ travel pattern, which reflects travel demand, trip generation, travel distribution, and so on. The traditional models of establishing OD matrix usually rely on travel behavior survey. In practice, household travel surveys are conducted in many countries [11]. However, the survey data contain many limitations and errors. For example, some metropolitans in Japan (Tokyo, Kyoto, Osaka, etc.) survey residents’ travel behavior every 10 years. Since the cities are growing rapidly, the survey data will definitely be out of date [11]. Furthermore, the sampling rate is usually very low, which brings sampling errors [12]. Meanwhile, many human factors may affect the accuracy of the OD matrix, such as, willfully filtering some trips, forgetfulness, and other related factors [13].Compared with traditional survey data, GPS data and smart card data exhibit wider coverage, lower cost, and higher accuracy. With the rapid development of data-based technology, various intelligent transportation systems are widely applied in public transit system. These systems could collect residents’ mobility data every day, including longitude, latitude, boarding time, and dropping off time [14]. In the last decade, various researches based on these data have been carried out, for example, mining urban recurrent congestion evolution patterns from GPS-equipped vehicle mobility data [14], comparing accessibility in urban slums using smart card and bus GPS data [15, 16], discovering functional zones using bus smart card data [17], and partitioning bus operating hours into time of day intervals based on bus GPS data [18], which makes the data based transportation research to be a hot spot of transportation field [19].GPS data and smart card data are usually collected from different subsystems of one whole intelligent transportation system, or even from different systems. This is because bus, taxi, and rail transit usually belong to different public transit companies [3]. Therefore, most previous data based researches into transit traveler behaviors utilize smart data [20–22] or GPS data [23, 24], respectively. Accordingly, bus and rail transit riders’ travel pattern can be generated [20–22] or travel behavior of taxi riders [23, 24], respectively. To the best of my knowledge, public transit riders’ travel pattern researches are barely utilizing both smart card and GPS data. However, smart card data only provide riders’ boarding and dropping off information and lack locating information [21, 22]. As a result, only approximate location can be acquired, which causes inaccuracy of origin and destination inference. Moreover, on the microlevel, GPS-equipped public transit vehicle riders’ trips are far less than other public transit vehicles [3]. And the insufficient sampling number will definitely lead to inaccuracy of trip distribution. In this light, the smart data and GPS data should be integrated to discover public transit riders’ travel pattern.The aim of this paper is to propose an effective method to explore the public transit riders’ travel pattern in an urban area. There are two subgoals identified: detecting urban origin and destination areas at a cell level and measuring public transit riders’ travel pattern.This paper is organized as follows. Section2 discusses the definition of cells and locating points that will be applied to this research. Section 3 describes the proposed urban public transit riders’ travel pattern measuring method. Section 4 applies the proposed public transit riders’ travel pattern measuring methodology using taxi and bus GPS data and the urban road network of Harbin. Section 5 provides conclusions and recommendations for future research.


## 2. Definition of Cells and Locating Points

In this part, we are going to define some parameters of cells and locating points. Firstly, the urban area is divided inton×m small cells with same size. Cx,y is one of these cells, where x=1,2,…,n and y=1,2,…,m.According to the taxi GPS dataset, there are four types of occupation status. In this paper, we are going to study the public transit riders’ travel pattern, so we defined two types of locating points according to the boarding and dropping off status, which are described as follows:Type 1. Ptaxi.boarding represents taxi vehicles’ locating points whose occupation status value is 768 (i.e. the occupation status is boarding).Type 2. Ptaxi.droppingoff represents taxi vehicles’ locating points whose occupation status value is 16640 (i.e. the occupation status is dropping off).The locating points whose occupation status value is 256 (represents the taxi vehicles being vacant) are of no use to the public transit riders’ travel pattern, so they are not utilized in our research.P
t
a
x
i
I
D
,
t
s represents a specific locating point of a taxi vehicle, where ID and ts are the Taxi ID and Timestamp from the taxi GPS datasets.Passenger only needs to touch smart card once while boarding bus in many cities of China, such as Guangzhou, Xi’an, and Harbin, whereas, in Beijing, passenger needs to touch smart card once again while dropping off. It is easy to know both boarding locating points and dropping off locating points in Beijing. In Harbin, the dropping off locating points can be inferred based on boarding points and time period. In order to simplify the method, in this paper, we use the most periodic commuters’ data from the bus datasets. That is to say, if a specific bus rider gets to work every morning and comes back home every evening, only this type of rider’s locating points will be included in our research. This type of bus rider generates two trips each day. In this light, we defined two types of bus locating points as follows.Type 1. Pbus.boarding represents boarding passengers’ locating points.Type 2. Pbus.droppingoff represents dropping off passengers’ locating points.P
b
u
s
I
D
,
t
s represents a specific locating point of a bus vehicle, where ID and ts are the Card ID and Timestamp from the bus GPS datasets.How to identify the boarding locating points and dropping off locating points from smart card data is described as follows.Step 1.
Extract the locating points with same card ID (i.e.,PbusIDi,ts1,PbusIDi,ts2,…,PbusIDi,tsn) from the whole smart dataset, where Pbus_iIDi,ts is the locating point whose card ID is IDi and n is the number of this kind of locating points.Step 2.
Extract the locating points whosen=2 and ts1∈07:00,10:00 and ts2∈15:00,18:00.Step 3.
LetPbus.boarding=PbusID,ts1 and Pbus.droppingoff=PbusID,ts2.


## 3. Methodology

The methodology for discovering public transit riders’ travel pattern is described in this section. Two stepwise methods are proposed to achieve the main goal, including detecting Origin and Destination areas, measuring the commuter pattern between each OD pair.


### 3.1. Detecting Origin and Destination Areas

In this part, we are going to detect Origin and Destination areas from public transit riders’ GPS data (i.e., taxi and bus GPS data). Because of the reasonable urban planning and constructions in recent decades, the urban area is usually divided into many different function zones. Therefore, all the origin points and destination points of passengers’ trips will be clustered into several origin areas and destination areas. In this light, we use cluster algorithm to detect origin areas and destination areas in urban area. Moreover, we do not know how many clusters are in advance, and this type of cluster is usually not spherically-shaped. Therefore, we apply a customized Density-Based Spatial Clustering of Applications with Noise (DBSCAN) algorithm to solve this problem.For a specific cellCx,y, we define a parameter Nx,y to present the locating point (i.e. Ptaxi and Pbus) number in this cell during a specific period P. Taking a 3×3 cells, for example (as shown in Figure 1(a)), in one period all cells’ parameters Nx,y can be easily calculated. For two specific cells Cx1,y1 and Cx2,y2, the distance between them distCx1,y1,Cx2,y2 is calculated as follows:(1)distCx1,y1,Cx2,y2=l·x1-x22+y1-y22,where l is the length of the cell. As shown in Figure 1(b), the distance between Ci-1,j and Ci-1,j+1 is l, and the distance between Ci-1,j-1 and Ci,j is l·2.Figure 1
Example of cells withN value.
(a)
(b)Some relative parameters are defined as follows:Objecto represents the cell.Core objectco: the cell satisfies N≥δN, where δN is the threshold of parameter N.ε-Neighborhood of a core object is the space within a radius ε centered at co.Figure2 illustrates the flow chart of the customized DBSCAN algorithm in this research. The Original DBSCAN algorithm is based on density. Comparatively speaking, the density reachable points in our customized algorithm are defined by N≥δN. And we define theminPts value as 1 in our paper.Figure 2
The flow chart of the customized DBSCAN algorithm.


### 3.2. Measuring Public Transit Riders’ Travel Pattern

The public transit riders’ travel pattern can be reflected by some indicators, like trip number between each OD pair, the proportion of different transit, path between each OD pair, travel time of the path, and so on. In order to show these indicators, we defined some preliminary terms as follows.T
r
m
o
d
e
C
L
i
→
C
L
j is the trip number of a specific transit mode (i.e. taxi or bus) from ith cluster to jth cluster.Based on clusters of origin and destination area, given a specific periodP, the public transit OD matrix can be calculated.p
m
o
d
e
C
L
i
→
C
L
j is the proportion of a specific transit mode (i.e. taxi or bus) from ith cluster to jth cluster, which is calculated as follows:(2)pmodeCLi⟶CLj=TrmodeCLi→CLjTrMODECLi→CLj×100%,where TrMODECLi→CLj is the trip number of all transit modes from ith cluster to jth cluster.T
t
m
o
d
e
C
L
i
→
C
L
j is the travel time by the specific transit mode (i.e. taxi or bus) from ith cluster to jth cluster.


## 3.1. Detecting Origin and Destination Areas

In this part, we are going to detect Origin and Destination areas from public transit riders’ GPS data (i.e., taxi and bus GPS data). Because of the reasonable urban planning and constructions in recent decades, the urban area is usually divided into many different function zones. Therefore, all the origin points and destination points of passengers’ trips will be clustered into several origin areas and destination areas. In this light, we use cluster algorithm to detect origin areas and destination areas in urban area. Moreover, we do not know how many clusters are in advance, and this type of cluster is usually not spherically-shaped. Therefore, we apply a customized Density-Based Spatial Clustering of Applications with Noise (DBSCAN) algorithm to solve this problem.For a specific cellCx,y, we define a parameter Nx,y to present the locating point (i.e. Ptaxi and Pbus) number in this cell during a specific period P. Taking a 3×3 cells, for example (as shown in Figure 1(a)), in one period all cells’ parameters Nx,y can be easily calculated. For two specific cells Cx1,y1 and Cx2,y2, the distance between them distCx1,y1,Cx2,y2 is calculated as follows:(1)distCx1,y1,Cx2,y2=l·x1-x22+y1-y22,where l is the length of the cell. As shown in Figure 1(b), the distance between Ci-1,j and Ci-1,j+1 is l, and the distance between Ci-1,j-1 and Ci,j is l·2.Figure 1
Example of cells withN value.
(a)
(b)Some relative parameters are defined as follows:Objecto represents the cell.Core objectco: the cell satisfies N≥δN, where δN is the threshold of parameter N.ε-Neighborhood of a core object is the space within a radius ε centered at co.Figure2 illustrates the flow chart of the customized DBSCAN algorithm in this research. The Original DBSCAN algorithm is based on density. Comparatively speaking, the density reachable points in our customized algorithm are defined by N≥δN. And we define theminPts value as 1 in our paper.Figure 2
The flow chart of the customized DBSCAN algorithm.


## 3.2. Measuring Public Transit Riders’ Travel Pattern

The public transit riders’ travel pattern can be reflected by some indicators, like trip number between each OD pair, the proportion of different transit, path between each OD pair, travel time of the path, and so on. In order to show these indicators, we defined some preliminary terms as follows.T
r
m
o
d
e
C
L
i
→
C
L
j is the trip number of a specific transit mode (i.e. taxi or bus) from ith cluster to jth cluster.Based on clusters of origin and destination area, given a specific periodP, the public transit OD matrix can be calculated.p
m
o
d
e
C
L
i
→
C
L
j is the proportion of a specific transit mode (i.e. taxi or bus) from ith cluster to jth cluster, which is calculated as follows:(2)pmodeCLi⟶CLj=TrmodeCLi→CLjTrMODECLi→CLj×100%,where TrMODECLi→CLj is the trip number of all transit modes from ith cluster to jth cluster.T
t
m
o
d
e
C
L
i
→
C
L
j is the travel time by the specific transit mode (i.e. taxi or bus) from ith cluster to jth cluster.


## 4. Case Study

We apply the proposed methods to Harbin city (China). First of all, the datasets used in this case study are described, and the stepwise methods are implemented one by one.


### 4.1. Datasets

The smart card data and GPS data are collected from all operating buses and taxis in Harbin. Harbin is the capital of Heilongjiang province in northeast of China. With 4.74 million populations and 7,086 km2 areas in urban area [25], Harbin is a typical developing city in China. The public transit system of Harbin consists of taxi, bus, metro, and other rail transit. In addition, all taxis are equipped with GPS device and all buses are equipped with smart card system, which make collecting GPS data and smart card data possible.Approximately 16,000 operating taxis equipped with GPS device are running around Harbin’s urban area day and night. The location information is uploaded to the management system every 30 s during the day and 2 min at night. The data are accumulated to 2G in size and around 25 million rows each day. The taxi GPS data collected from 3rd Aug. to 7th Aug. 2015 is used, consisting of taxi ID, timestamp, latitude, longitude, and status, as shown in Table1. There are 4 kinds of “status” in the table: 17152, 16640, 256, and 768, which represent occupation, dropping off, vacant, and boarding, respectively.Table 1
Typical GPS records of a taxi.
Taxi ID
Timestamp
Latitude
Longitude
Status
0100320311
2015-08-21 17:08:21
45.726025
126.69883
17152
0100320311
2015-08-21 17:08:51
45.725937
126.70046
17152
0100320311
2015-08-21 17:09:21
45.72592
126.70041
16640
0100320311
2015-08-21 17:09:51
45.72424
126.70028
256
0100320311
2015-08-21 17:10:21
45.727768
126.678665
256
0100320311
2015-08-21 17:10:51
45.72415
126.70134
768
0100320311
2015-08-21 17:11:21
45.7243
126.70027
17152
Note. The columns from the left to right, respectively, represent the taxi identification; the taxi sampling time; the latitude and longitude of the position; and the occupation status.There are nearly 1500 buses traveling around urban area of Harbin, and these buses belong to 100 routes. The sampling frequency is 30 s. The bus IC records collected from 3rd Aug. to 7th Aug. in 2015 are used, consisting of Route ID, Bus ID, Card ID, Timestamp, Latitude, and Longitude, as shown in Table2.Table 2
Typical Bus GPS records of a bus.
Route ID
Bus ID
Card ID
Timestamp
Latitude
Longitude
120
0100320311
313011090
2015-08-11 08:08:21
45.739658
126.643119
120
0100320311
313013116
2015-08-11 08:08:24
45.739658
126.643119
120
0100320311
313008101
2015-08-11 08:08:28
45.739658
126.643119
120
0100320311
313191231
2015-08-11 08:08:33
45.73966
126.643119
Note. The columns from the left to right, respectively, represent the bus route number; bus ID; smart card ID; the card sampling time; the latitude and longitude of the position.In our research, another dataset is digital map of Harbin, which consists of most urban areas of Harbin. The research area is nearly 100 square kilometers and just covers the range of 2nd Ring Road of Harbin, as shown in Figure3. About 80 percent of GPS points locate in this research area. We divided this area into 250 square cells with same size. And each of the cell is 200×200 meters square, which is shown in Figure 3.Figure 3
Illustration of divided cells.The base map is OpenStreetMap.


### 4.2. Detecting Origin and Destination Areas

In this paper, we only measure the public transit riders’ travel pattern during working days. According to Section3.1, we can calculate the Nx,y values of all cells. Taking one week’s dropping off data (consisting of Ptaxi.droppingoff and Pbus.droppingoff) for example, that is, from 3rd Aug. to 7th Aug. in 2015, the relative data are illustrated in Table 3. And the sampling time is from 07:00 am to 19:00 pm.Table 3
The dropping off data in sampling 5 days.
Date
Total
Average
Minimum
Maximum
20150803
1,121,198
448.5
0
26,371
20150804
1,112,459
445.0
0
24,508
20150805
1,112,058
444.8
0
25,100
20150806
1,111,987
444.8
0
24,766
20150807
1,110,522
444.2
0
24,958
Sum
1,113,645
445.5
0
25,140
Note. Total is the number of dropping off points in all cells, average is the average number of dropping off points in one cell, minimum is the minimum number of dropping off points in one cell, and maximum is the maximum number of dropping off points in one cell.And then we apply the customized DBSCAN algorithm to measure the origin and the destination clusters. We setε value as 282.8 meters (i.e. 2·l, where l is 200 meters). We apply different δN values to our experiment, in order to find out the optimal value. Table 4 illustrates the clustering results by computing with different δN values. According to the results, when the value is larger than 1500, the cluster number is going to be stable. So we set the δN value as 2,000 each day and 10,000 for 5 days. In this light, we can measure the origin clusters (i.e. boarding clusters) and destination clusters (i.e. dropping off clusters) for the 5 days, as shown in Figure 4.Table 4
Dropping off clustering results by differentδN variously.
Date
Cluster number
δ
N
=
500
δ
N
=
1000
δ
N
=
1500
δ
N
=
2000
δ
N
=
2500
δ
N
=
3000
20150803
38
71
49
30
24
19
20150804
36
75
44
27
20
21
20150805
32
70
47
29
18
19
20150806
35
77
46
29
21
17
20150807
41
75
50
28
20
18
Average
36.4
73.6
47.2
28.6
20.6
18.8
Variance
11.3
8.8
5.7
1.3
4.8
2.2Figure 4
Origin and destination clusters results.
(a)
Origin clusters (boarding clusters)
(b)
Destination clusters (dropping off clusters)


### 4.3. Measuring Public Transit Riders’ Travel Pattern

For a specific OD pair, we can measure the travel pattern between them, that is,TrmodeCLi→CLj, pmodeCLi→CLj, and TtmodeCLi→CLj, as mentioned in Section 3.2. Taking a 20×17 cells area, for example, as shown in Figure 5, there are one typical origin cluster and one typical destination cluster. And the three indicators’ values in the survey 5 days can be calculated, which are shown in Table 5.Table 5
Results of the three indicators.
Taxi
Bus
T
r
367
509
p
41.9%
58.1%
T
t
16 min
40 minFigure 5
Illustration of a20×17 cells area.The Origin Cluster in this case study is the cellC(48,22), and the Destination Cluster consists of five cells (i.e. C(35,33), C(35,34), C(36,34), C(36,35), and C(36,36)). The cell C(48,22) is located nearby Harbin Hongqi Resident District, Heping Resident District, Yuanda Central Park Resident District, and so on. And there is a bus station locating in this cell, which consists of 12 bus routes (e.g. 14, 25, 31, and 44). And the Destination Cluster is located nearby Harbin International Golf Club, Harbin Wanda Plaza, several banks, and so on. This cluster area locates in an important commercial district in Harbin and definitely attracts many residents’ travels. In this Cluster, there are several bus stations consisting of more than 20 bus routes (e.g. 17, 27, 34, and 71). From the Origin Cluster to Destination Cluster, there are at least three direct bus routes (i.e. 71, 82, and 209).


## 4.1. Datasets

The smart card data and GPS data are collected from all operating buses and taxis in Harbin. Harbin is the capital of Heilongjiang province in northeast of China. With 4.74 million populations and 7,086 km2 areas in urban area [25], Harbin is a typical developing city in China. The public transit system of Harbin consists of taxi, bus, metro, and other rail transit. In addition, all taxis are equipped with GPS device and all buses are equipped with smart card system, which make collecting GPS data and smart card data possible.Approximately 16,000 operating taxis equipped with GPS device are running around Harbin’s urban area day and night. The location information is uploaded to the management system every 30 s during the day and 2 min at night. The data are accumulated to 2G in size and around 25 million rows each day. The taxi GPS data collected from 3rd Aug. to 7th Aug. 2015 is used, consisting of taxi ID, timestamp, latitude, longitude, and status, as shown in Table1. There are 4 kinds of “status” in the table: 17152, 16640, 256, and 768, which represent occupation, dropping off, vacant, and boarding, respectively.Table 1
Typical GPS records of a taxi.
Taxi ID
Timestamp
Latitude
Longitude
Status
0100320311
2015-08-21 17:08:21
45.726025
126.69883
17152
0100320311
2015-08-21 17:08:51
45.725937
126.70046
17152
0100320311
2015-08-21 17:09:21
45.72592
126.70041
16640
0100320311
2015-08-21 17:09:51
45.72424
126.70028
256
0100320311
2015-08-21 17:10:21
45.727768
126.678665
256
0100320311
2015-08-21 17:10:51
45.72415
126.70134
768
0100320311
2015-08-21 17:11:21
45.7243
126.70027
17152
Note. The columns from the left to right, respectively, represent the taxi identification; the taxi sampling time; the latitude and longitude of the position; and the occupation status.There are nearly 1500 buses traveling around urban area of Harbin, and these buses belong to 100 routes. The sampling frequency is 30 s. The bus IC records collected from 3rd Aug. to 7th Aug. in 2015 are used, consisting of Route ID, Bus ID, Card ID, Timestamp, Latitude, and Longitude, as shown in Table2.Table 2
Typical Bus GPS records of a bus.
Route ID
Bus ID
Card ID
Timestamp
Latitude
Longitude
120
0100320311
313011090
2015-08-11 08:08:21
45.739658
126.643119
120
0100320311
313013116
2015-08-11 08:08:24
45.739658
126.643119
120
0100320311
313008101
2015-08-11 08:08:28
45.739658
126.643119
120
0100320311
313191231
2015-08-11 08:08:33
45.73966
126.643119
Note. The columns from the left to right, respectively, represent the bus route number; bus ID; smart card ID; the card sampling time; the latitude and longitude of the position.In our research, another dataset is digital map of Harbin, which consists of most urban areas of Harbin. The research area is nearly 100 square kilometers and just covers the range of 2nd Ring Road of Harbin, as shown in Figure3. About 80 percent of GPS points locate in this research area. We divided this area into 250 square cells with same size. And each of the cell is 200×200 meters square, which is shown in Figure 3.Figure 3
Illustration of divided cells.The base map is OpenStreetMap.


## 4.2. Detecting Origin and Destination Areas

In this paper, we only measure the public transit riders’ travel pattern during working days. According to Section3.1, we can calculate the Nx,y values of all cells. Taking one week’s dropping off data (consisting of Ptaxi.droppingoff and Pbus.droppingoff) for example, that is, from 3rd Aug. to 7th Aug. in 2015, the relative data are illustrated in Table 3. And the sampling time is from 07:00 am to 19:00 pm.Table 3
The dropping off data in sampling 5 days.
Date
Total
Average
Minimum
Maximum
20150803
1,121,198
448.5
0
26,371
20150804
1,112,459
445.0
0
24,508
20150805
1,112,058
444.8
0
25,100
20150806
1,111,987
444.8
0
24,766
20150807
1,110,522
444.2
0
24,958
Sum
1,113,645
445.5
0
25,140
Note. Total is the number of dropping off points in all cells, average is the average number of dropping off points in one cell, minimum is the minimum number of dropping off points in one cell, and maximum is the maximum number of dropping off points in one cell.And then we apply the customized DBSCAN algorithm to measure the origin and the destination clusters. We setε value as 282.8 meters (i.e. 2·l, where l is 200 meters). We apply different δN values to our experiment, in order to find out the optimal value. Table 4 illustrates the clustering results by computing with different δN values. According to the results, when the value is larger than 1500, the cluster number is going to be stable. So we set the δN value as 2,000 each day and 10,000 for 5 days. In this light, we can measure the origin clusters (i.e. boarding clusters) and destination clusters (i.e. dropping off clusters) for the 5 days, as shown in Figure 4.Table 4
Dropping off clustering results by differentδN variously.
Date
Cluster number
δ
N
=
500
δ
N
=
1000
δ
N
=
1500
δ
N
=
2000
δ
N
=
2500
δ
N
=
3000
20150803
38
71
49
30
24
19
20150804
36
75
44
27
20
21
20150805
32
70
47
29
18
19
20150806
35
77
46
29
21
17
20150807
41
75
50
28
20
18
Average
36.4
73.6
47.2
28.6
20.6
18.8
Variance
11.3
8.8
5.7
1.3
4.8
2.2Figure 4
Origin and destination clusters results.
(a)
Origin clusters (boarding clusters)
(b)
Destination clusters (dropping off clusters)


## 4.3. Measuring Public Transit Riders’ Travel Pattern

For a specific OD pair, we can measure the travel pattern between them, that is,TrmodeCLi→CLj, pmodeCLi→CLj, and TtmodeCLi→CLj, as mentioned in Section 3.2. Taking a 20×17 cells area, for example, as shown in Figure 5, there are one typical origin cluster and one typical destination cluster. And the three indicators’ values in the survey 5 days can be calculated, which are shown in Table 5.Table 5
Results of the three indicators.
Taxi
Bus
T
r
367
509
p
41.9%
58.1%
T
t
16 min
40 minFigure 5
Illustration of a20×17 cells area.The Origin Cluster in this case study is the cellC(48,22), and the Destination Cluster consists of five cells (i.e. C(35,33), C(35,34), C(36,34), C(36,35), and C(36,36)). The cell C(48,22) is located nearby Harbin Hongqi Resident District, Heping Resident District, Yuanda Central Park Resident District, and so on. And there is a bus station locating in this cell, which consists of 12 bus routes (e.g. 14, 25, 31, and 44). And the Destination Cluster is located nearby Harbin International Golf Club, Harbin Wanda Plaza, several banks, and so on. This cluster area locates in an important commercial district in Harbin and definitely attracts many residents’ travels. In this Cluster, there are several bus stations consisting of more than 20 bus routes (e.g. 17, 27, 34, and 71). From the Origin Cluster to Destination Cluster, there are at least three direct bus routes (i.e. 71, 82, and 209).


## 5. Conclusion

This paper presented a cell-based urban public transit riders’ travel pattern measurement method. The method used GPS-equipped public transit vehicle’s locating data, which is more realistic and easy to collect. We proposed a customized DBSCAN algorithm to detect the origin and destination areas. We computed three indicators for each OD pair, which can reflect the relationship between the origin area and destination area. We carried out a numerical case study to evaluate our method, which uses taxi and bus GPS data in Harbin, China. The results can reflect OD pairs and relationship between each OD pair.In the future, we will improve the customized DBSCAN algorithm by accuracy and efficiency. The travel pattern is a little simple in this paper, so we will further discover some other indicators in depth. And we plan to extend this research by utilizing more kinds of GPS data from different public transit vehicles.



---

*Source: 5290795-2017-07-02.xml*