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 PDFInfo
- 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
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
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.
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)
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)
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 |
-
2018
- 2018-07-02 CN CN201810719377.0A patent/CN109189868A/en active Pending
Patent Citations (4)
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)
Title |
---|
张成校等: "格网内插法坐标转换", 《全球定位系统》 * |
段金梅等: "基于格网的坐标转换方法问题研究", 《测绘技术装备》 * |
王志杰等: "大范围区域内地心坐标系与参心坐标系的格网转换方法", 《测绘通报》 * |
Cited By (5)
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 |