An ecological protection method for identifying poaching based on IMSI and location tracking

Abstract

The present application relates to the technical field of ecological protection, and discloses an ecological protection method for identifying poaching based on IMSI and location tracking, which comprises the following steps: deploying a mobile terminal tracking device to realize identity information collection of IMSI in a covered area; confirming IMSI with a staying behavior through a space-time algorithm and marking the IMSI as a staying IMSI; when the staying IMSI moves along a track, calculating the accompanying IMSI at the current moving time of the staying IMSI and the accompanying IMSI before the staying through the space-time algorithm; marking the accompanying IMSI with a secondary accompanying result with the staying IMSI as a poaching device; marking the staying IMSI with the poaching device as a suspect object; comparing the accompanying and separating relationships of the poaching device and the suspect object to confirm the poaching behavior. The above technical solution compares the accompanying and separating relationships between the IMSI information of the animal catching clamp and the mobile terminal carried by the poacher, effectively identifies the poaching behavior and locks the illegal personnel, and improves the intelligent and accurate level of anti-poaching monitoring.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of ecological protection, and in particular to an ecological protection method for identifying poaching based on IMSI and location tracking. BACKGROUND

[0002] In the field of wild animal protection, illegal poaching behavior seriously threatens the survival of endangered species. Traditional monitoring methods mainly rely on devices such as cameras and infrared cameras for deployment and control, but have many limitations. For example, camera monitoring has limited coverage and cannot effectively monitor a large area. It is also greatly affected by weather, light, and terrain, and its monitoring effect is significantly reduced at night or in bad weather. In addition, traditional monitoring devices cannot actively identify and track the identity of poachers, and can only record the event process, lacking the ability to immediately locate and identify the identity of the violator.

[0003] Chinese patent document CN113486810B discloses an "intelligent identification method for park poaching birds". The steps include: step one, using a photographing device in the park to take pictures of people and trees inside the park; step two, using a tree detection module to process the obtained photos, then using the identification function in the tree detection module to identify the trees in the image, and determining whether there are birds on the trees in the image, and the tree detection module also identifies whether there are humans around the trees; step three, using a personnel detection module to detect and process the obtained photos. The above technical solution uses image recognition, which has limited protection range and is limited by the on-site poaching scene. SUMMARY

[0004] The present application mainly solves the technical problems of the original technical solution, which has limited protection range and is limited by the on-site poaching scene. The present application provides an ecological protection method for identifying poaching based on IMSI and location tracking. By collecting the IMSI information of the SIM card-equipped animal trapping device used by poachers and the mobile terminal carried by them, and combining space-time trajectory analysis, the present application realizes modeling of the accompanying relationship of poachers and identity binding. By comparing the accompanying and separating relationship between the IMSI information of the animal trapping device and the mobile terminal carried by poachers, the present application effectively identifies poaching behavior and locks down illegal personnel, improving the intelligence and precision level of anti-poaching monitoring.

[0005] The above technical problems of the present application are mainly solved by the following technical solution: the present application includes the following steps:

[0006] Deploying a mobile terminal tracking device to realize the collection of identity information of IMSI in the covered area;

[0007] Using a space-time algorithm to confirm the IMSI that has a stay behavior and marking it as a stay IMSI;

[0008] When the trajectory of the stay IMSI moves, the space-time algorithm is used to calculate the accompanying IMSI when the stay IMSI moves currently and the accompanying IMSI before the stay IMSI moves respectively;

[0009] The accompanying IMSI which has the second accompanying result with the stay IMSI is marked as a poaching device;

[0010] The stay IMSI which has the poaching device is marked as a suspect object;

[0011] The accompanying and separation relationship of the poaching device and the suspect object is compared, and then the poaching behavior is confirmed.

[0012] By comparing the IMSI, location information and timestamp data of the trap and the poaching terminal, the space-time correlation between the two is analyzed. Before the poacher installs the trap, the space-time is used to determine that the two are in an accompanying relationship. When the trap is stationary and the poacher leaves, it is determined that the two are in a separation relationship. When the trap sends a notification, the poacher is close to the trap and takes away the trap, it can be determined that the two are in an accompanying relationship, and then the poaching behavior is confirmed. According to the changes of the IMSI and location information of the trap and the poaching terminal, the state of the trap is automatically detected.

