CN113587900A - Geographic altitude difference algorithm and system - Google Patents

Geographic altitude difference algorithm and system Download PDF

Info

Publication number
CN113587900A
CN113587900A CN202110888874.5A CN202110888874A CN113587900A CN 113587900 A CN113587900 A CN 113587900A CN 202110888874 A CN202110888874 A CN 202110888874A CN 113587900 A CN113587900 A CN 113587900A
Authority
CN
China
Prior art keywords
elevation
measurement point
point
measuring point
abnormal
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.)
Granted
Application number
CN202110888874.5A
Other languages
Chinese (zh)
Other versions
CN113587900B (en
Inventor
曲维荣
傅文祥
赵君毅
张永生
徐玉章
李纯洁
张涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qingdao Jielida Geographic Information Group Co ltd
Original Assignee
Qingdao Jielida Geographic Information Group Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Qingdao Jielida Geographic Information Group Co ltd filed Critical Qingdao Jielida Geographic Information Group Co ltd
Priority to CN202110888874.5A priority Critical patent/CN113587900B/en
Publication of CN113587900A publication Critical patent/CN113587900A/en
Application granted granted Critical
Publication of CN113587900B publication Critical patent/CN113587900B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C5/00—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Remote Sensing (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Algebra (AREA)
  • Evolutionary Biology (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)

Abstract

According to the geographic altitude difference algorithm and the geographic altitude difference system, firstly, the elevation abnormal change rate alpha in the direction from the fourth measuring point to the second measuring point is calculatedGII04‑GIIO2And the abnormal elevation change rate alpha from the fifth measuring point to the sixth measuring pointGII05‑GIIO6Then according to the abnormal change rate alpha of the elevationGII04‑GIIO2And the elevation anomaly rate of change alphaGII05‑GIIO6Calculating the abnormal elevation change rate alpha from the second measuring point to the fifth measuring pointGII02‑GIIO5. Finally, according to the difference of the geodetic heights of the second measuring point and the fifth measuring point and the abnormal elevation change rate alphaGII02‑GIIO5And the distance between the second measurement point and the fifth measurement point, and obtaining the difference in height Δ H between the second measurement point and the fifth measurement point6. The calculation method is simple, the calculation precision is high, and the work difficulty is greatly reduced.

Description

Geographic altitude difference algorithm and system
Technical Field
The present application relates to the field of geographic mapping, and in particular, to a geographic height difference algorithm and system.
Background
In the field of geographical mapping, the geographical height difference is an important geographical parameter. The geographical height difference is an important data index for construction and geographical drawing. However, the existing surveying method has poor precision, which affects the surveying quality of the geographic height.
Disclosure of Invention
In view of the above, it is necessary to provide a geographic height difference algorithm and system for solving the above technical problems.
A geographic elevation difference algorithm, sequentially distributing a fourth measurement point, a second measurement point, a fifth measurement point and a sixth measurement point along a first direction, the elevation difference algorithm being configured to measure an elevation difference between the second measurement point and the fifth measurement point, the method comprising:
calculating the abnormal elevation change rate alpha from the fourth measuring point to the second measuring pointGII04-GIIO2And the abnormal elevation change rate alpha from the fifth measuring point to the sixth measuring pointGII05-GIIO6；
According to the abnormal change rate alpha of the elevationGII04-GIIO2And the elevation anomaly rate of change alphaGII05-GIIO6Calculating the abnormal elevation change rate alpha from the second measuring point to the fifth measuring pointGII02-GIIO5；
According to the difference of the geodetic heights of the second measuring point and the fifth measuring point and the abnormal elevation change rate alphaGII02-GIIO5And the distance between the second measurement point and the fifth measurement point, and obtaining the difference in height Δ H between the second measurement point and the fifth measurement point6。
In one embodiment, the calculation of the abnormal elevation change rate α from the fourth measurement point to the second measurement point is performedGII04-GIIO2And in the direction from the fifth measuring point to the sixth measuring pointElevation anomaly rate of change alphaGII05-GIIO6Includes calculating the abnormal elevation change rate alpha by the following formulaGII04-GIIO2：
αGII04-GIIO2＝(⊿H1-⊿H2)/SGII04-GIIO2；
Wherein S isGII04-GIIO2For a pitch between said fourth measurement point and said second measurement point, Δ H1 being a geodetic height difference, Δ H, between said fourth measurement point and said second measurement point2Is the normal height difference from the fourth measuring point to the second measuring point.
In one embodiment, the calculation of the abnormal elevation change rate α from the fourth measurement point to the second measurement point is performedGII04-GIIO2And the abnormal elevation change rate alpha from the fifth measuring point to the sixth measuring pointGII05-GIIO6Includes calculating the abnormal elevation change rate alphaGII05-GIIO6：
αGII05-GIIO6＝(⊿H3-⊿H4)/SGII05-GIIO6；
Wherein S isGII05-GIIO6A pitch, Δ H, from said fifth measurement point to said sixth measurement point3A difference in geodetic height, Δ H, between said fifth and sixth measurement points4Is the normal height difference from the fifth measuring point to the sixth measuring point.
In one embodiment, the abnormal rate of change a according to the elevation isGII04-GIIO2And the elevation anomaly rate of change alphaGII05-GIIO6Calculating the abnormal elevation change rate alpha from the second measuring point to the fifth measuring pointGII02-GIIO5The method comprises the following steps:
for the abnormal change rate alpha of the elevationGII04-GIIO2And the elevation anomaly rate of change alphaGII05-GIIO6Averaging to obtain the abnormal elevation change rate alpha from the second measuring point to the fifth measuring pointGII02-GIIO5。
In one embodiment, the altitude abnormal change rate α is determined according to the difference between the geodetic heights of the second measurement point and the fifth measurement pointGII02-GIIO5And the distance between the second measurement point and the fifth measurement point, and obtaining the difference in height Δ H between the second measurement point and the fifth measurement point6In the first measurement point, a height difference Δ H between the second measurement point and the fifth measurement point6Calculated by the following method:
⊿H6＝⊿H5-αGII02-GIIO5﹡SGII02-GIIO5；
wherein (delta H)6The normal height difference from the second measuring point to the fifth measuring point;
⊿H5the geodetic height difference from the second measuring point to the fifth measuring point;
αGII02-GIIO5the elevation abnormal change rate from the second measuring point to the fifth measuring point is obtained;
SGII02-GIIO5is the straight distance from the second measuring point to the fifth measuring point.
In one embodiment, the elevation abnormal change rate alpha is obtained according to the difference of the earth heights of the second measuring point and the fifth measuring pointGII02-GIIO5And the distance between the second measurement point and the fifth measurement point, and obtaining the difference in height Δ H between the second measurement point and the fifth measurement point6And then, S40, the step of verifying by the distance measurement triangle elevation leveling measurement is also included.
In one embodiment, the elevation abnormal change rate alpha is obtained according to the difference of the earth heights of the second measuring point and the fifth measuring pointGII02-GIIO5And the distance between the second measurement point and the fifth measurement point, and obtaining the difference in height Δ H between the second measurement point and the fifth measurement point6And then, the step of verifying by a second-class leveling method is also included.
A geographic elevation difference algorithm system comprising:
a first calculating module, configured to calculate an elevation anomaly change rate α in a direction from the fourth measurement point to the second measurement pointGII04-GIIO2And the abnormal elevation change rate alpha from the fifth measuring point to the sixth measuring pointGII05-GIIO6；
Second calculationA module for determining the abnormal change rate alpha of the elevationGII04-GIIO2And the elevation anomaly rate of change alphaGII05-GIIO6Calculating the abnormal elevation change rate alpha from the second measuring point to the fifth measuring pointGII02-GIIO；
A height difference module for calculating the height difference between the second measurement point and the fifth measurement point and the abnormal elevation change rate alphaGII02-GIIO5And the distance between the second measurement point and the fifth measurement point, and obtaining the difference in height Δ H between the second measurement point and the fifth measurement point6。
A computer device comprising a memory and a processor, the memory storing a computer program, characterized in that the processor implements the steps of the method when executing the computer program.
A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the steps of the method.
In the geographic altitude difference algorithm provided by the embodiment of the application, the elevation abnormal change rate alpha in the direction from the fourth measuring point to the second measuring point is calculatedGII04-GIIO2And the abnormal elevation change rate alpha from the fifth measuring point to the sixth measuring pointGII05-GIIO6Then according to the abnormal change rate alpha of the elevationGII04-GIIO2And the elevation anomaly rate of change alphaGII05-GIIO6Calculating the abnormal elevation change rate alpha from the second measuring point to the fifth measuring pointGII02-GIIO5. Finally, according to the difference of the geodetic heights of the second measuring point and the fifth measuring point and the abnormal elevation change rate alphaGII02-GIIO5And the distance between the second measurement point and the fifth measurement point, and obtaining the difference in height Δ H between the second measurement point and the fifth measurement point6. The calculation method is simple, the calculation precision is high, and the work difficulty is greatly reduced.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present application, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
Fig. 1 is a flowchart of a geographic altitude difference algorithm provided in an embodiment of the present application;
FIG. 2 is a schematic diagram of a geographic elevation difference algorithm location provided by an embodiment of the present application;
fig. 3 is a schematic location diagram of a geographic altitude difference algorithm provided in an embodiment of the present application.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application.
Referring to fig. 1-2, an embodiment of the present application provides a geographic height difference algorithm, where a fourth measurement point, a second measurement point, a fifth measurement point, and a sixth measurement point are sequentially distributed along a first direction, and the height difference algorithm is configured to measure a height difference between the second measurement point and the fifth measurement point. The method comprises the following steps:
s10, calculating the abnormal elevation change rate alpha from the fourth measuring point to the second measuring pointGII04-GIIO2And the abnormal elevation change rate alpha from the fifth measuring point to the sixth measuring pointGII05-GIIO6；
S20, according to the abnormal elevation change rate alphaGII04-GIIO2And the elevation anomaly rate of change alphaGII05-GIIO6Calculating the abnormal elevation change rate alpha from the second measuring point to the fifth measuring pointGII02-GIIO5；
S30, according to the difference between the ground height of the second measuring point and the ground height of the fifth measuring point and the abnormal elevation change rate alphaGII02-GIIO5And the distance between the second measurement point and the fifth measurement point, and obtaining the difference in height Δ H between the second measurement point and the fifth measurement point6。
In this embodiment, the fourth measurement point gii 04, the second measurement point gii 02, the fifth measurement point gii 05, and the sixth measurement point gii 06 may be approximately distributed along a certain direction. Among them, there may be obstacles such as a river between the second measurement point gii 02 and the fifth measurement point gii 05. The height difference between the second and fifth measuring points gii 02, gii 05 can thus be measured by calculation. Therefore, the method provided by the embodiment can be widely applied to altitude difference measurement in the fields of bridge crossing, sea crossing and the like.
In the geographic altitude difference algorithm provided by the embodiment of the application, the elevation abnormal change rate alpha in the direction from the fourth measuring point to the second measuring point is calculatedGII04-GIIO2And the abnormal elevation change rate alpha from the fifth measuring point to the sixth measuring pointGII05-GIIO6Then according to the abnormal change rate alpha of the elevationGII04-GIIO2And the elevation anomaly rate of change alphaGII05-GIIO6Calculating the abnormal elevation change rate alpha from the second measuring point to the fifth measuring pointGII02-GIIO5. Finally, according to the difference of the geodetic heights of the second measuring point and the fifth measuring point and the abnormal elevation change rate alphaGII02-GIIO5And the distance between the second measurement point and the fifth measurement point, and obtaining the difference in height Δ H between the second measurement point and the fifth measurement point6. The calculation method is simple, the calculation precision is high, and the work difficulty is greatly reduced.
Referring to fig. 2, a fourth measuring point and a second measuring point may be set on the west bank of a river. The fifth and sixth measurement points are on the east bank of the river. The west fourth measurement point gii 04 to the second measurement point gii 02 and the east fifth measurement point gii 05 to the sixth measurement point gii 06 are as long as possible.
In one embodiment, the distance from the west land fourth measurement point gii 04 to the second measurement point gii 02 and the east land fifth measurement point gii 05 to the sixth measurement point gii 06 is preferably as long as the width distance of the river.
In one embodiment, the west bank fourth measurement point gii 04, the second measurement point gii 02, the east bank fifth measurement point gii 05, and the sixth measurement point gii 06 are on a straight line.
In one embodiment, the step S10 includes calculating the abnormal elevation change rate α according to the following formulaGII04-GIIO2：
αGII04-GIIO2＝(⊿H1-⊿H2)/SGII04-GIIO2；
Wherein S isGII04-GIIO2For a pitch between said fourth measurement point and said second measurement point, Δ H1 being a geodetic height difference, Δ H, between said fourth measurement point and said second measurement point2Is the normal height difference from the fourth measuring point to the second measuring point.
In one embodiment, the step S10 includes calculating the abnormal change rate α of elevation byGII05-GIIO6：
αGII05-GIIO6＝(⊿H3-⊿H4)/SGII05-GIIO6；
Wherein S isGII05-GIIO6A pitch, Δ H, from said fifth measurement point to said sixth measurement point3A difference in geodetic height, Δ H, between said fifth and sixth measurement points4Is the normal height difference from the fifth measuring point to the sixth measuring point.
That is, αGII04-GIIO2、αGII05-GIIO6The elevation abnormal change rates are respectively the elevation abnormal change rates in the directions from GII 04 to GII 02 and from GII 05 to GII 06, and the unit is m/km.
SGII04-GIIO2、SGII05-GIIO6The respective distances in km for GII 04 to GII 02 and GII 05 to GII 06.
⊿H1、⊿H3The large height differences of GII 04 to GII 02 and GII 05 to GII 06, respectively, are given in m.
⊿H2、⊿H4The normal height differences in m for GII 04 to GII 02 and GII 05 to GII 06, respectively.
Respectively calculating alpha values of two banks of the bridge according to the formula, and finally taking the average value of alpha values of the two banks as the elevation abnormal change rate alpha from GII 02 to GII 05GII02-GIIO5。
In one embodiment, the S20 includes:
for the abnormal change rate alpha of the elevationGII04-GIIO2And the elevation anomaly rate of change alphaGII05-GIIO6Averaging to obtain the abnormal elevation change rate alpha from the second measuring point to the fifth measuring pointGII02-GIIO5。
In one embodiment, in S30, the difference in height Δ H between the second measurement point and the fifth measurement point6Calculated by the following method:
⊿H6＝⊿H5-αGII02-GIIO5﹡SGII02-GIIO5；
wherein (delta H)6The normal height difference from the second measuring point to the fifth measuring point;
⊿H5the geodetic height difference from the second measuring point to the fifth measuring point;
αGII02-GIIO5the elevation abnormal change rate from the second measuring point to the fifth measuring point is obtained;
SGII02-GIIO5is the straight distance from the second measuring point to the fifth measuring point.
That is, Δ H6The normal height difference from GII 02 to GII 05; delta H5The ground height difference is from GII 02 to GII 05; alpha (a)GII02-GIIO5The elevation abnormal change rate from GII 02 to GII 05; sGII02-GIIO5The distances from GII 02 to GII 05.
In one embodiment, after the step of S30, the method further includes a step of S40 of performing verification through ranging triangulation elevation leveling.
Referring to fig. 3, the present embodiment is illustrated by using cross-sea measurement points: namely, the sea-crossing leveling measurement is carried out by adopting the distance measuring triangle elevation.
Firstly 1, selecting and laying a site, and 1) carrying out sea crossing leveling and surveying four times due to the long length of the bridge.
2) And (3) laying observation graphs according to the established priority pier layout, wherein the cross-sea graph layout is shown in figure 3.
In the figure, KH1 and KH2 are shore points on 118# priority piers (bridgeheads), KH3, KH4, KH5 and KH6 are points on 110# and 80# priority piers respectively, and KH7, KH8, KH9 and KH10 are points on 47# and 46# pier working platforms respectively; the GII 02 (IIwater 02) and the GII 05 (IIwater 07) are primary and secondary leveling points (also secondary GPS points).
Because of the influence of the size of the bridge pier (the size of the top of the bridge pier is about 7m × 3m), the lengths between KH3 and KH4 and between KH5 and KH6 can only be set to be 6.1m to 6.2m, and therefore the height near-shore points of the sea-crossing distance measuring triangle are uniformly set to be 6.1m to 6.2 m.
2. Setting of auxiliary points
KH1, KH2, KH3, KH4, KH5 and KH6 are all arranged on the preferential pier, a round hole with the diameter of 1.8cm is drilled by impact firstly, then a twisted steel bar with the diameter of 2cm is driven into the round hole, and then the twisted steel bar is reinforced by concrete, and the top of the twisted steel bar is carved with a cross. KH7, KH8 and KH10 are arranged on the working platform, and the setting method of the point positions is as follows: a reinforcing steel bar with the diameter of about 2cm and the length of about 10cm is welded on a main cross beam of the platform, and the top of the reinforcing steel bar is carved with a cross. KH9 utilizes the three equal GPS points (north beacon piers). The auxiliary point is stable and reliable.
3 Cross-sea level observation
1) Two NET05 total stations (angle measurement accuracy 0.5 second, distance measurement accuracy 0.8mm +1ppm d 10) were used for cross-sea level observation-6) And (6) carrying out observation.
2) The instrument is placed under the open-air shadow 30 minutes before observation, so that the temperature of the instrument is consistent with the outside temperature, and the umbrella is used for shielding the sunlight during observation.
3) After the observation of one test-loop is finished, the interval is 15-20 minutes, and then the observation of the next test-loop is carried out.
4) The elevation difference between two points on the shore is observed back and forth by using a second-class leveling method.
5) The target board on the opposite bank is manufactured according to a diagram C.2 in appendix C1 of national first and second-class leveling regulations.
4. Measured return number and limit difference
1) The measured data and the group data of the distance-measuring triangulation elevation measurement are shown in Table 4
TABLE 4 observed number of test returns and number of groups
Figure BDA0003195122480000091
2) The difference between the heights of the two measured loops is smaller than the limit difference calculated according to the following formula.
Figure BDA0003195122480000092
In the formula: n-the number of double-survey, and S-the sea-crossing sight length, with km as a unit.
5 Observation method
1) Measurement of height difference between the present bank survey stations
The ZDL700 electronic level is used, KH 1-KH 2, KH 3-KH 4, KH 5-KH 6, KH 7-KH 8 and KH 9-KH 10 are respectively used as a measuring station, the backward measurement is carried out according to the requirements of the second-class leveling, the backward measurement is 0.5mm at the maximum, and 0.6mm is allowed. And taking the number of round trip measurements as a formal result of the height difference of the measuring station points, and taking the result as a reference for detecting and calculating the height of the measuring station instrument.
2) Distance measurement
The measurement of the distance between the shore and the sea is carried out by using a NET05-10302 total station (the distance measurement precision is 0.8mm +1ppm d 10)-6) And 4, measuring back by using a round-trip observation plain gauge. The temperature is read to 0.2 ℃, the air pressure is read to 50 Pa, and the temperature is directly input into a total station to be automatically corrected by an instrument. The reading range of one measured distance is less than or equal to 10mm, the distance between measured distances is less than or equal to 15mm, the distance between the measured distances is 1.0mm after various corrections, the maximum distance is 1.54mm, and the allowable distance is 1.60 mm.
3) Vertical angle observation
The observation procedure and observation method of the vertical angle were performed according to the national first, second class leveling Specification, page 23, items 8.9.4.1 and 8.9.4.2.
The error of each group of vertical angle observation is within the tolerance range.
6. Side length correction and checking
1) Side length correction calculation
The length of the observation side is projected to a reference ellipsoid after being corrected by an addition and multiplication constant of an instrument and is reduced to a Gaussian plane, and each correction calculation formula is as follows:
correcting addition and multiplication constants of the instrument:
D＝D′*k+d
secondly, projection correction:
S′0＝D{1-(Hm+hm)/R+(Hm+hm)2/R2}；
③ reducing and correcting:
S0＝S′0{1+Ym2/(2R2)+(⊿Y)2/(24R2)}。
in the formula: d', observing the flat pitch;
k: instrument multiplier constant (k as 0.99999776)
d: instrument constant (d ═ 0.00058)
Hm: average elevation of two end points of the distance measuring edge;
hm: height of geodesic to reference ellipsoid (hm 67 m);
ym: measuring the average natural value of the horizontal coordinates of the edges;
r: average radius of curvature (R-6366255 m) of the measurement area.
2) Checking calculation
The maximum difference between the heights of the two double measured returns is 8.03mm, and the allowable range is +/-10.34 mm; three independent closed rings are formed by the earth quadrangle, and the height difference of each side in the same time period is used for calculating the closing difference, wherein the maximum value is 3mm and 11.1mm, and the allowance is +/-21.9 mm. The accuracy statistics of each span are shown in table 1.
TABLE 1 accuracy statistics for each span-range triangulation
Figure BDA0003195122480000111
The distance measurement triangle elevation cross-sea leveling measurement is finally checked and calculated, 4 spans are carried out simultaneously, the height difference number calculated by each measuring loop is adopted, the height difference number is calculated on a microcomputer by using 'Qinghuashan Wei NASEW95 engineering measurement control net adjustment software', the full-mean error is out-of-limit (the full-mean error is +/-1.1 mm, and is allowed +/-1.0 mm) due to the fact that the observation height difference error from KH5 to KH6 is large (the ring closure difference is 4.3mm and is allowed to be 6.6mm) in the checking and calculating process, and the full-mean error is removed in the checking and calculating process. After eliminating the height difference from KH5 to KH6, the maximum ring closure difference is-2.3 mm, and the allowance is + -6.7 mm. The error of the whole is +/-0.96 mm.
In one embodiment, the method further comprises step S50, verifying by the second-class leveling method:
the height difference between a connecting distance-measuring triangle elevation cross-sea leveling starting point and a leveling point (second-class GPS point) is measured by a ZDL700 electronic level, and the joint measuring points are KH1 to GII 05 (II water 05) and KH9 to GII 02 (II water 02) and respectively form two closed rings.
1. Field observation
1) The second-class leveling measurement is observed by using a ZDL700(DS03 type) electronic level and a bar code indium tile leveling staff, and one-time reciprocating observation is adopted by adopting a single line. The operation process is strictly carried out according to the relevant requirements of national first-class and second-class leveling regulations. Before observation, the observation route number and the test section sequence number are arranged according to the net layout sketch map in the design.
2) During observation, the temperature is low, the sunlight is mild, and the sunlight is often cloudy, the sunlight is not very strong around noon and the temperature difference change is not very obvious, so that the scale is clearly divided, and the imaging is stable. The forward measurement and the backward measurement of the same measurement section are generally performed symmetrically in the morning or afternoon, respectively. The total station number of the asymmetric observation does not exceed 30 percent of the total station number of the observation,
3) the visual distance between the front and the back of the survey station is generally controlled to be 5m to 10m on the trestle, the maximum is 15m, and the visual distance on the land is generally controlled to be about 30 m. The visual distance is controlled by adopting the pulling force of the measuring rope, the visual distance difference of a single station is controlled to be 1.50 meters, and the accumulated difference of the visual distances of the front and the back of any measuring station is less than 6 meters.
4) During observation, the position correctness of the collimation axis is checked once every day (i-angle inspection), and automatic correction is carried out by using a self-contained program of the instrument.
5) Leveling observations are made with the instrument and staff as in-line as possible. The observation adopts a 5 kg ruler platform, and the ruler platform is stabilized by a foot-operated pile before the ruler is not placed.
6) The observation sequence and method of the observation station adopt the own program of the instrument, and meet the requirements of the specification and the design book.
7) The forward test line or the backward test line are even stations.
2. Recording
1) The observation data is automatically stored in the memory of the electronic level, and the automatically recorded original data is transmitted to a microcomputer and automatically formed into a field observation file to be output through a printer.
2) In the observation process, factors such as date, temperature, weather, wind power and the like are recorded according to requirements.
3. Checking calculation
The observation height difference of each section of the second-class leveling is corrected by the true length per meter of the leveling staff, the temperature of the staff and the unparallel correction of the normal leveling surface, and a leveling field height difference and an approximate height table are compiled. The formula for calculating the correction of each term is as follows:
the correction number of the height difference leveling staff per meter is calculated as follows: δ is f × h
The scale temperature correction calculation formula is as follows: α ═ T (T-T)0)×α×h
The normal horizontal plane non-parallelism correction is: e ═ y — (γ)i+1-γi)×Hm/γm
γm＝(γi+γi+1)/2-0.1543×Hm
γ＝978032×(1+0.005302sin2Φ-0.0000058sin22Φ)
In the formula: f is the true length per meter of the leveling staff (f is 0.006 mm/m); t is the temperature of the scale, (. degree.C.); t is0Temperature (T) for scale length identification020 ℃); alpha is the expansion coefficient of invar tape (alpha is 0.0015 mm/(m.deg.C); gamma isi、γi+1＝Normal gravity values on ellipsoids of point i and point i +1, respectively (10)-5m/s2) (ii) a Hm is the approximate height average (m) of two level points; phi is the horizontal point latitude; gammam is the average value of normal gravity of two level points (10)-5m/s2)。
The precision of the second-class level after each correction is shown in Table 2
TABLE 2 statistical table of the accuracy of the second-class level
Figure BDA0003195122480000131
Figure BDA0003195122480000141
The difference between the altitude difference from GII 02 (IIwater 02) to GII 05 (IIwater 05) obtained by the geographical altitude difference algorithm and the altitude difference from GII 02 (IIwater 02) to GII 05 (IIwater 05) obtained by the conventional method is 0.009 m, and the allowable difference between the altitude differences of the detected sections detected according to the second-class level detection is equal to
Figure BDA0003195122480000142
The allowable error of the calculation (L is the route length) is +/-0.014 m, which indicates that the altitude difference between GII 02 (IIwater 02) and GII 05 (IIwater 05) obtained by a geographical altitude difference algorithm completely meets the specification requirement.
The embodiment of the present application further provides a geographic height difference algorithm system, including:
a first calculating module, configured to calculate an elevation anomaly change rate α in a direction from the fourth measurement point to the second measurement pointGII04-GIIO2And the abnormal elevation change rate alpha from the fifth measuring point to the sixth measuring pointGII05-GIIO6；
A second calculation module for calculating the abnormal change rate alpha of the elevationGII04-GIIO2And the elevation anomaly rate of change alphaGII05-GIIO6Calculating the abnormal elevation change rate alpha from the second measuring point to the fifth measuring pointGII02-GIIO；
A height difference module for calculating the height difference between the second measurement point and the fifth measurement point and the abnormal elevation change rate alphaGII02-GIIO5And the distance between the second measurement point and the fifth measurement point, and obtaining the difference in height Δ H between the second measurement point and the fifth measurement point6。
The embodiment of the present application further provides a computer device, which includes a memory and a processor, where the memory stores a computer program, and the processor implements the steps of the method according to any one of the above embodiments when executing the computer program.
A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the steps of the method of any of the above embodiments.
Those skilled in the art will appreciate that the implementation of the modules such as the program and the terminal in the above embodiments can be accomplished by instructing related hardware through a computer program, where the computer program can be stored in a non-volatile computer-readable storage medium, and when executed, the computer program can include the flow of the above embodiments of the methods. Any reference to memory, storage, database or other medium used in the embodiments provided herein can include at least one of non-volatile and volatile memory. Non-volatile memory may include Read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical storage, or the like. Volatile memory can include Random Access Memory (RAM) or external cache memory. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM), among others.
The technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, but should be considered as the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present application, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent shall be subject to the appended claims.

Claims (10)

1. A geographic elevation difference algorithm, which sequentially distributes a fourth measurement point, a second measurement point, a fifth measurement point and a sixth measurement point along a first direction, the elevation difference algorithm being configured to measure an elevation difference between the second measurement point and the fifth measurement point, the method comprising:
calculating the abnormal elevation change rate alpha from the fourth measuring point to the second measuring pointGII04-GIIO2And the abnormal elevation change rate alpha from the fifth measuring point to the sixth measuring pointGII05-GIIO6；
According to the abnormal change rate alpha of the elevationGII04-GIIO2And the elevation anomaly rate of change alphaGII05-GIIO6Calculating the abnormal elevation change rate alpha from the second measuring point to the fifth measuring pointGII02-GIIO5；
According to the difference of the geodetic heights of the second measuring point and the fifth measuring point and the abnormal elevation change rate alphaGII02-GIIO5And the distance between the second measurement point and the fifth measurement point, and obtaining the difference in height Δ H between the second measurement point and the fifth measurement point6。
2. The geographic height difference algorithm of claim 1, wherein the calculating of the abnormal rate of change α in elevation from the fourth measurement point to the second measurement point is performed at a rate of change α in elevationGII04-GIIO2And the abnormal elevation change rate alpha from the fifth measuring point to the sixth measuring pointGII05-GIIO6Includes calculating the abnormal elevation change rate alpha by the following formulaGII04-GIIO2：
αGII04-GIIO2＝(⊿H1-⊿H2)/SGII04-GIIO2；
Wherein S isGII04-GIIO2For a pitch between said fourth measurement point and said second measurement point, Δ H1 being a geodetic height difference, Δ H, between said fourth measurement point and said second measurement point2Is the normal height difference from the fourth measuring point to the second measuring point.
3. The geographic height difference algorithm of claim 2, wherein the calculating of the abnormal rate of change α in elevation from the fourth measurement point to the second measurement point is performed at a rate of change α in elevationGII04-GIIO2And the abnormal elevation change rate alpha from the fifth measuring point to the sixth measuring pointGII05-GIIO6Comprises the steps ofCalculating the abnormal change rate alpha of the elevation in the following mannerGII05-GIIO6：
αGII05-GIIO6＝(⊿H3-⊿H4)/SGII05-GIIO6；
Wherein S isGII05-GIIO6A pitch, Δ H, from said fifth measurement point to said sixth measurement point3A difference in geodetic height, Δ H, between said fifth and sixth measurement points4Is the normal height difference from the fifth measuring point to the sixth measuring point.
4. The geographic height difference algorithm of claim 3, wherein the rate of change α is a function of the elevation anomaly rateGII04-GIIO2And the elevation anomaly rate of change alphaGII05-GIIO6Calculating the abnormal elevation change rate alpha from the second measuring point to the fifth measuring pointGII02-GIIO5The method comprises the following steps:
for the abnormal change rate alpha of the elevationGII04-GIIO2And the elevation anomaly rate of change alphaGII05-GIIO6Averaging to obtain the abnormal elevation change rate alpha from the second measuring point to the fifth measuring pointGII02-GIIO5。
5. The geographic elevation difference algorithm of claim 4 wherein the elevation anomaly rate of change α is based on the geodetic elevation difference between the second and fifth measurement pointsGII02-GIIO5And the distance between the second measurement point and the fifth measurement point, and obtaining the difference in height Δ H between the second measurement point and the fifth measurement point6In the first measurement point, a height difference Δ H between the second measurement point and the fifth measurement point6Calculated by the following method:
⊿H6＝⊿H5-αGII02-GIIO5﹡SGII02-GIIO5；
wherein (delta H)6The normal height difference from the second measuring point to the fifth measuring point;
⊿H5the geodetic height difference from the second measuring point to the fifth measuring point;
αGII02-GIIO5is the second measurement point toThe elevation abnormal change rate of a fifth measuring point;
SGII02-GIIO5is the straight distance from the second measuring point to the fifth measuring point.
6. The geographic elevation difference algorithm of claim 1 wherein the elevation anomaly rate of change α is based on the geodetic elevation difference between the second measurement point and the fifth measurement pointGII02-GIIO5And the distance between the second measurement point and the fifth measurement point, and obtaining the difference in height Δ H between the second measurement point and the fifth measurement point6And then, S40, the step of verifying by the distance measurement triangle elevation leveling measurement is also included.
7. The geographic elevation difference algorithm of claim 1 wherein the elevation anomaly rate of change α is based on the geodetic elevation difference between the second measurement point and the fifth measurement pointGII02-GIIO5And the distance between the second measurement point and the fifth measurement point, and obtaining the difference in height Δ H between the second measurement point and the fifth measurement point6And then, the step of verifying by a second-class leveling method is also included.
8. A geographic elevation difference algorithm system, comprising:
a first calculating module, configured to calculate an elevation anomaly change rate α in a direction from the fourth measurement point to the second measurement pointGII04-GIIO2And the abnormal elevation change rate alpha from the fifth measuring point to the sixth measuring pointGII05-GIIO6；
A second calculation module for calculating the abnormal change rate alpha of the elevationGII04-GIIO2And the elevation anomaly rate of change alphaGII05-GIIO6Calculating the abnormal elevation change rate alpha from the second measuring point to the fifth measuring pointGII02-GIIO；
A height difference module for calculating the height difference between the second measurement point and the fifth measurement point and the abnormal elevation change rate alphaGII02-GIIO5And the second measuring point and theThe straight distance of the fifth measurement point is obtained, and the height difference delta H between the second measurement point and the fifth measurement point is obtained6。
9. A computer device comprising a memory and a processor, the memory storing a computer program, wherein the processor implements the steps of the method of any one of claims 1 to 7 when executing the computer program.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method of any one of claims 1 to 7.
CN202110888874.5A 2021-08-04 2021-08-04 Geographic altitude difference algorithm and system Active CN113587900B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110888874.5A CN113587900B (en) 2021-08-04 2021-08-04 Geographic altitude difference algorithm and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110888874.5A CN113587900B (en) 2021-08-04 2021-08-04 Geographic altitude difference algorithm and system

Publications (2)

Publication Number Publication Date
CN113587900A true CN113587900A (en) 2021-11-02
CN113587900B CN113587900B (en) 2023-06-06

Family

ID=78254672

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110888874.5A Active CN113587900B (en) 2021-08-04 2021-08-04 Geographic altitude difference algorithm and system

Country Status (1)

Country Link
CN (1) CN113587900B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102175213A (en) * 2011-01-30 2011-09-07 中国测绘科学研究院 GPS leveling height transfer method
US20130204567A1 (en) * 2012-02-03 2013-08-08 Suunto Oy Method and arrangement for determining altitude
WO2016050871A1 (en) * 2014-10-03 2016-04-07 Universite Pierre Et Marie Curie (Paris 6) Method for characterizing the cracking mechanism of a material from the fracture surface thereof
CN106989717A (en) * 2017-03-28 2017-07-28 湖北省测绘工程院 A kind of quasigeoid detection method and device
CN108761510A (en) * 2018-04-09 2018-11-06 中铁二局集团有限公司 The method for carrying out height differences by leveling measurement using the gravity field model based on topographic correction
CN108917706A (en) * 2018-06-26 2018-11-30 中铁大桥勘测设计院集团有限公司 A kind of GNSS leveling measuring method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102175213A (en) * 2011-01-30 2011-09-07 中国测绘科学研究院 GPS leveling height transfer method
US20130204567A1 (en) * 2012-02-03 2013-08-08 Suunto Oy Method and arrangement for determining altitude
WO2016050871A1 (en) * 2014-10-03 2016-04-07 Universite Pierre Et Marie Curie (Paris 6) Method for characterizing the cracking mechanism of a material from the fracture surface thereof
CN106989717A (en) * 2017-03-28 2017-07-28 湖北省测绘工程院 A kind of quasigeoid detection method and device
CN108761510A (en) * 2018-04-09 2018-11-06 中铁二局集团有限公司 The method for carrying out height differences by leveling measurement using the gravity field model based on topographic correction
CN108917706A (en) * 2018-06-26 2018-11-30 中铁大桥勘测设计院集团有限公司 A kind of GNSS leveling measuring method

