Engineering ellipsoid establishing method and device suitable for construction measurement control network

Abstract

The invention relates to the technical field of construction measurement, in particular to an engineering ellipsoid establishing method and device suitable for a construction measurement control network. The method comprises the following steps: establishing a reference ellipsoid, and obtaining geodetic coordinates and normal heights of a position reference point, a plane and an elevation control point of an engineering area; changing the semi-major axis and eccentricity of the reference ellipsoid based on geodetic coordinates and normal height, and establishing a new ellipsoid; converting the geodetic coordinates of the position reference point, the plane and the elevation control point into geodetic coordinates under a new ellipsoid; determining an inclination angle between the ellipsoid surface of the new ellipsoid engineering area and the engineering projection surface based on the deviation of the geodetic coordinate variation of the plane and the elevation control point when the new ellipsoid rotates; and rotating the plane and the elevation control point based on the inclination angle, and converting the geodetic coordinates of the rotated plane and elevation control point under the new ellipsoid to obtain the coordinates under the engineering independent coordinate system, thereby completing the construction of the engineering ellipsoid. The method solves the problem that the reference ellipsoid and the projection plane of the engineering area cannot coincide well.

Description

Technical Field

[0001] The invention relates to the technical field of construction measurement, and in particular to a method and a device for establishing an engineering ellipsoid suitable for a construction measurement control network. Background Art

[0002] The construction survey control network is a survey control network set up to meet the needs of engineering construction. It is the starting data for construction layout of the project area, project completion, building settlement observation, and subsequent building reconstruction and expansion. The accuracy requirements of the construction control network are determined according to the nature of the project construction. In order to meet the accuracy of the layout of the construction site and the coordination of the overall layout of the project, a construction survey control network must be established before construction.

[0003] The construction survey control network is divided into a plane control network and an elevation control network. The plane control network is mainly constructed using methods such as the corner network and GNSS (Global Navigation Satellite System) measurement. Among them, the GNSS measurement method has the advantages of flexible control network construction and convenient observation, but its calculation and adjustment are performed on the reference ellipsoid, and the construction survey control network usually selects a specified elevation surface as the reference surface for side length reduction and projection. When the two are inconsistent, the coordinate inverse calculated side length of the GNSS control network will be different from the actual side length, especially for areas with a large inclination angle between the reference ellipsoid and the reference surface, which will produce large errors. Summary of the invention

[0004] In view of this, the present invention provides a method and device for establishing an engineering ellipsoid suitable for a construction survey control network, so as to solve the technical problem that the use of a GNSS measurement method to construct a plane control network easily leads to errors in the back-calculated side lengths.

[0005] In the first aspect, the present invention provides a method for establishing an engineering ellipsoid suitable for a construction measurement control network, the method comprising: establishing a reference ellipsoid, obtaining the geodetic coordinates and normal heights of the position reference points, plane and elevation control points of the engineering area; changing the major semi-axis and eccentricity of the reference ellipsoid based on the geodetic coordinates and normal height to establish a new ellipsoid; converting the geodetic coordinates of the position reference points, plane and elevation control points into geodetic coordinates under the new ellipsoid; determining the inclination angle between the ellipsoidal surface of the new ellipsoid engineering area and the engineering projection surface based on the deviation of the change in the geodetic coordinates of the plane and elevation control points when the new ellipsoid is rotated; rotating the plane and elevation control points based on the inclination angle, and converting the geodetic coordinates of the rotated plane and elevation control points under the new ellipsoid to obtain the coordinates under the engineering independent coordinate system, thereby completing the establishment of the engineering ellipsoid.

[0006] The method for establishing an engineering ellipsoid applicable to a construction measurement control network provided by an embodiment of the present invention adjusts the parameters of a reference ellipsoid of an engineering area so that it is as close to the engineering projection surface as possible, and then further rotates the ellipsoid to eliminate the deviation between the normal direction of the ellipsoid surface and the vertical direction of the engineering projection surface, and based on the adjustment and rotation of the ellipsoid parameters, the coordinates of each control point are converted to obtain the coordinates in the independent engineering coordinate system, and the engineering ellipsoid is established. In this way, the inclination angle between the ellipsoid surface of the engineering area and the engineering projection surface can be eliminated, so that the inverse distance of the projection coordinates under the new ellipsoid is consistent with the actual electromagnetic wave ranging, and it is easy to use.

[0007] In an optional implementation, the geodetic coordinates and normal heights of the position reference points, plane and elevation control points of the project area are obtained, including: determining the geodetic longitude, geodetic latitude, geodetic height and normal height of the position reference points based on the project area, and determining the normal height of the project projection surface; using the GNSS method to calculate the geodetic longitude, geodetic latitude and geodetic height of each plane and elevation control point; and using the leveling method to calculate the leveling height of each plane and elevation control point.

[0008] In this embodiment, by calculating the geodetic longitude, geodetic latitude, geodetic height and normal height of the position reference point, the normal height of the projection surface, and the geodetic longitude, geodetic latitude, geodetic height and normal height of each plane and elevation control point under the reference ellipsoid, a data basis is provided for subsequent coordinate transformation and ellipsoid parameter adjustment and rotation.

[0009] In an optional embodiment, the major semi-axis and eccentricity of the reference ellipsoid are changed to establish a new ellipsoid, including: determining the change of the major semi-axis and eccentricity of the reference ellipsoid based on the change of the geoid height of the position reference point to the difference between the normal height of the position reference point and the normal height of the engineering projection surface; and obtaining the major semi-axis and eccentricity of the new ellipsoid based on the change of the major semi-axis and eccentricity of the reference ellipsoid.

