CN108197347B - method for detecting geographic space distribution uniformity of continuous irregular area sampling points - Google Patents
method for detecting geographic space distribution uniformity of continuous irregular area sampling points Download PDFInfo
- Publication number
- CN108197347B CN108197347B CN201711293026.XA CN201711293026A CN108197347B CN 108197347 B CN108197347 B CN 108197347B CN 201711293026 A CN201711293026 A CN 201711293026A CN 108197347 B CN108197347 B CN 108197347B
- Authority
- CN
- China
- Prior art keywords
- sampling
- area
- sampling points
- grid
- points
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Computational Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Mathematical Analysis (AREA)
- Remote Sensing (AREA)
- Data Mining & Analysis (AREA)
- Image Analysis (AREA)
- Complex Calculations (AREA)
Abstract
The invention provides a method for detecting the geographic spatial distribution uniformity of sampling points in a continuous irregular area, which comprises the steps of coding each sampling point in a sampling area to obtain a coding value of each sampling point; and calculating the geographic spatial distribution uniformity of all sampling points in the sampling area according to the coding values of all sampling points. The geospatial distribution uniformity detection method of the continuous irregular area sampling points, provided by the invention, comprises the steps of coding each sampling point, calculating the geospatial distribution uniformity of all the sampling points in the sampling area according to the coding values of all the sampling points, generating a uniform coding graph of all the sampling points, evaluating the uniformity degree of the continuous irregular area sampling points according to the uniformity degree, determining an abnormal unit grid through the uniform coding graph, ensuring the usability of the sampling point data, being beneficial to the redundancy removal refinement and mining analysis of the sampling point data in the next step, and ensuring the accuracy and reliability of the analysis and the practical application.
Description
Technical Field
The invention relates to the technical field of spatial mapping, in particular to a method for detecting geographic spatial distribution uniformity of continuous irregular area sampling points.
Background
the continuous irregular sampling area is generally an area which is continuous inside and irregular in boundary, such as a certain village and a certain town. Whether the distribution of the sampling points of the continuous irregular area is uniform is a key factor for later-stage sampling point data detection and evaluation and is also an important basis for judging whether the later-stage analysis result is accurate, so that the detection of the geographic space distribution uniformity of the sampling points of the continuous irregular area is very important.
in the prior art, the following method is generally adopted to detect the geospatial distribution uniformity of sampling points: assuming that N points are uniformly distributed in a certain sampling area, constructing a Thiessen polygon network of the points, and then calculating the change accumulation of the Thiessen polygon area to construct an index of the uniformity of the point distribution in a planewherein N is the number of sampling points, SiAnd (3) the area of the Thiessen polygon corresponding to the ith sampling point, and S is the nominal Thiessen polygon area, namely, the area of the whole sampling area is averagely divided into N points, and the area of the Thiessen polygon corresponding to each point is obtained. The index reflects the uniformity of the distribution of the points in a plane space, N sampling points are distributed in the same area, different distributions may bring different uniformities, the uniformity of the sampling points can be measured by the dimensionless E, and the higher the E value is, the more uniform the corresponding distribution is.
The method in the prior art does not distinguish the specific types of fragmentation, continuous irregularity, continuous regulation and the like of a sampling area, and does not aim at the specific sampling area of the continuous irregular area, and an index system of the method is single and lacks systematicness, so that the detection result of the geographic space distribution uniformity of sampling points in the continuous irregular area is easy to deviate, local uniform details of the sampling point data cannot be reflected, and the redundancy removal refinement and mining analysis of the sampling point data are not facilitated in the next step.
Disclosure of Invention
technical problem to be solved
the invention aims to provide a method for detecting the geographic spatial distribution uniformity of continuous irregular area sampling points, which solves the technical problems that the geographic spatial distribution uniformity of the continuous irregular area sampling points is low in detection accuracy and reliability, and sampling points with non-uniform geographic spatial distribution of all the sampling points in the continuous irregular area cannot be positioned in the prior art.
(II) technical scheme
in order to solve the above technical problem, in one aspect, the present invention provides a method for detecting a geospatial distribution uniformity of continuous irregular area sampling points, including:
generating a first-level grid, and covering a sampling region with the first-level grid, wherein the first-level grid comprises a plurality of first-level unit grids, the first-level unit grids are square grids with the area equal to the average sampling area, the sampling region is a continuous irregular region, and the average sampling area is the ratio of the total area of the sampling region to the number of sampling points;
If the number of the sampling points in the first-level unit grid is judged to be larger than one, dividing the first-level unit grid into four second-level unit grids with equal areas, wherein the second-level unit grids are square grids;
If the number of the sampling points in the second-level unit grid is judged to be less than or equal to one, coding each sampling point in the sampling area to obtain a coding value of each sampling point;
And calculating the geographic spatial distribution uniformity of all sampling points in the sampling region according to the code values of all sampling points in the sampling region.
further, after the calculating the geospatial distribution uniformity of all sampling points in the sampling region according to the coded values of all sampling points in the sampling region, the method further includes:
generating a uniform coding graph of all sampling points in the sampling area according to the coding values of all sampling points in the sampling area;
And determining an abnormal unit grid in the sampling area according to the uniform coding pattern, wherein the abnormal unit grid is a unit grid with a coding value not equal to 1.
Further, the encoding each sampling point in the sampling region specifically includes:
wherein, Ciis the coded value of the ith sampling point, and t is the level of the unit grid where the ith sampling point is located.
Further, the calculating the geospatial distribution uniformity of all sampling points in the sampling region according to the code values of all sampling points in the sampling region specifically includes:
Wherein E is the geospatial distribution uniformity of all sampling points, CiIs the code value of the ith sampling point, and N is the number of sampling points in the sampling area.
Further, the generating a uniform code pattern of all sampling points in the sampling region according to the code values of all sampling points in the sampling region specifically includes:
Acquiring the coding value of a unit grid corresponding to each sampling point according to the coding value of the sampling point, wherein the coding value of each unit grid is equal to the coding value of the sampling point in the unit grid, and the coding value of the unit grid without the sampling point in the unit grid is 0;
And marking the code value of each unit grid into the corresponding unit grid, and generating a uniform code pattern of all sampling points in the sampling area.
In another aspect, the present invention provides a device for detecting geospatial distribution uniformity of continuous irregular area sampling points, comprising:
The grid generating module is used for generating a first-level grid and covering a sampling region by the first-level grid, wherein the first-level grid comprises a plurality of first-level unit grids, the first-level unit grids are square grids with the area equal to the average sampling area, the sampling region is a continuous irregular region, and the average sampling area is the ratio of the total area of the sampling region to the number of sampling points;
the dividing module is used for dividing the first-stage unit grid into four second-stage unit grids with equal areas if the number of the sampling points in the first-stage unit grid is judged to be larger than one, wherein the second-stage unit grids are square grids;
the coding module is used for coding each sampling point in the sampling area to obtain a coding value of each sampling point if the number of the sampling points in the second-level unit grid is judged to be less than or equal to one;
And the calculation module is used for calculating the geographic space distribution uniformity of all sampling points in the sampling area according to the code values of all sampling points in the sampling area.
further, the method also comprises a code pattern generating module and a determining module, wherein,
The code pattern generating module is used for generating a uniform code pattern of all sampling points in the sampling area according to the code values of all sampling points in the sampling area;
And the determining module is used for determining an abnormal unit grid in the sampling area according to the uniform code pattern, wherein the abnormal unit grid is a unit grid with a code value not equal to 1.
In yet another aspect, the present invention provides an electronic device for geospatial distribution uniformity detection of consecutive irregular area sampling points, comprising:
The processor and the memory are communicated with each other through a bus; the memory stores program instructions executable by the processor, which when called by the processor are capable of performing the methods described above.
in yet another aspect, the invention provides a computer program product comprising a computer program stored on a non-transitory computer readable storage medium, the computer program comprising program instructions which, when executed by a computer, cause the computer to perform the method described above.
in a further aspect, the invention provides a computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements the method described above.
(III) advantageous effects
the method for detecting the geographic spatial distribution uniformity of the continuous irregular area sampling points, provided by the invention, has the advantages that each sampling point is coded, the geographic spatial distribution uniformity of all the sampling points in the sampling area is calculated according to the coding values of all the sampling points, the uniformity degree of the sampling points is evaluated according to the uniformity degree, the evaluation accuracy and reliability are high, and the availability of sampling data is evaluated according to the geographic spatial distribution uniformity of the sampling points, so that the accuracy and reliability of the specific analysis and the practical application of the sampling point data are ensured. And the geospatial distribution uniformity of all sampling points in the sampling region can be determined by combining the uniform coding pattern, and the abnormal unit grid in the sampling region is determined according to the uniform coding pattern, so that the abnormal unit grid can be further processed conveniently, and redundancy removal, refinement, mining and analysis of the sampling point data in the next step are facilitated.
drawings
FIG. 1 is a schematic diagram illustrating a method for detecting a geospatial distribution uniformity of continuous irregular area sampling points according to an embodiment of the present invention;
FIG. 2 is a logic flow diagram of a method for detecting geospatial distribution uniformity of consecutive irregular area sampling points in accordance with an embodiment of the present invention;
FIG. 3 is a schematic diagram of a uniform code pattern for all sampling points in a sampling region according to an embodiment of the present invention;
FIG. 4 is a schematic diagram of an apparatus for detecting the geospatial distribution uniformity of consecutive irregular area sampling points according to an embodiment of the present invention;
fig. 5 is a schematic structural diagram of an electronic device for detecting the geospatial distribution uniformity of consecutive irregular area sampling points according to an embodiment of the present invention.
Detailed Description
in order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments of the present invention without making any creative effort, shall fall within the protection scope of the present invention.
Example 1:
Fig. 1 is a schematic diagram illustrating a method for detecting a geospatial distribution uniformity of consecutive irregular area sampling points according to an embodiment of the present invention, and as shown in fig. 1, an embodiment of the present invention provides a method for detecting a geospatial distribution uniformity of consecutive irregular area sampling points, including:
step S10, generating a first-level grid, and covering a sampling region with the first-level grid, wherein the first-level grid comprises a plurality of first-level unit grids, the first-level unit grids are square grids with the area equal to the average sampling area, the sampling region is a continuous irregular region, and the average sampling area is the ratio of the total area of the sampling region to the number of sampling points;
Step S20, if the number of sampling points in the first-level unit grid is judged and known to be more than one, dividing the first-level unit grid into four second-level unit grids with equal areas, wherein the second-level unit grids are square grids;
step S30, if the number of the sampling points in the second-level unit grid is judged and obtained to be less than or equal to one, coding each sampling point in the sampling area to obtain a coding value of each sampling point;
And step S40, calculating the geospatial distribution uniformity of all sampling points in the sampling area according to the code values of all sampling points in the sampling area.
further, after the calculating the geospatial distribution uniformity of all sampling points in the sampling region according to the coded values of all sampling points in the sampling region, the method further includes:
generating a uniform coding graph of all sampling points in the sampling area according to the coding values of all sampling points in the sampling area;
And determining an abnormal unit grid in the sampling area according to the uniform coding pattern, wherein the abnormal unit grid is a unit grid with a coding value not equal to 1.
further, the encoding each sampling point in the sampling region specifically includes:
wherein, CiIs the coded value of the ith sampling point, and t is the level of the unit grid where the ith sampling point is located.
Further, the calculating the geospatial distribution uniformity of all sampling points in the sampling region according to the code values of all sampling points in the sampling region specifically includes:
wherein E is the geospatial distribution uniformity of all sampling points, Ciis the code value of the ith sampling point, and N is the number of sampling points in the sampling area.
Further, the generating a uniform code pattern of all sampling points in the sampling region according to the code values of all sampling points in the sampling region specifically includes:
Acquiring the coding value of a unit grid corresponding to each sampling point according to the coding value of the sampling point, wherein the coding value of each unit grid is equal to the coding value of the sampling point in the unit grid, and the coding value of the unit grid without the sampling point in the unit grid is 0;
And marking the code value of each unit grid into the corresponding unit grid, and generating a uniform code pattern of all sampling points in the sampling area.
Specifically, fig. 2 is a logic flow diagram of a geospatial distribution uniformity detection method of continuous irregular area sampling points according to an embodiment of the present invention, and as shown in fig. 2, the method of the embodiment of the present invention is further described below by using the specific logic flow diagram.
firstly, acquiring distribution data of a sampling region and sampling points to be evaluated, counting the total area S of the sampling region and the total number N of the sampling points, and calculating the average sampling area S of the sampling points with uniformly distributed sampling points under an ideal condition0. Wherein the average sampling area S0the calculation formula of (a) is as follows:
wherein S is0is the average sampling area, S is the total area of the sampling region, and N is the number of sampling points.
In this embodiment, taking the distribution data of sampling points for pollution monitoring in large villages and towns in great prosperous area in beijing as an example, the geospatial uniformity of the pollution monitoring sampling points is detected. The sampling area is a continuous irregular area, and the total area S of the sampling area is calculated to be 109.384km2the total number of sampling points N is 27 according to the average sampling area S0calculating the average sampling area S0is 4.051km2；
then, a unit grid area S is generated0The square grid is used as a first-level grid, the first-level grid is used for covering the whole sampling region, and M first-level grids containing the sampling region are extracted and obtained based on space superposition analysis. Wherein the spatial overlay analysis can be implemented based on ArcGIS software; the total number of the extracted first-stage grids is M, wherein M is larger than or equal to N.
In this example, the unit grid area is 4.051km2the square grids cover numerous village and town sampling areas, and 39 square grids containing the sampling areas are extracted based on the spatial superposition analysis of ArcGIS software. Generating a square grid number M which is 39 and is larger than the total number of the sampling points which is 27;
and then, judging whether the number of the sampling points in the first-level unit grid is more than one.
and if the number of the sampling points in the first-level unit grid is judged to be larger than one, the first-level unit grid is divided into four second-level unit grids with equal areas, the second-level unit grids are square grids, the second-level grids divided by the first-level unit grid are sub-grids of the first-level unit grid, and the divided second-level grids comprise a plurality of second-level unit grids.
and if the maximum number of the sampling points in the second-level unit grid is judged and obtained, the grid division is finished.
And if the number of the sampling points in the second-level unit grid is judged and obtained to be still larger than one, performing secondary division by taking the second-level unit grid as an object.
recursion operation is carried out in the mode until the maximum number of sampling points in each unit grid is met, the grid is divided, the unit grids comprise first-level unit grids, second-level unit grids, third-level unit grids and the like, and the unit grids are generated after division.
in this embodiment, a grid division termination condition that the number of sampling points in each unit grid is at most one is satisfied, and a grid division scheme of this time is counted as follows: the first grid division is 39, namely the first level unit grids are 39; 4 secondary dividing grids, namely 4 first-level unit grids are subjected to secondary division; the three-time division grids are 1, namely 1 second-level unit grid is subjected to three-time division.
And then, carrying out coding operation on the divided unit grids, and coding each sampling point to obtain a coding value of each sampling point, and finally generating a uniform coding graph of all sampling points in the sampling area.
Fig. 3 is a schematic diagram of a uniform code pattern of all sampling points in a sampling region according to an embodiment of the present invention, as shown in fig. 3, wherein grids with decreasing areas are a first-level unit grid, a second-level unit grid, a third-level unit grid, …, and a t-th-level unit grid in sequence. The second-level mesh, the third-level mesh, … and the t-level mesh are collectively referred to as sub-meshes, and each first-level unit mesh, each second-level unit mesh, each third-level unit mesh, … and each t-level unit mesh are collectively referred to as unit meshes.
uniform coding gauge for sampling pointsthen it is: for a unit grid with sampling points, firstly dividing the codes into 1, namely, the first-level unit grid codes are 1, and the coding values of the sampling points in the first-level unit grid are 1; the secondary partition code is 1/4, that is, the second-level unit grid code is 1/4, and the code values of the sampling points in the second-level unit grid are 1/4; the third-level unit grid code is 1/16, namely 1/16, and the code value of the sampling point in the third-level unit grid is 1/16; by analogy, the t-th partition is coded as 1/4t-1I.e. the t-th level unit grid is coded as 1/4t-1the coding value of the sampling point in the t-th level unit grid is 1/4t-1(ii) a Is uniformly coded to 0 for a unit grid without sampling points. And marking the coding value of each sampling point in the unit grid corresponding to the sampling point, and marking the unit grid without the sampling points by using the coding value 0 to generate a uniform coding graph of all the sampling points in the sampling area.
In this embodiment, the grid division is divided into three levels. The first level mesh encodes 35, which contains 18 "1" s and 17 "0" s; the second level mesh encodes 15, of which 7 "1/4" and 8 "0" are included; the third level mesh encodes 4 cells, which contain 2 "1/16" and 2 "0" s. The unit grids whose code values are not equal to 1 are abnormal unit grids, that is, the unit grids whose code values are equal to 0, 1/4, and 1/16 are abnormal unit grids. Sampling points in the abnormal unit grid are abnormal sampling points, and the abnormal sampling points can cause that the geographic space distribution of all the sampling points in the sampling area is not uniform. The abnormal unit grids and the abnormal sampling points can be rapidly and visually determined from the uniform coding graph, and the abnormal unit grids and the abnormal sampling points are processed according to actual needs, so that redundancy removal, refinement, mining and analysis of sampling point data in the next step are facilitated.
Then, according to the code values of all sampling points in the sampling region, calculating the geographic spatial distribution uniformity of all sampling points in the sampling region, wherein the specific calculation formula is as follows:
Wherein E is the geospatial distribution uniformity of all sampling points, Ciis the code value of the ith sampling point, and N is the number of sampling points in the sampling area. The larger the uniformity value, the more uniform the geospatial distribution of the sampling point data, and on the contrary, the more non-uniform the geospatial distribution of the sampling point data.
In this embodiment, the geospatial distribution uniformity of the sampling points in each village and town sampling area of the great-interest area of beijing city is calculated to be 0.736 based on the specific calculation formula of the geospatial distribution uniformity E of all the sampling points.
Uniformity E can be used to characterize the uniformity of the distribution of the detected sampling points in geographic space. The uniformity of the distribution of the sampling points in the geographic space can also be characterized by the uniformity code pattern and the uniformity E of the sampling points. Wherein, the local uniform details of sampling points of each village and town of the great happy area of Beijing City are represented by the uniform coding chart in FIG. 3; the overall uniformity of sampling points of each village and town of the great prosperity area of Beijing is represented by uniformity E, and the uniformity of the sampling points in the sampling area is 0.736. The detection results of the uniform code pattern and the uniformity of the sampling points jointly represent that the detection method for the sampling points in the continuous irregular area, which is provided by the invention, is very accurate, reliable and effective.
the method for detecting the geographic spatial distribution uniformity of the continuous irregular area sampling points, provided by the invention, has the advantages that each sampling point is coded, the geographic spatial distribution uniformity of all the sampling points in the sampling area is calculated according to the coding values of all the sampling points, the uniformity degree of the sampling points is evaluated according to the uniformity degree, the evaluation accuracy and reliability are high, and the availability of sampling data is evaluated according to the geographic spatial distribution uniformity of the sampling points, so that the accuracy and reliability of the specific analysis and the practical application of the sampling point data are ensured. And the geospatial distribution uniformity of all sampling points in the sampling region can be determined by combining the uniform coding pattern, and the abnormal unit grid in the sampling region is determined according to the uniform coding pattern, so that the abnormal unit grid can be further processed conveniently, and redundancy removal, refinement, mining and analysis of the sampling point data in the next step are facilitated.
example 2:
fig. 4 is a schematic diagram of an apparatus for detecting geospatial distribution uniformity of consecutive irregular area sampling points according to an embodiment of the present invention, and as shown in fig. 4, an embodiment of the present invention provides an apparatus for detecting geospatial distribution uniformity of consecutive irregular area sampling points, which is used for performing the method described in the above embodiment, and specifically includes a grid generating module 10, a dividing module 20, an encoding module 30 and a calculating module 40, wherein,
The grid generating module 10 is configured to generate a first-level grid, and cover a sampling region with the first-level grid, where the first-level grid includes a plurality of first-level unit grids, the first-level unit grid is a square grid having an area equal to an average sampling area, the sampling region is a continuous irregular region, and the average sampling area is a ratio between a total area of the sampling region and the number of sampling points;
The dividing module 20 is configured to, if it is determined that the number of sampling points in the first-level unit grid is greater than one, divide the first-level unit grid into four second-level unit grids with equal areas, where each second-level unit grid is a square grid;
The encoding module 30 is configured to encode each sampling point in the sampling region to obtain an encoded value of each sampling point if it is determined that the number of sampling points in the second-level unit grid is less than or equal to one;
The calculating module 40 is configured to calculate a geospatial distribution uniformity of all sampling points in the sampling region according to the code values of all sampling points in the sampling region.
further, the method also comprises a code pattern generating module and a determining module, wherein,
The code pattern generating module is used for generating a uniform code pattern of all sampling points in the sampling area according to the code values of all sampling points in the sampling area;
And the determining module is used for determining an abnormal unit grid in the sampling area according to the uniform code pattern, wherein the abnormal unit grid is a unit grid with a code value not equal to 1.
the geospatial distribution uniformity detection device for the continuous irregular area sampling points, provided by the invention, has the advantages that each sampling point is coded, the geospatial distribution uniformity of all the sampling points in the sampling area is calculated according to the coding values of all the sampling points, the uniformity degree of the sampling points is evaluated according to the uniformity degree, the evaluation accuracy and reliability are high, and the availability of sampling data is evaluated according to the geospatial distribution uniformity of the sampling points, so that the accuracy and reliability of the specific analysis and the practical application of the sampling point data are ensured. And the geospatial distribution uniformity of all sampling points in the sampling region can be determined by combining the uniform coding pattern, and the abnormal unit grid in the sampling region is determined according to the uniform coding pattern, so that the abnormal unit grid can be further processed conveniently, and redundancy removal, refinement, mining and analysis of the sampling point data in the next step are facilitated.
Example 3:
Fig. 5 is a schematic structural diagram of an electronic device for detecting geospatial distribution uniformity of consecutive irregular area sampling points according to an embodiment of the present invention, and as shown in fig. 5, the device includes: a processor (processor)801, a memory (memory)802, and a bus (bus) 803;
wherein, the processor 801 and the memory 802 complete the communication with each other through the bus 803;
The processor 801 is configured to call program instructions in the memory 802 to perform the methods provided by the above-described method embodiments, including, for example:
generating a first-level grid, and covering a sampling region with the first-level grid, wherein the first-level grid comprises a plurality of first-level unit grids, the first-level unit grids are square grids with the area equal to the average sampling area, the sampling region is a continuous irregular region, and the average sampling area is the ratio of the total area of the sampling region to the number of sampling points;
if the number of the sampling points in the first-level unit grid is judged to be larger than one, dividing the first-level unit grid into four second-level unit grids with equal areas, wherein the second-level unit grids are square grids;
if the number of the sampling points in the second-level unit grid is judged to be less than or equal to one, coding each sampling point in the sampling area to obtain a coding value of each sampling point;
And calculating the geographic spatial distribution uniformity of all sampling points in the sampling region according to the code values of all sampling points in the sampling region.
Example 4:
An embodiment of the present invention discloses a computer program product, which includes a computer program stored on a non-transitory computer readable storage medium, the computer program including program instructions, when the program instructions are executed by a computer, the computer can execute the methods provided by the above method embodiments, for example, the method includes:
generating a first-level grid, and covering a sampling region with the first-level grid, wherein the first-level grid comprises a plurality of first-level unit grids, the first-level unit grids are square grids with the area equal to the average sampling area, the sampling region is a continuous irregular region, and the average sampling area is the ratio of the total area of the sampling region to the number of sampling points;
If the number of the sampling points in the first-level unit grid is judged to be larger than one, dividing the first-level unit grid into four second-level unit grids with equal areas, wherein the second-level unit grids are square grids;
if the number of the sampling points in the second-level unit grid is judged to be less than or equal to one, coding each sampling point in the sampling area to obtain a coding value of each sampling point;
And calculating the geographic spatial distribution uniformity of all sampling points in the sampling region according to the code values of all sampling points in the sampling region.
example 5:
embodiments of the present invention provide a non-transitory computer-readable storage medium, which stores computer instructions, where the computer instructions cause the computer to perform the methods provided by the above method embodiments, for example, the methods include:
generating a first-level grid, and covering a sampling region with the first-level grid, wherein the first-level grid comprises a plurality of first-level unit grids, the first-level unit grids are square grids with the area equal to the average sampling area, the sampling region is a continuous irregular region, and the average sampling area is the ratio of the total area of the sampling region to the number of sampling points;
If the number of the sampling points in the first-level unit grid is judged to be larger than one, dividing the first-level unit grid into four second-level unit grids with equal areas, wherein the second-level unit grids are square grids;
If the number of the sampling points in the second-level unit grid is judged to be less than or equal to one, coding each sampling point in the sampling area to obtain a coding value of each sampling point;
and calculating the geographic spatial distribution uniformity of all sampling points in the sampling region according to the code values of all sampling points in the sampling region.
Those of ordinary skill in the art will understand that: all or part of the steps for implementing the method embodiments may be implemented by hardware related to program instructions, and the program may be stored in a computer readable storage medium, and when executed, the program performs the steps including the method embodiments; and the aforementioned storage medium includes: various media that can store program codes, such as ROM, RAM, magnetic or optical disks.
The above-described embodiments of the apparatuses and devices are merely illustrative, where the units described as separate parts may or may not be physically separate, and the parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods described in the embodiments or some parts of the embodiments.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (7)
1. A geospatial distribution uniformity detection method for continuous irregular area sampling points is characterized by comprising the following steps:
Generating a first-level grid, and covering a sampling region with the first-level grid, wherein the first-level grid comprises a plurality of first-level unit grids, the first-level unit grids are square grids with the area equal to the average sampling area, the sampling region is a continuous irregular region, and the average sampling area is the ratio of the total area of the sampling region to the number of sampling points;
if the number of the sampling points in the first-level unit grid is judged to be larger than one, dividing the first-level unit grid into four second-level unit grids with equal areas, wherein the second-level unit grids are square grids;
If the number of the sampling points in the second-level unit grid is judged to be less than or equal to one, coding each sampling point in the sampling area to obtain a coding value of each sampling point;
calculating the geographic spatial distribution uniformity of all sampling points in the sampling area according to the coding values of all sampling points in the sampling area;
the encoding of each sampling point in the sampling region specifically includes:
wherein, Cithe coding value of the ith sampling point is, and t is the level of the unit grid where the ith sampling point is located;
The calculating the geospatial distribution uniformity of all sampling points in the sampling region according to the code values of all sampling points in the sampling region specifically comprises:
Wherein E is the geospatial distribution uniformity of all sampling points, CiIs the code value of the ith sampling point, and N is the number of sampling points in the sampling area.
2. the method according to claim 1, wherein the calculating the geospatial distribution uniformity of all the sampling points in the sampling region according to the coded values of all the sampling points in the sampling region further comprises:
Generating a uniform coding graph of all sampling points in the sampling area according to the coding values of all sampling points in the sampling area;
And determining an abnormal unit grid in the sampling area according to the uniform coding pattern, wherein the abnormal unit grid is a unit grid with a coding value not equal to 1.
3. the method according to claim 2, wherein the generating a uniform code pattern of all sampling points in the sampling region according to the code values of all sampling points in the sampling region specifically comprises:
Acquiring the coding value of a unit grid corresponding to each sampling point according to the coding value of the sampling point, wherein the coding value of each unit grid is equal to the coding value of the sampling point in the unit grid, and the coding value of the unit grid without the sampling point in the unit grid is 0;
and marking the code value of each unit grid into the corresponding unit grid, and generating a uniform code pattern of all sampling points in the sampling area.
4. a geospatial distribution uniformity detection device for continuous irregular area sampling points, comprising:
The grid generating module is used for generating a first-level grid and covering a sampling region by the first-level grid, wherein the first-level grid comprises a plurality of first-level unit grids, the first-level unit grids are square grids with the area equal to the average sampling area, the sampling region is a continuous irregular region, and the average sampling area is the ratio of the total area of the sampling region to the number of sampling points;
The dividing module is used for dividing the first-stage unit grid into four second-stage unit grids with equal areas if the number of the sampling points in the first-stage unit grid is judged to be larger than one, wherein the second-stage unit grids are square grids;
The coding module is used for coding each sampling point in the sampling area to obtain a coding value of each sampling point if the number of the sampling points in the second-level unit grid is judged to be less than or equal to one;
The calculation module is used for calculating the geographic space distribution uniformity of all sampling points in the sampling area according to the code values of all the sampling points in the sampling area;
the encoding of each sampling point in the sampling region specifically includes:
Wherein, CiThe coding value of the ith sampling point is, and t is the level of the unit grid where the ith sampling point is located;
The calculating the geospatial distribution uniformity of all sampling points in the sampling region according to the code values of all sampling points in the sampling region specifically comprises:
wherein E is the geospatial distribution uniformity of all sampling points, CiIs the code value of the ith sampling point, and N is the number of sampling points in the sampling area.
5. The apparatus of claim 4, further comprising a code pattern generation module and a determination module, wherein,
the code pattern generating module is used for generating a uniform code pattern of all sampling points in the sampling area according to the code values of all sampling points in the sampling area;
and the determining module is used for determining an abnormal unit grid in the sampling area according to the uniform code pattern, wherein the abnormal unit grid is a unit grid with a code value not equal to 1.
6. an electronic device for geospatial distribution uniformity detection of continuous irregular area sampling points, comprising:
the processor and the memory are communicated with each other through a bus; the memory stores program instructions executable by the processor, the processor invoking the program instructions to perform the method of any of claims 1 to 3.
7. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the method according to any one of claims 1 to 3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711293026.XA CN108197347B (en) | 2017-12-08 | 2017-12-08 | method for detecting geographic space distribution uniformity of continuous irregular area sampling points |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711293026.XA CN108197347B (en) | 2017-12-08 | 2017-12-08 | method for detecting geographic space distribution uniformity of continuous irregular area sampling points |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108197347A CN108197347A (en) | 2018-06-22 |
| CN108197347B true CN108197347B (en) | 2019-12-17 |
Family
ID=62573712
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711293026.XA Active CN108197347B (en) | 2017-12-08 | 2017-12-08 | method for detecting geographic space distribution uniformity of continuous irregular area sampling points |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108197347B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110426230A (en) * | 2019-08-08 | 2019-11-08 | 中山市疾病预防控制中心 | A kind of appraisal procedure of Food Monitoring sampled point spatial distribution |
| CN112785910B (en) * | 2019-11-07 | 2023-07-25 | 中国石油天然气集团有限公司 | Large dynamic range nonlinear geophysical contour drawing method and device |
| CN113516731B (en) * | 2021-05-14 | 2022-03-01 | 北京市农林科学院信息技术研究中心 | Sample point data refinement method and system |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106021480A (en) * | 2016-05-19 | 2016-10-12 | 福建农林大学 | Parallel space declustering method and system based on mesh declustering |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102750361B (en) * | 2012-06-13 | 2015-02-18 | 中国矿业大学 | Discovery method of track data hot spot based on local multilayer grids |
| GB2532948B (en) * | 2014-12-02 | 2021-04-14 | Vivo Mobile Communication Co Ltd | Object Recognition in a 3D scene |
| CN105701250A (en) * | 2016-03-04 | 2016-06-22 | 农业部环境保护科研监测所 | Agricultural producing area soil monitoring point arrangement method based on step-by-step grid refinement |
| CN107092680B (en) * | 2017-04-21 | 2019-12-10 | 中国测绘科学研究院 | Government affair information resource integration method based on geographic grids |
-
2017
- 2017-12-08 CN CN201711293026.XA patent/CN108197347B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106021480A (en) * | 2016-05-19 | 2016-10-12 | 福建农林大学 | Parallel space declustering method and system based on mesh declustering |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108197347A (en) | 2018-06-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108197347B (en) | method for detecting geographic space distribution uniformity of continuous irregular area sampling points | |
| CN110689084B (en) | Abnormal user identification method and device | |
| CN111314173A (en) | Monitoring information abnormity positioning method and device, computer equipment and storage medium | |
| CN108805174A (en) | clustering method and device | |
| CN113992340B (en) | User abnormal behavior identification method, device, equipment and storage medium | |
| CN115329904B (en) | Outlier detection method, system and storage medium for discriminating population data | |
| CN112381272A (en) | Power grid load prediction method, system, terminal and storage medium | |
| CN114722091A (en) | Data processing method, data processing device, storage medium and processor | |
| CN110210774B (en) | Landslide risk evaluation method and system | |
| CN113657644B (en) | False data detection method, device, equipment and storage medium in smart power grid | |
| CN114446019A (en) | Alarm information processing method, device, equipment, storage medium and product | |
| CN116610821B (en) | Knowledge graph-based enterprise risk analysis method, system and storage medium | |
| CN112398226A (en) | Power supply system electricity stealing prevention method, system, terminal and storage medium | |
| CN106528048B (en) | Method and apparatus for evaluating quality of random number generator | |
| CN110390160B (en) | Method and device for detecting period of time sequence signal and related equipment | |
| CN108133234B (en) | Sparse subset selection algorithm-based community detection method, device and equipment | |
| CN111783883A (en) | Abnormal data detection method and device | |
| CN115618714A (en) | Solar radiation probability forecasting method and related equipment | |
| CN109059989B (en) | Method, system and equipment for evaluating remaining life of instrument | |
| CN114679335A (en) | Network security risk assessment training and assessment method and equipment for power monitoring system | |
| CN114124456A (en) | Safety detection method and device for comprehensive energy system, electronic equipment and storage medium | |
| CN111125685A (en) | Method and device for predicting network security situation | |
| CN114219680A (en) | Method, device, equipment and storage medium for checking topology of transformer area | |
| CN106295026B (en) | Traffic similarity analysis method and device | |
| CN116932255A (en) | Scheduling method and device in cloud environment, electronic equipment and storage medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |