CN109189868A - 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
CN109189868A
CN109189868A CN201810719377.0A CN201810719377A CN109189868A CN 109189868 A CN109189868 A CN 109189868A CN 201810719377 A CN201810719377 A CN 201810719377A CN 109189868 A CN109189868 A CN 109189868A
Authority
CN
China
Prior art keywords
coordinate
grid
coordinate system
point
evaluation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201810719377.0A
Other languages
Chinese (zh)
Inventor
符彦
王剑辉
杨映新
韩菲
伍素贞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
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 CN201810719377.0A priority Critical patent/CN109189868A/en
Publication of CN109189868A publication Critical patent/CN109189868A/en
Pending legal-status Critical Current

Links

Landscapes

  • Processing Or Creating Images (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 evaluation distance between the grid points in multiple known points and the grid pre-established;Coordinate evaluation coefficient is determined according to evaluation distance, and the assessment coordinate of grid points in the second coordinate system is determined according to the coordinate of coordinate evaluation coefficient and the multiple known point in the second coordinate system;The grid reduction of the grid points is determined according to coordinate of the grid points in the first coordinate system and assessment coordinate;Transformation model is established according to 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.Spatial coherence is accounted for, high-precision coordinate conversion is realized.

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 estimate sky Between correlation, obtain accurately converting coordinate.
The content of the embodiment of the present invention is as follows:
A kind of method of coordinate conversion, comprising the following steps: determine the lattice in multiple known points and the grid pre-established Evaluation distance between site;The grid points have the coordinate of the first coordinate system, and the known point has the second coordinate system Coordinate;Coordinate evaluation coefficient is determined according to the evaluation distance, is existed according to the coordinate evaluation coefficient and the multiple known point Coordinate in second coordinate system determines the assessment coordinate of the grid points in the second coordinate system;According to the grid points first Coordinate and the assessment coordinate in coordinate system determine the grid reduction of the grid points;It is established according to the grid reduction Transformation model;The transformation model is used to determine that the point to be converted in the first coordinate system to be transformed into the seat in second coordinate system Mark.
Described the step of coordinate evaluation coefficient is determined according to the evaluation distance in one of the embodiments, comprising: root Coordinate evaluation coefficient is determined by Kriging regression method according to the evaluation distance.
It is described in one of the embodiments, that coordinate evaluation system is determined by Kriging regression method according to the evaluation distance Before several steps, further includes: obtain the reference distance between known point two-by-two and refer to semivariance, determine the reference distance With the matched curve with reference to semivariance.
It is described in one of the embodiments, that coordinate evaluation system is determined by Kriging regression method according to the evaluation distance Several step, comprising: the multiple known point and the grid points are determined according to the evaluation distance and the matched curve Candidate semivariance;The coordinate evaluation coefficient is determined according to the candidate semivariance and the semivariance that refers to.
It is described in one of the embodiments, that the coordinate is determined according to the candidate semivariance and the semivariance that refers to The step of evaluation coefficient, comprising: the candidate semivariance and the reference semivariance are substituted into Kriging regression equation group, asked The optimal weights coefficient for solving the Kriging regression equation group obtains the coordinate evaluation coefficient.
In one of the embodiments, it is described according to the coordinate evaluation coefficient and the multiple known point in the second coordinate Coordinate in system determines the step of assessment coordinate of the grid points in the second coordinate system, comprising: the coordinate is evaluated system The coordinate of several and the multiple known point in the second coordinate system substitutes into the Kriging regression equation group, and solution obtains described The assessment coordinate of grid points in the second coordinate system.
It in one of the embodiments, include the grid reduction of multiple grid points in the transformation model;Described After the step of establishing transformation model according to the grid reduction, further includes: determining and described wait turn from the transformation model Change the grid reduction of a little corresponding target mesh point;According to the grid reduction of the target mesh point, by bilinearity Slotting method determines the reduction of the point to be converted;According to described to be converted o'clock coordinate in the first coordinate system and described wait turn It changes reduction a little and determines coordinate of the point to be converted in second coordinate system.
Correspondingly, the embodiment of the present invention provides a kind of device of coordinate conversion, comprising: apart from determining module, for determining Evaluation distance between grid points in multiple known points and the grid pre-established;The grid points have the first coordinate system Coordinate, the known point have the coordinate of the second coordinate system;Coordinate determining module, for determining coordinate according to the evaluation distance Evaluation coefficient determines the grid according to the coordinate evaluation coefficient and the coordinate of the multiple known point in the second coordinate system The assessment coordinate of point in the second coordinate system;Reduction determining module, for according to the grid points in the first coordinate system Coordinate and the assessment coordinate determine the grid reduction of the grid points;And model building module, for according to the lattice Net reduction establishes transformation model;The transformation model is for determining that the point to be converted in the first coordinate system is transformed into described second Coordinate in coordinate system.
The method and device of above-mentioned coordinate conversion determines coordinate evaluation according to the distance between multiple known points and grid points Coefficient, and then determine grid points assessment coordinate in the second coordinate system, according to coordinate of the grid points in the first coordinate system and Assessment coordinate just can determine that the grid reduction of grid points, and establish for carrying out between the first coordinate system and the second coordinate system The transformation model of coordinate conversion.Grid reduction is solved on the basis of considering the evaluation distance between known point and grid points, Spatial coherence is accounted for, high-precision coordinate conversion is realized.
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 multiple known points and builds in advance Evaluation distance between grid points in vertical grid;The grid points have the coordinate of the first coordinate system, the known point tool There is the coordinate of the second coordinate system;Coordinate evaluation coefficient is determined according to the evaluation distance, according to the coordinate evaluation coefficient and institute It states the coordinate of multiple known points in the second coordinate system and determines the assessment coordinate of the grid points in the second coordinate system;According to institute State the grid reduction that coordinate and the assessment coordinate of the grid points in the first coordinate system determine the grid points;According to described 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.
Above-mentioned computer equipment solves grid correction on the basis of considering the evaluation distance between known point and grid points Amount, accounts for spatial coherence, 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 evaluation distance determined between the grid points in multiple known points and the grid pre-established is performed the steps of when row;It is described Grid points have the coordinate of the first coordinate system, and the known point has the coordinate of the second coordinate system;Really according to the evaluation distance Position fixing evaluation coefficient determines institute according to the coordinate evaluation coefficient and the coordinate of the multiple known point in the second coordinate system State the assessment coordinate of grid points in the second coordinate system;According to coordinate of the grid points in the first coordinate system and the assessment Coordinate determines the grid reduction of the grid points;Transformation model is established according to the grid reduction;The transformation model is used In determining that the point to be converted in the first coordinate system is transformed into the coordinate in second coordinate system.
Above-mentioned computer readable storage medium solves on the basis of considering the evaluation distance between known point and grid points Grid reduction, accounts for spatial coherence, realizes 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 schematic diagram for determining grid reduction in one embodiment according to known point;
Fig. 5 is the three dimensional structure diagram of transformation model in one embodiment;
Fig. 6 projects to the schematic diagram of grid for point to be converted 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 for store coordinate evaluation coefficient, coordinate, grid reduction, The data such as transformation model.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, evaluation distance between grid points in multiple known points and the grid pre-established is determined;The grid Point has the coordinate of the first coordinate system, and the known point has the coordinate of the second coordinate system.
In this step, it is known that put the coordinate of the coordinate and the second coordinate system with the first coordinate system, grid points have the The coordinate of one coordinate system, therefore, according to the coordinate of known point and grid points in the first coordinate system it become convenient substantially that calculating The distance between they, i.e. evaluation distance.
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.
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.
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.
S202, coordinate evaluation coefficient is determined according to evaluation distance, according to the coordinate evaluation coefficient and the multiple The coordinate of known point in the second coordinate system determines the assessment coordinate of the grid points in the second coordinate system.
This step evaluates distance according to determined by previous step and determines coordinate evaluation coefficient, and determines grid points second Assessment coordinate in coordinate system.
Wherein, coordinate evaluation coefficient is referred to for evaluating the transformational relation between the first coordinate system and the second coordinate system Coefficient can refer to weight etc..In this step, the determination of coordinate evaluation coefficient is realized according to evaluation distance, can be passed through Kriging regression method or other calculation method carry out.
S203, the grid points are determined according to coordinate of the grid points in the first coordinate system and the assessment coordinate Grid reduction.
Assessment coordinate refer to according to preceding method the coordinate of calculated grid points in the second coordinate system, therefore, In this step the coordinate by grid points in the first coordinate system and the assessment coordinate be compared (such as: seeking difference) can be true Make the grid reduction of each site.
Grid reduction refers to the reduction of some grid points, grid points be distributed in region to be converted eachly Side, each point to be converted can determine the reduction of oneself according to neighbouring grid points.
S204, transformation model is established according to the grid reduction;The transformation model is for determining in the first coordinate system Point to be converted be transformed into the coordinate in second coordinate system.
In this step, transformation model is established according to grid reduction, according to the transformation model and to be converted o'clock Coordinate in one coordinate system just can determine that the coordinate of point to be converted in the 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.
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 solves grid reduction on the basis of considering the evaluation distance between known point and grid points, cares for And spatial coherence, realize high-precision coordinate conversion.
In one embodiment, described the step of coordinate evaluation coefficient is determined according to the evaluation distance, comprising: according to institute Commentary valence distance determines coordinate evaluation coefficient by Kriging regression method.
Wherein, Kriging regression method is a kind of interpolation method based on statistics, is that one kind is based on auto-correlation (that is, measurement point Between statistical relationship) statistical model.Therefore, which not only has the function of predicting surface, but also can be to pre- The certainty or accuracy of survey provide certain measurement.Kriging regression method assumes that the distance between sampled point or direction can reflect The spatial coherence of surface variation.Kriging regression method can will own in the point or specified radius of mathematical function and specified quantity Point is fitted, with the output valve of each point of determination.Meanwhile Kriging regression method is multistep process: the spy including data Statistical analysis, variation function modeling, creation surface and research variance surface etc. without hesitation.
In gram in golden differential technique any point interpolation resultSee following formula (golden difference value equation in gram):
Wherein, Z (si) indicate ith measurement value;λiIndicate the weight of ith measurement value;s0Indicate predicted position;N is indicated The total number of measured value.
The present embodiment determines coordinate evaluation coefficient by Kriging regression method, can fully consider that grid points and surrounding are multiple Know the distance between a little, and then can determine the grid reduction of grid points according to required coordinate evaluation coefficient.
In one embodiment, described that coordinate evaluation coefficient is determined by Kriging regression method according to the evaluation distance Before step, further includes: obtain the reference distance between known point two-by-two and refer to semivariance, determine the reference distance and institute State the matched curve with reference to semivariance.
Wherein it is possible to calculate the reference distance two-by-two of all known points in region to be converted and refer to semivariance, and according to These reference distances and reference semivariance are fitted, and find out a matched curve.So that in the case where known any distance, Corresponding semivariance can be calculated according to the matched curve.
The present embodiment is accomplished that the determination of matched curve, includes multiple reference distances and reference half in the matched curve Therefore variance just can determine that out corresponding candidate semivariance, in turn after knowing the evaluation distance between known point and grid points Coordinate evaluation coefficient can be determined according to candidate's semivariance.
In one embodiment, described that coordinate evaluation coefficient is determined by Kriging regression method according to the evaluation distance Step, comprising: the time of the multiple known point Yu the grid points is determined according to the evaluation distance and the matched curve Select semivariance;The coordinate evaluation coefficient is determined according to the candidate semivariance and the semivariance that refers to.
According to the description of a upper embodiment it is found that it is known evaluation apart from when corresponding candidate can be found according to matched curve Semivariance.The present embodiment determines coordinate evaluation coefficient according to these candidate semivariances and with reference to semivariance.
Wherein, according to candidate semivariance and with reference to semivariance determine coordinate evaluation coefficient can by Kriging regression method or Person's others algorithms is realized.
The present embodiment determines coordinate evaluation coefficient by Kriging regression method, can fully consider that grid points and surrounding are multiple Know the distance between a little, and then can determine the grid reduction of grid points according to required coordinate evaluation coefficient.
In one embodiment, described to determine that the coordinate is evaluated according to the candidate semivariance and the semivariance that refers to The step of coefficient, comprising: the candidate semivariance and the reference semivariance are substituted into Kriging regression equation group, institute is solved The optimal weights coefficient for stating Kriging regression equation group obtains the coordinate evaluation coefficient.
According to the formula of Kriging regression method it is found that substituting into Kriging regression equation by candidate semivariance and with reference to semivariance After in group, the weight in formula be it is unknown, solve system of equation just can determine that out optimal weight coefficient, and then obtain coordinate Evaluation coefficient.The present embodiment solves coordinate evaluation coefficient by solving the formula of Kriging regression method, has fully considered coordinate Relationship between evaluation coefficient and evaluation distance, semivariance.Therefore, the evaluation distance between known point and grid points is being determined When, coordinate evaluation coefficient can be used to accurately evaluate the coordinate relationship between known point and grid points.
In one embodiment, it is described according to the coordinate evaluation coefficient and the multiple known point in the second coordinate system Coordinate the step of determining the assessment coordinate of the grid points in the second coordinate system, comprising: by the coordinate evaluation coefficient and The coordinate of the multiple known point in the second coordinate system substitutes into the Kriging regression equation group, and solution obtains the grid The assessment coordinate of point in the second coordinate system.
Coordinate evaluation coefficient can be used to accurately evaluate the coordinate relationship between known point and grid points, therefore the present embodiment Assessment coordinate is determined according to the coordinate of coordinate evaluation coefficient and known point in the second coordinate system.Method of determination is equally that will count According to substituting into Kriging regression formula, the assessment coordinate of grid points in the second coordinate system is determined by solving equation.
In one embodiment, as shown in figure 4, determining grid points (being connected with the point of a plurality of straight line in figure) and surrounding respectively Evaluation distance between all known points determines coordinate evaluation coefficient according to these evaluation distances, according to coordinate evaluation coefficient It determines the assessment coordinate of grid points in the second coordinate system and then determines grid reduction.
It in one embodiment, include the grid reduction of multiple grid points in the transformation model;It is described according to institute After stating the step of grid reduction establishes transformation model, further includes: determining with the point to be converted from the transformation model The grid reduction of corresponding target mesh point;According to the grid reduction of the target mesh point, pass through bilinear interpolation method Determine the reduction of the point to be converted;According to the described to be converted o'clock coordinate and the point to be converted in the first coordinate system Reduction determine 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 5, in Fig. 5 different coordinates Grid reduction may be different.Grid reduction in Fig. 5 is bigger to illustrate that deviation is bigger, and Fig. 5 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 be converted o'clock in the first coordinate system and point to be converted determines point to be converted The realization process of coordinate in the second coordinate system can be with are as follows: by the to be converted o'clock coordinate in the first coordinate system plus to be converted The coordinate of point to be converted in the second coordinate system can be obtained in the reduction of point.
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. 6) as shown in Figure 6, 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, as shown in fig. 7, providing a kind of method of coordinate conversion, comprising the following steps:
S701, evaluation distance between grid points in multiple known points and the grid pre-established is determined.
S702, it obtains the reference distance between known point two-by-two and refers to semivariance, determine the reference distance and described With reference to the matched curve of semivariance.
S703, the multiple known point and the grid points are determined according to the evaluation distance and the matched curve Candidate semivariance;The candidate semivariance and the reference semivariance are substituted into Kriging regression equation group, solve described gram In golden interpolation equation group optimal weights coefficient, obtain coordinate evaluation coefficient.
S704, the coordinate of the coordinate evaluation coefficient and the multiple known point in the second coordinate system is substituted into described gram In in golden interpolation equation group, solution obtains the assessment coordinate of the grid points in the second coordinate system.
S705, the grid points are determined according to coordinate of the grid points in the first coordinate system and the assessment coordinate Grid reduction.
S706, transformation model is established according to the grid reduction;From the transformation model determine with it is described to be converted The grid reduction of the corresponding target mesh point of point.
S707, 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.
The present embodiment solves grid reduction on the basis of considering the evaluation distance between known point and grid points, cares for And spatial coherence, realize high-precision coordinate conversion.
The application of the above method in order to better understand, the method for a coordinate conversion of the present invention detailed below is real Example.The application example is accomplished that the coordinate between 1980 coordinate of Xi'an and CGCS2000 coordinate is converted, wherein region to be converted Inside there is known 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 in area map, calculate the reference distance between known point two-by-two and refer to semi-square Difference determines these reference distances and the matched curve with reference to semivariance;
S32, the evaluation distance for determining each known point Yu grid points A are determined each according to the evaluation distance and matched curve The candidate semivariance of a known point and grid points A;
S33, it is substituted into Kriging regression equation group by candidate semivariance and with reference to semivariance, the optimal power of solve system of equation Weight coefficient, obtains coordinate evaluation coefficient;
S34, the coordinate of required coordinate evaluation coefficient and each known point in the second coordinate system is substituted into the Ke Lijin In interpolation equation group, solution obtains the assessment coordinate of grid points A in the second coordinate system;
S35, the grid that grid points A is determined according to coordinate of the grid points A in the first coordinate system and required assessment coordinate Reduction;
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 true apart from determining module 801, coordinate determining module 802, reduction Cover half block 803 and model building module 804, detailed description are as follows:
Apart from determining module 801, for determining commenting between the grid points in multiple known points and the grid pre-established Valence distance;The grid points have the coordinate of the first coordinate system, and the known point has the coordinate of the second coordinate system.
Coordinate determining module 802 is evaluated for determining coordinate evaluation coefficient according to the evaluation distance according to the coordinate The coordinate of coefficient and the multiple known point in the second coordinate system determines that the assessment of the grid points in the second coordinate system is sat Mark.
Reduction determining module 803, for the coordinate and assessment seat according to the grid points in the first coordinate system Mark determines the grid reduction of the grid points.
And model building module 804, for establishing 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 present embodiment solves grid reduction on the basis of considering the evaluation distance between known point and grid points, cares for And spatial coherence, realize high-precision coordinate conversion.
In one embodiment, coordinate evaluation module 802 is also used to pass through Kriging regression method according to the evaluation distance Determine coordinate evaluation coefficient.
In one embodiment, further includes: curve fitting module, for obtain the reference distance between known point two-by-two and With reference to semivariance, the reference distance and the matched curve with reference to semivariance are determined.
In one embodiment, coordinate evaluation module 802, comprising: semivariance determines submodule, for according to the evaluation Distance and the matched curve determine the candidate semivariance of the multiple known point Yu the grid points;Evaluation coefficient determines son Module, for determining the coordinate evaluation coefficient according to the candidate semivariance and the semivariance that refers to.
In one embodiment, evaluation coefficient determines submodule, is also used to the candidate semivariance and the reference half Variance substitutes into Kriging regression equation group, solves the optimal weights coefficient of the Kriging regression equation group, obtains the seat Mark evaluation coefficient.
In one embodiment, coordinate evaluation module 802 is also used to the coordinate evaluation coefficient and the multiple known The coordinate of point in the second coordinate system substitutes into the Kriging regression equation group, and solution obtains the grid points in the second coordinate Assessment coordinate in system.
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.
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 performs the steps of determination when executing computer program multiple Know the evaluation distance a little between the grid points in the grid pre-established;The grid points have the coordinate of the first coordinate system, The known point has the coordinate of the second coordinate system;Coordinate evaluation coefficient is determined according to the evaluation distance, according to the coordinate The coordinate of evaluation coefficient and the multiple known point in the second coordinate system determines the grid points commenting in the second coordinate system Estimate coordinate;Determine that the grid of the grid points changes according to coordinate of the grid points in the first coordinate system and the assessment coordinate Positive quantity;Transformation model is established according to the grid reduction;The transformation model is to be converted in the first coordinate system for determining The coordinate that point is transformed into second coordinate system.
In one embodiment, it is also performed the steps of when processor executes computer program described according to the evaluation Distance determines the step of coordinate evaluation coefficient, comprising: determines that coordinate is evaluated by Kriging regression method according to the evaluation distance Coefficient.
In one embodiment, it is also performed the steps of when processor executes computer program described according to the evaluation Before the step of distance determines coordinate evaluation coefficient by Kriging regression method, further includes: obtain the ginseng between known point two-by-two It examines distance and with reference to semivariance, determines the reference distance and the matched curve with reference to semivariance.
In one embodiment, it is also performed the steps of when processor executes computer program described according to the evaluation The step of distance determines coordinate evaluation coefficient by Kriging regression method, comprising: according to the evaluation distance and the fitting Curve determines the candidate semivariance of the multiple known point Yu the grid points;According to the candidate semivariance and the reference half Variance determines the coordinate evaluation coefficient.
In one embodiment, it is also performed the steps of when processor executes computer program described according to the candidate Semivariance and described the step of determining the coordinate evaluation coefficient with reference to semivariance, comprising: by the candidate semivariance and described It is substituted into Kriging regression equation group with reference to semivariance, solves the optimal weights coefficient of the Kriging regression equation group, obtain The coordinate evaluation coefficient.
In one embodiment, it is also performed the steps of when processor executes computer program described according to the coordinate The coordinate of evaluation coefficient and the multiple known point in the second coordinate system determines the grid points commenting in the second coordinate system The step of estimating coordinate, comprising: substitute into the coordinate of the coordinate evaluation coefficient and the multiple known point in the second coordinate system In the Kriging regression equation group, solution obtains the assessment coordinate of the grid points in the second coordinate system.
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;After described the step of establishing transformation model according to the grid reduction, also wrap It includes: determining the grid reduction of target mesh point corresponding with the point to be converted from the transformation model;According to the mesh The grid reduction of style site determines the reduction of the point to be converted by bilinear interpolation method;According to described to be converted The reduction of coordinate and the point to be converted o'clock in the first coordinate system determines the point to be converted in second coordinate Coordinate in system.
In one embodiment, a kind of computer readable storage medium is provided, computer program is stored thereon with, is calculated Machine program is performed the steps of when being executed by processor between the grid points determined in multiple known points and the grid pre-established Evaluation distance;The grid points have the coordinate of the first coordinate system, and the known point has the coordinate of the second coordinate system;According to The evaluation distance determines coordinate evaluation coefficient, according to the coordinate evaluation coefficient and the multiple known point in the second coordinate system In coordinate determine the assessment coordinate of the grid points in the second coordinate system;According to the grid points in the first coordinate system Coordinate and the assessment coordinate determine the grid reduction of the grid points;Transformation model is established 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.
In one embodiment, it is also performed the steps of when computer program is executed by processor described according to institute's commentary The step of valence distance determines coordinate evaluation coefficient, comprising: determine that coordinate is commented by Kriging regression method according to the evaluation distance Valence coefficient.
In one embodiment, it is also performed the steps of when computer program is executed by processor described according to institute's commentary Before the step of valence distance determines coordinate evaluation coefficient by Kriging regression method, further includes: obtain between known point two-by-two Reference distance and refer to semivariance, determine the reference distance and the matched curve with reference to semivariance.
In one embodiment, it is also performed the steps of when computer program is executed by processor described according to institute's commentary Valence distance the step of coordinate evaluation coefficient is determined by Kriging regression method, comprising: according to the evaluation apart from and described intend Close the candidate semivariance that curve determines the multiple known point Yu the grid points;According to the candidate semivariance and the reference Semivariance determines the coordinate evaluation coefficient.
In one embodiment, it is also performed the steps of when computer program is executed by processor described according to the time Select semivariance and described the step of determining the coordinate evaluation coefficient with reference to semivariance, comprising: by the candidate semivariance and institute It states and is substituted into Kriging regression equation group with reference to semivariance, solve the optimal weights coefficient of the Kriging regression equation group, obtain To the coordinate evaluation coefficient.
In one embodiment, it is also performed the steps of when computer program is executed by processor described according to the seat Mark evaluation coefficient and the multiple known point coordinate in the second coordinate system determines the grid points in the second coordinate system The step of assessing coordinate, comprising: by the coordinate generation of the coordinate evaluation coefficient and the multiple known point in the second coordinate system Enter in the Kriging regression equation group, solution obtains the assessment coordinate of the grid points in the second coordinate system.
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;After described the step of establishing transformation model according to the grid reduction, also It include: the grid reduction that target mesh point corresponding with the point to be converted is determined from the transformation model;According to described The grid reduction of target mesh point determines the reduction of the point to be converted by bilinear interpolation method;According to described wait turn The reduction for changing coordinate and the point to be converted o'clock in the first coordinate system determines that the point to be converted is sat described second Coordinate in mark system.
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:
Determine the evaluation distance between the grid points in multiple known points and the grid pre-established;The grid points have first The coordinate of coordinate system, the known point have the coordinate of the second coordinate system;
Coordinate evaluation coefficient is determined according to evaluation distance, according to the coordinate evaluation coefficient and the multiple known point the Coordinate in two coordinate systems determines the assessment coordinate of the grid points in the second coordinate system;
The grid correction of the grid points is determined according to coordinate of the grid points in the first coordinate system and the assessment coordinate Amount;
Transformation model is established according to the grid reduction;The transformation model is used to determine the point to be converted in the first coordinate system The coordinate being transformed into second coordinate system.
2. the method for coordinate conversion according to claim 1, which is characterized in that described determined according to the evaluation distance is sat The step of marking evaluation coefficient, comprising:
Coordinate evaluation coefficient is determined by Kriging regression method according to the evaluation distance.
3. the method for coordinate according to claim 2 conversion, which is characterized in that described to be passed through gram according to the evaluation distance In before golden interpolation method the step of determining coordinate evaluation coefficient, further includes:
It obtains the reference distance between known point two-by-two and refers to semivariance, determine the reference distance and described with reference to semivariance Matched curve.
4. the method for coordinate according to claim 3 conversion, which is characterized in that described to be passed through gram according to the evaluation distance In golden interpolation method the step of determining coordinate evaluation coefficient, comprising:
The candidate semi-square of the multiple known point Yu the grid points is determined according to the evaluation distance and the matched curve Difference;
The coordinate evaluation coefficient is determined according to the candidate semivariance and the semivariance that refers to.
5. the method for coordinate conversion according to claim 4, which is characterized in that described according to the candidate semivariance and institute State the step of determining the coordinate evaluation coefficient with reference to semivariance, comprising:
The candidate semivariance and the reference semivariance are substituted into Kriging regression equation group, the Kriging regression is solved The optimal weights coefficient of equation group obtains the coordinate evaluation coefficient.
6. the method for coordinate according to claim 5 conversion, which is characterized in that it is described according to the coordinate evaluation coefficient and The coordinate of the multiple known point in the second coordinate system determines the step of the assessment coordinate of the grid points in the second coordinate system Suddenly, comprising:
The coordinate of the coordinate evaluation coefficient and the multiple known point in the second coordinate system is substituted into the Kriging regression In equation group, solution obtains the assessment coordinate of the grid points in the second coordinate system.
7. the method for coordinate conversion according to any one of claims 1 to 6, which is characterized in that wrapped in the transformation model Grid reduction containing multiple grid points;
After described the step of establishing transformation model according to the grid reduction, 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.
8. a kind of device of coordinate conversion characterized by comprising
Apart from determining module, for determining the evaluation distance between the grid points in multiple known points and the grid pre-established; The grid points have the coordinate of the first coordinate system, and the known point has the coordinate of the second coordinate system;
Coordinate determining module, for determining coordinate evaluation coefficient according to evaluation distance, according to the coordinate evaluation coefficient with The coordinate of the multiple known point in the second coordinate system determines the assessment coordinate of the grid points in the second coordinate system;
Reduction determining module, for determining institute according to coordinate of the grid points in the first coordinate system and the assessment coordinate State the grid reduction of grid points;
And model building module, for establishing 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.
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.
CN201810719377.0A 2018-07-02 2018-07-02 Method, apparatus, computer equipment and the storage medium of coordinate conversion Pending CN109189868A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810719377.0A CN109189868A (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
CN201810719377.0A CN109189868A (en) 2018-07-02 2018-07-02 Method, apparatus, computer equipment and the storage medium of coordinate conversion

Publications (1)

Publication Number Publication Date
CN109189868A true CN109189868A (en) 2019-01-11

Family

ID=64948860

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810719377.0A Pending CN109189868A (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) CN109189868A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111723167A (en) * 2019-05-10 2020-09-29 腾讯科技(深圳)有限公司 Map coordinate conversion method, map coordinate conversion device, readable storage medium and computer equipment
CN114297324A (en) * 2021-11-08 2022-04-08 广东省国土资源测绘院 Coordinate conversion method, equipment, medium and product of geographic data
CN114399541A (en) * 2021-12-29 2022-04-26 北京师范大学 Regional coordinate conversion method and device
CN117708960A (en) * 2024-02-04 2024-03-15 武汉大学 Real-time conversion method, device, equipment and medium for plane coordinates and normal height

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103646395A (en) * 2013-11-28 2014-03-19 中南大学 A high-precision image registering method based on a grid method
CN106484997A (en) * 2016-10-10 2017-03-08 南京林业大学 A kind of waterfront band silt depth based on Kriging regression calculates and goes out drawing method
CN107391578A (en) * 2017-06-20 2017-11-24 国家测绘地理信息局海南基础地理信息中心 A kind of Map Service of Network dynamic coordinate conversion method based on grid method
US20180046652A1 (en) * 2016-02-02 2018-02-15 Southeast University A definition method for urban dynamic spatial structure circle

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103646395A (en) * 2013-11-28 2014-03-19 中南大学 A high-precision image registering method based on a grid method
US20180046652A1 (en) * 2016-02-02 2018-02-15 Southeast University A definition method for urban dynamic spatial structure circle
CN106484997A (en) * 2016-10-10 2017-03-08 南京林业大学 A kind of waterfront band silt depth based on Kriging regression calculates and goes out drawing method
CN107391578A (en) * 2017-06-20 2017-11-24 国家测绘地理信息局海南基础地理信息中心 A kind of Map Service of Network dynamic coordinate conversion method based on grid method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
张成校等: "格网内插法坐标转换", 《全球定位系统》 *
段金梅等: "基于格网的坐标转换方法问题研究", 《测绘技术装备》 *
王志杰等: "大范围区域内地心坐标系与参心坐标系的格网转换方法", 《测绘通报》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111723167A (en) * 2019-05-10 2020-09-29 腾讯科技(深圳)有限公司 Map coordinate conversion method, map coordinate conversion device, readable storage medium and computer equipment
CN114297324A (en) * 2021-11-08 2022-04-08 广东省国土资源测绘院 Coordinate conversion method, equipment, medium and product of geographic data
CN114399541A (en) * 2021-12-29 2022-04-26 北京师范大学 Regional coordinate conversion method and device
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

Similar Documents

Publication Publication Date Title
CN109189868A (en) Method, apparatus, computer equipment and the storage medium of coordinate conversion
Zamzam et al. Physics-aware neural networks for distribution system state estimation
Gyulassy et al. Efficient computation of Morse-Smale complexes for three-dimensional scalar functions
CN103139907B (en) A kind of indoor wireless positioning method utilizing fingerprint technique
CN110880717B (en) Automatic cable laying method based on equipment and bridge positioning digitization
CN102819568B (en) Submarine topography data establishment method based on topographical sampling point positions
CN107093206A (en) The method modeled using 3D laser scanner technique fast Bs IM
CN108763825B (en) Numerical simulation method for simulating wind field of complex terrain
CN109191579B (en) Method, apparatus, computer equipment and the storage medium of coordinate conversion
CN113532397B (en) Regional elevation anomaly fitting method based on expansion algorithm
CN104252549A (en) Analytical well spacing method based on kriging interpolation
CN104573234B (en) A kind of basin storage capacity density space method based on ArcGIS
CN108242962B (en) Indoor signal propagation loss calculation method and device based on measurement report
CN109887073A (en) The three-dimensional digital model construction method and device of rock core
Chaturvedi et al. Solar potential analysis and integration of the time-dependent simulation results for semantic 3D city models using dynamizers
CN117077922A (en) BP neural network-based intelligent power grid power distribution method
CN108769996A (en) Beaconing nodes layout method and localization method used for positioning
Tavelli et al. A high‐order parallel Eulerian‐Lagrangian algorithm for advection‐diffusion problems on unstructured meshes
CN106658538B (en) Mobile phone base station signal coverage area simulation method based on Thiessen polygon
Gosciewski Reduction of deformations of the digital terrain model by merging interpolation algorithms
Vázquez‐Delgado et al. Real‐time multi‐window stereo matching algorithm with fuzzy logic
Molinero et al. How the geometry of cities explains urban scaling laws and determines their exponents
CN113051637A (en) BIM modeling method for public medical building
CN108563915A (en) Vehicle digitizes emulation testing model construction system and method, computer program
CN108053481A (en) Generation method, device and the storage medium of three-dimensional point cloud normal vector

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20190111

RJ01 Rejection of invention patent application after publication