[0010] In this embodiment, when establishing a new ellipsoid, in order to make the ellipsoidal surface of the new ellipsoid as close as possible to the engineering projection surface, the major semi-axis and eccentricity of the reference ellipsoid are changed so that the geoid height of the position reference point is changed to the difference between the normal height of the position reference point and the normal height of the engineering projection surface, thereby determining the major semi-axis and eccentricity of the new ellipsoid and realizing the establishment of the new ellipsoid.

[0011] In an optional embodiment, the geodetic coordinates of the position reference point, plane and elevation control point are converted into geodetic coordinates under the new ellipsoid, including: converting the geodetic coordinates of the position reference point, plane and elevation control point into spatial rectangular coordinates based on the major semi-axis and eccentricity of the reference ellipsoid; converting the spatial rectangular coordinates into geodetic coordinates under the new ellipsoid based on the major semi-axis and eccentricity of the new ellipsoid.

[0012] In this embodiment, when performing geodetic coordinate conversion, the geodetic coordinates may be first converted into spatial rectangular coordinates, and then the spatial rectangular coordinates may be converted into geodetic coordinates under the new ellipsoid, thereby realizing the conversion of geodetic coordinates under the reference ellipsoid and the new ellipsoid.

[0013] In an optional embodiment, the inclination angle between the ellipsoidal surface of the new ellipsoid engineering area and the engineering projection surface is determined based on the deviation of the geodetic coordinate change of the plane and elevation control points when the new ellipsoid is rotated, including: determining the geodetic height change based on the relationship between the rotation angle when the new ellipsoid is rotated and the geodetic height change of the plane and elevation control points; constructing an indirect adjustment model based on the deviation of the geodetic height change; solving the indirect adjustment model using the least squares method to obtain the rotation angle of the new ellipsoid, where the rotation angle is the inclination angle between the ellipsoidal surface of the new ellipsoid engineering area and the engineering projection surface.

[0014] In an optional implementation, the rotation angle of the new ellipsoid is determined by the rotation angle of the new ellipsoid around the equator and the rotation angle around the north and south poles with the position reference point as the center.

[0015] In this embodiment, when determining the rotation angle of the new ellipsoid, the rotation angle is determined based on the principle of minimizing the deviation of the change in the earth height after rotation, thereby improving the accuracy of the determined rotation angle.

[0016] In an optional implementation, the plane and elevation control points are rotated based on the inclination angle, and the geodetic coordinates of the rotated plane and elevation control points under the new ellipsoid are converted to obtain coordinates under the engineering independent coordinate system, thereby completing the establishment of the engineering ellipsoid, including: rotating the spatial rectangular coordinates of the plane and elevation control points based on the inclination angle to obtain the rotated spatial rectangular coordinates; converting the rotated spatial rectangular coordinates into the geodetic coordinates of the rotated plane and elevation control points under the new ellipsoid based on the major semi-axis and eccentricity of the new ellipsoid; performing projection conversion and one-point orientation conversion on the geodetic coordinates of the rotated plane and elevation control points under the new ellipsoid to obtain coordinates under the engineering independent coordinate system, thereby completing the establishment of the engineering ellipsoid.

[0017] In this embodiment, when establishing the engineering ellipsoid, the coordinates of the plane and elevation control points are rotated according to the rotation angle of the new ellipsoid, and then projection transformation and one-point orientation transformation are performed based on the rotated coordinates to obtain the coordinates in the engineering independent coordinate system, providing an accurate data basis for the subsequent calculation of the construction measurement control network.

[0018] In a second aspect, the present invention provides an engineering ellipsoid establishment device suitable for a construction measurement control network, the device comprising: a coordinate acquisition module, used to establish a reference ellipsoid, obtain the geodetic coordinates and normal heights of the position reference points, plane and elevation control points of the engineering area; a new ellipsoid establishment module, used to change the major semi-axis and eccentricity of the reference ellipsoid based on the geodetic coordinates and normal height to establish a new ellipsoid; a conversion module, used to convert the geodetic coordinates of the position reference points, plane and elevation control points into geodetic coordinates under the new ellipsoid; an inclination angle determination module, used to determine the inclination angle between the ellipsoidal surface of the new ellipsoid engineering area and the engineering projection surface based on the deviation of the geodetic coordinate change of the plane and elevation control points when the new ellipsoid is rotated; an engineering ellipsoid establishment module, used to rotate the plane and elevation control points based on the inclination angle, and convert the geodetic coordinates of the rotated plane and elevation control points under the new ellipsoid to obtain the coordinates under the engineering independent coordinate system, thereby completing the establishment of the engineering ellipsoid.

[0019] In a third aspect, the present invention provides a computer device, comprising: a memory and a processor, the memory and the processor being communicatively connected to each other, the memory storing computer instructions, and the processor executing the method for establishing an engineering ellipsoid applicable to a construction surveying control network of the above-mentioned first aspect or any corresponding embodiment thereof by executing the computer instructions.

[0020] In a fourth aspect, the present invention provides a computer-readable storage medium having computer instructions stored thereon, the computer instructions being used to enable a computer to execute the method for establishing an engineering ellipsoid applicable to a construction surveying control network of the above-mentioned first aspect or any corresponding embodiment thereof.

[0021] In a fifth aspect, the present invention provides a computer program product, comprising computer instructions for causing a computer to execute the method for establishing an engineering ellipsoid applicable to a construction surveying control network according to the first aspect or any corresponding embodiment thereof. BRIEF DESCRIPTION OF THE DRAWINGS

[0022] In order to more clearly illustrate the specific implementation methods of the present invention or the technical solutions in the prior art, the drawings required for use in the specific implementation methods or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some implementation methods of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0023] Figure 1 It is a flow chart of a method for establishing an engineering ellipsoid applicable to a construction measurement control network according to an embodiment of the present invention;

[0024] Figure 2 It is a flowchart of another method for establishing an engineering ellipsoid applicable to a construction measurement control network according to an embodiment of the present invention;

[0025] Figure 3 It is a structural block diagram of an engineering ellipsoid establishment device applicable to a construction measurement control network according to an embodiment of the present invention;

[0026] Figure 4 It is a schematic diagram of the hardware structure of a computer device according to an embodiment of the present invention. DETAILED DESCRIPTION

[0027] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative work are within the scope of protection of the present invention.

[0028] According to an embodiment of the present invention, an embodiment of a method for establishing an engineering ellipsoid suitable for a construction measurement control network is provided. It should be noted that the steps shown in the flowchart of the accompanying drawings can be executed in a computer system such as a set of computer executable instructions, and although a logical order is shown in the flowchart, in some cases, the steps shown or described can be executed in an order different from that shown here.

[0029] In this embodiment, a method for establishing an engineering ellipsoid suitable for a construction measurement control network is provided, which can be used in electronic devices such as computers, mobile phones, tablet computers, etc. Figure 1 is a flow chart of a method for establishing an engineering ellipsoid applicable to a construction measurement control network according to an embodiment of the present invention. Figure 1 As shown, the process includes the following steps:

[0030] Step S101, establish a reference ellipsoid, obtain the geodetic coordinates and normal height of the position reference point, plane and elevation control point of the project area. Specifically, the project area is a specific geographical space range defined when carrying out project construction, planning, design or construction activities. That is, before establishing the project ellipsoid, the project area must first be determined according to the actual project situation and the geographical boundaries of the project area must be divided.

[0031] Among them, the position reference point is the starting point or starting point for constructing the construction measurement control network (which can be referred to as the control network). For example, the center of the project area can be used as the position reference point. Multiple plane and elevation control points are arranged according to the preset design plan. That is, the control point includes plane coordinate information and elevation information.

[0032] After the location reference point, plane and elevation control point are arranged, it is also necessary to determine its actual coordinates in the engineering area. That is, the coordinates of the location reference point, plane and elevation control point in the design drawing are converted into real coordinates. When converting, it is necessary to determine the geographic coordinate system of the engineering area. In the present embodiment, the geographic coordinates of the engineering area are characterized by constructing a reference ellipsoid, which is a mathematical model for approximately describing the shape of the earth. Thus, the reference ellipsoid can be used to more accurately represent the surface of the earth. When establishing the reference ellipsoid, it is necessary to consider the engineering area situation to determine the parameters of the reference ellipsoid such as the major semi-axis, the minor semi-axis and the eccentricity, etc. At the same time, it is also necessary to select a suitable geographic coordinate system, for example, WGS84 (World Geodetic System-1984 Coordinate System, World Geodetic System 1984), CGCS2000 (China Geodetic Coordinate System 2000, China 2000 National Geodetic Coordinate System) can be selected as the geographic coordinate system. In the present embodiment, WGS84 is selected to construct the reference ellipsoid, so the constructed reference ellipsoid can also be called WGS84 ellipsoid.

[0033] After constructing the reference ellipsoid, the geodetic coordinates and normal heights of the position reference points, plane and elevation control points can be calculated and determined.

[0034] Step S102, based on the geodetic coordinates and the normal height, the major semi-axis and eccentricity of the reference ellipsoid are changed to establish a new ellipsoid. Specifically, when establishing a construction measurement control network, it includes two aspects: determining the ellipsoid and selecting the projection, that is, first of all, it is necessary to establish an engineering ellipsoid that conforms to the reality. In this embodiment, when the engineering ellipsoid is used, it is achieved by transforming the reference ellipsoid. Since in the construction project, the ground observation results need to be projected onto a certain engineering projection surface, therefore, the ellipsoidal surface of the obtained new ellipsoid can be made as close to the engineering projection surface as possible by changing the parameters of the reference ellipsoid, such as the major semi-axis and eccentricity. Among them, it should be noted that the engineering projection surface can be predetermined before the engineering ellipsoid is established.

[0035] Step S103, converting the geodetic coordinates of the position reference point, plane and elevation control point into geodetic coordinates under the new ellipsoid. Specifically, after the parameters of the new ellipsoid are adjusted, the geodetic coordinates of the original position reference point, plane and elevation control point will also change, so it is necessary to convert their geodetic coordinates into the new ellipsoid. During the conversion, the geodetic coordinates under the new ellipsoid can be determined according to the parameters of the new ellipsoid.

[0036] Step S104, based on the deviation of the geodetic coordinate change amount of the plane and elevation control points when the new ellipsoid rotates, the inclination angle between the ellipsoidal surface of the new ellipsoid engineering area and the engineering projection surface is determined. Specifically, although the ellipsoidal surface of the new ellipsoid engineering area is close to the engineering projection surface by changing the major semi-axis and eccentricity of the reference ellipsoid, there is still a certain angle deviation between the two, that is, the two are not parallel. Among them, the angular deviation (i.e., the inclination angle) between the ellipsoidal surface of the new ellipsoid engineering area and the engineering projection surface is expressed as the deviation between the normal direction of the ellipsoidal surface of the new ellipsoid engineering area and the vertical direction of the engineering projection surface, and the deviation can be eliminated specifically by rotating the new ellipsoid. The angle of rotation of the new ellipsoid is the inclination angle between the ellipsoidal surface of the new ellipsoid engineering area and the engineering projection surface. In the present embodiment, when determining the inclination angle, it is determined with the minimum deviation of the geodetic height change amount when the new ellipsoid rotates, thereby ensuring the accuracy of rotation.

