CN107270922B - Traffic accident spatial positioning method based on POI (Point of interest) index - Google Patents

Traffic accident spatial positioning method based on POI (Point of interest) index Download PDF

Info

Publication number
CN107270922B
CN107270922B CN201710451696.3A CN201710451696A CN107270922B CN 107270922 B CN107270922 B CN 107270922B CN 201710451696 A CN201710451696 A CN 201710451696A CN 107270922 B CN107270922 B CN 107270922B
Authority
CN
China
Prior art keywords
traffic
road
information
accident
link
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
Application number
CN201710451696.3A
Other languages
Chinese (zh)
Other versions
CN107270922A (en
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.)
BEIJING TRANSPORTATION INFORMATION CENTER
Beihang University
Original Assignee
BEIJING TRANSPORTATION INFORMATION CENTER
Beihang University
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 BEIJING TRANSPORTATION INFORMATION CENTER, Beihang University filed Critical BEIJING TRANSPORTATION INFORMATION CENTER
Priority to CN201710451696.3A priority Critical patent/CN107270922B/en
Publication of CN107270922A publication Critical patent/CN107270922A/en
Application granted granted Critical
Publication of CN107270922B publication Critical patent/CN107270922B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/26Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network
    • G01C21/34Route searching; Route guidance
    • G01C21/3453Special cost functions, i.e. other than distance or default speed limit of road segments
    • G01C21/3492Special cost functions, i.e. other than distance or default speed limit of road segments employing speed data or traffic data, e.g. real-time or historical
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/09Arrangements for giving variable traffic instructions
    • G08G1/091Traffic information broadcasting

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Multimedia (AREA)
  • Traffic Control Systems (AREA)
  • Navigation (AREA)

Abstract

The invention provides a traffic accident spatial positioning method based on POI (point of interest) index, which is characterized in that POI information is screened according to location information in all accident text broadcast information to obtain the center coordinates of accident locations; performing map matching on the center coordinates of the accident site, and calculating the vertical distance from the center coordinates to the road link to obtain the road link within a fixed radius range; calculating the traffic flow direction of all the road chains running in one direction in the electronic map; acquiring the traffic direction and road grade of the POI to form a POI index table; extracting accident locations and traffic flow directions for the text broadcast information of the traffic accidents; and searching accident places in the constructed index table, and then screening according to the traffic flow direction and the road grade to finally obtain a matching road chain. The invention can position the accident location according to the accident location, the traffic flow direction and the road grade, and has small search range and high matching speed.

Description

Traffic accident spatial positioning method based on POI (Point of interest) index
Technical Field
The invention belongs to the technical field of intelligent traffic, and particularly relates to a traffic accident spatial positioning method based on POI (point of interest) index.
Background
The accident text broadcast information is that when an accident occurs at a certain place, the traffic administration office receives alarm information, and broadcasts accident information on a broadcast station in a manual broadcast mode to inform a driver, so that the traffic jam pressure is relieved. However, the broadcasted text information has ambiguity of human factors because of no fixed broadcasting frame, and the text information cannot be directly reflected on the electronic map, so that the influence of the accident on the road condition is difficult to analyze. The broadcasted accident text information is spatially positioned, so that the accident site can be positioned on the electronic map, the road link number of the accident site is obtained, the spatial range is obtained according to the topological structure of the road link number, and the method is the basis for researching the road condition change caused by the traffic accident.
In the PGIS-based traffic violation accident location rapid positioning method (application number/patent number: 201210200518.0), road space geographic information is set based on PGIS (police geographic information system), and electronic coordinates of the traffic violation accident location are obtained through complicated calculation according to specific kilometer number and meter number input by a user.
In the automatic recording and positioning method for semantic expression information of traffic accident sites (application number/patent number: 201210297477.1), site recording and positioning are realized according to the selection of a user on road sections, side directions, distances, reference points and the like, and the interaction between the user and a system is more, so that the method is not suitable for positioning broadcast text information.
The conventional positioning method for the accident site has more interaction with users, depends on input information of the users, and the input information is not comprehensive enough in manual text broadcast information.
The data used includes POI location data and electronic map information. POI location data: the city-wide city system comprises information of all interest points in a city range, wherein the information comprises traffic place information such as road entrances and exits and overpasses, public facility information such as a gas station, a hospital and a school, and GPS information of life-related places such as restaurants, living halls and entertainment halls, and the total number of the information is 131 ten thousand. Screening is required according to accident location information. Electronic map information: the method comprises all road link information in a city, and each road link is provided with information such as sampling point information, map large grid numbers, road link numbers, road directions and road grades. According to the sampling point data, the traffic flow direction of the one-way road can be calculated.
The method used comprises large grid numbers located into the map according to the GPS coordinates of the traffic location and a fast map matching algorithm based on the organization of the small grid network (application/patent No.: 200610112433.1): the traffic location GPS coordinates can be quickly located into 9 small grids.
The method for positioning the GPS coordinates of the traffic place to a large grid number comprises the following steps: the map large grid number adopts a six-bit coding mode, for example 595673 is a geographical large grid code. The six-bit code may be denoted abcdef. The meaning of the different bits is shown in table 1 below.
TABLE 1 Large trellis encoding scheme
Figure BDA0001322631670000021
The calculation mode of calculating the geographical grid code according to the longitude and latitude is as follows:
LatCode ═ latitude 1.5, then rounding
Longitude code LonCode longitude-60 and rounding
The large grid is divided by taking one eighth of the longitude and latitude code as a unit.
y1 (left decimal part of 8) when rounding is coded and then rounded
x1 is the fractional part left by 8 when longitude code is rounded
Disclosure of Invention
The invention solves the problems: in order to solve the problems of dependence on user input, complex calculation and the like, a spatial positioning method for accident sites based on a POI index table is provided, the accident sites can be positioned according to the accident sites, the traffic flow direction and the road grade, the search range is small, and the matching speed is high.
The technical scheme of the invention is as follows: a traffic accident spatial positioning method based on POI (point of interest) index is characterized in that traffic accident text broadcast information issued by a traffic management bureau is used for acquiring accident occurrence places and traffic flow directions of the occurrence places, an electronic map data and POI place data are combined to construct a POI index table, the road link number of the accident places is searched, and basic support is provided for subsequent researches on space range, duration and the like influenced by traffic accidents. A method flow diagram is shown in fig. 1.
The method comprises the following steps:
the first step is as follows: screening POI (Point of interest) location information, and extracting a GPS coordinate of a specific traffic location;
the second step is that: carrying out map matching on the GPS coordinates of the specific traffic places to obtain a road link set within a range of 30m from the GPS coordinates;
the third step: acquiring road grade according to the map file, and calculating the traffic flow direction of the road;
the fourth step: according to the traffic place and the road link set obtained in the second step, the traffic flow direction and the road grade corresponding to each road link number in the road link set are searched in the third step, and a POI index table is constructed and comprises four fields: traffic location, link number, road grade and traffic direction;
the fifth step: extracting traffic accident information according to the accident broadcasting information;
and a sixth step: matching accident sites and traffic flow directions in the POI index table constructed in the fourth step through the traffic accident information obtained in the fifth step; and finally, obtaining the final space positioning according to the road grade screening result.
The second step is specifically realized as follows:
(21) traversing the GPS point coordinate information of all traffic points in the first step, firstly positioning to a large grid number in a map, and then positioning to a small grid number according to a map matching algorithm of a small grid network organization structure to obtain all link sets in the small grid;
(22) each link in the link set is a line and is formed by connecting a plurality of sampling points. Traversing the road links in the road link set, calculating the spherical distance between the GPS coordinates of the traffic place and each sampling point on the road links, selecting the shortest distance as the distance between the GPS coordinates of the traffic place and the road links, and generating the distance set between the GPS coordinates of the traffic place and the road links;
(23) and if the distance from the GPS coordinates of the traffic place to the road links in the road link set is less than 30m, outputting the traffic place and the corresponding road link numbers.
The third step is specifically realized as follows:
(31) the road links on the electronic map information comprise grid numbers, road link numbers, road grades and passing direction information, the passing direction is represented by 4 numbers which are respectively 0, 1, 2 and 3, 0 represents no investigation and can pass in both directions by default; 1 represents that the road is a bidirectional road, and both directions can pass; 2, representing one-way passing, wherein the passing direction is from the starting point to the end point of the road chain; 3, representing one-way passing, wherein the passing direction is from the end point to the starting point of the road chain, and if the passing direction of the road chain is 2, sequentially acquiring the longitude and latitude of the sampling point of the road chain; if the traffic direction of the road chain is 3, acquiring the longitude and latitude of the sampling point on the road chain number in a reverse order, and if the traffic direction of the road is 0 or 1, not considering the traffic direction;
(32) and after the longitude and latitude of the sampling point on the road link number are obtained, the longitude and latitude information of the first sampling point and the longitude and latitude information of the last sampling point on the road link are taken, and the longitude difference and the latitude difference between two points are calculated. Because the longitude and the latitude of the sampling point are accurate to 5 bits after the decimal point, the default road chain head and tail sampling points cannot have the condition of the same longitude and latitude;
(33) if the absolute value of the longitude difference is larger than the absolute value of the latitude difference, the trend of the road is considered to be the east-west direction, and if the longitude difference is larger than 0, the traffic flow direction of the road link is the east-west direction; if the longitude difference is less than 0, the traffic flow direction of the road link is west to east;
(34) if the absolute value of the latitude difference is greater than the absolute value of the longitude difference, the trend of the road is considered to be the north-south direction, and if the latitude difference is greater than 0, the traffic flow direction of the road chain is the north-south direction; if the latitude difference is less than 0, the traffic flow direction of the road link is from south to north;
in the fifth step, the concrete implementation is as follows:
(51) the text information broadcasted by the accident contains the location information and the traffic flow direction information, and the broadcasted text information is subjected to word segmentation by using an open source word segmentation packet and a traffic inherent word bank;
(52) sequentially traversing the obtained word list, and finding out the position of the accident keywords in the index list, wherein the accident keywords comprise accidents, scratch and rub, faults, scratch and collision, spontaneous combustion, rear-end collision and scattering;
(53) after the accident keyword index position is determined, searching all words before the accident keyword index position in a traffic inherent word bank according to a text habit, if the words exist, obtaining accident information, and reading the next accident text; if the broadcast message does not exist, searching after the index position, if the broadcast message does not exist, acquiring an accident site, and if the broadcast message does not exist, indicating that the broadcast message does not exist in a traffic word bank, and updating the traffic word bank after manual correction;
(54) the obtained accident information includes: location information, traffic direction, lane information.
In the step (22), a calculation formula of the spherical distance between the GPS points is as follows:
the traffic location GPS coordinates are: c (x)0,y0) The coordinates of the sampling point of the road link are A (x)1,y1)
Firstly, converting the longitude and latitude of the GPS into radian:
Figure BDA0001322631670000051
the traffic place radian coordinates are: c (radX)0,radY0) The coordinates of the sampling point of the road chain are A (radX)1,radY1)
Radius of the earth R6378137
α=radY1-radY2
b=radX1-radX2
Figure BDA0001322631670000052
The sixth step is specifically realized as follows:
(61) matching the traffic points in the POI index table established in the fourth step through the traffic accident information acquired in the fifth step, finding a POI index set of the same traffic points, and screening the POI index set according to traffic flow information in the traffic accident information;
(62) and determining road grade according to the occurrence place of the traffic accident information, and screening a result according to the road grade to obtain final space positioning.
Compared with the prior art, the invention has the advantages that:
(1) the input information of the interactive accident site is not needed to be carried out with the broadcasting personnel.
(2) And acquiring the road link number of the accident occurrence place without carrying out complex calculation.
(3) Only need be based on electronic map and POI data information, need not resort to other platforms.
Drawings
FIG. 1 is a flow chart of a method for spatial location of traffic accidents based on POI indexing;
FIG. 2 is a schematic diagram of POI point location to a small grid;
FIG. 3 is a schematic diagram illustrating a distance between a POI point and a road link sampling point;
FIG. 4 is a flow chart of a traffic direction acquisition method;
fig. 5 is a flow chart for acquiring traffic accident information.
Detailed Description
As shown in fig. 1, the spatial positioning method of traffic accidents based on POI index of the present invention is as follows:
the first step is as follows: and acquiring GPS coordinates of the traffic place center. And screening POI location information according to the road information, and extracting the GPS point coordinates of the specific traffic location.
The second step is that: traffic location center coordinate map matching
(1) Traversing the GPS information of all the locations obtained in the first step, performing map large grid positioning on the GPS location, and positioning to 9 small grids related to the GPS location according to the small grid network organization structure (see 200610112433.1), as shown in fig. 2.
(2) And calculating the distance between the point and all the links in the nine small grids. Each road link is provided with GPS sampling point information, all GPS sampling points are obtained according to the road link number, the spherical distance between the traffic place center GPS and each GPS sampling point is calculated, and the shortest distance is selected as the distance from the traffic place GPS to the road link. The spherical distance between the GPS points is calculated as follows:
(3) the distance between the point and all the links in the 9 small grids is calculated. Each road link has sampling point information, all sampling points are obtained according to the road link number, the spherical distance between the GPS coordinate of the traffic place and each sampling point is calculated, and the shortest distance is selected as the distance from the GPS coordinate of the traffic place to the road link. The calculation formula of the spherical distance between the GPS coordinates of the traffic points and the sampling points is as follows:
as shown in fig. 3, the traffic location GPS coordinates are: c (x)0,y0) And the coordinates of the road link sampling points are respectively as follows:
A1(x1,y1),A2(x2,y2),A3(x3,y3),A4(x4,y4) The longitude and latitude of the GPS are converted into radian first when the radius R of the earth is 6378137:
Figure BDA0001322631670000071
the traffic place radian coordinates are: c (radX)0,radY0) Road link sampling point A1The radian coordinates of (a) are: (radX)1,radY1),
a=radY1-radY2
b=radY1-radY2
Figure BDA0001322631670000072
Because the road link contains a plurality of sampling points, a plurality of distance values exist between the GPS coordinates of the traffic point and the sampling points, and the distance between the traffic point C and the road link A is defined as: distaneCA. Traffic site C and road link sampling point A1,A2,A3,A4The distances between are respectively: d1,d2,d3,d4. Among these distances d2The shortest distance, so d is selected2As the distance between the traffic point C and the link a.
DISTANCECA=d2
(3) And if the distance from the traffic place center GPS to the road link is less than 30m, outputting the traffic place and the corresponding road link number.
The third step: map link traffic direction acquisition
(1) The link on the electronic map information comprises the following attribute information: the information of grid number, road link number, road grade, traffic direction and the like. If the passing direction of the road is 2 or 3, the road is a one-way road, and if the passing direction of the road is 2, the longitude and latitude of the GPS sampling point on the road link number are sequentially acquired; and if the road passing direction is 3, reversely acquiring the longitude and latitude of the GPS sampling point on the road link number. If the passing direction of the road is 0 or 1, the road is a bidirectional road, and the traffic direction is not considered. Some road trends are: northwest, northeast, southwest, southeast, these directions do not consider, only consider the east-west trend and the south-north trend.
(2) And sequentially acquiring information of all GPS sampling points on the road link, and acquiring longitude and latitude information of a first sampling point and longitude and latitude information of a last sampling point on the road link. The longitude difference and the latitude difference between the two points are calculated. The GPS longitude and latitude are accurate to 5 bits after decimal point, and the GPS sampling points at the head and the tail of the road link do not have the same longitude and latitude.
(3) If the absolute value of the longitude difference is larger than the absolute value of the latitude difference, the general trend of the road is considered to be the east-west direction, and if the longitude difference is larger than 0, the traffic flow direction of the road link is the east-west direction; if the longitude difference is less than 0, the traffic direction of the road link is west to east.
(4) If the absolute value of the latitude difference is larger than the absolute value of the longitude difference, the general trend of the road is considered to be the north-south direction, and if the latitude difference is larger than 0, the traffic flow direction of the road chain is the north-south direction; if the latitude difference is less than 0, the traffic direction of the road link is north-south.
The flow chart for calculating the traffic flow direction is shown in fig. 4.
The fourth step: building POI index table
According to the traffic place and the road link number obtained in the second step, the traffic direction and the road grade corresponding to the road link number are searched in the map electronic information of the traffic direction obtained in the third step, and a POI index table is constructed and comprises four fields: traffic location, link number, road grade, and traffic direction.
The fifth step: and selecting traffic accident information from the broadcast information. And segmenting the text information according to the word segmentation packet and a word bank of the traffic place to obtain a word list. Searching the keyword position related to the accident in the word list, wherein the keyword related to the accident is represented by the following statistics on the text information of the accident: accidents, faults, scraping and collision, spontaneous combustion, rear-end collision and scattering. After the position of the keyword is obtained, searching the place and the traffic flow direction before the position of the keyword, and if the accident place and the traffic flow direction are obtained, stopping searching; otherwise, searching is carried out after the position of the keyword, if the result is searched, output is obtained, if the text information is not clear enough, the traffic place or the traffic flow direction cannot be extracted, the accident cannot be spatially positioned, and the next accident information is read. The flow chart is shown in fig. 5.
And a sixth step: and searching a corresponding road link number meeting the conditions in the POI index table according to the accident location, the traffic flow direction and the road grade. If the matching road link is searched, the positioning is successful, if the appointed road link number is not searched in the POI index table, the place information of the accident is lost in the POI, and the next accident information is read.
The present invention will be described in detail with reference to specific examples.
1. Screening in POI data according to the urban traffic accident location to obtain a POI location; for example, the meadow bridge: the GPS point coordinates are acquired as follows: (116.3804,39.96817).
2. And positioning to a small grid number according to a GPS coordinate of a meadow bridge center and a map matching algorithm based on a small grid road network organization structure, and acquiring all road chain sets in the small grid.
The link in the set is as follows (truncated) from the central GPS coordinate:
road and chain number Distance between two adjacent plates
59567350173911 4.9115593
59567300681052 5.5821333
59567350173914 6.9895267
59567300677000 7.1791806
59567312676374 7.7923665
59567300666673 8.9055573
59567350173912 12.966893
59567350173913 14.924163
59567300671462 19.464306
59567300676159 19.533117
59567387757507 20.257863
59567300680551 20.774972
59567300680554 22.328519
59567314095212 23.377091
59567300671454 28.050505
59567387757517 28.050505
59567314095210 30.649508
59567312676373 30.860383
59567300671461 33.042658
59567387757516 34.571205
59567300676160 36.465611
59567300687289 36.816878
59567300680583 42.199368
The link number with a distance of less than 30m is selected as shown in the following table:
location of a site Road and chain number
Meadow bridge 59567300666673
Meadow bridge 59567300671454
Meadow bridge 59567300671462
Meadow bridge 59567300676159
Meadow bridge 59567300677000
Meadow bridge 59567300680551
Meadow bridge 59567300680554
Meadow bridge 59567300681052
Meadow bridge 59567312676374
Meadow bridge 59567314095212
Meadow bridge 59567350173911
Meadow bridge 59567350173912
Meadow bridge 59567350173913
Meadow bridge 59567350173914
Meadow bridge 59567387757507
Meadow bridge 59567387757517
3. And acquiring the road traffic flow direction. For example, the link number: 59567300680551, the traffic direction is 3, and the sampling points are (116.38062,39.96809), (116.38123,39.96811), (116.38187,39.96813) obtained in reverse order in the electronic map file.
(1) Calculating the difference in longitude between the first and last sample points
Δlng=116.38062-116.38187=-0.00125
(2) Calculating the latitude difference between the first point and the last sampling point
Δlat=39.96809-39.96813=-0.00004
(3) So that the link runs east-west
4. And (3) extracting the traffic flow direction and the road grade corresponding to the electronic map information according to the road link number acquired in the step (2) to form a POI index table, wherein the POI index table comprises: location, link number, road grade, traffic direction. The Madian bridge index Table is as follows:
location of a site Road and chain number Road grade Direction of traffic flow
Meadow bridge 59567300666673 0102|010c W-E
Meadow bridge 59567300671454 030a S-N
Meadow bridge 59567300671462 030a N-S
Horse countyBridge with a bridge body 59567300676159 030a S-N
Meadow bridge 59567300677000 0102|010c S-N
Meadow bridge 59567300680551 0102|010c W-E
Meadow bridge 59567300680554 0102|010c E-W
Meadow bridge 59567300681052 0102|010c N-S
Meadow bridge 59567312676374 0102|010c E-W
Meadow bridge 59567314095212 030a W-E
Meadow bridge 59567350173911 0002|0003|000c N-S
Meadow bridge 59567350173912 0002|000c N-S
Meadow bridge 59567350173913 0002|000c S-N
Meadow bridge 59567350173914 0002|0003|000c S-N
Meadow bridge 59567387757507 030a E-W
Meadow bridge 59567387757517 030a W-E
5. Given the accident information: 2016/5/1015: 04:16, northeast of the tricyclic Madian bridge, accident of the inner lane, and tail of the rear fleet of vehicles to the Jimen bridge.
(1) Firstly, performing word segmentation processing to obtain a word segmentation list: tricyclic north, meadow bridge, west to east, direction, inner lane, accident, rear car, team tail, row, arrival, jimen bridge.
(2) Finding the accident keyword "accident", searching the place information in front of the word, which respectively comprises: tricyclic north, meadow bridge, west east, inboard lane.
(3) Searching location information in all POI index tables to obtain matched locations: a meadow bridge;
the direction information table is (east to west, west to east, south to north, north to south), and the matched traffic direction is obtained: west to east.
6. Space positioning:
(1) firstly, according to the traffic flow direction, all west-east road chains are obtained
Figure BDA0001322631670000111
Figure BDA0001322631670000121
(2) Then, the main link chains 59567300666673 and 59567300680551 are selected because the text message does not include the sub road information according to the road grade.
The above examples are provided only for the purpose of describing the present invention, and are not intended to limit the scope of the present invention. The scope of the invention is defined by the appended claims. Various equivalent substitutions and modifications can be made without departing from the spirit and principles of the invention, and are intended to be within the scope of the invention.