Cited By (1)

* Cited by examiner, † Cited by third party
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

Also Published As

Publication number Publication date
CN113587900B (en) 2023-06-06

Similar Documents

Publication Publication Date Title
CN102607506A (en) Free stationing transformation monitoring method of high-fill airport side slope unit set total station
CN109470222B (en) Supervision control method for super high-rise building engineering measurement
CN106840210B (en) High-precision clinometer calibration method
CN105388494A (en) Laser ranging positioning method for RTK receiver
CN114046770A (en) Marine settlement observation and measurement method
CN111609833B (en) Settlement observation method for high-rise building
CN113587900B (en) Geographic altitude difference algorithm and system
CN108917706A (en) A kind of GNSS leveling measuring method
Abduraufovich SIMPLE AND ACCURATE METHODS OF SYOMKAS PERFORMED IN THE FORMATION OF THE EARTH
CN112964237B (en) Measurement control system, method and device for construction engineering and computer equipment
Pirti et al. Evaluation of some levelling techniques in surveying application
Banasik et al. The use of quasigeoid in leveling through terrain obstacles
Chukwuocha et al. Solving control reference azimuth problems of traversing using reorientation traversing
Shin et al. A study on airborne LiDAR calibration and operation techniques for bathymetric survey
Idoko et al. Comparison of Orthometric Heights Obtained Using Total Station and Differential Global Positioning Systems (DGPS) with Precise Levels Instruments
Bikbulatova et al. Geodetic control in construction of steel roof structures
CN103954302A (en) Flat track design method of rocket sled test
CN108955635A (en) Pair observations unit composite level bilateral observation leveling measuring method
Lehmuskoski Precise levelling of the Olkiluoto GPS Network in 2005
Plopeanu et al. SPECIAL GEODETIC MONITORING OF ROCKFILL DAMS-RÂUŞOR DAM CASE
Підгурський Methodical instructions for “Engineering geodesy” course work
Bryś et al. Precise integrated hybrid networks in 21st-century tunnel surveying
Zrinjski et al. Geodetic Basis of the Longest Tunnel in the Republic of Croatia
Audu et al. Comparative Assessment of the Accuracy of the Elevation differences obtained from different Geomatics Techniques and Instruments
Umar et al. Determination of Vertical Movement for Deformation Monitoring of a Bridge Using GPS Technique

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