CN108731648B - 2000 independent coordinate system parameter obtaining method, device and computer readable storage medium - Google Patents
2000 independent coordinate system parameter obtaining method, device and computer readable storage medium Download PDFInfo
- Publication number
- CN108731648B CN108731648B CN201810215665.2A CN201810215665A CN108731648B CN 108731648 B CN108731648 B CN 108731648B CN 201810215665 A CN201810215665 A CN 201810215665A CN 108731648 B CN108731648 B CN 108731648B
- Authority
- CN
- China
- Prior art keywords
- parameter
- projection
- coordinate
- plane
- coordinate system
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C15/00—Surveying instruments or accessories not provided for in groups G01C1/00 - G01C13/00
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Processing Or Creating Images (AREA)
Abstract
The invention provides a 2000 independent coordinate system parameter acquisition method, a device and a computer readable storage medium, wherein the method comprises the steps of setting a first projection parameter according to a central meridian and a projection surface elevation of a target city coordinate system; according to the first projection parameter, carrying out Gaussian projection on CGCS2000 geodetic coordinates of coincident control points in a plurality of control networks of a target city, which are acquired in advance, so as to acquire a first Gaussian projection coordinate; calculating a first plane conversion parameter according to the first Gaussian projection coordinate, the pre-acquired existing plane coordinate of the target city coincidence control point and a preset parameter conversion model; and determining the optimal translation parameter and the optimal projection parameter according to the rotation angle parameter and the scale parameter of the first plane conversion parameter. By the method, the projection parameters and the translation parameters of the target city 2000 independent coordinate system can be accurately obtained, the target city 2000 independent coordinate system and the existing coordinate system are in smooth transition, and the coordinate consistency between the target city 2000 independent coordinate system and the existing coordinate system of the city is improved.
Description
Technical Field
The invention relates to the technical field of coordinate system parameter acquisition, in particular to a 2000 independent coordinate system parameter acquisition method, a device and a computer readable storage medium.
Background
At present, under the requirement of establishing a 2000 independent coordinate system of a city by the national mapping bureau, a 2000 independent coordinate system is established in a part of domestic cities, for example: wuhan 2000 city independent coordinate system, Ningbo 2000 city independent coordinate system and Nanchang 2000 city independent coordinate system; however, there are many problems in establishing 2000 city independent coordinate systems, and most typically, there is no guarantee of consistency between old and new coordinates, for example: the Wuhan 2000 city independent coordinate system has certain difference in coordinates because of the out-of-limit projection deformation of east and west edge areas of the city and areas with the ground height of more than 145m and no consideration of the connection problem of new and old data; the Ningbo 2000 city independent coordinate system only needs translation between new and old coordinates because the projection deformation of east and west sides of the city does not meet the requirement of less than 1/4 ten thousand (2.5cm/km), but a big number A with translation parameters of 100km plus a small number B, but different areas B are inconsistent and are not convenient to use; the independent coordinate system of the Nanchang 2000 city has large coordinate difference between a new coordinate system and an old coordinate system due to the fact that projection deformation of east and west marginal areas does not meet requirements, topographic map data of a planning area 1:500 can be continuously used after coordinate conversion, and a large amount of manpower, material resources and financial resources are required to be invested to develop map conversion work. Therefore, how to realize meeting the standard requirement that the projection length deformation does not exceed 1/4 ten thousands in the process of newly building the city 2000 independent coordinate system based on the CGCS2000 reference ellipsoid and meeting the seamless connection with historical data enables the achievement under the existing coordinate system to be continuously used, thereby avoiding huge conversion workload and cost caused by the inconsistency of the two coordinate systems, and becoming the problem to be urgently solved for newly building the city 2000 independent coordinate system based on the CGCS2000 reference ellipsoid.
Disclosure of Invention
The invention aims to provide a 2000 independent coordinate system parameter acquisition method, a device and a computer readable storage medium, which can accurately acquire projection parameters and translation parameters of a 2000 coordinate system of a target city based on a 2000 ellipsoid in the city, realize the stable transition of the 2000 independent coordinate system of the target city and the existing coordinate system, improve the coordinate consistency between the 2000 independent coordinate system of the target city and the existing coordinate system of the city, enable new and old data to be seamlessly connected, simultaneously enable achievements in the existing coordinate system to be continuously used, and reduce the conversion workload and cost.
The embodiment of the invention provides a 2000 independent coordinate system parameter acquisition method, which comprises the following steps:
the method comprises the steps of obtaining CGCS2000 geodetic coordinates of coincident control points in a plurality of control networks of a target city and plane coordinates of a target city coordinate system corresponding to the coincident control points;
setting a first projection parameter according to a central meridian and a projection surface elevation of the target city coordinate system;
according to the first projection parameter, carrying out Gaussian projection on the CGCS2000 geodetic coordinate of the superposition control point, and calculating a first Gaussian projection coordinate of the superposition control point;
calculating a first plane conversion parameter between the first Gaussian projection coordinate of the coincident control point and the plane coordinate according to the first Gaussian projection coordinate of the coincident control point, the plane coordinate and a preset parameter conversion model;
and determining an optimal translation parameter and an optimal projection parameter according to the rotation angle parameter and the scale parameter of the first plane transformation parameter.
Preferably, the determining an optimal translation parameter and an optimal projection parameter according to a rotation angle parameter and a scale parameter of the first planar transformation parameter specifically includes;
setting a second projection parameter according to the first plane conversion parameter;
according to the second projection parameter, carrying out Gaussian projection on the CGCS2000 geodetic coordinate of the superposition control point, and calculating a second Gaussian projection coordinate of the superposition control point;
calculating a second plane conversion parameter between the second Gaussian projection coordinate of the coincident control point and the plane coordinate of the coincident control point according to the second Gaussian projection coordinate of the coincident control point, the plane coordinate and the preset parameter conversion model;
and determining the optimal translation parameter and the optimal projection parameter according to the rotation angle parameter and the scale parameter of the second plane conversion parameter.
Preferably, the setting of the first projection parameter according to the central meridian and the elevation of the projection plane of the target city coordinate system specifically includes:
the first projection parameters comprise a first central meridian and a first projection surface elevation;
determining a numerical interval of the first central meridian according to a central meridian of the target city coordinate system, and setting the first central meridian according to a set first interval;
and determining a numerical interval of the first projection elevation surface according to the projection elevation surface of the target city coordinate system, and setting the first projection elevation surface according to a set second interval.
Preferably, the setting a second projection parameter according to the first plane conversion parameter specifically includes:
the second projection parameters comprise a second central meridian and a second projection surface elevation;
determining a numerical interval of the second central meridian and a numerical interval of the second projection surface elevation according to the first plane conversion parameter;
the second central meridian is set at a set third interval and the second projected elevation surface is set at a set fourth interval.
Preferably, the setting of the second central meridian at a set third interval and the setting of the second projected elevation surface at a set fourth interval specifically include;
when the second central meridian is L, setting the first projection elevation surface according to the set fourth interval, and obtaining N groups of parameter combinations of the second central meridian and the second projection elevation surface; wherein, L is a central meridian which belongs to the numerical interval of the second central meridian and is arranged according to a set third interval;
when the second projection elevation surface is H, setting the second central meridian according to a set third interval to obtain M groups of parameter combinations of the second central meridian and the second projection elevation surface; h is a projection elevation surface which belongs to the numerical interval of the second projection elevation surface and is arranged according to a set fourth interval.
Preferably, the determining an optimal translation parameter and an optimal projection parameter according to the rotation angle parameter and the scale parameter of the second planar transformation parameter specifically includes:
acquiring an ith second plane conversion parameter which accords with the rotation angle parameter alpha of the second plane conversion parameter and the scale parameter M of the second plane conversion parameter, wherein the rotation angle parameter alpha is approximately equal to 0, and the scale parameter M of the second plane conversion parameter is approximately equal to 1 from the M + N second plane conversion parameters; wherein M + N second plane transformation parameters are obtained according to N sets of parameter combinations of the second central meridian and the second projection elevation surface and M sets of parameter combinations of the second central meridian and the second projection elevation surface;
acquiring a translation parameter of the ith second plane conversion parameter; wherein, the translation parameter of the ith second plane conversion parameter is the optimal translation parameter;
the optimal projection parameters comprise an optimal central meridian and an optimal projection elevation surface;
fitting a first fitting function of the rotation angle parameter and the second central meridian according to the rotation angle parameter alpha of the M + N second plane transformation parameters and the second central meridian L;
calculating the corresponding optimal central meridian when the rotation angle parameter alpha is 0 according to the first fitting function;
fitting a second fitting function of the scale parameters and the second projection elevation surface according to the scale parameters M of the M + N second planar transformation parameters and the second projection elevation surface H;
and calculating the corresponding optimal projection elevation surface when the scale parameter m is 1 according to the second fitting function.
Preferably, the 2000 independent coordinate system parameter obtaining method further includes:
calculating the internal precision of the coincident control point according to the first Gaussian projection coordinate of the coincident control point, the plane coordinate and a preset parameter conversion model;
obtaining external check points of the plurality of control networks; wherein the external check point is a control point in the plurality of control nets that is not coincident with the coincident control point;
calculating the external precision of the coincident control point according to the coincident control point and the external check point;
determining an abnormal control point in the superposition control points according to the internal precision and the external precision of the superposition control points;
and eliminating abnormal control points in the overlapped control points, and recalculating the first plane conversion parameters corresponding to the overlapped control points after the abnormal control points are eliminated.
The embodiment of the present invention further provides a 2000 independent coordinate system parameter obtaining apparatus, including:
the system comprises a coordinate acquisition module, a data processing module and a data processing module, wherein the coordinate acquisition module is used for acquiring CGCS2000 geodetic coordinates of coincident control points in a plurality of control networks of a target city and plane coordinates of a target city coordinate system corresponding to the coincident control points;
the first projection parameter setting module is used for setting first projection parameters according to a central meridian and projection surface elevation of the target city coordinate system;
the first Gaussian projection coordinate calculation module is used for performing Gaussian projection on the CGCS2000 geodetic coordinates of the superposition control points according to the first projection parameters and calculating first Gaussian projection coordinates of the superposition control points;
the first plane conversion parameter calculation module is used for calculating a first plane conversion parameter between the first Gaussian projection coordinate of the superposition control point and the plane coordinate according to the first Gaussian projection coordinate of the superposition control point, the plane coordinate and a preset parameter conversion model;
and the coordinate system parameter acquisition module is used for determining an optimal translation parameter and an optimal projection parameter according to the rotation angle parameter and the scale parameter of the first plane conversion parameter.
The embodiment of the present invention further provides a 2000 independent coordinate system parameter acquiring apparatus, which includes a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, and when the processor executes the computer program, the 2000 independent coordinate system parameter acquiring method is implemented.
An embodiment of the present invention further provides a computer-readable storage medium, where the computer-readable storage medium includes a stored computer program, and when the computer program runs, the apparatus where the computer-readable storage medium is located is controlled to execute the above 2000 independent coordinate system parameter obtaining method.
Compared with the prior art, the method for acquiring the parameters of the 2000 independent coordinate system has the beneficial effects that: the 2000 independent coordinate system parameter acquisition method comprises the following steps: the method comprises the steps of obtaining CGCS2000 geodetic coordinates of coincident control points in a plurality of control networks of a target city and plane coordinates of a target city coordinate system corresponding to the coincident control points; setting a first projection parameter according to a central meridian and a projection surface elevation of the target city coordinate system; according to the first projection parameter, carrying out Gaussian projection on the CGCS2000 geodetic coordinate of the superposition control point, and calculating a first Gaussian projection coordinate of the superposition control point; calculating a first plane conversion parameter between the first Gaussian projection coordinate of the coincident control point and the plane coordinate according to the first Gaussian projection coordinate of the coincident control point, the plane coordinate and a preset parameter conversion model; and determining an optimal translation parameter and an optimal projection parameter according to the rotation angle parameter and the scale parameter of the first plane transformation parameter. The 2000 independent coordinate system parameter acquisition method can accurately acquire the projection parameters of the 2000 independent coordinate system of the target city, thereby improving the coordinate consistency between the 2000 independent coordinate system of the target city and the old coordinate system, enabling new and old data to be seamlessly connected, and simultaneously, achievements under the existing coordinate system can be continuously used, thereby reducing the conversion workload and the cost. The embodiment of the invention also provides a 2000 independent coordinate system parameter acquisition device and a computer readable storage medium.
Drawings
Fig. 1 is a flowchart of a 2000 independent coordinate system parameter obtaining method according to an embodiment of the present invention;
fig. 2 is a schematic diagram of a 2000 independent coordinate system parameter obtaining apparatus according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Please refer to fig. 1, which is a flowchart illustrating a 2000 independent coordinate system parameter obtaining method according to an embodiment of the present invention; the 2000 independent coordinate system parameter acquisition method comprises the following steps:
s100: the method comprises the steps of obtaining CGCS2000 geodetic coordinates of coincident control points in a plurality of control networks of a target city and plane coordinates of a target city coordinate system corresponding to the coincident control points;
s200: setting a first projection parameter according to a central meridian and a projection surface elevation of the target city coordinate system;
s300: according to the first projection parameter, carrying out Gaussian projection on the CGCS2000 geodetic coordinate of the superposition control point, and calculating a first Gaussian projection coordinate of the superposition control point;
s400: calculating a first plane conversion parameter between the first Gaussian projection coordinate of the coincident control point and the plane coordinate according to the first Gaussian projection coordinate of the coincident control point, the plane coordinate and a preset parameter conversion model;
s500: and determining an optimal translation parameter and an optimal projection parameter according to the rotation angle parameter and the scale parameter of the first plane transformation parameter.
In this embodiment, the plurality of control networks include a guangzhou CORS reference station main network, a guangzhou quasi-geoid refinement frame network, a guangzhou quasi-geoid refinement high-precision control network, and a guangzhou GPS primary (secondary) plane measurement control network in 2000 years; the target city coordinate system is a Guangzhou coordinate system, and the CGCS2000 geodetic coordinate is a CGCS2000 geodetic coordinate.
Before the step S100, CGCS2000 geodetic coordinates of control points of a plurality of control networks need to be calculated in advance, including: taking a control point of a province CORS system of a target city as a reference, performing adjustment calculation on a CORS control point main network of the target city to obtain a CGCS2000 geodetic coordinate corresponding to the control point of the CORS control point main network; and respectively carrying out adjustment calculation processing on a first control network, a second control network and a third control network of the target city by taking the CORS control point main network as a reference, and respectively obtaining the CGCS2000 geodetic coordinates corresponding to the control points of the first control network, the second control network and the third control network. The first control network, the second control network and the third control network are respectively the Guangzhou city geoid level refinement frame network, the Guangzhou city geoid level refinement high-precision control network and the Guangzhou city 2000-year GPS first-level (second-level) plane measurement control network; and the main network of the CORS control point is a main network of a Guangzhou CORS reference station.
Specifically, the method comprises the steps of obtaining a CGCS2000 geodetic coordinate corresponding to a control point of an existing CGCS2000 geodetic coordinate of a main network of a Guangzhou CORS reference station, taking a CORS control point of Guangdong province as a starting point according to the CGCS2000 geodetic coordinate of the control point, adopting POWERNET adjustment software of Wuhan university to perform compatibility analysis on the starting point, eliminating incompatible control points, and performing adjustment calculation on CGCS2000 geodetic coordinates of the remaining control points of the main network of the Guangzhou CORS reference station. Then, the control point of the major Guangzhou CORS reference station network is taken as a starting point, the same method is adopted to sequentially carry out adjustment on the Guangzhou similar geoid precision frame network, the Guangzhou similar geoid precision high-precision control network and the Guangzhou 2000-year GPS first-level (second-class) plane measurement control network, and the CGCS2000 geodetic coordinates of the Guangzhou similar geoid precision frame network, the Guangzhou similar geoid precision high-precision control network and the Guangzhou 2000-year GPS first-level (second-class) plane measurement control network are calculated.
In step S300, according to the first projection parameter, performing gaussian projection on the CGCS2000 geodetic coordinates of the coincidence control point under a 2000 reference ellipsoid, and calculating a first gaussian projection coordinate of the coincidence control point;
wherein, the Gaussian projection model is as follows:
wherein, a2The long half shaft is formed after the ellipsoid expands (contracts); e.g. of the type2A square of the ellipsoidal first eccentricity; n is a radical of1The curvature radius of the prime fourth unitary ring before ellipsoid expansion; b is2The latitude after the ellipsoid expansion; b is1The latitude before the expansion of the ellipsoid, B1≠B2；I＝L-L0L is the central meridian; h is projection surface elevation; ρ is a constant; η ═ e' cos B; n is a/W; x represents the arc length of the central meridian; b represents the latitude of a control point。
The preset parameter conversion model is as follows:
wherein x is1，y1The plane coordinates of a target city coordinate system corresponding to the coincidence control points; x is the number of2，y2The Gaussian projection coordinates of the coincident control points are obtained; Δ x, Δ y are translation parameters; alpha is a rotation parameter; and m is a scale parameter.
By the method, the projection parameters and the translation parameters of the 2000 coordinate system of the target city based on the 2000 ellipsoid in the city can be accurately acquired, the stable transition between the 2000 independent coordinate system of the target city and the existing coordinate system is realized, the coordinate consistency between the 2000 independent coordinate system of the target city and the existing coordinate system of the city is improved, new and old data are seamlessly connected, meanwhile, the result under the existing coordinate system can be continuously used, and the conversion workload and cost are reduced.
In an alternative embodiment, S500: determining an optimal translation parameter and an optimal projection parameter according to a rotation angle parameter and a scale parameter of the first plane conversion parameter, wherein the optimal translation parameter and the optimal projection parameter specifically comprise;
setting a second projection parameter according to the first plane conversion parameter;
according to the second projection parameter, carrying out Gaussian projection on the CGCS2000 geodetic coordinate of the superposition control point, and calculating a second Gaussian projection coordinate of the superposition control point;
calculating a second plane conversion parameter between the second Gaussian projection coordinate of the coincident control point and the plane coordinate of the coincident control point according to the second Gaussian projection coordinate of the coincident control point, the plane coordinate and the preset parameter conversion model;
and determining the optimal translation parameter and the optimal projection parameter according to the rotation angle parameter and the scale parameter of the second plane conversion parameter.
In an optional embodiment, the setting a first projection parameter according to the central meridian and the elevation of the projection plane of the target city coordinate system specifically includes:
the first projection parameters comprise a first central meridian and a first projection surface elevation;
determining a numerical interval of the first central meridian according to a central meridian of the target city coordinate system, and setting the first central meridian according to a set first interval;
and determining a numerical interval of the first projection elevation surface according to the projection elevation surface of the target city coordinate system, and setting the first projection elevation surface according to a set second interval.
In an optional embodiment, the setting a second projection parameter according to the first plane transformation parameter specifically includes:
the second projection parameters comprise a second central meridian and a second projection surface elevation;
determining a numerical interval of the second central meridian and a numerical interval of the second projection surface elevation according to the first plane conversion parameter;
the second central meridian is set at a set third interval and the second projected elevation surface is set at a set fourth interval.
In an alternative embodiment, the setting the second central meridian at a set third interval and the setting the second projected elevation surface at a set fourth interval specifically include;
when the second central meridian is L, setting the first projection elevation surface according to the set fourth interval, and obtaining N groups of parameter combinations of the second central meridian and the second projection elevation surface; wherein, L is a central meridian which belongs to the numerical interval of the second central meridian and is arranged according to a set third interval;
when the second projection elevation surface is H, setting the second central meridian according to a set third interval to obtain M groups of parameter combinations of the second central meridian and the second projection elevation surface; h is a projection elevation surface which belongs to the numerical interval of the second projection elevation surface and is arranged according to a set fourth interval.
In an optional embodiment, the determining an optimal translation parameter and an optimal projection parameter according to the rotation angle parameter and the scale parameter of the second planar transformation parameter specifically includes:
acquiring the ith second plane conversion parameter when the rotation angle parameter alpha of the second plane conversion parameter is approximately equal to a first threshold value and the scale parameter M of the second plane conversion parameter is approximately equal to a second threshold value from the M + N second plane conversion parameters; wherein M + N second plane transformation parameters are obtained according to N sets of parameter combinations of the second central meridian and the second projection elevation surface and M sets of parameter combinations of the second central meridian and the second projection elevation surface;
acquiring a translation parameter of the ith second plane conversion parameter; wherein, the translation parameter of the ith second plane conversion parameter is the optimal translation parameter;
the optimal projection parameters comprise an optimal central meridian and an optimal projection elevation surface;
fitting a first fitting function of the rotation angle parameter and the second central meridian according to the rotation angle parameter alpha of the M + N second plane transformation parameters and the second central meridian L;
calculating a corresponding optimal central meridian when the rotation angle parameter alpha is equal to the first threshold value according to the first fitting function;
fitting a second fitting function of the scale parameters and the second projection elevation surface according to the scale parameters M of the M + N second planar transformation parameters and the second projection elevation surface H;
and calculating the corresponding optimal projection elevation surface when the scale parameter m is equal to the second threshold value according to the second fitting function.
Wherein the first threshold is 0 and the second threshold is 1. Calculating first Gaussian projection coordinates which are in one-to-one correspondence with parameter combinations of the first central meridian and the first projection elevation surface through the Gaussian model, then calculating first plane conversion parameters which are in one-to-one correspondence with the first Gaussian projection coordinates through the parameter conversion model, and determining the numerical ranges of the projection elevation surface and the central meridian which are corresponding to the situation that m is approximately equal to 1 and alpha is approximately equal to 0 according to the rotation angle parameter alpha and the scale parameter m in the first plane conversion parameters, namely the numerical range of the second central meridian and the numerical range of the second projection elevation surface. Specifically, the numerical range of the second central meridian and the numerical range of the second projection elevation surface can be obtained according to the occurrence probabilities of m ≈ 1 and α ≈ 0 by calculating the occurrence probabilities of m ≈ 1 and α ≈ 0 in the numerical range of the first central meridian and the numerical range of the first projection elevation surface.
Similarly, second Gaussian projection coordinates which are in one-to-one correspondence with parameter combinations of the second central meridian and the second projection elevation surface are calculated through the Gaussian model, then second plane conversion parameters which are in one-to-one correspondence with the second Gaussian projection coordinates are calculated through the parameter conversion model, and according to a rotation angle parameter alpha and a scale parameter m in the second plane conversion parameters, an ith second plane conversion parameter corresponding to the situation that m is approximately equal to 1 and alpha is approximately equal to 0 is searched, wherein a translation parameter in the ith second plane conversion parameter is the optimal translation parameter.
The first fitting function is: α ═ a1*L-b1(ii) a The second fitting function is: a is2*H+b2(ii) a Wherein, a1、b1Two constants for the first fitting function; a is2、b2Are two constants of the second fitting function.
In an optional embodiment, the 2000 independent coordinate system parameter obtaining method further includes:
calculating the internal precision of the coincident control point according to the first Gaussian projection coordinate of the coincident control point, the plane coordinate and a preset parameter conversion model;
obtaining external check points of the plurality of control networks; wherein the external check point is a control point in the plurality of control nets that is not coincident with the coincident control point;
calculating the external precision of the coincident control point according to the coincident control point and the external check point;
determining an abnormal control point in the superposition control points according to the internal precision and the external precision of the superposition control points;
and eliminating abnormal control points in the overlapped control points, and recalculating the first plane conversion parameters corresponding to the overlapped control points after the abnormal control points are eliminated.
Please refer to fig. 2, which is a schematic diagram of a 2000 independent coordinate system parameter obtaining apparatus according to an embodiment of the present invention, wherein the 2000 independent coordinate system parameter obtaining apparatus includes:
the system comprises a coordinate acquisition module 1, a coordinate processing module and a control module, wherein the coordinate acquisition module is used for acquiring CGCS2000 geodetic coordinates of coincident control points in a plurality of control networks of a target city and plane coordinates of a target city coordinate system corresponding to the coincident control points;
the first projection parameter setting module 2 is used for setting first projection parameters according to a central meridian and a projection surface elevation of the target city coordinate system;
the first Gaussian projection coordinate calculating module 3 is configured to perform Gaussian projection on the CGCS2000 geodetic coordinates of the coincidence control points according to the first projection parameter, and calculate first Gaussian projection coordinates of the coincidence control points;
the first plane conversion parameter calculation module 4 is configured to calculate a first plane conversion parameter between the first gaussian projection coordinate of the coincidence control point and the plane coordinate according to the first gaussian projection coordinate of the coincidence control point, the plane coordinate, and a preset parameter conversion model;
and the coordinate system parameter acquisition module 5 is configured to determine an optimal translation parameter and an optimal projection parameter according to the rotation angle parameter and the scale parameter of the first plane transformation parameter.
In this embodiment, the plurality of control networks include a guangzhou CORS reference station main network, a guangzhou quasi-geoid refinement frame network, a guangzhou quasi-geoid refinement high-precision control network, and a guangzhou GPS primary (secondary) plane measurement control network in 2000 years; the target city coordinate system is a Guangzhou coordinate system, and the CGCS2000 geodetic coordinate is a CGCS2000 geodetic coordinate.
The 2000 independent coordinate system parameter acquisition device further comprises a CGCS2000 geodetic coordinate calculation module, wherein the CGCS2000 geodetic coordinate calculation module is used for performing adjustment calculation on a CORS control point main network of a target city by taking a control point of a province CORS system of the target city as a reference to obtain a CGCS2000 geodetic coordinate corresponding to the control point of the CORS control point main network; and respectively carrying out adjustment calculation processing on a first control network, a second control network and a third control network of the target city by taking the CORS control point main network as a reference, and respectively obtaining the CGCS2000 geodetic coordinates corresponding to the control points of the first control network, the second control network and the third control network. The first control network, the second control network and the third control network are respectively the Guangzhou city geoid level refinement frame network, the Guangzhou city geoid level refinement high-precision control network and the Guangzhou city 2000-year GPS first-level (second-level) plane measurement control network; and the main network of the CORS control point is a main network of a Guangzhou CORS reference station.
Specifically, the method comprises the steps of obtaining a CGCS2000 geodetic coordinate corresponding to a control point of an existing CGCS2000 geodetic coordinate of a main network of a Guangzhou CORS reference station, taking a CORS control point of Guangdong province as a starting point according to the CGCS2000 geodetic coordinate of the control point, adopting POWERNET adjustment software of Wuhan university to perform compatibility analysis on the starting point, eliminating incompatible control points, and performing adjustment calculation on CGCS2000 geodetic coordinates of the remaining control points of the main network of the Guangzhou CORS reference station. Then, the control point of the major Guangzhou CORS reference station network is taken as a starting point, the same method is adopted to sequentially carry out adjustment on the Guangzhou similar geoid precision frame network, the Guangzhou similar geoid precision high-precision control network and the Guangzhou 2000-year GPS first-level (second-class) plane measurement control network, and the CGCS2000 geodetic coordinates of the Guangzhou similar geoid precision frame network, the Guangzhou similar geoid precision high-precision control network and the Guangzhou 2000-year GPS first-level (second-class) plane measurement control network are calculated.
The first gaussian projection coordinate calculation module 3 is configured to perform gaussian projection on the CGCS2000 geodetic coordinate of the coincidence control point under a 2000 reference ellipsoid according to the first projection parameter, and calculate a first gaussian projection coordinate of the coincidence control point;
wherein, the Gaussian projection model is as follows:
wherein, a2The long half shaft is formed after the ellipsoid expands (contracts); e.g. of the type2A square of the ellipsoidal first eccentricity; n is a radical of1The curvature radius of the prime fourth unitary ring before ellipsoid expansion; b is2The latitude after the ellipsoid expansion; b is1The latitude before the expansion of the ellipsoid, B1≠B2；I＝L-L0L is the central meridian; h is projection surface elevation; ρ is a constant; η ═ e' cosB; n is a/W; x represents the arc length of the central meridian; b represents the latitude of a control point.
The preset parameter conversion model is as follows:
wherein x is1，y1The plane coordinates of a target city coordinate system corresponding to the coincidence control points; x is the number of2，y2The Gaussian projection coordinates of the coincident control points are obtained; Δ x, Δ y are translation parameters; alpha is a rotation parameter; and m is a scale parameter.
The projection parameters and the translation parameters of the 2000 coordinate system of the target city based on the 2000 ellipsoid in the city can be accurately acquired through the device, the stable transition between the 2000 independent coordinate system of the target city and the existing coordinate system is realized, the coordinate consistency between the 2000 independent coordinate system of the target city and the existing coordinate system of the city is improved, new and old data are in seamless connection, meanwhile, the achievement can be continuously used under the existing coordinate system, and the conversion workload and cost are reduced.
In an alternative embodiment, the coordinate system parameter obtaining module 5 includes;
the second projection parameter setting module is used for setting second projection parameters according to the first plane conversion parameters;
the second Gaussian projection coordinate calculation module is used for performing Gaussian projection on the CGCS2000 geodetic coordinates of the superposition control points according to the second projection parameters and calculating second Gaussian projection coordinates of the superposition control points;
the second plane conversion parameter calculation module is used for calculating a second plane conversion parameter between the second Gaussian projection coordinate of the superposition control point and the plane coordinate of the superposition control point according to the second Gaussian projection coordinate of the superposition control point, the plane coordinate and the preset parameter conversion model;
and the coordinate system parameter determining module is used for determining the optimal translation parameter and the optimal projection parameter according to the rotation angle parameter and the scale parameter of the second plane transformation parameter.
In an alternative embodiment, the first projection parameters include a first central meridian and a first projection surface elevation;
the first projection parameter setting module is further configured to determine a numerical interval of the first central meridian according to the central meridian of the target city coordinate system, and set the first central meridian according to a set first interval;
the first projection parameter setting module is further configured to determine a numerical interval of the first projection elevation surface according to the projection elevation surface of the target city coordinate system, and set the first projection elevation surface according to a set second interval.
In an alternative embodiment, the second projection parameters include a second central meridian and a second projection surface elevation;
the second projection parameter setting module is further configured to determine a numerical interval of the second central meridian and a numerical interval of the second projection surface elevation according to the first plane conversion parameter;
the second projection parameter setting module is further configured to set the second central meridian according to a set third interval and set the second projection elevation surface according to a set fourth interval.
In an optional embodiment, the second projection parameter setting module is further configured to set the first projection elevation surface according to the set fourth interval when the second central meridian is L, so as to obtain N sets of parameter combinations of the second central meridian and the second projection elevation surface; wherein, L is a central meridian which belongs to the numerical interval of the second central meridian and is arranged according to a set third interval;
the second projection parameter setting module is further configured to set the second central meridian according to a set third interval when the second projection elevation surface is H, so as to obtain M sets of parameter combinations of the second central meridian and the second projection elevation surface; h is a projection elevation surface which belongs to the numerical interval of the second projection elevation surface and is arranged according to a set fourth interval.
In an alternative embodiment, the coordinate system parameter determination module comprises:
a second planar conversion parameter obtaining unit, configured to obtain, from M + N second planar conversion parameters, an ith second planar conversion parameter when a rotation angle parameter α of the second planar conversion parameter is approximately equal to a first threshold and a scale parameter M of the second planar conversion parameter is approximately equal to a second threshold; wherein M + N second plane transformation parameters are obtained according to N sets of parameter combinations of the second central meridian and the second projection elevation surface and M sets of parameter combinations of the second central meridian and the second projection elevation surface;
a translation parameter obtaining unit, configured to obtain a translation parameter of the ith second plane conversion parameter; wherein, the translation parameter of the ith second plane conversion parameter is the optimal translation parameter;
the optimal projection parameters comprise an optimal central meridian and an optimal projection elevation surface;
a first function fitting unit, configured to fit a first fitting function of the rotation angle parameter and the second central meridian according to the rotation angle parameter α and the second central meridian L of the M + N second plane transformation parameters;
a central meridian acquisition unit for calculating, according to the first fitting function, an optimal central meridian corresponding to when a rotation angle parameter α is equal to the first threshold;
a second function fitting unit, configured to fit a second fitting function of the scale parameter and the second projection elevation surface according to the scale parameter M of the M + N second plane transformation parameters and the second projection elevation surface H;
and the projection elevation surface acquisition unit is used for calculating the corresponding optimal projection elevation surface when the scale parameter m is equal to the second threshold according to the second fitting function.
Wherein the first threshold is 0 and the second threshold is 1. Calculating first Gaussian projection coordinates which are in one-to-one correspondence with parameter combinations of the first central meridian and the first projection elevation surface through the Gaussian model, then calculating first plane conversion parameters which are in one-to-one correspondence with the first Gaussian projection coordinates through the parameter conversion model, and determining the numerical ranges of the projection elevation surface and the central meridian which are corresponding to the situation that m is approximately equal to 1 and alpha is approximately equal to 0 according to the rotation angle parameter alpha and the scale parameter m in the first plane conversion parameters, namely the numerical range of the second central meridian and the numerical range of the second projection elevation surface. Specifically, the numerical range of the second central meridian and the numerical range of the second projection elevation surface can be obtained according to the occurrence probabilities of m ≈ 1 and α ≈ 0 by calculating the occurrence probabilities of m ≈ 1 and α ≈ 0 in the numerical range of the first central meridian and the numerical range of the first projection elevation surface.
Similarly, second Gaussian projection coordinates which are in one-to-one correspondence with parameter combinations of the second central meridian and the second projection elevation surface are calculated through the Gaussian model, then second plane conversion parameters which are in one-to-one correspondence with the second Gaussian projection coordinates are calculated through the parameter conversion model, and according to a rotation angle parameter alpha and a scale parameter m in the second plane conversion parameters, an ith second plane conversion parameter corresponding to the situation that m is approximately equal to 1 and alpha is approximately equal to 0 is searched, wherein a translation parameter in the ith second plane conversion parameter is the optimal translation parameter.
The first fitting function is: α ═ a1*L-b1(ii) a The second fitting function is: a is2*H+b2(ii) a Wherein, a1、b1Two constants for the first fitting function; a is2、b2Are two constants of the second fitting function.
In an alternative embodiment, the 2000 independent coordinate system parameter acquiring apparatus further comprises:
the internal precision calculation module is used for calculating the internal precision of the superposition control point according to the first Gaussian projection coordinate of the superposition control point, the plane coordinate and a preset parameter conversion model;
the external check point acquisition module is used for acquiring external check points of the plurality of control networks; wherein the external check point is a control point in the plurality of control nets that is not coincident with the coincident control point;
the external precision calculation unit is used for calculating the external precision of the superposition control point according to the superposition control point and the external check point;
the abnormal control point confirming module is used for determining an abnormal control point in the superposition control points according to the internal precision and the external precision of the superposition control points;
and the abnormal control point removing module is used for removing the abnormal control points in the overlapped control points and recalculating the first plane conversion parameters corresponding to the overlapped control points from which the abnormal control points are removed.
The embodiment of the present invention further provides a 2000 independent coordinate system parameter acquiring apparatus, which includes a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, and when the processor executes the computer program, the 2000 independent coordinate system parameter acquiring method is implemented.
Illustratively, the computer program may be partitioned into one or more modules/units that are stored in the memory and executed by the processor to implement the invention. The one or more modules/units may be a series of computer program instruction segments capable of performing specific functions, which are used for describing the execution process of the computer program in the 2000 independent coordinate system parameter acquisition device. For example, the computer program may be divided into a coordinate acquisition module for acquiring CGCS2000 geodetic coordinates of coincident control points in a plurality of control nets of a target city and planar coordinates of a target city coordinate system corresponding to the coincident control points; the first projection parameter setting module is used for setting first projection parameters according to a central meridian and projection surface elevation of the target city coordinate system; the first Gaussian projection coordinate calculation module is used for performing Gaussian projection on the CGCS2000 geodetic coordinates of the superposition control points according to the first projection parameters and calculating first Gaussian projection coordinates of the superposition control points; the first plane conversion parameter calculation module is used for calculating a first plane conversion parameter between the first Gaussian projection coordinate of the superposition control point and the plane coordinate according to the first Gaussian projection coordinate of the superposition control point, the plane coordinate and a preset parameter conversion model; and the coordinate system parameter acquisition module is used for determining an optimal translation parameter and an optimal projection parameter according to the rotation angle parameter and the scale parameter of the first plane conversion parameter.
The 2000 independent coordinate system parameter acquiring device can be a desktop computer, a notebook, a palm computer, a cloud server and other computing equipment. The 2000 independent coordinate system parameter acquisition device may include, but is not limited to, a processor, a memory. Those skilled in the art will appreciate that the schematic diagram is merely an example of the 2000-independent coordinate system parameter acquiring apparatus, and does not constitute a limitation to the 2000-independent coordinate system parameter acquiring apparatus, and may include more or less components than those shown in the drawings, or combine some components, or different components, for example, the 2000-independent coordinate system parameter acquiring apparatus may further include an input-output device, a network access device, a bus, and the like.
The Processor may be a Central Processing Unit (CPU), other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), an off-the-shelf Programmable Gate Array (FPGA) or other Programmable logic device, discrete Gate or transistor logic, discrete hardware components, etc. The general purpose processor may be a microprocessor or the processor may be any conventional processor or the like, the processor is a control center of the 2000 independent coordinate system parameter acquisition device, and various interfaces and lines are used to connect various parts of the entire 2000 independent coordinate system parameter acquisition device.
The memory may be used to store the computer programs and/or modules, and the processor may implement the various functions of the 2000 independent coordinate system parameter acquisition apparatus by executing or executing the computer programs and/or modules stored in the memory and calling the data stored in the memory. The memory may mainly include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required by at least one function (such as a sound playing function, an image playing function, etc.), and the like; the storage data area may store data (such as audio data, a phonebook, etc.) created according to the use of the cellular phone, and the like. In addition, the memory may include high speed random access memory, and may also include non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), at least one magnetic disk storage device, a Flash memory device, or other volatile solid state storage device.
Wherein, the module/unit integrated with the 2000 independent coordinate system parameter acquiring device can be stored in a computer readable storage medium if it is implemented in the form of software functional unit and sold or used as an independent product. Based on such understanding, all or part of the flow of the method according to the embodiments of the present invention may also be implemented by a computer program, which may be stored in a computer-readable storage medium, and when the computer program is executed by a processor, the steps of the method embodiments may be implemented. Wherein the computer program comprises computer program code, which may be in the form of source code, object code, an executable file or some intermediate form, etc. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, usb disk, removable hard disk, magnetic disk, optical disk, computer Memory, Read-Only Memory (ROM), Random Access Memory (RAM), electrical carrier wave signals, telecommunications signals, software distribution medium, and the like. It should be noted that the computer readable medium may contain content that is subject to appropriate increase or decrease as required by legislation and patent practice in jurisdictions, for example, in some jurisdictions, computer readable media does not include electrical carrier signals and telecommunications signals as is required by legislation and patent practice.
An embodiment of the present invention further provides a computer-readable storage medium, where the computer-readable storage medium includes a stored computer program, and when the computer program runs, the apparatus where the computer-readable storage medium is located is controlled to execute the above 2000 independent coordinate system parameter obtaining method.
Compared with the prior art, the method for acquiring the parameters of the 2000 independent coordinate system has the beneficial effects that: the 2000 independent coordinate system parameter acquisition method comprises the following steps: the method comprises the steps of obtaining CGCS2000 geodetic coordinates of coincident control points in a plurality of control networks of a target city and plane coordinates of a target city coordinate system corresponding to the coincident control points; setting a first projection parameter according to a central meridian and a projection surface elevation of the target city coordinate system; according to the first projection parameter, carrying out Gaussian projection on the CGCS2000 geodetic coordinate of the superposition control point, and calculating a first Gaussian projection coordinate of the superposition control point; calculating a first plane conversion parameter between the first Gaussian projection coordinate of the coincident control point and the plane coordinate according to the first Gaussian projection coordinate of the coincident control point, the plane coordinate and a preset parameter conversion model; and determining an optimal translation parameter and an optimal projection parameter according to the rotation angle parameter and the scale parameter of the first plane transformation parameter. The 2000 independent coordinate system parameter acquisition method can accurately acquire the projection parameters of the 2000 independent coordinate system of the target city, thereby improving the coordinate consistency between the 2000 independent coordinate system of the target city and the old coordinate system, enabling new and old data to be seamlessly connected, and simultaneously, achievements under the existing coordinate system can be continuously used, thereby reducing the conversion workload and the cost. The embodiment of the invention also provides a 2000 independent coordinate system parameter acquisition device and a computer readable storage medium.
While the foregoing is directed to the preferred embodiment of the present invention, it will be understood by those skilled in the art that various changes and modifications may be made without departing from the spirit and scope of the invention.
Claims (7)
1. A 2000 independent coordinate system parameter acquisition method, comprising:
the method comprises the steps of obtaining CGCS2000 geodetic coordinates of coincident control points in a plurality of control networks of a target city and plane coordinates of a target city coordinate system corresponding to the coincident control points;
setting a first projection parameter according to a central meridian and a projection surface elevation of the target city coordinate system;
according to the first projection parameter, carrying out Gaussian projection on the CGCS2000 geodetic coordinate of the superposition control point, and calculating a first Gaussian projection coordinate of the superposition control point;
calculating a first plane conversion parameter between the first Gaussian projection coordinate of the coincident control point and the plane coordinate according to the first Gaussian projection coordinate of the coincident control point, the plane coordinate and a preset parameter conversion model;
determining an optimal translation parameter and an optimal projection parameter according to the rotation angle parameter and the scale parameter of the first plane transformation parameter, specifically comprising:
setting a second projection parameter according to the first plane conversion parameter;
according to the second projection parameter, carrying out Gaussian projection on the CGCS2000 geodetic coordinate of the superposition control point, and calculating a second Gaussian projection coordinate of the superposition control point;
calculating a second plane conversion parameter between the second Gaussian projection coordinate of the coincident control point and the plane coordinate of the coincident control point according to the second Gaussian projection coordinate of the coincident control point, the plane coordinate and the preset parameter conversion model;
determining the optimal translation parameter and the optimal projection parameter according to a rotation angle parameter and a scale parameter of the second plane conversion parameter;
wherein, according to the first plane conversion parameter, setting a second projection parameter specifically comprises:
the second projection parameters comprise a second central meridian and a second projection surface elevation;
determining a numerical interval of the second central meridian and a numerical interval of the second projection surface elevation according to the first plane conversion parameter;
setting the second central meridian according to a set third interval and setting the second projected elevation surface according to a set fourth interval, specifically including: when the second central meridian is L, setting the second projection elevation surface according to the set fourth interval, and obtaining N groups of parameter combinations of the second central meridian and the second projection elevation surface; wherein, L is a central meridian which belongs to the numerical interval of the second central meridian and is arranged according to a set third interval;
when the second projection elevation surface is H, setting the second central meridian according to a set third interval to obtain M groups of parameter combinations of the second central meridian and the second projection elevation surface; h is a projection elevation surface which belongs to the numerical interval of the second projection elevation surface and is arranged according to a set fourth interval.
2. The method of claim 1, wherein the setting of the first projection parameters according to the central meridian and the elevation of the projection plane of the target city coordinate system comprises:
the first projection parameters comprise a first central meridian and a first projection surface elevation;
determining a numerical interval of the first central meridian according to a central meridian of the target city coordinate system, and setting the first central meridian according to a set first interval;
and determining a numerical interval of the first projection elevation surface according to the projection elevation surface of the target city coordinate system, and setting the first projection elevation surface according to a set second interval.
3. The method of 2000 independent coordinate system parameter acquisition as claimed in claim 1, wherein the determining an optimal translation parameter and an optimal projection parameter according to a rotation angle parameter and a scale parameter of the second planar transformation parameter specifically comprises:
acquiring an ith second plane conversion parameter which accords with a rotation angle parameter alpha of the second plane conversion parameter and a scale parameter M of the second plane conversion parameter, wherein the rotation angle parameter alpha is approximately equal to a first threshold value, and the scale parameter M of the second plane conversion parameter is approximately equal to a second threshold value from the M + N second plane conversion parameters;
wherein M + N second plane transformation parameters are obtained according to N sets of parameter combinations of the second central meridian and the second projection elevation surface and M sets of parameter combinations of the second central meridian and the second projection elevation surface;
acquiring a translation parameter of the ith second plane conversion parameter; wherein, the translation parameter of the ith second plane conversion parameter is the optimal translation parameter;
the optimal projection parameters comprise an optimal central meridian and an optimal projection elevation surface;
fitting a first fitting function of the rotation angle parameter and the second central meridian according to the rotation angle parameter alpha of the M + N second plane transformation parameters and the second central meridian L;
calculating a corresponding optimal central meridian when the rotation angle parameter alpha is equal to the first threshold value according to the first fitting function;
fitting a second fitting function of the scale parameters and the second projection elevation surface according to the scale parameters M of the M + N second planar transformation parameters and the second projection elevation surface H;
and calculating the corresponding optimal projection elevation surface when the scale parameter m is equal to a second threshold value according to the second fitting function.
4. The 2000 independent coordinate system parameter acquisition method of claim 1, wherein the 2000 independent coordinate system parameter acquisition method further comprises:
calculating the internal precision of the coincident control point according to the first Gaussian projection coordinate of the coincident control point, the plane coordinate and a preset parameter conversion model;
obtaining external check points of the plurality of control networks; wherein the external check point is a control point in the plurality of control nets that is not coincident with the coincident control point;
calculating the external precision of the coincident control point according to the coincident control point and the external check point;
determining an abnormal control point in the superposition control points according to the internal precision and the external precision of the superposition control points;
and eliminating abnormal control points in the overlapped control points, and recalculating the first plane conversion parameters corresponding to the overlapped control points after the abnormal control points are eliminated.
5. A 2000 independent coordinate system parameter acquisition apparatus, comprising:
the system comprises a coordinate acquisition module, a data processing module and a data processing module, wherein the coordinate acquisition module is used for acquiring CGCS2000 geodetic coordinates of coincident control points in a plurality of control networks of a target city and plane coordinates of a target city coordinate system corresponding to the coincident control points;
the first projection parameter setting module is used for setting first projection parameters according to a central meridian and projection surface elevation of the target city coordinate system;
the first Gaussian projection coordinate calculation module is used for performing Gaussian projection on the CGCS2000 geodetic coordinates of the superposition control points according to the first projection parameters and calculating first Gaussian projection coordinates of the superposition control points;
the first plane conversion parameter calculation module is used for calculating a first plane conversion parameter between the first Gaussian projection coordinate of the superposition control point and the plane coordinate according to the first Gaussian projection coordinate of the superposition control point, the plane coordinate and a preset parameter conversion model;
a coordinate system parameter obtaining module, configured to determine an optimal translation parameter and an optimal projection parameter according to a rotation angle parameter and a scale parameter of the first plane transformation parameter;
more specifically, the method is used for setting a second projection parameter according to the first plane conversion parameter;
according to the second projection parameter, carrying out Gaussian projection on the CGCS2000 geodetic coordinate of the superposition control point, and calculating a second Gaussian projection coordinate of the superposition control point;
calculating a second plane conversion parameter between the second Gaussian projection coordinate of the coincident control point and the plane coordinate of the coincident control point according to the second Gaussian projection coordinate of the coincident control point, the plane coordinate and the preset parameter conversion model;
determining the optimal translation parameter and the optimal projection parameter according to a rotation angle parameter and a scale parameter of the second plane conversion parameter;
wherein, according to the first plane conversion parameter, setting a second projection parameter specifically comprises:
the second projection parameters comprise a second central meridian and a second projection surface elevation;
determining a numerical interval of the second central meridian and a numerical interval of the second projection surface elevation according to the first plane conversion parameter;
setting the second central meridian according to a set third interval and setting the second projected elevation surface according to a set fourth interval, specifically including: when the second central meridian is L, setting the second projection elevation surface according to the set fourth interval, and obtaining N groups of parameter combinations of the second central meridian and the second projection elevation surface; wherein, L is a central meridian which belongs to the numerical interval of the second central meridian and is arranged according to a set third interval;
when the second projection elevation surface is H, setting the second central meridian according to a set third interval to obtain M groups of parameter combinations of the second central meridian and the second projection elevation surface; h is a projection elevation surface which belongs to the numerical interval of the second projection elevation surface and is arranged according to a set fourth interval.
6. A 2000 independent coordinate system parameter acquisition apparatus comprising a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, the processor implementing the 2000 independent coordinate system parameter acquisition method of any one of claims 1 to 4 when executing the computer program.
7. A computer-readable storage medium, comprising a stored computer program, wherein the computer program, when executed, controls an apparatus in which the computer-readable storage medium is located to perform the 2000 independent coordinate system parameter acquisition method according to any one of claims 1 to 4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810215665.2A CN108731648B (en) | 2018-03-15 | 2018-03-15 | 2000 independent coordinate system parameter obtaining method, device and computer readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810215665.2A CN108731648B (en) | 2018-03-15 | 2018-03-15 | 2000 independent coordinate system parameter obtaining method, device and computer readable storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108731648A CN108731648A (en) | 2018-11-02 |
| CN108731648B true CN108731648B (en) | 2020-12-22 |
Family
ID=63940492
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810215665.2A Active CN108731648B (en) | 2018-03-15 | 2018-03-15 | 2000 independent coordinate system parameter obtaining method, device and computer readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108731648B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114001650B (en) * | 2021-09-16 | 2023-09-29 | 北京市测绘设计研究院 | Encryption method for conversion parameters of local coordinate system and arbitrary plane coordinate system |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3595451B2 (en) * | 1998-07-24 | 2004-12-02 | ペンタックス株式会社 | Laser surveying equipment |
| CN103389080A (en) * | 2013-08-19 | 2013-11-13 | 重庆市地理信息中心 | Method for obtaining independent city coordinate system parameter based on geographic information application |
| CN103413352A (en) * | 2013-07-29 | 2013-11-27 | 西北工业大学 | Scene three-dimensional reconstruction method based on RGBD multi-sensor fusion |
| CN103827631A (en) * | 2011-09-27 | 2014-05-28 | 莱卡地球系统公开股份有限公司 | Measuring system and method for marking a known target point in a coordinate system |
| CN104123695A (en) * | 2014-06-23 | 2014-10-29 | 唐文兴 | Method for realizing coordinate conversion |
| CN104567840A (en) * | 2015-01-27 | 2015-04-29 | 国家测绘地理信息局大地测量数据处理中心 | New and old urban independent coordinate system fusion method |
| CN104598746A (en) * | 2015-01-27 | 2015-05-06 | 国家测绘地理信息局大地测量数据处理中心 | Applicability discrimination method for coordinate transformation model |
| CN104766302A (en) * | 2015-02-05 | 2015-07-08 | 武汉大势智慧科技有限公司 | Method and system for optimizing laser scanning point cloud data by means of unmanned aerial vehicle images |
| CN105388507A (en) * | 2015-11-24 | 2016-03-09 | 中国电建集团西北勘测设计研究院有限公司 | Method for determining ellipsoid parameters based on area GNSS and precise distance measurement scale ratio estimation method |
| CN106202000A (en) * | 2016-07-22 | 2016-12-07 | 武汉大学 | Seven-parameter transformation method between country's three-dimensional system of coordinate and anywhere plane coordinate system |
| CN106289195A (en) * | 2016-08-29 | 2017-01-04 | 长江空间信息技术工程有限公司(武汉) | The method for building up of plateau distance heavy construction Measurement and Control System |
| CN107167119A (en) * | 2017-03-29 | 2017-09-15 | 中交烟台环保疏浚有限公司 | The data processing method of distortion of projection |
-
2018
- 2018-03-15 CN CN201810215665.2A patent/CN108731648B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3595451B2 (en) * | 1998-07-24 | 2004-12-02 | ペンタックス株式会社 | Laser surveying equipment |
| CN103827631A (en) * | 2011-09-27 | 2014-05-28 | 莱卡地球系统公开股份有限公司 | Measuring system and method for marking a known target point in a coordinate system |
| CN103413352A (en) * | 2013-07-29 | 2013-11-27 | 西北工业大学 | Scene three-dimensional reconstruction method based on RGBD multi-sensor fusion |
| CN103389080A (en) * | 2013-08-19 | 2013-11-13 | 重庆市地理信息中心 | Method for obtaining independent city coordinate system parameter based on geographic information application |
| CN104123695A (en) * | 2014-06-23 | 2014-10-29 | 唐文兴 | Method for realizing coordinate conversion |
| CN104567840A (en) * | 2015-01-27 | 2015-04-29 | 国家测绘地理信息局大地测量数据处理中心 | New and old urban independent coordinate system fusion method |
| CN104598746A (en) * | 2015-01-27 | 2015-05-06 | 国家测绘地理信息局大地测量数据处理中心 | Applicability discrimination method for coordinate transformation model |
| CN104766302A (en) * | 2015-02-05 | 2015-07-08 | 武汉大势智慧科技有限公司 | Method and system for optimizing laser scanning point cloud data by means of unmanned aerial vehicle images |
| CN105388507A (en) * | 2015-11-24 | 2016-03-09 | 中国电建集团西北勘测设计研究院有限公司 | Method for determining ellipsoid parameters based on area GNSS and precise distance measurement scale ratio estimation method |
| CN106202000A (en) * | 2016-07-22 | 2016-12-07 | 武汉大学 | Seven-parameter transformation method between country's three-dimensional system of coordinate and anywhere plane coordinate system |
| CN106289195A (en) * | 2016-08-29 | 2017-01-04 | 长江空间信息技术工程有限公司(武汉) | The method for building up of plateau distance heavy construction Measurement and Control System |
| CN107167119A (en) * | 2017-03-29 | 2017-09-15 | 中交烟台环保疏浚有限公司 | The data processing method of distortion of projection |
Non-Patent Citations (4)
| Title |
|---|
| Combined Adjustment Project of National Astronomical Geodetic Networks and 2000 National GPS Control Network;YANG Y 等;《Progress in National Science》;20151231;第15卷(第4期);第435-441页 * |
| 参心独立坐标系向地心独立坐标系转换地块面积变化分析和计算方法研究;王维 等;《测绘与空间地理信息》;20130228;第36卷(第2期);第145-147页 * |
| 地方独立坐标系向2000国家大地坐标系转换研究;耿晓燕;《中国优秀硕士学位论文全文数据库基础科学辑》;20110515(第5期);第A008-41页 * |
| 宁波市2000坐标系建立及推广应用研究;蔡元波 等;《城市勘探》;20140430(第2期);第119-121页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108731648A (en) | 2018-11-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR102068419B1 (en) | Method, apparatus and computer readable medium for adjusting point cloud data collection trajectory | |
| CN109269472B (en) | Method and device for extracting characteristic line of oblique photogrammetry building and storage medium | |
| CN109186551B (en) | Method and device for extracting characteristic points of oblique photogrammetry building and storage medium | |
| WO2020107326A1 (en) | Lane line detection method, device and computer readale storage medium | |
| CN112233240A (en) | Three-dimensional vector data slicing method and device of three-dimensional vector map and electronic equipment | |
| CN108731649B (en) | 2000 mapping reference frame unifying method, device and computer readable storage medium | |
| WO2020233083A1 (en) | Image restoration method and apparatus, storage medium, and terminal device | |
| WO2022110862A1 (en) | Method and apparatus for constructing road direction arrow, electronic device, and storage medium | |
| CN108731648B (en) | 2000 independent coordinate system parameter obtaining method, device and computer readable storage medium | |
| CN113709006B (en) | Flow determination method and device, storage medium and electronic device | |
| CN110910446A (en) | Method and device for determining building removal area and method and device for determining indoor area of building | |
| JP2024511018A (en) | Methods, devices, computer devices and storage media for determining spatial relationships | |
| CN114627206A (en) | Grid drawing method and device, electronic equipment and computer readable storage medium | |
| CN109341704B (en) | Map precision determination method and device | |
| CN111723607A (en) | Antenna engineering parameter determination method and device | |
| CN113139218B (en) | Method and device for drawing outer diameter side line of shield segment and computer equipment | |
| CN116086411A (en) | Digital topography generation method, device, equipment and readable storage medium | |
| CN110489510B (en) | Road data processing method and device, readable storage medium and computer equipment | |
| CN109064393B (en) | Face special effect processing method and device | |
| CN113015117A (en) | User positioning method and device, electronic equipment and storage medium | |
| CN110704897A (en) | Method for placing connecting node between wall keel model and bottom guide beam model and product | |
| CN111191395A (en) | Nested model modeling method and equipment | |
| CN109813255A (en) | The area computation method of big trapezoidal segment on a kind of earth ellipsoid face | |
| CN113990066B (en) | Road information matching and intersection identification method and device | |
| CN111210500B (en) | Three-dimensional point cloud processing method and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |