CN102324096B - Method and device for splicing characteristic point correction partitions - Google Patents

Method and device for splicing characteristic point correction partitions Download PDF

Info

Publication number
CN102324096B
CN102324096B CN201110201465.XA CN201110201465A CN102324096B CN 102324096 B CN102324096 B CN 102324096B CN 201110201465 A CN201110201465 A CN 201110201465A CN 102324096 B CN102324096 B CN 102324096B
Authority
CN
China
Prior art keywords
area
square error
image
unique point
mean square
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201110201465.XA
Other languages
Chinese (zh)
Other versions
CN102324096A (en
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.)
State Grid Corp of China SGCC
State Grid Economic and Technological Research Institute
Beijing Electric Power Corp
Original Assignee
State Grid Corp of China SGCC
State Grid Economic and Technological Research Institute
Beijing Electric Power Corp
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 State Grid Corp of China SGCC, State Grid Economic and Technological Research Institute, Beijing Electric Power Corp filed Critical State Grid Corp of China SGCC
Priority to CN201110201465.XA priority Critical patent/CN102324096B/en
Publication of CN102324096A publication Critical patent/CN102324096A/en
Priority to US14/117,405 priority patent/US9256981B2/en
Priority to EP12785279.6A priority patent/EP2711893B1/en
Priority to PCT/CN2012/070646 priority patent/WO2012155540A1/en
Application granted granted Critical
Publication of CN102324096B publication Critical patent/CN102324096B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention provides a method for splicing characteristic point correction partitions and is used for solving the problem of digital convergence of historical precious drawings in the prior art. The method comprises the following steps of: 2, acquiring first image data, and selecting a characteristic point group of a first area on a first image according to a preset coordinate; 4, correcting each characteristic point in the characteristic point group, and thus forming a second area; and 6, splicing the second area at a position, corresponding to the first area, on a second image. By adoption of the technical scheme of the invention, the problem of the digital convergence of various kinds of paper drawings which are accumulated over the years and include paper geologic maps with different time, different dryness and humidity and different scales can be effectively solved, the paper drawings can be accurately superposed and converged with geologic, geographic and grid vector data, and use requirements of stages such as scanning, designing, construction, running and maintenance on geologic information are met.

Description

Characteristic point correction partitions joining method and device
Technical field
The present invention relates to a kind of characteristic point correction partitions joining method and device.
Background technology
Along with urban construction, geologic condition becomes increasingly complex, and expropriation of land removal cost is more and more higher, and power engineering addressing pressure is huge.And traditional figure is board-like, distributing geologic information management mode, work efficiency is low, and data utilization factor is low, is more and more not suitable with the demand of Informatization Development.Sweep vectorization be exactly by geological information drawing through overscanning, and lines in figure and print hand writing are identified, and are converted into vector lines and vector word, so that with electronic chart convergent-divergent.Because the precious drawing of history is in preserving for a long time, due to folding, extruding, cause the deformation of drawing, after carrying out drawing scanning, produce certain error.
There is in the prior art the problem that vector scan is existed to the risk of information distortion, and for this problem, not yet propose at present effective solution.
Summary of the invention
Fundamental purpose of the present invention is to provide a kind of characteristic point correction partitions joining method and device, to solve the problem that aligns vector scan in prior art and exist the risk of information distortion.
To achieve these goals, according to an aspect of the present invention, provide a kind of characteristic point correction partitions joining method.
Characteristic point correction partitions of the present invention joining method comprises: S2: gather the first view data, according to the characteristic point group of the first area on coordinate selection first image of existing electric power facility; S4: each unique point in characteristic point group is proofreaied and correct, form second area; And S6: second area is spliced on the position of corresponding first area on the second image.
Further, step S4 comprises: S41: the Minimum Mean Square Error that calculates each unique point; And S43: each unique point is proofreaied and correct according to Minimum Mean Square Error.
Further, in step S41, the computing method of Minimum Mean Square Error are: the basic point coordinate of each unique point is obtained to the Minimum Mean Square Error of each unique point by rigid transformation, the formula that Minimum Mean Square Error adopts is:
Figure BDA0000076743620000011
wherein, x i, y ifor the basic point coordinate of unique point;
Figure BDA0000076743620000012
for least mean-square error; F (x i) be rigid transformation, in piece image, the distance between any two points remains unchanged before and after conversion.
Further, in step S41, the computing method of Minimum Mean Square Error are: the basic point coordinate of each unique point is obtained to the Minimum Mean Square Error of each unique point by rotational transform, the formula that Minimum Mean Square Error adopts is
Figure BDA0000076743620000021
wherein, x i, y ifor the basic point coordinate of unique point;
Figure BDA0000076743620000022
for least mean-square error; R (x i) be rotational transform, allow every bit P determine angle around one of a point of fixity rotation, become another P ', the conversion so producing is called the rotational transform in plane; T is translation transformation, between image along x, the translational movement in y direction.
Further, the method in step S6, second area being spliced on the second image is: second area is projected on the position of corresponding first area on the second image; Or second area is added on the position of corresponding first area on the second image.
According to a further aspect in the invention, provide a kind of characteristic point correction partitions splicing apparatus.
Characteristic point correction partitions of the present invention splicing apparatus comprises: acquisition module, for gathering the first view data, according to the characteristic point group of the first area on coordinate selection first image of existing electric power facility; Correction module, proofreaies and correct for each unique point to characteristic point group, forms second area; And concatenation module, for second area being spliced to the position of corresponding first area on the second image.
Further, correction module comprises: computing module, for calculating the Minimum Mean Square Error of each unique point; And syndrome module, for each unique point being proofreaied and correct according to Minimum Mean Square Error.
Further, computing module comprises: rigid transformation module, obtain the Minimum Mean Square Error of each unique point by rigid transformation for the basic point coordinate to each unique point.
Further, computing module comprises: rotation transform module, obtain the Minimum Mean Square Error of each unique point by rotational transform for the basic point coordinate to each unique point.
Further, concatenation module comprises: projection module, for second area being projected to the position of corresponding first area on the second image; Or laminating module, for the position of corresponding first area on the second image that second area is added to.
According to technical scheme of the present invention, effectively solve the digital convergence of accumulative various paper drawings over the years, the papery geologic map that comprises different time, different dry humidity, different scale closes, and can superpose, merge with geology, geography, electrical network vector data accurately.Meet the request for utilizations of each stage to geological information such as planning, design, construction, operation, maintenance.
Brief description of the drawings
Figure of description is used to provide a further understanding of the present invention, forms the application's a part, and schematic description and description of the present invention is used for explaining the present invention, does not form inappropriate limitation of the present invention.In the accompanying drawings:
Fig. 1 is according to the process flow diagram of the key step of the characteristic point correction partitions joining method of the embodiment of the present invention;
Fig. 2 is the concrete splicing schematic diagram according to the characteristic point correction partitions joining method of the embodiment of the present invention; And
Fig. 3 is according to the structural representation of the characteristic point correction partitions splicing apparatus of the embodiment of the present invention.
Embodiment
It should be noted that, in the situation that not conflicting, the feature in embodiment and embodiment in the application can combine mutually.Describe below with reference to the accompanying drawings and in conjunction with the embodiments the present invention in detail.
Fig. 1 is according to the process flow diagram of the key step of the characteristic point correction partitions joining method of the embodiment of the present invention.
Shown in Figure 1, characteristic point correction partitions of the present invention joining method comprises: S2: gather the first view data, and according to the characteristic point group of the first area on coordinate selection first image of existing electric power facility; S4: each unique point in characteristic point group is proofreaied and correct, form second area; And S6: second area is spliced on the position of corresponding first area on the second image.
By the technical scheme of the present embodiment, characteristic point correction partitions splicing can realize drawing error and drop to meter level by hundred meter levels.
Particularly, adopting characteristic point correction partitions joining method is subarea-scanning by figure such as " Beijing geologic maps " of large format, be characteristic point group by selecting the existing coordinate of relevant multiple power engineering, corresponding scintigram Sub-region and hierarchical is projected in electric power geography information space, Beijing, improve the stack precision of raster graphic and geographical space, electrical network is distributed and is associated with corresponding landforms, geological space data, be designed to object with power engineering, by planning, just establish, arrange, move, transform as timing node, dynamic management geology and power engineering information.
Preferably, step S4 comprises: S41: the Minimum Mean Square Error that calculates each unique point; And S43: each unique point is proofreaied and correct according to Minimum Mean Square Error.
Alternatively, in step S41, the computing method of Minimum Mean Square Error are: the basic point coordinate of each unique point is obtained to the Minimum Mean Square Error of each unique point by rigid transformation, the formula that Minimum Mean Square Error adopts is:
Figure BDA0000076743620000031
wherein, x i, y ifor the basic point coordinate of unique point;
Figure BDA0000076743620000032
for least mean-square error; F (x i) be rigid transformation, in piece image, the distance between any two points remains unchanged before and after conversion.
Above-mentioned algorithm is cost function, uses conventional methods using iterative method, and the time overhead of this method is larger, and needs more registration point.
Preferably, in step S41, the computing method of Minimum Mean Square Error are: the basic point coordinate of each unique point is obtained to the Minimum Mean Square Error of each unique point by rotational transform, the formula that Minimum Mean Square Error adopts is
Figure BDA0000076743620000033
wherein, x i, y ifor the basic point coordinate of unique point;
Figure BDA0000076743620000034
for least mean-square error; R (x i) be rotational transform, allow every bit P determine angle around one of a point of fixity rotation, become another P ', the conversion so producing is called the rotational transform in plane; T is translation transformation, between image along x, the translational movement in y direction.
Above-mentioned algorithm is the least-squares algorithm based on svd (SVD).This algorithm only needs less registration point just can calculate fast rotational transform matrix and calculate translation vector simultaneously, then this formula is minimized, and obtains least mean-square error.
Fig. 2 is the concrete splicing schematic diagram according to the characteristic point correction partitions joining method of the embodiment of the present invention.
Shown in Figure 2, by above-mentioned two kinds of algorithms, Minimum Mean Square Error calculating is carried out in the region of Fig. 2, and according to the Minimum Mean Square Error calculating, the region of Fig. 2 is adjusted.For example figure is upper local or global figure layer, and for bending, the broken figure layer that waits, by hinting obliquely at or projection, finally by crossing the adjustment of characteristic point correction partitions of the present invention joining method, obtains meeting the image of expection imagination.
Preferably, the method in step S6, second area being spliced on the second image is: second area is projected on the position of corresponding first area on the second image; Or second area is added on the position of corresponding first area on the second image.
By the technical scheme of the present embodiment, make large scale drawing scanning, can meet permissible accuracy requirement in engineering early stage, because the precision of paper figure is subject to the Accuracy of printing equipment, be subject to the impact of archiving mode, be subject to scanning the impact of on-the-spot humidity, temperature, scanning device precision, the figure that size is larger, accuracy error is larger, in order to improve precision, meet the basic demand of previous work, after the characteristic point group of presumptive area is proofreaied and correct, be added on corresponding position.
According to a further aspect in the invention, provide a kind of characteristic point correction partitions splicing apparatus.
Shown in Figure 3, characteristic point correction partitions of the present invention splicing apparatus comprises: acquisition module 20, and for gathering the first view data, and according to the characteristic point group of the first area on coordinate selection first image of existing electric power facility; Correction module 40, proofreaies and correct for each unique point to characteristic point group, forms second area; And concatenation module 60, for second area being spliced to the position of corresponding first area on the second image.
Preferably, correction module 40 comprises: computing module, for calculating the Minimum Mean Square Error of each unique point; And syndrome module, for each unique point being proofreaied and correct according to Minimum Mean Square Error.
Preferably, computing module comprises: rigid transformation module, obtain the Minimum Mean Square Error of each unique point by rigid transformation for the basic point coordinate to each unique point.
Preferably, computing module comprises: rotation transform module, obtain the Minimum Mean Square Error of each unique point by rotational transform for the basic point coordinate to each unique point.
Preferably, concatenation module 60 comprises: projection module, for second area being projected to the position of corresponding first area on the second image; Or laminating module, for the position of corresponding first area on the second image that second area is added to.
Obviously, those skilled in the art should be understood that, according to technical scheme of the present invention, effectively solve the digital convergence of accumulative various paper drawings over the years, the papery geologic map that comprises different time, different dry humidity, different scale closes, and can superpose, merge with geology, geography, electrical network vector data accurately.Meet the request for utilizations of each stage to geological information such as planning, design, construction, operation, maintenance.
The foregoing is only the preferred embodiments of the present invention, be not limited to the present invention, for a person skilled in the art, the present invention can have various modifications and variations.Within the spirit and principles in the present invention all, any amendment of doing, be equal to replacement, improvement etc., within all should being included in protection scope of the present invention.

Claims (4)

1. a characteristic point correction partitions joining method, is characterized in that, comprising:
S2: gather the first view data, according to the characteristic point group of the first area on the first image described in the coordinate selection of existing electric power facility;
S4: each unique point in described characteristic point group is proofreaied and correct, form second area; And
S6: described second area is spliced on the position of corresponding described first area on the second image,
Described step S4 comprises:
S41: the Minimum Mean Square Error that calculates described each unique point; And
S43: according to described Minimum Mean Square Error, described each unique point is proofreaied and correct,
In described step S41, the computing method of Minimum Mean Square Error are:
The Minimum Mean Square Error that the basic point coordinate of described each unique point is obtained to described each unique point by rigid transformation, the formula that calculates described Minimum Mean Square Error employing is: min F E = min F 1 N Σ i = 0 N | | y i - F ( x i ) | | 2 ;
Wherein, x i, y ifor the basic point coordinate of unique point;
Figure FDA0000456233160000014
e is least mean-square error; F (x i) be rigid transformation, in piece image, the distance between any two points remains unchanged before and after conversion,
Or in described step S41, the computing method of Minimum Mean Square Error are:
The Minimum Mean Square Error that the basic point coordinate of described each unique point is obtained to described each unique point by rotational transform, the formula that calculates described Minimum Mean Square Error employing is min F E = 1 N Σ i = 0 N | | y i - R ( x i ) - T | | 2 ;
Wherein, x i, y ifor the basic point coordinate of unique point;
Figure FDA0000456233160000013
e is least mean-square error; R (x i) be rotational transform, allow every bit P determine angle around one of a point of fixity rotation, become another P ', the conversion so producing is called the rotational transform in plane; T is translation transformation, between image along x, the translational movement in y direction.
2. method according to claim 1, is characterized in that, the method in described step S6, described second area being spliced on the second image is:
Described second area is projected on the position of corresponding described first area on described the second image; Or
Described second area is added on the position of corresponding described first area on described the second image.
3. a characteristic point correction partitions splicing apparatus, is characterized in that, comprising:
Acquisition module, for gathering the first view data, according to the characteristic point group of the first area on the first image described in the coordinate selection of existing electric power facility;
Correction module, proofreaies and correct for each unique point to described characteristic point group, forms second area; And
Concatenation module, for described second area being spliced to the position of corresponding described first area on the second image,
Described correction module comprises:
Computing module, for calculating the Minimum Mean Square Error of described each unique point; And
Syndrome module, for according to described Minimum Mean Square Error, described each unique point being proofreaied and correct,
Described computing module comprises: rigid transformation module, the Minimum Mean Square Error of each unique point described in obtaining by rigid transformation for the basic point coordinate to described each unique point.
Or described computing module comprises: rotation transform module, the Minimum Mean Square Error of each unique point described in obtaining by rotational transform for the basic point coordinate to described each unique point.
4. device according to claim 3, is characterized in that, described concatenation module comprises:
Projection module, for projecting to described second area the position of corresponding described first area on described the second image; Or
Laminating module, for the position of corresponding described first area on described the second image that described second area is added to.
CN201110201465.XA 2011-05-13 2011-07-19 Method and device for splicing characteristic point correction partitions Active CN102324096B (en)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN201110201465.XA CN102324096B (en) 2011-07-19 2011-07-19 Method and device for splicing characteristic point correction partitions
US14/117,405 US9256981B2 (en) 2011-05-13 2012-01-20 Method and device for processing geological information
EP12785279.6A EP2711893B1 (en) 2011-05-13 2012-01-20 Method and device for processing geological information
PCT/CN2012/070646 WO2012155540A1 (en) 2011-05-13 2012-01-20 Method and device for processing geological information

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110201465.XA CN102324096B (en) 2011-07-19 2011-07-19 Method and device for splicing characteristic point correction partitions

Publications (2)

Publication Number Publication Date
CN102324096A CN102324096A (en) 2012-01-18
CN102324096B true CN102324096B (en) 2014-07-02

Family

ID=45451835

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110201465.XA Active CN102324096B (en) 2011-05-13 2011-07-19 Method and device for splicing characteristic point correction partitions

Country Status (1)

Country Link
CN (1) CN102324096B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104599257B (en) * 2013-10-30 2018-11-13 重庆伟踱医疗设备股份有限公司 A kind of image co-registration display methods
CN107036628B (en) * 2017-04-10 2020-09-01 中国船舶重工集团公司第七0七研究所 Method for correcting paper chart
CN110472538B (en) * 2019-07-31 2023-06-06 河南冠图信息科技有限公司 Image recognition method and storage medium of electronic drawing

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002298141A (en) * 2001-03-29 2002-10-11 Nec Corp Pattern collating device, pattern collating method thereof, and pattern collating program
JP4946730B2 (en) * 2007-08-27 2012-06-06 ソニー株式会社 Face image processing apparatus, face image processing method, and computer program
JP5406705B2 (en) * 2009-12-28 2014-02-05 キヤノン株式会社 Data correction apparatus and method
CN101833644B (en) * 2010-06-01 2012-06-06 福建新大陆电脑股份有限公司 Correction graph searching method based on dynamic template

Also Published As

Publication number Publication date
CN102324096A (en) 2012-01-18

Similar Documents

Publication Publication Date Title
US9256981B2 (en) Method and device for processing geological information
CN102645209B (en) Joint positioning method for spatial points by means of onboard LiDAR point cloud and high resolution images
CN103196426A (en) Building surveying method utilizing total station and three-dimensional laser scanner
CN103823981A (en) DEM (Digital Elevation Model)-assisted satellite image block adjustment method
CN101882320A (en) Method for digital restoration of urban spatial pattern
CN102779165A (en) Building method of grid map picture base
US20220414982A1 (en) Methods for Correcting and Encrypting Space Coordinates of Three-Dimensional Model
CN110298778B (en) Traffic cell demographic method based on raster data and area ratio correction
CN115796712B (en) Regional land ecosystem carbon reserve estimation method and device and electronic equipment
KR100686287B1 (en) Distorting Modeling method for Transforming the Presize Position of Partial/Positional information
CN102324096B (en) Method and device for splicing characteristic point correction partitions
Tuno et al. Improving the positional accuracy of digital cadastral maps through optimal geometric transformation
CN107067454B (en) Radar force diagram three-dimensional display method based on hypergraph component development platform
CN111696156A (en) Control point-free remote sensing image coordinate conversion method
CN113536412A (en) BIM and GIS-based high-precision fusion method for three-dimensional live-action model
KR101385111B1 (en) Calculation method of coordinate transformation coefficient of linear form for notify of topographic drawing in framework act on restrictions on land utilization [farlu]
CN105608714A (en) Framing image generation method based on CAD and GIS cooperation and system thereof
CN103017749B (en) Method, apparatus and navigator for converting narrow and long water system surface element into line element
CN107504959A (en) Utilize the method for oblique aerial radiographic measurement house wall base profile
CN103235290B (en) Geographic space lattice-based radar detection data processing method
CN109829254B (en) Method for drawing positioning line of diaphragm plate after top and bottom plates of curved steel box bridge are unfolded
Tong et al. Geometric integration of aerial and QuickBird imagery for high accuracy geopositioning and mapping application: A case study in Shanghai
CN109947877A (en) A kind of method and system improving GIS mobile terminal Orientation on map precision
CN105468693B (en) A kind of efficient topological relationship calculation method
CN116643290B (en) Metering method and system for double-platform motion compensation of irregular contour

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
ASS Succession or assignment of patent right

Owner name: STATE ELECTRIC NET CROP.

Free format text: FORMER OWNER: BEIJING ELECTRIC POWER CORP.

Effective date: 20130105

Owner name: BEIJING ELECTRIC POWER CORP. BEIJING POWER ECONOMI

Free format text: FORMER OWNER: BEIJING POWER ECONOMIC RESEARCH INSTITUTE

Effective date: 20130105

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 100031 XICHENG, BEIJING TO: 100031 DONGCHENG, BEIJING

TA01 Transfer of patent application right

Effective date of registration: 20130105

Address after: 100031 West Chang'an Avenue, Beijing, No. 86

Applicant after: State Grid Corporation of China

Applicant after: Beijing Electric Power Corp.

Applicant after: State Power Economic Research Institute

Address before: 100031 Beijing Qianmen West Street, Xicheng District, No. 41

Applicant before: Beijing Electric Power Corp.

Applicant before: State Power Economic Research Institute

C14 Grant of patent or utility model
GR01 Patent grant