[0037] Step S105, rotate the plane and elevation control points based on the inclination angle, and convert the geodetic coordinates of the rotated plane and elevation control points under the new ellipsoid to obtain the coordinates under the engineering independent coordinate system, and complete the establishment of the engineering ellipsoid. Specifically, after determining the rotation angle of the new ellipsoid, i.e., the inclination angle, the geodetic coordinates of the plane and elevation control points are also rotated, and then converted based on the rotated geodetic coordinates with reference to relevant technologies, so as to obtain the coordinates under the engineering independent coordinate system, and complete the establishment of the engineering ellipsoid. As a result, the back-calculated distance of the projection coordinates under the engineering ellipsoid is consistent with the actual electromagnetic wave ranging. This solves the problem that the reference ellipsoid and the projection surface of the engineering area cannot be well overlapped when the existing GNSS method is used to carry out construction measurement control network.

[0038] The method for establishing an engineering ellipsoid applicable to a construction measurement control network provided by an embodiment of the present invention adjusts the parameters of a reference ellipsoid of an engineering area so that it is as close to the engineering projection surface as possible, and then further rotates the ellipsoid to eliminate the deviation between the normal direction of the ellipsoid surface and the vertical direction of the engineering projection surface, and based on the adjustment and rotation of the ellipsoid parameters, the coordinates of each control point are converted to obtain the coordinates in the engineering independent coordinate system, and the engineering ellipsoid is established. In this way, the inclination angle between the ellipsoid surface and the engineering projection surface can be eliminated, so that the inverse distance of the projection coordinates under the new ellipsoid is consistent with the actual electromagnetic wave ranging, and it is easy to use.

[0039] In this embodiment, a method for establishing an engineering ellipsoid applicable to a construction measurement control network is provided, and the process includes the following steps:

[0040] Step S201, establish a reference ellipsoid, obtain the geodetic coordinates and normal heights of the location reference points, plane and elevation control points of the project area.

[0041] Specifically, the above step S201 includes:

[0042] Step S2011, establish a reference ellipsoid, determine the geodetic longitude, geodetic latitude, geodetic height and normal height of the location reference point based on the project area, and determine the normal height of the project projection surface.

[0043] Specifically, when determining the geodetic longitude and latitude, the geodetic longitude B0, geodetic latitude L0 and geodetic height H0 of the position reference point in the reference ellipsoid coordinate system can be obtained by using equipment such as GNSS. The normal height h0 is obtained by the leveling method.

[0044] In addition, when establishing the engineering ellipsoid, it is also necessary to determine the engineering projection surface for the projection of ground observation results. Specifically, when determining the normal height of the engineering projection surface, it is usually the average elevation surface of the engineering area, or an elevation surface specified according to the characteristics of the project.

[0045] Step S2012, using the GNSS method to calculate the geodetic longitude, latitude and height of each plane and elevation control point; specifically, during the calculation, a GNSS receiver can be set up at each plane and elevation control point, the receiver receives satellite observation signals through observation, and then uses relevant processing software to pre-process the signals, perform baseline calculations, and perform network adjustment, etc., to obtain the geodetic longitude, latitude and height of the corresponding control point. Among them, the plane and elevation control point P i The geodetic longitude, geodetic longitude and geodetic altitude (i=1,2,3,...) can be expressed as: geodetic latitude B=[B1B2…B i …B n- 1B n ] T 、Geod longitude L=[L1L2…L i …L n-1 L n ] T 、Earth height H=[H1H2…H i …H n-1 H n ] T .

[0046] Step S2013, using the leveling method to calculate the leveling elevation of each plane and elevation control point. Specifically, when using the leveling method to calculate the leveling elevation, it can be based on equipment such as a level and a leveling rod, and a point with a known elevation can be selected as the starting point, and then the measurement is performed according to the pre-set leveling route, and the height difference between two adjacent points is calculated based on the measured value, and the leveling elevation of each plane and elevation control point is obtained in combination with the elevation of the starting point. Among them, the plane and elevation control point P i The level elevation of (i=1,2,3,……) can be expressed as:

[0047] [h1h2…h i …h n-1 hn ] T .

[0048] Step S202, the major semi-axis and eccentricity of the reference ellipsoid are changed based on the geodetic coordinates and the normal height to create a new ellipsoid.

[0049] Specifically, the above step S202 includes:

[0050] Step S2021, when the geoid height of the position reference point is changed to the difference between the normal height of the position reference point and the normal height of the engineering projection surface, the changes of the major semi-axis and the eccentricity of the reference ellipsoid are determined.

[0051] Step S2022, based on the changes of the semi-major axis and eccentricity of the reference ellipsoid, obtain the semi-major axis and eccentricity of the new ellipsoid.

[0052] Specifically, in order to make the ellipsoid surface of the new ellipsoid engineering area as close to the engineering projection surface as possible, the major semi-axis a and eccentricity e of the reference ellipsoid can be changed without changing the positioning and orientation of the original ellipsoid (i.e., the reference ellipsoid). 2 , so that the geodetic height of the position reference point P0 changes from H0 to h0-h, while keeping the normal direction of the position reference point P0 under the new ellipsoid and the geodetic longitude and latitude B0, L0 unchanged. Therefore, we can first calculate the changes in the semi-major axis and eccentricity of the reference ellipsoid when the geodetic height of the position reference point P0 changes from H0 to the difference h0-h between the normal height of the position reference point and the normal height of the engineering projection surface, and then obtain the semi-major axis and eccentricity of the new ellipsoid based on the change and the semi-major axis and eccentricity of the reference ellipsoid. Specifically, the semi-major axis a1 and eccentricity of the new ellipsoid They are:

