CN108737951A - One kind monitoring stream of people's movement track system based on Mobile Location Technology - Google Patents

One kind monitoring stream of people's movement track system based on Mobile Location Technology Download PDF

Info

Publication number
CN108737951A
CN108737951A CN201710270333.XA CN201710270333A CN108737951A CN 108737951 A CN108737951 A CN 108737951A CN 201710270333 A CN201710270333 A CN 201710270333A CN 108737951 A CN108737951 A CN 108737951A
Authority
CN
China
Prior art keywords
stay
candidate
point
stop
time
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.)
Pending
Application number
CN201710270333.XA
Other languages
Chinese (zh)
Inventor
王振宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wenzhou Lucheng District New Research Institute Of Advanced Technology
Original Assignee
Wenzhou Lucheng District New Research Institute Of Advanced Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wenzhou Lucheng District New Research Institute Of Advanced Technology filed Critical Wenzhou Lucheng District New Research Institute Of Advanced Technology
Priority to CN201710270333.XA priority Critical patent/CN108737951A/en
Publication of CN108737951A publication Critical patent/CN108737951A/en
Pending legal-status Critical Current

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W4/00Services specially adapted for wireless communication networks; Facilities therefor
    • H04W4/02Services making use of location information
    • H04W4/025Services making use of location information using location based information parameters
    • H04W4/027Services making use of location information using location based information parameters using movement velocity, acceleration information

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Navigation (AREA)
  • Traffic Control Systems (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

One kind monitoring stream of people's movement track system based on Mobile Location Technology, its tracing point based on stream of people's speed merges stroke identification, it is calculated including speed, candidate stop place generates, dwell point identification, the present invention solves the problems, such as the fix drift of mobile phone location data and shake, stroke accuracy of identification is high, recognition result, which can be traffic programme work, provides places number and the data in terms of the travel time, and has cost more lower than conventional traffic investigation method and shorter data update period.

Description

One kind monitoring stream of people's movement track system based on Mobile Location Technology
Technical field
The present invention relates to a kind of trip stroke identifying systems based on mobile phone positioning track data.
Background technology
Foreign countries achieve certain achievement in research in terms of the acquisition of GPS track data and processing at present, can be according to rail Mark data carry out stroke and identify to obtain relevant information of going on a journey, and with trip daily record comparative analysis result accuracy.
Domestic scholar has also carried out the research that GPS track data information automatically extracts aspect successively.Know specific to stroke Mainly there are heuristic approaches and clustering procedure two major classes in terms of other algorithm.Heuristic approaches consider that stop activity occurs for individual first When space-time performance or positioning device feature, using these as known experience, algorithm for design handles initial data.Tool Body algorithm includes the method based on interrecord gap (IGP), the method based on rest point and the method based on missing point, is based on direction character Method etc., clustering procedure mainly considers that track data is generally acquired in a manner of constant duration, therefore when stop activity occurs A large amount of tracing points will be had to be gathered near a certain position, therefore the method that can use cluster is identified.Specific algorithm includes It is clustered based on K- mean clusters, based on DBSCAN.Since GPS data and mobile phone location data are in positioning accuracy and the quality of data There is difference in aspect, the stroke identification method for GPS data is not necessarily suitable mobile phone location data.
Invention content
The purpose of the present invention is to provide one kind monitoring stream of people's movement track system based on Mobile Location Technology, to overcoming The a wide range of drift that mobile phone positioning method is generated in quiescing activities has good result, solves the positioning of mobile phone location data The problem of drift and shake, stroke accuracy of identification is high.
Solution above-mentioned technical problem, the present invention use following technical scheme:One kind monitoring people based on Mobile Location Technology Movement track system is flowed, the tracing point based on stream of people's speed merges stroke identification, including speed calculates, candidate stop place life At the identification of, dwell point,
Speed calculates:Velocity information is not included in the initial trace data of acquisition, according to the longitude of tracing point record, latitude Traveler is calculated in the speed of each tracing point with temporal information, and instantaneous velocity is using being averaged on one section of track where tracing point Speed replaces;For mobile phone positioning track data, use the air line distance between tracing point that trajectory path distance is replaced to participate in Speed calculates;
Candidate stop place generates:Merge including rest point, stop center calculation, dwell point identification;
Rest point merges:The division of rest point and transfer point relies primarily on the threshold speed of setting, and threshold speed generally takes Minimum speed lower limit in resident trip mode, the i.e. lower velocity limit of manner of walking, the walking speed of normal person generally 3~ Between 6km/h, the most jogging speed of manner of walking is about 0.8m/s, considers that speed calculates influence and the pre-stage test of error Basis takes 1m/s that tracing point is classified as threshold speed, and 2 or more continuous rest point ps are then merged into candidate Stop place Stay ', candidate stop place are that traveler is likely to occur the movable position of stop;
Stop center calculation:The coordinate of candidate stop place is calculated using time-weighted mode, first, calculating successively connects Continuous two candidate dwell point psi, psi+1Average coordinates (x(i,i+1),y(i,i+1)),
Ps in formula (1)i.x it is candidate dwell point psiLongitude coordinate;psi.y it is candidate dwell point psiLatitude coordinate;
Secondly, by the time interval Δ t between 2 candidate dwell points(i, i+1)With the stay time of entire candidate stop place Weight of the ratio of Stay ' Δs t as the average coordinates
Finally, the coordinate of candidate stop place is calculated,
In formula (2):Stay ' .x are the longitude coordinate of candidate stop place Stay ';Stay ' .y are candidate stop place The latitude coordinate of Stay ';
Dwell point identifies:By considering that two away from discrete time factors are accepted or rejected and merged to candidate stop place, Steps are as follows:
Step 1 reads first candidate stop place Stay '1, it puts it into and stops sequence Sq, it will
Stay′1Coordinate as stop sequence Sq centre coordinate;
Step 2 judges whether the candidate stop place that do not read also, if so, reading next candidate stop place Stay′i, calculate Stay 'iCentre coordinate and Sq centre coordinate distances D(i, Sq), go to step 3;If not, going to step 4;
Step 3D(i, Sq)Whether the distance threshold Td of setting is less than, if so, by Stay 'iIt is put into and stops sequence Sq, again The centre coordinate for stopping sequence Sq is calculated, step 2 is gone to;If not, going to step 4;
Step 4 calculates the time interval Sq. Δ t for stopping start time Sq.st and stopping finish time Sq.st of Sq, this The time interval Sq. Δs t at place is to stop in sequence to be carved into last at the beginning of the stay time of first candidate stop place Time interval between the finish time of the stay time of a candidate stop place;
Whether step 5Sq. Δs t is more than the time threshold Tt of setting, if so, the candidate stop place in Sq is merged into and stopped Stationary point:(1) by included in Sq track point deletion at the time of between Sq.st and moment Sq.et;(2) by moment Sq.st and moment The coordinate of all tracing points uses the centre coordinate of Sq to replace between Sq.et, goes to step 6;If not, the candidate in Sq stops Position is stayed not constitute dwell point, then tracing point all between moment Sq.st and moment Sq.et is determined as transfer point, judges Whether comprising the last one candidate stop place in Sq, if so, terminating the judgement of this dwell point, if not, going to step 6;
Step 6 empties the candidate stop place in Sq, by Stay 'iIt is put into Sq, by Stay 'iStop centre coordinate conduct The centre coordinate for stopping sequence Sq, goes to step 2, when comprising the last one candidate stop place, then also continuing in Sq Judgement to candidate stop place, to identify new dwell point, then in Stay 'i-1With Stay 'iBetween candidate transfer point It is judged as transfer point, when recalculating the centre coordinate for stopping sequence Sq in above-mentioned steps 2, using time-weighted method, Calculation formula such as following formula (3):
Sq.x=wighti·Stay′i.x+(1-wighti)·Sq.x
Sq.y=wighti·Stay′i.y+(1-wighti)·Sq.y (3)
In formula:wightiFor candidate stop place Stay 'iWeighted value;When Sq.st is that the stop of stop sequence Sq starts Between i.e. first candidate stop place stay time at the beginning of carve;Sq.et is to stop stop end time of sequence Sq to be The finish time of the stay time of the last one candidate stop place;Sq.x is to stop sequence Sq central point longitude coordinates;Sq.y To stop sequence Sq central point latitude coordinates.
The beneficial effects of the invention are as follows:The invention solves the problems, such as the fix drift of mobile phone location data and shake, row Journey accuracy of identification is high, and recognition result, which can be traffic programme work, provides places number and the data in terms of the travel time, and has Cost more lower than conventional traffic investigation method and shorter data update period.
Specific implementation mode
In order to enable the technological means that the present invention realizes, character of innovation are apparent to function, it is further illustrated below The present invention.
It is of the invention a kind of based on Mobile Location Technology monitoring stream of people's movement track system, for acquired mobile phone positioning rail Mark data characteristics uses for reference the heuristic approaches of GPS track data identification, proposes and realizes a kind of based on Mobile Location Technology prison Survey stream of people's movement track system.Tracing point is divided by speed, certain speed threshold value tracing point below is will be less than and closes And be candidate stop place, recycle distance threshold and time threshold to merge candidate stop place, so that it is determined that going out true Positive dwell point.This system includes that speed calculates, candidate stop place generates, dwell point identification.
Speed calculates:Velocity information is not included in the initial trace data of acquisition, according to the longitude of tracing point record, latitude Traveler is calculated in the speed of each tracing point with temporal information, and instantaneous velocity is using being averaged on one section of track where tracing point Speed replaces;For mobile phone positioning track data, use the air line distance between tracing point that trajectory path distance is replaced to participate in Speed calculates;
Candidate stop place generates:Merge including rest point, stop center calculation, dwell point identification;
Rest point merges:The division of rest point and transfer point relies primarily on the threshold speed of setting, and threshold speed generally takes Minimum speed lower limit in resident trip mode, the i.e. lower velocity limit of manner of walking, the walking speed of normal person generally 3~ Between 6km/h, the most jogging speed of manner of walking is about 0.8m/s, considers that speed calculates influence and the pre-stage test of error Basis takes 1m/s that tracing point is classified as threshold speed, and 2 or more continuous rest point ps are then merged into candidate Stop place Stay ', candidate stop place are that traveler is likely to occur the movable position of stop;
Stop center calculation:The coordinate of candidate stop place is calculated using time-weighted mode, first, calculating successively connects Continuous two candidate dwell point psi, psi+1Average coordinates (x(i,i+1),y(i,i+1)),
Ps in formula (1)i.x it is candidate dwell point psiLongitude coordinate;psi.y it is candidate dwell point psiLatitude coordinate;
Secondly, by the time interval Δ t between 2 candidate dwell points(i, i+1)With the stay time of entire candidate stop place Weight of the ratio of Stay ' Δs t as the average coordinates
Finally, the coordinate of candidate stop place is calculated,
In formula (2):Stay ' .x are the longitude coordinate of candidate stop place Stay ';Stay ' .y are candidate stop place The latitude coordinate of Stay ';
Dwell point identifies:By considering that two away from discrete time factors are accepted or rejected and merged to candidate stop place, Steps are as follows:
Step 1 reads first candidate stop place Stay '1, put it into and stop sequence Sq, by Stay '1Coordinate make To stop the centre coordinate of sequence Sq;
Step 2 judges whether the candidate stop place that do not read also, if so, reading next candidate stop place Stay′i, calculate Stay 'iCentre coordinate and Sq centre coordinate distances D(i, Sq), go to step 3;If not, going to step 4;
Step 3D(i, Sq)Whether the distance threshold Td of setting is less than, if so, by Stay 'iIt is put into and stops sequence Sq, again The centre coordinate for stopping sequence Sq is calculated, step 2 is gone to;If not, going to step 4;
Step 4 calculates the time interval Sq. Δ t for stopping start time Sq.st and stopping finish time Sq.st of Sq, this The time interval Sq. Δs t at place is to stop in sequence to be carved into last at the beginning of the stay time of first candidate stop place Time interval between the finish time of the stay time of a candidate stop place;
Whether step 5Sq. Δs t is more than the time threshold Tt of setting, if so, the candidate stop place in Sq is merged into and stopped Stationary point:(1) by included in Sq track point deletion at the time of between Sq.st and moment Sq.et;(2) by moment Sq.st and moment The coordinate of all tracing points uses the centre coordinate of Sq to replace between Sq.et, goes to step 6;If not, the candidate in Sq stops Position is stayed not constitute dwell point, then tracing point all between moment Sq.st and moment Sq.et is determined as transfer point, judges Whether comprising the last one candidate stop place in Sq, if so, terminating the judgement of this dwell point, if not, going to step 6;
Step 6 empties the candidate stop place in Sq, by Stay 'iIt is put into Sq, by Stay 'iStop centre coordinate conduct The centre coordinate for stopping sequence Sq, goes to step 2, when comprising the last one candidate stop place, then also continuing in Sq Judgement to candidate stop place, to identify new dwell point, then in Stay 'i-1With Stay 'iBetween candidate transfer point It is judged as transfer point, when recalculating the centre coordinate for stopping sequence Sq in above-mentioned steps 2, using time-weighted method, Calculation formula such as following formula (3):
Sq.x=wighti·Stay′i.x+(1-wighti)·Sq.x
Sq.y=wighti·Stay′i.y+(1-wighti)·Sq.y (3)
In formula:wightiFor candidate stop place Stay 'iWeighted value;When Sq.st is that the stop of stop sequence Sq starts Between i.e. first candidate stop place stay time at the beginning of carve;Sq.et is to stop stop end time of sequence Sq to be The finish time of the stay time of the last one candidate stop place;Sq.x is to stop sequence Sq central point longitude coordinates;Sq.y To stop sequence Sq central point latitude coordinates.
The problem of system solves the problem the fix drift of mobile phone location data and shakes, stroke accuracy of identification is high, identification knot Fruit, which can be traffic programme work, provides places number and the data in terms of the travel time, and with more than conventional traffic investigation method Low cost and shorter data update period.The present invention is a wide range of to overcoming mobile phone positioning method to be generated in quiescing activities Drift has good result, solves the problems, such as the fix drift of mobile phone location data and shake, and stroke accuracy of identification is high, identification As a result be further mining analysis traffic analysis region in inhabitant go on a journey chain information, including trip number, trip distance, go out Row time etc. has established reliable basis.

Claims (1)

1. one kind monitoring stream of people's movement track system based on Mobile Location Technology, it is characterised in that:Its rail based on stream of people's speed Mark point merges stroke identification, including speed calculates, candidate stop place generates, dwell point identification,
Speed calculates:In the initial trace data of acquisition do not include velocity information, according to tracing point record longitude, latitude and when Between information calculate traveler each tracing point speed, instantaneous velocity use tracing point where one section of track on average speed To replace;For mobile phone positioning track data, use the air line distance between tracing point that trajectory path distance is replaced to participate in speed It calculates;
Candidate stop place generates:Merge including rest point, stop center calculation, dwell point identification;
Rest point merges:The division of rest point and transfer point relies primarily on the threshold speed of setting, and threshold speed generally takes resident Minimum speed lower limit in trip mode, the i.e. lower velocity limit of manner of walking, the walking speed of normal person is generally in 3~6km/h Between, the most jogging speed of manner of walking is about 0.8m/s, considers that speed calculates the influence of error and the basis of pre-stage test, It takes 1m/s that tracing point is classified as threshold speed, 2 or more continuous rest point ps is then merged into candidate stop Position Stay ', candidate stop place are that traveler is likely to occur the movable position of stop;
Stop center calculation:The coordinate of candidate stop place is calculated using time-weighted mode, first, calculates continuous two successively Candidate dwell point psi, psi+1Average coordinates (x(i,i+1),y(i,i+1)),
Ps in formula (1)i.x it is candidate dwell point psiLongitude coordinate;psi.y it is candidate dwell point psiLatitude coordinate;
Secondly, by the time interval Δ t between 2 candidate dwell points(i, i+1)With the stay time Stay ' Δs of entire candidate stop place Weight of the ratio of t as the average coordinates
Finally, the coordinate of candidate stop place is calculated,
In formula (2):Stay ' .x are the longitude coordinate of candidate stop place Stay ';Stay ' .y are candidate stop place Stay's ' Latitude coordinate;
Dwell point identifies:By considering that two away from discrete time factors are accepted or rejected and merged to candidate stop place, step It is as follows:
Step 1 reads first candidate stop place Stay '1, put it into and stop sequence Sq, by Stay '1Coordinate as stopping Stay the centre coordinate of sequence Sq;
Step 2 judges whether the candidate stop place that do not read also, if so, reading next candidate stop place Stay ′i, calculate Stay 'iCentre coordinate and Sq centre coordinate distances D(i, Sq), go to step 3;If not, going to step 4;
Step 3 D(i, Sq)Whether the distance threshold Td of setting is less than, if so, by Stay 'iIt is put into and stops sequence Sq, recalculate The centre coordinate for stopping sequence Sq, goes to step 2;If not, going to step 4;
Step 4 calculates the time interval Sq. Δ t for stopping start time Sq.st and stopping finish time Sq.st of Sq, herein Time interval Sq. Δs t is to stop in sequence to be carved into the last one time at the beginning of the stay time of first candidate stop place Select the time interval between the finish time of the stay time of stop place;
Whether step 5 Sq. Δs t is more than the time threshold Tt of setting, if so, stop is merged into the candidate stop place in Sq Point:(1) by included in Sq track point deletion at the time of between Sq.st and moment Sq.et;(2) by moment Sq.st and moment The coordinate of all tracing points uses the centre coordinate of Sq to replace between Sq.et, goes to step 6;If not, the candidate in Sq stops Position is stayed not constitute dwell point, then tracing point all between moment Sq.st and moment Sq.et is determined as transfer point, judges Whether comprising the last one candidate stop place in Sq, if so, terminating the judgement of this dwell point, if not, going to step 6;
Step 6 empties the candidate stop place in Sq, by Stay 'iIt is put into Sq, by Stay 'iStop centre coordinate as stop The centre coordinate of sequence Sq, goes to step 2, when not comprising the last one candidate stop place, then also continuing to waiting in Sq The judgement for selecting stop place, to identify new dwell point, then in Stay 'i-1With Stay 'iBetween candidate transfer point be also judged to It is set to transfer point, when recalculating the centre coordinate for stopping sequence Sq in above-mentioned steps 2, using time-weighted method, calculates Formula such as following formula (3):
Sq.x=wighti·Stay′i.x+(1-wighti)·Sq.x
Sq.y=wighti·Stay′i·y+(1-wighti)·Sq·y (3)
In formula:wightiFor candidate stop place Stay 'iWeighted value;Sq.st is to stop stop time started of sequence Sq to be It is carved at the beginning of the stay time of first candidate stop place;The stop end time that Sq.et is stop sequence Sq, that is, last The finish time of the stay time of one candidate stop place;Sq.x is to stop sequence Sq central point longitude coordinates;Sq.y is to stop Stay sequence Sq central point latitude coordinates.
CN201710270333.XA 2017-04-24 2017-04-24 One kind monitoring stream of people's movement track system based on Mobile Location Technology Pending CN108737951A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710270333.XA CN108737951A (en) 2017-04-24 2017-04-24 One kind monitoring stream of people's movement track system based on Mobile Location Technology

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710270333.XA CN108737951A (en) 2017-04-24 2017-04-24 One kind monitoring stream of people's movement track system based on Mobile Location Technology

Publications (1)

Publication Number Publication Date
CN108737951A true CN108737951A (en) 2018-11-02

Family

ID=63933799

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710270333.XA Pending CN108737951A (en) 2017-04-24 2017-04-24 One kind monitoring stream of people's movement track system based on Mobile Location Technology

Country Status (1)

Country Link
CN (1) CN108737951A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110769572A (en) * 2019-11-14 2020-02-07 安徽节源环保科技有限公司 Light control system and method based on GIS and mobile phone positioning
CN111508228A (en) * 2020-04-01 2020-08-07 佛山市城市规划设计研究院 Method for acquiring public transport trip chain by using mobile phone GPS and electronic map data
CN111812689A (en) * 2020-07-23 2020-10-23 中国平安财产保险股份有限公司 User behavior analysis method and device based on GPS track, electronic equipment and medium

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110769572A (en) * 2019-11-14 2020-02-07 安徽节源环保科技有限公司 Light control system and method based on GIS and mobile phone positioning
CN110769572B (en) * 2019-11-14 2021-08-13 安徽节源环保科技有限公司 Light control system and method based on GIS and mobile phone positioning
CN111508228A (en) * 2020-04-01 2020-08-07 佛山市城市规划设计研究院 Method for acquiring public transport trip chain by using mobile phone GPS and electronic map data
CN111508228B (en) * 2020-04-01 2021-01-01 佛山市城市规划设计研究院 Method for acquiring public transport trip chain by using mobile phone GPS and electronic map data
CN111812689A (en) * 2020-07-23 2020-10-23 中国平安财产保险股份有限公司 User behavior analysis method and device based on GPS track, electronic equipment and medium

Similar Documents

Publication Publication Date Title
CN110047277B (en) Urban road traffic jam ranking method and system based on signaling data
US9900747B1 (en) Using telematics data to identify a type of a trip
US9763055B2 (en) Travel and activity capturing
CN109636829A (en) A kind of multi-object tracking method based on semantic information and scene information
CN106096662B (en) Human motion state identification based on acceleration transducer
CN110505583A (en) A kind of path matching algorithm based on bayonet data and signaling data
CN108413968A (en) A kind of method and system of movement identification
CN108737951A (en) One kind monitoring stream of people's movement track system based on Mobile Location Technology
US20150285639A1 (en) Method and system for crowd sensing to be used for automatic semantic identification
CN103477241A (en) Device position estimates from motion and ambient light classifiers
CN105279483A (en) Fall-down behavior real-time detection method based on depth image
CN112543427A (en) Method and system for analyzing and identifying urban traffic corridor based on signaling track and big data
CN106448173B (en) A kind of long range trip traffic modal splitting method based on data in mobile phone
CN110287825A (en) It is a kind of that motion detection method is fallen down based on crucial skeleton point trajectory analysis
CN112862156A (en) Ship path planning method based on ship track and ant colony algorithm
Zhu et al. Learning transportation annotated mobility profiles from GPS data for context-aware mobile services
CN108510011B (en) User travel mode analysis method based on mobile phone multi-sensor
CN109034187A (en) A kind of subscriber household work address excavation process
CN108538054A (en) A kind of method and system obtaining traffic information based on mobile phone signaling data
CN110310477A (en) Bus passenger flow detection method based on public transport GPS Yu mobile phone signaling data
CN112163568B (en) Scenic spot person searching system based on video detection
CN109297489A (en) A kind of indoor navigation method based on user characteristics, electronic equipment and storage medium
CN109002189A (en) A kind of motion recognition method, device, equipment and computer storage medium
CN113079463A (en) Tourist attraction tourist travel activity identification method based on mobile phone signaling data
CN115239839A (en) Indoor swimming track data processing method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20181102