CN109598925B - Taxi gathering alarm method, terminal equipment and storage medium - Google Patents

Taxi gathering alarm method, terminal equipment and storage medium Download PDF

Info

Publication number
CN109598925B
CN109598925B CN201710916836.XA CN201710916836A CN109598925B CN 109598925 B CN109598925 B CN 109598925B CN 201710916836 A CN201710916836 A CN 201710916836A CN 109598925 B CN109598925 B CN 109598925B
Authority
CN
China
Prior art keywords
taxi
longitude
latitude
data values
data value
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
CN201710916836.XA
Other languages
Chinese (zh)
Other versions
CN109598925A (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.)
Xiamen Yaxon Networks Co Ltd
Original Assignee
Xiamen Yaxon Networks Co Ltd
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 Xiamen Yaxon Networks Co Ltd filed Critical Xiamen Yaxon Networks Co Ltd
Priority to CN201710916836.XA priority Critical patent/CN109598925B/en
Publication of CN109598925A publication Critical patent/CN109598925A/en
Application granted granted Critical
Publication of CN109598925B publication Critical patent/CN109598925B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0125Traffic data processing
    • G08G1/0133Traffic data processing for classifying traffic situation

Abstract

The invention discloses a taxi vehicle gathering alarm method based on road network and station information, which comprises the following steps of S1-S4: determining a mapping method from latitude and longitude data values to reference point data values, S5-S6: mapping longitude and latitude data values of road nodes in an operation area into reference point data values, recording and storing the reference point data values and the occurrence times thereof to form a road network information database DBwayS7-S8: mapping the longitude and latitude data values of taxi stations in the operation area into reference point data values, recording and storing the reference point data values, and forming a station information database DBsiteAnd S9: binding DBwayAnd DBsiteAnd (4) giving an alarm for the abnormal gathering condition of the taxi vehicles according to the data value. The latitude and longitude data values are mapped into the datum point data values, so that the fault tolerance of the supervision platform is improved, and the processing efficiency of the platform is guaranteed. The system is integrated with road network and site information, is associated with historical aggregation conditions, strengthens the triggering conditions of aggregation alarm and effectively reduces the false alarm rate.

Description

Taxi gathering alarm method, terminal equipment and storage medium
Technical Field
The invention relates to the technical field of taxi monitoring, in particular to a taxi vehicle gathering alarm method based on road network and site information, terminal equipment and a storage medium.
Background
The taxi is a transportation tool for providing individual door-to-door convenient service for social public with special travel demands in modes of tour, site, reservation and the like. In China, the vehicle mainly has tour service, so that certain negative externality exists. If necessary measures are not taken for effective supervision, once the negative externality is maliciously amplified, the normal social order can be seriously disturbed, for example, the vicious events of 'road congestion caused by taxi aggregation' frequently occurring in a plurality of cities in China in recent years are the like.
In the prior art, chinese patent publication No. CN102117532A, entitled "method for implementing illegal aggregation pre-alarm of taxis based on GPS system", uniformly divides the area under jurisdiction into grids, regularly counts the number of vehicles whose traveling speeds meet the requirements in each grid, and sends out an alarm when the number exceeds a preset value. The method ignores that the taxi waiting at the station such as an airport, a railway station and the like generally has lower running speed and a plurality of vehicles, which can generate false alarm of legal vehicle aggregation; meanwhile, the set running speed preset value is used as a preposed judgment condition for the vehicle to participate in the gathering event, and if the preset running speed preset value is mastered in advance, the preset running speed preset value is easy to become a vulnerability for avoiding supervision.
Disclosure of Invention
In order to solve the problems, the invention provides a taxi vehicle gathering alarm method, terminal equipment and a storage medium based on road network and station information, wherein phi and psi are used for respectively representing longitude values and latitude values of any positioning points on the earth surface, wherein east longitude is positive, west longitude is negative, north latitude is positive, south latitude is negative, and longitude and latitude of a positioning position are used as spherical data for modeling.
The invention relates to a taxi vehicle gathering alarm method based on road network and station information, which comprises the following steps:
s1: longitude and latitude data value of arbitrary positioning point A in taxi operation area
Figure BDA0001425954170000021
According to the longitude and latitude data value of the positioning point A, the longitude and latitude value of a point B1 with the length distance h along the longitude line from the positioning point A is calculated
Figure BDA0001425954170000022
The longitude and latitude value of a point B2 with the length distance h from the positioning point A along the latitude line
Figure BDA0001425954170000023
Wherein h is a fixed distance step length, and the step S2 is entered;
s2: selecting longitude and latitude data value (phi) of a positioning point C on the earth surfaceanyany) Then, againDetermining an integer N greater than 1, and proceeding to step S3;
s3: according to phi and psi obtained by calculation in S1, a group of integers p and q are taken on the basis of S2, wherein p is more than or equal to 0 and is less than N, and q is more than or equal to 0 and is less than N, so that
Figure BDA0001425954170000024
Thereby obtaining a plurality of reference points with data values of (phi)0+kφ,Ψ0+ l ψ) where k and l are each satisfied
Figure BDA0001425954170000025
Step (3), the process proceeds to step S4;
s4: for any set of k and l meeting the requirements in S3, the latitude and longitude ranges will be satisfied
Figure BDA0001425954170000026
The latitude and longitude data values (phi, psi) of the positioning point D of (1) are all mapped into a data value (phi) represented by the data value of the reference point through a many-to-one correspondence0+kφ,Ψ0+ l ψ), wherein the many-to-one correspondence formula is
Figure BDA0001425954170000027
Wherein [ ] is the operator of rounding, go to step S5;
s5: according to the urban map data of the operation area, acquiring longitude and latitude data values corresponding to road nodes, performing spherical equal-division interpolation between adjacent nodes of each road to enable the distance between adjacent nodes in the interpolated road to be smaller than a preset value r, storing all the longitude and latitude data values of the interpolated road nodes, and entering the step S6;
s6: go through all possible combination values of p and q, according toMapping the longitude and latitude data values of the road nodes stored in the S5 into corresponding reference point data values according to the corresponding relation in the S4, aiming at the combined value of each p and q, counting the occurrence times of each reference point data value after all the mapping is finished, recording and storing each reference point data value and the occurrence times thereof, and forming a road network information database DBwayThen, the process proceeds to step S7;
s7: obtaining the position of each taxi station according to the city map data of the operation area, storing the longitude and latitude data value corresponding to the station position, and entering the step S8;
s8: traversing all possible combination values of p and q, respectively mapping the latitude and longitude data values stored in S7 to corresponding reference point data values according to the corresponding relation in S4 aiming at each combination value of p and q, recording and storing the reference point data values after all mappings are finished, and forming a site information database DBsiteThen, the process proceeds to step S9;
s9: according to the road network information database DBwayAnd a site information database DBsiteAnd (4) giving an alarm for the abnormal gathering condition of the taxi vehicles according to the data value.
Further, in S1, the calculation formula of (ψ, Φ) is:
Figure BDA0001425954170000031
where π is the circumference ratio and R is the earth's mean radius.
Further, in S5, performing equal spherical interpolation between adjacent nodes of each road to make the distance between adjacent nodes in the interpolated road smaller than a preset value r, and storing all longitude and latitude data values of the interpolated road nodes, specifically:
s51, using formula
Figure BDA0001425954170000032
Calculating road adjacent node (phi)AA) And (phi)BB) Spherical distance D between0Wherein Δ ═ arccos (sin (Ψ)A)sin(ΨB)+cos(ΨA)cos(ΨB)cos(ΦAB));
S52, spherical distance D obtained according to S510Calculating the number of equal parts of the equal-part interpolation
Figure BDA0001425954170000033
Wherein
Figure BDA0001425954170000041
Operator for rounding in infinite direction;
s53, using formula
Figure BDA0001425954170000042
Converting the longitude and latitude data value of the road node spherical coordinate into rectangular coordinate, and then converting the longitude and latitude data value (phi) of the road adjacent nodeAA) And (phi)BB) Can be respectively converted into (x)A,yA,zA) And (x)B,yB,zB);
S54, using formula
Figure BDA0001425954170000043
The road adjacent node (x) obtained in S53 is calculatedA,yA,zA) And (x)B,yB,zB) The linear distance therebetween;
s55, using formula
Figure BDA0001425954170000044
Calculating the spatial weight of spherical equal-division interpolation between adjacent nodes, wherein i is 1,2, …, n-1;
s56, using formula
Figure BDA0001425954170000045
(x) obtained in S54A,yA,zA) And (x)B,yB,zB) All spatial linear interpolation points in between
Figure BDA0001425954170000046
S57, using formula
Figure BDA0001425954170000047
The spatial linear interpolation point obtained in S56
Figure BDA0001425954170000048
Projecting the spherical surface along the center of the sphere to the direction of the spherical surface to obtain the rectangular coordinate (x) of the spherical interpolation pointAB,i,yAB,i,zAB,i) Wherein
Figure BDA0001425954170000049
S58, using formula
Figure BDA0001425954170000051
Wherein R is the mean radius of the earth, R1=Rcos(Ψ)
Converting the rectangular coordinates of the road nodes into longitude and latitude data values of spherical coordinates, then obtaining (x) in S57AB,i,yAB,i,zAB,i) Can be converted into (phi)AB,iAB,i);
And S59, storing all longitude and latitude data values of the road nodes after interpolation.
Further, in S7, each taxi station includes, but is not limited to: airports, train stations and bus stations.
Further, in S9, the road network information database DB is usedwayAnd a site information database DBsiteThe data value of (2) gives an alarm to abnormal gathering conditions of taxi vehicles, and specifically comprises the following steps:
s91, acquiring the real-time positioning position of the taxi operated at the jth moment (wherein j is 1,2, …), and storing the corresponding taxi ID number and the longitude and latitude data value of the real-time positioning position;
s92, traversing all possible combination values of p and q, mapping the latitude and longitude data values stored in S91 to corresponding datum point data values according to the corresponding relation in S4 aiming at each combination value of p and q, and counting the occurrence times of each datum point data value after all mappings are finished;
s93, if the times of the datum point data values in the S92 are all smaller than the non-negative preset value X, enabling the vehicle at the jth moment to gather the pre-warning database DBtaxi,jIf the position is empty, increasing j by 1 and returning to S91, otherwise, recording the reference point data value with the occurrence frequency larger than X, and mapping the occurrence frequency and the positioning position into a taxi ID number set corresponding to the reference point data value;
s94, for each reference point data value recorded in S93, respectively in the road network information database DBwayAnd a site information database DBsiteIf it is in DBwayIf the data value exists in the database, the beta is taken as the number of times of the data value of the corresponding datum point, otherwise the data value is taken as zero, if the data value exists in the database, the data value of the corresponding datum point is taken as the number of times of the data value of the corresponding datum point, and ifsiteIf yes, the value of gamma is 1, otherwise the value of gamma is 0;
s95, if the occurrence times of the datum point data values recorded in S93 are all smaller than the corresponding (X + beta Y + gamma Z) values searched in S94 (wherein Y, Z is a non-negative preset value), making DB have non-negative valuestaxi,jIf the value is null, increasing j by 1 and returning to S91, otherwise, recording the datum point data value with the occurrence frequency greater than the corresponding (X + beta Y + gamma Z) value, and mapping the occurrence frequency and the positioning position of the datum point data value into a taxi ID number set corresponding to the datum point data value;
s96, the reference point data value recorded in S95 is used
Figure BDA0001425954170000061
Converting the datum data values (phi ', psi') into latitude and longitude data value (phi, psi) ranges, making the taxi ID number sets recorded in S95 all correspond to a latitude and longitude range, and comparingWhether the same taxi ID numbers are contained in different taxi ID number sets or not at the same time, if so, respectively combining the corresponding longitude and latitude range and the taxi ID number set, so that each taxi ID number after combination is uniquely contained in one taxi ID number set and corresponds to one longitude and latitude range, and simultaneously recording and storing the corresponding taxi ID number set and the corresponding longitude and latitude range to form a vehicle aggregation pre-alarm database DBtaxi,j
S97, if j is less than t (wherein t is a preset value of a positive integer), or if the vehicle aggregation pre-alarm database at a certain moment from max {1, j-t +1} to the jth moment is empty, self-increment j by 1 and return to S91, otherwise, executing S98 to process the vehicle aggregation pre-alarm database at the consecutive t moments from the jth-t +1 to the jth moment;
s98, merging the vehicle gathering pre-alarm database DB on the basis of keeping the corresponding relation between the taxi ID number set and the latitude and longitude rangetaxi,j-t+1,DBtaxi,j-t+2,…,DBtaxi,jThen, whether the same taxi ID number is contained in different taxi ID number sets or not is compared, if the same taxi ID number is contained in different taxi ID number sets, the corresponding longitude and latitude ranges and the corresponding taxi ID number sets are respectively merged, and then u taxi ID number sets S are respectively obtained1,S2,…,SuMaking each combined taxi ID number only be contained in a taxi ID number set and corresponding to a longitude and latitude range, and making statistics on the combined taxi ID number set SeNumber m of taxi vehiclese,1And SeThe number m of taxi vehicles simultaneously contained in the t pre-alarm databasese,2Where e ═ 1,2, …, u;
s99, if e is 1,2, …, u, the ratio of the number of vehicles counted in S98 satisfies me,2/me,1More than or equal to w (wherein w is a proportional preset value and satisfies that w is more than 0 and less than 1), a vehicle aggregation alarm is sent out, and DB is output simultaneouslytaxi,jIs contained in SeThe taxi ID number set and the corresponding latitude and longitude range, then increase j by 1 and return to S91.
The invention relates to taxi vehicle gathering alarm terminal equipment based on road network and site information, which comprises a memory, a processor and a computer program stored in the memory and capable of running on the processor, wherein the processor realizes the steps of a taxi vehicle gathering alarm method based on road network and site information when executing the computer program.
The invention relates to a computer-readable storage medium, which stores a computer program that, when executed by a processor, implements the steps of a taxi vehicle aggregation alarm method based on road network and site information.
The invention has the beneficial effects that:
according to the invention, the longitude and latitude data values are mapped into the datum point data values, so that the fault tolerance of the supervision platform is improved, and the processing efficiency of the platform is ensured. Meanwhile, the information of the road network and the station of the taxi operation area is integrated, and is associated with the historical gathering condition, so that the triggering condition of gathering alarm is strengthened, and the false alarm rate is effectively reduced. And the spherical model is used for carrying out data processing on the longitude and latitude, and the accuracy of aggregation alarm is also improved to a certain extent. Furthermore, a tracking function for aggregation targets in supervision may also be implemented in association with historical aggregation situations.
Drawings
Fig. 1 is a flowchart of a method according to a first embodiment of the present invention.
Detailed Description
To further illustrate the various embodiments, the invention provides the accompanying drawings. The accompanying drawings, which are incorporated in and constitute a part of this disclosure, illustrate embodiments of the invention and, together with the description, serve to explain the principles of the embodiments. Those skilled in the art will appreciate still other possible embodiments and advantages of the present invention with reference to these figures. Elements in the figures are not drawn to scale and like reference numerals are generally used to indicate like elements.
The invention will now be further described with reference to the accompanying drawings and detailed description.
The first embodiment is as follows:
the same latitude difference value on any warp circle corresponds to equal spherical distance, but the spherical distances corresponding to the same longitude difference value on different weft circles are unequal. The prior art basically treats longitude and latitude as plane data for processing. The invention respectively represents longitude values and latitude values of any positioning points on the earth surface by using phi and psi, wherein east longitude is positive, west longitude is negative, north latitude is positive, and south latitude is negative, and then longitude and latitude of the positioning position are used as spherical data for modeling, thereby providing a taxi vehicle gathering alarm method based on road network and station information.
The method specifically comprises the following steps:
s1: longitude and latitude data value of arbitrary positioning point A in taxi operation area
Figure BDA0001425954170000081
According to the longitude and latitude data value of the positioning point A, the longitude and latitude value of a point B1 with the length distance h along the longitude line from the positioning point A is calculated
Figure BDA0001425954170000082
The longitude and latitude value of a point B2 with the length distance h from the positioning point A along the latitude line
Figure BDA0001425954170000083
Wherein h is a fixed distance step, and the calculation formula of (psi, phi) is as follows:
Figure BDA0001425954170000084
where pi is the circumferential ratio, R is the mean radius of the earth, and the trigonometric functions are all calculated in an angular system.
Proceeding to step S2;
s2: selecting longitude and latitude data value (phi) of a positioning point C on the earth surfaceanyany) Usually, it is directly (0 ° ), and then an integer N greater than 1 is determined, and the process proceeds to step S3;
s3: according to phi and psi obtained by calculation in S1, a group of integers p and q are taken on the basis of S2, wherein p is more than or equal to 0 and is less than N, and q is more than or equal to 0 and is less than N, so that
Figure BDA0001425954170000085
Thereby obtaining a plurality of reference points with data values of (phi)0+kφ,Ψ0+ l ψ) where k and l are each satisfied
Figure BDA0001425954170000086
Step (3), the process proceeds to step S4;
s4: for any set of k and l meeting the requirements in S3, the latitude and longitude ranges will be satisfied
Figure BDA0001425954170000087
The latitude and longitude data values (phi, psi) of the positioning point D of (1) are all mapped into a data value (phi) represented by the data value of the reference point through a many-to-one correspondence0+kφ,Ψ0+ l ψ), wherein the many-to-one correspondence formula is
Figure BDA0001425954170000091
Wherein [ ] is the operator of rounding, go to step S5;
s5: according to the urban map data of the operation area, acquiring longitude and latitude data values corresponding to road nodes, carrying out spherical equal-division interpolation between adjacent nodes of each road to enable the distance between adjacent nodes in the interpolated road to be smaller than a preset value r, then storing all the longitude and latitude data values of the interpolated road nodes,
the method comprises the following specific steps:
s51, using formula
Figure BDA0001425954170000092
Calculating road adjacent node (phi)AA) And (phi)BB) Spherical distance D between0Wherein Δ ═ arccos (sin (Ψ)A)sin(ΨB)+cos(ΨA)cos(ΨB)cos(ΦAB));
S52, spherical distance D obtained according to S510Calculating the number of equal parts of the equal-part interpolation
Figure BDA0001425954170000093
Wherein
Figure BDA0001425954170000094
Operator for rounding in infinite direction;
s53, using formula
Figure BDA0001425954170000095
Converting the longitude and latitude data value of the road node spherical coordinate into rectangular coordinate, and then converting the longitude and latitude data value (phi) of the road adjacent nodeAA) And (phi)BB) Can be respectively converted into (x)A,yA,zA) And (x)B,yB,zB);
S54, using formula
Figure BDA0001425954170000096
The road adjacent node (x) obtained in S53 is calculatedA,yA,zA) And (x)B,yB,zB) The linear distance therebetween;
s55, using formula
Figure BDA0001425954170000097
Calculating the spatial weight of spherical equal-division interpolation between adjacent nodes, wherein i is 1,2, …, n-1;
s56, using formula
Figure BDA0001425954170000101
(x) obtained in S54A,yA,zA) And (x)B,yB,zB) All spatial linear interpolation points in between
Figure BDA0001425954170000102
S57, using formula
Figure BDA0001425954170000103
The spatial linear interpolation point obtained in S56
Figure BDA0001425954170000104
Projecting the spherical surface along the center of the sphere to the direction of the spherical surface to obtain the rectangular coordinate (x) of the spherical interpolation pointAB,i,yAB,i,zAB,i) Wherein
Figure BDA0001425954170000105
S58, using formula
Figure BDA0001425954170000106
Wherein R is the mean radius of the earth, R1=Rcos(Ψ)
Converting the rectangular coordinates of the road nodes into longitude and latitude data values of spherical coordinates, then obtaining (x) in S57AB,i,yAB,i,zAB,i) Can be converted into (phi)AB,iAB,i);
And S59, storing all longitude and latitude data values of the road nodes after interpolation.
Proceeding to step S6;
s6: traversing all possible combination values of p and q, and mapping the latitude and longitude data values of the road nodes stored in S5 to corresponding reference point data values for each combination value of p and q according to the correspondence in S4Counting the occurrence times of each datum point data value after all the mappings are finished, recording and storing each datum point data value and the occurrence times thereof, and forming a road network information database DBwayThen, the process proceeds to step S7;
s7: obtaining the position of each taxi station according to the city map data of the operation area, wherein each taxi station comprises but is not limited to: airport, railway station and bus station, and storing the longitude and latitude data value corresponding to the station position, and entering the step S8;
s8: traversing all possible combination values of p and q, respectively mapping the latitude and longitude data values stored in S7 to corresponding reference point data values according to the corresponding relation in S4 aiming at each combination value of p and q, recording and storing the reference point data values after all mappings are finished, and forming a site information database DBsiteThen, the process proceeds to step S9;
s9: according to the road network information database DBwayAnd a site information database DBsiteAnd (4) giving an alarm for the abnormal gathering condition of the taxi vehicles according to the data value.
The method comprises the following specific steps:
s91, acquiring the real-time positioning position of the taxi operated at the jth moment (wherein j is 1,2, …), and storing the corresponding taxi ID number and the longitude and latitude data value of the real-time positioning position;
s92, traversing all possible combination values of k and l, mapping the latitude and longitude data values stored in the S91 to corresponding datum point data values according to the corresponding relation in the S4 aiming at the combination values of k and l, and counting the occurrence times of each datum point data value after all mappings are finished;
s93, if the times of the datum point data values in the S92 are all smaller than the non-negative preset value X, enabling the vehicle at the jth moment to gather the pre-warning database DBtaxi,jIf the position is empty, increasing j by 1 and returning to S91, otherwise, recording the reference point data value with the occurrence frequency larger than X, and mapping the occurrence frequency and the positioning position into a taxi ID number set corresponding to the reference point data value;
s94, for each reference point data value recorded in S93, respectively in the road network information dataLibrary DBwayAnd a site information database DBsiteIf it is in DBwayIf the data value exists in the database, the beta is taken as the number of times of the data value of the corresponding datum point, otherwise the data value is taken as zero, if the data value exists in the database, the data value of the corresponding datum point is taken as the number of times of the data value of the corresponding datum point, and ifsiteIf yes, the value of gamma is 1, otherwise the value of gamma is 0;
s95, if the occurrence times of the datum point data values recorded in S93 are all smaller than the corresponding (X + beta Y + gamma Z) values searched in S94 (wherein Y, Z is a non-negative preset value), making DB have non-negative valuestaxi,jIf the value is null, increasing j by 1 and returning to S91, otherwise, recording the datum point data value with the occurrence frequency greater than the corresponding (X + beta Y + gamma Z) value, and mapping the occurrence frequency and the positioning position of the datum point data value into a taxi ID number set corresponding to the datum point data value;
s96, the reference point data value recorded in S95 is used
Figure BDA0001425954170000121
Converting the datum data values (phi ', psi') into ranges of longitude and latitude data values (phi, psi), enabling taxi ID number sets recorded in S95 to correspond to a longitude and latitude range, comparing whether the same taxi ID numbers are contained in different taxi ID number sets at the same time, respectively combining the corresponding longitude and latitude range and taxi ID number set if the same taxi ID number set is contained, enabling each merged taxi ID number to be uniquely contained in one taxi ID number set and correspond to one longitude and latitude range, and simultaneously recording and storing the corresponding taxi ID number set and the longitude and latitude range to form a vehicle aggregation pre-alarm database DBtaxi,j
S97, if j is less than t (wherein t is a preset value of a positive integer), or if the vehicle aggregation pre-alarm database at a certain moment from max {1, j-t +1} to the jth moment is empty, self-increment j by 1 and return to S91, otherwise, executing S98 to process the vehicle aggregation pre-alarm database at the consecutive t moments from the jth-t +1 to the jth moment;
s98, merging the vehicle gathering pre-alarm database DB on the basis of keeping the corresponding relation between the taxi ID number set and the latitude and longitude rangetaxi,j-t+1,DBtaxi,j-t+2,…,DBtaxi,jThen, whether the same taxi ID number is contained in different taxi ID number sets or not is compared, if the same taxi ID number is contained in different taxi ID number sets, the corresponding longitude and latitude ranges and the corresponding taxi ID number sets are respectively merged, and then u taxi ID number sets S are respectively obtained1,S2,…,SuMaking each combined taxi ID number only be contained in a taxi ID number set and corresponding to a longitude and latitude range, and making statistics on the combined taxi ID number set SeNumber m of taxi vehiclese,1And SeThe number m of taxi vehicles simultaneously contained in the t pre-alarm databasese,2Where e ═ 1,2, …, u;
s99, if e is 1,2, …, u, the ratio of the number of vehicles counted in S98 satisfies me,2/me,1More than or equal to w (wherein w is a proportional preset value and satisfies that w is more than 0 and less than 1), a vehicle aggregation alarm is sent out, and DB is output simultaneouslytaxi,jIs contained in SeThe taxi ID number set and the corresponding latitude and longitude range, then increase j by 1 and return to S91.
Example two:
the invention also provides taxi vehicle gathering alarm terminal equipment based on road network and site information, which comprises a memory, a processor and a computer program stored in the memory and capable of running on the processor, wherein the steps in the method embodiments of the invention are realized when the processor executes the computer program.
Further, as an executable scheme, the taxi vehicle gathering alarm terminal device based on the road network and the station information may be a desktop computer, a notebook, a palm computer, a cloud server and other computing devices. The taxi vehicle gathering alarm terminal device based on the road network and station information can comprise, but is not limited to, a processor and a memory. It is understood by those skilled in the art that the above-mentioned structure of the taxi vehicle aggregation alarm terminal device based on the road network and the station information is only an example of the taxi vehicle aggregation alarm terminal device based on the road network and the station information, and does not constitute a limitation on the taxi vehicle aggregation alarm terminal device based on the road network and the station information, and may include more or less components than the above-mentioned structure, or combine some components, or different components, for example, the taxi vehicle aggregation alarm terminal device based on the road network and the station information may further include an input and output device, a network access device, a bus, etc., which is not limited by the embodiments of the present invention.
Further, as an executable solution, the Processor may be a Central Processing Unit (CPU), other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), an off-the-shelf Programmable Gate Array (FPGA) or other Programmable logic device, a discrete Gate or transistor logic device, a discrete hardware component, and the like. The general processor can be a microprocessor or the processor can be any conventional processor and the like, the processor is a control center of the taxi vehicle aggregation alarm terminal device based on the road network and site information, and various interfaces and lines are used for connecting various parts of the whole taxi vehicle aggregation alarm terminal device based on the road network and site information.
The memory can be used for storing the computer program and/or module, and the processor realizes various functions of the taxi vehicle gathering alarm terminal device based on road network and site information by operating or executing the computer program and/or module stored in the memory and calling data stored in the memory. The memory can mainly comprise a program storage area and a data storage area, wherein the program storage area can store an operating system and an application program required by at least one function; the storage data area may store data created according to the use of the mobile phone, and the like. In addition, the memory may include high speed random access memory, and may also include non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), at least one magnetic disk storage device, a Flash memory device, or other volatile solid state storage device.
The invention also provides a computer-readable storage medium, in which a computer program is stored, which, when being executed by a processor, carries out the steps of the above-mentioned method of an embodiment of the invention.
The integrated module/unit of the taxi vehicle gathering alarm terminal device based on the road network and site information can be stored in a computer readable storage medium if the module/unit is implemented in the form of a software functional unit and sold or used as an independent product. Based on such understanding, all or part of the flow of the method according to the embodiments of the present invention may also be implemented by a computer program, which may be stored in a computer-readable storage medium, and when the computer program is executed by a processor, the steps of the method embodiments may be implemented. Wherein the computer program comprises computer program code, which may be in the form of source code, object code, an executable file or some intermediate form, etc. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, usb disk, removable hard disk, magnetic disk, optical disk, computer Memory, Read-Only Memory (ROM), Random Access Memory (RAM), electrical carrier wave signals, telecommunications signals, software distribution medium, and the like. It should be noted that the computer readable medium may contain content that is subject to appropriate increase or decrease as required by legislation and patent practice in jurisdictions, for example, in some jurisdictions, computer readable media does not include electrical carrier signals and telecommunications signals as is required by legislation and patent practice.
According to the taxi vehicle gathering alarm method based on the road network and site information, the terminal device and the storage medium, the longitude and latitude data values are mapped into the datum point data values, so that the fault tolerance of the supervision platform is improved, and the processing efficiency of the platform is ensured. Meanwhile, the information of the road network and the station of the taxi operation area is integrated, and is associated with the historical gathering condition, so that the triggering condition of gathering alarm is strengthened, and the false alarm rate is effectively reduced. And the spherical model is used for carrying out data processing on the longitude and latitude, and the accuracy of aggregation alarm is also improved to a certain extent. Furthermore, a tracking function for aggregation targets in supervision may also be implemented in association with historical aggregation situations.
While the invention has been particularly shown and described with reference to a preferred embodiment, it will be understood by those skilled in the art that various changes in form and detail may be made therein without departing from the spirit and scope of the invention as defined by the appended claims.

Claims (7)

1. A taxi vehicle gathering alarm method based on road network and station information is characterized in that: the method comprises the following steps:
s1: longitude and latitude data value of arbitrary positioning point A in taxi operation area
Figure FDA0002817076420000011
According to the longitude and latitude data value of the positioning point A, the longitude and latitude value of a point B1 with the length distance h along the longitude line from the positioning point A is calculated
Figure FDA0002817076420000012
The longitude and latitude value of a point B2 with the length distance h from the positioning point A along the latitude line
Figure FDA0002817076420000013
Wherein h is a fixed distance step length, and the step S2 is entered;
s2: selecting longitude and latitude data value (phi) of a positioning point C on the earth surfaceanyany) Then determining an integer N greater than 1, and entering the step S3;
s3: according to phi and psi obtained by calculation in S1, a group of integers p and q are taken on the basis of S2, wherein p is more than or equal to 0 and less than N, q is more than or equal to 0 and less than N, so that
Figure FDA0002817076420000014
Thereby obtaining a plurality of reference points with data values of (phi)0+kφ,Ψ0+ l ψ) where k and l are each satisfied
Figure FDA0002817076420000015
Step (3), the process proceeds to step S4;
s4: for any set of k and l meeting the requirements in S3, the latitude and longitude ranges will be satisfied
Figure FDA0002817076420000016
The latitude and longitude data values (phi, psi) of the positioning point D of (1) are all mapped into a data value (phi) represented by the data value of the reference point through a many-to-one correspondence0+kφ,Ψ0+ l ψ), wherein the many-to-one correspondence formula is
Figure FDA0002817076420000017
Wherein [ ] is the operator of rounding, go to step S5;
s5: according to the urban map data of the operation area, acquiring longitude and latitude data values corresponding to road nodes, performing spherical equal-division interpolation between adjacent nodes of each road to enable the distance between adjacent nodes in the interpolated road to be smaller than a preset value r, storing all the longitude and latitude data values of the interpolated road nodes, and entering the step S6;
s6: traversing all possible combination values of p and q, mapping the longitude and latitude data values of the road nodes stored in the S5 into corresponding reference point data values according to the corresponding relation in the S4 aiming at each combination value of p and q, counting the occurrence times of each reference point data value after all the mappings are finished, recording and storing each reference point data value and the occurrence times thereof, and recording and storing the shape and the shape of each reference point data valueInformation database DB of adult road networkwayThen, the process proceeds to step S7;
s7: obtaining the position of each taxi station according to the city map data of the operation area, storing the longitude and latitude data value corresponding to the station position, and entering the step S8;
s8: traversing all possible combination values of p and q, respectively mapping the latitude and longitude data values stored in S7 to corresponding reference point data values according to the corresponding relation in S4 aiming at each combination value of p and q, recording and storing the reference point data values after all mappings are finished, and forming a site information database DBsiteThen, the process proceeds to step S9;
s9: according to the road network information database DBwayAnd a site information database DBsiteAnd (4) giving an alarm for the abnormal gathering condition of the taxi vehicles according to the data value.
2. The taxi vehicle gathering alarm method based on road network and site information as claimed in claim 1, wherein: in S1, the formula for calculating (ψ, Φ) is:
Figure FDA0002817076420000021
where π is the circumference ratio and R is the earth's mean radius.
3. The taxi vehicle gathering alarm method based on road network and site information as claimed in claim 1, wherein: in S5, performing equal spherical interpolation between adjacent nodes of each road to make the distance between adjacent nodes inside the interpolated road smaller than a preset value r, and then storing all longitude and latitude data values of the interpolated road nodes, specifically:
s51, using formula
Figure FDA0002817076420000031
Calculating road adjacent node (phi)AA) And (phi)BB) Spherical distance D between0Wherein Δ ═arccos(sin(ΨA)sin(ΨB)+cos(ΨA)cos(ΨB)cos(ΦAB));
S52, spherical distance D obtained according to S510Calculating the number of equal parts of the equal-part interpolation
Figure FDA0002817076420000032
Wherein
Figure FDA0002817076420000033
Operator for rounding in infinite direction;
s53, using formula
Figure FDA0002817076420000034
Converting the longitude and latitude data value of the road node spherical coordinate into rectangular coordinate, and then converting the longitude and latitude data value (phi) of the road adjacent nodeAA) And (phi)BB) Can be respectively converted into (x)A,yA,zA) And (x)B,yB,zB);
S54, using formula
Figure FDA0002817076420000035
The road adjacent node (x) obtained in S53 is calculatedA,yA,zA) And (x)B,yB,zB) The linear distance therebetween;
s55, using formula
Figure FDA0002817076420000036
Calculating the spatial weight of spherical equal-division interpolation between adjacent nodes, wherein i is 1,2, …, n-1;
s56, using formula
Figure FDA0002817076420000037
(x) obtained in S54A,yA,zA) And (x)B,yB,zB) All spatial linear interpolation points in between
Figure FDA0002817076420000038
S57, using formula
Figure FDA0002817076420000039
The spatial linear interpolation point obtained in S56
Figure FDA0002817076420000041
Projecting the spherical surface along the center of the sphere to the direction of the spherical surface to obtain the rectangular coordinate (x) of the spherical interpolation pointAB,i,yAB,i,zAB,i) Wherein
Figure FDA0002817076420000042
S58, using formula
Figure FDA0002817076420000043
Wherein R is the mean radius of the earth, R1=Rcos(Ψ)
Converting the rectangular coordinates of the road nodes into longitude and latitude data values of spherical coordinates, then obtaining (x) in S57AB,i,yAB,i,zAB,i) Can be converted into (phi)AB,iAB,i);
And S59, storing all longitude and latitude data values of the road nodes after interpolation.
4. The taxi vehicle gathering alarm method based on road network and site information as claimed in claim 1, wherein: in S7, each taxi station includes, but is not limited to: airports, train stations and bus stations.
5. The taxi vehicle gathering alarm method based on road network and site information as claimed in claim 1, wherein: in S9, according to road network information database DBwayAnd a site information database DBsiteThe data value of (2) gives an alarm to abnormal gathering conditions of taxi vehicles, and specifically comprises the following steps:
s91, acquiring a real-time positioning position of the taxi operated at the jth moment, and storing a corresponding taxi ID number and a longitude and latitude data value of the real-time positioning position, wherein j is 1,2, …;
s92, traversing all possible combination values of p and q, mapping the latitude and longitude data values stored in S91 to corresponding datum point data values according to the corresponding relation in S4 aiming at each combination value of p and q, and counting the occurrence times of each datum point data value after all mappings are finished;
s93, if the times of the datum point data values in the S92 are all smaller than the non-negative preset value X, enabling the vehicle at the jth moment to gather the pre-warning database DBtaxi,jIf the position is empty, increasing j by 1 and returning to S91, otherwise, recording the reference point data value with the occurrence frequency larger than X, and mapping the occurrence frequency and the positioning position into a taxi ID number set corresponding to the reference point data value;
s94, for each reference point data value recorded in S93, respectively in the road network information database DBwayAnd a site information database DBsiteIf it is in DBwayIf the data value exists in the database, the beta is taken as the number of times of the data value of the corresponding datum point, otherwise the data value is taken as zero, if the data value exists in the database, the data value of the corresponding datum point is taken as the number of times of the data value of the corresponding datum point, and ifsiteIf yes, the value of gamma is 1, otherwise the value of gamma is 0;
s95, if the occurrence times of the datum point data values recorded in S93 are all smaller than the corresponding X + beta Y + gamma Z values searched in S94, making DBtaxi,jIf the value is empty, increasing j by 1 and returning to S91, otherwise, recording a reference point data value with the occurrence frequency larger than the corresponding X + beta Y + gamma Z value, and mapping the occurrence frequency and the positioning position of the reference point data value into a taxi ID number set corresponding to the reference point data value, wherein Y, Z is a non-negative preset value;
s96, the reference point data value recorded in S95 is used
Figure FDA0002817076420000051
Converting the datum data values (phi ', psi') into ranges of longitude and latitude data values (phi, psi), enabling taxi ID number sets recorded in S95 to correspond to a longitude and latitude range, comparing whether the same taxi ID numbers are contained in different taxi ID number sets at the same time, respectively combining the corresponding longitude and latitude range and taxi ID number set if the same taxi ID number set is contained, enabling each merged taxi ID number to be uniquely contained in one taxi ID number set and correspond to one longitude and latitude range, and simultaneously recording and storing the corresponding taxi ID number set and the longitude and latitude range to form a vehicle aggregation pre-alarm database DBtaxi,j
S97, if j < t or the vehicle aggregation pre-alarm database at a certain moment from max {1, j-t +1} to j is empty, enabling j to be increased by 1 and returning to S91, otherwise executing S98 to process the vehicle aggregation pre-alarm database at t continuous moments from j-t +1 to j, wherein t is a preset value of a positive integer;
s98, merging the vehicle gathering pre-alarm database DB on the basis of keeping the corresponding relation between the taxi ID number set and the latitude and longitude rangetaxi,j-t+1,DBtaxi,j-t+2,…,DBtaxi,jThen, whether the same taxi ID number is contained in different taxi ID number sets or not is compared, if the same taxi ID number is contained in different taxi ID number sets, the corresponding longitude and latitude ranges and the corresponding taxi ID number sets are respectively merged, and then u taxi ID number sets S are respectively obtained1,S2,…,SuMaking each combined taxi ID number only be contained in a taxi ID number set and corresponding to a longitude and latitude range, and making statistics on the combined taxi ID number set SeNumber m of taxi vehiclese,1And SeThe number m of taxi vehicles simultaneously contained in the t pre-alarm databasese,2Where e ═ 1,2, …, u;
s99, if e is 1,2, …, u, the ratio of the number of vehicles counted in S98 satisfies me,2/me,1If not less than w, the vehicle aggregation alarm is sent out, and DB is output simultaneouslytaxi,jIs contained in SeThe taxi ID number set and the corresponding latitude and longitude range thereof, then, j is increased by 1 and returns to S91, wherein w is a proportion preset value and satisfies 0<w<1。
6. A taxi vehicle aggregation alarm terminal device based on road network and site information, comprising a memory, a processor and a computer program stored in the memory and operable on the processor, characterized in that: the processor, when executing the computer program, realizes the steps of the method according to claims 1-5.
7. A computer-readable storage medium storing a computer program, characterized in that: the computer program realizing the steps of the method as claimed in claims 1-5 when executed by a processor.
CN201710916836.XA 2017-09-30 2017-09-30 Taxi gathering alarm method, terminal equipment and storage medium Active CN109598925B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710916836.XA CN109598925B (en) 2017-09-30 2017-09-30 Taxi gathering alarm method, terminal equipment and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710916836.XA CN109598925B (en) 2017-09-30 2017-09-30 Taxi gathering alarm method, terminal equipment and storage medium

Publications (2)

Publication Number Publication Date
CN109598925A CN109598925A (en) 2019-04-09
CN109598925B true CN109598925B (en) 2021-03-02

Family

ID=65955831

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710916836.XA Active CN109598925B (en) 2017-09-30 2017-09-30 Taxi gathering alarm method, terminal equipment and storage medium

Country Status (1)

Country Link
CN (1) CN109598925B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110197291A (en) * 2019-04-09 2019-09-03 安徽四创电子股份有限公司 Net based on Beidou about vehicle information exchange and supervisory systems
CN114863670B (en) * 2022-01-17 2022-12-02 深圳市慧视通科技股份有限公司 Method and device for early warning of taxi gathering condition

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090051510A1 (en) * 2007-08-21 2009-02-26 Todd Follmer System and Method for Detecting and Reporting Vehicle Damage
CN102117532A (en) * 2009-12-31 2011-07-06 厦门雅迅网络股份有限公司 Method for pre-alarming illegal gathering of taxis based on GPS (global positioning system)
CN104021668A (en) * 2014-06-26 2014-09-03 中国科学院自动化研究所 Public traffic supply-demand detection and prediction system and method
CN104715608A (en) * 2015-03-26 2015-06-17 杭州电子科技大学 Around-the-clock all-area taxi gathering real-time monitoring method based on HBase
US20170243485A1 (en) * 2012-04-24 2017-08-24 Zetta Research and Development LLC, ForC series V2v safety system using learned signal timing

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20150072471A (en) * 2013-12-19 2015-06-30 건국대학교 산학협력단 Traffic flow prediction system using spatiotemporal stochastic model

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090051510A1 (en) * 2007-08-21 2009-02-26 Todd Follmer System and Method for Detecting and Reporting Vehicle Damage
CN102117532A (en) * 2009-12-31 2011-07-06 厦门雅迅网络股份有限公司 Method for pre-alarming illegal gathering of taxis based on GPS (global positioning system)
US20170243485A1 (en) * 2012-04-24 2017-08-24 Zetta Research and Development LLC, ForC series V2v safety system using learned signal timing
CN104021668A (en) * 2014-06-26 2014-09-03 中国科学院自动化研究所 Public traffic supply-demand detection and prediction system and method
CN104715608A (en) * 2015-03-26 2015-06-17 杭州电子科技大学 Around-the-clock all-area taxi gathering real-time monitoring method based on HBase

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于动态道路映射的基站定位方法;陆萍等;《计算机与数字工程》;20130830(第8期);全文 *