[0013] As preferred, the mobile terminal tracking device includes a wireless communication base station and a phased array radar. After the wireless communication base station is started, a wireless cell coverage is formed. The mobile terminal tracking device is composed of two main components, the wireless communication base station and the phased array radar. After the wireless communication base station module is started, a wireless cell coverage of 90 degrees and a radius of 13 kilometers is formed.

[0014] As preferred, when the mobile terminal enters the coverage area of the tracking device, the wireless communication base station collects the identity information of the mobile terminal and controls the mobile terminal to continuously report the uplink signal through signaling.

[0015] As preferred, the phased array radar measures the direction of arrival of the uplink signal and calculates the distance to realize the confirmation of the position of the mobile terminal. After the position is confirmed, the wireless base station sends a rejection signaling to the mobile terminal, and the mobile terminal returns to the operator network.

[0016] As preferred, the wireless communication base station periodically updates the wireless cell parameters. After the update, the mobile terminal will be collected again. The periodic collection of the SIM card identity information of the mobile terminal realizes positioning and identification.

[0017] As preferred, the judgment of the stay behavior specifically includes:

[0018] The system periodically locates the collected IMSI and records the historical location information thereof;

[0019] For each historical trajectory of IMSI, extract N_samples consecutive location points;

[0020] If the N_samples points meet the stay condition, mark as "stay".

[0021] T_min_stay: minimum stay time threshold (unit: seconds), default setting is 300 seconds;

[0022] D_radius: stay area radius (unit: meters), default setting is 200 meters;

[0023] N_samples: the number of sampling points for stay judgment, usually take no less than 5 consecutive valid location points;

[0024] Δt_max_interval: the maximum allowed time interval between adjacent points (unit: seconds), default is 300 seconds;

[0025] Δd_max_distance: the maximum allowed displacement distance between adjacent points (unit: meters), default is 200 meters;

[0026] As a preferred, the stay condition includes: the time interval between all adjacent points is less than or equal to the adjacent point time interval threshold; the maximum enclosing circle radius formed by all points is less than or equal to the stay area radius; the total stay time is greater than or equal to the minimum stay time threshold.

[0027] As a preferred, the judgment of the accompanying IMSI in the current movement specifically includes:

[0028] When a "stay IMSI" starts to move, extract the trajectory data of the "stay IMSI" in the moving process;

[0029] For each other IMSI, judge whether it appears in the associated distance threshold range of the "stay IMSI" in the same time window;

[0030] If the trajectories of the two meet the minimum coexistence time and are close in space, mark the other IMSI as the "current accompanying IMSI" of the "stay IMSI";

[0031] All "current accompanying IMSIs" that meet the conditions form an accompanying set S1.

[0032] Δt_window: time window size, default ±60 seconds;

[0033] Δd_associate: associated distance threshold, default 100 meters;

[0034] T_min_assoc: minimum coexistence time, default 30 seconds.

[0035] As preferred, the judgment of the accompanying IMSI before the stay specifically includes:

[0036] Backtracking the trajectory data within the time window before the stay IMSI is marked;

[0037] For each other IMSI, judging whether it appears within the pre-marking distance threshold range of the pre-marking associated distance of the stay IMSI;

[0038] If the trajectories of both exist a pre-marking minimum coexistence time and are kept close in space, the other IMSI is marked as the pre-marking accompanying IMSI of the stay IMSI;

[0039] All the pre-marking accompanying IMSIs satisfying the conditions form an accompanying set S2.

[0040] Δt_before: the backtracking time window before the marking 300 seconds (5 minutes);

[0041] Δd_before: the pre-marking associated distance 150 meters;

[0042] T_min_before: the pre-marking minimum coexistence time 60 seconds.

[0043] As preferred, the intersection of the accompanying IMSI of the stay IMSI at the current movement and the accompanying IMSI before the stay is taken as the accompanying IMSI of the secondary accompanying result, and is marked as a poaching device, and the more the accompanying IMSI appears, the higher the suspicion level is.

[0044] For each IMSI in the final output result set containing all the IMSIs satisfying the secondary accompanying condition, the following operations are performed: the accompanying frequency between the IMSI and multiple stay IMSIs is counted; if an IMSI appears in the secondary accompanying set of different stay IMSIs for multiple times, the suspicion level of the IMSI is improved; finally, the IMSI is marked as a suspect IMSI according to a preset rule (such as appearing for more than 3 times, or the total accompanying time is more than 5 minutes, etc.).

[0045] The beneficial effects of the present application are: by collecting the IMSI information of the SIM card trap used by the poacher and the mobile terminal carried by the poacher, and combining the space-time trajectory analysis, the accompanying relationship modeling and identity binding of the poacher are realized, the accompanying and separation relationship between the IMSI information of the trap and the mobile terminal carried by the poacher is compared, the poaching behavior is effectively identified and the illegal personnel is locked, and the intelligent and accurate level of anti-poaching monitoring is improved. BRIEF DESCRIPTION OF DRAWINGS

[0046] Figure 1 is a flow chart of the present application. DETAILED DESCRIPTION

[0047] In order to make the objects, technical solutions and advantages of the present application clearer, the technical solutions of the present application are further specifically described below in conjunction with the embodiments and the accompanying drawings. It should be understood that the specific embodiments described herein are only the best mode of the present application, which are used to explain the present application and do not limit the protection scope of the present application. All other embodiments obtained by those skilled in the art without creative labor fall within the protection scope of the present application.

[0048] Before the example embodiments are discussed in more detail, it should be mentioned that some example embodiments are described as processes or methods depicted as flow charts. Although the flow charts describe the operations (or steps) as sequential processes, many of the operations (or steps) can be implemented in parallel, concurrently or simultaneously. In addition, the order of the operations can be rearranged. The processes can be terminated when their operations are completed, but can also have additional steps not included in the drawings; the processes can correspond to methods, functions, procedures, subroutines, subprograms, etc.

[0049] The technical solutions of the present application are further specifically described below in conjunction with the embodiments and the accompanying drawings. EMBODIMENT

[0050] An ecological protection method for identifying poaching based on IMSI and location tracking according to the present embodiment, as shown in Figure 1 includes the following steps:

[0051] S1. Deploy a mobile terminal tracking device, which includes a wireless communication base station and a phased array radar. After the wireless communication base station is started, a wireless cell coverage is formed. The deployment of the mobile terminal tracking device is composed of two main components, the wireless communication base station and the phased array radar. After the wireless communication base station module is started, a wireless cell coverage of 90 degrees and a radius of 13 kilometers is formed.

[0052] Identity information collection of IMSI in the coverage area is implemented. When the mobile terminal enters the tracking device coverage area, the wireless communication base station collects the identity information of the mobile terminal and controls the mobile terminal to continuously report the uplink signal through signaling. The phased array radar measures the direction of arrival of the uplink signal and calculates the distance to realize the confirmation of the mobile terminal position. After the position is confirmed, the wireless base station sends a rejection signaling to the mobile terminal, and the mobile terminal returns to the operator network.

[0053] Wireless communication base station periodically updates wireless cell parameters, after updating, mobile terminal will be collected again, periodically collect SIM card identity information of mobile terminal to realize positioning and identification.

[0054] S2. Confirm the IMSI with staying behavior by space-time algorithm, and mark it as staying IMSI.

[0055] The judgment of staying behavior specifically includes:

[0056] The system periodically locates the collected IMSI, and records its historical position information;

[0057] For the historical trajectory of each IMSI, extract N_samples consecutive position points;

[0058] If the N_samples points meet the staying condition, mark it as "staying".

[0059] T_min_stay: minimum staying time threshold (unit: second), default setting is 300 seconds;

[0060] D_radius: staying area radius (unit: meter), default setting is 200 meters;

[0061] N_samples: sampling point number for judging staying, usually take no less than 5 consecutive valid position points;

[0062] Δt_max_interval: maximum allowed time interval between adjacent two position points (unit: second), default is 300 seconds;

[0063] Δd_max_distance: maximum allowed displacement distance between adjacent two position points (unit: meter), default is 200 meters;

[0064] The staying condition includes: the time interval between all adjacent points is less than or equal to the adjacent point time interval threshold; the maximum enclosing circle radius formed by all points is less than or equal to the staying area radius; the total staying time is greater than or equal to the minimum staying time threshold.

[0065] S3. When the staying IMSI moves, calculate the accompanying IMSI at the current moving time and the accompanying IMSI before staying of the staying IMSI respectively by space-time algorithm.

[0066] The judgment of the accompanying IMSI at the current moving time specifically includes:

[0067] When a "staying IMSI" starts to move, extract the trajectory data of the "staying IMSI" in the moving process;

[0068] For each other IMSI, judge whether it appears in the range of the associated distance threshold of the "stay IMSI" in the same time window;

[0069] If the trajectories of the two exist and meet the minimum coexistence time, and are close in space, mark the other IMSI as the "current accompanying IMSI" of the "stay IMSI";

[0070] All the "current accompanying IMSIs" that meet the condition form an accompanying set S1.

[0071] Δt_window: time window size, default ±60 seconds;

[0072] Δd_associate: associated distance threshold, default 100 meters;

[0073] T_min_assoc: minimum coexistence time, default 30 seconds.

[0074] The judgment of the accompanying IMSI before staying specifically includes:

[0075] Trace back the trajectory data in the time window before the "stay IMSI" is marked;

[0076] For each other IMSI, judge whether it appears in the range of the associated distance threshold of the "stay IMSI" in the same time window;

[0077] If the trajectories of the two exist and meet the minimum coexistence time, and are close in space, mark the other IMSI as the "current accompanying IMSI" of the "stay IMSI";

[0078] All the "current accompanying IMSIs" that meet the condition form an accompanying set S2.

[0079] Δt_before: backtracking time window before marking 300 seconds (5 minutes);

[0080] Δd_before: associated distance before marking 150 meters;

[0081] T_min_before: minimum coexistence time before marking 60 seconds.

[0082] S4. Mark the accompanying IMSI that has a secondary accompanying result with the stay IMSI as a poaching device. Take the intersection of the accompanying IMSI when the stay IMSI is currently moving and the accompanying IMSI before staying as the accompanying IMSI that has a secondary accompanying result, and mark it as a poaching device. The more the accompanying IMSI appears, the higher the suspicion level.

[0083] S5. Mark the stay IMSI existing poaching device as a suspect. Perform the following operations on each IMSI in the final output result set containing all the IMSIs satisfying the secondary accompanying condition: count the accompanying frequency between the IMSI and multiple "stay IMSIs"; if an IMSI appears in the secondary accompanying set of different "stay IMSIs" multiple times, increase its suspicion level; and finally mark it as a "suspect" according to a preset rule (such as the number of occurrences ≥ 3 times, or the total sum of accompanying time ≥ 5 minutes, etc.).

[0084] S6. Compare the accompanying and separation relationship between the poaching device and the suspect, and then confirm the poaching behavior.

[0085] By comparing the IMSI, location information and timestamp data of the trap and the poacher's mobile terminal, the spatio-temporal correlation between the two is analyzed. Before the poacher installs the trap, the spatio-temporal relationship between the two is determined. When the trap is installed and stationary, the poacher leaves, and the separation relationship between the two is determined. When the trap sends a notification, the poacher approaches the trap and removes the trap, and the accompanying relationship between the two is determined, and the poaching behavior is confirmed. According to the changes of the IMSI and location information of the trap and the poacher's mobile terminal, the state of the trap is automatically detected. Embodiment

[0086] A specific embodiment of an ecological protection method for identifying poaching based on IMSI and location tracking:

[0087] 1. Deploy a mobile terminal tracking device, which consists of two main components: a wireless communication base station and a phased array radar. After the wireless communication base station module is started, it forms a wireless cell coverage of 90 degrees and a radius of 13 kilometers. When a mobile terminal enters the coverage area, the wireless communication base station collects the identity information of the mobile terminal and controls it to continuously report uplink signals through signaling. Then, the phased array radar measures the direction of arrival of the uplink signal and calculates the distance to achieve the confirmation of the mobile terminal's location. After the location is confirmed, the wireless base station sends a rejection signaling to the mobile terminal, and the mobile terminal returns to the operator network. The wireless communication base station can periodically change its cell parameters. After the cell parameters are updated, the mobile terminal will be collected again.

[0088] 2. According to the collected data, the IMSI with stay behavior is confirmed through a spatio-temporal algorithm and is marked.

[0089] 3. When the IMSI marked with stay behavior moves along the trajectory, the accompanying IMSI at the current moving time is calculated through a spatio-temporal algorithm.

[0090] 4. When the IMSI marked with stay behavior moves along the trajectory, the accompanying IMSI before being marked with stay behavior is calculated through a spatio-temporal algorithm.

[0091] 5. Perform an intersection operation on the result sets of steps 3 and 4 to determine the IMSIs that have a secondary association with the IMSI marked as the stationary IMSI.

[0092] 6. Mark the secondary associated IMSI result calculated in step 5 as a suspect.

[0093] I. Criteria for Determining Stay Behavior

[0094] "Stay" refers to the behavior of a mobile terminal (IMSI) remaining in a geographical area for more than a set time threshold. To accurately identify stay behavior, the system defines the following stay determination criteria:

[0095] 1. Dwell Time Determination Parameters:

[0096] T_min_stay: Minimum stay time threshold (unit: seconds), default setting is 300 seconds.

[0097] D_radius: Radius of the dwell area (unit: meters), the default setting is 200 meters.

[0098] N_samples: The number of sampling points used to determine the dwell time, usually no less than 5 consecutive valid location points.

[0099] Δt_max_interval: The maximum allowed time interval (in seconds) between two adjacent locations, with a default value of 300 seconds.

[0100] Δd_max_distance: The maximum allowable displacement distance between two adjacent locations (in meters), the default is 200 meters.

[0101] 2. Decision logic flow:

[0102] The system periodically locates the collected IMSIs and records their historical location information (including latitude, longitude, and timestamp).

[0103] For each IMSI's historical trajectory, extract N_samples consecutive location points;

[0104] If these N_samples points satisfy the following condition, they are marked as "staying":

[0105] The time interval between all adjacent points is less than or equal to Δt_max_interval;

[0106] The radius of the largest enclosing circle formed by all points is less than or equal to D_radius;

[0107] Total stay time ≥ T_min_stay;

[0108] IMSI satisfying the above conditions will be marked as "stay IMSI" and enter the subsequent associated analysis process.

[0109] II. Spatiotemporal algorithm description

[0110] Spatiotemporal algorithm is used to identify other IMSIs that have spatial and temporal correlation with "stay IMSI". It is divided into two scenarios: current moving associated IMSI identification and pre-marking associated IMSI identification.

[0111] (I) Current moving associated IMSI identification algorithm

[0112] Parameter definition

[0113] Δt_window: Time window size, default ±60 seconds.

[0114] Δd_associate: Association distance threshold, default 100 meters.

[0115] T_min_assoc: Minimum coexistence time, default 30 seconds.

[0116] Algorithm flow:

[0117] 1. When a "stay IMSI" starts to move, the system extracts the trajectory data of the IMSI during the moving process;

[0118] 2. For each other IMSI, determine whether it appears near the "stay IMSI" (distance ≤ Δd_associate) in the same time period (Δt_window);

[0119] 3. If the trajectories of the two have at least T_min_assoc coexistence time and are close in space, mark the IMSI as "current associated IMSI";

[0120] 4. All IMSIs satisfying the conditions form an associated set S1.

[0121] Δt_before: Marking back time window, 300 seconds (5 minutes).

[0122] Δd_before; Marking before association distance, 150 meters.

[0123] T_min_before; Minimum coexistence time before marking, default 60 seconds.

[0124] (II) Pre-marking associated IMSI identification algorithm flow:

[0125] 1. Backtrack the trajectory of "stay IMSI" within Δt_before time before being marked;

[0126] 2. In this time period, check if there is any spatial overlap (distance ≤ Δd_before) with other IMSIs;

[0127] 3. If both trajectories coexist for at least T_min_before and remain close in space, mark this IMSI as "pre-marked companion IMSI";

[0128] 4. All IMSIs that satisfy the condition form another companion set S2.

[0129] Three, secondary companion determination mechanism

[0130] Operation flow:

[0131] 1. Perform intersection operation between the current companion IMSI set S1 obtained in (I) and the pre-marked companion IMSI set S2 obtained in (II);

[0132] 2. The result obtained is the IMSI that simultaneously satisfies the spatial and temporal association with "stay IMSI" before and after marking;

[0133] 3. These IMSIs are called "secondary companion IMSIs".

[0134] Parameter description:

[0135] Intersection (S1, S2): represents the intersection operation of sets S1 and S2;

[0136] Output_set: the final output result set, containing all IMSIs that satisfy the secondary companion condition.

[0137] Four, suspect object marking

[0138] Perform the following operations on each IMSI in Output_set:

[0139] 1. Count the number of times it appears in the secondary companion set of different "stay IMSIs";

[0140] 2. If an IMSI appears multiple times in the secondary companion set of different "stay IMSIs", increase its suspicion level;

[0141] 3. Finally, according to the preset rules (such as appearing ≥ 3 times, or the total time of companionship ≥ 5 minutes, etc.), mark it as "suspect IMSI".

[0142] The specific embodiments described herein are merely illustrative of the principles of this application. These embodiments are presented by way of example only, and are not intended to limit the scope of the application. Numerous alternatives to the embodiments described herein will be readily apparent to those skilled in the art. The principles of the application described herein can be employed in any of various embodiments or examples not expressly described herein. Accordingly, the scope of the application is to be construed as not being limited to the embodiments described herein, but rather is to be accorded the full scope consistent with the language of the claims, and any equivalents thereof.

Claims

1. An ecological protection method for identifying poaching based on IMSI and location tracking, characterized in that, The method comprises the following steps: deploying a mobile terminal tracking device to collect the identity information of IMSIs in the coverage area; confirming the IMSIs with staying behavior through a space-time algorithm and marking them as staying IMSIs; when the staying IMSIs move, calculating the accompanying IMSIs at the current time and the accompanying IMSIs before the staying through a space-time algorithm; marking the accompanying IMSIs with secondary accompanying results with the staying IMSIs as poaching devices; marking the staying IMSIs with poaching devices as suspects; comparing the accompanying and separating relationships of the poaching devices and the suspects to confirm the poaching behavior.

2. The ecological protection method for identifying poaching based on IMSI and location tracking according to claim 1, characterized in that, The mobile terminal tracking device comprises a wireless communication base station and a phased array radar. After the wireless communication base station is started, a wireless cell coverage is formed.

3. The ecological protection method for identifying poaching based on IMSI and location tracking according to claim 2, characterized in that, When a mobile terminal enters the coverage area of the tracking device, the wireless communication base station collects the identity information of the mobile terminal and controls the mobile terminal to continuously report uplink signals through signaling.

4. The ecological protection method for identifying poaching based on IMSI and location tracking according to claim 3, characterized in that, The phased array radar measures the direction of arrival of the uplink signals and calculates the distance to confirm the position of the mobile terminal. After the position is confirmed, the wireless base station sends a rejection signaling to the mobile terminal, and the mobile terminal returns to the operator network.

5. The method of claim 2, wherein the method further comprises: determining whether the IMSI is associated with a location tracking device. The wireless communication base station periodically updates the wireless cell parameters. After the update, the mobile terminal will be collected again. The periodic collection of the SIM card identity information of the mobile terminal realizes positioning and identification.

6. The ecological protection method for identifying poaching based on IMSI and location tracking according to claim 1, characterized in that, The judgment of the staying behavior specifically comprises: the system periodically locates the collected IMSIs and records their historical position information; for the historical trajectory of each IMSI, N_samples consecutive position points are extracted; if the N_samples points meet the staying condition, they are marked as "staying".

7. The ecological protection method for identifying poaching based on IMSI and location tracking according to claim 6, characterized in that, The staying condition comprises: the time interval between all adjacent points is less than or equal to the adjacent point time interval threshold; the maximum enclosing circle radius formed by all points is less than or equal to the staying area radius; the total staying time is greater than or equal to the minimum staying time threshold.

8. The ecological protection method for identifying poaching based on IMSI and location tracking according to claim 1, characterized in that, The judgment of the accompanying IMSIs at the current time specifically comprises: when a "staying IMSI" starts to move, the trajectory data of the "staying IMSI" during the movement is extracted; for each other IMSI, it is judged whether it appears within the associated distance threshold range of the "staying IMSI" in the same time window; if the trajectories of the two exist and meet the minimum coexistence time, and are close in space, the other IMSI is marked as the "current accompanying IMSI" of the "staying IMSI"; all the "current accompanying IMSIs" that meet the condition form an accompanying set S1.

9. The ecological protection method for identifying poaching based on IMSI and location tracking according to claim 1, characterized in that, The judgment of the accompanying IMSIs before the staying specifically comprises: backtracking the trajectory data in the time window before the "staying IMSI" is marked; for each other IMSI, it is judged whether it appears within the associated distance threshold range of the "staying IMSI" before being marked in the time window before being marked; if the trajectories of the two exist and meet the minimum coexistence time before being marked, and are close in space, the other IMSI is marked as the "pre-marking accompanying IMSI" of the "staying IMSI". All the "tagged pre-IMSI" that meet the conditions form a set of accompanying S2.

10. The ecological protection method for identifying poaching based on IMSI and location tracking according to claim 1 or 8 or 9, characterized in that, Take the intersection of the accompanying IMSI when the current mobile stays and the accompanying IMSI before staying as the accompanying IMSI of the secondary accompanying result, and mark it as a stolen device. The more the accompanying IMSI appears, the higher the suspicion level.

Images