CN112462401B - Urban canyon rapid detection method and device based on floating vehicle track data - Google Patents
Urban canyon rapid detection method and device based on floating vehicle track data Download PDFInfo
- Publication number
- CN112462401B CN112462401B CN202110147584.5A CN202110147584A CN112462401B CN 112462401 B CN112462401 B CN 112462401B CN 202110147584 A CN202110147584 A CN 202110147584A CN 112462401 B CN112462401 B CN 112462401B
- Authority
- CN
- China
- Prior art keywords
- track
- points
- grid
- data
- point
- 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.)
- Active
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 20
- 238000007667 floating Methods 0.000 title claims abstract description 18
- 238000005070 sampling Methods 0.000 claims abstract description 43
- 230000002159 abnormal effect Effects 0.000 claims abstract description 32
- 238000004458 analytical method Methods 0.000 claims abstract description 18
- 239000004576 sand Substances 0.000 claims description 12
- BULVZWIRKLYCBC-UHFFFAOYSA-N phorate Chemical compound CCOP(=S)(OCC)SCSCC BULVZWIRKLYCBC-UHFFFAOYSA-N 0.000 claims description 10
- 230000004931 aggregating effect Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000013075 data extraction Methods 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 18
- 230000008859 change Effects 0.000 abstract description 7
- 238000010276 construction Methods 0.000 abstract description 3
- 238000004220 aggregation Methods 0.000 abstract 1
- 230000008569 process Effects 0.000 description 7
- 230000000694 effects Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000004807 localization Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000005206 flow analysis Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 238000012732 spatial analysis Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/393—Trajectory determination or predictive tracking, e.g. Kalman filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/50—Determining position whereby the position solution is constrained to lie upon a particular curve or surface, e.g. for locomotives on railway tracks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
-
- 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
- G08G1/0112—Measuring and analyzing of parameters relative to traffic conditions based on the source of data from the vehicle, e.g. floating car data [FCD]
-
- 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
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Theoretical Computer Science (AREA)
- Databases & Information Systems (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Traffic Control Systems (AREA)
- Navigation (AREA)
Abstract
The invention relates to the technical field of satellite navigation and positioning, and provides an urban canyon rapid detection method and device based on floating car track data, which comprises the following steps: selecting a sampling area to obtain a track data setP(ii) a Establishing a buffer area for the urban road in the sampling area, and generating a grid set corresponding to the urban roadR(ii) a Extracting a set of trajectory dataPTrack data of a certain vehicleP K (ii) a Extracting trajectory dataP K The abnormal points in (1) are set into the gridsRPerforming superposition analysis with the abnormal points and marking the grid setRA grid intersecting the outlier; from a set of trajectory dataPRemoving track dataP K (ii) a Determining a set of trajectory dataPWhether it is empty; statistical grid aggregationRThe number of times each grid is marked exceeds a thresholdHThe grid output of (a) is an urban canyon region. The track data used by the method is high in updating speed, the timeliness of the detection result is better, and the urban canyon region change in urban construction change can be reflected more truly.
Description
Technical Field
The invention relates to the technical field of satellite navigation and positioning, in particular to a method and a device for rapidly detecting urban canyons based on floating car track data.
Background
A Global Navigation Satellite System (GNSS) is a space-based wireless Navigation System that can provide information on three-dimensional coordinates, motion trajectories, etc. to people at any position on the earth's surface or in the near-earth space. The GNSS is widely applied to various fields of national defense and economic society, such as military affairs, communication, electric power, finance, transportation, basic mapping, disaster prevention and reduction and the like. With the continuous development of GNSS, the positioning accuracy is also improved. The performance of satellite positioning is closely related to the scene, including the number of visible satellites, geometric accuracy factors, signal interruption frequency, signal strength, multipath effect, deception source and other factors, which all affect the accuracy and usability of the satellite positioning result. In an open environment, the quality of a GNSS signal is good, in an urban environment, the GNSS signal is often shielded by high buildings, terrains, dense forests or reflected for multiple times to cause the situation that the positioning precision is poor or even the positioning cannot be performed, and the urban environment with the problems is called as an urban canyon. When the positioning terminal is located in an urban canyon, the satellite signal positioning performance is remarkably reduced, so that the related application is difficult to meet the actual production and living requirements. If the range of the urban canyon can be detected in advance, corresponding algorithm design and road traffic infrastructure transformation can be carried out, and the method is favorable for improving the positioning accuracy and the positioning effect in the urban canyon region.
In an urban complex environment, two main factors influencing the performance of the GNSS are occlusion and multipath effects, which are reflected by data loss or position drift at a certain positioning point on trajectory data. The traditional urban detection method mainly utilizes an urban three-dimensional model or a street view image, analyzes building height and position information, and combines a satellite altitude angle to calculate urban canyon distribution. However, the acquisition of the three-dimensional model data of the urban building has the defects of high cost, long time consumption and slow updating. The city is a constantly changing subject, and the position of the urban canyon is constantly changing along with the change of ground surface coverage data such as buildings, vegetation, water bodies and the like, so that the requirement of urban canyon change detection is difficult to meet by using urban three-dimensional data.
The above is only for the purpose of assisting understanding of the technical aspects of the present invention, and does not represent an admission that the above is prior art.
Disclosure of Invention
The invention mainly aims to solve the technical problem that the urban canyon change detection requirement is difficult to meet by using urban three-dimensional data in the prior art.
In order to achieve the purpose, the invention provides an urban canyon rapid detection method based on floating car track data, which comprises the following steps:
s1: selecting a sampling area, and carrying out track data acquisition on a plurality of vehicles in the sampling area to obtain a track data setP;
S2: establishing a buffer area for the urban road in the sampling area, and generating a grid set corresponding to the urban roadR;
S3: extracting the set of trajectory dataPTrack data of a certain vehicleP K ,KIs a positive integer and has an initial value of 1; proceeding to step S4;
s4: extracting the trajectory dataP K The outlier of (1), the outlier comprising: data missing points and space positioning drift points; aggregating the gridsRPerforming superposition analysis with the abnormal points and marking the grid setRThe grid intersecting the abnormal point is entered to step S5;
s5: from the set of trajectory dataPRemoving the trajectory dataP K To makeKAdding 1 to the value of (c); determining the trajectory data setPWhether it is empty; if yes, go to step S6; if not, returning to the step S3;
s6: counting the grid setRThe number of times each grid is marked is made to exceed a threshold valueHThe grid output of (a) is an urban canyon region.
Preferably, in step S2, the grid setREach grid in (a) is specifically represented as:
wherein,ja number representing a grid;CNT j representing the number of outliers in the grid;xl j represents the longitude of the lower left corner of the grid;yl j representing the latitude of the lower left corner of the grid;xR j represents the longitude of the upper right corner of the grid;yR j representing the latitude of the upper right corner of the grid; h denotes the number of the mesh set.
Preferably, in step S3, the trajectory dataP K The concrete expression is as follows:
wherein,Tthe maximum sampling times of the track;Knumbering the vehicles;P i is the first in the trackiTracing points;t i is shown asiSampling time of each trace point;v i indicating that the vehicle ist i The speed of travel at that moment;x i to representt i Longitude of the time track point;y i to representt i The latitude of the time trace point.
Preferably, in step S4, the step of extracting the data missing point in the outlier specifically includes:
s401: from the trajectory dataP K Middle pass sampling timet n Two adjacent track points are taken out in sequenceP n AndP n+1 ,nis from 1 toT-A positive integer of 1, or a mixture thereof,nis 1; if the track pointP n And track pointP n+1 The sampling time interval is greater than the preset timet m If the data missing point exists, the data missing point and the track point existP n And track pointP n+1 To (c) to (d);
s402: by means of tracing pointsP n And track pointP n+1 Obtaining the coordinates of the data missing point by the coordinate calculation of (2) so thatnThe value of (a) is added to 1, and the formula is specifically as follows:
wherein (A), (B), (C), (D), (C), (B), (C)x,y) Coordinates of the data missing point; (x n ,y n ) Is a track pointP n The coordinates of (a); (x n+1 ,y n+1 ) Tracing pointP n+1 The coordinates of (a);representing points of trackP n And track pointP n+1 The time interval between two adjacent points;t m is a preset time;
s403: repeating step S401 and step S402T-1 time, calculating said trajectory dataP K All of the data missing points in (1).
Preferably, in step S4, the step of extracting the spatial localization drift point in the outlier specifically includes:
s411: from the trajectory dataP K Middle pass sampling timet n Sequentially taking out three adjacent track pointsP n 、P n+1 AndP n+2 ,nis from 1 toT-A positive integer of 2 is a positive integer of,nis 1;
s412: judging the track pointP n 、P n+1 AndP n+2 the distance between and the average distanceSThe magnitude relationship of (1);
if it isP n AndP n+1 is less than or equal toSAnd is andP n+1 andP n+2 is greater thanSThen point of trackP n+2 For spatially locating drift points and track pointsP n+2 The coordinates of (a) are:
if it isP n AndP n+1 is greater thanSAnd is andP n+1 andP n+2 if the distance between the two points is less than or equal to S, the track points areP n For spatially locating drift points and track pointsP n The coordinates of (a) are:
if it isP n AndP n+1 is greater thanSAnd is andP n+1 andP n+2 if the distance between the two points is greater than S, the track points areP n+1 For spatially locating drift points and track pointsP n+1 The coordinates of (a) are:
wherein (A), (B), (C), (D), (C), (B), (C)x n ,y n ) Is a track pointP n The coordinates of (a); (x n+1 ,y n+1 ) Tracing pointP n+1 The coordinates of (a); (x n+2 ,y n+2 ) Is a track pointP n+2 The coordinates of (a);
s413: repeating the steps S411 and S412 for T-2 times, and calculating the track dataP K All of the spatially located drift points.
Preferably, step S6 is specifically:
when the grid is assembledRWhen a grid in (2) is marked once, the gridCNT j Adding 1 to the value of (c); when the grid is inCNT j Value greater than thresholdHThe grid is then output as an urban canyon region.
An urban canyon rapid detection device based on floating car track data comprises the following modules:
a track data acquisition module for selecting a sampling area, acquiring track data of a plurality of vehicles in the sampling area to obtain a track data setP;
A grid set generating module for establishing a buffer area for the urban road in the sampling area and generating a grid set corresponding to the urban roadR;
A track data extraction module for extracting the track data setPTrack data of a certain vehicleP K ,KIs a positive integer and has an initial value of 1; entering a superposition analysis module;
a superposition analysis module for extracting the trajectory dataP K The outlier of (1), the outlier comprising: data missing points and space positioning drift points; aggregating the gridsRPerforming superposition analysis with the abnormal points and marking the grid setREntering a grid intersected with the abnormal point into an updating track data set module;
an update track data set module for generating a track data set from the track data setPRemoving the trajectory dataP K To makeKAdding 1 to the value of (c); determining the trajectory data setPWhether it is empty; if yes, entering an urban canyon region judgment module; if not, returning to the track data extraction module;
the urban canyon region judgment module is used for counting the grid setRAnd outputting the grids with the marked times exceeding a threshold value as the urban canyon region.
The invention has the following beneficial effects:
1. the spatial distribution of the urban canyon region is detected by using the abnormal points existing in the large-scale track data, and the track data is simple to obtain and easy to realize;
2. track sampling time abnormity and space position abnormity are fully considered in the abnormal point extraction process, and the extraction result is more reliable and accurate;
3. according to the method, only the processing and analysis of the trajectory data need to be considered, the calculation of complex factors such as satellite height, angle, position, time period and the like related to urban canyon detection based on urban three-dimensional modeling data is avoided, and the calculation speed is higher;
4. the track data used by the method is high in updating speed, the timeliness of the detection result is better, and the urban canyon region change in urban construction change can be reflected more truly.
Drawings
FIG. 1 is a flow chart of a method for rapid urban canyon detection based on floating car trajectory data;
FIG. 2 is a diagram of an urban canyon region in Wuhan City;
FIG. 3 is a schematic diagram of an urban canyon rapid detection device based on floating car trajectory data;
the implementation, functional features and advantages of the objects of the present invention will be further explained with reference to the accompanying drawings.
Detailed Description
It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
With the increase of urban traffic track data in recent years, particularly the increase of GPS positioning devices installed on taxies and buses, floating car track data are acquired in a large scale, the data comprise geographic coordinates, speed, direction, time stamps and other information, the driving process of a vehicle can be completely recorded at a high frequency by virtue of a standard data structure, and a large amount of researches are carried out on the aspects of crowd flow analysis, trip hot spot area extraction, urban space analysis, urban traffic condition prediction and the like by adopting the urban floating car track data; the floating car track data inevitably has abnormal points, part of the abnormal points are generated due to random errors of the positioning terminal, more of the abnormal points are due to systematic influences existing in the driving process of the car, and most of the abnormal points are caused by the phenomena of signal shielding, reflection, refraction and the like in the driving process, and the phenomena are the concentrated reflection of the influence of urban canyon areas on GNSS signals; therefore, when the abnormal points in the trajectory data reach a certain scale, the distribution of the urban canyon region can be extracted completely by performing spatial analysis and processing on the abnormal points in the trajectory data.
Due to the fact that the track data are easy to obtain and fast to update, by exploring the relation between track positioning abnormal points and GNSS positioning abnormal areas, urban canyon areas around roads can be obtained fast in a large range, effective judgment means is provided for accurate updating of road infrastructures, construction and perfection of urban road infrastructures are facilitated, and effective data support is provided for intelligent transportation, unmanned driving, GNSS chip testing and the like;
referring to fig. 1, in view of the above analysis, the present invention provides a method for rapidly detecting an urban canyon based on floating car trajectory data, which obtains abnormal points that cannot be matched to a road network in the trajectory data by matching coordinates of trajectory points in the trajectory data with the road network, and performs spatial clustering on the abnormal points, thereby obtaining an urban canyon region along the road network, and specifically includes the following steps:
s1: selecting a sampling area, and carrying out track data acquisition on a plurality of vehicles in the sampling area to obtain a track data setP;
S2: establishing a buffer area for the urban road in the sampling area, and generating a grid set corresponding to the urban roadR;
S3: extracting the set of trajectory dataPTrack data of a certain vehicleP K ,KIs a positive integer and has an initial value of 1; proceeding to step S4;
s4: extracting the trajectory dataP K The outlier of (1), the outlier comprising: data missing points and space positioning drift points; aggregating the gridsRPerforming superposition analysis with the abnormal points and marking the grid setRThe grid intersecting the abnormal point is entered to step S5; the superposition analysis specifically comprises the following steps: comparing the coordinates of the abnormal points with the grid set R one by one, and when the coordinates of the abnormal points are positioned in the grid setRThe mesh in (1) is marked when inside or at the boundary of the mesh.
S5: from the set of trajectory dataPRemoving the trajectory dataP K To makeKAdding 1 to the value of (c); determining the trajectory data setPWhether it is empty; if yes, go to step S6; if not, returning to the step S3;
s6: counting the grid setRThe number of times each grid is marked is made to exceed a threshold valueHThe grid output of (a) is an urban canyon region.
In this embodiment, wuhan city is used as a sampling region, a grid set of a road buffer area is established on the basis of wuhan city road network data, 9447362 taxi track data in the wuhan city are collected, and an urban canyon region is extracted by extracting abnormal points in the track data set and performing superposition analysis and statistics with a road grid.
Further, in step S2, the grid setREach grid in (a) is specifically represented as:
wherein,ja number representing a grid;CNT j representing the number of outliers in the grid;xl j represents the longitude of the lower left corner of the grid;yl j representing the latitude of the lower left corner of the grid;xR j represents the longitude of the upper right corner of the grid;yR j representing the latitude of the upper right corner of the grid;ha number representing a set of grids.
In the embodiment, urban roads in Wuhan City are downloaded from OpenStreetMap, contain road grade information, and are specifically divided into four levels of urban expressways, main roads, secondary roads and branches; the road data covers the whole Wuhan city; establishing a buffer area for urban roads in Wuhan city, considering the positioning accuracy of GNSS, setting the distance of the buffer area as 50 meters, dividing the buffer area into grids with the width of 25 meters, and endowing each grid with an attribute fieldCNTThis field stores the long shaping variable.
Further, in step S3, the trajectory dataP K The concrete expression is as follows:
wherein,Tthe maximum sampling times of the track;Knumbering the vehicles;P i is the first in the trackiTracing points;t i is shown asiSampling time of each trace point;v i indicating that the vehicle ist i The speed of travel at that moment;x i to representt i Longitude of the time track point;y i to representt i The latitude of the time trace point.
In the embodiment, the track data is derived from the vehicle history record of the management system of the passenger taxi in Wuhan city, and the track data comprises information such as a vehicle number, longitude and latitude coordinates, data uploading time, whether the vehicle carries passengers or not, a power supply state, a GPS (global positioning system) positioning state and the like; for track data setPThe vehicles are grouped according to the vehicle numbers, and after the vehicles are sequenced according to the data uploading time, different vehicles can be sequencedkThe track data of (1) is recorded as track data increasing according to timeP K 。
Further, in step S4, the step of extracting the data missing point in the outlier specifically includes:
s401: from the trajectory dataP K Middle pass sampling timet n Two adjacent track points are taken out in sequenceP n AndP n+1 ,nis from 1 toT-A positive integer of 1, or a mixture thereof,nis 1; if the track pointP n And track pointP n+1 The sampling time interval is greater than the preset timet m If the data missing point exists, the data missing point and the track point existP n And track pointP n+1 To (c) to (d);
in this embodiment, the sampling time average of the vehicle trajectory data in wuhan cityThe interval is 60 seconds, so that the preset time for judging whether the trace point has the data missing point or not is sett m Set to 60 seconds;
s402: by means of tracing pointsP n And track pointP n+1 Obtaining the coordinates of the data missing point by the coordinate calculation of (2) so thatnThe value of (a) is added to 1, and the formula is specifically as follows:
wherein (A), (B), (C), (D), (C), (B), (C)x,y) Coordinates of the data missing point; (x n ,y n ) Is a track pointP n The coordinates of (a); (x n+1 ,y n+1 ) Tracing pointP n+1 The coordinates of (a);representing points of trackP n And track pointP n+1 The time interval between two adjacent points;t m is a preset time;
s403: repeating step S401 and step S402T-1 time, calculating said trajectory dataP K All of the data missing points in (1).
Further, when judging whether the track points exist in the space positioning drift points, considering that the distance of a running road network between two adjacent track points on one track should be close to the average distance of running on the road where the track points are located in the same sampling intervalSThus, when the distance between two tracing points is greater thanSIn time, the two trace points can be considered to have abnormal drift points; average distanceSThe calculation formula of (2) is as follows:
In this embodiment, different speed thresholds are set for different road grades, and the urban expressway is usedSet to 80km/h, of the trunkSet to 60km/h, secondary trunkSet at 50km/h, branchedSet to 40 km/h;
in step S4, the step of extracting the spatial localization drift point in the outlier specifically includes:
s411: from the trajectory dataP K Middle pass sampling timet n Sequentially taking out three adjacent track pointsP n 、P n+1 AndP n+2 ,nis from 1 toT-A positive integer of 2 is a positive integer of,nis 1;
s412: judging the track pointP n 、P n+1 AndP n+2 the distance between and the average distanceSThe magnitude relationship of (1);
if it isP n AndP n+1 is less than or equal toSAnd is andP n+1 andP n+2 is greater thanSThen point of trackP n+2 For spatially locating drift points and track pointsP n+2 The coordinates of (a) are:
if it isP n AndP n+1 is greater thanSAnd is andP n+1 andP n+2 if the distance between the two points is less than or equal to S, the track points areP n For spatially locating drift points and track pointsP n The coordinates of (a) are:
if it isP n AndP n+1 is greater thanSAnd is andP n+1 andP n+2 if the distance between the two points is greater than S, the track points areP n+1 For spatially locating drift points and track pointsP n+1 The coordinates of (a) are:
wherein (A), (B), (C), (D), (C), (B), (C)x n ,y n ) Is a track pointP n The coordinates of (a); (x n+1 ,y n+1 ) Tracing pointP n+1 The coordinates of (a); (x n+2 ,y n+2 ) Is a track pointP n+2 The coordinates of (a);
s413: repeating the steps S411 and S412 for T-2 times, and calculating the track dataP K All of the spatially located drift points.
Referring to fig. 2, further, step S6 specifically includes:
when the grid is assembledRWhen a grid in (2) is marked once, the gridCNT j Adding 1 to the value of (c); when the grid is inCNT j Value greater than thresholdHThen, outputting the grid as an urban canyon region; in this embodiment, the threshold value is setHIs set to 50, the urban canyon region in wuhan city is finally obtained.
Referring to fig. 3, an urban canyon rapid detection device based on floating car track data includes the following modules:
a track data acquisition module 10, configured to select a sampling area, acquire track data of multiple vehicles in the sampling area, and obtain a track data setP;
A grid set generating module 20, configured to establish a buffer area for the urban road in the sampling area, and generate a grid set corresponding to the urban roadR;
A trajectory data extraction module 30 for extracting the trajectory data setPTrack data of a certain vehicleP K ,KIs a positive integer and has an initial value of 1; entering a superposition analysis module;
an overlay analysis module 40 for extracting the trajectory dataP K The outlier of (1), the outlier comprising: data missing points and space positioning drift points; aggregating the gridsRPerforming superposition analysis with the abnormal points and marking the grid setREntering a grid intersected with the abnormal point into an updating track data set module;
an update track data set module 50 for generating a track data set from the track data setPRemoving the trajectory dataP K To makeKAdding 1 to the value of (c); determining the trajectory data setPWhether it is empty; if yes, entering an urban canyon region judgment module; if not, returning to the track data extraction module;
an urban canyon region determination module 60 for counting the grid setRThe number of times each grid is marked is made to exceed a threshold valueHThe grid output of (a) is an urban canyon region.
Other embodiments or specific implementation manners of the target area detection apparatus according to the present invention may refer to the above method embodiments, and are not described herein again.
It should be noted that, in this document, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or system that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or system. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other like elements in a process, method, article, or system that comprises the element.
The above-mentioned serial numbers of the embodiments of the present invention are merely for description and do not represent the merits of the embodiments. In the unit claims enumerating several means, several of these means may be embodied by one and the same item of hardware. The use of the words first, second, third and the like do not denote any order, but rather the words first, second and the like may be interpreted as indicating any order.
The above description is only a preferred embodiment of the present invention, and not intended to limit the scope of the present invention, and all modifications of equivalent structures and equivalent processes, which are made by using the contents of the present specification and the accompanying drawings, or directly or indirectly applied to other related technical fields, are included in the scope of the present invention.
Claims (5)
1. A floating vehicle track data-based urban canyon rapid detection method is characterized by comprising the following steps:
s1: selecting a sampling area, and carrying out track data acquisition on a plurality of vehicles in the sampling area to obtain a track data setP;
S2: establishing a buffer area for the urban road in the sampling area, and generating a grid set corresponding to the urban roadR;
S3: extracting the set of trajectory dataPTrack data of a certain vehicleP K ,KIs a positive integer and has an initial value of 1; proceeding to step S4;
the trackData ofP K The concrete expression is as follows:
wherein,Tthe maximum sampling times of the track;Knumbering the vehicles;p i is the first in the trackiTracing points;t i is shown asiSampling time of each trace point;v i indicating that the vehicle ist i The speed of travel at that moment;x i to representt i Longitude of the time track point;y i to representt i Latitude of the moment track point;
s4: extracting the trajectory dataP K The outlier of (1), the outlier comprising: data missing points and space positioning drift points; aggregating the gridsRPerforming superposition analysis with the abnormal points and marking the grid setRThe grid intersecting the abnormal point is entered to step S5;
the step of extracting the space positioning drift point in the abnormal points specifically comprises the following steps:
s411: from the trajectory dataP K Middle pass sampling timet n Sequentially taking out three adjacent track pointsP n 、P n+1 AndP n+2 ,nis from 1 toT-A positive integer of 2 is a positive integer of,nis 1;
s412: judging the track pointP n 、P n+1 AndP n+2 the distance between and the average distanceSThe magnitude relationship of (1);
if it isP n AndP n+1 is less than or equal toSAnd is andP n+1 andP n+2 is greater thanSThen point of trackP n+2 For positioning in spaceDrift points, track pointsP n+2 The coordinates of (a) are:
if it isP n AndP n+1 is greater thanSAnd is andP n+1 andP n+2 if the distance between the two points is less than or equal to S, the track points areP n For spatially locating drift points and track pointsP n The coordinates of (a) are:
if it isP n AndP n+1 is greater thanSAnd is andP n+1 andP n+2 if the distance between the two points is greater than S, the track points areP n+1 For spatially locating drift points and track pointsP n+1 The coordinates of (a) are:
wherein (A), (B), (C), (D), (C), (B), (C)x n ,y n ) Is a track pointP n The coordinates of (a); (x n+1 ,y n+1 ) Tracing pointP n+1 The coordinates of (a); (x n+2 ,y n+2 ) Is a track pointP n+2 The coordinates of (a);
s413: repeating the steps S411 and S412 for T-2 times, and calculating the track dataP K All of the spatially located drift points;
s5: from the set of trajectory dataPRemoving the trajectory dataP K To makeKAdding 1 to the value of (c); determining the trajectory data setPWhether it is empty; if yes, go to step S6; if not, returning to the step S3;
s6: counting the grid setRThe number of times each grid is marked is made to exceed a threshold valueHThe grid output of (a) is an urban canyon region.
2. The floating vehicle trajectory data-based urban canyon rapid detection method according to claim 1, wherein in step S2, said grid assemblyREach grid in (a) is specifically represented as:
wherein,ja number representing a grid;CNT j representing the number of outliers in the grid;xl j represents the longitude of the lower left corner of the grid;yl j representing the latitude of the lower left corner of the grid;xR j represents the longitude of the upper right corner of the grid;yR j representing the latitude of the upper right corner of the grid;ha number representing a set of grids.
3. The floating vehicle trajectory data-based urban canyon rapid detection method according to claim 1, wherein in step S4, the step of extracting the data missing point in the outlier specifically comprises:
s401: from the trajectory dataP K Middle pass sampling timet n Two adjacent track points are taken out in sequenceP n AndP n+1 ,nis from 1 toT-A positive integer of 1, or a mixture thereof,nis 1; if the track pointP n And track pointP n+1 The sampling time interval is greater than the preset timet m If the data missing point exists, the data missing point and the track point existP n And track pointP n+1 To (c) to (d);
s402: by means of tracing pointsP n And track pointP n+1 Obtaining the coordinates of the data missing point by the coordinate calculation of (2) so thatnThe value of (a) is added to 1, and the formula is specifically as follows:
wherein (A), (B), (C), (D), (C), (B), (C)x,y) Coordinates of the data missing point; (x n ,y n ) Is a track pointP n The coordinates of (a); (x n+1 ,y n+1 ) Is a track pointP n+1 The coordinates of (a);representing points of trackP n And track pointP n+1 The time interval between two adjacent points;t m is a preset time;
s403: repeating step S401 and step S402T-1 time, calculating said trajectory dataP K All of the data missing points in (1).
4. The floating vehicle trajectory data-based urban canyon rapid detection method according to claim 2, wherein step S6 specifically comprises:
when the grid is assembledRWhen a grid in (2) is marked once, the gridCNT j Adding 1 to the value of (c); when the grid is inCNT j Value greater than thresholdHThe grid is then output as an urban canyon region.
5. The utility model provides a quick detection device of urban canyon based on floating car orbit data which characterized in that includes following module:
a track data acquisition module for selecting a sampling region and comparing the selected sampling region with a corresponding placeAcquiring track data of a plurality of vehicles in the sampling area to obtain a track data setP;
A grid set generating module for establishing a buffer area for the urban road in the sampling area and generating a grid set corresponding to the urban roadR;
A track data extraction module for extracting the track data setPTrack data of a certain vehicleP K ,KIs a positive integer and has an initial value of 1; entering a superposition analysis module;
the trajectory dataP K The concrete expression is as follows:
wherein,Tthe maximum sampling times of the track;Knumbering the vehicles;p i is the first in the trackiTracing points;t i is shown asiSampling time of each trace point;v i indicating that the vehicle ist i The speed of travel at that moment;x i to representt i Longitude of the time track point;y i to representt i Latitude of the moment track point;
a superposition analysis module for extracting the trajectory dataP K The outlier of (1), the outlier comprising: data missing points and space positioning drift points; aggregating the gridsRPerforming superposition analysis with the abnormal points and marking the grid setREntering a grid intersected with the abnormal point into an updating track data set module;
the step of extracting the space positioning drift point in the abnormal points specifically comprises the following steps:
s411: from the trajectory dataP K Middle pass sampling timet n Sequentially taking out three adjacent track pointsP n 、P n+1 AndP n+2 ,nis from 1 toT-A positive integer of 2 is a positive integer of,nis 1;
s412: judging the track pointP n 、P n+1 AndP n+2 the distance between and the average distanceSThe magnitude relationship of (1);
if it isP n AndP n+1 is less than or equal toSAnd is andP n+1 andP n+2 is greater thanSThen point of trackP n+2 For spatially locating drift points and track pointsP n+2 The coordinates of (a) are:
if it isP n AndP n+1 is greater thanSAnd is andP n+1 andP n+2 if the distance between the two points is less than or equal to S, the track points areP n For spatially locating drift points and track pointsP n The coordinates of (a) are:
if it isP n AndP n+1 is greater thanSAnd is andP n+1 andP n+2 if the distance between the two points is greater than S, the track points areP n+1 For spatially locating drift points and track pointsP n+1 The coordinates of (a) are:
wherein (A), (B), (C), (D), (C), (B), (C)x n ,y n ) Is a track pointP n The coordinates of (a); (x n+1 ,y n+1 ) Tracing pointP n+1 The coordinates of (a); (x n+2 ,y n+2 ) Is a track pointP n+2 The coordinates of (a);
s413: repeating the steps S411 and S412 for T-2 times, and calculating the track dataP K All of the spatially located drift points;
an update track data set module for generating a track data set from the track data setPRemoving the trajectory dataP K To makeKAdding 1 to the value of (c); determining the trajectory data setPWhether it is empty; if yes, entering an urban canyon region judgment module; if not, returning to the track data extraction module;
the urban canyon region judgment module is used for counting the grid setRAnd outputting the grids with the marked times exceeding a threshold value as the urban canyon region.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110147584.5A CN112462401B (en) | 2021-02-03 | 2021-02-03 | Urban canyon rapid detection method and device based on floating vehicle track data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110147584.5A CN112462401B (en) | 2021-02-03 | 2021-02-03 | Urban canyon rapid detection method and device based on floating vehicle track data |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112462401A CN112462401A (en) | 2021-03-09 |
CN112462401B true CN112462401B (en) | 2022-02-08 |
Family
ID=74802462
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110147584.5A Active CN112462401B (en) | 2021-02-03 | 2021-02-03 | Urban canyon rapid detection method and device based on floating vehicle track data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112462401B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112685685B (en) * | 2021-03-17 | 2021-06-04 | 中国人民解放军国防科技大学 | Analytic calculation method for discontinuous coverage of regressive orbit satellite constellation to ground |
CN114898559B (en) * | 2022-05-26 | 2023-08-15 | 中国科学院深圳先进技术研究院 | Urban vehicle movement perception capability measuring method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111160480A (en) * | 2019-12-31 | 2020-05-15 | 斑马网络技术有限公司 | Underground parking lot entrance and exit excavation method and device and electronic equipment |
WO2020103558A1 (en) * | 2018-11-19 | 2020-05-28 | 华为技术有限公司 | Positioning method and electronic device |
CN111983655A (en) * | 2020-07-22 | 2020-11-24 | 汉海信息技术(上海)有限公司 | Method and device for determining urban canyon region, electronic device and storage medium |
-
2021
- 2021-02-03 CN CN202110147584.5A patent/CN112462401B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020103558A1 (en) * | 2018-11-19 | 2020-05-28 | 华为技术有限公司 | Positioning method and electronic device |
CN111160480A (en) * | 2019-12-31 | 2020-05-15 | 斑马网络技术有限公司 | Underground parking lot entrance and exit excavation method and device and electronic equipment |
CN111983655A (en) * | 2020-07-22 | 2020-11-24 | 汉海信息技术(上海)有限公司 | Method and device for determining urban canyon region, electronic device and storage medium |
Non-Patent Citations (2)
Title |
---|
基于HMM模型的老年人出行异常检测研究;刘翔;《中国优秀硕士学位论文全文数据库 信息科技辑》;20151015(第10期);正文第4,9-13页 * |
基于改进滑动窗口的渔船AIS轨迹数据压缩算法;董婉婷 等;《大连海洋大学学报》;20200630;第35卷(第3期);第462-468页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112462401A (en) | 2021-03-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110850439B (en) | High-precision three-dimensional point cloud map construction method | |
CN109946731B (en) | Vehicle high-reliability fusion positioning method based on fuzzy self-adaptive unscented Kalman filtering | |
CN103743383B (en) | A kind of Automatic extraction method for road information based on a cloud | |
CN103745018B (en) | Multi-platform point cloud data fusion method | |
JP5162849B2 (en) | Fixed point position recorder | |
CN106197460B (en) | A method of it is predicted with carrying out trip purpose using GPS trip data | |
CN102208013A (en) | Scene matching reference data generation system and position measurement system | |
CN103453913A (en) | Integrated positioning method of vehicle in parking lot | |
CN112462401B (en) | Urban canyon rapid detection method and device based on floating vehicle track data | |
CN104677361A (en) | Comprehensive positioning method | |
CN105740505A (en) | GPS-RTK technology based road space line shape recovery method | |
CN107907134A (en) | A kind of mileage information aids in the matched Vehicle positioning system of earth magnetism and method | |
Santos et al. | A roadside unit-based localization scheme to improve positioning for vehicular networks | |
KR100510835B1 (en) | Method for constituting geographic information system applied digital map using real time measuring systems | |
Gustafsson et al. | Navigation and tracking of road-bound vehicles | |
KR100448543B1 (en) | Method for Preparing Geographical Information System | |
CN1455268A (en) | Satellite locating error processing method based on map matching | |
Hsu et al. | Intelligent viaduct recognition and driving altitude determination using GPS data | |
Liu et al. | Determination of routing velocity with GPS floating car data and webGIS-based instantaneous traffic information dissemination | |
CN117741717A (en) | Global positioning system bias detection and reduction | |
CN105427739A (en) | Kalman filter based road slope enhanced digital map making method | |
Magrath et al. | Evaluating different methods for determining road grade best suited to advanced bus transportation systems | |
Wang et al. | Intersection and stop bar position extraction from crowdsourced gps trajectories | |
CN113360849B (en) | Sea ice density estimation method and system based on ICESat-2 data | |
KR100448054B1 (en) | Method for Preparing Geographical Information System Employing the Amended Value as Road Data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20230419 Address after: 430072 Hubei Province, Wuhan city Wuchang District of Wuhan University Luojiashan Patentee after: WUHAN University Patentee after: Wuhan zhunwang Technology Co.,Ltd. Address before: No.10, 14 / F, building A-2, National Geospatial Information Industry base zone II (phase 6), Wuda science and Technology Park, Donghu New Technology Development Zone, Wuhan City, Hubei Province, 430000 Patentee before: Wuhan zhunwang Technology Co.,Ltd. |