CN104199937A - Multi-website POI position mapping method and device - Google Patents
Multi-website POI position mapping method and device Download PDFInfo
- Publication number
- CN104199937A CN104199937A CN201410455721.1A CN201410455721A CN104199937A CN 104199937 A CN104199937 A CN 104199937A CN 201410455721 A CN201410455721 A CN 201410455721A CN 104199937 A CN104199937 A CN 104199937A
- Authority
- CN
- China
- Prior art keywords
- mtd
- mtr
- msub
- msup
- poi
- 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.)
- Pending
Links
- 238000013507 mapping Methods 0.000 title claims abstract description 151
- 238000000034 method Methods 0.000 title claims abstract description 68
- 238000012937 correction Methods 0.000 claims abstract description 131
- 230000009466 transformation Effects 0.000 claims description 24
- 230000011218 segmentation Effects 0.000 claims description 19
- 239000011159 matrix material Substances 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 15
- 238000009827 uniform distribution Methods 0.000 claims description 12
- 238000011156 evaluation Methods 0.000 claims description 10
- 108700041286 delta Proteins 0.000 claims description 4
- 230000000694 effects Effects 0.000 description 13
- 238000005516 engineering process Methods 0.000 description 6
- 238000011426 transformation method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000006073 displacement reaction Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 230000003993 interaction Effects 0.000 description 4
- 238000007796 conventional method Methods 0.000 description 3
- 238000013506 data mapping Methods 0.000 description 3
- 238000011549 displacement method Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000002452 interceptive effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/95—Retrieval from the web
- G06F16/953—Querying, e.g. by the use of web search engines
- G06F16/9537—Spatial or temporal dependent retrieval, e.g. spatiotemporal queries
-
- 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)
- Databases & Information Systems (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
Abstract
The invention relates to a multi-website POI position mapping method and device. The multi-website POI position mapping method includes the following specific steps: (S110) building a correction unit set; (S120) conducting area splitting in each correction unit of the correction unit set; (S130) conducting control point matching based on semantics; (S140) automatically and evenly obtaining a multi-website POI same-name control point set and an examination and checking point set; (S150) resolving coefficients of a position mapping model; (S160) evaluating the position mapping accuracy. By means of the multi-website POI position mapping method and device, POI positions are evenly divided into a plurality of units, control points are evenly selected, the control points are completely and automatically selected by comprehensively using the semantic matching and the least-square method, and the optimal position mapping model is built according to the dividing-and-conquering strategy.
Description
Technical Field
The invention relates to the field of spatial information, in particular to a spatial position mapping method and a spatial position mapping device, which are used for establishing a position mapping relation among a plurality of website POIs.
Background
POI (Point Of Interest) is important content Of geographic information Service, and plays an important role in the fields Of emergency command, intelligent transportation, public safety, logistics management, electronic commerce, various other LBS (Location Based Service) services and the like. With the development of internet technology, internet POI information resources are increasingly abundant, and rich PO.I data contained in the internet is used as an important data source for expanding and updating a POI database. However, the multi-website POI data source has a problem of inconsistent position information, and the establishment of the position mapping relationship in the conventional manual/human-computer interaction manner is time-consuming, labor-consuming and low in precision, so that a method for automatically establishing the multi-website POI position mapping relationship is urgently needed.
1. Control point selecting method
The method comprises the steps of establishing a position mapping relation of the POI of the multi-website, and reasonably selecting a control point. The number of control points and the point location distribution will seriously affect the accuracy of the position mapping. An improper point location distribution directly leads to local deformation. In order to reduce adverse effects caused by inaccurate control point selection and uneven distribution, a large number of redundant control points are needed, and simultaneously, the selected POI control points should be uniformly distributed in the geographical range where the POI point set is located. The number, distribution and accuracy of the control points will directly affect the effect of the position mapping.
In the prior art, there are three main methods for selecting control points: the method comprises the steps of purely manually selecting control points, selecting control points in a man-machine interactive mode and matching and automatically generating control points.
(1) Manually selecting control points: the selection mode of the control points of the POI of different websites can be used for manually interpreting the attribute information of the POI according to the experience of an operator. And (3) the POI coordinate information to be mapped is (r, c), the POI coordinate information of the reference is (x, y), and a corresponding coordinate point pair is established, so that the control point is generated. The method has larger dependence on the workload and experience of an operator, time cost is large, and the selection precision of the control point is difficult to ensure.
(2) Selecting a ground control point in a man-machine interactive mode: the method is characterized in that an operator selects a plurality of control points through a man-machine interaction mode according to experience, rough conversion parameters are determined by utilizing the control points, and then sufficient control points are obtained in the geographic range of the POI through an identification technology, so that the conversion parameters are determined more accurately. The method realizes man-machine combination, reduces the time consumption to a certain extent, and has better objectivity for the result. In the selection process, the accuracy of the control point selected by the man-machine interaction and the reliability of the identification technology are used as the determining factors.
(3) Automatically selecting control points: the process of selecting the control points does not need human-computer interaction and manual participation, and the control points are found out completely by means of automatic identification of a computer, so that a position mapping model is established. At present, the automatic selection method of the control point is more researched in the field of remote sensing images, and is less involved in the field of multi-website POI. Aiming at multi-website POI data, control points are mainly selected automatically through an identification technology, but the positions of the control points matched by the method are often unevenly distributed, so that the local deformation is caused by the pathological condition of a law equation coefficient matrix in the position mapping calculation process; in a large geographic range, the position mapping is difficult to achieve the global optimal effect only by means of identification technology without combination of methods such as region division and the like, and the requirement of practical application on POI precision cannot be met.
2. Position mapping method
Common position mapping methods in the prior art include an average displacement method, a similarity transformation method, a polynomial method, and the like:
(1) average displacement method: the position mapping relationship is established by adding displacement amounts to the geographical positions of different data sources, and the displacement amounts generally take the average displacement of a common point. The mean displacement model is expressed in matrix form as follows:
wherein, r, c and x, y are the position information of different source POIs.
The average displacement method needs at least 1 control point to map the position, is suitable for data with a smaller geographic range of a data source and more obvious system error relative to accidental error, but has poor data mapping effect on large geographic range, accidental error and system error, and the effect of increasing the number of control points on improving the position mapping effect is not great.
(2) Similarity transformation method: transformations in euclidean space that transform a graph into a similar graph, such as translation, rotation, and scaling. The similarity transformation model is expressed in matrix form as follows:
wherein, r, c and x, y are the position information of different source POIs.
The similarity transformation method needs at least 2 control points to carry out position mapping, has the advantages of definite geometric relationship, regular transformation formula and proper extrapolation, can partially improve the position mapping effect by increasing the data quantity of the control points, and has poor data mapping effect on large geographic range, accidental errors and system errors.
(3) Polynomial transformation: and expressing the position mapping relation between different geographic positions by using a polynomial.
The second order polynomial transformation model is as follows:
wherein, r, c and x, y are the position information of different source POIs.
The second order polynomial transformation method requires at least 6 control points to perform the position mapping. The position mapping effect can be improved by increasing the number of control points, the data mapping effect which is not obvious to accidental errors and system errors is good, and the global optimal mapping effect cannot be achieved by only adopting a polynomial transformation method for a large geographic range.
In summary, the existing method for establishing the multi-website POI position mapping relationship still stays in a pure manual or semi-automatic stage, the method is long in time consumption, low in efficiency and high in cost, and for a POI point set with a large geographic range, only a common position mapping method is adopted without being combined with a matching technology, so that the position mapping cannot achieve the global optimal effect. Therefore, how to automatically and uniformly select control points and establish a position mapping relationship of POIs between network stations becomes a technical problem to be solved in the prior art.
Disclosure of Invention
The invention aims to solve the problem that the establishment of the prior multi-website POI position mapping relation cannot be completely automated, and the common position mapping method cannot achieve the global optimal effect on a POI point set with a large geographic range. And establishing a position mapping relation of the POI of the multi-website by combining a semantic matching method and a least square principle.
In order to achieve the purpose, the invention adopts the following technical scheme:
a multi-website POI position mapping method comprises the following steps:
s110, establishing a correction unit set:
respectively calculating minimum and maximum geographic coordinates according to a geographic coordinate set of the POI point set, wherein the obtained minimum and maximum geographic coordinates are the geographic range of the POI point set, dividing the geographic range of the POI point set into N x N small grids by utilizing a certain geographic interval, and forming a correction unit set by the N x N grid full set;
s120, performing region splitting in each correction unit of the correction unit set:
after the correction unit set is established, in order to ensure the uniform distribution of control points in each correction unit of the correction unit set, uniformly splitting each correction unit into a plurality of coordinate sub-grids to form a coordinate sub-grid of the control points, and in order to ensure the uniform distribution of check points in each correction unit of the correction unit set, further uniformly splitting each correction unit into a plurality of coordinate sub-grids to form a coordinate sub-grid of the check points;
s130, semantic-based control point matching:
carrying out word segmentation on POI names from different websites by using a Chinese word segmentation tool, and judging whether the POI names are the same control point or not according to word group strings after word segmentation;
s140, automatically and uniformly acquiring a multi-website POI homonymous control point set and a check point set:
in each control point sub-grid range, utilizing a semantic matching method to automatically acquire control points and add a control point set, and in each check point sub-grid range, utilizing a semantic matching method to automatically acquire at least 1 check point and add a check point set;
s150, resolving the position mapping model coefficient
And (4) utilizing the control points and the check point set obtained in the steps, selecting a second-order polynomial transformation model, and respectively resolving the position mapping model coefficients for each correction unit.
Wherein the second order polynomial transformation model is expressed in a matrix form by equation (1) as:
Where r, c is the geographic coordinate of a certain website POI, x, y is the geographic coordinate of another website POI, m0、m1、m2、m3、m4、m5、n0、n1、n2、n3、n4、n5Is the location mapping coefficient in a second order polynomial model, the geographic coordinates are longitude and latitude.
Preferably, in the step of establishing a set of correction units, the geographical interval is 4 °, 2 °, 1 °, 30 'or 15'.
Preferably, in the step of splitting the region in each correction unit of the correction unit set, each correction unit is split into 6 × 6 coordinate sub-grids to form coordinate sub-grids of the control points, and each correction unit is split into 12 × 12 coordinate sub-grids to form coordinate sub-grids of the check points.
Preferably, in the step of solving the coefficients of the position mapping model, the error equations are firstly listed according to a second-order polynomial transformation model, and then the position mapping coefficients in the error equations are calculated according to an indirect adjustment method.
Preferably, after S150, there is also S160: location mapping accuracy assessment
And traversing each correction unit in the correction unit set, calculating the median error, the maximum error and the minimum error of the check point set by using the position mapping coefficient delta obtained by resolving each correction unit, and finishing the evaluation of the multi-website POI position mapping precision.
The invention also discloses a multi-website POI position mapping device, which comprises the following units:
correction unit set creation unit (210):
respectively calculating minimum and maximum geographic coordinates according to a geographic coordinate set of the POI point set, wherein the obtained minimum and maximum geographic coordinates are the geographic range of the POI point set, dividing the geographic range of the POI point set into N x N small grids by utilizing a certain geographic interval, and forming a correction unit set by the N x N grid full set;
region splitting unit (220):
after the correction unit set is established, in order to ensure the uniform distribution of control points in each correction unit of the correction unit set, uniformly splitting each correction unit into a plurality of coordinate sub-grids to form a coordinate sub-grid of the control points, and in order to ensure the uniform distribution of check points in each correction unit of the correction unit set, further uniformly splitting each correction unit into a plurality of coordinate sub-grids to form a coordinate sub-grid of the check points;
semantic-based control point matching unit (230):
carrying out word segmentation on POI names from different websites by using a Chinese word segmentation tool, and judging whether the POI names are the same control point or not according to word group strings after word segmentation;
a control point set and check point set acquisition unit (240):
in each control point sub-grid range, utilizing a semantic matching method to automatically acquire control points and add a control point set, and in each check point sub-grid range, utilizing a semantic matching method to automatically acquire at least 1 check point and add a check point set;
position mapping model coefficient calculating unit (250)
And respectively resolving a position mapping model coefficient for each correction unit by using the control point and the check point set obtained by the units and selecting a second-order polynomial transformation model.
Wherein the second order polynomial transformation model is expressed in a matrix form by equation (1) as:
Where r, c is the geographic coordinate of a certain website POI, x, y is the geographic coordinate of another website POI, m0、m1、m2、m3、m4、m5、n0、n1、n2、n3、n4、n5Is the location mapping coefficient in a second order polynomial model, the geographic coordinates are longitude and latitude.
Preferably, in the establishing of the correction unit set unit, the geographical interval is 4 °, 2 °, 1 °, 30 'or 15'.
Preferably, in the region splitting unit, each correction unit is split into 6 × 6 coordinate sub-grids to form coordinate sub-grids of the control points, and each correction unit is split into 12 × 12 coordinate sub-grids to form coordinate sub-grids of the check points.
Preferably, in the position mapping model coefficient solving unit, the error equations are first listed according to a second-order polynomial transformation model, and then the position mapping coefficients in the error equations are calculated according to an indirect adjustment method.
Preferably, the method further comprises the following steps: position mapping accuracy evaluation unit (260)
And traversing each correction unit in the correction unit set, calculating the median error, the maximum error and the minimum error of the check point set by using the position mapping coefficient delta obtained by resolving each correction unit, and finishing the evaluation of the multi-website POI position mapping precision.
The method uniformly divides the POI position into a plurality of units, uniformly selects the control points, fully automatically selects the control points by comprehensively utilizing a semantic matching and least square method, and establishes an optimal position mapping model by applying a divide-and-conquer strategy.
Drawings
FIG. 1 is a flow diagram of a multi-site POI location mapping method in accordance with an embodiment of the present invention;
FIG. 2 is a diagram illustrating a correction unit set and a split of control point and check point regions in the correction units according to an embodiment of the present invention;
FIG. 3 is an automatic route map obtained by a multi-website POI homonymous control point set;
FIG. 4 is an automatic route map obtained by a multi-website POI homonymous check point set;
FIG. 5 is a block diagram of a multi-site POI location mapping device in accordance with an embodiment of the present invention;
FIG. 6 is a diagram of the result of location mapping for a multi-site POI.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings and examples. It is to be understood that the specific embodiments described herein are merely illustrative of the invention and are not limiting of the invention. It should be further noted that, for the convenience of description, only some of the structures related to the present invention are shown in the drawings, not all of the structures.
The principle of the invention is as follows: firstly, a correction unit set is established according to the geographical range of the multi-website POI. And obtaining a control point set and a check point set by a semantic matching method through region splitting in each correction unit. And resolving a position mapping coefficient by using a control point set according to a least square principle, and finally evaluating the position mapping precision.
Referring to fig. 1, a flow chart of a multi-site POI location mapping method according to the present invention is disclosed.
Specifically, the method comprises the following steps:
s110: establishing a set of correction units
And respectively calculating minimum and maximum geographic coordinates according to the geographic coordinate set of the POI point set, wherein the obtained minimum and maximum geographic coordinates are the geographic range of the POI point set, dividing the geographic range of the POI point set into N x N small grids by utilizing a certain geographic interval, and forming a correction unit set by the N x N grid full set.
Each POI contains geographic coordinates, and the whole set of the geographic coordinates of the POI in the POI point set forms the geographic coordinate set of the POI point set.
And respectively calculating the minimum and maximum geographic coordinates according to the geographic coordinate set of the POI point set, wherein the obtained minimum and maximum geographic coordinates are the geographic range of the POI point set.
In this embodiment, the beijing area (39.26 ° N, 115.25 ° E, 41.03 ° N, 117.30 ° E) is used as a test area, and the test data are two different sets of POI data sets of websites.
The Beijing 39.26 degrees and 41.03 degrees are minimum and maximum geographic coordinates of the Beijing POI point set in the latitude direction, and form a geographic range of the POI point set in the latitude direction; the east longitude 115.25 ° and 117.30 ° are minimum and maximum geographic coordinates of the beijing POI point set in the longitude direction, and constitute the geographic range of the POI point set in the longitude direction.
The geographical range of the set of POI points is divided into a small grid of N x N with a certain geographical interval. The smaller the geographical interval, the larger the value of N. Each small grid corresponds to one correction unit, and the full grid set of N x N forms a correction unit set.
In this embodiment, the correction unit sets can be respectively established at five geographical intervals of 4 °, 2 °, 1 °, 30 ', and 15'.
A typical geogrid grid spacing is shown in table 1.
TABLE 1 geography grid spacing
S120: performing region splitting within each correction unit of a set of correction units
After the correction unit set is established, in order to ensure the uniform distribution of the control points in each correction unit of the correction unit set, each correction unit is uniformly split into a plurality of coordinate sub-grids to form the coordinate sub-grids of the control points, and in order to ensure the uniform distribution of the check points in each correction unit of the correction unit set, each correction unit is also uniformly split into a plurality of coordinate sub-grids to form the coordinate sub-grids of the check points.
In a preferred embodiment, each correction unit is uniformly split into 6 × 6 coordinate sub-grids, which form coordinate sub-grids of control points. Each correction unit is evenly split into 12 × 12 coordinate subgrids, which form coordinate subgrids of the check points, as shown in fig. 2.
S130: semantic-based control point matching
And (3) segmenting the POI names from different websites by using a Chinese word segmentation tool, and judging whether the POI names are the same control point or not according to the word group string after word segmentation.
For example, in this embodiment, a POI name a of a certain website is the national mapping information bureau, and a POI name b of another website is the national mapping bureau. The phrase strings after word segmentation of the POI name a of one website are { { nation }, { mapping }, { information }, { office } }, and the phrase strings after word segmentation of the POI name b of another website are { { nation }, { mapping }, { office } }, both refer to the State geographic information office of surveying and mapping, and only the names are inconsistent, 2 can be considered as the same location and serve as a control point object.
S140: multi-website POI homonymous control point set and check point set are automatically and uniformly acquired
In the sub-grid range of each control point, a semantic matching method is used to automatically acquire the control points and add the control point set, as shown in fig. 4.
Similarly, in each checking point sub-grid range, at least 1 checking point is automatically obtained by using a semantic matching method, and is added into the checking point set, as shown in fig. 5.
This step is to use the method of steps S110-S130 to traverse all the control point sets and check point sets, and ensure that each grid obtains the corresponding control point and check point.
Those skilled in the art will appreciate that step S130 may be directly merged into step S140, and each grid obtains the corresponding control point and check point by using a matching method directly in each control point sub-grid and each check point sub-grid.
S150: resolving location mapping model coefficients
The position model coefficients are solved using the control points.
Establishing a position mapping model among different website POIs is equivalent to establishing a position mapping relation, and the position mapping model is expressed by position mapping model coefficients. Therefore, the position mapping model coefficients are determined through calculation, and the position mapping model among different website POIs can be well established.
The control points are used for calculating position mapping coefficients, and the check points are used for evaluating the model reliability of the position mapping relation.
The calculation of the position mapping model coefficients relies on sufficient control points for calculating the position mapping coefficients and check points for evaluating the accuracy of the position mapping.
In a specific embodiment, the control point and the check point set obtained in the above steps are adopted, a second-order polynomial transformation model is selected, and the position mapping model coefficients are respectively solved for each correction unit.
Wherein the second order polynomial transformation model is expressed in a matrix form by equation (1) as:
Where r, c is the geographic coordinate of a certain website POI, i.e. longitude and latitude, x, y is the geographic coordinate of another website POI, i.e. longitude and latitude, m0、m1、m2、m3、m4、m5、n0、n1、n2、n3、n4、n5Is the position mapping coefficient in a second order polynomial model. The position mapping relation between two different websites can be established through the transformation model.
Specifically, the basic principle of solving the position mapping model is equivalent to an indirect adjustment method of measurement, firstly, an error equation is listed according to a second-order polynomial transformation model, and then undetermined parameters (position mapping coefficients) in the error equation are calculated according to the indirect adjustment method.
The error equation is relative to the homogeneous equation, the number of equations is more than the number of unknowns, and the most probable value of the unknowns is obtained according to the least square principle.
In this embodiment, a second-order polynomial model is used to obtain the most probable value of the position mapping coefficient, and the error equation is listed in formula (1) as follows:
Namely: v ═ B Δ -1
Wherein, <math>
<mrow>
<mi>B</mi>
<mo>=</mo>
<mfenced open='(' close=')'>
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mi>r</mi>
</mtd>
<mtd>
<mi>c</mi>
</mtd>
<mtd>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
</mtd>
<mtd>
<mi>r</mi>
<mo>×</mo>
<mi>c</mi>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mi>r</mi>
</mtd>
<mtd>
<mi>c</mi>
</mtd>
<mtd>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
</mtd>
<mtd>
<mi>r</mi>
<mo>×</mo>
<mi>c</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
</math>
mi=mi initial+Δmi
Wherein V comprises Vx、VyThe correction of the geographical coordinates of another website POI, namely longitude x and latitude y is represented, B is a coefficient matrix of an error equation, and B is obtained by calculation through the geographical coordinates, longitude r and latitude c of a certain website POI; l includes x and y, and is the geographical coordinates of another site POI, longitude x and latitude y. Position mapping coefficient miIs a solution of the position mapping coefficienti initialAnd the number of corrections Δ miAnd (4) forming. In general, the solution initial value m of the position mapping coefficienti initialAre all given 0, in which case mi=ΔmiTherefore,. DELTA.m is used hereinafteriTo represent the position mapping coefficient mi。
According to the calculation method of indirect adjustment, matrix operation is carried out:
Δ=(BTB)-1BT1
the position mapping coefficient Δ is solved.
The basic idea of this step is to establish an error equation based on the least squares adjustment method. The coordinates of the known control points are substituted in the equation to solve for the position mapping coefficients.
After the position mapping coefficient is solved, the POI coordinates are directly input, and the coordinates after the position mapping can be solved. For example, after the position mapping coefficient between two different websites has been solved, the geographical coordinates (i.e., longitude r and latitude c) of a certain website POI are input, and then the longitude coordinate x and the latitude coordinate y of another website POI are obtained by matrix operation using formula 1.
S160: location mapping accuracy assessment
And traversing each correction unit in the correction unit set, calculating the median error, the maximum error and the minimum error of the check point set by using the position mapping coefficient delta obtained by resolving each correction unit, and finishing the evaluation of the multi-website POI position mapping precision.
That is, each correction unit calculates a corresponding position coefficient, and the calculation error of the correction unit is calculated through the check point in each correction unit using the position coefficient of the correction unit. Those skilled in the art will appreciate that the method of calculating the error using the check point is any other conventional method in the art, for example, a conventional method of measuring the adjustment such as: conditional adjustment with parameters and indirect adjustment with limiting conditions; conventional methods for accuracy assessment are as follows: and calculating the error of the model by introducing a control point set.
In this embodiment, the accuracy is shown in table 2, the unit of error is km, and the result of the position mapping accuracy evaluation is shown in fig. 6.
TABLE 2 precision statistics table
In the embodiment, compared with the manual and semi-automatic mode of selecting the control points required by the current multi-website POI position mapping model establishment, the automatic and uniform selection of the control points is completely realized, and the working efficiency is greatly improved.
In this embodiment, as shown in table 2, a traditional method is adopted to establish a mapping relationship of POI positions of different websites, where a median error in a longitude direction is 1.25km, and a median error in a latitude direction is 0.074 km; by adopting the method (taking fifteen as geographical intervals as an example), the median error in the longitude direction is 0.88m, the median error in the latitude direction is 0.01m, the accuracy of position mapping is improved from kilometer level to decimeter level, and the requirement of practical application on the POI position mapping accuracy is met.
Referring to fig. 5, the present invention also discloses a position mapping apparatus for a multi-website POI, comprising the following units
The correction unit set creating unit 210:
respectively calculating minimum and maximum geographic coordinates according to a geographic coordinate set of the POI point set, wherein the obtained minimum and maximum geographic coordinates are the geographic range of the POI point set, dividing the geographic range of the POI point set into N x N small grids by utilizing a certain geographic interval, and forming a correction unit set by the N x N grid full set;
the region splitting unit 220:
after the correction unit set is established, in order to ensure the uniform distribution of control points in each correction unit of the correction unit set, uniformly splitting each correction unit into a plurality of coordinate sub-grids to form a coordinate sub-grid of the control points, and in order to ensure the uniform distribution of check points in each correction unit of the correction unit set, further uniformly splitting each correction unit into a plurality of coordinate sub-grids to form a coordinate sub-grid of the check points;
semantic-based control point matching unit 230:
carrying out word segmentation on POI names from different websites by using a Chinese word segmentation tool, and judging whether the POI names are the same control point or not according to word group strings after word segmentation;
control point set/check point set acquisition unit 240:
in each control point sub-grid range, utilizing a semantic matching method to automatically acquire control points and add a control point set, and in each check point sub-grid range, utilizing a semantic matching method to automatically acquire at least 1 check point and add a check point set;
the position mapping model coefficient solving unit 250:
and respectively resolving a position mapping model coefficient for each correction unit by using the control point and the check point set obtained by the units and selecting a second-order polynomial transformation model.
Wherein the second order polynomial transformation model is expressed in a matrix form by equation (1) as:
Where r, c is the geographic coordinate of a certain website POI, x, y is the geographic coordinate of another website POI, m0、m1、m2、m3、m4、m5、n0、n1、n2、n3、n4、n5Is the location mapping coefficient in a second order polynomial model, the geographic coordinates are longitude and latitude.
Wherein, in the unit of establishing the correction unit set, the geographic interval is 4 °, 2 °, 1 °, 30 'or 15'.
In the region splitting unit, each correction unit is uniformly split into 6 × 6 coordinate sub-grids to form coordinate sub-grids of control points, and each correction unit is uniformly split into 12 × 12 coordinate sub-grids to form coordinate sub-grids of check points.
In the position mapping model coefficient resolving unit, firstly listing an error equation according to a second-order polynomial transformation model, and then calculating a position mapping coefficient in the error equation according to an indirect adjustment method;
wherein, the most probable value of the position mapping coefficient is obtained by a second-order polynomial model, and the error equation is listed by the formula (1) as follows
Namely: v ═ B Δ -1
Wherein, <math>
<mrow>
<mi>B</mi>
<mo>=</mo>
<mfenced open='(' close=')'>
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mi>r</mi>
</mtd>
<mtd>
<mi>c</mi>
</mtd>
<mtd>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
</mtd>
<mtd>
<mi>r</mi>
<mo>×</mo>
<mi>c</mi>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mi>r</mi>
</mtd>
<mtd>
<mi>c</mi>
</mtd>
<mtd>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
</mtd>
<mtd>
<mi>r</mi>
<mo>×</mo>
<mi>c</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
</math>
mi=mi initial+Δmi
Wherein V comprises Vx、VyThe correction of the geographical coordinates of another website POI, namely longitude x and latitude y is represented, B is a coefficient matrix of an error equation, and B is obtained by calculation through the geographical coordinates, longitude r and latitude c of a certain website POI; i includes x and y, which are geographical coordinates of POI, longitude x and latitude y of another website. Position mapping coefficient miIs a solution of the position mapping coefficienti initialAnd the number of corrections Δ miAnd (4) forming. In general, the solution initial value m of the position mapping coefficienti initialAre all given 0, in which case mi=ΔmiTherefore,. DELTA.m is used hereinafteriTo represent the position mapping coefficient mi。
According to the calculation method of indirect adjustment, matrix operation is carried out:
Δ=(BTB)-1BT1
the position mapping coefficient Δ is solved.
Wherein, the device also comprises a position mapping precision evaluation unit (260)
And traversing each correction unit in the correction unit set, calculating the median error, the maximum error and the minimum error of the check point set by using the position mapping coefficient delta obtained by resolving each correction unit, and finishing the evaluation of the multi-website POI position mapping precision.
That is, each correction unit has a set of position coefficients, and the calculation error of the correction unit is calculated in each correction unit using the position coefficients of the correction unit and the check point.
The method uniformly divides the POI position into a plurality of units, uniformly selects the control points, fully automatically selects the control points by comprehensively utilizing a semantic matching and least square method, and establishes an optimal position mapping model by applying a divide-and-conquer strategy.
It will be apparent to those skilled in the art that the various elements or steps of the invention described above may be implemented using a general purpose computing device, they may be centralized on a single computing device, or alternatively, they may be implemented using program code that is executable by a computing device, such that they may be stored in a memory device and executed by a computing device, or they may be separately fabricated into various integrated circuit modules, or multiple ones of them may be fabricated into a single integrated circuit module. Thus, the present invention is not limited to any specific combination of hardware and software.
While the invention has been described in further detail with reference to specific preferred embodiments thereof, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention as defined by the appended claims.
Claims (10)
1. A multi-website POI position mapping method comprises the following steps:
s110, establishing a correction unit set:
respectively calculating minimum and maximum geographic coordinates according to a geographic coordinate set of the POI point set, wherein the obtained minimum and maximum geographic coordinates are the geographic range of the POI point set, dividing the geographic range of the POI point set into N x N small grids by utilizing a certain geographic interval, and forming a correction unit set by the N x N grid full set;
s120, performing region splitting in each correction unit of the correction unit set:
after the correction unit set is established, in order to ensure the uniform distribution of control points in each correction unit of the correction unit set, uniformly splitting each correction unit into a plurality of coordinate sub-grids to form a coordinate sub-grid of the control points, and in order to ensure the uniform distribution of check points in each correction unit of the correction unit set, further uniformly splitting each correction unit into a plurality of coordinate sub-grids to form a coordinate sub-grid of the check points;
s130, semantic-based control point matching:
carrying out word segmentation on POI names from different websites by using a Chinese word segmentation tool, and judging whether the POI names are the same control point or not according to word group strings after word segmentation;
s140, automatically and uniformly acquiring a multi-website POI homonymous control point set and a check point set:
in each control point sub-grid range, utilizing a semantic matching method to automatically acquire control points and add a control point set, and in each check point sub-grid range, utilizing a semantic matching method to automatically acquire at least 1 check point and add a check point set;
s150, resolving the position mapping model coefficient
Using the control points and the check point set obtained in the above steps, selecting a second-order polynomial transformation model, and respectively resolving a position mapping model coefficient for each correction unit;
wherein the second order polynomial transformation model is expressed in a matrix form by equation (1) as:
Where r, c is the geographic coordinate of a certain website POI, x, y is the geographic coordinate of another website POI, m0、m1、m2、m3、m4、m5、n0、n1、n2、n3、n4、n5Is the location mapping coefficient in a second order polynomial model, the geographic coordinates are longitude and latitude.
2. The multi-site POI location mapping method of claim 1, wherein:
in the step of establishing a correction unit set, the geographical interval is 4 °, 2 °, 1 °, 30 'or 15'.
3. The multi-site POI location mapping method of claim 1, wherein:
in the step of carrying out region splitting in each correction unit of the correction unit set, uniformly splitting each correction unit into 6 × 6 coordinate sub-grids to form coordinate sub-grids of the control points, and uniformly splitting each correction unit into 12 × 12 coordinate sub-grids to form coordinate sub-grids of the check points.
4. The multi-site POI location mapping method of claim 1, wherein:
in the step of calculating the coefficients of the position mapping model, firstly, listing an error equation according to a second-order polynomial transformation model, and then calculating the position mapping coefficients in the error equation according to an indirect adjustment method;
wherein, the most probable value of the position mapping coefficient is obtained by a second-order polynomial model, and the error equation is listed by the formula (1) as follows
Namely: v ═ B Δ -1
Wherein, <math>
<mrow>
<mi>B</mi>
<mo>=</mo>
<mfenced open='(' close=')'>
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mi>r</mi>
</mtd>
<mtd>
<mi>c</mi>
</mtd>
<mtd>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
</mtd>
<mtd>
<mi>r</mi>
<mo>×</mo>
<mi>c</mi>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mi>r</mi>
</mtd>
<mtd>
<mi>c</mi>
</mtd>
<mtd>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
</mtd>
<mtd>
<mi>r</mi>
<mo>×</mo>
<mi>c</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
</math>
mi=mi initial+Δmi
V comprises Vx、VyThe correction of the geographical coordinates of another website POI is expressed, B is a coefficient matrix of an error equation, and B is obtained by calculation of the geographical coordinates of a certain website POI; 1 includes x and y, is the geographical coordinates of another site POI, mi initialIndicating the solved initial value of the position mapping coefficient, Δ miA calculation correction number representing a position mapping coefficient;
according to the calculation method of indirect adjustment, matrix operation is carried out:
Δ=(BTB)-1BT1
the position mapping coefficient Δ is solved.
5. The multi-site POI location mapping method of claim 1, wherein:
after S150, there is also S160: and (4) evaluating the precision of the position mapping,
and traversing each correction unit in the correction unit set, calculating the median error, the maximum error and the minimum error of the check point set by using the position mapping coefficient delta obtained by resolving each correction unit, and finishing the evaluation of the multi-website POI position mapping precision.
6. A multi-website POI position mapping device comprises the following units:
correction unit set creation unit (210):
respectively calculating minimum and maximum geographic coordinates according to a geographic coordinate set of the POI point set, wherein the obtained minimum and maximum geographic coordinates are the geographic range of the POI point set, dividing the geographic range of the POI point set into N x N small grids by utilizing a certain geographic interval, and forming a correction unit set by the N x N grid full set;
region splitting unit (220):
after the correction unit set is established, in order to ensure the uniform distribution of control points in each correction unit of the correction unit set, uniformly splitting each correction unit into a plurality of coordinate sub-grids to form a coordinate sub-grid of the control points, and in order to ensure the uniform distribution of check points in each correction unit of the correction unit set, further uniformly splitting each correction unit into a plurality of coordinate sub-grids to form a coordinate sub-grid of the check points;
semantic-based control point matching unit (230):
carrying out word segmentation on POI names from different websites by using a Chinese word segmentation tool, and judging whether the POI names are the same control point or not according to word group strings after word segmentation;
a control point set and check point set acquisition unit (240):
in each control point sub-grid range, utilizing a semantic matching method to automatically acquire control points and add a control point set, and in each check point sub-grid range, utilizing a semantic matching method to automatically acquire at least 1 check point and add a check point set;
position mapping model coefficient calculating unit (250)
Using the control points and the check point sets obtained by the units, selecting a second-order polynomial transformation model, and respectively resolving a position mapping model coefficient for each correction unit;
wherein the second order polynomial transformation model is expressed in a matrix form by equation (1) as:
Where r, c is the geographic coordinate of a certain website POI, x, y is the geographic coordinate of another website POI, m0、m1、m2、m3、m4、m5、n0、n1、n2、n3、n4、n5Is the location mapping coefficient in a second order polynomial model, the geographic coordinates are longitude and latitude.
7. The multi-website POI location mapping apparatus of claim 6, wherein:
in the cells for establishing the correction cell set, the geographical interval is 4 °, 2 °, 1 °, 30 'or 15'.
8. The multi-website POI location mapping apparatus of claim 6, wherein:
in the region splitting unit, each correction unit is uniformly split into 6 × 6 coordinate sub-grids to form coordinate sub-grids of control points, and each correction unit is uniformly split into 12 × 12 coordinate sub-grids to form coordinate sub-grids of check points.
9. The multi-website POI location mapping apparatus of claim 6, wherein:
in a position mapping model coefficient resolving unit, firstly, listing an error equation according to a second-order polynomial transformation model, and then calculating a position mapping coefficient in the error equation according to an indirect adjustment method;
wherein, the most probable value of the position mapping coefficient is obtained by a second-order polynomial model, and the error equation is listed by the formula (1) as follows
Namely: v ═ B Δ -1
Wherein, <math>
<mrow>
<mi>B</mi>
<mo>=</mo>
<mfenced open='(' close=')'>
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mi>r</mi>
</mtd>
<mtd>
<mi>c</mi>
</mtd>
<mtd>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
</mtd>
<mtd>
<mi>r</mi>
<mo>×</mo>
<mi>c</mi>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mi>r</mi>
</mtd>
<mtd>
<mi>c</mi>
</mtd>
<mtd>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
</mtd>
<mtd>
<mi>r</mi>
<mo>×</mo>
<mi>c</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
</math>
mi=mi initial+Δmi
V comprises Vx、VyThe correction of the geographical coordinates of another website POI is expressed, B is a coefficient matrix of an error equation, and B is obtained by calculation of the geographical coordinates of a certain website POI; 1 includes x and y, is the geographical coordinates of another site POI, mi initialIndicating the solved initial value of the position mapping coefficient, Δ miA calculation correction number representing a position mapping coefficient;
according to the calculation method of indirect adjustment, matrix operation is carried out:
Δ=(BTB)-1BT1
the position mapping coefficient Δ is solved.
10. The multi-website POI location mapping apparatus of claim 1, wherein:
further comprising: position mapping accuracy evaluation unit (260)
And traversing each correction unit in the correction unit set, calculating the median error, the maximum error and the minimum error of the check point set by using the position mapping coefficient delta obtained by resolving each correction unit, and finishing the evaluation of the multi-website POI position mapping precision.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410455721.1A CN104199937A (en) | 2014-09-10 | 2014-09-10 | Multi-website POI position mapping method and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410455721.1A CN104199937A (en) | 2014-09-10 | 2014-09-10 | Multi-website POI position mapping method and device |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104199937A true CN104199937A (en) | 2014-12-10 |
Family
ID=52085230
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410455721.1A Pending CN104199937A (en) | 2014-09-10 | 2014-09-10 | Multi-website POI position mapping method and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104199937A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104572902A (en) * | 2014-12-26 | 2015-04-29 | 北京中交兴路车联网科技有限公司 | Information point matching method and device |
CN104850649A (en) * | 2015-05-29 | 2015-08-19 | 苏州大学张家港工业技术研究院 | Method and system for sampling points of interest on map |
CN110096540A (en) * | 2019-04-16 | 2019-08-06 | 湖北地信科技集团股份有限公司 | Surveying and mapping data conversion method, equipment, storage medium and device |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103124966A (en) * | 2010-07-09 | 2013-05-29 | 诺基亚公司 | Method and apparatus for aggregating and linking place data |
CN103150309A (en) * | 2011-12-07 | 2013-06-12 | 清华大学 | Method and system for searching POI (Point of Interest) points of awareness map in space direction |
-
2014
- 2014-09-10 CN CN201410455721.1A patent/CN104199937A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103124966A (en) * | 2010-07-09 | 2013-05-29 | 诺基亚公司 | Method and apparatus for aggregating and linking place data |
CN103150309A (en) * | 2011-12-07 | 2013-06-12 | 清华大学 | Method and system for searching POI (Point of Interest) points of awareness map in space direction |
Non-Patent Citations (1)
Title |
---|
陈瑞等: "一种自动提取控制点的多源兴趣点位置纠正方法", 《测绘科学》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104572902A (en) * | 2014-12-26 | 2015-04-29 | 北京中交兴路车联网科技有限公司 | Information point matching method and device |
CN104572902B (en) * | 2014-12-26 | 2018-01-23 | 北京中交兴路车联网科技有限公司 | A kind of method and device of information Point matching |
CN104850649A (en) * | 2015-05-29 | 2015-08-19 | 苏州大学张家港工业技术研究院 | Method and system for sampling points of interest on map |
CN104850649B (en) * | 2015-05-29 | 2018-06-19 | 苏州大学张家港工业技术研究院 | A kind of method and system that point of interest sampling is carried out on map |
CN110096540A (en) * | 2019-04-16 | 2019-08-06 | 湖北地信科技集团股份有限公司 | Surveying and mapping data conversion method, equipment, storage medium and device |
CN110096540B (en) * | 2019-04-16 | 2022-02-18 | 湖北地信科技集团股份有限公司 | Mapping data conversion method, device, storage medium and device |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109376761B (en) | Address identification and longitude and latitude mining method and device thereof | |
US11512970B2 (en) | Method and apparatus for generating information | |
CN106919567B (en) | Method and device for processing POI (Point of interest) address | |
JP2017016625A (en) | Method and apparatus for determining location of new point of interest | |
Sehra et al. | A systematic study of OpenStreetMap data quality assessment | |
KR100686287B1 (en) | Distorting Modeling method for Transforming the Presize Position of Partial/Positional information | |
US20180285371A1 (en) | Method and apparatus for automatically discovering gas station poi, storage medium and device | |
CN115796712B (en) | Regional land ecosystem carbon reserve estimation method and device and electronic equipment | |
CN106323266B (en) | Method and device for processing information of interest point | |
JP2010128898A (en) | Content position estimating device | |
CN104199937A (en) | Multi-website POI position mapping method and device | |
CN111177289A (en) | Method and system for extracting and checking related information of data space of multi-source network | |
CN103177189B (en) | Register Data Quality Analysis method in a kind of many source positions | |
CN103064953B (en) | Position-based is registered the POI update method of data | |
CN105426387B (en) | Map aggregation method based on K-means algorithm | |
ESTOQUE et al. | < Original Papers> Validating ALOS PRISM DSM-derived surface feature height: Implications for urban volume estimation | |
CN111986552B (en) | Map data quality freshness acquisition method, device and storage medium | |
CN106291756B (en) | The construction method of near space air virtual environment resource | |
CN107203547B (en) | Method and device for auditing point of interest information | |
CN117197639A (en) | True value acquisition method and device, electronic equipment and storage medium | |
CN116629265A (en) | Address verification method and device, electronic equipment and storage medium | |
CN103902633A (en) | Method, device and system for generating POI hotspot data | |
Krdžalić et al. | A precise geoid model of Bosnia and Herzegovina by the KTH method and its validation | |
CN104809336A (en) | Method for sampling region factor by considering spatial correlation | |
CN113590544B (en) | Open-air land cover sample sampling auxiliary system based on photo EXIF information reads |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20141210 |
|
RJ01 | Rejection of invention patent application after publication |