[0053] a1=a+da1 formula (1)

[0054]

[0055] In formula (1),

[0056]

[0057] In formula (2),

[0058]

[0059] Among them, da1 represents the change of the major semi-axis of the reference ellipsoid, de 2 represents the change of the eccentricity of the reference ellipsoid. Therefore, the semi-major axis is a1 and the square of the eccentricity is The new ellipsoid.

[0060] Step S203, converting the geodetic coordinates of the position reference point, plane and elevation control point into geodetic coordinates under the new ellipsoid; specifically, during the conversion, the geodetic coordinates may be first converted into spatial rectangular coordinates, and then converted into geodetic coordinates under the new ellipsoid.

[0061] Specifically, the above step S203 includes:

[0062] Step S2031, based on the major semi-axis and eccentricity of the reference ellipsoid, the geodetic coordinates of the position reference point, the plane control point and the elevation control point are converted into spatial rectangular coordinates; specifically, the conversion process can be implemented by referring to the following formula:

[0063]

[0064] In formula (5), for the obtained X i , Y i , Z i When i takes different values ​​such as 0, 1, 2, 3...n, it represents the spatial rectangular coordinates corresponding to the position reference point, plane control point and elevation control point, N i is the radius of the truncated circle of the reference ellipsoid, and its calculation formula is:

[0065]

[0066] Step S2032: based on the major semi-axis and eccentricity of the new ellipsoid, the spatial rectangular coordinates are converted into geodetic coordinates under the new ellipsoid. Specifically, the conversion process is implemented using the following formula:

[0067]

[0068] In formula (7), for the obtained B i ′, L′ i , H i ′, when i takes different values ​​such as 0, 1, 2, 3...n, it represents the geodetic coordinates of the position reference point, plane and elevation control point under the new ellipsoid, where B i ′ needs to be calculated iteratively, that is, substitute formula (8) into the second formula in formula (7), and iterate to make B i The difference between the two is less than 1e10 -10 Stop iteration when N i ' is the radius of the umbo circle of each point in the new ellipsoid, and its calculation formula is:

[0069]

[0070] Step S204, determining the inclination angle between the new ellipsoid and the engineering projection surface based on the deviation of the geodetic coordinate change of the plane and elevation control points when the new ellipsoid rotates.

[0071] Specifically, the above step S204 includes:

[0072] Step S2041, determining the geodetic height change amount based on the relationship between the rotation angle of the new ellipsoid when rotating and the geodetic height change amount of the plane and the elevation control point. The rotation angle of the new ellipsoid when rotating is determined by the rotation angle of the new ellipsoid when rotating around the equator with the position reference point as the center and the rotation angle around the north and south poles.

[0073] Specifically, for the new ellipsoid obtained by changing the parameters, since the geodetic longitude and latitude of the position reference point remain unchanged when adjusting the parameters, and the orientation and positioning of the ellipsoid are not changed, the normal of the new ellipsoid is consistent with the normal of the reference ellipsoid, but the ellipsoid surface of the new ellipsoid engineering area is not parallel to the engineering projection surface. Therefore, in this embodiment, the new ellipsoid is rotated around the equator and the north and south poles with the position reference point as the center, and the rotation angle θ of the new ellipsoid is determined by the rotation angle θ1 around the equator and the rotation angle θ2 around the north and south poles, that is, θ=[θ1θ2] T .

[0074] In addition, when determining the rotation angle, it is also necessary to further consider the geodetic height change of the plane and elevation control points caused by the rotation of the new ellipsoid. Among them, the rotation angle of the new ellipsoid and the geodetic height change of the plane control point and the elevation control point dH i The relationship between them is expressed by the following formula:

[0075] dH i = = F2θ2 + F1θ1 Formula (9)

[0076] in,

[0077]

[0078] Step S2042, constructing an indirect adjustment model based on the deviation of the geoid height variation. Specifically, the deviation of the geoid height variation obtained based on the geoid height variation is expressed by the following formula:

[0079] v i =dH i +H i -h i +h formula (12)

[0080] Substituting formula (9) into formula (12), we can get the standard form of the indirect adjustment model:

[0081] v i =a i θ+l i Formula (13)

[0082] In formula (13),

[0083] ai =(F2, F1) Formula (14)

[0084] θ=[θ1θ2] T Formula (15)

[0085] l i =H i -h i +h formula (16)

[0086] Step S2043, using the least square method to solve the indirect adjustment model, to obtain the rotation angle of the new ellipsoid, the rotation angle is the inclination angle between the ellipsoid surface of the new ellipsoid project area and the project projection surface. Specifically, the least square method can be used to solve the indirect adjustment model to obtain the rotation angle of the new ellipsoid, the rotation angle can be expressed by the following formula:

[0087] θ=(A T A) -1 A T L-type (17)

[0088] In formula (17), the matrix A is a i is the i-th row vector, and the column vector L is l i is the ith component.

[0089] Step S205, rotating the plane and the elevation control points based on the inclination angle, and converting the geodetic coordinates of the rotated plane and the elevation control points under the new ellipsoid to obtain the coordinates under the engineering independent coordinate system, thereby completing the establishment of the engineering ellipsoid.

[0090] Specifically, the above step S205 includes:

[0091] Step S2051, based on the inclination angle, the spatial rectangular coordinates of the plane and the elevation control point are rotated to obtain the rotated spatial rectangular coordinates; specifically, the rotated spatial rectangular coordinates (X′ i , Y′ i , Z′ i ) is expressed by the following formula:

[0092]

[0093] In formula (18),

[0094]

[0095] Among them, R1 and R2 represent the rotation matrix during the spatial coordinate rotation transformation.

[0096] Step S2052, based on the major semi-axis and eccentricity of the new ellipsoid, the rotated spatial rectangular coordinates are converted into the geodetic coordinates of the rotated plane and the elevation control point under the new ellipsoid; specifically, after obtaining the rotated spatial rectangular coordinates, the geodetic coordinates of the rotated spatial rectangular coordinates under the new ellipsoid can be calculated according to the above formula 7 ( Bi新椭球 , L i新椭球 , H i新椭球 ).

[0097] Step S2053, the geodetic coordinates of the rotated plane and the elevation control point under the new ellipsoid are projected and transformed, and the coordinates under the independent coordinate system of the project are obtained, and the establishment of the project ellipsoid is completed. Specifically, during the conversion, the rectangular coordinates of the Gaussian plane can be calculated according to the central meridian and the elevation of the projection surface and other parameters formulated in the project area, and the "one point one direction" conversion is performed according to the fixed coordinates and the fixed coordinate azimuth of the project starting point to obtain the coordinates under the independent coordinate system of the project.

[0098] When performing the conversion, first select a suitable central meridian according to the location of the project area. Then calculate the rectangular coordinates of the Gaussian plane. When calculating, you can use the Gaussian projection formula to convert the longitude and latitude coordinates into rectangular coordinates (x, y) on the Gaussian plane. When converting, you need to consider the new ellipsoid parameters (such as the major semi-axis, minor semi-axis, eccentricity, etc.), as well as the elevation of the projection surface. After calculating the rectangular coordinates of the Gaussian plane, perform a point-to-azimuth conversion. When converting, first select a known point with accurate coordinates in the national coordinate system as the starting point of the project. Then translate the rectangular coordinates of the Gaussian plane to the Gaussian plane coordinates corresponding to the starting point of the project, and then rotate the Gaussian plane coordinates according to the fixed coordinate azimuth to make it conform to the definition of the engineering coordinate system.

[0099] As a specific application example of the embodiment of the present invention, Figure 2 As shown, the method for establishing the engineering ellipsoid applicable to the construction measurement control network can be implemented by the following process:

[0100] Step 1: Select the location reference point P0 in the center of the project area (survey area), set the normal height of the project projection surface as h, the geodetic longitude and latitude and geodetic height of the location reference point P0 as B0, L0 and H0, and the normal height as h0.

[0101] Step 2: Change the semi-major axis a and eccentricity e of the reference ellipsoid without changing the positioning and orientation of the reference ellipsoid. 2 , so that the geodetic height of the position reference point P0 changes from H0 to h0-h, while keeping the normal direction and geodetic longitude and latitude B0 and L0 of the position reference point P0 under the new ellipsoid unchanged, thus obtaining the major semi-axis a1 and eccentricity of the new ellipsoid

[0102] Step 3: Evenly distribute several plane and elevation control points P in the project area i (i=1,2,3,……), and the geodetic latitude of each control point is obtained by GNSS method B=[B1B2…B i …B n-c B n ] T , geodetic longitude L = [L ′ L2…L i …L n-1 L n ] T 、Earth height H=[H1H2…H i …H n-1 H n ] T , the leveling method is used to obtain the leveling height of each control point h = [h1h2…h i …h n-1 h n ] T .

[0103] Step 4: Position reference point P0 and each plane and elevation control point P i The longitude and latitude of the geodetic height B i , L i and H i According to formula (5), it is converted into spatial rectangular coordinates.

[0104] Step 5: Position reference point P0 and each plane and elevation control point P i The spatial rectangular coordinate X i , Y i , Z i Convert the longitude, latitude and geodetic height under the new ellipsoid according to the parameters of the new ellipsoid.

[0105] Step 6: To eliminate the deviation between the normal direction of the new ellipsoid and the vertical direction of the project projection surface, that is, the inclination angle θ between the new ellipsoid and the projection surface, take the position reference point P0 as the center, rotate around the equator by an angle θ1, and rotate around the north and south poles by an angle θ2, that is, θ = [θ1θ2] T .

[0106] Step 7: At each control point P i At the point, the relationship between the change in geodetic height and the rotation angle component is determined. Then, based on the longitude and latitude of each control point and the normal height, the error equations (indirect adjustment model) shown in equations (13), (14), (15), and (16) are constructed with the principle of minimizing the change in geodetic height after rotation.

[0107] Step 8: Solve the error equation using the least squares method to obtain the rotation angle.

[0108] Step 9: Rotate the spatial coordinates of each control point according to the obtained rotation angle to obtain the spatial rectangular coordinates after the rotation transformation.

[0109] Step 10: The spatial rectangular coordinates of each control point after rotation transformation are calculated according to the new ellipsoid parameters a1, Convert to geodetic coordinates.

[0110] Step 11: The geodetic coordinates of each control point under the new ellipsoid after rotation transformation ( Bi新椭球 , L i新椭球 , H i新椭球 ) Calculate the rectangular coordinates of the Gaussian plane according to the central meridian, projection surface elevation and other parameters set for the project area, and perform a "one point one azimuth" conversion according to the fixed coordinates and fixed coordinate azimuth of the project starting point to obtain the coordinates in the independent coordinate system of the project.

[0111] In this embodiment, a device for establishing an engineering ellipsoid suitable for a construction measurement control network is also provided, and the device is used to implement the above-mentioned embodiments and preferred implementation modes, and the descriptions that have been made will not be repeated. As used below, the term "module" can implement a combination of software and / or hardware for a predetermined function. Although the device described in the following embodiments is preferably implemented in software, the implementation of hardware, or a combination of software and hardware, is also possible and conceivable.

[0112] This embodiment provides an engineering ellipsoid establishment device suitable for construction measurement control network, such as Figure 3 As shown, including:

[0113] The coordinate acquisition module 31 is used to establish a reference ellipsoid and obtain the geodetic coordinates and normal height of the position reference points, plane and elevation control points of the project area;

[0114] A new ellipsoid establishment module 32 is used to change the major semi-axis and eccentricity of the reference ellipsoid based on the geodetic coordinates and the normal height to establish a new ellipsoid;

[0115] A conversion module 33, used for converting the geodetic coordinates of the position reference point, plane and elevation control points into geodetic coordinates under the new ellipsoid;

[0116] The tilt angle determination module 34 is used to determine the tilt angle between the ellipsoid surface of the new ellipsoid engineering area and the engineering projection surface based on the deviation of the geodetic coordinate change amount of the plane and elevation control points when the new ellipsoid rotates;

[0117] The engineering ellipsoid establishment module 35 is used to rotate the plane and elevation control points based on the inclination angle, and convert the geodetic coordinates of the rotated plane and elevation control points under the new ellipsoid to obtain the coordinates under the engineering independent coordinate system, thereby completing the establishment of the engineering ellipsoid.

[0118] The further functional description of each of the above modules and units is the same as that of the above corresponding embodiments and will not be repeated here.

[0119] The embodiment of the present invention also provides a computer device having the above Figure 3 The engineering ellipsoid establishment device suitable for construction measurement control network is shown.

[0120] See also Figure 4 , Figure 4 is a schematic diagram of the structure of a computer device provided by an optional embodiment of the present invention, such as Figure 4 As shown, the computer device includes: one or more processors 10, a memory 20, and interfaces for connecting various components, including high-speed interfaces and low-speed interfaces. Various components are connected to each other using different buses for communication, and can be installed on a common mainboard or installed in other ways as needed. The processor can process the instructions executed in the computer device, including instructions stored in or on the memory to display the graphical information of the GUI on an external input / output device (such as, a display device coupled to the interface). In some optional embodiments, if necessary, multiple processors and / or multiple buses can be used together with multiple memories and multiple memories. Similarly, multiple computer devices can be connected, and each device provides some necessary operations (for example, as a server array, a group of blade servers, or a multi-processor system). Figure 4 A processor 10 is taken as an example.

[0121] The processor 10 may be a central processing unit, a network processor or a combination thereof. The processor 10 may further include a hardware chip. The hardware chip may be a dedicated integrated circuit, a programmable logic device or a combination thereof. The programmable logic device may be a complex programmable logic device, a field programmable gate array, a general purpose array logic or any combination thereof.

[0122] The memory 20 stores instructions executable by at least one processor 10, so that at least one processor 10 executes the method shown in the above embodiment.

[0123] The memory 20 may include a program storage area and a data storage area, wherein the program storage area may store an operating system, an application required for at least one function; the data storage area may store data created by the use of a computer device based on the presentation of a small program landing page, etc. In addition, the memory 20 may include a high-speed random access memory, and may also include a non-transient memory, such as at least one disk storage device, a flash memory device, or other non-transient solid-state storage device. In some optional embodiments, the memory 20 may optionally include a memory remotely arranged relative to the processor 10, and these remote memories may be connected to the computer device via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0124] The memory 20 may include a volatile memory, such as a random access memory; the memory may also include a non-volatile memory, such as a flash memory, a hard disk or a solid state drive; the memory 20 may also include a combination of the above types of memory.

[0125] The computer device further comprises a communication interface 30 for the computer device to communicate with other devices or a communication network.

[0126] The embodiment of the present invention also provides a computer-readable storage medium. The method according to the embodiment of the present invention can be implemented in hardware, firmware, or can be implemented as a computer code that can be recorded in a storage medium, or can be implemented as a computer code that is originally stored in a remote storage medium or a non-temporary machine-readable storage medium and will be stored in a local storage medium through a network download, so that the method described herein can be stored in such software processing on a storage medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware. Among them, the storage medium can be a magnetic disk, an optical disk, a read-only storage memory, a random access memory, a flash memory, a hard disk or a solid-state hard disk, etc.; further, the storage medium can also include a combination of the above types of memories. It can be understood that a computer, a processor, a microprocessor controller, or programmable hardware includes a storage component that can store or receive software or computer code. When the software or computer code is accessed and executed by a computer, a processor, or hardware, the method shown in the above embodiment is implemented.

[0127] A part of the present invention may be applied as a computer program product, such as a computer program instruction, which, when executed by a computer, can call or provide the method and / or technical solution according to the present invention through the operation of the computer. Those skilled in the art should understand that the existence of the computer program instruction in a computer-readable medium includes, but is not limited to, a source file, an executable file, an installation package file, etc., and accordingly, the way in which the computer program instruction is executed by the computer includes, but is not limited to: the computer directly executes the instruction, or the computer compiles the instruction and then executes the corresponding compiled program, or the computer reads and executes the instruction, or the computer reads and installs the instruction and then executes the corresponding installed program. Here, the computer-readable medium may be any available computer-readable storage medium or communication medium accessible to the computer.

[0128] Although the embodiments of the present invention have been described in conjunction with the accompanying drawings, those skilled in the art may make various modifications and variations without departing from the spirit and scope of the present invention, and such modifications and variations are all within the scope defined by the appended claims.

Claims

1. A method for establishing an engineering ellipsoid suitable for a construction measurement control network, characterized in that: The method comprises: Establish a reference ellipsoid to obtain the geodetic coordinates and normal height of the location reference points, plane and elevation control points in the project area; Changing the major semi-axis and eccentricity of the reference ellipsoid based on the geodetic coordinates and the normal height to establish a new ellipsoid; Convert the geodetic coordinates of the position reference points, plane and elevation control points into geodetic coordinates under the new ellipsoid; Determine the inclination angle between the ellipsoid surface of the new ellipsoid project area and the project projection surface based on the deviation of the geodetic coordinate change of the plane control point and the elevation control point when the new ellipsoid rotates; The plane and the elevation control points are rotated based on the inclination angle, and the geodetic coordinates of the rotated plane and the elevation control points under the new ellipsoid are converted to obtain the coordinates under the engineering independent coordinate system, thereby completing the establishment of the engineering ellipsoid.

2. The method according to claim 1, characterized in that Obtain the geodetic coordinates and normal heights of the location benchmarks, plane and elevation control points of the project area, including: Determine the geodetic longitude, geodetic latitude, geodetic height and normal height of the location reference point based on the project area, and determine the normal height of the project projection surface; The GNSS method is used to calculate the geodetic longitude, latitude and height of each plane and elevation control point; The leveling method is used to calculate the level elevation of each plane and elevation control point.

3. The method according to claim 2, characterized in that The major semi-axis and eccentricity of the reference ellipsoid are changed based on the geodetic coordinates and the normal height to establish a new ellipsoid, including: When the geodetic height of the position reference point is changed to the difference between the normal height of the position reference point and the normal height of the engineering projection surface, the change of the major semi-axis and eccentricity of the reference ellipsoid is determined; Based on the changes of the semi-major axis and eccentricity of the reference ellipsoid, the semi-major axis and eccentricity of the new ellipsoid are obtained.

4. The method according to claim 1, characterized in that: Convert the geodetic coordinates of the position reference points, plane and elevation control points to the geodetic coordinates under the new ellipsoid, including: Based on the major semi-axis and eccentricity of the reference ellipsoid, the geodetic coordinates of the position reference points, plane and elevation control points are converted into spatial rectangular coordinates; The spatial rectangular coordinates are converted into geodetic coordinates under the new ellipsoid based on the major semi-axis and eccentricity of the new ellipsoid.

5. The method according to claim 2, characterized in that: The inclination angle between the ellipsoid surface of the new ellipsoid project area and the project projection surface is determined based on the deviation of the geodetic coordinate change of the plane and elevation control points when the new ellipsoid rotates, including: Determine the geodetic height variation based on the relationship between the rotation angle of the new ellipsoid and the geodetic height variation of the plane and elevation control points; Constructing an indirect adjustment model based on the deviation of the geodetic height change; The indirect adjustment model is solved by the least square method to obtain the rotation angle of the new ellipsoid, which is the inclination angle between the ellipsoid surface of the new ellipsoid engineering area and the engineering projection surface.

6. The method according to claim 5, characterized in that The rotation angle of the new ellipsoid is determined by the rotation angle of the new ellipsoid around the equator with the position reference point as the center and the rotation angle around the north and south poles.

7. The method according to claim 4, characterized in that The plane and the elevation control point are rotated based on the inclination angle, and the geodetic coordinates of the rotated plane and the elevation control point under the new ellipsoid are converted to obtain the coordinates under the engineering independent coordinate system, and the engineering ellipsoid is established, including: Rotating the spatial rectangular coordinates of the plane and the elevation control point based on the inclination angle to obtain the rotated spatial rectangular coordinates; Based on the major semi-axis and eccentricity of the new ellipsoid, the rotated spatial rectangular coordinates are converted into the geodetic coordinates of the rotated plane and elevation control points under the new ellipsoid; The geodetic coordinates of the rotated plane and the elevation control point under the new ellipsoid are projected and transformed, and one-point orientation is transformed to obtain the coordinates under the engineering independent coordinate system, thus completing the establishment of the engineering ellipsoid.

8. An engineering ellipsoid establishment device suitable for construction measurement control network, characterized in that: The device comprises: The coordinate acquisition module is used to establish a reference ellipsoid and obtain the geodetic coordinates and normal height of the location reference points, plane and elevation control points of the project area; A new ellipsoid establishment module, used for changing the major semi-axis and eccentricity of the reference ellipsoid based on the geodetic coordinates and normal height to establish a new ellipsoid; A conversion module, used to convert the geodetic coordinates of the position reference point, plane and elevation control point into the geodetic coordinates under the new ellipsoid; The inclination angle determination module is used to determine the inclination angle between the ellipsoid surface of the new ellipsoid engineering area and the engineering projection surface based on the deviation of the geodetic coordinate change amount of the plane and elevation control points when the new ellipsoid rotates; The engineering ellipsoid establishment module is used to rotate the plane and elevation control points based on the inclination angle, and convert the geodetic coordinates of the rotated plane and elevation control points under the new ellipsoid to obtain the coordinates under the engineering independent coordinate system, thereby completing the establishment of the engineering ellipsoid.

9. A computer device, characterized in that: include: A memory and a processor, wherein the memory and the processor are communicatively connected to each other, the memory stores computer instructions, and the processor executes the method for establishing an engineering ellipsoid applicable to a construction surveying control network as described in any one of claims 1 to 7 by executing the computer instructions.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a computer to execute the method for establishing an engineering ellipsoid applicable to a construction survey control network as described in any one of claims 1 to 7.

Images