Method for constructing electromagnetic interference simulation model for shared iron tower

Abstract

The invention relates to the technical field of data processing, and provides a shared iron tower-oriented electromagnetic interference simulation model construction method, which comprises the steps of establishing a three-dimensional terrain model of a to-be-constructed region, and collecting electromagnetic data of all sampling positions of the three-dimensional terrain model; determining geographic information complexity and sampling selection goodness of sampling positions, and screening optimal sampling positions; acquiring an equipotential line and an electromagnetic fluctuation degree of the optimal sampling position; and according to the electromagnetic data, the geographic information complexity and the electromagnetic fluctuation of all the optimal sampling positions, dividing the region to be constructed into different grid regions, and obtaining an electromagnetic interference simulation model of the shared iron tower. The invention aims to adaptively divide the grid region and improve the quality of the electromagnetic interference simulation model.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of data processing, in particular to a method for constructing an electromagnetic interference simulation model for a shared tower. BACKGROUND

[0002] A shared tower refers to a tower on which optical cables, communication base stations, mobile antennas and other communication equipment are installed to realize the co-construction and sharing of power and communication infrastructure. The various communication equipment installed on the power tower is susceptible to mutual influence and may produce electromagnetic interference, resulting in signal attenuation, increased bit error rate and other problems. At the same time, the terrain will affect the propagation path of electromagnetic waves and also change the intensity, reflection, refraction and other behaviors of the waves. Therefore, when establishing an electromagnetic interference simulation model for a shared tower, the terrain and electromagnetic data at different locations need to be considered. Further, in order to improve the accuracy of the electromagnetic interference simulation model, the electric field region needs to be divided into different grid regions during electromagnetic signal simulation.

[0003] When dividing the grid region, the size of the grid region will affect the quality and simulation speed of the electromagnetic interference simulation model. A larger grid region can improve the calculation efficiency, but it may affect the authenticity of the simulation. A smaller grid region can improve the simulation effect, but the dramatic increase in data volume may affect the normal operation of the equipment. Therefore, the specific division method of the grid region needs to be determined according to the specific conditions of the terrain and electromagnetic data at different locations in the electric field region during electromagnetic signal simulation. SUMMARY

[0004] The present application provides a method for constructing an electromagnetic interference simulation model for a shared tower to solve the problem of unreasonable division of grid regions affecting the quality and simulation speed of the electromagnetic interference simulation model. The technical solution adopted is as follows: One embodiment of the present application provides a method for constructing an electromagnetic interference simulation model for a shared tower, which includes the following steps: establishing a three-dimensional terrain model of the to-be-constructed region of the electromagnetic interference simulation model, and collecting electromagnetic data at all sampling locations of the three-dimensional terrain model multiple times; determining the geographical information complexity of each sampling location according to the difference between the heights of adjacent sampling locations and the difference in height between each sampling location and the shared tower, determining the sampling selection preference of each sampling location according to the geographical information complexity of all sampling locations between the sampling location and the shared tower and the distance between different sampling locations, and selecting all preferred sampling locations according to the sampling selection preference; According to all the electromagnetic data collected at all the preferred sampling positions in the three-dimensional terrain model, equipotential lines of each preferred sampling position are fitted respectively, and the electromagnetic fluctuation degree of each preferred sampling position is calculated according to the differences between the equipotential lines of different preferred sampling positions, the distances between different preferred sampling positions, and the differences between all the electromagnetic data collected at the same preferred sampling position. According to the electromagnetic data, the geographical information complexity, and the electromagnetic fluctuation degree of all the preferred sampling positions, the region to be constructed is divided into different grid regions, and the electromagnetic interference simulation model of the shared tower is obtained by processing all the grid regions divided from the region to be constructed, the electromagnetic data of all the sampling positions in all the grid regions, and the three-dimensional terrain model of the region to be constructed using electromagnetic wave propagation simulation software.

[0005] Further, the method for obtaining the geographical information complexity of the sampling position is as follows: Any one sampling position is recorded as a target sampling position, all the sampling positions in a window centered on the target sampling position are recorded as the neighboring sampling positions of the target sampling position, the standard deviation of the heights of all the neighboring sampling positions of the target sampling position is recorded as the first standard deviation of the target sampling position. The absolute value of the difference between the height of the target sampling position and the height of the shared tower is recorded as the first absolute value of the target sampling position. The geographical information complexity of the target sampling position is determined according to the first standard deviation and the first absolute value of the target sampling position.

[0006] Further, the method for determining the geographical information complexity of the target sampling position according to the first standard deviation and the first absolute value of the target sampling position includes the following specific method: The sum of the first standard deviation and the first absolute value of the target sampling position is recorded as the geographical information complexity of the target sampling position.

[0007] Further, the method for determining the sampling selection priority of the sampling position is as follows: The cumulative sum of the geographical information complexity of all the sampling positions contained in the line segment determined by the target sampling position and the shared tower as the end points is recorded as the first cumulative sum of the target sampling position. The distance between the target sampling position and the nearest sampling position is recorded as the first distance of the target sampling position. The ratio of the product of the geographical information complexity of the target sampling position and the first distance to the first cumulative sum is recorded as the sampling selection priority of the target sampling position.

[0008] Further, the method for screening the preferred sampling position is as follows: The sampling positions with the sampling selection priority greater than the preset preferred threshold value are all recorded as the preferred sampling positions.​

[0009] Further, the method for obtaining the equipotential line of the preferred sampling position is: The average electromagnetic data of the preferred sampling position is recorded as the average electromagnetic data of the preferred sampling position. Any preferred sampling position is recorded as a target preferred sampling position. The average electromagnetic data of all preferred sampling positions is used to interpolate the electromagnetic data of different positions to obtain the electromagnetic data interpolation of all positions equal to the average electromagnetic data of the target preferred sampling position. The average electromagnetic data of the target preferred sampling position and the electromagnetic data interpolation of all positions equal to the average electromagnetic data of the target preferred sampling position are used for curve fitting to obtain the equipotential line of the target preferred sampling position.

[0010] Further, the electromagnetic fluctuation degree of each preferred sampling position is calculated according to the difference between the equipotential lines of different preferred sampling positions, the distance between different preferred sampling positions, and the difference between all electromagnetic data collected at the same preferred sampling position, including the specific method: The adjacent equipotential line difference of the target preferred sampling position is determined according to the difference between the equipotential lines of different preferred sampling positions. The distance between the two preferred sampling positions corresponding to the adjacent equipotential line difference of the target preferred sampling position is recorded as the adjacent equipotential line distance of the target preferred sampling position. The ratio of the adjacent equipotential line difference of the target preferred sampling position to the adjacent equipotential line distance is recorded as the plane affected degree of the target preferred sampling position. The normalized value of the product of the variance of all electromagnetic data collected at the target preferred sampling position and the plane affected degree of the target preferred sampling position is recorded as the electromagnetic fluctuation degree of the target preferred sampling position.

[0011] Further, the adjacent equipotential line difference of the target preferred sampling position is determined according to the difference between the equipotential lines of different preferred sampling positions, including the specific method: The Frechet distance between the equipotential line of the target preferred sampling position and the equipotential line of each preferred sampling position is calculated respectively. The minimum value of the Frechet distance between the equipotential line of the target preferred sampling position is recorded as the adjacent equipotential line difference of the target preferred sampling position.

[0012] Further, the method for dividing the region to be constructed into different grid regions according to the electromagnetic data, geographical information complexity, and electromagnetic fluctuation degree of all preferred sampling positions includes the specific method: The wavelength corresponding to the electromagnetic data at the same preferred sampling position is taken as the x-axis value, the frequency corresponding to the electromagnetic data at the same preferred sampling position is taken as the y-axis value, and the electromagnetic fluctuation at the preferred sampling position is taken as the z-axis value. Surface fitting is performed on the data points corresponding to different frequencies and wavelengths of the electromagnetic data at the same preferred sampling position to obtain the fitted surface of the same preferred sampling position. Based on the electromagnetic fluctuation degree and geographical information complexity of all preferred sampling locations, calculate the stopping feature value corresponding to the target preferred sampling location when it extends to different numbers of preferred sampling locations in different directions, until the stopping feature value is greater than or equal to the preset stopping threshold, and then stop calculating the stopping feature value in the corresponding direction. The preferred sampling positions corresponding to all stopping feature values ​​less than the stopping threshold corresponding to the target preferred sampling position are denoted as the sampling positions of the same sub-region; the convex hull of all preferred sampling positions of the same sub-region sampling positions is denoted as the grid region; and the region composed of different grid regions that overlap is considered as a single grid region. All regions in the area to be constructed that have not been assigned to a grid region are denoted as grid regions.

[0013] Furthermore, the formula for calculating the stopping feature value is: in, Indicates the first The preferred sampling position is oriented in the direction extend The stopping feature value corresponding to the preferred sampling position; Indicates the first The preferred sampling position is oriented in the direction The extended first The geographical information complexity of each preferred sampling location; Indicates the first The preferred sampling position is oriented in the direction The extended first Electromagnetic fluctuations at preferred sampling locations.

[0014] The beneficial effects of this application are: This application first considers that when electromagnetic information is affected by electromagnetic interference and ground-type interference, the electromagnetic information at adjacent locations is quite complex. When constructing an electromagnetic interference simulation model based on this electromagnetic information, the electromagnetic information at different sampling locations is relatively dense. The complexity of the electromagnetic information at each sampling location is evaluated, the geographical information complexity of the sampling location is obtained, and the sampling location's selection quality is obtained in conjunction with the sampling location's effectiveness. Based on the sampling selection quality, all preferred sampling locations are selected, i.e., reasonable sampling points are selected for constructing the electromagnetic interference simulation model. Furthermore, considering that signal transmission on shared towers is radial, and that signals are affected by buildings and ground-type interference during transmission... When interference occurs due to various reasons such as shape influence, transmitted signals at equal distances from the source within the same two-dimensional plane will exhibit differences. Equipotential lines are constructed for each preferred sampling location, and the electromagnetic fluctuation degree at each preferred sampling location is calculated based on these lines. Finally, based on the electromagnetic data, geographical information complexity, and electromagnetic fluctuation degree of all preferred sampling locations, the area to be constructed is adaptively divided into grid regions. This addresses the problem of unreasonable grid region division affecting the quality and simulation speed of the electromagnetic interference simulation model. Based on the divided grid regions, electromagnetic wave propagation simulation software is used to obtain the electromagnetic interference simulation model of the shared tower, thereby improving the quality of the electromagnetic interference simulation model. Attached Figure Description

[0015] To more clearly illustrate the technical solutions in the embodiments of this application or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0016] Figure 1 This is a schematic flowchart illustrating a method for constructing an electromagnetic interference simulation model for shared iron towers, provided in one embodiment of this application. Figure 2 This is a schematic diagram of the equipotential lines of a preferred sampling location provided in one embodiment of this application. Detailed Implementation

[0017] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0018] Please see Figure 1The diagram illustrates a flowchart of a method for constructing an electromagnetic interference simulation model for shared iron towers, according to an embodiment of this application. The method includes the following steps: Step S001: Establish a three-dimensional terrain model of the area to be constructed for the electromagnetic interference simulation model, and collect electromagnetic data from all sampling locations of the three-dimensional terrain model multiple times.

[0019] The DEM data of the area to be constructed in the electromagnetic interference simulation model was acquired using remote sensing satellites. A LiDAR scanner mounted on a drone was used to acquire LiDAR point cloud data of the area. A 3D terrain model was then constructed using the DEM and LiDAR point cloud data.

[0020] Among them, DEM data, or Digital Elevation Model data, is usually raster data, while each point in LiDAR point cloud data has corresponding three-dimensional coordinates. Constructing a three-dimensional terrain model using DEM data and LiDAR point cloud data is a well-known technique and will not be elaborated further.

[0021] Testers used a handheld signal receiver to acquire electromagnetic data at different sampling locations on a 3D terrain model, and also acquired electromagnetic data multiple times at the same sampling location.

[0022] Preferably, in one embodiment of this application, electromagnetic data is acquired five times at the same sampling location. In practical applications, as other implementation methods, the implementer can decide the value of the number of times electromagnetic data is acquired at the same sampling location according to the actual situation. This application does not impose any special restrictions.

[0023] Thus, electromagnetic data from multiple sampling locations on the 3D terrain model were obtained.

[0024] Step S002: Based on the height difference between adjacent sampling locations and the height difference between each sampling location and the shared tower, determine the geographical information complexity of each sampling location. Based on the geographical information complexity of all sampling locations between the sampling location and the shared tower, and the distance between different sampling locations, determine the sampling selection quality of each sampling location. Based on the sampling selection quality, select all preferred sampling locations.

[0025] The various communication devices installed on the tower are susceptible to mutual interference, potentially causing electromagnetic interference that leads to signal attenuation and increased bit error rate. When the electromagnetic wave source is located at a low altitude, obstacles such as tall buildings and mountains can cause shadows to appear on the waves, while when the electromagnetic wave source is located at a high altitude, the electromagnetic waves can propagate further. Different types of ground surfaces, such as water bodies, forests, and cities, have different characteristics in terms of reflection, diffraction, and absorption. For example, water bodies reflect electromagnetic waves strongly, while forests attenuate signals significantly. Therefore, different types of ground surfaces also affect the propagation of electromagnetic waves. Thus, terrain affects the propagation path of electromagnetic waves and can also change the wave's intensity, reflection, and refraction. When establishing an electromagnetic interference simulation model for a shared tower, it is necessary to consider both the terrain and electromagnetic data from different locations simultaneously.

[0026] When electromagnetic information is affected by electromagnetic interference and ground-type interference, the electromagnetic information at adjacent locations is quite complex. Therefore, when constructing an electromagnetic interference simulation model based on this electromagnetic information, the electromagnetic information at different sampling locations should also be relatively dense. At the same time, it is also necessary to consider the effectiveness of each sampling location and select reasonable sampling points for constructing the electromagnetic interference simulation model.

[0027] FDTD (Finite-Difference Time-Domain) simulation is a numerical simulation method used to calculate the distribution and propagation of electromagnetic fields in the time domain. The sampling locations used in FDTD simulations are the appropriate sampling points for constructing electromagnetic interference simulation models. The setting of the sampling locations in FDTD simulations has a significant impact on the effectiveness of the established electromagnetic interference simulation model; therefore, it is necessary to set the sampling locations appropriately. Ideally, the sampling locations in FDTD simulations should be located at key locations in the region to be constructed, such as near the wave source, the center of the region of interest, or areas where the field distribution changes drastically.

[0028] The geographic information complexity of each sampling location is determined based on the height difference between adjacent sampling locations and the height difference between each sampling location and the shared tower.

[0029] Let any sampling location be designated as the target sampling location, and then define the sampling location as the center of the target sampling location. All sampling locations within the window are denoted as neighboring sampling locations of the target sampling location. The standard deviation of the height of all neighboring sampling locations of the target sampling location is denoted as the first standard deviation of the target sampling location. The absolute value of the difference between the height of the target sampling location and the height of the shared tower is denoted as the first absolute value of the target sampling location. The sum of the first standard deviation and the first absolute value of the target sampling location is denoted as the geographic information complexity of the target sampling location.

[0030] The same method can be used to obtain the geographic information complexity of each sampling location in the 3D terrain model.

[0031] The greater the geographical complexity of the sampling location, the more complex the electromagnetic information at the sampling location, and the denser the electromagnetic information provided by the sampling location. At the same time, the greater the likelihood that the sampling location is located at a lower position, the greater the likelihood of a wave shadow area appearing when the sampling location is used as a reasonable sampling point for FDTD simulation.

[0032] Based on the geographical information complexity of all sampling locations between the sampling location and the shared tower, and the distance between different sampling locations, the sampling selection quality of each sampling location is determined.

[0033] The sum of the geographic information complexity of all sampling locations included in the line segment determined by the target sampling location and the shared tower as endpoints is denoted as the first sum of the target sampling location; the distance between the target sampling location and the nearest sampling location is denoted as the first distance of the target sampling location; the ratio of the product of the geographic information complexity of the target sampling location and the first distance to the first sum is denoted as the sampling selection quality of the target sampling location.

[0034] The first accumulated sum of the target sampling locations is used to evaluate the complexity of the electromagnetic information at each sampling location along the path from the target sampling location to the shared tower. The greater the complexity of the electromagnetic information at each sampling location along the path, the higher the sampling quality of the target sampling location. When selecting the target sampling location as a reasonable sampling point for FDTD simulation, more dense electromagnetic information can be extracted. In calculating the sampling quality of the target sampling location, the first distance of the target sampling location is used to avoid setting the reasonable sampling points for FDTD simulation to be too dense.

[0035] The same method can be used to obtain the sampling quality of each sampling location in the 3D terrain model.

[0036] Sampling positions with a sampling quality greater than the preset preferred threshold are all recorded as preferred sampling positions.

[0037] It is understandable that the preferred sampling location is the reasonable sampling point for FDTD simulation, which is also the reasonable sampling point used to construct the electromagnetic interference simulation model.

[0038] At this point, all preferred sampling locations have been selected.

[0039] Step S003: Based on all electromagnetic data collected from all preferred sampling locations in the three-dimensional terrain model, fit the equipotential line of each preferred sampling location. Calculate the electromagnetic fluctuation degree of each preferred sampling location based on the differences between equipotential lines at different preferred sampling locations, the distance between different preferred sampling locations, and the differences between all electromagnetic data collected at the same preferred sampling location.

[0040] Signal transmission on shared towers is radial. Transmitted signals at equal distances from the source within the same two-dimensional plane should be consistent. However, during transmission, signals are affected by various factors such as buildings and terrain, resulting in differences in transmitted signals at equal distances from the source within the same two-dimensional plane. Therefore, equipotential lines can be constructed to determine the complexity of interference affecting the electromagnetic signal.

[0041] Based on all electromagnetic data collected from all preferred sampling locations in the 3D terrain model, fit the equipotential lines for each preferred sampling location.

[0042] Preferably, as an embodiment of this application, the average value of all electromagnetic data collected at the same preferred sampling location of the three-dimensional terrain model is recorded as the average electromagnetic data of the same preferred sampling location. Any preferred sampling location is designated as the target preferred sampling location. Based on the average electromagnetic data collected at all preferred sampling locations of the three-dimensional terrain model, interpolation calculations are performed on the electromagnetic data at different locations to obtain the electromagnetic data interpolation values ​​for all locations whose average electromagnetic data is equal to that of the target preferred sampling location. Curve fitting is then performed based on the average electromagnetic data of the target preferred sampling location and the electromagnetic data interpolation values ​​of the locations whose average electromagnetic data is equal to that of the target preferred sampling location to obtain the equipotential line of the target preferred sampling location.

[0043] Interpolating electromagnetic data and curve fitting data are well-known techniques and will not be elaborated further. Specifically, this embodiment uses cubic spline interpolation to interpolate electromagnetic data and least squares method to perform curve fitting.

[0044] When the electromagnetic data at each preferred sampling location is undisturbed, the equipotential lines at all preferred sampling locations are circular. However, when electromagnetic data is interfered with, the affected positions on the equipotential lines at the preferred sampling locations will be distorted. Due to interference from various communication devices installed on the power transmission tower and the type of ground, the equipotential lines at each preferred sampling location are distorted circles. A schematic diagram of the equipotential lines at the preferred sampling locations is shown in Figure 2. Figure 2 In the diagram, 1 represents the preferred sampling location, and 2 represents the equipotential line of the preferred sampling location.

[0045] Calculate the Frechet distance between the equipotential line of the target preferred sampling position and the equipotential line of each preferred sampling position. The minimum Frechet distance between the target preferred sampling position and the equipotential line of the target preferred sampling position is recorded as the difference between adjacent equipotential lines of the target preferred sampling position. The distance between the two preferred sampling positions corresponding to the difference between adjacent equipotential lines of the target preferred sampling position is recorded as the distance between adjacent equipotential lines of the target preferred sampling position. The ratio of the difference between adjacent equipotential lines of the target preferred sampling position to the distance between adjacent equipotential lines is recorded as the degree of planar influence of the target preferred sampling position.

[0046] Among them, the equipotential lines are curves, and the calculation of Frechet distance between different curves is a well-known technique, which will not be elaborated here.

[0047] The smaller the plane influence of the target preferred sampling location, the smaller the grid area should be divided in the electric field region where the target preferred sampling location is located, so as to improve the quality of the electromagnetic interference simulation model.

[0048] The normalized value of the product of the variance of all electromagnetic data collected at the target preferred sampling location and the plane influence degree of the target preferred sampling location is denoted as the electromagnetic fluctuation degree of the target preferred sampling location.

[0049] It should be noted that this embodiment uses the Z-Score standard normalization method to calculate the normalized value. In practical applications, implementers may use other methods of existing technology, such as the maximum-minimum normalization method or the sigmoid function, to calculate the normalized value, and no limitation is made here.

[0050] The electromagnetic fluctuations at each preferred sampling location of the three-dimensional terrain model can be obtained using the same method.

[0051] At this point, the electromagnetic fluctuations at all preferred sampling locations of the 3D terrain model have been obtained.

[0052] Step S004: Based on the electromagnetic data, geographic information complexity, and electromagnetic fluctuation of all preferred sampling locations, the area to be constructed is divided into different grid areas. Electromagnetic wave propagation simulation software is used to process all grid areas divided into the area to be constructed, the electromagnetic data of all sampling locations in all grid areas, and the three-dimensional terrain model of the area to be constructed to obtain the electromagnetic interference simulation model of the shared tower.

[0053] The wavelength corresponding to the electromagnetic data at the same preferred sampling location is taken as the x-axis value, the frequency corresponding to the electromagnetic data at the same preferred sampling location is taken as the y-axis value, and the electromagnetic fluctuation at the preferred sampling location is taken as the z-axis value. The least squares method is used to perform surface fitting on the data points corresponding to different frequencies and wavelengths of the electromagnetic data at the same preferred sampling location to obtain the fitted surface of the same preferred sampling location.

[0054] Surface fitting for different data points to obtain the fitted surface is a well-known technique and will not be elaborated further.

[0055] Based on the electromagnetic fluctuation degree and geographical information complexity of all preferred sampling locations, the stopping characteristic value corresponding to the extension of the target preferred sampling location to different numbers of preferred sampling locations in different directions is calculated. The calculation formula is as follows: in, Indicates the first The preferred sampling position is oriented in the direction extend The stopping feature value corresponding to the preferred sampling position; Indicates the first The preferred sampling position is oriented in the direction The extended first The geographical information complexity of each preferred sampling location; Indicates the first The preferred sampling position is oriented in the direction The extended first Electromagnetic fluctuations at preferred sampling locations.

[0056] in, The value starts from 1 and increases gradually in increments of 1 until the stopping feature value is greater than or equal to a preset stopping threshold; direction It can be any direction within the plane.

[0057] The stopping threshold is set to 0.75. The preferred sampling positions corresponding to all stopping feature values ​​less than the stopping threshold at the target preferred sampling position are denoted as the same sub-region sampling positions.

[0058] The convex hull of all preferred sampling locations that are determined to be sampling locations in the same sub-region is denoted as the grid region.

[0059] It is understandable that when different grid regions overlap, the region consisting of the overlapping different grid regions is considered as a single grid region.

[0060] All regions in the area to be constructed that have not been assigned to a grid region are denoted as grid regions.

[0061] At this point, the area to be constructed is divided into different grid areas.

[0062] Using electromagnetic wave propagation simulation software, an electromagnetic interference simulation model of the shared tower is obtained based on the electromagnetic data of all grid areas divided into the area to be constructed, the electromagnetic data of all sampling locations within all grid areas, and the three-dimensional terrain model of the area to be constructed.

[0063] Electromagnetic wave propagation simulation software can be used, such as CST Studio Suite and COMSOL Multiphysics. Obtaining electromagnetic interference simulation models through electromagnetic wave propagation simulation software is a well-known technique and will not be elaborated further.

[0064] This completes the construction of the electromagnetic interference simulation model for the shared tower.

[0065] The above description is only a preferred embodiment of this application and is not intended to limit this application. Any modifications, equivalent substitutions, improvements, etc., made within the principles of this application should be included within the protection scope of this application.

Claims

1. A method for constructing an electromagnetic interference simulation model for shared iron towers, characterized in that, The method includes the following steps: A three-dimensional terrain model of the area to be constructed is established for the electromagnetic interference simulation model, and electromagnetic data of all sampling locations of the three-dimensional terrain model are collected multiple times. Based on the height differences between adjacent sampling locations and the height differences between each sampling location and the shared tower, the geographic information complexity of each sampling location is determined. Any sampling location is designated as the target sampling location. The sum of the geographic information complexities of all sampling locations included in the line segment defined by the target sampling location and the shared tower as endpoints is designated as the first sum of the target sampling location. The distance between the target sampling location and the nearest sampling location is designated as the first distance of the target sampling location. The ratio of the product of the geographic information complexity of the target sampling location and the first distance to the first sum is designated as the sampling selection quality of the target sampling location. Based on the sampling selection quality, all preferred sampling locations are selected. Based on all electromagnetic data collected from all preferred sampling locations in the 3D terrain model, equipotential lines are fitted for each preferred sampling location. The differences between adjacent equipotential lines at different preferred sampling locations are used to determine the differences between adjacent equipotential lines at the target preferred sampling location. The distance between two preferred sampling locations corresponding to the differences between adjacent equipotential lines at the target preferred sampling location is denoted as the adjacent equipotential line distance of the target preferred sampling location. The ratio of the adjacent equipotential line difference to the adjacent equipotential line distance of the target preferred sampling location is denoted as the planar influence degree of the target preferred sampling location. The normalized value of the product of the variance of all electromagnetic data collected at the target preferred sampling location and the planar influence degree of the target preferred sampling location is denoted as the electromagnetic fluctuation degree of the target preferred sampling location. Using the wavelength corresponding to the electromagnetic data at the same preferred sampling location as the x-axis value, the frequency corresponding to the electromagnetic data at the same preferred sampling location as the y-axis value, and the electromagnetic fluctuation degree of the preferred sampling location as the z-axis value, surface fitting is performed on the data points corresponding to different frequencies and wavelengths of the electromagnetic data at the same preferred sampling location to obtain the fitted surface of the same preferred sampling location. Based on the electromagnetic fluctuation degree and geographical information complexity of all preferred sampling locations, the stopping feature value corresponding to the target preferred sampling location extended in different directions by different numbers of preferred sampling locations is calculated, until the stopping feature value is greater than or equal to a preset stopping threshold, at which point the calculation of the corresponding direction is stopped. Stop feature value; record the preferred sampling positions corresponding to all stop feature values ​​less than the stop threshold corresponding to the target preferred sampling position as the same sub-region sampling positions; record the convex hull of all preferred sampling positions of the same sub-region sampling positions as the grid region; record the region composed of different overlapping grid regions as a single grid region; record all regions in the region to be constructed that have not been divided into grid regions as grid regions; use electromagnetic wave propagation simulation software to process the electromagnetic data of all grid regions divided in the region to be constructed, all sampling positions in all grid regions, and the three-dimensional terrain model of the region to be constructed to obtain the electromagnetic interference simulation model of the shared tower.

2. The method for constructing an electromagnetic interference simulation model for shared towers according to claim 1, characterized in that, The method for obtaining the geographical information complexity of the sampling location is as follows: Centered on the target sampling location All sampling positions within the window are recorded as neighboring sampling positions of the target sampling position, and the standard deviation of the height of all neighboring sampling positions of the target sampling position is recorded as the first standard deviation of the target sampling position. The absolute value of the difference between the target sampling location and the height of the shared tower is recorded as the first absolute value of the target sampling location; The geographic information complexity of the target sampling location is determined based on the first standard deviation and the first absolute value of the target sampling location.

3. The method for constructing an electromagnetic interference simulation model for shared towers according to claim 2, characterized in that, The method for determining the geographic information complexity of the target sampling location based on the first standard deviation and the first absolute value of the target sampling location includes the following specific methods: The sum of the first standard deviation and the first absolute value of the target sampling location is denoted as the geographic information complexity of the target sampling location.

4. The method for constructing an electromagnetic interference simulation model for shared iron towers according to claim 1, characterized in that, The method for selecting the preferred sampling location is as follows: Sampling positions with a sampling quality greater than the preset preferred threshold are all recorded as preferred sampling positions.

5. The method for constructing an electromagnetic interference simulation model for shared iron towers according to claim 1, characterized in that, The method for obtaining the equipotential lines at the preferred sampling locations is as follows: The mean value of all electromagnetic data collected at the same preferred sampling location of the three-dimensional terrain model is denoted as the average electromagnetic data at the same preferred sampling location. Any preferred sampling location is designated as the target preferred sampling location. Based on the average electromagnetic data collected from all preferred sampling locations in the 3D terrain model, the electromagnetic data at different locations are interpolated to obtain the electromagnetic data interpolation values ​​of all locations with the same average electromagnetic data as the target preferred sampling location. Based on the average electromagnetic data of the target preferred sampling location and the electromagnetic data interpolation values ​​of the locations with the same average electromagnetic data as the target preferred sampling location, curve fitting is performed to obtain the equipotential line of the target preferred sampling location.

6. The method for constructing an electromagnetic interference simulation model for shared towers according to claim 1, characterized in that, The method for determining the difference between adjacent equipotential lines at the target preferred sampling position based on the difference between equipotential lines at different preferred sampling positions includes: Calculate the Frechet distance between the equipotential line of the target preferred sampling position and the equipotential line of each preferred sampling position. The minimum Frechet distance between the target preferred sampling position and the equipotential line of the target preferred sampling position is recorded as the difference between adjacent equipotential lines of the target preferred sampling position.

7. The method for constructing an electromagnetic interference simulation model for shared towers according to claim 1, characterized in that, The formula for calculating the stopping characteristic value is: in, Indicates the first The preferred sampling position is oriented in the direction extend The stopping feature value corresponding to the preferred sampling position; Indicates the first The preferred sampling position is oriented in the direction The extended first The geographical information complexity of each preferred sampling location; Indicates the first The preferred sampling position is oriented in the direction The extended first Electromagnetic fluctuations at preferred sampling locations.

Images