CN104199937A - Multi-website POI position mapping method and device - Google Patents

Multi-website POI position mapping method and device Download PDF

Info

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
Application number
CN201410455721.1A
Other languages
Chinese (zh)
Inventor
刘纪平
王勇
杨毅
罗安
张福浩
王亮
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chinese Academy of Surveying and Mapping
Original Assignee
Chinese Academy of Surveying and Mapping
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chinese Academy of Surveying and Mapping filed Critical Chinese Academy of Surveying and Mapping
Priority to CN201410455721.1A priority Critical patent/CN104199937A/en
Publication of CN104199937A publication Critical patent/CN104199937A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/90Details of database functions independent of the retrieved data types
    • G06F16/95Retrieval from the web
    • G06F16/953Querying, e.g. by the use of web search engines
    • G06F16/9537Spatial or temporal dependent retrieval, e.g. spatiotemporal queries
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/29Geographical 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

Multi-website POI position mapping method and device
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:
x = r + m 0 y = c + n 0
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:
<math> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mi>x</mi> <mo>=</mo> <msub> <mi>m</mi> <mn>1</mn> </msub> <mo>&times;</mo> <mi>r</mi> <mo>+</mo> <msub> <mi>m</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mi>y</mi> <mo>=</mo> <msub> <mi>n</mi> <mn>1</mn> </msub> <mo>&times;</mo> <mi>c</mi> <mo>+</mo> <msub> <mi>n</mi> <mn>0</mn> </msub> </mtd> </mtr> </mtable> </mfenced> </math>
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:
<math> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mi>x</mi> <mo>=</mo> <msub> <mi>m</mi> <mn>1</mn> </msub> <mo>&times;</mo> <mi>r</mi> <mo>+</mo> <msub> <mi>m</mi> <mn>2</mn> </msub> <mo>&times;</mo> <mi>c</mi> <mo>+</mo> <msub> <mi>m</mi> <mn>3</mn> </msub> <mo>&times;</mo> <msup> <mi>r</mi> <mn>2</mn> </msup> <mo>+</mo> <msub> <mi>m</mi> <mn>4</mn> </msub> <mo>&times;</mo> <msup> <mi>c</mi> <mn>2</mn> </msup> <mo>+</mo> <msub> <mi>m</mi> <mn>5</mn> </msub> <mo>&times;</mo> <mi>r</mi> <mo>&times;</mo> <mi>c</mi> <mo>+</mo> <msub> <mi>m</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mi>y</mi> <mo>=</mo> <msub> <mi>n</mi> <mn>1</mn> </msub> <mo>&times;</mo> <mi>r</mi> <mo>+</mo> <msub> <mi>n</mi> <mn>2</mn> </msub> <mo>&times;</mo> <mi>c</mi> <mo>+</mo> <msub> <mi>n</mi> <mn>3</mn> </msub> <mo>&times;</mo> <msup> <mi>r</mi> <mn>2</mn> </msup> <mo>+</mo> <msub> <mi>n</mi> <mn>4</mn> </msub> <mo>&times;</mo> <msup> <mi>c</mi> <mn>2</mn> </msup> <mo>+</mo> <msub> <mi>n</mi> <mn>5</mn> </msub> <mo>&times;</mo> <mi>r</mi> <mo>&times;</mo> <mi>c</mi> <mo>+</mo> <msub> <mi>n</mi> <mn>0</mn> </msub> </mtd> </mtr> </mtable> </mfenced> </math>
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:
<math> <mrow> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>x</mi> </mtd> </mtr> <mtr> <mtd> <mi>y</mi> </mtd> </mtr> </mtable> </mfenced> <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>&times;</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>&times;</mo> <mi>c</mi> </mtd> </mtr> </mtable> </mfenced> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>m</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>4</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>5</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>4</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>5</mn> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow> </math> formula (1)
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:
<math> <mrow> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>x</mi> </mtd> </mtr> <mtr> <mtd> <mi>y</mi> </mtd> </mtr> </mtable> </mfenced> <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>&times;</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>&times;</mo> <mi>c</mi> </mtd> </mtr> </mtable> </mfenced> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>m</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>4</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>5</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>4</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>5</mn> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow> </math> formula (1)
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:
<math> <mrow> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>x</mi> </mtd> </mtr> <mtr> <mtd> <mi>y</mi> </mtd> </mtr> </mtable> </mfenced> <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>&times;</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>&times;</mo> <mi>c</mi> </mtd> </mtr> </mtable> </mfenced> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>m</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>4</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>5</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>4</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>5</mn> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow> </math> formula (1)
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:
<math> <mrow> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>v</mi> <mi>x</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>v</mi> <mi>y</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <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>&times;</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>&times;</mo> <mi>c</mi> </mtd> </mtr> </mtable> </mfenced> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>4</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>5</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>4</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>5</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>x</mi> </mtd> </mtr> <mtr> <mtd> <mi>y</mi> </mtd> </mtr> </mtable> </mfenced> </mrow> </math> formula (2)
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>&times;</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>&times;</mo> <mi>c</mi> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow> </math>
l = x y
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:
<math> <mrow> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>x</mi> </mtd> </mtr> <mtr> <mtd> <mi>y</mi> </mtd> </mtr> </mtable> </mfenced> <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>&times;</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>&times;</mo> <mi>c</mi> </mtd> </mtr> </mtable> </mfenced> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>m</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>4</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>5</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>4</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>5</mn> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow> </math> formula (1)
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
<math> <mrow> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>v</mi> <mi>x</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>v</mi> <mi>y</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <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>&times;</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>&times;</mo> <mi>c</mi> </mtd> </mtr> </mtable> </mfenced> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>4</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>5</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>4</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>5</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>x</mi> </mtd> </mtr> <mtr> <mtd> <mi>y</mi> </mtd> </mtr> </mtable> </mfenced> </mrow> </math> Formula (2)
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>&times;</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>&times;</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:
<math> <mrow> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>x</mi> </mtd> </mtr> <mtr> <mtd> <mi>y</mi> </mtd> </mtr> </mtable> </mfenced> <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>&times;</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>&times;</mo> <mi>c</mi> </mtd> </mtr> </mtable> </mfenced> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>m</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>4</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>5</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>4</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>5</mn> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow> </math> formula (1)
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
<math> <mrow> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>v</mi> <mi>x</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>v</mi> <mi>y</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <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>&times;</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>&times;</mo> <mi>c</mi> </mtd> </mtr> </mtable> </mfenced> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>4</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>5</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>4</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>5</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>x</mi> </mtd> </mtr> <mtr> <mtd> <mi>y</mi> </mtd> </mtr> </mtable> </mfenced> </mrow> </math> Formula (2)
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>&times;</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>&times;</mo> <mi>c</mi> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow> </math>
l = x y
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:
<math> <mrow> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>x</mi> </mtd> </mtr> <mtr> <mtd> <mi>y</mi> </mtd> </mtr> </mtable> </mfenced> <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>&times;</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>&times;</mo> <mi>c</mi> </mtd> </mtr> </mtable> </mfenced> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>m</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>4</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>m</mi> <mn>5</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>4</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mn>5</mn> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow> </math> formula (1)
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
<math> <mrow> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>v</mi> <mi>x</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>v</mi> <mi>y</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <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>&times;</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>&times;</mo> <mi>c</mi> </mtd> </mtr> </mtable> </mfenced> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>4</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;m</mi> <mn>5</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>4</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&Delta;n</mi> <mn>5</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>x</mi> </mtd> </mtr> <mtr> <mtd> <mi>y</mi> </mtd> </mtr> </mtable> </mfenced> </mrow> </math> Formula (2)
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>&times;</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>&times;</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.
CN201410455721.1A 2014-09-10 2014-09-10 Multi-website POI position mapping method and device Pending CN104199937A (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
陈瑞等: "一种自动提取控制点的多源兴趣点位置纠正方法", 《测绘科学》 *

Cited By (6)

* Cited by examiner, † Cited by third party
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