CN113473398A - Mobile phone signaling data stop point identification method, device and storage medium - Google Patents

Mobile phone signaling data stop point identification method, device and storage medium Download PDF

Info

Publication number
CN113473398A
CN113473398A CN202110613302.6A CN202110613302A CN113473398A CN 113473398 A CN113473398 A CN 113473398A CN 202110613302 A CN202110613302 A CN 202110613302A CN 113473398 A CN113473398 A CN 113473398A
Authority
CN
China
Prior art keywords
base station
time
distance
longitude
user
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.)
Granted
Application number
CN202110613302.6A
Other languages
Chinese (zh)
Other versions
CN113473398B (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.)
Sun Yat Sen University
Original Assignee
Sun Yat Sen 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 Sun Yat Sen University filed Critical Sun Yat Sen University
Priority to CN202110613302.6A priority Critical patent/CN113473398B/en
Publication of CN113473398A publication Critical patent/CN113473398A/en
Application granted granted Critical
Publication of CN113473398B publication Critical patent/CN113473398B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W4/00Services specially adapted for wireless communication networks; Facilities therefor
    • H04W4/20Services signaling; Auxiliary data signalling, i.e. transmitting data via a non-traffic channel
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

The invention discloses a method, a device and a storage medium for identifying a mobile phone signaling data stop point; the method comprises the steps of obtaining mobile phone signaling data; sequencing the signaling data of the same user mobile phone according to the generation time of the signaling record to obtain a sequencing list; calculating the continuous residence time of the user in the first base station according to the sorted list; if the continuous stay time is longer than the first threshold time, determining the first base station as a stay point; and if the continuous residence time is not greater than the first threshold value, acquiring a first distance group and a first time sum, and determining a residence point according to the first distance group and the first time sum. The method can accurately identify the user staying points, so that the method is suitable for extracting the trip chain of the mobile phone user, and mobile phone signaling data can be better served for the research in the fields of urban traffic planning, population distribution, privacy protection and the like; the invention can be widely applied to the technical field of mobile phone signaling data processing.

Description

