Method and device for judging intercity transportation mode based on mobile phone data
A method and device for judging intercity transportation modes based on mobile phone data solves the problems of difficult data acquisition and low data processing efficiency in the existing solutions. The method includes: obtaining the trajectory of each user from the mobile phone data; extracting an intercity traffic sub-trajectory in the trajectory of each user; and judging the transportation mode of the intercity traffic sub-trajectory. The solution has the advantages of easy data acquisition, high data processing efficiency and capability of achieving long-term monitoring.
1. A method for judging intercity transportation modes based on mobile phone data, the method comprising:
obtaining a trajectory of each user from the mobile phone data such that m trajectories are obtained and being represented as:
S={s 1 ,S 2 . . . s m },S i ∈S;
extracting an intercity traffic sub-trajectory in the trajectory of each user including:
judging whether a city C 3 and a city C 4 corresponding to the sub-trajectory s i are adjacent,
if the city C 3 and the city C 4 are adjacent, adding the sub-trajectory s i to a sub-trajectory set R, and
if the city C 3 and the city C 4 are not adjacent, judging whether the sub-trajectory s i is a sub-trajectory with a transit city,
if the sub-trajectory has the transit city, correspondingly dividing the sub-trajectory s i into two segments of sub-trajectories according to the transit city, and adding the two segments of sub-trajectories to the sub-trajectory set R,
if the sub-trajectory has no transit city, adding the sub-trajectory s i to the sub-trajectory set R, and
processing all sub-trajectories in the sub-trajectory set R based on a method of processing the sub-trajectory s i until all sub-trajectories in the sub-trajectory set R do not include a transit city; and
judging the transportation mode of the intercity traffic sub-trajectory as being a flight travel, a train travel, or a road transportation, wherein:
the sub-trajectory r i is judged as being the flight travel based on:
if there is a flight between a city C 5 and a city C 6 corresponding to the sub-trajectory r i , if the time of the sub-trajectory r i passing A 1 or A 2 exceeds an airport stay threshold, and at least the A 1 or the A 2 has no railway station or bus station, the transportation mode corresponding to the sub-trajectory r i is the flight travel, or
if the maximum speed and the average speed of the sub-trajectory r i is within a preset aircraft navigation speed range, the transportation mode corresponding to the sub-trajectory r i is the flight travel;
the sub-trajectory r i is judged as being the train travel based on:
if there is a train between the city C 5 and the city C 6 corresponding to the sub-trajectory r i , the time of the sub-trajectory r i passing T 1 and the time of the sub-trajectory r i passing T 2 both exceed a railway station stay threshold, the maximum speed or the average speed of the sub-trajectory r i exceeds a preset motor vehicle speed range, and the time of the sub-trajectory r i passing the T 1 or T 2 meets a train number in a railway timetable, the transportation mode corresponding to the sub-trajectory r i is the train travel,
if one of the time of the sub-trajectory r i passing the T 1 and the time of the sub-trajectory r i passing the T 2 exceeds the railway station stay threshold, and the maximum speed or the average speed of the sub-trajectory r i is within a preset first train speed range, wherein the lower boundary of the first train speed range exceeds the upper boundary of the motor vehicle speed range, and the difference value between the lower boundary of the first train speed range and the upper boundary of the motor vehicle speed range is greater than a speed difference preset value, the transportation mode corresponding to the sub-trajectory r i is the train travel,
if one of the time of the sub-trajectory r i passing the T 1 and the time of the sub-trajectory r i passing the T 2 exceeds the railway station stay threshold, the maximum speed or the average speed of the sub-trajectory r i exceeds the preset motor vehicle speed range, and the difference value between the movement direction of the sub-trajectory r i within an outbound time after passing the T 1 and the movement direction of the sub-trajectory r i within an inbound time before entering the T 2 is smaller than a direction deviation preset value, the transportation mode corresponding to the sub-trajectory r i is the train travel, or
if neither of the time of the sub-trajectory r i passing the T 1 and the time of the sub-trajectory r i passing the T 2 exceeds the railway station stay threshold, and the maximum speed or the average speed of the sub-trajectory r i is within a preset second train speed range, then the transportation mode corresponding to the sub-trajectory r i is the train travel; and
if the sub-trajectory r i fails to be determined as the flight travel or the train travel, the sub-trajectory r i is judged as being the road transportation, the sub-trajectory r i corresponds to the trajectory between the city C 5 and the city C 6 ; and
A 1 represents an airport area corresponding to the city C 5 , A 2 represents the airport area corresponding to the city C 6 , T 1 represents a railway station area corresponding to the city C 5 , T 2 represents the railway station area corresponding to the city C 6 , B 1 represents a bus station area corresponding to the city C 5 , and B 2 represents the bus station area corresponding to the city C 6 , the sub-trajectory r i ∈R.
2. The method according to claim 1 , wherein
after the transportation mode corresponding to the sub-trajectory r i is the train travel, the method further comprises the following steps:
if the maximum speed or the average speed of the sub-trajectory r i is within a preset high-speed rail range, then the transportation mode corresponding to the sub-trajectory r i is high-speed rail travel, or otherwise, the transportation mode corresponding to the sub-trajectory r i is ordinary train travel.
3. The method according to claim 1 , wherein
after the transportation mode corresponding to the sub-trajectory r i is the road transportation, the method further comprises the following steps:
if the sub-trajectory r i passes a highway toll station, or a highway entrance, exit or service area, and the average speed of the sub-trajectory r i is greater than or equal to 100 km/h, then the transportation mode corresponding to the sub-trajectory r i is highway travel; and if the sub-trajectory r i does not pass a highway toll station, or a highway entrance, exit or service area, and the average speed of the sub-trajectory r i is less than 100 km/h, then the transportation mode corresponding to the sub-trajectory r i is ordinary road travel.
4. The method according to claim 1 , wherein
after the step of judging whether the sub-trajectory r i is the road transportation, the method further comprises the following steps:
if the transportation mode corresponding to the sub-trajectory r i fails to be determined, obtaining a corresponding pre-judgment base station sequence set according to the sub-trajectory r i , separately obtaining corresponding typical base station sequence sets according to a typical aviation trajectory, a typical high-speed rail trajectory, a typical ordinary train trajectory and a typical highway trajectory, calculating the similarities between the pre-judgment base station sequence set and the typical base station sequence sets respectively, if one similarity S 1 is greater than the other similarities, the difference value between the similarity S 1 and the closest similarity S 2 is greater than a preset similarity difference value, and the similarity is greater than a similarity threshold, then the transportation mode of the sub-trajectory r i is the same as the transportation mode corresponding to the base station sequence set for calculating the similarity.
5. The method according to claim 1 , wherein
after the step of judging whether the sub-trajectory r i is the road transportation, the method further comprises the following steps:
if the transportation mode corresponding to the sub-trajectory r i fails to be determined, the transportation mode corresponding to the sub-trajectory r i is a transportation mode other than the flight travel, the train travel and the road transportation.
6. A non-transitory computer readable storage medium storing a computer program that causes a processor to perform steps comprising:
obtaining a trajectory of each user from the mobile phone data such that m trajectories are obtained and being represented as:
S={s 1 ,s 2 . . . s m },s i ∈S;
extracting an intercity traffic sub-trajectory in the trajectory of each user including:
judging whether a city C 3 and a city C 4 corresponding to the sub-trajectory s i are adjacent,
if the city C 3 and the city C 4 are adjacent, adding the sub-trajectory s i to a sub-trajectory set R, and
if the city C 3 and the city C 4 are not adjacent, judging whether the sub-trajectory s i is a sub-trajectory with a transit city,
if the sub-trajectory has the transit city, correspondingly dividing the sub-trajectory s i into two segments of sub-trajectories according to the transit city, and adding the two segments of sub-trajectories to the sub-trajectory set R,
if the sub-trajectory has no transit city, adding the sub-trajectory s i to the sub-trajectory set R, and
processing all sub-trajectories in the sub-trajectory set R based on a method of processing the sub-trajectory s i until all sub-trajectories in the sub-trajectory set R do not include a transit city; and
judging the transportation mode of the intercity traffic sub-trajectory as being a flight travel, a train travel, or a road transportation, wherein:
the sub-trajectory r i is judged as being the flight travel based on:
if there is a flight between a city C 5 and a city C 6 corresponding to the sub-trajectory r i , if the time of the sub-trajectory r i passing A 1 or A 2 exceeds an airport stay threshold, and at least the A 1 or the A 2 has no railway station or bus station, the transportation mode corresponding to the sub-trajectory r i is the flight travel, or
if the maximum speed and the average speed of the sub-trajectory r i is within a preset aircraft navigation speed range, the transportation mode corresponding to the sub-trajectory r i is the flight travel;
the sub-trajectory r i is judged as being the train travel based on:
if there is a train between the city C 5 and the city C 6 corresponding to the sub-trajectory r i , the time of the sub-trajectory r i passing T 1 and the time of the sub-trajectory r i passing T 2 both exceed a railway station stay threshold, the maximum speed or the average speed of the sub-trajectory r i exceeds a preset motor vehicle speed range, and the time of the sub-trajectory r i passing the T 1 or T 2 meets a train number in a railway timetable, the transportation mode corresponding to the sub-trajectory r i is the train travel,
if one of the time of the sub-trajectory r i passing the T 1 and the time of the sub-trajectory r i passing the T 2 exceeds the railway station stay threshold, and the maximum speed or the average speed of the sub-trajectory r i is within a preset first train speed range, wherein the lower boundary of the first train speed range exceeds the upper boundary of the motor vehicle speed range, and the difference value between the lower boundary of the first train speed range and the upper boundary of the motor vehicle speed range is greater than a speed difference preset value, the transportation mode corresponding to the sub-trajectory r i is the train travel,
if one of the time of the sub-trajectory r i passing the T 1 and the time of the sub-trajectory r i passing the T 2 exceeds the railway station stay threshold, the maximum speed or the average speed of the sub-trajectory r i exceeds the preset motor vehicle speed range, and the difference value between the movement direction of the sub-trajectory r i within an outbound time after passing the T 1 and the movement direction of the sub-trajectory r i within an inbound time before entering the T 2 is smaller than a direction deviation preset value, the transportation mode corresponding to the sub-trajectory r i is the train travel, or
if neither of the time of the sub-trajectory r i passing the T 1 and the time of the sub-trajectory r i passing the T 2 exceeds the railway station stay threshold, and the maximum speed or the average speed of the sub-trajectory r i is within a preset second train speed range, then the transportation mode corresponding to the sub-trajectory r i is the train travel; and
if the sub-trajectory r i fails to be determined as the flight travel or the train travel, the sub-trajectory r i is judged as being the road transportation, the sub-trajectory r i corresponds to the trajectory between the city C 5 and the city C 6 ; and
A 1 represents an airport area corresponding to the city C 5 , A 2 represents the airport area corresponding to the city C 6 , T 1 represents a railway station area corresponding to the city C 5 , T 2 represents the railway station area corresponding to the city C 6 , B 1 represents a bus station area corresponding to the city C 5 , and B 2 represents the bus station area corresponding to the city C 6 , the sub-trajectory r i ∈R.
7. A computer device, comprising:
a processor programmed to:
obtaining a trajectory of each user from the mobile phone data such that m trajectories are obtained and being represented as:
S={s 1 ,s 2 . . . s m },s i ∈S;
extracting an intercity traffic sub-trajectory in the trajectory of each user including:
judging whether a city C 3 and a city C 4 corresponding to the sub-trajectory s i are adjacent,
if the city C 3 and the city C 4 are adjacent, adding the sub-trajectory s i to a sub-trajectory set R, and
if the city C 3 and the city C 4 are not adjacent, judging whether the sub-trajectory s i is a sub-trajectory with a transit city,
if the sub-trajectory has the transit city, correspondingly dividing the sub-trajectory s i into two segments of sub-trajectories according to the transit city, and adding the two segments of sub-trajectories to the sub-trajectory set R,
if the sub-trajectory has no transit city, adding the sub-trajectory s i to the sub-trajectory set R, and
processing all sub-trajectories in the sub-trajectory set R based on a method of processing the sub-trajectory s i until all sub-trajectories in the sub-trajectory set R do not include a transit city; and
judging the transportation mode of the intercity traffic sub-trajectory as being a flight travel, a train travel, or a road transportation, wherein:
the sub-trajectory r i is judged as being the flight travel based on:
if there is a flight between a city C 5 and a city C 6 corresponding to the sub-trajectory r i , if the time of the sub-trajectory r i passing A 1 or A 2 exceeds an airport stay threshold, and at least the A 1 or the A 2 has no railway station or bus station, the transportation mode corresponding to the sub-trajectory r i is the flight travel, or
if the maximum speed and the average speed of the sub-trajectory r i is within a preset aircraft navigation speed range, the transportation mode corresponding to the sub-trajectory r i is the flight travel;
the sub-trajectory r i is judged as being the train travel based on:
if there is a train between the city C 5 and the city C 6 corresponding to the sub-trajectory r i , the time of the sub-trajectory r i passing T 1 and the time of the sub-trajectory r i passing T 2 both exceed a railway station stay threshold, the maximum speed or the average speed of the sub-trajectory r i exceeds a preset motor vehicle speed range, and the time of the sub-trajectory r i passing the T 1 or T 2 meets a train number in a railway timetable, the transportation mode corresponding to the sub-trajectory r i is the train travel,
if one of the time of the sub-trajectory r i passing the T 1 and the time of the sub-trajectory r i passing the T 2 exceeds the railway station stay threshold, and the maximum speed or the average speed of the sub-trajectory r i is within a preset first train speed range, wherein the lower boundary of the first train speed range exceeds the upper boundary of the motor vehicle speed range, and the difference value between the lower boundary of the first train speed range and the upper boundary of the motor vehicle speed range is greater than a speed difference preset value, the transportation mode corresponding to the sub-trajectory r i is the train travel,
if one of the time of the sub-trajectory r i passing the T 1 and the time of the sub-trajectory r i passing the T 2 exceeds the railway station stay threshold, the maximum speed or the average speed of the sub-trajectory r i exceeds the preset motor vehicle speed range, and the difference value between the movement direction of the sub-trajectory r i within an outbound time after passing the T 1 and the movement direction of the sub-trajectory r i within an inbound time before entering the T 2 is smaller than a direction deviation preset value, the transportation mode corresponding to the sub-trajectory r i is the train travel, or
if neither of the time of the sub-trajectory r i passing the T 1 and the time of the sub-trajectory r i passing the T 2 exceeds the railway station stay threshold, and the maximum speed or the average speed of the sub-trajectory r i is within a preset second train speed range, then the transportation mode corresponding to the sub-trajectory r i is the train travel; and
if the sub-trajectory r i fails to be determined as the flight travel or the train travel, the sub-trajectory r i is judged as being the road transportation, the sub-trajectory r i corresponds to the trajectory between the city C 5 and the city C 6 ; and
A 1 represents an airport area corresponding to the city C 5 , A 2 represents the airport area corresponding to the city C 6 , T 1 represents a railway station area corresponding to the city C 5 , T 2 represents the railway station area corresponding to the city C 6 , B 1 represents a bus station area corresponding to the city C 5 , and B 2 represents the bus station area corresponding to the city C 6 , the sub-trajectory r i ∈R.