Also Published As

Publication number Publication date
CN109598925A (en) 2019-04-09

Similar Documents

Publication Publication Date Title
CN107657637B (en) Method for acquiring operation area of agricultural machine
WO2019205020A1 (en) Road condition recognition method, apparatus and device
CN113421432B (en) Traffic restriction information detection method and device, electronic equipment and storage medium
CN105551082A (en) Method and device of pavement identification on the basis of laser-point cloud
CN110111599B (en) Parking guidance method based on big data, terminal equipment and storage medium
CN104281709A (en) Method and system for generating traffic information tiled map
CN107993001A (en) A kind of method for visualizing of risk assessment, device and storage medium
CN109598925B (en) Taxi gathering alarm method, terminal equipment and storage medium
CN105489011B (en) Extensive vehicle road network real-time speed limit computational methods based on grid cutting
US10295650B2 (en) Method and device for determining vehicle site location
CN110008921B (en) Road boundary generation method and device, electronic equipment and storage medium
CN113050660B (en) Error compensation method, error compensation device, computer equipment and storage medium
CN111220147B (en) Positioning correction method and device for terminal
CN110853412B (en) Method and device for identifying abnormal track point
CN104850657A (en) Holographic position map superposing method
CN113160589A (en) Vehicle early warning method and device, electronic equipment and computer readable storage medium
CN111710187A (en) Method, device, medium and electronic equipment for determining vehicle driving risk
CN103675763B (en) A kind of based on hardware-accelerated vector quantization VTS radar video rendering method
CN112304281A (en) Road slope measuring method, terminal equipment and storage medium
CN108810036B (en) Method, device and system for processing geographical position information
CN111223293A (en) System and method for analyzing traffic congestion
US11823572B2 (en) Method, electronic device, and system for predicting future overspeeding
CN110807954B (en) Method and device for identifying abnormal track point
CN111767295B (en) Map data processing method, device, computing equipment and medium
WO2022081083A1 (en) Method, electronic device, and system for detecting overspeeding

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