Mobile phone signaling data stop point identification method, device and storage medium
Technical Field
The invention relates to the technical field of mobile phone signaling data processing, in particular to a mobile phone signaling data stop point identification method, a mobile phone signaling data stop point identification device and a storage medium.
Background
The mobile phone signaling data is a communication record captured by a communication base station of a mobile phone operator when a mobile phone user makes a call, sends a short message or moves a position to cause base station switching. The mobile phone signaling data is communication record data between the mobile phone and the base station. Once the mobile phone is started, the signaling data starts to be generated. After desensitization processing, the mobile phone signaling data can be used for the research on the aspects of resident behavior preference, resident movement track analysis, city planning, privacy protection and the like. The mobile phone signaling data is a novel big data resource, and has the advantages that the mobile phone signaling data has real-time performance, integrity, space-time full coverage and other data resources, and the like, compared with other types of data, the mobile phone signaling data is a novel big data resource, and has the characteristics of large sample size, long observation period, short sampling period, strong trackability and the like, so that the mobile phone signaling data has the advantages that other data resources do not have, and the mobile phone signaling data has the advantages of real-time performance, integrity, space-time full coverage and the like, and has attracted the wide attention of researchers.
The application research of the mobile phone signaling data is based on the identification of a stop point, and a user travel chain can be constructed only by identifying the stop position of a user through the mobile phone signaling data, so that the research on the aspects of resident behavior preference, resident movement track analysis, city planning, privacy protection and the like can be carried out. Due to the generation mechanism of the mobile phone signaling data, when a user does not use the mobile phone, the mobile phone can communicate with a nearby base station irregularly to generate signaling records, so that the mobile phone signaling data has partial redundancy. Meanwhile, the mobile phone may be switched back and forth with several nearby base stations due to the slow movement of the mobile phone user in one place, and the mode is considered as a potential stay mode of the user; resulting in certain difficulties in identifying the location where the user is staying.
Disclosure of Invention
The present invention is directed to solving at least one of the problems of the prior art. Therefore, the invention provides a method, a device and a storage medium for processing mobile phone signaling data.
The technical scheme adopted by the invention is as follows:
in one aspect, an embodiment of the present invention includes a method for identifying a mobile phone signaling data dwell point, including:
acquiring mobile phone signaling data, wherein the mobile phone signaling data comprises a signaling record number, a signaling user number, a base station longitude, a base station latitude and signaling record generation time;
sequencing the signaling data of the same user mobile phone according to the generation time of the signaling record to obtain a sequencing list;
calculating the continuous stay time of the user in a first base station according to the sorted list, wherein the first base station is any base station where the user stays;
if the continuous stay time is longer than a first threshold time, determining the first base station as a stay point;
if the continuous staying time is not more than the first threshold time, acquiring a first distance group and a first time sum, wherein the first distance group comprises the distance between the first base station and a first target base station and the distance between a second base station and the first target base station, the first time sum is the sum of the staying time of a user at the first base station and the staying time of the user at the second base station, the second base station is a base station which is not in addition to the first base station and has the continuous staying time not more than the first threshold time, and the first target base station is a possible base station calculated according to the first base station and the second base station;
and determining a stopping point according to the first distance group and the first time sum.
Further, the step of calculating the continuous staying time of the user at the first base station according to the sorted list includes:
acquiring a first time and a second time according to the sorted list, wherein the first time is the earliest signaling record generation time corresponding to the first base station, and the second time is the latest signaling record generation time corresponding to the first base station;
and calculating the continuous staying time of the user in the first base station according to the first time and the second time.
Further, the first distance group is obtained by the following steps:
acquiring a first information group according to the sorted list, wherein the first information group comprises longitude and latitude coordinates of the first base station, longitude and latitude coordinates of the second base station, earliest signaling record generation time corresponding to the first base station, latest signaling record generation time corresponding to the first base station, earliest signaling record generation time corresponding to the second base station and latest signaling record generation time corresponding to the second base station;
according to the first information group, calculating to obtain longitude and latitude coordinates of the first target base station;
and calculating to obtain the first distance group according to the longitude and latitude coordinates of the first base station, the longitude and latitude coordinates of the second base station and the longitude and latitude coordinates of the first target base station.
Further, the first time sum is obtained by:
calculating the staying time of the user in the first base station according to the earliest signaling record generation time corresponding to the first base station and the latest signaling record generation time corresponding to the first base station;
calculating the staying time of the user in the second base station according to the earliest signaling record generation time corresponding to the second base station and the latest signaling record generation time corresponding to the second base station;
and adding the stay time of the user in the first base station and the stay time of the user in the second base station to obtain a first time sum.
Further, the longitude and latitude coordinates of the first target base station are calculated by a first formula, where the first formula is:
Figure BDA0003096482760000021
in the formula, longnewIs the longitude, lat, of the first target base stationnewIs the latitude, t, of the first target base station1Generating a time, t, for the earliest signalling record of a base station inGenerating a time, t, for the latest signalling record of base station iiFor the stay time of the user at base station i, lngiIs the longitude, lat, of base station iiIs the latitude of base station i, i denotes the base station.
Further, the step of determining a stopping point according to the first distance group and the first time sum comprises:
if the distance between the first base station and the first target base station and the distance between the second base station and the first target base station are both smaller than a first threshold distance, extracting a third base station, wherein the third base station is a base station except the first base station and the second base station;
acquiring a second information group according to the ordered list, wherein the second information group comprises the first information group, longitude and latitude coordinates of the third base station, earliest signaling record generation time corresponding to the third base station and latest signaling record generation time corresponding to the third base station;
according to the second information group, calculating to obtain the longitude and latitude coordinates of a second target base station through the first formula;
calculating to obtain a second distance group according to the longitude and latitude coordinates of the first base station, the longitude and latitude coordinates of the second base station, the longitude and latitude coordinates of the third base station and the longitude and latitude coordinates of the second target base station, wherein the second distance group comprises the distance from the first base station to the second target base station, the distance from the second base station to the second target base station and the distance from the third base station to the second target base station;
if at least one distance in the second distance group is greater than the first threshold distance, judging whether the first time sum is greater than the first threshold time;
and if the first time sum is greater than the first threshold time, determining that the first target base station is a stop point.
Further, if all the distances in the second distance group are smaller than the first threshold distance, extracting a fourth base station, where the fourth base station is a base station other than the first base station, the second base station, and the third base station;
calculating to obtain a second time sum according to the second information group, wherein the second time sum is the sum of the staying time of the user in the first base station, the second base station and the third base station;
acquiring a third information group according to the ordered list, wherein the second information group comprises the second information group, longitude and latitude coordinates of the fourth base station, earliest signaling record generation time corresponding to the fourth base station and latest signaling record generation time corresponding to the fourth base station;
according to the third information group, calculating to obtain longitude and latitude coordinates of a third target base station through the first formula;
calculating a third distance group according to the longitude and latitude coordinates of the first base station, the longitude and latitude coordinates of the second base station, the longitude and latitude coordinates of the third base station, the longitude and latitude coordinates of the fourth base station and the longitude and latitude coordinates of the third target base station, wherein the third distance group comprises the distance from the first base station to the third target base station, the distance from the second base station to the third target base station, the distance from the third base station to the third target base station and the distance from the fourth base station to the third target base station;
if at least one distance in the third distance group is greater than the first threshold distance, judging whether the second time sum is greater than the first threshold time;
and if the second time sum is greater than the first threshold time, determining that the second target base station is a stop point.
Further, the first base station is marked as a moving point and recalculated from the second base station to obtain a stopping point if at least one of the following conditions exists:
at least one distance in the first distance group is greater than the first threshold distance;
alternatively, there is at least one distance in the second distance group that is greater than the first threshold distance but the first time sum is not greater than the first threshold time.
The invention has the beneficial effects that:
the method comprises the steps of obtaining mobile phone signaling data, and sequencing the same user mobile phone signaling data according to the signaling record generation time to obtain a sequence list; calculating the continuous stay time of the user in the first base station according to the sorted list, and determining the first base station as a stay point if the continuous stay time is greater than a first threshold time; if the continuous staying time is not more than the first threshold value, acquiring a first distance group and a first time sum, and determining a staying point according to the first distance group and the first time sum; therefore, the user staying points are accurately identified, the method is suitable for extracting the trip chain of the mobile phone user, and the mobile phone signaling data can be better used for the research in the fields of urban traffic planning, population distribution, privacy protection and the like.
Additional aspects and advantages of the invention will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the invention.
Drawings
The above and/or additional aspects and advantages of the present invention will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings of which:
fig. 1 is a flowchart illustrating steps of a method for identifying a mobile phone signaling data stop point according to an embodiment of the present invention;
fig. 2 is a schematic structural diagram of a device for identifying a mobile phone signaling data stop point according to an embodiment of the present invention.
Detailed Description
Reference will now be made in detail to embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to the same or similar elements or elements having the same or similar function throughout. The embodiments described below with reference to the accompanying drawings are illustrative only for the purpose of explaining the present invention, and are not to be construed as limiting the present invention.
In the description of the present invention, it should be understood that the orientation or positional relationship referred to in the description of the orientation, such as the upper, lower, front, rear, left, right, etc., is based on the orientation or positional relationship shown in the drawings, and is only for convenience of description and simplification of description, and does not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and thus, should not be construed as limiting the present invention.
In the description of the present invention, the meaning of a plurality of means is one or more, the meaning of a plurality of means is two or more, and larger, smaller, larger, etc. are understood as excluding the number, and larger, smaller, inner, etc. are understood as including the number. If the first and second are described for the purpose of distinguishing technical features, they are not to be understood as indicating or implying relative importance or implicitly indicating the number of technical features indicated or implicitly indicating the precedence of the technical features indicated.
In the description of the present invention, unless otherwise explicitly limited, terms such as arrangement, installation, connection and the like should be understood in a broad sense, and those skilled in the art can reasonably determine the specific meanings of the above terms in the present invention in combination with the specific contents of the technical solutions.
The embodiments of the present application will be further explained with reference to the drawings.
Referring to fig. 1, an embodiment of the present invention provides a method for identifying a mobile phone signaling data dwell point, including but not limited to the following steps:
s100, acquiring mobile phone signaling data, wherein the mobile phone signaling data comprise a signaling record number, a signaling user number, a base station longitude, a base station latitude and signaling record generation time;
s200, sequencing the mobile phone signaling data of the same user according to the generation time of the signaling record to obtain a sequencing list;
s300, calculating the continuous stay time of the user in a first base station according to the sorted list, wherein the first base station is any base station where the user stays;
s400, if the continuous stay time is longer than a first threshold time, determining that the first base station is a stay point;
s500, if the continuous staying time is not more than a first threshold value, obtaining a first distance group and a first time sum, wherein the first distance group comprises the distance between a first base station and a first target base station and the distance between a second base station and the first target base station, the first time sum is the sum of the staying time of a user at the first base station and the staying time of the user at the second base station, the second base station is a base station which is except the first base station and the continuous staying time of which is not more than the first threshold value time, and the first target base station is a possible base station obtained through calculation according to the first base station and the second base station;
s600, determining a stopping point according to the first distance group and the first time sum.
Regarding step S100, in this embodiment, the signaling record number is denoted as id, which represents a unique identifier of each signaling record; the number of the signaling user is marked as isdn which represents the unique identification of the signaling user; the longitude of the base station is recorded as lng, which represents the longitude of the base station connected by the user; the latitude of the base station is recorded as lat, which represents the latitude of the user connected with the base station; the generation time of the signaling record is recorded as date and represents the time stamp of the generation of the signaling record; recording the gender of the user as the gender, and representing the gender of the mobile phone user; the user age is recorded as age and represents the age of the mobile phone user.
In this embodiment, the mobile phone signaling data of the same user is obtained first, and then sorted according to the signaling record generation time, where the mobile phone signaling data of the same user can be extracted through the signaling user number isdn, and then sorted according to the signaling record generation time for the mobile phone signaling data of the same user obtained by extraction, in this embodiment, referring to table 1, in order to facilitate calculation of the residence time of the user in each base station, the mobile phone signaling data of the same user can be sorted according to the sequence of the signaling record generation time to obtain a sorted list, and different base stations can be named as O according to the time sequence1,O2,…,On
TABLE 1 Signaling data schematic table for the same user's handset
Figure BDA0003096482760000061
Regarding step S300, in this embodiment, after sequencing is performed according to the sequence of the generation time of the signaling record and a sequence table is obtained, step S300 may be executed, where the step specifically includes:
s301, acquiring a first time and a second time according to the sorted list, wherein the first time is the earliest signaling record generation time corresponding to the first base station, and the second time is the latest signaling record generation time corresponding to the first base station;
s302, calculating and obtaining the continuous staying time of the user in the first base station according to the first time and the second time.
Specifically, referring to table 1, the user is at base station O1(113.2603472, 20.1489583) the earliest signaling record generation time is 2021-05-22-12:14:39, namely the first time is 2021-05-22-12:14:39, and the user is at the base station O1(113.2603472,20.1489583) The latest signaling record generation time is 2021-05-22-12:46:39, namely the second time 2021-05-22-12:46:39, the first time is subtracted from the second time, and the user can be obtained at the base station O1(113.2603472, 20.1489583) has a residence time of 32 minutes; in the same way, the user in the base station O in Table 1 can be obtained2(113.2803274, 20.1468958) has a residence time of 1 minute; user at base station O3(113.2813473, 20.1468959) has a residence time of 1 minute 20 seconds; user at base station O4(113.2913473, 20.1468969) has a residence time of 25 minutes; user at base station O5The residence time of (113.2913474, 20.1568979) was 12 minutes 20 seconds.
Regarding step S400, in this embodiment, if the continuous staying time is greater than the first threshold time, it is determined that the first base station is the staying point; specifically, for example, if the first threshold time is set to 30 minutes, in table 1, the user is at the base station O1(113.2603472, 20.1489583) is 32 minutes, then base station O can be determined1The point is left for the user.
Regarding step S500, in this embodiment, if the continuous staying time is not greater than the first threshold, a first distance group and a first time sum are obtained, where the first distance group includes a distance from the first base station to the first target base station and a distance from the second base station to the first target base station, the first time sum is a sum of staying times of the user at the first base station and the second base station, the second base station is a base station that is different from the first base station and has a continuous staying time that is not greater than the first threshold time, and the first target base station is a possible base station calculated according to the first base station and the second base station; the first distance group is obtained through the following steps:
p1, acquiring a first information group according to the sorted list, wherein the first information group comprises longitude and latitude coordinates of a first base station, longitude and latitude coordinates of a second base station, earliest signaling record generation time corresponding to the first base station, latest signaling record generation time corresponding to the first base station, earliest signaling record generation time corresponding to the second base station and latest signaling record generation time corresponding to the second base station;
p2, calculating to obtain longitude and latitude coordinates of the first target base station according to the first information group;
and P3, calculating to obtain a first distance group according to the longitude and latitude coordinates of the first base station, the longitude and latitude coordinates of the second base station and the longitude and latitude coordinates of the first target base station.
Specifically, in table 1, since the user is at base station O2(113.2803274, 20.1468958) has a dwell time of 1 minute, i.e., 30 minutes less than the first threshold time, and thus, base station O can be obtained from Table 12Has longitude and latitude coordinates of (113.2803274, 20.1468958), and is base station O3Has longitude and latitude coordinates of (113.2813473, 20.1468959), and is base station O2The corresponding earliest signaling record generation time is 2021-05-22-12:37:15, and the base station O2The corresponding latest signaling record generation time is 2021-05-22-12:38:15, base station O3The corresponding earliest signaling record generation time is 2021-05-22-12:39:39, and the base station O3The corresponding latest signaling record generation time is 2021-05-22-12:40: 59; according to the information, the longitude and latitude coordinates of the first target base station can be further calculated and obtained, and then the base station O is used2Latitude and longitude coordinates, base station O3And calculating the longitude and latitude coordinates of the first target base station to obtain a first distance group.
In this embodiment, the longitude and latitude coordinates of the first target base station are calculated by a first formula, where the first formula is:
Figure BDA0003096482760000071
in the formula, longnewIs the longitude, lat, of the first target base stationnewIs the latitude, t, of the first target base station1Generating a time, t, for the earliest signalling record of a base station inGenerating a time, t, for the latest signalling record of base station iiFor the stay time of the user at base station i, lngiIs the longitude, lat, of base station iiIs the latitude of base station i, i denotes the base station.
In particular, base station O for example2And base station O3The corresponding information is substituted into the first formula to obtain the longitude and latitude of the first target base stationLabeled (113.2809012, 20.1468959); and then further according to base station O2Longitude and latitude coordinates (113.2803274, 20.1468958), base station O3And the longitude and latitude coordinates (113.2813473, 20.1468959) of the first target base station and the longitude and latitude coordinates (113.2809012, 20.1468959) of the first target base station are calculated to obtain O2The distance to the first target base station is 59 meters, O3The distance to the first target base station is 46 meters.
In this embodiment, the first time sum is obtained by the following steps:
D1. calculating the stay time of the user in the first base station according to the earliest signaling record generation time corresponding to the first base station and the latest signaling record generation time corresponding to the first base station;
D2. calculating the stay time of the user in the second base station according to the earliest signaling record generation time corresponding to the second base station and the latest signaling record generation time corresponding to the second base station;
D3. and adding the stay time of the user in the first base station and the stay time of the user in the second base station to obtain a first time sum.
Specifically, in this embodiment, the first time sum is the sum of the stay times of the user at the first base station and the second base station, and referring to table 1, since the user is at the base station O2(113.2803274, 20.1468958) has a residence time of 1 minute; user at base station O3The residence time of (113.2813473, 20.1468959) was 1 minute 20 seconds, and it can be seen that the first time sum was 2 minutes 20 seconds.
In this embodiment, the step S600, namely, the step of determining the stopping point according to the first distance group and the first time sum, includes:
s601, if the distance between the first base station and the first target base station and the distance between the second base station and the first target base station are both smaller than a first threshold distance, extracting a third base station, wherein the third base station is a base station except the first base station and the second base station;
s602, acquiring a second information group according to the sorted list, wherein the second information group comprises the first information group, the longitude and latitude coordinates of the third base station, the earliest signaling record generation time corresponding to the third base station and the latest signaling record generation time corresponding to the third base station;
s603, according to the second information group, calculating to obtain longitude and latitude coordinates of the second target base station through a first formula;
s604, calculating to obtain a second distance group according to the longitude and latitude coordinates of the first base station, the longitude and latitude coordinates of the second base station, the longitude and latitude coordinates of the third base station and the longitude and latitude coordinates of the second target base station, wherein the second distance group comprises the distance from the first base station to the second target base station, the distance from the second base station to the second target base station and the distance from the third base station to the second target base station;
s605, if at least one distance in the second distance group is larger than a first threshold distance, judging whether the first time sum is larger than the first threshold time;
and S606, if the first time sum is larger than the first threshold time, determining the first target base station as a stop point.
In this embodiment, if all the distances in the second distance group are smaller than the first threshold distance, a fourth base station is extracted, where the fourth base station is a base station other than the first base station, the second base station, and the third base station;
according to the second information group, calculating to obtain a second time sum, wherein the second time sum is the sum of the residence time of the user in the first base station, the second base station and the third base station;
acquiring a third information group according to the sorted list, wherein the second information group comprises the second information group, the longitude and latitude coordinates of a fourth base station, earliest signaling record generation time corresponding to the fourth base station and latest signaling record generation time corresponding to the fourth base station;
according to the third information group, calculating to obtain longitude and latitude coordinates of a third target base station through a first formula;
calculating to obtain a third distance group according to the longitude and latitude coordinates of the first base station, the longitude and latitude coordinates of the second base station, the longitude and latitude coordinates of the third base station, the longitude and latitude coordinates of the fourth base station and the longitude and latitude coordinates of the third target base station, wherein the third distance group comprises the distance from the first base station to the third target base station, the distance from the second base station to the third target base station, the distance from the third base station to the third target base station and the distance from the fourth base station to the third target base station;
if at least one distance in the third distance group is greater than the first threshold distance, judging whether the second time sum is greater than the first threshold time;
and if the second time sum is larger than the first threshold time, determining the second target base station as a stop point.
In this embodiment, if at least one of the following conditions exists, the first base station is marked as a moving point, and recalculation is started from the second base station to obtain a staying point:
at least one distance in the first distance group is greater than a first threshold distance;
alternatively, there is at least one distance in the second distance group that is greater than the first threshold distance but not greater than the first threshold time.
Specifically, referring to table 1, the user is at base station O2(113.2803274, 20.1468958) has a dwell time of 1 minute with the user at base station O3(113.2813473, 20.1468959) is 1 minute and 20 seconds, and the first target base station is calculated to be O according to the first formulanew(113.2809012, 20.1468959), then O2To OnewAt a distance of 59 m, O3To OnewIs 46 meters, and the first threshold distance is set to 1000 meters, since base station O2To Onew1Distance and base station O3To Onew1If the distances are all smaller than the first threshold distance, the next base station needs to be added for calculation, that is, the base station O is extracted from the table 14The information of (a); in this embodiment, if the base station O is obtained by calculation2To Onew1Distance and base station O3To Onew1If any one of the distances is greater than the first threshold distance, it indicates that the base station O is a base station O2And base station O3Does not form a cluster of stop points, at this time, the base station O is set2Marked as mobile point, slave base station O3Starting to recalculate; from table 1, it can be seen that the user is at base station O4(113.2913473, 20.1568969) a residence time of 25 minutes; by a first formula
Figure BDA0003096482760000091
Recalculating to obtain the second target base station as Onew2(113.290445, 20.156042), then base station O2To Onew2Is 1466 meters, base station O3To Onew2Distance of 1391 m, base station O4To Onew2Is 134 m, at this time, due to the presence of the base station O2And base station O3To Onew2Is greater than a first distance threshold, thus will be O2Marked as mobile point, slave base station O3Starting to recalculate; according to user at base station O3And base station O4By a first formula
Figure BDA0003096482760000092
Recalculating to obtain the possible coordinates of the stop point as Onew3(113.290842, 20.156390), base station O3To Onew3At a distance of 1448 m, base station O4To Onew3Is 77 meters, likewise, due to the base station O3To Onew3Is greater than the first threshold distance, so base station O3Marked as the moving point, continue from O4Starting to calculate; user at base station O4The residence time of (2) is 25 minutes, no stagnation point can be formed, so O is added5(113.2913474, 20.1568969) continue the calculation with the user at base station O5(113.2913474, 20.1568979) a residence time of 12 minutes 20 seconds; calculating a possible coordinate of a stopping point according to a first formula to be Onew4(113.291347,20.1568963), base station O4To Onew4At a distance of 4 m, base station O5To Onew4Is 7 meters, due to the base station O4To Onew4Distance and base station O5To Onew4Are all less than a first threshold distance, so base station O4And base station O5May be a potential stop point, at which point base station O in Table 1 continues to be added6(113.3013475, 20.2568979) further calculate, i.e., join base station O6Further calculating by a first formula to obtain a possible coordinate of the stopping point as Onew5Then calculates base station O4To Onew5Distance of, base stationO5To Onew5Distance and base station O6To Onew5The distance of (d); finding a distance in which there is more than a first threshold value; at this time, as can be seen from table 1, the user is at base station O4(113.2913473, 20.1468969) has a dwell time of 25 minutes, user at base station O5(113.2913474, 20.1568979) is 12 minutes 20 seconds, so that the user can be further calculated to be in the base station O4And base station O5Is 37 minutes 20 seconds, is greater than a first threshold time, and thus, the dwell point can be determined to be Onew4(113.291347,20.1568963)。
The method for identifying the mobile phone signaling data dwell point has the following technical effects:
the embodiment of the invention sorts the mobile phone signaling data of the same user according to the generation time of the signaling record by acquiring the mobile phone signaling data to obtain a sorted list; calculating the continuous stay time of the user in the first base station according to the sorted list, and determining the first base station as a stay point if the continuous stay time is greater than a first threshold time; if the continuous staying time is not more than the first threshold value, acquiring a first distance group and a first time sum, and determining a staying point according to the first distance group and the first time sum; therefore, the user staying points are accurately identified, the method is suitable for extracting the trip chain of the mobile phone user, and the mobile phone signaling data can be better used for the research in the fields of urban traffic planning, population distribution, privacy protection and the like.
Referring to fig. 2, an embodiment of the present invention further provides a device 200 for processing signaling data of a mobile phone, which specifically includes:
at least one processor 210;
at least one memory 220 for storing at least one program;
when executed by the at least one processor 210, causes the at least one processor 210 to implement the method as shown in fig. 1.
The memory 220, which is a non-transitory computer-readable storage medium, may be used to store non-transitory software programs and non-transitory computer-executable programs. The memory 220 may include high speed random access memory and may also include non-transitory memory, such as at least one magnetic disk storage device, flash memory device, or other non-transitory solid state storage device. In some embodiments, memory 220 may optionally include remote memory located remotely from processor 210, and such remote memory may be connected to processor 210 via a network. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.
It will be understood that the device configuration shown in fig. 2 is not intended to be limiting of device 200, and may include more or fewer components than shown, or some components may be combined, or a different arrangement of components.
In the apparatus 200 shown in fig. 2, the processor 210 may retrieve the program stored in the memory 220 and execute, but is not limited to, the steps of the embodiment shown in fig. 1.
The above-described embodiments of the apparatus 200 are merely illustrative, and the units illustrated as separate components may or may not be physically separate, may be located in one place, or may be distributed over a plurality of network units. Some or all of the modules can be selected according to actual needs to achieve the purposes of the embodiments.
Embodiments of the present invention also provide a computer-readable storage medium, which stores a program executable by a processor, and the program executable by the processor is used for implementing the method shown in fig. 1 when being executed by the processor.
The embodiment of the application also discloses a computer program product or a computer program, which comprises computer instructions, and the computer instructions are stored in a computer readable storage medium. The computer instructions may be read by a processor of a computer device from a computer-readable storage medium, and executed by the processor to cause the computer device to perform the method illustrated in fig. 1.
It will be understood that all or some of the steps, systems of methods disclosed above may be implemented as software, firmware, hardware, and suitable combinations thereof. Some or all of the physical components may be implemented as software executed by a processor, such as a central processing unit, digital signal processor, or microprocessor, or as hardware, or as an integrated circuit, such as an application specific integrated circuit. Such software may be distributed on computer readable media, which may include computer storage media (or non-transitory media) and communication media (or transitory media). The term computer storage media includes volatile and nonvolatile, removable and non-removable media implemented in any method or technology for storage of information such as computer readable instructions, data structures, program modules or other data, as is well known to those of ordinary skill in the art. Computer storage media includes, but is not limited to, RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, Digital Versatile Disks (DVD) or other optical disk storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired information and which can accessed by a computer. In addition, communication media typically embodies computer readable instructions, data structures, program modules or other data in a modulated data signal such as a carrier wave or other transport mechanism and includes any information delivery media as known to those skilled in the art.
The embodiments of the present invention have been described in detail with reference to the accompanying drawings, but the present invention is not limited to the above embodiments, and various changes can be made within the knowledge of those skilled in the art without departing from the gist of the present invention.

