A method of personnel trajectory backtracking based on communication data

Abstract

The present application relates to a kind of personnel trajectory backtracking method based on communication data, belong to personnel trajectory backtracking technical field.The personnel trajectory backtracking method based on communication data of the present application considers the trajectory backtracking problem of personnel who simultaneously possesses multiple phone cards, multiple communication terminals, this method is based on the association between the user identification number and the user equipment number generated by device card, with database software as a tool, the longitude and latitude coordinates of the stop point corresponding to multiple user identification numbers and equipment number are fused in time sequence, compared with the personnel trajectory backtracking method based on only communication terminal device geographic location information, the personnel trajectory backtracked is more complete, continuous.

Description

Technical Field

[0001] This invention relates to a method for tracing personnel trajectories based on communication data, belonging to the field of personnel trajectory tracing technology. Background Technology

[0002] Mobile phones have become an indispensable part of people's lives, and extracting their location information can directly reflect a person's movement trajectory. Tracing these trajectories has profound implications for research on human travel behavior and the construction of "electronic fences." Existing technologies utilize signaling data from mobile communication networks to track the unique IMEI device number of a communication terminal, collect the geographical coordinates (latitude and longitude) of all communication base stations communicating with that device number, and sort them by time to form the communication terminal's movement trajectory, which is generally considered to be the person's movement trajectory. Another approach uses the source addresses of base stations communicating with the terminal and map data to match geographical coordinates with roads, tracing back to more accurate human trajectories. Data extracted based on these methods, including human trajectories and device activation / deactivation behavior, can be used for personnel monitoring and other applications.

[0003] Currently, with the support of dual-SIM dual-standby technology in communication terminals, people can seamlessly switch between using data cards and secondary SIM cards from different operators. They often purchase multiple communication terminals and use them at different times, leading to an increasing trend of using multiple SIM cards on multiple communication terminals. Having multiple SIM cards and terminals makes it extremely difficult to trace a complete person's trajectory using their International Mobile Equipment Identity (IMEI). If only a single IMEI is used for trajectory tracing, the traced trajectory cannot reflect the trajectories of other communication terminals used by the same person, resulting in gaps and breaks in the trajectory. These gaps and breaks in the trajectory lose many features, leading to a decrease in the accuracy of applications based on this trajectory. Summary of the Invention

[0004] The technical problem solved by this invention is to overcome the shortcomings of the prior art and propose a method for tracing personnel trajectories based on communication data. This method takes into account the behavior of personnel with multiple phone cards and multiple communication terminals, and uses communication data such as International Mobile Subscriber Identification Number (IMSI), IMEI device number, and communication relationships to trace the complete trajectory of personnel.

[0005] The technical solution of this invention is:

[0006] A method for tracing personnel trajectories based on communication data, the method comprising the following steps:

[0007] The first step is to collect signaling data between the communication terminal and the communication base station, and then construct a database based on the collected signaling data. The constructed database includes a communication database and a communication relationship database. The signaling data includes at least the user identification number, the device number corresponding to the user identification number, the latitude and longitude coordinates of the communication base station's stopping point, the communication termination time of the device number, and the communication date. The communication relationship database includes at least the calling user identification number, the calling user identification number's home location, the called user identification number, the communication termination time of the device number, and the communication date.

[0008] The second step, within the set number of days, extracts all device numbers (International Mobile Equipment Identity, IMEI) associated with any user identification number (International Mobile Subscriber Identification Number, IMSI) used by each person whose trajectory is to be traced from the communication database built in the first step, as the device number to be queried;

[0009] The third step is to extract all signaling data containing any of the device numbers to be queried obtained in the second step from the communication database constructed in the first step within a set number of days, and use this data as the signaling dataset to be queried; and to extract all user identification numbers appearing in the signaling dataset to be queried, and use this data as the user identification numbers to be queried.

[0010] The fourth step is to use the signaling dataset extracted in the third step to count the number of different user identification numbers associated with each device number to be queried in the second step. When the number of user identification numbers associated with any device number to be queried is more than six, it is considered an abnormal situation and will not be processed, i.e., no personnel trajectory backtracking will be performed.

[0011] When the number of user identification numbers associated with all the device numbers in the signaling data extracted in the third step is less than six, the signaling data extracted in the third step is sorted according to the communication termination time of the device number. Then, the latitude and longitude coordinates of the stopping point of the communication base station included in the sorted signaling data are the preliminary tracing trajectory of the personnel's SIM card.

[0012] The fifth step involves using the communication database established in the first step to extract all appearing user identification numbers within the set geographical and time range of the initial SIM card tracing trajectory endpoint, thus creating a user set for SIM card and device replacement screening.

[0013] The sixth step is to extract signaling data from the communication relationship database built in the first step, within the set number of days, for users whose calling user identification number belongs to the user set for SIM card and device replacement screening established in the fifth step. This includes the calling user identification number, the location of the calling user identification number, and the called user identification number, and is used as the general database of personnel communication characteristics.

[0014] Step 7: Using the total database of personnel communication characteristics obtained in Step 6, extract the data whose caller ID is the same as the user ID to be queried obtained in Step 3, and denote it as set A. Extract the data whose caller ID is the same as the user set of users who have changed cards or devices obtained in Step 5, and denote it as set B.

[0015] Extract the calling user identification numbers from set B that share the same called user identification number and the same home location as set A, and use them as set C. Then, using the communication database constructed in the first step, extract the signaling data from the communication database within a set number of days whose user identification numbers are the same as those in set C, as the preliminary communication backtracking trajectory of the personnel. Merge the preliminary communication backtracking trajectory of the personnel with the preliminary SIM card backtracking trajectory of the personnel in the fourth step, and sort them according to the communication termination time of the device number. The latitude and longitude coordinates of the stopping point of the communication base station included in the signaling data are the final trajectory of the personnel to be traced.

[0016] When set A and set B do not share a common called user identification number or a common home location, it means that the communication relationship database has not detected the behavior of the person to be traced changing both the SIM card and the device. The trajectory tracing ends, and the final trajectory of the person to be traced is the preliminary SIM card tracing trajectory.

[0017] In the first step, the data center collects signaling data from communication terminals and communication base stations. Using data processing software capable of storing, querying, and analyzing this data, including Oracle databases, it constructs a communication database and a communication relationship database. The communication database, represented by COM, contains the following fields: IMSI, IMEI, Lng, Lat, Time, and Date. These represent: user identification number, location of the user identification number, device number corresponding to the user identification number, latitude and longitude coordinates of the communication base station's stopping point, communication termination time of the device number, and communication date, respectively. The communication relationship database, represented by TEL, contains the following fields: IMSI, Area, Call_IMSI, Time, and Date. These represent: calling user identification number, location of the calling user identification number, called user identification number, communication termination time of the device number, and communication date, respectively.

[0018] In the second step, using database query statements, including but not limited to Structured Query Language (SQL), all device numbers associated with any user identification number used by the person whose trajectory is being traced are extracted from the communication database constructed in the first step. The query object is limited to COM, the query start time is limited to date T, the query time range is limited to T ≤ COM.Date < T + N, and the query condition is that COM.IMSI is the same as any user identification number used by the person whose trajectory is being traced. The IMEI field is set as the query field. After deduplicating the obtained query results, the device number J1 to be queried is obtained.

[0019] In the third step, using database query statements, including but not limited to Structured Query Language, all signaling data containing any of the device numbers to be queried obtained in the second step are extracted from the communication database constructed in the first step. The query object is limited to COM, the query time start point is limited to date T, the query time range is limited to T≤COM.Date<T+N, and the query condition is limited to COM.IMEI=J1.IMEI. The query fields are set to IMSI, IMEI, Lng, Lat, Time, Date, resulting in the signaling dataset J2 to be queried. The J2.IMSI field is then deduplicated to obtain the user identification number J3 to be queried.

[0020] In the fourth step, the grouping function of the database is used to count the number N of different user identification numbers associated with the i-th IMEI field data in J2. i i = 1, 2, ..., N cn , where N cn This indicates the number of distinct IMEI field data in J2.

[0021] Step 4.1 If any N i The following conditions must be met:

[0022] N i ≥6

[0023] This is an abnormal situation and will not be handled; that is, no personnel trajectory will be traced back.

[0024] Step 4.2 If the conditions in Step 4.1 are not met, add a Name field to the signaling dataset J2 to be queried as a unique identifier to distinguish the personnel whose trajectory is to be traced back. The Name field is determined by taking the minimum value of all IMSIs in J2. Using the database sorting function, limit the query object to J2, limit the sorting field to Time, and limit the query fields to IMSI, Lng, Lat, Time, Date, Name to obtain the preliminary machine card tracing trajectory J5 of the personnel.

[0025] In the fifth step, the query object is limited to the preliminary personnel card retrospective trajectory J5. The data with the largest corresponding Time value for each Name field is extracted to obtain the intermediate result J9. J9.Time represents the time of the endpoint of the preliminary personnel card retrospective trajectory, and J9.Lng and J9.Lat represent the latitude and longitude coordinates of the endpoint of the preliminary personnel card retrospective trajectory. The query object is limited to the communication database COM, and the query time range is limited to T. R Limit the geographical coordinate range G LAT With G LNG Querying the COM.IMSI field yields the user set J10 for SIM card / device replacement screening, with the following query conditions:

[0026]

[0027] In the sixth step, the query time start point is limited to date T, the number of days is limited to M, the query object is limited to the communication database TEL constructed in the first step, the query conditions are limited to TEL.IMSI = J10.IMSI and T≤TEL.Date<T+M, and the query fields are IMSI, Area, Call_IMSI, Time, Date, to obtain the total database of personnel communication characteristics J11.

[0028] In the seventh step, the query object is limited to the personnel communication feature database J11, and the query condition is J11.IMSI = J3.IMSI. J11.IMSI, J11.Area, J11.Call_IMSI, J11.Time, and J11.Date are extracted and denoted as set A. The query condition is J11.IMSI ! = J3.IMSI. J11.IMSI, J11.Area, J11.Call_IMSI, J11.Time, and J11.Date are extracted and denoted as set B. Set A and set B are combined to obtain set C, with the joint condition A.Call_IMSI = B.Call_IMSI and A.Area = B.Area.

[0029] Step 7.1 If set C is not empty, limit the query time starting point to date T, limit the number of days to M, limit the query object to the communication database COM, and limit the query conditions to COM.IMSI = C.IMSI and T ≤ COM.Date < T + M. Extract COM.IMSI, COM.Area, COM.Call_IMSI, COM.Time, and COM.Date to obtain the preliminary communication tracing trajectory J12 of the personnel. Merge J12 and J5, and fill the Name field of the merged data with the Name field from J5. After sorting by the communication termination time of the device number, the latitude and longitude coordinates of the stopping point of the communication base station included in the signaling data are the final trajectory J13 of the personnel to be traced.

[0030] Step 7.2 If set C is empty, end the trajectory backtracking. The final trajectory of the person to be backtracked is the preliminary machine card backtracking trajectory.

[0031] Beneficial effects

[0032] 1. The present invention provides a method for tracing personnel trajectories based on communication data, which considers the problem of tracing the trajectories of personnel who simultaneously possess multiple phone cards and multiple communication terminals. This method is based on the association between the user identification number generated by the device card and the user device number, and uses database software as a tool to merge the latitude and longitude coordinates of the stopping points corresponding to multiple user identification numbers and device numbers in chronological order. Compared with the method of tracing personnel trajectories based solely on the geographical location information of communication terminal devices, the traced personnel trajectories are more complete and continuous.

[0033] 2. The method of the present invention takes into account the situation where the same device number is associated with too many user identification numbers, so as to avoid the personnel trajectory being repeated and overlapping due to the influence of abnormal situations.

[0034] 3. For certain individuals who change both their SIM cards and devices in a particular location, the association between the user identification number and the device number is no longer sufficient to trace their movements after the change. Therefore, this patent analyzes phone numbers within a specific geographical location and time range of the trajectory endpoint. Phone cards belonging to the same region and having dialed the same number over a period of time are considered to belong to the same individual. By fusing the latitude and longitude coordinates of these individuals' stops, a more complete trajectory can be traced back.

[0035] 4. Existing technologies do not consider the problem of discontinuous and interrupted tracing data caused by the same user constantly switching SIM cards and devices. This technology first associates the user's number with the IMSI, and then associates the IMSI with the IMEI. This allows for the integration of data before and after a person changes SIM cards but not devices, and vice versa, which helps to extract more complete trajectory data. If this trajectory is used for personnel analysis, the accuracy will be improved.

[0036] 5. In addition to considering the case where personnel change SIM cards without changing devices or change devices without changing SIM cards, this application also considers the case where personnel change both SIM cards and devices at the same time. By extracting communication relationships and location information near the endpoint of the trajectory, it determines whether the previous and subsequent communication trajectories belong to the same person. If so, the trajectories are merged, thus making the traced user trajectory more comprehensive. Detailed Implementation

[0037] The present invention will be further described below with reference to the embodiments.

[0038] Example

[0039] The first step is for the data center to collect signaling information from communication terminals and communication base stations, and construct a communication database and a communication relationship database, as shown in Table 1 and Table 2, respectively.

[0040] Table 1 Communication Database COM

[0041] IMSI IMEI Lng Lat Time Date 46**1 86**1 105.5 40.1 06010918 20230601 46**2 86**2 104.2 39.1 06022310 20230602 … … … … … …

[0042] Table 2. Connection Relationship Database TEL

[0043] IMSI Area Call_IMSI Time Date 46**1 0001 46**8 06010918 20230601 46**2 0002 46**9 06022310 20230602 … … … … …

[0044] The second step involves using Oracle database as the database software, limiting the query object to the communication database COM, and using the WHERE clause to specify the query start time T as "20230601", the query time range N as 7 days, and the query end time as "20230608". Using the SELECT clause, with the conditions "20230601" ≤ COM.Date < "20230608" and COM.IMSI matching any user identification number used by the person whose trajectory is being traced, the query field is set to IMEI. After deduplicating the query results using the DISTINCT keyword, the device number J1 to be queried is obtained, as shown in Table 3.

[0045] Table 3 Device Number to be Queryed J1

[0046] IMEI 86************1 86************2

[0047] The third step involves using an Oracle database, limiting the query object to the communication database COM, setting the query start time T to "20230601", the query time range N to 7 days, and the query end time to "20230608". A WHERE clause is used, with "20230601" ≤ COM.Date < "20230608" and COM.IMEI = J1.IMEI as query conditions. The extracted fields are IMSI, IMEI, Lng, Lat, Time, and Date, yielding the signaling dataset J2 to be queried, as shown in Table 4. The J2.IMSI field is then extracted and deduplicated to obtain the user identification number J3 to be queried, as shown in Table 5.

[0048] Table 4 Signaling Dataset J2 to be Queryed

[0049] IMSI IMEI Lng Lat Time Date 46**1 86**1 105.5 40.1 06010918 20230601 46**2 86**2 104.2 39.1 06022310 20230602 … … … … … … 46**3 86**1 105.2 38.5 06071105 20230607

[0050] Table 5 User ID J3 to be Queryed

[0051] IMSI 46************1 46************2 46************3

[0052] The fourth step involves using the database's grouping function to count the number of distinct user identification numbers associated with the i-th IMEI field in J2: the number of 86**1 ​​is 2, and the number of 86**2 is 1. Where N cn This represents the number of distinct IMEI field data in J2. Since 2 < 6 and 1 < 6, trajectory backtracking is performed.

[0053] Add a Name field to J2 as a unique identifier to distinguish the personnel whose trajectory is to be traced. The Name field is determined by taking the minimum value of all IMSI values ​​in J2. Using the database sorting function, limit the query object to J2, limit the sorting field to Time, and limit the query fields to IMSI, Lng, Lat, Time, Date. Name is used to obtain the preliminary machine card tracing trajectory J5 of the personnel, as shown in Table 6.

[0054] Table 6 Preliminary Machine Card Retrospective Trajectory J5

[0055] NAME IMSI IMEI Lng Lat Time Date 46**1 46**1 86**1 105.5 40.1 06010918 20230601 46**1 46**2 86**2 104.2 39.1 06022310 20230602 … … … … … … … 46**1 46**3 86**1 105.2 38.5 06071105 20230607

[0056] Fifth, limit the query object to the preliminary machine card retrospective trajectory J5 of personnel, extract the data with the largest Time corresponding to each Name field, and obtain the intermediate result J9, as shown in Table 7.

[0057] Table 7 Intermediate Results J9

[0058] NAME IMSI IMEI Lng Lat Time Date 46**1 46**3 86**1 105.2 38.5 06071105 20230607

[0059] The obtained J9.Time represents the time of the initial personnel SIM card tracing trajectory endpoint, and J9.Lng and J9.Lat represent the latitude and longitude coordinates of the endpoint of the initial personnel SIM card tracing trajectory endpoint. The query object is limited to the communication database COM, and the query time range is limited to T. R =30, limiting the geographical coordinate range G LAT =0.05 and G LNG =0.05. Using the SELECT statement to query the COM.IMSI field, the user set J10 for card and device replacement screening is obtained, as shown in Table 8:

[0060] Table 8 Screening User Set for Card / Equipment Replacement J10

[0061] IMSI IMEI Lng Lat Time Date 46**3 86**1 105.2 38.5 06071105 20230607 46**4 86**4 105.22 38.52 06071122 20230607 … … … … … … 46**6 86**6 105.23 38.47 06071133 20230607

[0062] The query conditions are:

[0063]

[0064] Step 6: Limit the query time to the date "20230601", limit the number of days to 10 days, limit the query object to the communication database TEL constructed in Step 1, and limit the query conditions to TEL.IMSI = J10.IMSI and "20230601" ≤ TEL.Date < "20230611". Use the SELECT statement to query the fields IMSI, Area, Call_IMSI, Time, and Date to obtain the total database of personnel communication characteristics J11, as shown in Table 9.

[0065] Table 9. Personnel Communication Characteristics Database J11

[0066] IMSI Area Call_IMSI Time Date 46**3 0001 46**8 06071105 20230607 46**4 0002 46**9 06080920 20230608 … … … … … 46**6 0001 46**8 06101502 20230610

[0067] Step 7: Limit the query object to the personnel communication feature database J11. Use SELECT and WHERE statements to limit the query condition to J11.IMSI = J3.IMSI, and extract J11.IMSI, J11.Area, J11.Call_IMSI, J11.Time, and J11.Date, denoted as set A. Limit the query condition to J11.IMSI ! = J3.IMSI, and extract J11.IMSI, J11.Area, J11.Call_IMSI, J11.Time, and J11.Date, denoted as set B. Use the JOIN method to combine set A and set B to obtain set C, with the combination condition A.Call_IMSI = B.Call_IMSI and A.Area = B.Area, retaining only the B.IMSI field as IMSI. The resulting set C is shown in Table 10.

[0068] Table 10 Set C

[0069] IMSI 46**6

[0070] If set C is not empty, as shown in Table 10, use SELECT and WHERE statements, limiting the query time start point to "20230601", the number of days to 10 days, the query object to the communication database COM, and the query conditions to COM.IMSI = C.IMSI and "20230601" ≤ COM.Date < "20230611". Extract COM.IMSI, COM.Area, COM.Call_IMSI, COM.Time, and COM.Date to obtain the preliminary communication tracing trajectory J12, as shown in Table 11. Merge J12 and J5, and fill the Name field of the merged data with the Name field from J5. After sorting by the communication termination time of the device number, the latitude and longitude coordinates of the stopping point of the communication base station included in the signaling data are the final trajectory of the person to be traced, resulting in J13, as shown in Table 12.

[0071] Table 11 Initial Contact Tracing Trajectory of Personnel J12

[0072] IMSI IMEI Lng Lat Time Date 46**6 86**6 105.29 38.5 06071133 20230607 46**6 86**6 105.90 38.65 06080955 20230608 … … … … … … 46**6 86**6 107.64 39.74 06101354 20230610

[0073] Table 12 Final Trajectory of Personnel to be Retrospectively Traced J13

[0074] NAME IMSI IMEI Lng Lat Time Date 46**1 46**1 86**1 105.5 40.1 06010918 20230601 46**1 46**2 86**2 104.2 39.1 06022310 20230602 … … … … … … … 46**1 46**3 86**1 105.2 38.5 06071105 20230607 46**1 46**6 86**6 105.29 38.5 06071133 20230607 46**1 46**6 86**6 105.90 38.65 06080955 20230608 46**1 … … … … … … 46**1 46**6 86**6 107.64 39.74 06101354 20230610

[0075] If set C is empty, the trajectory backtracking ends, and the final trajectory of the person to be backtracked is the preliminary machine card backtracking trajectory J5.

[0076] In summary, the above are merely preferred embodiments of the present invention and are not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A method for tracing personnel trajectories based on communication data, characterized in that... The steps of this method include: The first step is to collect signaling data between the communication terminal and the communication base station, and to build a database based on the collected signaling data. The database includes a communication database and a communication relationship database. The second step is to extract all device numbers associated with any user identification number used by each person whose trajectory is to be traced from the communication database built in the first step as the device numbers to be queried within the set number of days. The third step is to extract all signaling data containing any of the device numbers to be queried obtained in the second step from the communication database constructed in the first step within a set number of days, and use this data as the signaling dataset to be queried; then, extract all user identification numbers appearing in the signaling dataset to be queried, and use this data as the user identification number to be queried. The fourth step is to use the signaling dataset extracted in the third step to count the number of different user identification numbers associated with each device number to be queried in the second step. When the number of user identification numbers associated with any device number to be queried is more than six, it is considered an abnormal situation and will not be processed, i.e., no personnel trajectory backtracking will be performed. When the number of user identification numbers associated with all the device numbers to be queried in the signaling dataset extracted in the third step is less than six, add a Name field to the signaling dataset to be queried, and sort the signaling data after adding the Name field according to the communication termination time of the device number. Then the latitude and longitude coordinates of the stopping point of the communication base station in the sorted signaling data are the preliminary tracing trajectory of the personnel's SIM card. The fifth step involves using the communication database established in the first step to extract all appearing user identification numbers within the set geographical and time range of the initial SIM card tracing trajectory endpoint, as a user set for SIM card and device replacement screening. Step 6: Within the set number of days, extract the signaling data of the calling user identification number belonging to the user set of card and device replacement screening established in Step 5 from the communication relationship database built in Step 1. The calling user identification number, the location of the calling user identification number, and the called user identification number in the signaling data are used as the total database of personnel communication characteristics. Step 7: Using the total database of personnel communication characteristics obtained in Step 6, extract the data whose caller ID is the same as the user ID to be queried obtained in Step 3, and denote it as set A. Extract the data whose caller ID is the same as the user set of users who have changed cards or devices obtained in Step 5, and denote it as set B. Extract the calling user identification numbers from set B that share the same called user identification number and the same home location as set A, and use them as set C. Then, using the communication database constructed in the first step, extract the signaling data from the communication database within a set number of days whose user identification numbers are the same as those in set C, as the preliminary communication backtracking trajectory of the personnel. Merge the preliminary communication backtracking trajectory of the personnel with the preliminary SIM card backtracking trajectory of the personnel in the fourth step, and sort them according to the communication termination time of the device number. The latitude and longitude coordinates of the stopping point of the communication base station included in the signaling data are the final trajectory of the personnel to be traced. When set A and set B do not share a common called user identification number or a common location, it means that the communication relationship database has not detected the behavior of the person to be traced changing both the SIM card and the device. The trajectory tracing ends, and the final trajectory of the person to be traced is the preliminary SIM card tracing trajectory.

2. The method for tracing personnel trajectories based on communication data according to claim 1, characterized in that: In the first step, when collecting signaling data, an Oracle database is used. The communication database is represented by COM and contains the following fields: IMSI, IMEI, Lng, Lat, Time, Date, which have the following meanings: user identification number, location of user identification number, device number corresponding to user identification number, latitude and longitude coordinates of the communication base station's stopping point, communication termination time of device number, and communication date. The communication relationship database is represented by TEL and includes the following fields: IMSI, Area, Call_IMSI, Time, and Date, which mean: calling user identification number, the location of the calling user identification number, the called user identification number, the communication termination time of the device number, and the communication date, respectively.

3. A method for tracing personnel trajectories based on communication data according to claim 1 or 2, characterized in that: In the second step, the method for extracting all device numbers associated with any user identification number used by each person whose trajectory is to be traced is as follows: Using database query statements, including Structured Query Language, in the communication database constructed in the first step, extract all device numbers associated with any user identification number used by the person whose trajectory is to be traced. Limit the query object to COM, limit the query time start to date T, limit the query time range T≤COM.Date<T+N, where N is a set number of days, limit the query condition that COM.IMSI is the same as any user identification number used by the person whose trajectory is to be traced, set the query field to IMEI, and after deduplicating the obtained query results, obtain the device number J1 to be queried.

4. The method for tracing personnel trajectories based on communication data according to claim 3, characterized in that: In the third step, the method for extracting all signaling data containing any one of the device numbers to be queried obtained in the second step is as follows: Using database query statements, including but not limited to Structured Query Language, extract all signaling data containing any of the device numbers to be queried obtained in the second step from the communication database constructed in the first step. Limit the query object to COM, limit the query time start to date T, limit the query time range T≤COM.Date<T+N, limit the query condition to COM.IMEI=J1.IMEI, and set the query fields to IMSI, IMEI, Lng, Lat, Time, Date to obtain the signaling dataset J2 to be queried. The method for extracting all user identification numbers appearing in the signaling dataset to be queried is as follows: After removing duplicates from the J2.IMSI field, the user identification number J3 to be queried is obtained.

5. The method for tracing personnel trajectories based on communication data according to claim 4, characterized in that: In the fourth step, the method for counting the number of distinct user identification numbers associated with each queried device number in the second step is as follows: Using the grouping function of the database, count the number N distinct user identification numbers associated with the i-th IMEI field in J2. i i = 1, 2, ..., N cn , where N cn This indicates the number of distinct IMEI field data in J2; In the fourth step, the sorted signaling data includes IMSI, Lng, Lat, Time, Date, and Name; the Name field is determined by taking the minimum value of all IMSIs in J2.

6. A method for tracing personnel trajectories based on communication data according to claim 4 or 5, characterized in that: In the fifth step, the method for extracting all appearing user identification numbers is as follows: The query is limited to the preliminary personnel SIM card retrospective trajectory J5. The data with the largest corresponding Time value for each Name field is extracted to obtain intermediate result J9. J9.Time represents the time of the endpoint of the preliminary personnel SIM card retrospective trajectory, and J9.Lng and J9.Lat represent the latitude and longitude coordinates of the endpoint. The query is limited to the communication database COM, and the query time range is limited to T. R Limit the geographical coordinate range G LAT With G LNG Querying the COM.IMSI field yields the user set J10 for SIM card / device replacement screening. The specific query conditions are:

7. The method for tracing personnel trajectories based on communication data according to claim 6, characterized in that: In the sixth step, the method for extracting signaling data from the user set of users whose calling user identification number belongs to the SIM card / device change screening user set established in the fifth step is as follows: The query time is limited to date T, the number of days is limited to M, the query object is limited to the communication database TEL constructed in the first step, the query conditions are TEL.IMSI = J10.IMSI and T ≤ TEL.Date < T + M, and the query fields are IMSI, Area, Call_IMSI, Time, Date, resulting in the total database of personnel communication characteristics J11.

8. The method for tracing personnel trajectories based on communication data according to claim 6, characterized in that: In the seventh step, the method for obtaining set A is as follows: the query object is limited to the personnel communication feature database J11, the query condition is limited to J11.IMSI = J3.IMSI, and J11.IMSI, J11.Area, J11.Call_IMSI, J11.Time, and J11.Date are extracted and denoted as set A. The method to obtain set B is as follows: limit the query condition to J11.IMSI! = J3.IMSI, extract J11.IMSI, J11.Area, J11.Call_IMSI, J11.Time, and J11.Date, and denote them as set B.

9. A method for tracing personnel trajectories based on communication data according to claim 8, characterized in that: In the seventh step, the method to obtain set C is as follows: combine set A with set B to obtain set C, with the conditions being A.Call_IMSI = B.Call_IMSI and A.Area = B.Area.

10. A method for tracing personnel trajectories based on communication data according to claim 8, characterized in that: In the seventh step, the method for extracting signaling data in the communication database whose user identification number is the same as that in set C within a set number of days is as follows: the query time start point is limited to date T, the number of days is limited to M, the query object is limited to the communication database COM, and the query conditions are COM.IMSI = C.IMSI and T ≤ COM.Date < T + M. COM.IMSI, COM.Area, COM.Call_IMSI, COM.Time, and COM.Date are extracted.

Images