CN111721272A - Engineering surface measurement method based on ellipsoid calculation - Google Patents
Engineering surface measurement method based on ellipsoid calculation Download PDFInfo
- Publication number
- CN111721272A CN111721272A CN202010567367.7A CN202010567367A CN111721272A CN 111721272 A CN111721272 A CN 111721272A CN 202010567367 A CN202010567367 A CN 202010567367A CN 111721272 A CN111721272 A CN 111721272A
- Authority
- CN
- China
- Prior art keywords
- area
- engineering
- measured
- ellipsoid
- earth
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/2433—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures for measuring outlines by shadow casting
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/28—Measuring arrangements characterised by the use of optical techniques for measuring areas
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The invention discloses an engineering surface measurement method based on ellipsoid calculation, which comprises the following steps: s1, converting the space coordinates of each sampling point on the earth surface to an earth ellipsoid through Gaussian coordinates to obtain the earth coordinates of each sampling point; s2, determining the range of the engineering area to be tested on the earth ellipsoid through the point location relation according to the geodetic coordinates of the sampling points; s3, calculating the meridian curvature radius and the prime circle curvature radius of the engineering area to be measured according to the geodetic coordinates of the sampling points; s4, calculating the area of the engineering area to be measured according to the ellipsoidal meridian curvature radius and the prime circle curvature radius; the method is suitable for calculating the engineering surface area with large area or large topographic relief, and solves the problems that the prior projection method projects the measurement observation value onto a plane and then calculates the surface area, and the projection error is generated when the measurement observation value is projected onto the plane, so that the engineering surface area calculation is influenced.
Description
Technical Field
The invention belongs to the technical field of measurement, and particularly relates to an engineering surface measurement method based on ellipsoid calculation.
Background
At present, in actual engineering measurement operation, for an engineering surface area calculation method, a traditional projection method is generally adopted to project a measurement observation value onto a plane and then calculate the surface area, but any projection condition cannot avoid deformation in the earth ellipsoid projection process, so that various errors caused by the influence of the projection condition on the observation value in the projection process onto the plane cannot be avoided, and further errors are brought to the engineering surface area calculation.
The existing method for calculating the surface area of the irregular terrain is more, and the more common methods at home and abroad comprise the following steps: the method comprises the steps of generally projecting an observation value to a plane, and then calculating the area by using a plane analysis method, wherein in the process, the projection error is generated when the measurement observation value is projected to the plane, so that the calculation of the engineering surface area is influenced. For the actual engineering, the method generates larger projection errors for the ground surface measurement with large span range or the area with large elevation relief.
For example, in practical engineering, a calculation principle based on a helenlun formula is often adopted, and the helenlun formula is a method for solving an area by using three side lengths of a triangle. When the surface area of the irregular terrain is calculated, the irregular terrain is approximately divided into a plurality of small triangles, the area of each triangle is calculated based on the Helen formula and the Excel table and summed, and finally the approximate surface area of the terrain is obtained. However, this method is only suitable for small area surface calculation, and the difference in the control of the side length of the divided triangle may cause different surface area calculation results, so this method is not suitable for large area or large relief engineering surface area calculation.
Disclosure of Invention
Aiming at the defects in the prior art, the engineering surface measurement method based on ellipsoid calculation is suitable for large-area or large-relief engineering surface area calculation, and solves the problems that the existing projection method projects the measurement observation value onto a plane and then carries out surface area calculation, and projection errors are generated when the measurement observation value is projected onto the plane, so that the engineering surface area calculation is influenced.
In order to achieve the purpose of the invention, the invention adopts the technical scheme that: an engineering surface measurement method based on ellipsoid calculation comprises the following steps:
s1, performing inverse calculation on the spatial coordinates of the various sampling points on the earth to an ellipsoid of the earth through Gaussian coordinates to obtain the geodetic coordinates of the various sampling points;
s2, determining the range of the engineering area to be tested on the earth ellipsoid through the point location relation according to the geodetic coordinates of the sampling points;
s3, calculating the meridian curvature radius and the prime circle curvature radius of the engineering area to be measured according to the geodetic coordinates of the sampling points;
and S4, calculating the area of the engineering area to be measured according to the ellipsoidal meridian curvature radius and the prime circle curvature radius.
Further, the ellipsoidal meridian radius of curvature and the prime circle radius of curvature of the engineering area to be measured in step S3 satisfy the following formulas:
N=a(1-e2sin2(B2-B1))
wherein M is the radius of curvature of the meridian, N is the radius of curvature of the prime circle, a is the half-axis of the ellipsoid of the earth, e is the eccentricity of the engineering area to be measured, B2Upper boundary latitude, B, of the project area to be measured1The latitude of the lower boundary of the engineering area to be measured.
Further, the area of the engineering region to be measured in step S4 satisfies the following formula:
wherein P is the area of the engineering region to be measured, L2Is the right boundary longitude, L of the project area to be measured1The left boundary longitude of the project area to be tested.
The invention has the beneficial effects that:
1. the method does not need to adopt the traditional geometric analysis plane to calculate the surface area of the projection plane, and compared with the traditional plane analysis calculation method, the method has the advantages that the projection deformation is avoided or the deformation is minimized in the process of changing the measurement result into the ellipsoid, so that the error generated in the process of projecting the ground observation station to the plane in the traditional method is effectively avoided, and the purpose of no deformation or minimum deformation influence on the projection of the measurement calculation result is achieved.
2. The invention integrates the transformation and deformation calculation of various data on a reference ellipsoid surface and a Gaussian surface; the data calculated on the reference ellipsoid can obtain results on various projection surfaces according to various projection rules, so that the calculated data has reusability and universality.
3. In the actual engineering measurement, the geodetic coordinates of the measurement sampling points can be directly used for calculating the surface area of the engineering area to be measured, and the method is more suitable for the actual engineering measurement.
Drawings
FIG. 1 is a flow chart of an engineering surface measurement method based on ellipsoid calculations.
Fig. 2 is a schematic diagram of the transformation of the spatial coordinates of the sampling points onto the ellipsoid of the earth.
Detailed Description
The following description of the embodiments of the present invention is provided to facilitate the understanding of the present invention by those skilled in the art, but it should be understood that the present invention is not limited to the scope of the embodiments, and it will be apparent to those skilled in the art that various changes may be made without departing from the spirit and scope of the invention as defined and defined in the appended claims, and all matters produced by the invention using the inventive concept are protected.
As shown in fig. 1, an engineering surface measurement method based on ellipsoid calculation includes the following steps:
s1, performing inverse calculation on the spatial coordinates of the various sampling points on the earth to an ellipsoid of the earth through Gaussian coordinates to obtain the geodetic coordinates of the various sampling points;
in this embodiment, as shown in fig. 2, four sampling points Q1, N1, M1, and P1 are located on the earth surface, and are transformed to the earth ellipsoid through gaussian coordinates, so as to obtain respective geodetic coordinates a1, B1, C1, and D1 of the four sampling points.
S2, determining the range of the engineering area to be tested on the earth ellipsoid through the point location relation according to the geodetic coordinates of the sampling points;
s3, calculating the meridian curvature radius and the prime circle curvature radius of the engineering area to be measured according to the geodetic coordinates of the sampling points;
the ellipsoidal meridian radius of curvature and the unitary radius of curvature of the engineering area to be measured in the step S3 satisfy the following formulas:
N=a(1-e2sin2(B2-B1)) (2)
wherein M is the radius of curvature of the meridian, N is the radius of curvature of the prime circle, a is the half-axis of the ellipsoid of the earth, e is the eccentricity of the engineering area to be measured, B2Upper boundary latitude, B, of the project area to be measured1The latitude of the lower boundary of the engineering area to be measured.
And S4, calculating the area of the engineering area to be measured according to the ellipsoidal meridian curvature radius and the prime circle curvature radius.
The area of the engineering region to be measured in the step S4 satisfies the following formula:
after the project area range is determined through steps S2 and S3, equations (1) and (2) are taken into equation (3), and the equations are expanded as follows:
wherein P is the area of the engineering region to be measured, L2Is the right boundary longitude, L of the project area to be measured1The longitude of the left boundary of the engineering area to be measured, d curved surface infinitesimal, B ═ B2-B1,L=L2-L1And D is the range of the engineering area to be measured.
The invention has the beneficial effects that:
1. the method does not need to adopt the traditional geometric analysis plane to calculate the surface area of the projection plane, and compared with the traditional plane analysis calculation method, the method has the advantages that the projection deformation is avoided or the deformation is minimized in the process of changing the measurement result into the ellipsoid, so that the error generated in the process of projecting the ground observation station to the plane in the traditional method is effectively avoided, and the purpose of no deformation or minimum deformation influence on the projection of the measurement calculation result is achieved.
2. The invention integrates the transformation and deformation calculation of various data on a reference ellipsoid surface and a Gaussian surface; the data calculated on the reference ellipsoid can obtain results on various projection surfaces according to various projection rules, so that the calculated data has reusability and universality.
3. In the actual engineering measurement, the geodetic coordinates of the measurement sampling points can be directly used for calculating the surface area of the engineering area to be measured, and the method is more suitable for the actual engineering measurement.
Claims (3)
1. An engineering surface measurement method based on ellipsoid calculation is characterized by comprising the following steps:
s1, performing inverse calculation on the spatial coordinates of the various sampling points on the earth to an ellipsoid of the earth through Gaussian coordinates to obtain the geodetic coordinates of the various sampling points;
s2, determining the range of the engineering area to be tested on the earth ellipsoid through the point location relation according to the geodetic coordinates of the sampling points;
s3, calculating the meridian curvature radius and the prime circle curvature radius of the engineering area to be measured according to the geodetic coordinates of the sampling points;
and S4, calculating the area of the engineering area to be measured according to the ellipsoidal meridian curvature radius and the prime circle curvature radius.
2. The method for measuring an engineering surface based on ellipsoidal calculation according to claim 2, wherein the radius of curvature of the ellipsoidal meridian and the radius of curvature of the unitary mortise in the engineering area to be measured in step S3 satisfy the following formulas:
N=a(1-e2sin2(B2-B1))
wherein M is the radius of curvature of the meridian, N is the radius of curvature of the prime circle, a is the half-axis of the ellipsoid of the earth, e is the eccentricity of the engineering area to be measured, B2Upper boundary latitude, B, of the project area to be measured1The latitude of the lower boundary of the engineering area to be measured.
3. The method of claim 2, wherein the area of the engineering region to be measured in step S4 satisfies the following formula:
wherein P is the area of the engineering region to be measured, L2Is the right boundary longitude, L of the project area to be measured1The left boundary longitude of the project area to be tested.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010567367.7A CN111721272A (en) | 2020-06-19 | 2020-06-19 | Engineering surface measurement method based on ellipsoid calculation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010567367.7A CN111721272A (en) | 2020-06-19 | 2020-06-19 | Engineering surface measurement method based on ellipsoid calculation |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111721272A true CN111721272A (en) | 2020-09-29 |
Family
ID=72568140
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010567367.7A Pending CN111721272A (en) | 2020-06-19 | 2020-06-19 | Engineering surface measurement method based on ellipsoid calculation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111721272A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115098830A (en) * | 2022-06-28 | 2022-09-23 | 中国人民解放军战略支援部队信息工程大学 | Method and device for calculating difference between curvature radii of leveling surface and leveling ellipsoid |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101598786A (en) * | 2009-06-30 | 2009-12-09 | 武汉大学 | The rational function imaging model generation method of two/three grades of products of satellite image |
CN102288158A (en) * | 2011-08-02 | 2011-12-21 | 甘肃铁道综合工程勘察院有限公司 | Method for controlling Gaussian projection deformation based on normal section meridian ellipsoid |
CN102945570A (en) * | 2012-11-23 | 2013-02-27 | 华东师范大学 | Method for constructing full-space three-dimensional digital earth model |
CN106227237A (en) * | 2016-09-29 | 2016-12-14 | 广州极飞科技有限公司 | The distribution method of the aerial mission of unmanned plane and device |
CN106845035A (en) * | 2017-03-15 | 2017-06-13 | 武汉大学 | Sphere Terrain Niche influences tight method and system indirectly during the deviation of plumb line determines |
CN108088408A (en) * | 2016-11-23 | 2018-05-29 | 北京地林远景信息技术有限公司 | The method that a kind of smart mobile phone/tablet measures national territorial area |
-
2020
- 2020-06-19 CN CN202010567367.7A patent/CN111721272A/en active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101598786A (en) * | 2009-06-30 | 2009-12-09 | 武汉大学 | The rational function imaging model generation method of two/three grades of products of satellite image |
CN102288158A (en) * | 2011-08-02 | 2011-12-21 | 甘肃铁道综合工程勘察院有限公司 | Method for controlling Gaussian projection deformation based on normal section meridian ellipsoid |
CN102945570A (en) * | 2012-11-23 | 2013-02-27 | 华东师范大学 | Method for constructing full-space three-dimensional digital earth model |
CN106227237A (en) * | 2016-09-29 | 2016-12-14 | 广州极飞科技有限公司 | The distribution method of the aerial mission of unmanned plane and device |
CN108088408A (en) * | 2016-11-23 | 2018-05-29 | 北京地林远景信息技术有限公司 | The method that a kind of smart mobile phone/tablet measures national territorial area |
CN106845035A (en) * | 2017-03-15 | 2017-06-13 | 武汉大学 | Sphere Terrain Niche influences tight method and system indirectly during the deviation of plumb line determines |
Non-Patent Citations (2)
Title |
---|
吕朋一: "关于椭球区域面积计算问题的讨论", 《科技创新与应用》 * |
姜晨光: "我国陆地国土面积的初步测算", 《黑龙江测绘》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115098830A (en) * | 2022-06-28 | 2022-09-23 | 中国人民解放军战略支援部队信息工程大学 | Method and device for calculating difference between curvature radii of leveling surface and leveling ellipsoid |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106885576B (en) | AUV (autonomous Underwater vehicle) track deviation estimation method based on multipoint terrain matching positioning | |
CN103217688B (en) | Airborne laser radar point cloud adjustment computing method based on triangular irregular network | |
CN107504974B (en) | Terrain matching positioning method based on weighting of terrain blocks and terrain measuring points | |
CN108763825B (en) | Numerical simulation method for simulating wind field of complex terrain | |
CN107479078A (en) | Geodetic coordinates is converted to the method and system of separate planes coordinate in railroad survey | |
CN109141266B (en) | Steel structure measuring method and system | |
CN107167119B (en) | Data processing method for projection deformation | |
CN112233233B (en) | Three-dimensional linear recovery method for tunnel relative point cloud data | |
CN113378471B (en) | GPS elevation fitting method and system considering gravity terrain correction | |
CN112729258B (en) | River flow continuous measurement method based on satellite big data | |
CN112711027B (en) | Tunnel inner transverse positioning method based on laser radar point cloud data | |
CN114755661A (en) | Parameter calibration method and device for mobile laser scanning system | |
CN111854699A (en) | Unmanned aerial vehicle-based monitoring method for aerial survey river channel bank collapse process | |
CN108919319A (en) | Sea island reef satellite image Pillarless caving localization method and system | |
Du et al. | Cross‐section positioning based on a dynamic MLS tunnel monitoring system | |
CN112857315B (en) | Steel pipe column verticality measuring method based on three-dimensional laser scanning | |
CN108427741A (en) | A kind of DEM relative error evaluation methods based on a large amount of high-precision control points | |
CN113362468B (en) | Dimension measuring method for hub of train wheel | |
CN111721272A (en) | Engineering surface measurement method based on ellipsoid calculation | |
CN105204047A (en) | Detection and repair method for single gross error of observed quantity in satellite navigation system | |
CN110108287A (en) | A kind of unmanned vehicle high-precision map-matching method and system based on street lamp auxiliary | |
CN116412839A (en) | Amphibious integrated point cloud ranging system and calibration method | |
CN107727096B (en) | AUV terrain matching positioning method based on effective node screening | |
CN111765868A (en) | Earth measurement method based on oblique photography technology and divided according to different grids | |
Gao et al. | Accuracy comparison and analysis of interpolation methods in DEM generation with 3D laser point cloud data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200929 |