Claims (10)

1. A method for identifying a mobile phone signaling data stop point is characterized by comprising the following steps:
acquiring mobile phone signaling data, wherein the mobile phone signaling data comprises a signaling record number, a signaling user number, a base station longitude, a base station latitude and signaling record generation time;
sequencing the signaling data of the same user mobile phone according to the generation time of the signaling record to obtain a sequencing list;
calculating the continuous stay time of the user in a first base station according to the sorted list, wherein the first base station is any base station where the user stays;
if the continuous stay time is longer than a first threshold time, determining the first base station as a stay point;
if the continuous staying time is not more than the first threshold time, acquiring a first distance group and a first time sum, wherein the first distance group comprises the distance between the first base station and a first target base station and the distance between a second base station and the first target base station, the first time sum is the sum of the staying time of a user at the first base station and the staying time of the user at the second base station, the second base station is a base station which is not included in the first base station and has the continuous staying time not more than the first threshold time, and the first target base station is a possible base station calculated according to the first base station and the second base station;
and determining a stopping point according to the first distance group and the first time sum.
2. The method as claimed in claim 1, wherein the step of calculating the continuous staying time of the user at the first base station according to the sorted list comprises:
acquiring a first time and a second time according to the sorted list, wherein the first time is the earliest signaling record generation time corresponding to the first base station, and the second time is the latest signaling record generation time corresponding to the first base station;
and calculating the continuous staying time of the user in the first base station according to the first time and the second time.
3. The method for identifying a mobile phone signaling data dwell point according to claim 1, wherein the first distance group is obtained by the following steps:
acquiring a first information group according to the sorted list, wherein the first information group comprises longitude and latitude coordinates of the first base station, longitude and latitude coordinates of the second base station, earliest signaling record generation time corresponding to the first base station, latest signaling record generation time corresponding to the first base station, earliest signaling record generation time corresponding to the second base station and latest signaling record generation time corresponding to the second base station;
according to the first information group, calculating to obtain longitude and latitude coordinates of the first target base station;
and calculating to obtain the first distance group according to the longitude and latitude coordinates of the first base station, the longitude and latitude coordinates of the second base station and the longitude and latitude coordinates of the first target base station.
4. The method for identifying a mobile phone signaling data stop point according to claim 3, wherein the first time sum is obtained by the following steps:
calculating the staying time of the user in the first base station according to the earliest signaling record generation time corresponding to the first base station and the latest signaling record generation time corresponding to the first base station;
calculating the staying time of the user in the second base station according to the earliest signaling record generation time corresponding to the second base station and the latest signaling record generation time corresponding to the second base station;
and adding the stay time of the user in the first base station and the stay time of the user in the second base station to obtain a first time sum.
5. The method for identifying a mobile phone signaling data dwell point according to claim 3, wherein the latitude and longitude coordinates of the first target base station are calculated by a first formula, and the first formula is:
Figure FDA0003096482750000021
in the formula, longnewIs the channel of the first target base stationScale, latnewIs the latitude, t, of the first target base station1Generating a time, t, for the earliest signalling record of a base station inGenerating a time, t, for the latest signalling record of base station iiFor the stay time of the user at base station i, lngiIs the longitude, lat, of base station iiIs the latitude of base station i, i denotes the base station.
6. The method for identifying a stop point of signaling data of a mobile phone according to claim 5, wherein the step of determining a stop point according to the first distance group and the first time sum comprises:
if the distance between the first base station and the first target base station and the distance between the second base station and the first target base station are both smaller than a first threshold distance, extracting a third base station, wherein the third base station is a base station except the first base station and the second base station;
acquiring a second information group according to the ordered list, wherein the second information group comprises the first information group, longitude and latitude coordinates of the third base station, earliest signaling record generation time corresponding to the third base station and latest signaling record generation time corresponding to the third base station;
according to the second information group, calculating to obtain the longitude and latitude coordinates of a second target base station through the first formula;
calculating to obtain a second distance group according to the longitude and latitude coordinates of the first base station, the longitude and latitude coordinates of the second base station, the longitude and latitude coordinates of the third base station and the longitude and latitude coordinates of the second target base station, wherein the second distance group comprises the distance from the first base station to the second target base station, the distance from the second base station to the second target base station and the distance from the third base station to the second target base station;
if at least one distance in the second distance group is greater than the first threshold distance, judging whether the first time sum is greater than the first threshold time;
and if the first time sum is greater than the first threshold time, determining that the first target base station is a stop point.
7. The method for identifying the mobile phone signaling data stop point according to claim 6, wherein:
if all the distances in the second distance group are smaller than the first threshold distance, extracting a fourth base station, wherein the fourth base station is a base station except the first base station, the second base station and the third base station;
calculating to obtain a second time sum according to the second information group, wherein the second time sum is the sum of the staying time of the user in the first base station, the second base station and the third base station;
acquiring a third information group according to the ordered list, wherein the second information group comprises the second information group, longitude and latitude coordinates of the fourth base station, earliest signaling record generation time corresponding to the fourth base station and latest signaling record generation time corresponding to the fourth base station;
according to the third information group, calculating to obtain longitude and latitude coordinates of a third target base station through the first formula;
calculating a third distance group according to the longitude and latitude coordinates of the first base station, the longitude and latitude coordinates of the second base station, the longitude and latitude coordinates of the third base station, the longitude and latitude coordinates of the fourth base station and the longitude and latitude coordinates of the third target base station, wherein the third distance group comprises the distance from the first base station to the third target base station, the distance from the second base station to the third target base station, the distance from the third base station to the third target base station and the distance from the fourth base station to the third target base station;
if at least one distance in the third distance group is greater than the first threshold distance, judging whether the second time sum is greater than the first threshold time;
and if the second time sum is greater than the first threshold time, determining that the second target base station is a stop point.
8. The method as claimed in claim 6, wherein the first base station is marked as a moving point and recalculated from the second base station to obtain the staying point if at least one of the following conditions exists:
at least one distance in the first distance group is greater than the first threshold distance;
alternatively, there is at least one distance in the second distance group that is greater than the first threshold distance but the first time sum is not greater than the first threshold time.
9. A mobile phone signaling data stop point identification device is characterized by comprising:
at least one processor;
at least one memory for storing at least one program;
when executed by the at least one processor, cause the at least one processor to implement the method of any one of claims 1-8.
10. Computer-readable storage medium, on which a processor-executable program is stored, which, when being executed by a processor, is adapted to carry out the method according to any one of claims 1-8.
CN202110613302.6A 2021-06-02 2021-06-02 Mobile phone signaling data stop point identification method, device and storage medium Active CN113473398B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110613302.6A CN113473398B (en) 2021-06-02 2021-06-02 Mobile phone signaling data stop point identification method, device and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110613302.6A CN113473398B (en) 2021-06-02 2021-06-02 Mobile phone signaling data stop point identification method, device and storage medium

