CN113536559B - Wire data processing system and method in hydraulic engineering measurement - Google Patents
Wire data processing system and method in hydraulic engineering measurement Download PDFInfo
- Publication number
- CN113536559B CN113536559B CN202110762742.8A CN202110762742A CN113536559B CN 113536559 B CN113536559 B CN 113536559B CN 202110762742 A CN202110762742 A CN 202110762742A CN 113536559 B CN113536559 B CN 113536559B
- Authority
- CN
- China
- Prior art keywords
- point
- wire
- measurement
- station
- measurement data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Arrangements For Transmission Of Measured Signals (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
The invention discloses a wire data processing system and method in hydraulic engineering measurement, which consists of a server side and a plurality of user sides; the server side is used for carrying out data processing on the measurement data uploaded by each user side and feeding back the measurement data processing result to the corresponding user side in real time; the field measurement personnel guide the measurement data into a remote server through a user terminal, and the measurement data is processed, stored and fed back in real time through the server, so that the field personnel can obtain adjustment results, precision assessment and lead results in time, the problem of finding the measurement data in time is solved, the problem is solved in real time, reworking is avoided, the labor intensity and the cost are further reduced, and the lead layout speed and the working efficiency are improved; and the measurement data are prevented from being tidied again and the measurement result is prevented from being saved in the room.
Description
Technical Field
The invention relates to the field of field measurement in hydraulic engineering, in particular to a lead data processing system and method in hydraulic engineering measurement.
Background
The national surveying and mapping geographic information bureau is popularizing and applying unmanned aerial vehicle aerial photography technology for improving the surveying and mapping technical equipment in China. The unmanned aerial vehicle aerial photography technology plays a positive role in constructing digital water conservancy, detecting geographical national conditions, improving water conservancy management efficiency and the like. In hydraulic engineering measurement, a wire network needs to be laid before unmanned aerial vehicle aerial photography begins, so that the wire network layout needs to be realized rapidly on a mapping site, and current field operators need to return to the room to perform adjustment calculation and precision evaluation on observed data after recording the site observed data, then wire control network results are obtained, and then the wire network is laid on the return site. The process can not find the error of the observed data in time, and reworking is easy to cause, so that the working efficiency is low, the labor intensity is high, the wire laying speed is low, and the wire laying speed can not be matched with the efficiency of the unmanned aerial vehicle aerial photography technology.
Disclosure of Invention
The invention aims to provide a wire data processing system in hydraulic engineering measurement, and another aim of the invention is to provide a wire data processing method.
In order to achieve the above purpose, the present invention may adopt the following technical scheme:
the invention relates to a wire data processing system in hydraulic engineering measurement, which consists of a server side and a plurality of user sides; the server side is used for carrying out data processing on the measurement data uploaded by each user side and feeding back the measurement data processing results to the corresponding user side in real time; the field measurement personnel guide the measurement data into a remote server through a user terminal, and the measurement data is processed, stored and fed back in real time through the server, so that the field personnel can obtain adjustment results, precision assessment and lead results in time, the problem of finding the measurement data in time is solved, the problem is solved in real time, reworking is avoided, the labor intensity and the cost are further reduced, and the lead layout speed and the working efficiency are improved; and the measurement data are prevented from being tidied again and the measurement result is prevented from being saved in the room.
The user terminal is terminal equipment which has the capability of accessing the Internet and data input, and comprises a mobile phone, a mobile computer and a fixed computer.
The invention relates to a wire data processing method in hydraulic engineering measurement, which comprises data input, adjustment calculation and report output;
the data entry comprises the following steps:
s1.1, setting and determining observation parameters by the field measurement personnel through the user side, wherein the observation parameters comprise instrument types, observation grades, observation points, observation return numbers and weather conditions;
s1.2, carrying out wire observation of a 1 st station by field measurement personnel, and recording measurement data of left reading of a 1 st station horizontal dial, right reading of the horizontal dial, left reading of a vertical dial, right reading of the vertical dial, inclined distance, left length of the dial, right length of the dial, station measurement height and viewpoint height through a user terminal;
s1.3, the user side calculates a 2C value, an index difference and a side length difference of the 1 st station; the 2C value = horizontal dial left reading + 180-horizontal dial right reading; the index difference = (vertical scale left reading + vertical scale right degree-360)/2; the side length is poor = disc left side length-disc right side length;
s1.4, the user side judges whether the 2C value, the index difference and the side length difference of the 1 st station exceed preset limit values; if one of the 2C value, the index difference and the side length difference of the 1 st station exceeds a preset limit value, the 1 st station measurement data is invalidated, and the logging is observed again; if the preset limit value is not exceeded, continuing the next step;
s1.5, the user side calculates the horizontal angle, the vertical angle and the side length of the 1 st station; the horizontal angle= (horizontal dial left reading + horizontal dial right reading)/2; the vertical angle = (vertical scale left reading + vertical scale right reading)/2; the side length= (side length of disc + side length of disc)/2;
s1.6, the user side judges whether the horizontal angle, the vertical angle and the side length data of the 1 st station exceed preset limit values; if one of the horizontal angle, the vertical angle and the side length of the 1 st station exceeds a preset limit value, the 1 st station measurement data is invalidated, and the logging is observed again; if the 1 st station does not exceed the preset limit value, continuing the next step;
s1.7, after the measurement data of the 1 st measuring station is input, the user side firstly performs the warehouse entry inspection of the measurement data, timely feeds the inspection result back to field measuring staff, uploads the measurement data of the 1 st measuring station to the server side through a network after the inspection is correct, and inputs the measurement data into a database of the server side; the server side processes and stores the measurement data uploaded by the user side;
s1.8, a field measurer continuously observes the wire of the station 2, and inputs the measuring data of the left reading of the horizontal dial, the right reading of the horizontal dial, the left reading of the vertical dial, the right reading of the vertical dial, the inclined distance, the left length of the dial, the right length of the dial, the station height and the viewpoint height of the station 2 through a user terminal; then repeating the steps S1.3-S1.7, and sequentially completing the measurement data processing and storage of other measuring stations;
through the steps, the user side feeds back the measurement data condition of each measuring station to the field measuring personnel in real time, whether the measurement data of each measuring station has problems or not is found in time, if the problems are solved in real time on site, the field personnel are prevented from reworking, the labor intensity and the cost are further reduced, and the working efficiency is improved.
The adjustment calculation comprises the following steps:
s2.1, after the measurement data of all the measuring stations are put in storage, the service end calculates the starting azimuth angle and the final azimuth angle of the wire;
the azimuth angle of the starting edge is the azimuth angle of the known points A and B of the starting edge, which is called a in the following description AB The azimuth angle of the final edge is the azimuth angle of the known points C and D of the final edge, and is called a in the following description CD The method comprises the steps of carrying out a first treatment on the surface of the The known point is the abscissa and the ordinate of the point in the measurement coordinate system; the calculation is carried out according to a quadrant formula, and the calculation process is as follows:
，
，
，
，
，
，
the measuring coordinate system takes an x axis as an ordinate and a y axis as an abscissa, and is divided into four quadrants; wherein Deltay AB =Y B -Y A ，△x AB =X B -X A ；X A Represents the ordinate, Y, of the known point A A Represents the abscissa of the known point a; x is X B Represents the ordinate of the known point B, Y B Represents the abscissa of the known point B; deltax AB Representing the difference between the vertical coordinates of the known point B and the known point A; Δy AB Representing the difference between the horizontal coordinates of the known point B and the known point A;
s2.2, according to the advancing direction of the wire, sequentially calculating azimuth angles of every two adjacent points of the wire side;
the advancing direction of the wire is the direction indicated from the starting point A of the wire to the ending point D of the wire; the method comprises the steps that between the point A and the point D, besides the known point B and the known point C, n unknown points are encryption points arranged between the known points B and C according to design requirements, coordinate values of the encryption points are unknown, and the n unknown points correspond to n station measurement data;
wherein in step S2.1, the azimuth angle a between the known point A and the known point B of the wire edge is known AB And azimuth angle a between known point C and known point D on the wire side CD Marking any one of n unknown points as M i Wherein i is more than or equal to 0 and n is more than or equal to a natural number; point M according to the advancing direction of the wire i+1 At point M i Front of (C), point M i-1 At point M i Behind (a), calculate every two adjacent points M of the wire edge i And M is as follows i+1 Azimuth of (2); the calculation process is as follows:
f β = a AB +∑β left side -n*180- a CD ；a MiMi+1 =a Mi-1M +β Left side -f β /n-180；
Wherein: a, a MiMi+1 Representing adjacent points M i And M is as follows i+1 Azimuth angle between a Mi-1Mi Representing adjacent points M i-1 And M is as follows i Azimuth angles therebetween;
when i=0, M i-1 Is M 0-1 I.e. the known point a; m is M i Is M 0 I.e. the known point B;
when i=n, M i+1 Is M n+1 I.e. the known point C; thus, when i=0, the known point B and the unknown point M are calculated 1 The azimuth angle between the two is calculated as follows:
s2.2.1, calculate a CD’ =a AB Left-n 180 + Σβ;
s2.2.2 calculating f β = a CD’ - a CD ；
S2.2.3, calculate a BM1 =a AB +beta left-f β /n-180；
When i=1, it means calculating the known point M 1 With unknown point M 2 The azimuth angle between the two is calculated as follows:
calculation of a CD’ =a AB Left-n 180 + Σβ;
calculating f β = a CD’ - a CD ；
Calculation of a M1M2 =a BM1 +beta left-f β /n-180；
And by analogy, calculating azimuth angles of every two adjacent points of the wire side;
wherein:
∑β left side A sum of horizontal angles representing n unknown points; beta Left side Represents the sum of the horizontal angles of the unknown point Mi+1 and the points Mi and Mi-1 behind it;
s2.3, calculating the incremental coordinates DeltaX of any unknown point Mi on the wire Mi 、△Y Mi The method comprises the steps of carrying out a first treatment on the surface of the The calculation formula is as follows:
△X Mi =d*cosa MiMi+1 ，△Y Mi =d*sina MiMi+1 the method comprises the steps of carrying out a first treatment on the surface of the Wherein d is the side length between the adjacent point Mi and the point Mi+1, and the side length is obtained by measuring by field personnel;
s2.4, calculating an unknown point coordinate increment closing difference; the calculation formula is as follows:
Fx=∑△X –(x C -x B )，Fy=∑△Y –(y C -y B ) Wherein ΣΔx is Δx calculated in step S2.3 Mi Sigma DeltaY is DeltaY calculated in step S2.3 Mi Is the sum of (3); x is X B Representative takeThe starting edge is known as the ordinate of point B, Y B An abscissa representing the known point B of the starting edge; x is X C Representing the ordinate of the known point C of the starting edge, Y C An abscissa representing the known point C of the starting edge;
s2.5, calculating the relative closing difference F of the total length of the wire; the calculation formula is as follows:
F=（Fx 2 + Fy 2 ) k=1/(Σd/F); wherein Σd is the sum of all unknown point side lengths;
s2.6, evaluating whether the relative closing difference of the total length of the wire exceeds a preset limit value, and if so, searching the reason by the server side and returning to the user side;
s2.7, if the current value is not exceeded, the server calculates coordinates of all unknown points in the lead; the calculation formula is as follows:
X Mi+1 =X Mi +△X Mi - Fx*d Mi /∑d，Y Mi+1 =Y Mi +△Y Mi - Fy*d Mi a/Σd; wherein X is Mi+1 Is the ordinate of the point Mi+1, Y Mi+1 The abscissa of the point Mi+1; x is X Mi Is the ordinate, X, of the adjacent point Mi behind the point Mi+1 Mi Is the abscissa of the adjacent point Mi behind the point Mi+1; deltaX Mi And DeltaY Mi Incremental coordinates of the adjacent point Mi behind Mi+1; d, d Mi The side length of the adjacent point Mi behind the point Mi+1;
s2.8, after the calculation is completed, the server stores the calculated data into a database of the server;
the report output comprises the following steps:
s3.1, the user side requests report output to the server side;
s3.2, the server automatically draws a control measurement route map of the calculation result, and simultaneously generates a control point result table and returns the control point result table to the user side; through automatic calculation of the server and drawing of the control measurement route map, field workers can quickly and conveniently obtain measurement results, labor intensity of the field workers is reduced, returning of indoor repeated arrangement of measurement data and preservation of the measurement results are avoided, and working efficiency is improved.
The invention has the advantages that field workers can guide the measured data into a remote server through the user terminal, and the field workers can timely obtain adjustment results, precision evaluation and lead results through automatic processing, storage and real-time feedback of the measured data by the server, so that the problem of timely finding the measured data is realized, the problem is solved in real time, reworking is avoided, the labor intensity and the cost are further reduced, and the lead layout speed and the working efficiency are improved; and the measurement data are prevented from being tidied again and the measurement result is prevented from being saved in the room.
Drawings
Fig. 1 is a diagram of a wire data processing system in hydraulic engineering measurement according to the present invention.
Fig. 2 is a flow chart of the method of the present invention.
FIG. 3 is a flow chart of data entry according to the present invention.
Fig. 4 is a flowchart of the adjustment calculation according to the present invention.
Fig. 5 is a four-quadrant schematic of the present invention.
FIG. 6 is a flowchart of report output according to the present invention.
FIG. 7 is a table of control point achievements generated by the method of the present invention.
Fig. 8 is a control measurement roadmap generated by the method of the invention.
Detailed Description
The following description of the technical solutions in the embodiments of the present invention will be clear and complete, and it is obvious that the described embodiments are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
As shown in fig. 1, the wire data processing system in hydraulic engineering measurement of the present invention comprises a server side and a plurality of user sides; the server side is used for carrying out data processing on the measurement data uploaded by each user side and feeding back the data processing result to the corresponding user side in real time; the user terminal can be any device which can be accessed to the Internet and has data input capability, such as a mobile phone, a mobile computer, a fixed computer and the like; the field staff can guide the measured data into a remote server through the user side, and can timely obtain adjustment results, precision assessment and lead results through automatic processing, storage and real-time feedback of the measured data by the server, so that the problem of timely finding the measured data is solved, the problem is timely solved, reworking is avoided, the labor intensity and cost are further reduced, and the lead layout speed and the working efficiency are improved; and the measurement data are prevented from being tidied again and the measurement result is prevented from being saved in the room.
As shown in FIG. 2, the wire data processing method in hydraulic engineering measurement comprises data input, adjustment calculation and report output;
as shown in fig. 3, the data entry includes the steps of:
s1.1, field measurement personnel set and determine observation parameters through a user side, wherein the observation parameters comprise instrument types, observation grades, observation points, observation measurement back numbers and weather conditions;
s1.2, carrying out wire observation of a 1 st station by field measurement personnel, and recording measurement data such as left reading of a 1 st station horizontal dial, right reading of the horizontal dial, left reading of a vertical dial, right reading of the vertical dial, inclined distance, left length of the dial, right length of the dial, station measurement height, viewpoint height and the like through a user terminal; the user terminal can be a mobile phone, a mobile computer and a fixed computer;
s1.3, the user side calculates a 2C value, an index difference and a side length difference of the 1 st station; the 2C value = horizontal dial left reading + 180-horizontal dial right reading; the index difference = (vertical scale left reading + vertical scale right degree-360)/2; the side length is poor = disc left side length-disc right side length;
s1.4, the user side judges whether the 2C value, the index difference and the side length difference of the 1 st station exceed preset limit values; if one of the 2C value, the index difference and the side length difference of the 1 st station measurement data exceeds a preset limit value, the 1 st station measurement data is invalidated, and the logging is required to be observed again; if the preset limit value is not exceeded, continuing the next step;
s1.5, the user side calculates the horizontal angle, the vertical angle and the side length of the 1 st station; the horizontal angle= (horizontal dial left reading + horizontal dial right reading)/2; the vertical angle = (vertical scale left reading + vertical scale right reading)/2; the side length= (side length of disc + side length of disc)/2;
s1.6, the user side judges whether the horizontal angle, the vertical angle and the side length data of the 1 st station exceed preset limit values; if one of the horizontal angle, the vertical angle and the side length of the 1 st station measurement data exceeds a preset limit value, the 1 st station measurement data is invalidated, and the logging is required to be observed again; if the 1 st station does not exceed the preset limit value, continuing the next step;
s1.7, after the measurement data of the 1 st measuring station is input, the user side firstly performs the warehouse entry inspection of the measurement data, timely feeds the inspection result back to field measurement personnel, uploads the measurement data of the 1 st measuring station to the server side through a network after the inspection is correct, and inputs the measurement data into a database of the server side; the server side processes and stores the measurement data uploaded by the user side;
s1.8, the field survey staff continues to observe the wire of the station 2, and inputs the measurement data such as the left reading of the horizontal dial of the station 2, the right reading of the horizontal dial, the left reading of the vertical dial, the right reading of the vertical dial, the inclined distance, the left length of the dial, the right length of the dial, the station height, the viewpoint height and the like through a user terminal; repeating the steps S1.3-S1.7, and sequentially completing the measurement data processing and storage of other measuring stations;
through the steps, the user side feeds back the measurement data condition of each measuring station to the field measuring personnel in real time, whether the measurement data of each measuring station has problems or not is found in time, if the problems are solved in real time on site, the field personnel are prevented from reworking, the labor intensity and the cost are further reduced, and the working efficiency is improved.
As shown in fig. 4, the adjustment calculation includes the following steps:
s2.1, after the measurement data of all the measuring stations are put in storage, the service end calculates the starting azimuth angle and the final azimuth angle of the wire;
the azimuth angle of the starting edge is the azimuth angle of the known points A and B of the starting edge, and is as followsIn the description simply called a AB The azimuth angle of the final edge is the azimuth angle of the known points C and D of the final edge, and is called a in the following description CD The method comprises the steps of carrying out a first treatment on the surface of the The known point is the abscissa and the ordinate of the known point; the calculation is calculated according to a quadrant formula, and a is calculated according to a quadrant formula AB For example, the calculation process is as follows:
，
，
，
，
，
，
the measurement coordinate system is shown in fig. 5, and is divided into four quadrants by taking the x axis as an ordinate and the y axis as an abscissa; wherein Deltay AB =Y B -Y A ，△x AB =X B -X A ，X A Represents the ordinate, Y, of the known point A A Represents the abscissa of the known point a; x is X B Represents the ordinate of the known point B, Y B Represents the abscissa of the known point B; deltax AB Representing the difference between the vertical coordinates of the known point B and the known point A; Δy AB Representing the difference between the horizontal coordinates of the known point B and the known point A;
s2.2, according to the advancing direction of the wire, sequentially calculating azimuth angles of every two adjacent points of the wire side;
the advancing direction of the wire is the direction indicated from the starting point A of the wire to the ending point D of the wire; the method comprises the steps that between the point A and the point D, besides the known point B and the known point C, n unknown points are encryption points arranged between the known points B and C according to design requirements, coordinate values of the encryption points are unknown, and the n unknown points correspond to n station measurement data;
wherein in step S2.1, the azimuth angle a between the known point A and the known point B of the wire edge is known AB And azimuth angle a between known point C and known point D on the wire side CD Marking any one of n unknown points as M i Wherein i is more than or equal to 0 and n is more than or equal to a natural number; point M according to the advancing direction of the wire i+1 At point M i Front of (C), point M i-1 At point M i Behind (a), calculate every two adjacent points M of the wire edge i And M is as follows i+1 Azimuth of (2); the calculation process is as follows:
f β = a AB +∑β left side -n*180- a CD ；a MiMi+1 =a Mi-1M +β Left side -f β /n-180；
Wherein: a, a MiMi+1 Representing adjacent points M i And M is as follows i+1 Azimuth angle between a Mi-1Mi Representing adjacent points M i-1 And M is as follows i Azimuth angles therebetween;
when i=0, M i-1 Is M 0-1 I.e. the known point a; m is M i Is M 0 I.e. the known point B;
when i=n, M i+1 Is M n+1 I.e. the known point C; thus, when i=0, the known point B and the unknown point M are calculated 1 The azimuth angle between the two is calculated as follows:
s2.2.1, calculate a CD’ =a AB Left-n 180 + Σβ;
s2.2.2 calculating f β = a CD’ - a CD ；
S2.2.3, calculate a BM1 =a AB +beta left-f β /n-180；
When i=1, it means calculating the known point M 1 With unknown point M 2 The azimuth angle between the two is calculated as follows:
calculation of a CD’ =a AB Left-n 180 + Σβ;
calculating f β = a CD’ - a CD ；
Calculation of a M1M2 =a BM1 +beta left-f β /n-180；
And by analogy, calculating azimuth angles of every two adjacent points of the wire side;
wherein:
∑β left side Representing the sum of n unknown point horizontal angles; beta Left side Represents the sum of the horizontal angles of the unknown point Mi+1 and the points Mi and Mi-1 behind it; a, a MiMi+1 Represents the azimuth angle between the adjacent point Mi to be calculated and the point Mi+1, a Mi-1Mi Representing the azimuth angle between the adjacent point Mi-1 and the point Mi;
s2.3, calculating the increment coordinates of any unknown point Mi on the lead; the increment coordinates are abbreviated as DeltaX Mi 、△Y Mi The method comprises the steps of carrying out a first treatment on the surface of the The calculation process is as follows:
△X Mi =d*cosa MiMi+1 ，△Y Mi =d*sina MiMi+1 the method comprises the steps of carrying out a first treatment on the surface of the Wherein d is the side length between the adjacent point Mi and the point Mi+1, and is obtained by measuring by field personnel;
s2.4, calculating an unknown point coordinate increment closing difference; the calculation process is as follows:
Fx=∑△X –(x C -x B )，Fy=∑△Y –(y C -y B ) Wherein ΣΔx is Δx calculated in step S2.3 Mi Is the sum of (3); wherein ΣΔy is Δy calculated in step S2.3 Mi Is the sum of (3); x is X B Represents the ordinate of the known point B of the starting edge, Y B An abscissa representing the known point B of the starting edge; x is X C Representing the ordinate of the known point C of the starting edge, Y C An abscissa representing the known point C of the starting edge;
s2.5, calculating the relative closing difference F of the total length of the wire; f= (Fx 2 + Fy 2 ) k=1/(Σd/F); wherein Σd is the sum of all unknown point side lengths;
s2.6, evaluating whether the total length relative closing difference of the lead exceeds a preset limit value or not, and judging the precision grade; if the threshold is exceeded, the server side searches the reason and returns to the user side;
s2.7, if the current value is not exceeded, the server calculates coordinates of all unknown points in the lead; the calculation process is illustrated below by taking the unknown point Mi+1 as an example:
X Mi+1 =X Mi +△X Mi - Fx*d Mi /∑d，Y Mi+1 =Y Mi +△Y Mi - Fy*d Mi a/Σd; wherein X is Mi+1 Is the ordinate of the point Mi+1, Y Mi+1 The abscissa of the point Mi+1; x is X Mi Is the ordinate, X, of the adjacent point Mi behind the point Mi+1 Mi Is the abscissa of the adjacent point Mi behind the point Mi+1; deltaX Mi And DeltaY Mi Incremental coordinates of the adjacent point Mi behind Mi+1; d, d Mi The side length of the adjacent point Mi behind the point Mi+1;
s2.8, after the calculation is completed, the server stores the calculated data into a database of the server;
as shown in fig. 6, the report output includes the following steps:
s3.1, the user side requests report output to the server side;
s3.2, the server automatically draws a control measurement route diagram of the calculation result, as shown in fig. 7, and simultaneously generates a control point result table, as shown in fig. 8; and returns to the user side; the measuring route map is automatically calculated and drawn by the server, so that field workers can quickly and conveniently obtain the measuring result, the labor intensity of the field workers is reduced, the situation that the measuring data are returned to the room for finishing again and the measuring result is saved is avoided, and the working efficiency is improved;
s3.3, the user side can share the calculation result to related staff; and measurement data and achievements in the project group are shared, so that the working efficiency is improved.
Claims (3)
1. A wire data processing method in hydraulic engineering measurement is characterized in that: the method comprises data input, adjustment calculation and report output;
the data entry comprises the following steps:
s1.1, field measurement personnel set and determine observation parameters through a user side, wherein the observation parameters comprise instrument types, observation grades, observation points, observation measurement back numbers and weather conditions;
s1.2, carrying out wire observation of a 1 st station by field measurement personnel, and recording measurement data of left reading of a horizontal dial, right reading of the horizontal dial, left reading of a vertical dial, right reading of the vertical dial, inclined distance, left length of the dial, right length of the dial, station measurement height and viewpoint height of the 1 st station through a user terminal;
s1.3, the user side calculates a 2C value, an index difference and a side length difference of the 1 st station; the 2C value = horizontal dial left reading + 180-horizontal dial right reading; the index difference = (vertical scale left reading + vertical scale right degree-360)/2; the side length is poor = disc left side length-disc right side length;
s1.4, the user side judges whether the 2C value, the index difference and the side length difference of the 1 st station exceed preset limit values; if one of the 2C value, the index difference and the side length difference of the 1 st station exceeds a preset limit value, the 1 st station measurement data is invalidated, and the logging is observed again; if the preset limit value is not exceeded, continuing the next step;
s1.5, the user side calculates the horizontal angle, the vertical angle and the side length of the 1 st station; the horizontal angle= (horizontal dial left reading + horizontal dial right reading)/2; the vertical angle = (vertical scale left reading + vertical scale right reading)/2; the side length= (side length of disc + side length of disc)/2;
s1.6, the user side judges whether the horizontal angle, the vertical angle and the side length data of the 1 st station exceed preset limit values; if one of the horizontal angle, the vertical angle and the side length of the 1 st station exceeds a preset limit value, the 1 st station measurement data is invalidated, and the logging is observed again; if the preset limit value is not exceeded, continuing the next step;
s1.7, after the measurement data of the 1 st measuring station is input, the user side firstly performs the warehouse entry inspection of the measurement data, timely feeds the inspection result back to field measurement personnel, uploads the measurement data of the 1 st measuring station to the server side through a network after the inspection is correct, and inputs the measurement data into a database of the server side; the server side processes and stores the measurement data uploaded by the user side;
s1.8, a field measurer continuously observes the wire of the station 2, and inputs the measuring data of the left reading of the horizontal dial, the right reading of the horizontal dial, the left reading of the vertical dial, the right reading of the vertical dial, the inclined distance, the left length of the dial, the right length of the dial, the station height and the viewpoint height of the station 2 through a user terminal; then repeating the steps S1.3-S1.7, and sequentially completing the measurement data processing and storage of other measuring stations;
the adjustment calculation comprises the following steps:
s2.1, after the measurement data of all the measuring stations are put in storage, the service end calculates the starting azimuth angle and the final azimuth angle of the wire;
the azimuth angle of the starting edge is the azimuth angle of the known points A and B of the starting edge, which is called a in the following description AB The azimuth angle of the final edge is the azimuth angle of the known points C and D of the final edge, and is called a in the following description CD The method comprises the steps of carrying out a first treatment on the surface of the The known point is the abscissa and the ordinate of the point in the measurement coordinate system; the calculation is carried out according to a quadrant formula, and the calculation process is as follows:
，
，
；
，
，
；
the measuring coordinate system takes an x axis as an ordinate and a y axis as an abscissa, and is divided into four quadrants; wherein Deltay AB =Y B -Y A ，△x AB =X B -X A ；X A Represents the ordinate, Y, of the known point A A Represents the abscissa of the known point a; x is X B Represents the ordinate of the known point B, Y B Represents the abscissa of the known point B; deltax AB Representing the difference between the vertical coordinates of the known point B and the known point A; Δy AB Representing the difference between the horizontal coordinates of the known point B and the known point A;
s2.2, according to the advancing direction of the wire, sequentially calculating azimuth angles of every two adjacent points of the wire side;
the advancing direction of the wire is the direction indicated from the starting point A of the wire to the ending point D of the wire; the method comprises the steps that between the point A and the point D, besides the known point B and the known point C, n unknown points are encryption points arranged between the known points B and C according to design requirements, coordinate values of the encryption points are unknown, and the n unknown points correspond to n station measurement data; wherein in step S2.1, the azimuth angle a between the known point A and the known point B of the wire edge is known AB And azimuth angle a between known point C and known point D on the wire side CD Marking any one of n unknown points as M i Wherein i is more than or equal to 0 and n is more than or equal to a natural number; point M according to the advancing direction of the wire i+1 At point M i Front of (C), point M i-1 At point M i Behind (a), calculate every two adjacent points M of the wire edge i And M is as follows i+1 Azimuth of (2); the calculation process is as follows:
f β = a AB +∑β left side -n*180- a CD ；a MiMi+1 =a Mi-1M +β Left side -f β /n-180；
Wherein: a, a MiMi+1 Representing adjacent points M i And M is as follows i+1 Azimuth angle between a Mi-1Mi Representing adjacent points M i-1 And M is as follows i Azimuth angles therebetween;
when i=0, M i-1 Is M 0-1 I.e. the known point a; m is M i Is M 0 I.e. the known point B;
when i=n, M i+1 Is M n+1 I.e. the known point C; thus, when i=0, the known point B and the unknown point M are calculated 1 The azimuth angle between the two is calculated as follows:
s2.2.1, calculate a CD’ =a AB +∑β Left side -n*180；
S2.2.2 calculating f β = a CD’ -a CD ；
S2.2.3, calculate a BM1 =a AB +β Left side -f β /n-180；
When i=1, it means calculating the known point M 1 With unknown point M 2 The azimuth angle between the two is calculated as follows:
calculation of a CD’ =a AB +∑β Left side -n*180；
Calculating f β =a CD’ -a CD ；
Calculation of a M1M2 =a BM1 +β Left side -f β /n-180；
And by analogy, calculating azimuth angles of every two adjacent points of the wire side;
wherein:
∑β left side A sum of horizontal angles representing n unknown points; beta Left side Representing an unknown point M i+1 And its rear point M i Sum point M i-1 Is the sum of the horizontal angles of (2);
s2.3, calculating any unknown point M on the wire i Is a function of the incremental coordinates of (a); the increment coordinates are abbreviated as DeltaX Mi 、△Y Mi The method comprises the steps of carrying out a first treatment on the surface of the The calculation formula is as follows:
△X Mi =d*cosa MiMi+1 ，△Y Mi =d*sina MiMi+1 the method comprises the steps of carrying out a first treatment on the surface of the Wherein d is the adjacent point M i And point M i+1 The side length of the two parts is measured and obtained by field personnel;
s2.4, calculating an unknown point coordinate increment closing difference; the calculation formula is as follows:
Fx=∑△X –(X C -X B )，Fy=∑△Y –(Y C -Y B ) Wherein ΣΔx is Δx calculated in step S2.3 Mi Sigma DeltaY is DeltaY calculated in step S2.3 Mi Is the sum of (3); x is X B Represents the ordinate of the known point B of the starting edge, Y B An abscissa representing the known point B of the starting edge; x is X C Representing the ordinate of the known point C of the starting edge, Y C An abscissa representing the known point C of the starting edge;
s2.5, calculating the relative closing difference K of the total length of the wire; the calculation formula is as follows:
k=1/(Σd/F); wherein Σd is the sum of all unknown point side lengths;
s2.6, evaluating whether the relative closing difference of the total length of the wire exceeds a preset limit value, and if so, searching the reason by the server side and returning to the user side;
s2.7, if the current value is not exceeded, the server calculates coordinates of all unknown points in the lead; the calculation formula is as follows:
X Mi+1 =X Mi +△X Mi - Fx*d Mi /∑d，Y Mi+1 =Y Mi +△Y Mi - Fy*d Mi a/Σd; wherein X is Mi+1 For point M i+1 Y is the ordinate of (2) Mi+1 For point M i+1 Is the abscissa of (2); x is X Mi For point M i+ 1 rear adjacent point M i X, X Mi For point M i+1 Rear adjacent point M i Is the abscissa of (2); deltaX Mi And DeltaY Mi Is M i+1 Rear adjacent point M i Is a function of the incremental coordinates of (a); d, d Mi For point M i+1 Rear adjacent point M i Is a side length of (2);
s2.8, after the calculation is completed, the server stores the calculated data into a database of the server;
the report output comprises the following steps:
s3.1, the user side requests report output to the server side;
and S3.2, the server automatically draws the calculation result into a control measurement route map, simultaneously generates a control point result table and returns the control point result table to the user side.
2. A wire data processing system in hydraulic engineering measurement is used in the method of claim 1, and consists of a service end and a plurality of user ends; the method is characterized in that: the server side is used for carrying out data processing on the measurement data uploaded by each user side and feeding back the measurement data processing result to the corresponding user side in real time; the field measurement personnel can import the measurement data into a remote server through the user terminal, and the measurement data is processed, stored and fed back in real time through the server, so that the field measurement personnel can obtain adjustment results, precision assessment and lead results in time.
3. The system for processing wire data in hydraulic engineering measurement according to claim 2, wherein: the user terminal is terminal equipment which has the capability of accessing the Internet and data input, and comprises a mobile phone, a mobile computer and a fixed computer.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110762742.8A CN113536559B (en) | 2021-07-06 | 2021-07-06 | Wire data processing system and method in hydraulic engineering measurement |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110762742.8A CN113536559B (en) | 2021-07-06 | 2021-07-06 | Wire data processing system and method in hydraulic engineering measurement |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113536559A CN113536559A (en) | 2021-10-22 |
| CN113536559B true CN113536559B (en) | 2023-08-01 |
Family
ID=78126880
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110762742.8A Active CN113536559B (en) | 2021-07-06 | 2021-07-06 | Wire data processing system and method in hydraulic engineering measurement |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113536559B (en) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111866253A (en) * | 2020-07-24 | 2020-10-30 | 耿天真 | Mobile terminal software testing system based on big data sharing |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2014090310A1 (en) * | 2012-12-13 | 2014-06-19 | Abb Technology Ag | System and method for monitoring and/or diagnosing operation of a production line of an industrial plant |
| CN103926913A (en) * | 2014-05-08 | 2014-07-16 | 中国水利水电第七工程局有限公司 | Total station tunnel deformation remote monitoring system and establishing method thereof |
| CN106197335A (en) * | 2016-06-14 | 2016-12-07 | 武汉瑞得信息工程有限责任公司 | A kind of control measuring method and system |
| CN110069584B (en) * | 2017-11-29 | 2023-04-14 | 星际空间（天津）科技发展有限公司 | Information acquisition and feature matching method based on mobile terminal |
| CN108534766A (en) * | 2018-04-09 | 2018-09-14 | 辽宁林业职业技术学院 | A kind of improvement computational methods of back traverse |
| CN109947024A (en) * | 2019-03-29 | 2019-06-28 | 中铁隧道局集团有限公司 | A kind of rail traffic data monitoring method based on mobile terminal |
-
2021
- 2021-07-06 CN CN202110762742.8A patent/CN113536559B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111866253A (en) * | 2020-07-24 | 2020-10-30 | 耿天真 | Mobile terminal software testing system based on big data sharing |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113536559A (en) | 2021-10-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104703143A (en) | Indoor positioning method based on WIFI signal strength | |
| CN112767464A (en) | Ground laser scanning three-dimensional point cloud data registration method | |
| CN109826248A (en) | Method for applying 3D Laser Scanner to foundation pit deformation monitoring | |
| CN107860374A (en) | A kind of tunnel tunnel face excavation line intelligent positioning measuring method | |
| CN112857315B (en) | Steel pipe column verticality measuring method based on three-dimensional laser scanning | |
| CN105651311A (en) | Method for measuring satellite navigation autopilot accuracy of agricultural machinery operation | |
| CN113536559B (en) | Wire data processing system and method in hydraulic engineering measurement | |
| CN112884823B (en) | Guiding rule simulation point cloud selection and calculation method for actual measurement actual quantity application | |
| CN108415105B (en) | Method for inspecting relative humidity value observed by ground meteorological observation station | |
| CN113450076B (en) | Hydraulic engineering measures integrated intelligent processing service platform | |
| CN111550604A (en) | Pipeline detection method of directional director | |
| CN110648280A (en) | Data processing method for splicing large-scale karst cave mass point cloud data | |
| CN113873429B (en) | Underground tunnel-oriented ridge regression two-dimensional positioning method and system | |
| CN111787608B (en) | Anchor node laying method and device | |
| CN113536194B (en) | Leveling data processing system and method in hydraulic engineering measurement | |
| CN106683126A (en) | Method for quantitatively evaluating distribution quality of targets in point cloud registration | |
| CN113567966B (en) | Monte Carlo simulation-based airborne/vehicular laser point cloud precision prediction method | |
| Yeh et al. | Establishment of Taiwan’s standard operating procedure for network-based RTK in cadastral surveying | |
| CN111629321A (en) | Median Kalman filtering positioning method based on outlier stripping | |
| CN112946736A (en) | Three-dimensional observation system reconstruction method and system | |
| CN117935080B (en) | Power transmission and transformation project geological distribution map drawing processing method based on big data | |
| Hermawan et al. | Evaluation of Precision and Accuracy of Mobile Mapping System (MMS) Leica Pegasus Two Ultimate for Road Monitoring | |
| CN110207699A (en) | A kind of localization method and device | |
| CN115033835B (en) | Ballistic parameter calculation method and system based on multiple interaction least square estimation | |
| NL2030845B1 (en) | Positioning and correction method and system of video data using seabed topographic profile |
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 |