CN109191579A - Method, apparatus, computer equipment and the storage medium of coordinate conversion - Google Patents

Method, apparatus, computer equipment and the storage medium of coordinate conversion Download PDF

Info

Publication number
CN109191579A
CN109191579A CN201810706859.2A CN201810706859A CN109191579A CN 109191579 A CN109191579 A CN 109191579A CN 201810706859 A CN201810706859 A CN 201810706859A CN 109191579 A CN109191579 A CN 109191579A
Authority
CN
China
Prior art keywords
coordinate
grid
reduction
coordinate system
point
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.)
Granted
Application number
CN201810706859.2A
Other languages
Chinese (zh)
Other versions
CN109191579B (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.)
Guangdong Provincial Institute Of Geological Surveying And Mapping
Original Assignee
Guangdong Provincial Institute Of Geological 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 Guangdong Provincial Institute Of Geological Surveying And Mapping filed Critical Guangdong Provincial Institute Of Geological Surveying And Mapping
Priority to CN201810706859.2A priority Critical patent/CN109191579B/en
Publication of CN109191579A publication Critical patent/CN109191579A/en
Application granted granted Critical
Publication of CN109191579B publication Critical patent/CN109191579B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite element generation, e.g. wire-frame surface description, tesselation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Processing Or Creating Images (AREA)
  • Image Processing (AREA)

Abstract

The present invention relates to method, apparatus, computer equipment and the storage mediums of coordinate conversion, belong to technical field of geographic information.The described method includes: determining the characteristic parameter that the first coordinate system is converted to the second coordinate system according to known point;The conversion residual error of the known point and the coordinate correction amount of grid points are determined according to the characteristic parameter;The residual error reduction of the grid points is determined according to the conversion residual error of the known point;The grid reduction that the grid points are determined according to the coordinate correction amount and the residual error reduction establishes transformation model according to the grid reduction;The transformation model is used to determine that the point to be converted in the first coordinate system to be transformed into the coordinate in second coordinate system.Above-mentioned technical proposal, solve the problems, such as coordinate transform process ignore point spatial coherence make convert coordinate accuracy it is not high.Grid reduction is solved on the basis of considering known point reduction, accounts for spatial coherence, realizes high-precision coordinate conversion.

Description

Method, apparatus, computer equipment and the storage medium of coordinate conversion
Technical field
The present invention relates to technical field of geographic information, more particularly to the method, apparatus of coordinate conversion, computer equipment and Storage medium.
Background technique
Presently, there are a variety of coordinate systems, due to differences such as the ellipsoids, frame, epoch that are used between different coordinate systems, This just needs to carry out the conversion of coordinate, such as 1980 coordinate system of Xi'an is turned to CGCS2000 (2000 national earth coordinates) It changes.Traditional coordinate conversion is often using modes such as four parameters, seven parameter models.But four parameters are suitable only for difference of longitude, latitude Area coordinate conversion of the degree difference less than 2 degree, seven parameters are suitable only for the area coordinate conversion of provincial, and municipal level.In recent years, also there is research Personnel propose the conversion that coordinate is realized using grid correction model.In realizing process of the present invention, inventor has found existing skill At least there are the following problems in art: conventional grid correction model corrects in grid points only according to the anti-weighting of distance to realize The calculating of amount has ignored the spatial coherence of point in this way, and obtained conversion coordinate accuracy is not high.
Summary of the invention
Based on this, the present invention provides method, apparatus, computer equipment and the storage mediums of coordinate conversion, can obtain standard True conversion coordinate.
The content of the embodiment of the present invention is as follows:
A kind of method of coordinate conversion, comprising the following steps: determined according to known point and the first coordinate system is converted to second The characteristic parameter of coordinate system;The known point has the coordinate of first coordinate system and the coordinate of the second coordinate system;According to institute It states characteristic parameter and determines the conversion residual error of the known point and the coordinate correction amount of the grid points in the grid that pre-establishes;Root The residual error reduction of the grid points is determined according to the conversion residual error of the known point;According to the coordinate correction amount and the residual error Reduction determines the grid reduction of the grid points, establishes transformation model according to the grid reduction;The transformation model Point to be converted for determining in the first coordinate system is transformed into the coordinate in second coordinate system.
The characteristic parameter includes four parameters in one of the embodiments,;Described determined according to known point is sat first The step of mark system is converted to the characteristic parameter of the second coordinate system, comprising: determine first coordinate of the known point in the first coordinate system And the second coordinate in the second coordinate system, calculate four parameters that first coordinate is converted to the second coordinate.
The coordinate of the grid points in the grid pre-established is determined according to the characteristic parameter in one of the embodiments, The step of reduction, comprising: determine first grid coordinate of the grid points in the first coordinate system, it is true according to four parameter Fixed the second grid coordinate of the grid points in the second coordinate system;Institute is determined according to the first grid coordinate and the second grid coordinate State the coordinate correction amount of grid points.
The conversion residual error according to the known point determines that the residual error of the grid points changes in one of the embodiments, The step of positive quantity, comprising: determine the known point at a distance from the grid points;According to the conversion residual sum institute of the known point Distance is stated, the residual error reduction of the grid points is determined by distance weighting method.
It is described in one of the embodiments, to determine that grid corrects according to the coordinate correction amount and the residual error reduction The step of amount, comprising: the grid reduction is obtained according to the sum of the coordinate correction amount and the residual error reduction.
It in one of the embodiments, include the grid reduction of multiple grid points in the transformation model;Described The grid reduction that the grid points are determined according to the coordinate correction amount and the residual error reduction, according to the grid reduction After the step of establishing transformation model, further includes: determine target lattice corresponding with the point to be converted from the transformation model The grid reduction of site;According to the grid reduction of the target mesh point, determined by bilinear interpolation method described wait turn Change reduction a little;It is determined according to the reduction of the described to be converted o'clock coordinate in the first coordinate system and the point to be converted Coordinate of the point to be converted in second coordinate system.
It is described in one of the embodiments, that the spy that the first coordinate system is converted to the second coordinate system is determined according to known point Before the step of levying parameter, further includes: the region to be analyzed within the scope of the grid points pre-set radius is determined, from described to be analyzed The known point with the coordinate of first coordinate system and the second coordinate system is determined in region.
Correspondingly, the embodiment of the present invention provides a kind of device of coordinate conversion, comprising: parameter determination module is used for basis Known point determines the characteristic parameter that the first coordinate system is converted to the second coordinate system;The known point has first coordinate system Coordinate and the second coordinate system coordinate;Reduction determining module, for determining the known point according to the characteristic parameter The coordinate correction amount of conversion residual error and the grid points in the grid pre-established;It is determined according to the conversion residual error of the known point The residual error reduction of the grid points;And model building module, for being corrected according to the coordinate correction amount and the residual error Amount determines that the grid reduction of the grid points establishes transformation model;The transformation model be used for determine in the first coordinate system to Transfer point is transformed into the coordinate in second coordinate system.
The method and device of above-mentioned coordinate conversion is determined according to known point in the coordinate of the first coordinate system and the second coordinate system Characteristic parameter;The conversion residual error of known point is determined according to this feature parameter, and the residual error of grid points is determined according to the conversion residual error Reduction;The coordinate correction amount of grid points is determined according to required characteristic parameter;It is true according to residual error reduction and coordinate correction amount It fixes and net reduction and establishes transformation model, determined by transformation model and the point to be converted in the first coordinate system is transformed into Coordinate after two coordinate systems.Grid reduction is solved on the basis of considering known point reduction, accounts for spatial coherence, Realize high-precision coordinate conversion.
A kind of computer equipment can be run on a memory and on a processor including memory, processor and storage Computer program, the processor performs the steps of when executing the computer program to be determined according to known point by the first seat Mark system is converted to the characteristic parameter of the second coordinate system;The known point has the coordinate and the second coordinate system of first coordinate system Coordinate;According to the grid points in the grid that the characteristic parameter determines the conversion residual error of the known point and pre-establishes Coordinate correction amount;The residual error reduction of the grid points is determined according to the conversion residual error of the known point;Changed according to the coordinate Positive quantity and the residual error reduction determine the grid reduction of the grid points, establish modulus of conversion according to the grid reduction Type;The transformation model is used to determine that the point to be converted in the first coordinate system to be transformed into the coordinate in second coordinate system.
Above-mentioned computer equipment solves grid reduction on the basis of considering known point reduction, accounts for space Correlation realizes high-precision coordinate conversion.
A kind of computer readable storage medium, is stored thereon with computer program, and the computer program is held by processor The characteristic parameter for determining according to known point and the first coordinate system being converted to the second coordinate system is performed the steps of when row;It is described Know the coordinate of a little coordinate with first coordinate system and the second coordinate system;The known point is determined according to the characteristic parameter Conversion residual error and the grid points in the grid that pre-establishes coordinate correction amount;It is true according to the conversion residual error of the known point The residual error reduction of the fixed grid points;The lattice of the grid points are determined according to the coordinate correction amount and the residual error reduction Net reduction establishes transformation model according to the grid reduction;The transformation model be used for determine in the first coordinate system to Transfer point is transformed into the coordinate in second coordinate system.
Above-mentioned computer readable storage medium solves grid reduction on the basis of considering known point reduction, cares for And spatial coherence, realize high-precision coordinate conversion.
Detailed description of the invention
Fig. 1 is the applied environment figure of the method for coordinate conversion in one embodiment;
Fig. 2 is the flow diagram of the method for coordinate conversion in one embodiment;
Fig. 3 is the schematic diagram of the grid in region to be converted in one embodiment;
Fig. 4 is the three dimensional structure diagram of transformation model in one embodiment;
Fig. 5 projects to the schematic diagram of grid for point to be converted in one embodiment;
Fig. 6 is the schematic diagram that known point is chosen in one embodiment;
Fig. 7 is the flow diagram of the method for coordinate conversion in another embodiment;
Fig. 8 is the structural block diagram of the device of coordinate conversion in one embodiment.
Specific embodiment
In order to make the objectives, technical solutions, and advantages of the present invention clearer, with reference to the accompanying drawings and embodiments, right The present invention is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, and It is not used in the restriction present invention.
The method of coordinate conversion provided by the present application can be applied in computer equipment as shown in Figure 1.The computer Equipment can be server, and internal structure chart can be as shown in Figure 1.The computer equipment includes being connected by system bus Processor, memory, network interface and database.Wherein, the processor of the computer equipment calculates and controls energy for providing Power.The memory of the computer equipment includes non-volatile memory medium, built-in storage.The non-volatile memory medium is stored with Operating system, computer program and database.The built-in storage is operating system and computer in non-volatile memory medium The operation of program provides environment.The database of the computer equipment is for storing coordinate, coordinate correction amount, residual error reduction, lattice The data such as net reduction.The network interface of the computer equipment is used to communicate with external terminal by network connection.The calculating A kind of machine program realizes coordinate conversion method when being executed by processor.
It will be understood by those skilled in the art that structure shown in Fig. 1, only part relevant to application scheme is tied The block diagram of structure does not constitute the restriction for the computer equipment being applied thereon to application scheme, specific computer equipment It may include perhaps combining certain components or with different component layouts than more or fewer components as shown in the figure.
The embodiment of the present invention provides method, apparatus, computer equipment and the storage medium of a kind of coordinate conversion.Individually below It is described in detail.
In one embodiment, it as shown in Fig. 2, providing a kind of method of coordinate conversion, is able to achieve and treats transition region The coordinate of interior point to be converted is converted, method includes the following steps:
S201, the characteristic parameter that the first coordinate system is converted to the second coordinate system is determined according to known point;The known point The coordinate of coordinate and the second coordinate system with first coordinate system.
In this step, it is known that the coordinate with the first coordinate system and the second coordinate system is put, according to conventional calculation method The characteristic parameter of transformational relation between the first coordinate system of computational representation and the second coordinate system.
Wherein, the first coordinate system and the second coordinate system can refer to 1980 coordinate system of Xi'an, CGCS2000 coordinate system, Beijing 54 The coordinate systems such as coordinate system, WGS-84 coordinate system, the embodiment of the present invention do not limit the classification of the first coordinate system and the second coordinate system System can guarantee that the two is different and can carry out coordinate conversion.
Characteristic parameter refers to the parameter of the transformational relation between the first coordinate system of characterization and the second coordinate system, can be four The parameters such as parameter, seven parameters.
In addition, the known point can refer to all coordinates with the first coordinate system and the second coordinate system in region to be converted Point, be also possible to a part (such as: the point within the scope of pre-set radius) therein.Also, it is known that point can be area to be converted Random point (choosing as long as the coordinate that it has the first coordinate system and the second coordinate system) in domain, is also possible to according to one Certain points that fixed rule chooses from multiple known points in region to be converted.
Lattice in S202, the conversion residual error that the known point is determined according to the characteristic parameter and the grid pre-established The coordinate correction amount of site;The residual error reduction of the grid points is determined according to the conversion residual error of the known point.
This step is accomplished that the coordinate correction amount of grid points and the determination of residual error reduction.Wherein, coordinate correction amount It is the characteristic parameter that is obtained according to previous step to determine;And residual error reduction is the root after determining the conversion residual error of known point It is determined according to the conversion residual error and known point of known point and the relationship of grid points.
Wherein, the grid pre-established refers to the latice drawn on the first coordinate system according to certain mesh spacing Picture.Wherein mesh spacing refers to: the longitude and latitude difference of the grid adjacent in length and breadth of established grid.Mesh spacing is mainly known to The equispaced of point determines that general mesh spacing is less than the equispaced of known point;Can also according to transform graphical or to The size of transition region determines.Wherein transform graphical refers to geography that need to convert, according to certain standard progress framing Information graphic.
The grid established in a certain region to be converted can with as shown in figure 3, the grid in Fig. 3 by grid adjacent in length and breadth It staggeredly constitutes, each grid points has its coordinate on the first coordinate system.
S203, the grid reduction that the grid points are determined according to the coordinate correction amount and the residual error reduction, root Transformation model is established according to the grid reduction;The transformation model is for determining that the point to be converted in the first coordinate system is transformed into Coordinate in second coordinate system.
Wherein, point to be converted refers to the coordinate with the first coordinate system but does not know its coordinate in the second coordinate system Point.Grid reduction refers to the reduction of some grid points, and grid points are distributed in each place in region to be converted, each Point to be converted can determine the reduction of oneself according to neighbouring grid points.
In this step, the grid reduction of grid points is determined according to coordinate correction amount and residual error reduction, and according to lattice Net reduction establishes transformation model, according to the transformation model and the to be converted o'clock coordinate in the first coordinate system just can determine that The coordinate of transfer point in the second coordinate system.
Transformation model can be determined according to the grid reduction of one or more grid points in grid, such as can be It is determined according to the grid reduction of all grid points in region to be converted.
Further, the method for establishing transformation model, which can be, deposits identified grid reduction according to certain format It puts in one file, this document can serve as transformation model.Required mesh can be therefrom found out by reading this document The grid reduction of style site, and then point to be converted is determined in the second coordinate system according to the grid reduction of target mesh point Coordinate.Wherein, target mesh point can be multiple.
The present embodiment, the acquisition process of grid reduction using multiple known points of neighboring area reduction, i.e., these Known point has been involved in the solution of grid reduction, accounts for spatial coherence, and the transformation model obtained in this way is continuously smooth , so that the coordinate conversion in big region has continuity, to realize the seamless big area coordinate conversion of high-precision.
In one embodiment, the characteristic parameter includes four parameters;It is described to be determined according to known point by the first coordinate system The step of being converted to the characteristic parameter of the second coordinate system, comprising: determine first coordinate of the known point in the first coordinate system and The second coordinate in the second coordinate system calculates four parameters that first coordinate is converted to the second coordinate.
If known point be it is multiple, can therefrom select a known point, according to the first coordinate of selected known point and Second coordinate determines four parameters.Four parameters can also be determined according to the coordinate of all or part known point, if obtaining There is four parameter of multiple groups, then these four parameters can be integrated by certain operation.
It is four parameters embodiment defines characteristic parameter, the determination process of four parameters is using conventional method realization It can.Four parameters consider the deviation of the coordinate conversion in small range region, can effectively improve the precision of identified grid reduction. The determination process of four parameters is simply direct simultaneously, and the calculating for follow-up correction amount provides Data safeguard.
In one embodiment, the coordinate correction of the grid points in the grid pre-established is determined according to the characteristic parameter The step of amount, comprising: determine first grid coordinate of the grid points in the first coordinate system, institute is determined according to four parameter State the second grid coordinate of grid points in the second coordinate system;The lattice are determined according to the first grid coordinate and the second grid coordinate The coordinate correction amount of site.
Because grid is established on the basis of the first coordinate system, each grid points are in the first coordinate system First grid coordinate can be directly obtained.
Due to transformational relation of four parameter characterizations between the first coordinate system and the second coordinate system, according to four parameters It just can determine that out the second grid coordinate of grid in the second coordinate system, the first grid coordinate and the second grid coordinate compared Compared with the coordinate correction amount that (such as: seeking difference) just can determine that out grid points.
The present embodiment proposes to solve grid points using four parametric methods of part on the basis of generally conventional grid correction method Coordinate correction amount determines the second grid coordinate of grid according to four parameters, and according to the first grid coordinate and the second grid coordinate To determine the coordinate correction amount of grid.The coordinate correction amount had not only characterized grid in the coordinate of the first coordinate system, but also characterized lattice The coordinate of net in the second coordinate system, therefore can be used to determine the coordinate of point to be converted in the second coordinate system.
In one embodiment, the conversion residual error according to the known point determines the residual error reduction of the grid points The step of, comprising: determine the known point at a distance from the grid points;According to the conversion residual sum of the known point away from From determining the residual error reduction of the grid points by distance weighting method.
Wherein, determine that the process of the residual error reduction of grid points can be with by distance weighting method are as follows: pass through inverse distance-weighting Method determines weight (distance bigger weight smaller) of the conversion residual error for residual error reduction of each known point, combines conversion later The residual error reduction of residual error and its corresponding weight computing grid points.
The present embodiment is weighted place during solving residual error reduction, according to the conversion residual error of each known point Reason, so that the calculation result of grid points is more in line with true value, improves conversion accuracy.
In one embodiment, described that grid reduction is determined according to the coordinate correction amount and the residual error reduction Step, comprising: the grid reduction is obtained according to the sum of the coordinate correction amount and the residual error reduction.
The present embodiment determines that grid is corrected according to the residual error reduction and coordinate correction amount that acquire in previous embodiment Amount.Not only allow for the grid deviation of conversion process, it is also contemplated that conversion residual error, so that the considerations of deviation more fully, it can More accurate grid reduction is obtained, and then realizes more accurate coordinate conversion.
It in one embodiment, include the grid reduction of multiple grid points in the transformation model;It is described according to institute It states coordinate correction amount and the residual error reduction determines the grid reduction of the grid points, established according to the grid reduction After the step of transformation model, further includes: determine target mesh point corresponding with the point to be converted from the transformation model Grid reduction;According to the grid reduction of the target mesh point, the point to be converted is determined by bilinear interpolation method Reduction;According to described to be converted o'clock coordinate and the point to be converted in the first coordinate system reduction determine described in Coordinate of the point to be converted in second coordinate system.
It include that the grid reduction of multiple grid points can refer in transformation model;Transformation model includes in region to be converted The grid reduction of all grid points.The tomograph drawn according to transformation model as shown in figure 4, in Fig. 4 different coordinates Grid reduction may be different.Grid reduction in Fig. 4 is bigger to illustrate that deviation is bigger, and Fig. 4 has fully demonstrated grid correction The necessity of amount, because if not considering grid reduction, then biggish deviation will be had by converting later coordinate.
Wherein, the reduction of the coordinate according to described to be converted o'clock in the first coordinate system and the point to be converted determines The realization process of coordinate of the point to be converted in second coordinate system can be with are as follows: by be converted o'clock in the first coordinate system Coordinate coordinate of the point to be converted in second coordinate system can be obtained plus the reduction of point to be converted.
In the present embodiment, to carry out coordinate conversion to a certain point to be converted, it is first determined the lattice where the point to be converted Net determines the reduction of point to be converted (here, point to be converted claims by bilinear interpolation method according to the grid reduction of grid points For interpolated point).The reduction of point to be converted can be accurately determined out, and then determines that point to be converted converts later coordinate.
In one embodiment, determine that the process of the reduction of interpolated point can be with by bilinear interpolation method are as follows:
(five-pointed star indicates that interpolated point, soft dot indicate grid points in Fig. 5) as shown in Figure 5, the corresponding lattice of interpolated point There are four grid points in net.Interpolation point is projected respectively to the x-axis of place grid, y-axis direction, calculates separately projector distance, root Corresponding weight is determined according to the projector distance of x-axis, y-axis direction, is changed according to the grid of these weights and corresponding four grid points Positive quantity determines the reduction of interpolated point.
In one embodiment, the corresponding grid points of interpolated point can be with more than four, for example, can be 8,16 etc. Deng.
In one embodiment, described to determine that the feature that the first coordinate system is converted to the second coordinate system is joined according to known point Before several steps, further includes: the region to be analyzed within the scope of the grid points pre-set radius is determined, from the region to be analyzed Middle determination has the known point of the coordinate of first coordinate system and the second coordinate system.
Wherein, pre-set radius can be 1m, 10m, 10km etc., and the radius is changeable, and the embodiment of the present invention is to default half The size of diameter is with no restrictions.
The present embodiment determines region to be analyzed according to pre-set radius centered on grid points, and institute is determined from region to be analyzed There is known point, determines that the grid reduction of grid points, all known points have both participated in grid reduction according to these known points Solution, fully considered spatial coherence, accurate grid reduction can be obtained.
Wherein, the schematic diagram for known point being chosen in region to be analyzed is shown in Fig. 6, using the corresponding grid points in the center of circle as in The center of these circles and all known points are all chosen and (do not consider other grid points temporarily) by the heart.On the right side of Fig. 6 In figure, central point is grid points, other points are known point.
In one embodiment, as shown in fig. 7, providing a kind of method of coordinate conversion, comprising the following steps:
S701, region to be analyzed within the scope of grid points pre-set radius is determined, determination has from the region to be analyzed The known point of the coordinate of first coordinate system and the second coordinate system.
S702, first coordinate of the known point in the first coordinate system and the second coordinate in the second coordinate system are determined, Calculate four parameters that first coordinate is converted to the second coordinate.
S703, the conversion residual error that the known point is determined according to four parameter.
S704, determine the known point at a distance from the grid points;According to the conversion residual sum of the known point Distance determines the residual error reduction of the grid points by distance weighting method.
S705, it determines first grid coordinate of the grid points in the first coordinate system, institute is determined according to four parameter State the second grid coordinate of grid points in the second coordinate system;The lattice are determined according to the first grid coordinate and the second grid coordinate The coordinate correction amount of site.
S706, grid reduction is obtained according to the sum of the coordinate correction amount and the residual error reduction, according to the lattice Net reduction establishes transformation model.
S707, the grid reduction that target mesh point corresponding with point to be converted is determined from the transformation model;According to The grid reduction of the target mesh point determines the reduction of the point to be converted by bilinear interpolation method.
S708, it is determined according to the reduction of the described to be converted o'clock coordinate in the first coordinate system and the point to be converted Coordinate of the point to be converted in second coordinate system.
The present embodiment solves grid reduction on the basis of considering known point reduction, accounts for spatial coherence, Realize high-precision coordinate conversion.
The above method in order to better understand, the application example of the method for a coordinate conversion of the present invention set forth below.It should Application example is accomplished that the coordinate between 1980 coordinate of Xi'an and CGCS2000 coordinate is converted, wherein having in region to be converted Know point (each known point has 1980 coordinate of Xi'an and CGCS2000 coordinate).Detailed description are as follows:
One, transformation model is constructed:
1, in a known area map (region to be converted) known point 1980 coordinate system coordinate of Xi'an and CGCS2000 Coordinate system coordinate.
2, grid is made on this area map (grid coordinate is 1980 coordinate system coordinate of Xi'an).
3, by taking one of grid points A in grid as an example, grid points A reduction (D is determinedX802000, DY802000), wherein X indicate X-axis coordinate, Y indicate Y-axis coordinate, under be designated as 2000 expressions be CGCS2000 coordinate system X axis coordinate, subscript For 80 indicate be 1980 coordinate system of Xi'an X axis coordinate, 802000 indicate to be 1980 coordinate system of Xi'an sit to CGCS2000 The conversion of mark system,
The determination process of reduction is as follows:
S31, all known points are searched within the scope of the appropriate area of grid points A;
S32, coordinate (1980 coordinate system coordinate of Xi'an and CGCS2000 coordinate system seat according to each known point searched Mark) four parameters of conversion are calculated, and the conversion residual error of each known point is calculated based on four parameter;
S33, each known point is determined at a distance from grid points A, according to the conversion residual error of distance weighting method and known point Determine residual error the reduction ((d of grid points AX, dY);
S34, the CGCS2000 coordinate system coordinate that grid points A is calculated according to four parameters acquired, and then calculate grid points A's Coordinate correction amount (dX802000, dY802000);
S35, residual error reduction and the summation of coordinate correction amount obtain the grid reduction at grid points A:
(DX802000, DY802000)=(dX802000, dY802000)+(dX, dY);
Each grid reduction is placed on one according to certain format by S36, the grid reduction for determining other each grid points In a file, as transformation model.
Two, coordinate conversion is carried out according to transformation model:
4, determine which grid it is located within according to the 1980 coordinate system coordinate of Xi'an of point to be converted (B point).
5, determine point B's to be converted using the grid reduction of four grid points corresponding in bilinear interpolation method and grid Reduction.
6, the reduction of all points to be converted in 1980 coordinate system of Xi'an is determined, 1980 coordinate system coordinate of Xi'an is plus correction Amount just obtains coordinate of each point to be converted in CGCS2000, that is, completes 1980 coordinate system of Xi'an to CGCS2000 coordinate system Coordinate transform process.
The present embodiment solves grid reduction on the basis of considering known point reduction, accounts for spatial coherence, Realize high-precision coordinate conversion.
It should be noted that for the various method embodiments described above, describing for simplicity, it is all expressed as a series of Combination of actions, but those skilled in the art should understand that, the present invention is not limited by the sequence of acts described, because according to According to the present invention, certain steps can use other sequences or carry out simultaneously.
Based on thought identical with the method for coordinate conversion in above-described embodiment, the present invention also provides the dresses of coordinate conversion It sets, which can be used for executing the method for above-mentioned coordinate conversion.For ease of description, the structure of the Installation practice of coordinate conversion In schematic diagram, part related to the embodiment of the present invention illustrate only, it will be understood by those skilled in the art that schematic structure is simultaneously The not restriction of structure twin installation may include perhaps combining certain components or different than illustrating more or fewer components Component layout.
As described in Figure 8, the device of coordinate conversion includes that parameter determination module 801, reduction determining module 802 and model are built Formwork erection block 803, detailed description are as follows:
Parameter determination module 801, for determining the feature that the first coordinate system is converted to the second coordinate system according to known point Parameter;The known point has the coordinate of first coordinate system and the coordinate of the second coordinate system.
Reduction determining module 802, for determining the conversion residual error of the known point and pre- according to the characteristic parameter The coordinate correction amount for the grid points in grid first established;The residual of the grid points is determined according to the conversion residual error of the known point Poor reduction.
And model building module 803, for determining the lattice according to the coordinate correction amount and the residual error reduction The grid reduction of site establishes transformation model;The transformation model is for determining that the point to be converted in the first coordinate system is transformed into Coordinate in second coordinate system.
The present embodiment solves grid reduction on the basis of considering known point reduction, accounts for spatial coherence, Realize high-precision coordinate conversion.
In one embodiment, parameter determination module 801, for determining first coordinate of the known point in the first coordinate system And the second coordinate in the second coordinate system, calculate four parameters that first coordinate is converted to the second coordinate.
In one embodiment, reduction determining module 802, comprising: coordinate determines submodule, for determining the grid The first grid coordinate o'clock in the first coordinate system, determines of the grid points in the second coordinate system according to four parameter Two grid coordinates;Reduction determines submodule, for determining the grid points according to the first grid coordinate and the second grid coordinate Coordinate correction amount.
In one embodiment, reduction determining module 802, be also used to determine the known point and the grid points away from From;According to distance described in the conversion residual sum of the known point, determine that the residual error of the grid points corrects by distance weighting method Amount.
In one embodiment, model building module 803 are also used to be corrected according to the coordinate correction amount and the residual error The sum of amount obtains the grid reduction.
In one embodiment, further includes: target point reduction determining module, for from the transformation model determine with The grid reduction of the corresponding target mesh point of the point to be converted;To turning point reduction determining module, for according to the mesh The grid reduction of style site determines the reduction of the point to be converted by bilinear interpolation method;Coordinate transferring is used It is described to be converted in being determined according to the reduction of the described to be converted o'clock coordinate in the first coordinate system and the point to be converted Coordinate of the point in second coordinate system.
In one embodiment, further includes: known point determining module, for determining within the scope of the grid points pre-set radius Region to be analyzed, determined from the region to be analyzed known to the coordinate with first coordinate system and the second coordinate system Point.
It should be noted that the method that device and coordinate of the invention that coordinate of the invention is converted are converted corresponds, It is suitable for the device of coordinate conversion in the technical characteristic and its advantages of the embodiment elaboration of the method for above-mentioned coordinate conversion Embodiment in, for details, please refer to the description in the embodiment of the method for the present invention, and details are not described herein again, hereby give notice that.
In addition, the logical partitioning of each program module is only to lift in the embodiment for the device that the coordinate of above-mentioned example is converted Example explanation, can according to need in practical application, for example, for corresponding hardware configuration requirement or software realization convenience Consider, above-mentioned function distribution is completed by different program modules, i.e., the internal structure for the device converted the coordinate divides At different program modules, to complete all or part of the functions described above.
In one embodiment, a kind of computer equipment is provided, including memory, processor and storage are on a memory And the computer program that can be run on a processor, processor perform the steps of when executing computer program according to known point Determine the characteristic parameter that the first coordinate system is converted to the second coordinate system;The known point has the coordinate of first coordinate system With the coordinate of the second coordinate system;The grid for determining the conversion residual error of the known point according to the characteristic parameter and pre-establishing In grid points coordinate correction amount;The residual error reduction of the grid points is determined according to the conversion residual error of the known point;Root The grid reduction that the grid points are determined according to the coordinate correction amount and the residual error reduction, according to the grid reduction Establish transformation model;The transformation model is for determining that the point to be converted in the first coordinate system is transformed into second coordinate system Coordinate.
In one embodiment, also performing the steps of the characteristic parameter when processor executes computer program includes Four parameters;It is described that the step of first coordinate system is converted to the characteristic parameter of the second coordinate system is determined according to known point, comprising: really Determine first coordinate of the known point in the first coordinate system and the second coordinate in the second coordinate system, calculating is sat described first Mark is converted to four parameters of the second coordinate.
In one embodiment, it also performs the steps of when processor executes computer program according to the characteristic parameter The step of determining the coordinate correction amount of the grid points in the grid pre-established, comprising: determine the grid points in the first coordinate The first grid coordinate in system determines the second grid coordinate of the grid points in the second coordinate system according to four parameter; The coordinate correction amount of the grid points is determined according to the first grid coordinate and the second grid coordinate.
In one embodiment, it is also performed the steps of when processor executes computer program described according to described known The step of conversion residual error of point determines the residual error reduction of the grid points, comprising: determine the known point and the grid points Distance;According to distance described in the conversion residual sum of the known point, the residual error of the grid points is determined by distance weighting method Reduction.
In one embodiment, it is also performed the steps of when processor executes computer program described according to the coordinate The step of reduction and the residual error reduction determine grid reduction, comprising: according to the coordinate correction amount and the residual error The sum of reduction obtains the grid reduction.
In one embodiment, it also performs the steps of when processor executes computer program and is wrapped in the transformation model Grid reduction containing multiple grid points;It is described that the grid is determined according to the coordinate correction amount and the residual error reduction The grid reduction of point, after the step of establishing transformation model according to the grid reduction, further includes: from the transformation model The grid reduction of middle determination target mesh point corresponding with the point to be converted;It is corrected according to the grid of the target mesh point Amount, the reduction of the point to be converted is determined by bilinear interpolation method;According to described to be converted o'clock in the first coordinate system The reduction of coordinate and the point to be converted determines coordinate of the point to be converted in second coordinate system.
In one embodiment, it is also performed the steps of when processor executes computer program described true according to known point Before fixed the step of first coordinate system is converted to the characteristic parameter of the second coordinate system, further includes: determine that the grid points are default Region to be analyzed in radius determines there is first coordinate system and the second coordinate system from the region to be analyzed The known point of coordinate.
In one embodiment, a kind of computer readable storage medium is provided, computer program is stored thereon with, is calculated Machine program performs the steps of when being executed by processor is converted to the second coordinate system for the first coordinate system according to known point is determining Characteristic parameter;The known point has the coordinate of first coordinate system and the coordinate of the second coordinate system;Joined according to the feature The coordinate correction amount of the conversion residual error of the determining known point of number and the grid points in the grid pre-established;According to it is described The conversion residual error for knowing a little determines the residual error reduction of the grid points;It is true according to the coordinate correction amount and the residual error reduction The grid reduction of the fixed grid points, establishes transformation model according to the grid reduction;The transformation model is for determining Point to be converted in first coordinate system is transformed into the coordinate in second coordinate system.
In one embodiment, the characteristic parameter packet is also performed the steps of when computer program is executed by processor Include four parameters;It is described that the step of first coordinate system is converted to the characteristic parameter of the second coordinate system is determined according to known point, comprising: It determines first coordinate of the known point in the first coordinate system and the second coordinate in the second coordinate system, calculates described first Coordinate is converted to four parameters of the second coordinate.
In one embodiment, it also performs the steps of when computer program is executed by processor and is joined according to the feature Number determines the step of coordinate correction amount of the grid points in the grid pre-established, comprising: determines that the grid points are sat first The first grid coordinate in mark system determines that the second grid of the grid points in the second coordinate system is sat according to four parameter Mark;The coordinate correction amount of the grid points is determined according to the first grid coordinate and the second grid coordinate.
In one embodiment, also performed the steps of when computer program is executed by processor it is described according to The step of conversion residual error for knowing a little determines the residual error reduction of the grid points, comprising: determine the known point and the grid The distance of point;According to distance described in the conversion residual sum of the known point, the residual of the grid points is determined by distance weighting method Poor reduction.
In one embodiment, it is also performed the steps of when computer program is executed by processor described according to the seat Mark reduction and residual error reduction the step of determining grid reduction, comprising: according to the coordinate correction amount and described residual The sum of poor reduction obtains the grid reduction.
In one embodiment, it is also performed the steps of in the transformation model when computer program is executed by processor It include the grid reduction of multiple grid points;It is described that the lattice are determined according to the coordinate correction amount and the residual error reduction The grid reduction of site, after the step of establishing transformation model according to the grid reduction, further includes: from the modulus of conversion The grid reduction of target mesh point corresponding with the point to be converted is determined in type;Changed according to the grid of the target mesh point Positive quantity determines the reduction of the point to be converted by bilinear interpolation method;According to described to be converted o'clock in the first coordinate system Coordinate and the reduction of the point to be converted determine coordinate of the point to be converted in second coordinate system.
In one embodiment, it is also performed the steps of when computer program is executed by processor described according to known point Before determining the step of first coordinate system is converted to the characteristic parameter of the second coordinate system, further includes: determine that the grid points are pre- If the region to be analyzed in radius, determine that there is first coordinate system and the second coordinate system from the region to be analyzed Coordinate known point.
It will appreciated by the skilled person that realizing all or part of the process in above-described embodiment method, being can It is completed with instructing relevant hardware by computer program, the program can be stored in a computer-readable storage and be situated between In matter, sells or use as independent product.The more specific example (non-exhaustive list) of computer-readable medium includes Below: there is the electrical connection section (electronic device) of one or more wirings, portable computer diskette box (magnetic device), arbitrary access Memory (RAM), read-only memory (ROM), erasable edit read-only storage (EPROM or flash memory), optical fiber dress It sets and portable optic disk read-only storage (CDROM).In addition, computer-readable medium, which can even is that, to be printed on it The paper of described program or other suitable media, because can be for example by carrying out optical scanner to paper or other media, then It edited, interpreted or is handled when necessary with other suitable methods electronically to obtain described program, then by it Storage is in computer storage.
It should be appreciated that each section of the invention can be realized with hardware, software, firmware or their combination.Above-mentioned In embodiment, software that multiple steps or method can be executed in memory and by suitable instruction execution system with storage Or firmware is realized.It, and in another embodiment, can be under well known in the art for example, if realized with hardware Any one of column technology or their combination are realized: having a logic gates for realizing logic function to data-signal Discrete logic, with suitable combinational logic gate circuit specific integrated circuit, programmable gate array (PGA), scene Programmable gate array (FPGA) etc..
The term " includes " of the embodiment of the present invention and " having " and their any deformations, it is intended that cover non-exclusive Include.Such as contain series of steps or the process, method, system, product or equipment of (module) unit are not limited to The step of listing or unit, but optionally further comprising the step of not listing or unit, or optionally further comprising for these The intrinsic other step or units of process, method, product or equipment.
Each technical characteristic of embodiment described above can be combined arbitrarily, for simplicity of description, not to above-mentioned reality It applies all possible combination of each technical characteristic in example to be all described, as long as however, the combination of these technical characteristics is not deposited In contradiction, all should be considered as described in this specification.
The embodiments described above only express several embodiments of the present invention, should not be understood as to the invention patent range Limitation.It should be pointed out that for those of ordinary skill in the art, without departing from the inventive concept of the premise, Various modifications and improvements can be made, and these are all within the scope of protection of the present invention.Therefore, the scope of protection of the patent of the present invention It should be determined by the appended claims.

Claims (10)

1. a kind of method of coordinate conversion, which comprises the following steps:
The characteristic parameter that the first coordinate system is converted to the second coordinate system is determined according to known point;The known point has described the The coordinate of the coordinate of one coordinate system and the second coordinate system;
According to the seat of the grid points in the grid that the characteristic parameter determines the conversion residual error of the known point and pre-establishes Mark reduction;The residual error reduction of the grid points is determined according to the conversion residual error of the known point;
The grid reduction that the grid points are determined according to the coordinate correction amount and the residual error reduction, according to the grid Reduction establishes transformation model;The transformation model is used to determine that the point to be converted in the first coordinate system to be transformed into described second and sits Coordinate in mark system.
2. the method for coordinate conversion according to claim 1, which is characterized in that the characteristic parameter includes four parameters;
It is described that the step of first coordinate system is converted to the characteristic parameter of the second coordinate system is determined according to known point, comprising:
Determine first coordinate of the known point in the first coordinate system and the second coordinate in the second coordinate system, calculating will be described First coordinate is converted to four parameters of the second coordinate.
3. the method for coordinate conversion according to claim 2, which is characterized in that built in advance according to characteristic parameter determination The step of coordinate correction amount of grid points in vertical grid, comprising:
It determines first grid coordinate of the grid points in the first coordinate system, determines that the grid points exist according to four parameter The second grid coordinate in second coordinate system;
The coordinate correction amount of the grid points is determined according to the first grid coordinate and the second grid coordinate.
4. the method for coordinate conversion according to claim 1, which is characterized in that the conversion according to the known point is residual Difference determines the step of residual error reduction of the grid points, comprising:
Determine the known point at a distance from the grid points;According to distance described in the conversion residual sum of the known point, pass through Distance weighting method determines the residual error reduction of the grid points.
5. the method for coordinate conversion according to claim 1, which is characterized in that described according to the coordinate correction amount and institute State the step of residual error reduction determines grid reduction, comprising:
The grid reduction is obtained according to the sum of the coordinate correction amount and the residual error reduction.
6. the method for coordinate conversion according to claim 1, which is characterized in that include multiple lattice in the transformation model The grid reduction of site;
The grid reduction that the grid points are determined according to the coordinate correction amount and the residual error reduction, according to described After the step of grid reduction establishes transformation model, further includes:
The grid reduction of target mesh point corresponding with the point to be converted is determined from the transformation model;
According to the grid reduction of the target mesh point, the reduction of the point to be converted is determined by bilinear interpolation method;
It is determined according to the reduction of the described to be converted o'clock coordinate in the first coordinate system and the point to be converted described wait turn Change the coordinate a little in second coordinate system.
7. the method for coordinate conversion according to any one of claims 1 to 6, which is characterized in that described true according to known point Before fixed the step of first coordinate system is converted to the characteristic parameter of the second coordinate system, further includes:
It determines the region to be analyzed within the scope of the grid points pre-set radius, determines there is described the from the region to be analyzed The known point of the coordinate of one coordinate system and the second coordinate system.
8. a kind of device of coordinate conversion characterized by comprising
Parameter determination module, for determining the characteristic parameter that the first coordinate system is converted to the second coordinate system according to known point;Institute State the coordinate of coordinate and second coordinate system of the known point with first coordinate system;
Reduction determining module, for determining the conversion residual error of the known point according to the characteristic parameter and pre-establishing The coordinate correction amount of grid points in grid;The residual error correction of the grid points is determined according to the conversion residual error of the known point Amount;
And model building module, for determining the grid points according to the coordinate correction amount and the residual error reduction Grid reduction establishes transformation model;The point to be converted that the transformation model is used to determine in the first coordinate system is transformed into described the Coordinate in two coordinate systems.
9. a kind of computer equipment including memory, processor and stores the meter that can be run on a memory and on a processor Calculation machine program, which is characterized in that the processor is realized described in any one of claim 1 to 7 when executing the computer program Method the step of.
10. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the computer program The step of claim 1 to 7 described in any item methods are realized when being executed by processor.
CN201810706859.2A 2018-07-02 2018-07-02 Method, apparatus, computer equipment and the storage medium of coordinate conversion Active CN109191579B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810706859.2A CN109191579B (en) 2018-07-02 2018-07-02 Method, apparatus, computer equipment and the storage medium of coordinate conversion

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810706859.2A CN109191579B (en) 2018-07-02 2018-07-02 Method, apparatus, computer equipment and the storage medium of coordinate conversion

Publications (2)

Publication Number Publication Date
CN109191579A true CN109191579A (en) 2019-01-11
CN109191579B CN109191579B (en) 2019-06-28

Family

ID=64948775

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810706859.2A Active CN109191579B (en) 2018-07-02 2018-07-02 Method, apparatus, computer equipment and the storage medium of coordinate conversion

Country Status (1)

Country Link
CN (1) CN109191579B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114297324A (en) * 2021-11-08 2022-04-08 广东省国土资源测绘院 Coordinate conversion method, equipment, medium and product of geographic data
CN116859421A (en) * 2023-06-27 2023-10-10 国汽大有时空科技(安庆)有限公司 Positioning service method and device for multiple reference frames
CN117708960A (en) * 2024-02-04 2024-03-15 武汉大学 Real-time conversion method, device, equipment and medium for plane coordinates and normal height

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102339478A (en) * 2011-09-14 2012-02-01 北京地拓科技发展有限公司 Method and device for generating digital elevation model from contour map
CN102779231A (en) * 2012-06-18 2012-11-14 郑州大学 Coordinate conversion parameter calculating method based on adjacent control point sets
CN102853793A (en) * 2012-09-27 2013-01-02 中国科学院高能物理研究所 Coordinate transformation data processing method and coordinate transformation data processing device
CN103017653A (en) * 2012-11-27 2013-04-03 武汉海达数云技术有限公司 Registration and measurement method of spherical panoramic image and three-dimensional laser scanning point cloud
CN104537606A (en) * 2014-12-09 2015-04-22 中国农业大学 Geographical coordinate projection transformation method
CN105300373A (en) * 2015-09-16 2016-02-03 广州市城市规划勘测设计研究院 Three-dimensional coordinate conversion method and device

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102339478A (en) * 2011-09-14 2012-02-01 北京地拓科技发展有限公司 Method and device for generating digital elevation model from contour map
CN102779231A (en) * 2012-06-18 2012-11-14 郑州大学 Coordinate conversion parameter calculating method based on adjacent control point sets
CN102853793A (en) * 2012-09-27 2013-01-02 中国科学院高能物理研究所 Coordinate transformation data processing method and coordinate transformation data processing device
CN103017653A (en) * 2012-11-27 2013-04-03 武汉海达数云技术有限公司 Registration and measurement method of spherical panoramic image and three-dimensional laser scanning point cloud
CN104537606A (en) * 2014-12-09 2015-04-22 中国农业大学 Geographical coordinate projection transformation method
CN105300373A (en) * 2015-09-16 2016-02-03 广州市城市规划勘测设计研究院 Three-dimensional coordinate conversion method and device

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114297324A (en) * 2021-11-08 2022-04-08 广东省国土资源测绘院 Coordinate conversion method, equipment, medium and product of geographic data
CN116859421A (en) * 2023-06-27 2023-10-10 国汽大有时空科技(安庆)有限公司 Positioning service method and device for multiple reference frames
CN116859421B (en) * 2023-06-27 2024-03-15 国汽大有时空科技(安庆)有限公司 Positioning service method and device for multiple reference frames
CN117708960A (en) * 2024-02-04 2024-03-15 武汉大学 Real-time conversion method, device, equipment and medium for plane coordinates and normal height
CN117708960B (en) * 2024-02-04 2024-05-03 武汉大学 Real-time conversion method, device, equipment and medium for plane coordinates and normal height

Also Published As

Publication number Publication date
CN109191579B (en) 2019-06-28

Similar Documents

Publication Publication Date Title
CN109191579B (en) Method, apparatus, computer equipment and the storage medium of coordinate conversion
CN109814572A (en) Localization for Mobile Robot builds drawing method, device, mobile robot and storage medium
CN109189868A (en) Method, apparatus, computer equipment and the storage medium of coordinate conversion
CN107093206A (en) The method modeled using 3D laser scanner technique fast Bs IM
CN103971366B (en) A kind of solid matching method being polymerize based on double weights
CN110880717B (en) Automatic cable laying method based on equipment and bridge positioning digitization
CN103336894B (en) A kind of city plot ratio automatic Mesh Partition Method
CN106780262A (en) A kind of same bit pattern for considering urban road network constraint finds method and device
CN111028335B (en) Point cloud data block surface patch reconstruction method based on deep learning
CN110246095A (en) A kind of image repair method, device, storage medium and terminal device
CN109167348A (en) Method is determined based on the dominant eigenvalues feasible zone of multi-parametric programming theory
CN106023317B (en) A kind of weighted Voronoi diagrams drawing generating method for big data test
CN105741340A (en) Transmission line 3D scene simulation method and transmission line 3D scene simulation system for webpage display
CN109934866A (en) A kind of extracting method, copy method and the device of figure spot partition mode
CN113051637A (en) BIM modeling method for public medical building
CN117610354A (en) Method and device for mapping parallelepiped grid
CN108921448A (en) Power trade method, apparatus, equipment and storage medium
CN109785422A (en) The construction method and device of three-dimensional power grid scene
CN105809717B (en) A kind of depth estimation method, system and electronic equipment
CN117172134A (en) Moon surface multiscale DEM modeling method and system based on converged terrain features
CN104598678A (en) Generation method for logic geographic diagram of power transmission network
CN110489510B (en) Road data processing method and device, readable storage medium and computer equipment
CN102324096B (en) Method and device for splicing characteristic point correction partitions
CN102117495B (en) New data partitioning method and device in three-dimensional modeling
Galanda Optimization techniques for polygon generalization

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CB03 Change of inventor or designer information

Inventor after: Fu Yan

Inventor after: Wang Jianhui

Inventor after: Li Genger

Inventor after: Han Fei

Inventor after: Wu Lvchuan

Inventor before: Wang Jianhui

Inventor before: Fu Yan

Inventor before: Li Genger

Inventor before: Han Fei

Inventor before: Wu Lvchuan

CB03 Change of inventor or designer information