Publications (2)

Publication Number Publication Date
CN113473398A true CN113473398A (en) 2021-10-01
CN113473398B CN113473398B (en) 2022-08-30

Family

ID=77872240

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110613302.6A Active CN113473398B (en) 2021-06-02 2021-06-02 Mobile phone signaling data stop point identification method, device and storage medium

Country Status (1)

Country Link
CN (1) CN113473398B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113810992A (en) * 2021-11-15 2021-12-17 广州睿帆科技有限公司 Data processing method and device
CN114501419A (en) * 2021-12-30 2022-05-13 中国联合网络通信集团有限公司 Signaling data processing method, device and storage medium
CN115297499A (en) * 2022-09-26 2022-11-04 智慧足迹数据科技有限公司 Base station operating parameter determination method based on signaling data and related device

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107040894A (en) * 2017-04-21 2017-08-11 杭州市综合交通研究中心 A kind of resident trip OD acquisition methods based on mobile phone signaling data
CN109583640A (en) * 2018-11-23 2019-04-05 东南大学 A kind of Urban Traffic passenger flow attribute recognition approach based on multi-source location data
US20190124618A1 (en) * 2017-10-23 2019-04-25 Ubtech Robotics Corp Range finding base station selection method and apparatus
CN110113718A (en) * 2019-05-31 2019-08-09 南京瑞栖智能交通技术产业研究院有限公司 A kind of railway traffic hinge type of population recognition methods based on mobile phone signaling data
CN110377687A (en) * 2019-07-26 2019-10-25 智慧足迹数据科技有限公司 User's trip mode method of discrimination, device and server
CN111464950A (en) * 2020-03-31 2020-07-28 天津市城市规划设计研究院 Method for extracting travel stop point by using mobile phone signaling data
CN111681421A (en) * 2020-06-10 2020-09-18 南京瑞栖智能交通技术产业研究院有限公司 Mobile phone signaling data-based external passenger transport hub centralized-sparse space distribution analysis method
CN111770452A (en) * 2020-05-27 2020-10-13 中山大学 Mobile phone signaling stop point identification method based on personal travel track characteristics

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107040894A (en) * 2017-04-21 2017-08-11 杭州市综合交通研究中心 A kind of resident trip OD acquisition methods based on mobile phone signaling data
US20190124618A1 (en) * 2017-10-23 2019-04-25 Ubtech Robotics Corp Range finding base station selection method and apparatus
CN109583640A (en) * 2018-11-23 2019-04-05 东南大学 A kind of Urban Traffic passenger flow attribute recognition approach based on multi-source location data
CN110113718A (en) * 2019-05-31 2019-08-09 南京瑞栖智能交通技术产业研究院有限公司 A kind of railway traffic hinge type of population recognition methods based on mobile phone signaling data
CN110377687A (en) * 2019-07-26 2019-10-25 智慧足迹数据科技有限公司 User's trip mode method of discrimination, device and server
CN111464950A (en) * 2020-03-31 2020-07-28 天津市城市规划设计研究院 Method for extracting travel stop point by using mobile phone signaling data
CN111770452A (en) * 2020-05-27 2020-10-13 中山大学 Mobile phone signaling stop point identification method based on personal travel track characteristics
CN111681421A (en) * 2020-06-10 2020-09-18 南京瑞栖智能交通技术产业研究院有限公司 Mobile phone signaling data-based external passenger transport hub centralized-sparse space distribution analysis method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
钟舒琦 等: "基于兴趣点与导航数据的手机信令数据出行方式识别", 《中山大学学报(自然科学版)》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113810992A (en) * 2021-11-15 2021-12-17 广州睿帆科技有限公司 Data processing method and device
CN114501419A (en) * 2021-12-30 2022-05-13 中国联合网络通信集团有限公司 Signaling data processing method, device and storage medium
CN114501419B (en) * 2021-12-30 2023-05-12 中国联合网络通信集团有限公司 Signaling data processing method, apparatus and storage medium
CN115297499A (en) * 2022-09-26 2022-11-04 智慧足迹数据科技有限公司 Base station operating parameter determination method based on signaling data and related device
CN115297499B (en) * 2022-09-26 2022-12-16 智慧足迹数据科技有限公司 Base station parameter determining method based on signaling data and related device

Also Published As

Publication number Publication date
CN113473398B (en) 2022-08-30

Similar Documents

Publication Publication Date Title
CN113473398B (en) Mobile phone signaling data stop point identification method, device and storage medium
CN110166943B (en) Method for processing terminal position information
US9866997B2 (en) Systems and methods for geo-location based message streams
US7640027B2 (en) Mobile communication terminal capable of determining position and method of determining position in the mobile communication terminal
US20190301891A1 (en) Method and Apparatus for Obtaining and Displaying Map Data On a Mobile Device
US8542637B2 (en) Clustering crowd-sourced data for determining beacon positions
US7890120B2 (en) Method of locating a cellular terminal
JP5705110B2 (en) How to create a speed estimate
US8483962B2 (en) Method for providing route information and the system thereof
CN110213724B (en) Pseudo base station motion trajectory identification method
CN104380816A (en) Obtaining timing of lte wireless base stations using aggregated otdoa assistance data
CN109195113B (en) Method and device for identifying travel track of user and computer equipment
US20130166352A1 (en) Mobile categorization
CN107871400B (en) Road network information updating method and device
JP6564382B2 (en) Apparatus and method for probe data management
CN111222375A (en) Vehicle driving track generation method, system and server
CN111427983A (en) Service method, system, device and storage medium based on geographic information retrieval
CN111414558B (en) Navigation route sending method, navigation route displaying method, navigation route sending device, navigation route displaying server and navigation route displaying medium
CN109068282B (en) Method and device for identifying user travel scene and computer equipment
CN107274705A (en) The determination methods, devices and systems of position of bus
CN100387098C (en) Information providing method based on position and system
CN112153568B (en) Wi-Fi identification and binding method, device and equipment based on service scene
CN113271557A (en) Mobile phone signaling data processing method, device and storage medium
CN111050293A (en) Data processing method and device and electronic equipment
US20160157064A1 (en) Method and notifying unit for providing a notification about an event

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