CN105654721B - A kind of method for calculating public transport OD amounts using bus IC card and vehicle GPS data - Google Patents
A kind of method for calculating public transport OD amounts using bus IC card and vehicle GPS data Download PDFInfo
- Publication number
- CN105654721B CN105654721B CN201610040284.6A CN201610040284A CN105654721B CN 105654721 B CN105654721 B CN 105654721B CN 201610040284 A CN201610040284 A CN 201610040284A CN 105654721 B CN105654721 B CN 105654721B
- Authority
- CN
- China
- Prior art keywords
- trip
- passenger
- bus
- card
- ith
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000004364 calculation method Methods 0.000 claims description 20
- 239000011159 matrix material Substances 0.000 claims description 14
- 238000011835 investigation Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 2
- 238000013500 data storage Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003032 molecular docking Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0108—Measuring and analyzing of parameters relative to traffic conditions based on the source of data
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0125—Traffic data processing
Landscapes
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Traffic Control Systems (AREA)
- Devices For Checking Fares Or Tickets At Control Points (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Navigation (AREA)
Abstract
A kind of method for calculating public transport OD amounts using bus IC card and vehicle GPS data, the present invention relates to the method for calculating public transport OD amounts using bus IC card and vehicle GPS data.The present invention is to solve existing method precision it is low long with poll cycle the problem of.Establish rider history trip information table；Determine information；It is poor to calculate adjacent charge time twice；Judge computing unit CI, i+1In whether go on a journey twice be that " travel route is identical；Judge to utilize secondary website and the position relationship of ith travel route of swiping the card of passenger k i+1s；Search meets the history trip information of following condition in Table1；Judge computing unit CI, i+1Whether middle i+1 time trip is that same day last time is gone on a journey；Judge whether passenger k is last passenger of the same day；Establish passenger's bus trip OD tables；Estimated according to it with the get on the bus get-off stop of passenger of station.The present invention is applied to field of urban traffic.
Description
Technical Field
The invention relates to a method for calculating the OD (origin-destination) quantity of a bus by adopting a bus IC (integrated circuit) card and vehicle-mounted GPS (global positioning system) data.
Background
Urban traffic jam causes serious troubles to resident travelling, urban public transport is developed vigorously to be an effective measure for relieving urban traffic jam, and a reasonable urban public transport planning and bus operation scheduling system can improve urban public transport attractiveness and ensure safe and smooth urban traffic. The public transport trip OD data are basic data of urban public transport planning and public transport vehicle operation scheduling, timely and accurate acquisition of the public transport trip OD data can optimize an urban public transport system, can also guide a public transport company to reasonably schedule public transport vehicles, and meets public transport trip requirements of urban residents.
The bus trip OD data are mainly obtained by a manual sampling investigation method at present, the precision of the manual sampling investigation method is limited by the sampling method and the sample size, the investigation precision is low, and the investigation period is long; the manual survey needs huge manpower and material resources, and the survey method is not economical; under the restriction of factors such as cultural level of the respondents, the data obtained by manual investigation is often incomplete or unreal. Therefore, the reliability of the public transportation OD data obtained by the manual investigation method is poor, so that unreasonable public transportation planning and public transportation vehicle operation scheduling are generated, and the purpose of relieving urban traffic jam by public transportation cannot be achieved.
With the advent of public transportation IC cards and vehicle-mounted GPS, public transportation IC card systems and public transportation GPS systems have been established in almost all cities. The bus-mounted GPS and IC card data contain the running information of the bus and the card swiping information of the passengers, and the inventor finds that a plurality of trip information of the passengers swiping the bus can be obtained by processing and calculating the bus-mounted GPS and IC card data for a long time.
Disclosure of Invention
The invention provides a method for calculating the bus OD amount by adopting bus IC card and vehicle-mounted GPS data, aiming at solving the problems of low precision and long investigation period of the existing method.
A method for calculating the bus OD quantity by adopting a bus IC card and vehicle-mounted GPS data is realized according to the following steps:
the method comprises the following steps: establishing a historical travel information Table Table1 of the passenger, and storing and updating the bus travel information of the bus card holding passenger within the last 30 days;
wherein, bus trip information includes: the method comprises the following steps that a passenger card number is held, boarding card swiping time is obtained, boarding card swiping stations are obtained, a passenger takes a bus line number, and calculated getting-off card swiping time and getting-off card swiping stations are obtained;
step two: reading IC card swiping data and vehicle-mounted GPS data of all card swiping passengers in one day, and determining information: the time of swiping the card when the passenger gets on the bus, the station of swiping the card when the passenger gets on the bus, the serial number of the bus line taken by the passenger and the driving direction;
establishing a card swiping information Table Table2 of the passengers for the same day to store the information, starting from the 1 st passenger in Table2, marking the kth passenger as a passenger k, and executing a step three, wherein k =1,2,3 and …;
step three: the method comprises the following steps that 1 st card swiping on the same day is started from a passenger k, each card swiping corresponds to one public transport trip, the starting point of the trip is a corresponding card swiping station, and the time difference delta T between two adjacent card swiping times is calculated:
if the delta T is less than or equal to 30s, combining the two trips into one trip;
if Δ T>, 30s, i (i =1,2,3, …, m) th trip and i +1 th trip as a computing unit C i，i+1 Executing the step four; m is the total times of card swiping of the passenger k on the same day;
step four: judgment calculation unit C i，i+1 Whether the two trips are two trips of 'the same trip route and different traveling directions' or not is determined:
if the judgment result is yes, the ith trip and the (i + 1) th trip of the passenger k are round trip trips, the terminal point of the ith trip is the starting point of the (i + 1) th trip, and the terminal points of the ith trip and the (i + 1) th trip of the passenger k are recorded into Table1;
if the judgment result is negative, executing the step five;
step five: judging whether the terminal point of the ith trip of the passenger k can be determined by utilizing the position relationship between the terminal point of the ith +1 card swiping of the passenger k and the ith trip line:
if the judgment result is yes, recording the ith trip information of the passenger k into Table1; if the judgment result is negative, executing the step six;
step six: reading the passenger cardholder number C of passenger k in Table2 k Ith driving routeDirection of travelTime for getting on bus and swiping cardAnd getting-on card swiping stationSearching Table1 for historical travel information satisfying the following conditions:
1) The boarding card swiping station is
2) The ith driving route is
3) The running direction is
4) The time for passenger to get on the bus and to swipe the card is
If the trips meeting the conditions are searched, the starting point and the ending point of the searched trip and the ith trip of the passenger k are consistent, the ith trip ending point of the passenger k is obtained according to the searched trip ending point, and is stored in Table1, and the seventh step is executed;
if the travel meeting the conditions cannot be searched, marking the ith travel endpoint of the passenger k in Table1 as null; i = i +1, and returns to step four;
step seven: judgment calculation unit C i，i+1 Whether the middle (i + 1) th trip is the last trip on the day:
if so, taking the last trip end point of the day as the 1 st trip starting point of the day, obtaining the all-day trip information of the passenger k, and executing the step eight;
if the judgment result is negative, i = i +1, the calculation unit C is switched to i+1，i+2 Returning to the step four;
step eight: determine whether passenger k is the last passenger of the day:
if the judgment result is yes, executing the step nine;
if the judgment result is negative, k = k +1, returning to the step three, and calculating the travel OD of the next passenger;
step nine: establishing a current passenger bus trip OD Table, converting an getting-on card swiping station and a getting-off station of each trip of the passengers in Table1 into the current passenger bus trip OD Table, updating Table1, and calculating the trip OD of each passenger through the steps to obtainAnd (3) an initial OD matrix, wherein the bus trip OD quantity taking p as a trip starting point and q as a trip end point is recorded as Trips pq Executing the step ten;
step ten: for the passenger travel OD without getting off the station, estimating according to the off station of the passenger getting on the same station
Trips pq *＝Trips pq ·b p
Wherein,
b p -calculation coefficient of trip quantity with trip starting point p
-total bus trip amount with trip starting point p
-the travel starting point is p, and the travel end point is determined as the OD value of the bus travel
Trips pq * And the final bus travel OD amount of the whole day with the travel starting point being p and the travel end point being q.
The invention has the following effects:
the invention provides a method for calculating a bus OD (origin-destination) by adopting IC (integrated circuit) card data and bus-mounted GPS (global positioning system) data. The IC card data and the vehicle-mounted GPS data can be directly obtained through the vehicle-mounted GPS and the IC card data storage equipment, huge manpower and material resources are not needed, and the calculation method is economical and feasible.
The bus boarding stop and the boarding time of a bus card holder are determined through the existing bus IC card swiping data and the vehicle-mounted GPS data. The method has the advantages that the positions and the time of get-off stops of passengers during each trip in one day are determined according to the interrelation among the card swiping data of the IC cards of the passengers in each time in one day and the historical trip rules of the passengers who hold the cards in combination with the vehicle operation data provided by the vehicle-mounted GPS, the trip OD of the passengers who hold the cards in transit can be obtained in real time, basic data is provided for urban traffic planning and urban bus driving arrangement, and the method has important significance for promoting reasonable and rapid development of urban traffic.
Drawings
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2 is a schematic diagram of a round trip route;
fig. 3 is a schematic diagram of two travel routes of the i +1 th travel starting point on the ith travel route;
fig. 4 is a schematic diagram of the ith travel route downstream station searched within the range of the i +1 th travel starting point 500 m.
Detailed Description
The first embodiment is as follows: in the method for calculating the bus OD amount by using the bus IC card and the vehicle-mounted GPS data according to the embodiment, the method is specifically executed according to the following steps:
1. establishing a historical travel information Table Table1 of the passenger, storing and updating the bus travel information of the bus card holder in the last 30 days, wherein the information comprises: the bus system comprises a passenger card number, the time of swiping the card for the passenger to get on the bus on a trip, a station of swiping the card, the serial number of the bus line taken by the passenger, and the calculated time of getting off the bus and the station of getting off the bus. The table format corresponding to each card number is shown in table 1:
table1 passenger history trip information table
2. Reading the IC card data and the vehicle-mounted GPS data of all the card-swiping passengers in one day, and determining the boarding time, boarding stations, riding routes and traveling directions of the passengers. The card swiping information Table Table2 of the current day of the passenger is established to store the information. Starting with the 1 st passenger in Table2, the k-th passenger is denoted as passenger k, k =1,2,3, ….
The IC card data of passenger k read in this step includes a card number C k Time of card swipingAnd the bus number of the card swipingStored in the card swiping information table of passenger k on the same day. According to vehicle number in bus GPS databaseSearching the GPS data of the vehicle, and determining the time of card swipingThe vehicle taken by the passenger k is obtained by comparing the GPS data of the vehicleAt the time of card swipingAt a docking stationAnd direction of travelThe numbers of 0 and 1 are respectively numbered in two directions of the reciprocating of any bus line, the stop of the bus at the time of swiping the card by the passenger is the getting-on stop of the passenger, and through the processing, the getting-on stop of the passenger k going out every time isAnd direction of travelDetermined, stored in Table 2. Table2 format is shown in Table 2:
table2 passenger card information table for the same day
3. Starting from the 1 st card swiping of the passenger k on the same day, each card swiping is recorded as a bus trip, and the i (i =1,2,3, …, m) th card swiping and the i +1 th card swiping are used as a computing unit C i，i+1 . Calculating C i，i+1 The time difference Delta T of two card swiping times. If DeltaT is less than or equal to 30s, then C is added i，i+1 The middle and twice trips are combined into one trip and transferred to the next calculation unit, if delta T>, 30s, the step four is executed.
4. Using a computing unit C i，i+1 The bus route numbers and the automobile traveling directions of the buses with the cards of the middle and two times judge whether the ith and (i + 1) th trips of the passengers are to-and-fro trips, and the trip route relationship is shown in fig. 2. If it is not The i-th trip and the (i + 1) -th trip of the passenger are the round trip, and the end point of the i-th trip is the starting point of the (i + 1) -th trip. The ith trip end point of passenger k is recorded in Table1. If it is usedOrIt cannot be determined whether the i-th and i + 1-th trips of the passenger are round trips. And F, turning to the step five.
5. And determining the ith trip destination of the passenger k by utilizing the position relationship between the ith +1 card swiping destination of the passenger k and the ith trip route. If so, the ith trip information of the passenger k is recorded in Table1. If the determination cannot be made, step six is executed.
In this step, a specific method for determining the terminal point of the ith trip of the passenger k by using the position relationship between the station point of the ith +1 card swiping of the passenger k and the ith trip route is as follows:
1): station for judging i +1 th card swiping time of passenger kWhether to punch the card at the ith time for the automobile running lineIf yes, executing the next step, and if no, turning to 3).
2): judgment ofWhether the bus is in the ith card swiping bus driving directionIf the judgment result is yes, the relationship between the ith and (i + 1) th trip routes is shown in FIG. 3, thenIs the ith trip end point. If the judgment result is negative, the next step is carried out.
3): judging stationWhether or not a line exists within a range of 500 meters aroundEdge ofIf the judgment result of the station in the downstream direction is yes, the relationship between the ith and (i + 1) th trip routes is shown in fig. 4, the next step is executed, and if the judgment result is no, the terminal point of the ith trip of the passenger k cannot be determined by using the position relationship between the station point of the ith +1 card swiping of the passenger k and the ith trip route。
4): judging whether the station searched in the step 3) is unique or not, if so, the station is the get-off station for the ith card swiping trip, and if not, selecting the distance stationAnd the nearest station is used as a station for getting off the vehicle in the next trip, and the terminal point of the ith trip of the passenger is obtained.
6. Reading card number C of passenger k in card swiping information table of passenger k on the same day k Ith card swiping bus line numberDirection of travelTime of card swipingAnd card swiping stationSearching the historical travel information meeting the following conditions in the first table:
1) The card swiping station is
2) The riding line is numbered as
3) The running direction is
4) Get on the bus at the time
If the trips meeting the conditions are searched, the searched trip is consistent with the starting and ending point of the ith trip of the passenger k, the terminal point of the ith trip of the passenger k is obtained according to the searched terminal point of the trip, and the terminal point is stored in Table1, and the step 7 is executed. If the travel meeting the condition cannot be searched, the ith travel endpoint of the passenger k is marked as null in Table1. i = i +1, step seven is performed.
Through the above steps, the passenger k calculation unit C i，i+1 The determination work of the ith trip end point is finished, at this time, if the ith getting-off point value in Table1 is the station number, the get-off station of the ith bus trip of the passenger k is obtained, and if the ith getting-off station value in Table1 is null, the get-off station of the ith bus trip of the passenger k cannot be determined.
7. Judgment calculation unit C i，i+1 Whether the (i + 1) th trip is the last trip of the day. If the judgment result is yes, the last trip end point of the day is the 1 st trip starting point of the day, and the all-day trip information of the passenger k is obtained. Proceed to the next step. If the judgment result is negative, i = i +1, the next calculation unit is switched to, and the step four is returned.
8. It is determined whether passenger k is the last passenger of the day. And if so, turning to the ninth step. If the judgment result is negative, k = k +1, returning to the step three, and calculating the getting-off station of the next passenger for all-day bus trip.
9. Through the steps, the getting-on stations of all passengers for all-day bus trips and the getting-off stations of part of passengers for bus trips are determined, and the value of the getting-off stations for the trips of the getting-off stations in Table1 cannot be determined to be null. Counting Trips with an getting-on station i and a getting-off station j in Table1 to obtain the bus trip OD quantity with i as a trip starting point and j as a trip end point in the bus trip OD Table of the passenger, and recording the bus trip OD quantity as the Trips ij Table1 is updated. Step ten is performed.
10. Calculating the travel OD of each passenger through the steps to obtain an initial OD matrix, wherein the travel OD quantity of the passenger with the get-off station value null in Table1 is not counted in the initial OD matrix, and the passenger with the get-off station value null is estimated according to the get-off station of the passenger getting on the same station as the passenger:
Trips pq *＝Trips pq ·b p
wherein,
b p -calculation coefficient of trip quantity with trip starting point p
-total bus trip amount with trip starting point p
-the travel starting point is p, and the travel end point is determined as the OD value of the bus travel
Trips pq * And the final bus travel OD amount of the whole day with the travel starting point being p and the travel end point being q.
Through the steps, the bus card holding passengers travel OD amount Trips all day long pq * Thus obtaining the product.
The second embodiment is as follows: the first difference between the present embodiment and the specific embodiment is: the second step, the specific method for determining the time for swiping the card when the passenger gets on the bus, the station for swiping the card when the passenger gets on the bus, the serial number of the bus line taken by the passenger and the driving direction is as follows:
1) Card number C of passenger k after ith card swiping k Time for passenger to get on bus and to swipe cardNumber of bus line taken by passengerThe card number C of the passenger k is read by the system from the IC card data stored in the bus trip k Time of passenger getting on bus and swiping cardNumber of bus line taken by passenger
2) System search number isThe vehicle-mounted GPS data of which the reading number isAt the time of the automobileRunning route ofAnd direction of travelDeterminingTime of dayThe stop information of the stop of the automobile is used as the getting-on stop of the bus trip
Other steps and parameters are the same as those in the first embodiment.
The third concrete implementation mode: the present embodiment differs from the first or second embodiment in that: the fifth step is as follows: judging whether the terminal point of the ith trip of the passenger k can be determined by utilizing the position relationship between the terminal point of the ith +1 card swiping of the passenger k and the ith trip line:
1): judging the (i + 1) th card swiping stationWhether to punch the card at the ith time for the automobile running lineIf the judgment result is yes, the next step is executed, and if the judgment result is no, the current calculation unit C cannot be calculated i，i+1 Determining the bus trip end point of the ith card swiping of the passenger according to the twice card swiping information;
2): judgment ofWhether the bus is in the ith card swiping bus driving directionIf the judgment result is yes, thenAnd obtaining the get-off station of the ith card swiping bus for the get-off station of the ith card swiping bus. If the judgment result is negative, the next step is carried out;
3): judging stationWhether or not a line exists within a range of 500 meters aroundEdge ofIf the judgment result of the station downstream in the direction is yes, the next step is executed, and if the judgment result is no, the current calculation unit C cannot be calculated i，i+1 The bus trip terminal point of the ith card swiping of the passenger is determined by the twice card swiping information;
4): judging whether the station retrieved in the step 3) is unique or not, if so, the station is the get-off station for the ith card swiping trip, and if not, selecting the distance stationAnd the nearest station is used as a get-off station in the current trip, and the get-off station for the ith card swiping trip of the passenger is obtained.
Other steps and parameters are the same as those in the first or second embodiment.
Examples
156 public transport lines and 2000 stations are shared in a city, and card swiping data of 10 passengers randomly selected from a passenger card swiping database at a certain day in the city are shown in table 3:
TABLE 3 data summary table for 10 passengers swiping card in a day
The starting point and the terminal point of each trip of the 10 passengers on the day are obtained by calculating through the first step to the sixth step in the method, wherein the starting point and the terminal point of each trip of the 10 passengers on the day are obtained, and the terminal point of each trip is null and is used for representing that the terminal point of each trip cannot be obtained by calculating through the steps. As shown in table 4:
table 4 calculation table for getting on or off bus station
The bus station of the passenger in the table 4 is converted into the bus travel OD matrix of the bus card holding passenger, the travel record of which the station value of the passenger is not null in the travel record table is directly converted into the bus travel OD matrix, and the initial OD matrix of the bus travel of the 10 passengers is obtained as shown in the table 5. In the embodiment, bus trip card swiping data of 10 passengers in one day is only considered, and the total trip amount of the 10 passengers in one day is less, so that the OD amount between bus stops is mostly 0. If the travel of the bus card holding passengers in the whole city is calculated, most of the travel OD quantity between the bus stops in the obtained OD matrix is not 0.
TABLE 5 initial OD matrix for bus trip
Note: the passenger travel OD with a null travel station is not included in the table.
And calculating the travel starting point and the travel ending point of 10 card holders and converting the travel starting point and the travel ending point into a travel OD matrix to obtain the initial OD matrix of the bus travel on the same day. For the travel record with the lower stop point value null, the travel record needs to be estimated through the bus travel initial OD matrix and the calculation formulas (1) and (2). The total day final bus travel OD quantity with the travel starting point i and the travel end point j in the initial OD matrix of the bus travel is the Trips pq And calculating to obtain the total daily bus travel OD (origin-destination) quantity Trips of the bus card holding passenger pq * As shown in table 6.
TABLE 6 bus trip OD matrix
Note: the table contains the calculated passenger travel OD with null station as the get-off station
Table 6 shows the bus trip OD of 10 passengers in the city in one day, and the calculation of the bus trip OD of 10 passengers in the city in one day is completed.
Claims (2)
1. A method for calculating the bus OD quantity by adopting a bus IC card and vehicle-mounted GPS data is characterized by comprising the following steps:
the method comprises the following steps: establishing a passenger historical travel information Table (Table 1), and storing and updating bus travel information of the bus card holder in the last 30 days;
wherein, bus trip information includes: the method comprises the following steps of (1) calculating card-holding card numbers of passengers, boarding card-swiping time, boarding card-swiping stations, bus line numbers of passengers, and calculated getting-off card-swiping time and getting-off card-swiping stations;
step two: reading IC card swiping data and vehicle-mounted GPS data of all card swiping passengers in one day, and determining information: the time of swiping the card when the passenger gets on the bus, the station of swiping the card when the passenger gets on the bus, the serial number of the bus line taken by the passenger and the driving direction;
establishing a card swiping information Table Table2 of the passengers for the same day to store the information, starting from the 1 st passenger in Table2, marking the kth passenger as a passenger k, and executing a step three, wherein k =1,2,3 and …;
step three: starting from the 1 st card swiping of the passenger k on the day, the card swiping corresponds to one bus trip each time, the starting point of the trip is the corresponding card swiping station, and the time difference delta T between the two adjacent card swiping times is calculated:
if the delta T is less than or equal to 30s, combining the two trips into one trip;
if Δ T>, 30s, i (i =1,2,3, …, m) th trip and i +1 th trip as a computing unit C i，i+1 Executing the step four; m is the total times of card swiping of the passenger k on the same day;
step four: judgment calculation unit C i，i+1 Whether the middle and twice trips are in the same trip route and different travel directions"two trips:
if the judgment result is yes, the ith trip and the (i + 1) th trip of the passenger k are round trip trips, the terminal point of the ith trip is the starting point of the (i + 1) th trip, and the terminal points of the ith trip and the (i + 1) th trip of the passenger k are recorded into Table1;
if the judgment result is negative, executing the step five;
step five: judging whether the terminal point of the ith trip of the passenger k can be determined by utilizing the position relationship between the terminal point of the ith +1 card swiping of the passenger k and the ith trip line:
if the judgment result is yes, recording the ith trip information of the passenger k into Table1; if the judgment result is negative, executing the step six;
judging whether the terminal point of the ith trip of the passenger k can be determined by utilizing the position relationship between the terminal point of the ith +1 card swiping of the passenger k and the ith trip line:
1): judging the (i + 1) th card swiping stationWhether to punch the card at the ith time for the automobile running lineIf the judgment result is yes, the next step is executed, and if the judgment result is no, the current calculation unit C cannot be calculated i，i+1 Determining the bus trip end point of the ith card swiping of the passenger according to the twice card swiping information;
2): judgment ofWhether the bus is in the ith card swiping bus driving directionIf the judgment result is yes, thenObtaining an getting-off station for the ith card swiping bus for the getting-off station for the ith card swiping bus;if the judgment result is negative, the next step is carried out;
3): judging stationWhether or not a line exists within a range of 500 meters aroundEdge ofIf the judgment result of the station downstream in the direction is yes, the next step is executed, and if the judgment result is no, the current calculation unit C cannot be calculated i，i+1 Determining the bus trip end point of the ith card swiping of the passenger according to the twice card swiping information;
4): judging whether the station retrieved in the step 3) is unique or not, if so, the station is the get-off station for the ith card swiping trip, and if not, selecting the distance stationThe nearest station is used as a get-off station for the current trip, and a get-off station for the ith card swiping trip of the passenger is obtained;
step six: reading passenger card number C of passenger k in Table2 k Ith driving routeDirection of travelTime T for swiping card on boarding i k And getting-on card swiping stationSearching Table1 for historical trip information satisfying the following conditions:
1) The boarding card swiping station is
2) The ith driving route is
3) The running direction is
4) The time for passenger to get on the bus and to swipe the card is T i k -30min，T i k +30min]；
If the trips meeting the conditions are searched, the starting point and the ending point of the searched trip and the ith trip of the passenger k are consistent, the ith trip ending point of the passenger k is obtained according to the searched trip ending point, and is stored in Table1, and the seventh step is executed;
if the travel meeting the conditions cannot be searched, marking the ith travel endpoint of the passenger k in Table1 as null; i = i +1, and returns to step four;
step seven: judgment calculation unit C i，i+1 To middle i Whether +1 trip is the last trip on the day:
if so, taking the last trip end point of the day as the 1 st trip starting point of the day, obtaining the all-day trip information of the passenger k, and executing the step eight;
if the judgment result is negative, i = i +1, returning to the fourth step;
step eight: determine whether passenger k is the last passenger of the day:
if the judgment result is yes, executing the step nine;
if the judgment result is negative, k = k +1, returning to the step three, and calculating the travel OD of the next passenger;
step nine: establishing a current-day passenger bus trip OD Table, converting an getting-on card swiping station and a getting-off station of each trip of a passenger in Table1 into the current-day passenger bus trip OD Table, updating Table1, and calculating each passenger trip OD through the steps to obtain an initial OD matrix, wherein p is a trip starting point and q is a bus trip O of a trip terminal pointD quantity is recorded as Trips pq Executing the step ten;
step ten: for the passenger travel OD without getting off the station, estimating according to the off station of the passenger getting on the same station
Trips pq * ＝Trips pq ·b p
Wherein,
b p -calculation coefficient of trip quantity with trip starting point p
-total bus trip amount with trip starting point p
-the travel starting point is p, and the travel end point is determined as the OD value of the bus travel
Trips pq * And the final bus travel OD amount of the whole day with the travel starting point being p and the travel end point being q.
2. The method for calculating the OD (origin-destination) quantity of the bus by adopting the bus IC card and the vehicle-mounted GPS (global positioning system) data as claimed in claim 1, wherein the specific method for determining the card swiping time of the passenger getting on the bus, the card swiping stop of the passenger getting on the bus, the serial number of the passenger taking the bus line and the driving direction in the second step is as follows:
1) Card number C of passenger k after ith card swiping k Time T for passenger to get on bus and to swipe card i k Number of bus line taken by passengerThe card number C of the passenger k is read by the system from the IC card data stored in the bus trip k Time T of passenger getting on bus and swiping card i k Number of bus line taken by passenger
2) System search number ofThe vehicle-mounted GPS data of which the reading number isAt time T i k Running route ofAnd direction of travelDetermining T i k Time of dayThe stop information of the stop of the automobile is used as the getting-on stop of the bus trip
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610040284.6A CN105654721B (en) | 2016-01-21 | 2016-01-21 | A kind of method for calculating public transport OD amounts using bus IC card and vehicle GPS data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610040284.6A CN105654721B (en) | 2016-01-21 | 2016-01-21 | A kind of method for calculating public transport OD amounts using bus IC card and vehicle GPS data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105654721A CN105654721A (en) | 2016-06-08 |
| CN105654721B true CN105654721B (en) | 2018-01-30 |
Family
ID=56486942
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610040284.6A Expired - Fee Related CN105654721B (en) | 2016-01-21 | 2016-01-21 | A kind of method for calculating public transport OD amounts using bus IC card and vehicle GPS data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105654721B (en) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107545730B (en) * | 2017-09-08 | 2020-04-07 | 哈尔滨工业大学 | Bus IC card data-based station getting-on/off passenger number estimation method |
| CN108288321A (en) * | 2018-01-24 | 2018-07-17 | 哈尔滨工业大学 | Method is determined based on the volume of the flow of passengers above and below the bus station of IC card data and vehicle GPS information |
| CN108269399B (en) * | 2018-01-24 | 2020-12-25 | 哈尔滨工业大学 | High-speed rail passenger flow demand prediction method based on road network passenger flow OD reverse thrust technology |
| CN110580537B (en) * | 2018-06-08 | 2023-06-16 | 阿里巴巴集团控股有限公司 | Bus bearing capacity analysis method and bus configuration optimization method |
| CN109102114B (en) * | 2018-08-02 | 2022-04-05 | 中南大学 | Bus trip getting-off station estimation method based on data fusion |
| CN110188923B (en) * | 2019-05-06 | 2022-01-07 | 西南交通大学 | Multi-mode bus passenger flow calculation method based on big data technology |
| CN114691708B (en) * | 2020-12-31 | 2023-05-23 | 南京行者易智能交通科技有限公司 | Visualization method and device based on passenger flow OD detail |
| CN116994435B (en) * | 2023-08-09 | 2024-05-03 | 广州地铁设计研究院股份有限公司 | Bus OD acquisition method and device |
| CN117275274B (en) * | 2023-11-20 | 2024-02-02 | 河北省交通规划设计研究院有限公司 | Conventional bus trip information identification method, device and medium |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101080042B1 (en) * | 2009-09-02 | 2011-11-04 | (주)경봉 | Traffic information system and method using simultaneous measurement |
| JP4693935B1 (en) * | 2010-06-14 | 2011-06-01 | 株式会社野村総合研究所 | TRAVEL TIME PROVIDING DEVICE, TRAVEL TIME PROVIDING METHOD, AND TRAVEL TIME PROVIDING PROGRAM |
| CN102324128B (en) * | 2011-05-24 | 2013-07-31 | 北京交通大学 | Method for predicting OD (Origin-Destination) passenger flow among bus stations on basis of IC (Integrated Circuit)-card record and device |
| CN103730008A (en) * | 2014-01-15 | 2014-04-16 | 汪涛 | Bus congestion degree analysis method based on real-time data of bus GPS (Global Position System) and IC (Integrated Circuit) cards |
| CN104809344B (en) * | 2015-04-23 | 2017-11-07 | 中山大学 | A kind of interval passenger flow method of estimation in the bus station based on IC-card data |
| CN105185105B (en) * | 2015-06-26 | 2017-05-03 | 东南大学 | Bus transfer identification method based on vehicle GPS and bus IC card data |
-
2016
- 2016-01-21 CN CN201610040284.6A patent/CN105654721B/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN105654721A (en) | 2016-06-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105654721B (en) | A kind of method for calculating public transport OD amounts using bus IC card and vehicle GPS data | |
| CN106504577B (en) | A kind of park and shift traffic path method and device for planning | |
| CN108831149B (en) | Method and system for customizing bus route running based on historical OD information | |
| CN107818412B (en) | Large-scale bus passenger OD parallel computing method based on MapReduce | |
| WO2018219337A1 (en) | Systems and methods for route planning | |
| CN105185105B (en) | Bus transfer identification method based on vehicle GPS and bus IC card data | |
| WO2017193556A1 (en) | Speed prediction method and apparatus | |
| CN102324128A (en) | Method for predicting OD (Origin-Destination) passenger flow among bus stations on basis of IC (Integrated Circuit)-card record and device | |
| CN108279017B (en) | Method for calculating and adding via points in real time in navigation process | |
| CN105788260A (en) | Public transportation passenger OD calculation method based on intelligent public transportation system data | |
| CN106767873A (en) | A kind of map-matching method based on space-time | |
| CN107203945A (en) | Vehicle insurance grading evaluation method and device | |
| CN110310477B (en) | Bus passenger flow detection method based on bus GPS and mobile phone signaling data | |
| CN103177412A (en) | Method and system for calculating taxi passenger information | |
| CN104916133B (en) | Road altitude information extraction method and system based on traffic track data | |
| CN110853156B (en) | Passenger OD identification method integrating bus GPS track and IC card data | |
| CN106897955A (en) | A kind of Public Transport Transfer recognition methods based on public transport OD data | |
| CN112036757A (en) | Parking transfer parking lot site selection method based on mobile phone signaling and floating car data | |
| CN106251628B (en) | A kind of method and device of the traffic trip amount of determining motor vehicle | |
| CN107578619B (en) | Method for measuring public bicycle service range of subway station based on IC card data | |
| CN114358808A (en) | Public transport OD estimation and distribution method based on multi-source data fusion | |
| CN115655301A (en) | Vehicle navigation route selection method and device, electronic equipment and medium | |
| CN114912689A (en) | Map grid index and XGBOST-based over-limit vehicle destination prediction method and system | |
| CN109801384A (en) | A kind of expressway charge calculation method and system based on digital map navigation | |
| CN112750219B (en) | Vehicle travel expense calculation method and system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180130 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |