Summary of the invention
The purpose of the present invention is overcoming the above-mentioned prior art, provide it is a kind of meet accuracy is high, error is small,
The method and corresponding system that field strength distribution prediction processing is realized based on frequency spectrum drive test data applied widely.
To achieve the goals above, it is of the invention based on frequency spectrum drive test data realize field strength distribution prediction processing method and
Corresponding system is as follows:
The method for realizing field strength distribution prediction processing based on frequency spectrum drive test data, is mainly characterized by, the method
The following steps are included:
(1) the frequency spectrum drive test data of estimation range is obtained;
(2) rasterizing pretreatment is carried out to estimation range and frequency spectrum drive test data;
(3) pretreated frequency spectrum data is considered as sample point creation Thiessen polygon;
(4) according to the weight of each adjacent sample point of Natural neighbors interpolation calculation, the grid of unknown field intensity value is carried out slotting
Value processing;
(5) the field intensity prediction value of the grid is calculated according to the weight of each adjacent sample point of interpolation grid and field intensity value;
(6) it is divided according to color range and the field intensity value of whole region is converted into field intensity map.
Preferably, frequency spectrum drive test data includes latitude and longitude information, frequency information and field intensity value in the step (1).
Preferably, the step (2) specifically includes the following steps:
(2.1) grid division is carried out to estimation range according to designated size;
(2.2) frequency spectrum drive test data is filtered into cleaning, removes incomplete, mistake and duplicate data;
(2.3) longitude and latitude is located at a plurality of frequency spectrum drive test data arithmetic mean within the scope of a grid is a data,
And the longitude and latitude of the central point of grid where being set as the longitude and latitude of the data, grid where the field intensity value of the data represents
Field intensity value.
Preferably, pretreated frequency spectrum data is the frequency spectrum data after grid is average in the step (3), it is specific to wrap
Grid central point longitude and latitude, frequency information and corresponding field intensity value where including.
Preferably, the step (4) specifically includes the following steps:
(4.1) the grid central point of interpolation is considered as new sample point, generates a new Thiessen polygon;
(4.2) sample point in the Thiessen polygon adjacent with Thiessen polygon where interpolation grid is found out, it is described
Sample point is the adjoining sample point for participating in interpolation;
(4.3) weight is calculated according to the ratio of newly-generated Thiessen polygon and adjacent Thiessen polygon intersecting area.
Preferably, in the step (4.3) i-th of adjacent sample point weight w_i, specifically:
The weight w of i-th of adjacent sample point is calculated according to the following formulai:
Wherein, x is interpolation point, aiFor the area for participating in Thiessen polygon locating for the adjoining sample point of interpolation, a (x) is
The area of Thiessen polygon locating for interpolation point x, ai∩ a (x) is the area of the two intersection.
Preferably, adjoining sample point Tyson where interpolation grid of interpolation grid is more in the step (4.2)
Sample in the polygon of side shape all of its neighbor at.
Preferably, computation grid predicted value in the step (5), specifically:
Computation grid predicted value according to the following formula:
Wherein, f (x) is the field intensity prediction value of grid where interpolation point x, wiIt (x) is i-th (i=1 ... ..., n) a ginseng
Weight with the sample point of interpolation about interpolation point x, fiFor the field intensity value at sample point i.
This carries out the system that field strength distribution prediction is handled, main feature based on frequency spectrum drive test data using the above method
It is that the system includes:
Frequency spectrum data obtains module, for obtaining the frequency spectrum drive test data of estimation range;
Grid preprocessing module obtains module with the frequency spectrum data and is connected, for estimation range and Target area
The frequency spectrum drive test data in domain is pre-processed;
Field strength interpolation processing module is connected, for the grid to unknown field intensity value with the grid preprocessing module
Carry out interpolation processing;
Field intensity map generation module is connected with the field strength interpolation processing module, will for being divided according to color range
The field intensity prediction value of whole region is converted to field intensity map.
Preferably, the frequency spectrum data obtains the measurement that the frequency spectrum drive test data that module obtains includes each sampled data
The information such as time, longitude and latitude, frequency values, field intensity value.
Realize that field strength distribution predicts the method and corresponding system of processing based on frequency spectrum drive test data using of the invention,
The present invention can provide the field strength distribution situation of entire target area entirety, even if propagating in the very big different community of environmental difference
When, the error of prediction data also very little.The solution of the present invention does not survey the field intensity value in region according to the field intensity value prediction in actual measurement region,
To save cost of human and material resources.
Specific embodiment
It is further to carry out combined with specific embodiments below in order to more clearly describe technology contents of the invention
Description.
The method that field strength distribution prediction processing is realized based on frequency spectrum drive test data of the invention, including following steps:
(1) the frequency spectrum drive test data of estimation range is obtained;
(2) rasterizing pretreatment is carried out to estimation range and frequency spectrum drive test data;
(2.1) grid division is carried out to estimation range according to designated size;
(2.2) frequency spectrum drive test data is filtered into cleaning, removes incomplete, mistake and duplicate data;
(2.3) longitude and latitude is located at a plurality of frequency spectrum drive test data arithmetic mean within the scope of a grid is a data,
And the longitude and latitude of the central point of grid where being set as the longitude and latitude of the data, grid where the field intensity value of the data represents
Field intensity value;
(3) pretreated frequency spectrum data is considered as sample point creation Thiessen polygon;
(4) according to the weight of each adjacent sample point of Natural neighbors interpolation calculation, the grid of unknown field intensity value is carried out slotting
Value processing;
(4.1) the grid central point of interpolation is considered as new sample point, generates a new Thiessen polygon;
(4.2) sample point in the Thiessen polygon adjacent with Thiessen polygon where interpolation grid is found out, it is described
Sample point is the adjoining sample point for participating in interpolation;
(4.3) weight is calculated according to the ratio of newly-generated Thiessen polygon and adjacent Thiessen polygon intersecting area;
(5) the field intensity prediction value of the grid is calculated according to the weight of each adjacent sample point of interpolation grid and field intensity value;
(6) it is divided according to color range and the field intensity value of whole region is converted into field intensity map.
Preferably, frequency spectrum drive test data includes latitude and longitude information, frequency information and field intensity value in the step (1).
Preferably, pretreated frequency spectrum data is the frequency spectrum data after grid is average in the step (3), it is specific to wrap
Grid central point longitude and latitude, frequency information and corresponding field intensity value where including.
Preferably, in the step (4.3) i-th of adjacent sample point weight w_i, specifically:
The weight w of i-th of adjacent sample point is calculated according to the following formulai:
Wherein, x is interpolation point, aiFor the area for participating in Thiessen polygon locating for the adjoining sample point of interpolation, a (x) is
The area of Thiessen polygon locating for interpolation point x, ai∩ a (x) is the area of the two intersection.
Preferably, adjoining sample point Tyson where interpolation grid of interpolation grid is more in the step (4.2)
Sample in the polygon of side shape all of its neighbor at.
Preferably, computation grid predicted value in the step (5), specifically:
Computation grid predicted value according to the following formula:
Wherein, f (x) is the field intensity prediction value of grid where interpolation point x, wiIt (x) is i-th (i=1 ... ..., n) a ginseng
Weight with the sample point of interpolation about interpolation point x, fiFor the field intensity value at sample point i.
This of the invention carries out the system that field strength distribution prediction is handled based on frequency spectrum drive test data using the above method,
In include:
Frequency spectrum data obtains module, for obtaining the frequency spectrum drive test data of estimation range;
Grid preprocessing module obtains module with the frequency spectrum data and is connected, for estimation range and Target area
The frequency spectrum drive test data in domain is pre-processed;
Field strength interpolation processing module is connected, for the grid to unknown field intensity value with the grid preprocessing module
Carry out interpolation processing;
Field intensity map generation module is connected with the field strength interpolation processing module, will for being divided according to color range
The field intensity prediction value of whole region is converted to field intensity map.
As the preferred embodiment of the present invention, the frequency spectrum data obtains the frequency spectrum drive test data that module obtains and includes
The information such as time of measuring, longitude and latitude, frequency values, the field intensity value of each sampled data.
In a specific embodiment of the invention, the field strength distribution prediction based on frequency spectrum drive test data that the present invention provides a kind of
Method and device.Method includes:
Obtain the frequency spectrum drive test data of estimation range;Rasterizing pretreatment is carried out to estimation range and frequency spectrum drive test data;
Pretreated frequency spectrum data is considered as sample point creation Thiessen polygon;The grid central point of interpolation is considered as new sample
Point generates a new Thiessen polygon;It finds the adjoining sample point of interpolation grid and calculates the weight of each adjacent sample point;
The field intensity prediction value of the grid is calculated according to the weight of each adjacent sample point of interpolation grid and field intensity value;It is divided according to color range
The field intensity value of whole region is converted into field intensity map.The solution of the present invention predicts the area Wei Ce according to the field intensity value in actual measurement region
The field intensity value in domain, to save cost of human and material resources.
The present invention provides a kind of field strength distribution prediction technique based on frequency spectrum data, comprising the following steps:
(1) the frequency spectrum drive test data of estimation range is obtained;
(2) rasterizing pretreatment is carried out to estimation range and frequency spectrum drive test data;
(3) pretreated frequency spectrum data is considered as sample point creation Thiessen polygon;
(4) the grid central point of interpolation is considered as new sample point, generates a new Thiessen polygon;
(5) it finds the adjoining sample point of interpolation grid and calculates the weight of each adjacent sample point;
(6) the field intensity prediction value of the grid is calculated according to the weight of each adjacent sample point of interpolation grid and field intensity value;
(7) it is divided according to color range and the field intensity value of whole region is converted into field intensity map.
The frequency spectrum drive test data of estimation range at least contains following information in step (1): latitude and longitude information, frequency information and
Field intensity value.
Rasterizing pretreatment is carried out to estimation range in step (2) to specifically include according to designated size to estimation range progress
Grid division.
Rasterizing pretreatment is carried out to frequency spectrum drive test data in step (2) to specifically include to frequency spectrum drive test data into cleaning
Filter, removes incomplete, mistake and duplicate data;It further include a plurality of frequency spectrum being located at longitude and latitude within the scope of one grid
Drive test data arithmetic mean be a data, and by the longitude and latitude of the data be set as where grid central point longitude and latitude,
The field intensity value of grid where the field intensity value of the data represents.
Pretreated frequency spectrum data is considered as sample point creation Thiessen polygon in step (3), wherein pretreated
Frequency spectrum data is the frequency spectrum data after grid is average, grid central point longitude and latitude where specifically including, frequency information and corresponding
Field intensity value.
The adjoining sample point of interpolation grid is by the more of interpolation grid place Thiessen polygon all of its neighbor in step (5)
Sample in the shape of side at.
The weight of each adjacent sample point is calculated in step (5), specifically:
The weight Wi of i-th of adjacent sample point is calculated according to the following formula:
Wherein, x is interpolation point, aiFor the area for participating in Thiessen polygon locating for the adjoining sample point of interpolation, a (x) is
The area of Thiessen polygon locating for interpolation point x, ai∩ a (x) is the area of the two intersection.
The field strength for calculating the grid according to the weight of each adjacent sample point of interpolation grid and field intensity value in step (6) is pre-
Measured value, specifically:
Computation grid predicted value according to the following formula:
Wherein, f (x) is the field intensity prediction value of grid where interpolation point x, wi(x) (i=1 ..., n) a ginseng for i-th
Weight with the sample point of interpolation about interpolation point x, fiFor the field intensity value at sample point i.
Step (7) includes that the color that the grid field intensity value of whole region is corresponded to color range is divided according to color range, is exported whole
The field intensity map in a region.
In addition, the embodiment of the present invention also provides a kind of determining device for predicting field strength, comprising:
Frequency spectrum data obtains module, for obtaining the frequency spectrum drive test data of estimation range;Frequency spectrum drive test data includes each
The information such as time of measuring, longitude and latitude, frequency values, the field intensity value of sampled data;
Grid preprocessing module is pre-processed for the frequency spectrum drive test data to estimation range and estimation range;Including
Grid division is carried out according to specified granularity to estimation range;Road test data is filtered into cleaning, removes imperfect, mistake sum
Duplicate data;It is a data that longitude and latitude, which is located at a plurality of data arithmetic mean in grid, and by the data
The longitude and latitude of the central point of grid where longitude and latitude is set as;
Field strength interpolation processing module carries out interpolation processing for the grid to unknown field intensity value;
The field intensity prediction value of whole region is converted to field strength point for dividing according to color range by field intensity map generation module
Butut.
To keep the technical problem to be solved in the present invention, technical solution and advantage clearer, below in conjunction with attached drawing and implementation
Example, the present invention is described in detail.
A kind of field strength distribution prediction technique based on frequency spectrum drive test data is provided in specific embodiments of the present invention, comprising:
Step 1, the frequency spectrum drive test data of estimation range is obtained.When frequency spectrum drive test data includes the measurement of each sampled data
Between, the information such as longitude and latitude, frequency values, field intensity value;
Step 2, estimation range and frequency spectrum drive test data are pre-processed;
2.1, estimation range is subjected to grid division according to specified granularity;
2.2, frequency spectrum drive test data is filtered into cleaning, removes incomplete, mistake and duplicate data;
2.3, it is a data that longitude and latitude, which is located at a plurality of frequency spectrum drive test data arithmetic mean within the scope of a grid, and
The longitude and latitude of the central point of grid where the longitude and latitude of the data is set as;
Step 3, the data point after grid being managed is considered as sample point creation Thiessen polygon;
Step 4, according to the weight of each adjacent sample point of Natural neighbors interpolation calculation, the grid of unknown field intensity value is carried out
Interpolation processing:
4.1, the grid central point of interpolation is considered as new sample point, generates a new Thiessen polygon;
4.2, find the sample point in the Thiessen polygon adjacent with Thiessen polygon where interpolation grid, these samples
Point is the adjoining sample point for participating in interpolation;
4.3, weight is calculated according to the ratio of newly-generated Thiessen polygon and adjacent Thiessen polygon intersecting area;
Formula is as follows:
Wherein, f (x) is the interpolation result at interpolation point x, wi(x) (i=1 ..., n) a interpolation that participates in for i-th
Weight of the sample point about interpolation point x, fiFor the value at sample point i.
Weight is determined by following formula:
Wherein, aiFor the area for participating in Thiessen polygon locating for the sample point of interpolation, a (x) is Tyson locating for interpolation point x
The area of polygon, ai∩ a (x) is the area of the two intersection, as shown in Figure 3.
Step 5, the predicted value of the grid is calculated according to the weight of each adjacent sample point of interpolation grid and field intensity value;
Step 6, being divided according to color range by the prediction numerical value conversion in region is region field intensity map.
Realize that field strength distribution predicts the method and corresponding system of processing based on frequency spectrum drive test data using of the invention,
The present invention can provide the field strength distribution situation of entire target area entirety, even if propagating in the very big different community of environmental difference
When, the error of prediction data also very little.The solution of the present invention does not survey the field intensity value in region according to the field intensity value prediction in actual measurement region,
To save cost of human and material resources.
In this description, the present invention is described with reference to its specific embodiment.But it is clear that can still make
Various modifications and alterations are without departing from the spirit and scope of the invention.Therefore, the description and the appended drawings should be considered as illustrative
And not restrictive.