Claims (1)

1. A traffic accident space positioning method based on POI index is characterized by comprising the following steps:
the first step is as follows: screening POI location information, and extracting GPS coordinates of specific traffic locations;
the second step is that: carrying out map matching on the GPS coordinates of the specific traffic places to obtain a road link set within a range of 30m from the GPS coordinates;
the third step: acquiring road grade according to the map file, and calculating the traffic flow direction of the road;
the fourth step: according to the traffic place and the road link set obtained in the second step, the traffic flow direction and the road grade corresponding to each road link number in the road link set are searched in the third step, and a POI index table is constructed and comprises four fields: traffic location, link number, road grade and traffic direction;
the fifth step: extracting traffic accident information according to the accident broadcasting information;
and a sixth step: matching the traffic places and the traffic flow direction in the POI index table constructed in the fourth step through the traffic accident information obtained in the fifth step; finally, according to the screening result of the road grade, obtaining the final space positioning;
the second step is specifically realized as follows:
(21) traversing the GPS coordinate information of all traffic places in the first step, firstly positioning to a large grid number in a map, then positioning to a small grid number according to a map matching algorithm of a small grid network organization structure, and acquiring all link sets in the small grid;
(22) each link in the link set is a line and is formed by connecting a plurality of sampling points, the links in the link set are traversed, the spherical distance between the GPS coordinate of the traffic place and each sampling point on the link is calculated, the shortest distance is selected as the distance from the GPS coordinate of the traffic place to the link, and the distance set from the GPS coordinate of the traffic place to the link is generated;
(23) if the distance from the GPS coordinates of the traffic place to the road links in the road link set is less than 30m, outputting the traffic place and the corresponding road link numbers;
the third step is specifically realized as follows:
(31) the road links on the electronic map information comprise grid numbers, road link numbers, road grades and passing direction information, the passing direction is represented by 4 numbers which are respectively 0, 1, 2 and 3, 0 represents no investigation and can pass in both directions by default; 1 represents that the road is a bidirectional road, and both directions can pass; 2, representing one-way passing, wherein the passing direction is from the starting point to the end point of the road chain; 3, representing one-way passing, wherein the passing direction is from the end point to the starting point of the road chain, and if the passing direction of the road chain is 2, sequentially acquiring the longitude and latitude of the sampling point of the road chain; if the passing direction of the road link is 3, acquiring the longitude and latitude of the sampling point on the road link number in a reverse order; if the passing direction of the road is 0 or 1, the traffic flow direction is not considered;
(32) after the longitude and latitude of the sampling point on the road link number are obtained, the longitude and latitude information of the first sampling point and the longitude and latitude information of the last sampling point on the road link are taken, and the longitude difference and the latitude difference between two points are calculated; because the longitude and the latitude of the sampling point are accurate to 5 bits after the decimal point, the default road chain head and tail sampling points cannot have the condition of the same longitude and latitude;
(33) if the absolute value of the longitude difference is larger than the absolute value of the latitude difference, the trend of the road is considered to be the east-west direction, and if the longitude difference is larger than 0, the traffic flow direction of the road link is the east-west direction; if the longitude difference is less than 0, the traffic flow direction of the road link is west to east;
(34) if the absolute value of the latitude difference is greater than the absolute value of the longitude difference, the trend of the road is considered to be the north-south direction, and if the latitude difference is greater than 0, the traffic flow direction of the road chain is the north-south direction; if the latitude difference is less than 0, the traffic flow direction of the road link is from south to north;
in the fifth step, the concrete implementation is as follows:
(51) the text information broadcasted by the accident contains the location information and the traffic flow direction information, and the broadcasted text information is subjected to word segmentation by using an open source word segmentation packet and a traffic inherent word bank;
(52) sequentially traversing the obtained word list, and finding out the position of the accident keywords in the word list, wherein the accident keywords comprise accidents, scratch and rub, faults, scratch and collision, spontaneous combustion, rear-end collision and scattering;
(53) after the accident keyword index position is determined, searching all words before the accident keyword index position in a traffic inherent word bank according to a text habit, if the words exist, obtaining accident information, and reading the next broadcast text; if the broadcast message does not exist, searching after the index position, if the broadcast message does not exist, acquiring a traffic place, if the broadcast message does not exist, indicating that the position broadcasted by the drill message book does not exist in a traffic inherent word bank, and updating the traffic inherent word bank after manual correction;
(54) the obtained accident information includes: location information, traffic direction, lane information;
in the step (22), a calculation formula of the spherical distance between the GPS points is as follows:
the traffic location GPS coordinates are: c (x)0,y0) GPS coordinate of road link sampling point is A (x)1,y1)
Firstly, converting the longitude and latitude of the GPS into radian:
Figure FDA0002362351590000031
the traffic place radian coordinates are: c (radX)0,radY0) And the radian coordinate of the sampling point of the road chain is A (radX)1,radY1)
Radius of the earth R6378137 m
a=radY0-radY1
b=radX0-radX1
Figure FDA0002362351590000032
The sixth step is specifically realized as follows:
(61) matching the traffic points in the POI index table established in the fourth step through the traffic accident information acquired in the fifth step, finding a POI index set of the same traffic points, and screening the POI index set according to the traffic flow direction in the traffic accident information;
(62) and determining road grade according to the location information of the traffic accident information, and screening a result according to the road grade to obtain final space positioning.
CN201710451696.3A 2017-06-15 2017-06-15 Traffic accident spatial positioning method based on POI (Point of interest) index Active CN107270922B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710451696.3A CN107270922B (en) 2017-06-15 2017-06-15 Traffic accident spatial positioning method based on POI (Point of interest) index

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710451696.3A CN107270922B (en) 2017-06-15 2017-06-15 Traffic accident spatial positioning method based on POI (Point of interest) index

Publications (2)

Publication Number Publication Date
CN107270922A CN107270922A (en) 2017-10-20
CN107270922B true CN107270922B (en) 2020-04-14

Family

ID=60067724

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710451696.3A Active CN107270922B (en) 2017-06-15 2017-06-15 Traffic accident spatial positioning method based on POI (Point of interest) index

Country Status (1)

Country Link
CN (1) CN107270922B (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109990791B (en) * 2017-12-30 2021-04-13 中国移动通信集团辽宁有限公司 Method, apparatus, device and medium for road data extraction
CN110160535B (en) * 2018-02-11 2020-11-17 阿里巴巴(中国)有限公司 Road matching method and device
CN108320515B (en) * 2018-04-04 2021-07-02 公安部交通管理科学研究所 Automatic matching and checking method and system for road network of traffic accident site
CN109215338B (en) * 2018-06-27 2020-05-29 杭州叙简科技股份有限公司 Gridding-based Beidou gps vehicle track management system and method thereof
CN109841057A (en) * 2018-11-07 2019-06-04 厦门卫星定位应用股份有限公司 A kind of concern car tracing method based on GNSS
CN109993971B (en) * 2019-03-27 2021-09-17 江苏智通交通科技有限公司 Method for improving positioning accuracy of traffic accident site
CN110060472B (en) * 2019-04-09 2020-12-15 浙江大华技术股份有限公司 Road traffic event positioning method, system, readable storage medium and device
CN110515942B (en) * 2019-07-12 2023-08-04 同济大学 Storage and retrieval method of serialized lane line map
CN110375762B (en) * 2019-08-13 2021-06-11 福建工程学院 Method and device for assisting in planning navigation path
CN113701763B (en) * 2020-05-20 2024-06-21 百度在线网络技术(北京)有限公司 Method and device for generating information
CN113570107A (en) * 2021-06-08 2021-10-29 众能联合数字技术有限公司 Project address positioning method for project rental scene

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1908588A (en) * 2006-08-17 2007-02-07 北京航空航天大学 Fast map matching method based on small lattice road network organization and structure
JP2008157885A (en) * 2006-12-26 2008-07-10 Pioneer Electronic Corp Information guide device, navigation device, information guide method, navigation method, information guide program, navigation program, and recording medium
JP2009133666A (en) * 2007-11-29 2009-06-18 Aisin Aw Co Ltd Navigation device
CN101807196A (en) * 2010-01-08 2010-08-18 北京世纪高通科技有限公司 Method and device for searching information
CN102364545A (en) * 2011-10-12 2012-02-29 北京世纪高通科技有限公司 Method, device and system for inquiring road conditions
CN102521969A (en) * 2011-12-20 2012-06-27 北京世纪高通科技有限公司 Positioning method of traffic event information and apparatus thereof
CN104183144A (en) * 2013-05-20 2014-12-03 北京百度网讯科技有限公司 Real-time traffic condition information generating method and system thereof
CN104820721A (en) * 2015-05-26 2015-08-05 厦门雅迅网络股份有限公司 Method for broadcasting traffic emergency accidents based on landmarks and mileages
CN105787110A (en) * 2016-03-23 2016-07-20 苏州数字地图信息科技股份有限公司 Positioning method and device based on spatial search
CN106504746A (en) * 2016-10-28 2017-03-15 普强信息技术(北京)有限公司 A kind of method for extracting structuring traffic information from speech data

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1908588A (en) * 2006-08-17 2007-02-07 北京航空航天大学 Fast map matching method based on small lattice road network organization and structure
JP2008157885A (en) * 2006-12-26 2008-07-10 Pioneer Electronic Corp Information guide device, navigation device, information guide method, navigation method, information guide program, navigation program, and recording medium
JP2009133666A (en) * 2007-11-29 2009-06-18 Aisin Aw Co Ltd Navigation device
CN101807196A (en) * 2010-01-08 2010-08-18 北京世纪高通科技有限公司 Method and device for searching information
CN102364545A (en) * 2011-10-12 2012-02-29 北京世纪高通科技有限公司 Method, device and system for inquiring road conditions
CN102521969A (en) * 2011-12-20 2012-06-27 北京世纪高通科技有限公司 Positioning method of traffic event information and apparatus thereof
CN104183144A (en) * 2013-05-20 2014-12-03 北京百度网讯科技有限公司 Real-time traffic condition information generating method and system thereof
CN104820721A (en) * 2015-05-26 2015-08-05 厦门雅迅网络股份有限公司 Method for broadcasting traffic emergency accidents based on landmarks and mileages
CN105787110A (en) * 2016-03-23 2016-07-20 苏州数字地图信息科技股份有限公司 Positioning method and device based on spatial search
CN106504746A (en) * 2016-10-28 2017-03-15 普强信息技术(北京)有限公司 A kind of method for extracting structuring traffic information from speech data

Also Published As

Publication number Publication date
CN107270922A (en) 2017-10-20

Similar Documents

Publication Publication Date Title
CN107270922B (en) Traffic accident spatial positioning method based on POI (Point of interest) index
CN112133090A (en) Multi-mode traffic distribution model construction method based on mobile phone signaling data
CN110298500B (en) Urban traffic track data set generation method based on taxi data and urban road network
CN112015835B (en) Geohash compressed map matching method
CN105809962A (en) Traffic trip mode splitting method based on mobile phone data
CN109359162B (en) GIS-based school site selection method
CN106197460B (en) A method of it is predicted with carrying out trip purpose using GPS trip data
CN107133700A (en) Mobile phone signaling data road network method based on R* tree indexes
CN110796337B (en) System for evaluating service accessibility of urban bus stop
CN111931998B (en) Individual travel mode prediction method and system based on mobile positioning data
CN105512166A (en) Traffic parallel method with mapping between microblog public sentiments and city road conditions
CN110160538B (en) Map matching method based on mobile phone signaling data
CN105844031B (en) A kind of urban transportation gallery recognition methods based on mobile phone location data
CN111625732B (en) Address matching method and device
Sun et al. How the natural environment in downtown neighborhood affects physical activity and sentiment: Using social media data and machine learning
CN114331232B (en) Street space quality monitoring, evaluating and early warning method
CN113988659B (en) Three-dimensional compact digital city design method and device and computer equipment
Yao et al. Trip segmentation and mode detection for human mobility data
CN115129769A (en) Resident travel survey sample expansion method and device and storage medium
CN110442660B (en) Bus network length calculation method for intelligent bus system
CN108346299B (en) Vehicle speed evaluation method and device
Masoumi Transformation of urban form and the effects on travel behavior in Iran
Sekimoto et al. Digital archiving of people flow by recycling large-scale social survey data of developing cities
Rahımı et al. Monitoring the Long Term Urban Expansion of Zanjan City Using Remote Sensing and Geographic Information Systems
CN117557011B (en) Method and device for measuring and calculating average distance of travel of highway

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