CN113375644A - High-precision side slope occurrence measuring method - Google Patents
High-precision side slope occurrence measuring method Download PDFInfo
- Publication number
- CN113375644A CN113375644A CN202110783751.5A CN202110783751A CN113375644A CN 113375644 A CN113375644 A CN 113375644A CN 202110783751 A CN202110783751 A CN 202110783751A CN 113375644 A CN113375644 A CN 113375644A
- Authority
- CN
- China
- Prior art keywords
- slope
- points
- equation
- occurrence
- triangular
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
- G01C11/04—Interpretation of pictures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2200/00—Indexing scheme for image data processing or generation, in general
- G06T2200/08—Indexing scheme for image data processing or generation, in general involving all processing steps from image acquisition to 3D model generation
-
- 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
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/23—Dune restoration or creation; Cliff stabilisation
Abstract
The invention relates to a high-precision measurement method for slope occurrence, which belongs to the technical field of measurement and comprises the following steps: acquiring geological attribute information by using an unmanned aerial vehicle aerial photography technology; converting the acquired geological attribute information into a tilted photography model by using a graphic workstation, and then importing the tilted photography model into mobile equipment; determining a slope range within the oblique photography model; randomly generating n slope point data in the slope range, and then screening and filtering outlier slope points to obtain reserved slope points; fitting a slope surface equation by adopting a least square method, and solving the occurrence of the side slope through the slope surface equation; encrypting by using a Delaunay triangulation network, circularly traversing all the triangular surfaces, and solving the occurrence of all the triangular surfaces; and comparing the side slope attitude with the attitude of all the triangular faces, and determining the triangular face closest to the attitude elements of the side slope.
Description
Technical Field
The disclosure relates to the technical field of measurement, in particular to a high-precision measurement method for slope occurrence.
Background
The side slope occurrence measurement is a basic work in the field geological exploration process, and lays a foundation for the analysis of the stability of the side slope surface. The traditional measuring method is to measure by means of a geological compass, and has the disadvantages of complex process, large error, easy limitation to terrain conditions and poor implementation effect. With the continuous development of science and technology and the continuous enhancement of computer parallel computing capability, the technical field starts to utilize an unmanned aerial vehicle to collect surface data, then a series of datum points are selected to perform trial calculation, and finally a maximum apparent dip angle and a corresponding azimuth angle are selected as the attitude.
Disclosure of Invention
The present disclosure provides a high-precision measurement method for a side slope occurrence, which is used for solving the problems of calculation result deviation and incapability of in-depth analysis in the background art.
In view of the above problems, the present disclosure provides a method for measuring a slope occurrence with high accuracy, including:
carrying out aerial photography on the area to be measured by using an unmanned aerial vehicle aerial photography technology, and acquiring geological attribute information;
converting the acquired geological attribute information into a tilted photography model by using a graphic workstation, and then importing the tilted photography model into mobile equipment;
determining a natural side slope to be measured in the oblique photography model, randomly selecting at least three reference points on the natural side slope, and determining a side slope range;
randomly generating n slope point data within the slope range, and then screening and filtering the outlier slope points to obtain reserved slope points, wherein n is an integer greater than or equal to 1;
fitting a slope equation by the datum points by a least square method, and solving the occurrence of the side slope through the slope equation;
encrypting the reserved slope points by using a Delou internal triangular network, circularly traversing all the triangular surfaces, and solving the occurrence of all the triangular surfaces;
and comparing the side slope occurrence with the occurrences of all the triangular surfaces, and determining the triangular surface closest to the occurrence elements of the side slope.
According to the embodiment of the disclosure, the oblique photography model comprises a three-dimensional real scene model which is formed by superposing images acquired by aerial photography of the unmanned aerial vehicle and elevation data and can represent a field terrain environment.
According to an embodiment of the present disclosure, the method for determining a slope range includes: and clicking the screen through the mobile equipment to select the reference points, clicking three or more representative reference points, and determining the approximate range of the slope.
According to the embodiment of the present disclosure, the reference points are uniformly distributed.
According to the embodiment of the disclosure, the method for screening and filtering the outlier slope points comprises the following steps:
determining a range of a horizontal coordinate X, Y according to the selected reference points, randomly combining n groups of (X, Y) data point pairs in the range, and extracting elevation information on the oblique photography model according to the oblique photography model and the n groups of (X, Y) data point pairs to obtain a first set of three-dimensional coordinates of the n groups of (X, Y) data point pairs;
and then automatically screening the outlier slope points by using an isolated forest algorithm, and filtering to obtain the retained second set of three-dimensional coordinates of the slope points.
According to the embodiment of the disclosure, the method for solving the slope occurrence comprises the following steps:
selecting three reference points from the point selection, and preliminarily drawing up a slope surface equation S0：
A0X+B0Y+C0Z+1＝0 (1)
Wherein: a. the0、B0、C0Is the coefficient of the plane equation;
x, Y, Z is a spatial rectangular three-dimensional coordinate;
assuming that the final slope surface equation fitted by all the selected datum points is S:
AX+BY+CZ+1＝0 (2)
wherein: A. b, C are coefficients of a plane equation;
x, Y, Z is a spatial rectangular three-dimensional coordinate;
then S0And S satisfies the following relationship:
A＝A0+A*,B＝B0+B*,C＝C0+C* (3)
wherein: a. the*、B*、C*Respectively representing A, B, C corresponding error coefficients;
the reference point (X) is obtained by bringing the formula (3) into the formula (2)i,Yi,Zi) The fitted slope surface equation S is:
S＝A*Xi+B*Yi+C*Zi+A0Xi+B0Yi+C0Zi+1＝A*Xi+B*Yi+C*Zi+Ni (4)
wherein: n is a radical ofiThe error caused by the equation introduced by the ith slope point is represented as a constant term;
reference point is not less than 3, using least square method, using point (X)i,Yi,Zi) I-0-n-1, the slope equation should be fitted, so that equation (5) is minimized,
to make it
The system of equations is obtained as:
solving and calculating a linear equation set (6), solving to obtain an error coefficient, and finally obtaining a slope equation S:
AfX+BfY+CfZ+1＝0 (7)
wherein: a. thef、Bf、CfFinal coefficients representing slope equation X, Y, Z;
the trend of the slope is obtained, the intersection line of the slope S and the plane Z is 0, and the calculation formula is as follows:
(1) when B is presentfWhen the value is 0, the trend of the slope is as follows:
θ1＝90 (8)
(2) when B is presentfWhen not equal to 0, the trend calculation formula of the side slope is as follows:
θ1＝arctan(-Af/Bf) (9)
and (3) calculating the inclination angle of the side slope, wherein the calculation formula is as follows:
and (3) solving the slope inclination, considering the quadrant of the projection of the slope equation normal vector on the X-Y plane, and specifically calculating the following steps:
(1) when C is presentf＞0：
(2) When C is presentf＜0：
According to an embodiment of the present disclosure, a method of determining a triangular surface closest to a slope includes:
combining the reserved slope surface points and the reference points, and encrypting by using a Delaunay triangulation network;
solving the occurrence of each triangular surface through cyclic iteration;
based on all the triangular surface attitude and the side slope attitude, any attribute of the trend, the inclination and the dip angle is selected, and all the triangular surfaces and the side slopes are visualized in a thermodynamic diagram mode, so that the triangular surface closest to the side slope attitude elements is determined.
According to an embodiment of the present disclosure, a method of visualizing using thermodynamic diagrams includes: by calculating:
wherein, XfRepresenting any occurrence attribute on any triangle surface;
Xarearepresenting the corresponding occurrence attribute on the final slope surface equation;
According to an embodiment of the present disclosure, the mobile device includes a mobile tablet.
According to the embodiment of the disclosure, the system adopted by the mobile device comprises a Windows operating system.
The unmanned aerial vehicle aerial photography technology is suitable for field geological survey exploration, the three-dimensional occurrence of the side slope can be safely and accurately determined, automatic measurement of the three-dimensional occurrence of the slope is realized by the graphic workstation, the oblique photography model is displayed through the three-dimensional geological real scene platform, a large amount of labor force can be saved, the occurrence measurement can be carried out on regions with severe natural conditions, such as high side slopes and cliff, which are difficult for people to reach, and the field measurement mode is enriched.
The slope range is determined by selecting the datum points, then a plurality of slope point data are randomly generated in the slope range, and then the outlier slope points are screened and filtered, so that the deviation caused by selecting points with abnormal large elevation or points with concave slope is effectively avoided.
The occurrence of the side slope is solved by fitting the slope equation through the datum points, the reserved slope points are encrypted by utilizing a Delou internal triangular net, the occurrence of all triangular surfaces is solved, the occurrence of the side slope is compared with the occurrence of all the triangular surfaces, the triangular surface closest to the occurrence elements of the side slope is determined, and therefore geological prospecting personnel can conveniently conduct follow-up deep analysis, and the accuracy of geological exploration is improved.
Drawings
The foregoing and other objects, features and advantages of the disclosure will be apparent from the following description of embodiments of the disclosure, which proceeds with reference to the accompanying drawings, in which:
fig. 1 schematically shows an architectural diagram of the present disclosure.
Detailed Description
Hereinafter, embodiments of the present disclosure will be described with reference to the accompanying drawings. It should be understood that the description is illustrative only and is not intended to limit the scope of the present disclosure. In the following detailed description, for purposes of explanation, numerous specific details are set forth in order to provide a thorough understanding of the embodiments of the disclosure. It may be evident, however, that one or more embodiments may be practiced without these specific details. Moreover, in the following description, descriptions of well-known structures and techniques are omitted so as to not unnecessarily obscure the concepts of the present disclosure.
The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the disclosure. The terms "comprises," "comprising," and the like, as used herein, specify the presence of stated features, steps, operations, and/or components, but do not preclude the presence or addition of one or more other features, steps, operations, or components.
All terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art unless otherwise defined. It is noted that the terms used herein should be interpreted as having a meaning that is consistent with the context of this specification and should not be interpreted in an idealized or overly formal sense.
Where a convention analogous to "at least one of A, B and C, etc." is used, in general such a construction is intended in the sense one having skill in the art would understand the convention (e.g., "a system having at least one of A, B and C" would include but not be limited to systems that have a alone, B alone, C alone, a and B together, a and C together, B and C together, and/or A, B, C together, etc.).
Fig. 1 schematically shows an architectural diagram of the present disclosure. As shown in fig. 1, according to an embodiment of the present disclosure, a slope occurrence high-precision measurement system includes: unmanned aerial vehicle, figure workstation, mobile device. The unmanned aerial vehicle is used for aerial photography of an area to be measured, collecting geological attribute information, the geological attribute information collected by the unmanned aerial vehicle is converted into an oblique photography model by the graphic workstation, then the oblique photography model is guided into the mobile device, measurement of side slope occurrence is achieved through the mobile device, and finally the measured occurrence information is stored in the database.
As shown in fig. 1, according to an embodiment of the present disclosure, there is provided a high-precision side slope occurrence measurement method including:
carrying out aerial photography on the area to be measured by using an unmanned aerial vehicle aerial photography technology, and acquiring geological attribute information;
converting the acquired geological attribute information into a tilted photography model by using a graphic workstation, and then importing the tilted photography model into mobile equipment;
determining a natural side slope to be measured in an oblique photography model, randomly selecting at least three reference points on the natural side slope, and determining a side slope range;
randomly generating n slope point data within a side slope range, and then screening and filtering outlier slope points to obtain reserved slope points, wherein n is an integer greater than or equal to 1;
fitting a slope equation by the datum point by adopting a least square method, and solving the occurrence of the side slope through the slope equation;
the reserved slope points are encrypted by utilizing a Delaunay triangulation network, all the triangular surfaces are traversed circularly, and the occurrence of all the triangular surfaces is solved;
and comparing the occurrence of the side slope with the occurrence of all the triangular surfaces, and determining the triangular surface closest to the occurrence elements of the side slope.
The unmanned aerial vehicle aerial photography technology is suitable for field geological survey exploration, the three-dimensional occurrence of the side slope can be safely and accurately determined, automatic measurement of the three-dimensional occurrence of the slope is realized by the graphic workstation, the oblique photography model is displayed through the three-dimensional geological real scene platform, a large amount of labor force can be saved, the occurrence measurement can be carried out on regions with severe natural conditions, such as high side slopes and cliff, which are difficult for people to reach, and the field measurement mode is enriched.
The slope range is determined by selecting the datum points, then a plurality of slope point data are randomly generated in the slope range, and then the outlier slope points are screened and filtered, so that the deviation caused by selecting points with abnormal large elevation or points with concave slope is effectively avoided.
The occurrence of the side slope is solved by fitting the slope equation through the datum points, the reserved slope points are encrypted by utilizing a Delou internal triangular net, the occurrence of all triangular surfaces is solved, the occurrence of the side slope is compared with the occurrence of all the triangular surfaces, the triangular surface closest to the occurrence elements of the side slope is determined, and therefore geological prospecting personnel can conveniently conduct follow-up deep analysis, and the accuracy of geological exploration is improved.
As shown in fig. 1, according to an embodiment of the present disclosure, the oblique photography model includes a three-dimensional live-action model formed by superimposing an image acquired by unmanned aerial vehicle aerial photography and elevation data, and capable of representing a field terrain environment.
The high-precision unmanned aerial vehicle is used as a main tool for collecting geological attribute information, earth surface data are collected from a plurality of different angles, then a base map of aerial photography is rendered by utilizing a graphic workstation, elevation data are combined, an oblique photography model is manufactured, the oblique photography model not only contains terrain elevation information, but also contains specific earth surface information, the accuracy of the earth surface data is improved, and the oblique photography model can be applied to the fields of water conservancy, mineral products, resource exploration and the like.
As shown in fig. 1, according to an embodiment of the present disclosure, the mobile device includes a mobile tablet, and the system employed by the mobile device includes a Windows operating system.
The mobile flat plate is adopted as mobile equipment, so that the mobile flat plate is convenient to carry and operate, and a Windows operating system is adopted, so that the functions are complete, the compatibility is good, and the operation is convenient.
As shown in fig. 1, according to an embodiment of the present disclosure, a method of determining a side slope range includes: and clicking a screen through the mobile equipment to select the datum points, clicking three or more representative datum points, and determining the approximate range of the slope.
The method for determining the slope range by selecting the reference points is convenient to operate, convenient to implement and high in accuracy.
The selected reference points are uniformly distributed, so that the accuracy of the slope range can be improved.
As shown in fig. 1, according to an embodiment of the present disclosure, a method for screening and filtering outlier slope points includes:
determining the range of the horizontal and vertical coordinates X, Y according to the clicked reference points, randomly combining n groups of (X, Y) data point pairs in the range, and extracting elevation information on an oblique photography model according to the oblique photography model and the n groups of (X, Y) data point pairs to obtain a first set of three-dimensional coordinates of the n groups of (X, Y) data point pairs; preferred values of n may be 150, 200, 250.
And then automatically screening out outlier slope points by using an isolated forest algorithm, and filtering to obtain a second set of retained three-dimensional coordinates of the slope points.
The unsupervised algorithm-the isolated forest algorithm realizes automatic separation of outliers by training iTrees and then integrating the iTrees into a forest, and the specific implementation mode is shown in the following algorithm 1 and algorithm 2.
TABLE 1 Algorithm 1
TABLE 2 Algorithm 2
The isolated forest algorithm is adopted to automatically screen out the outlier slope points, so that the screening accuracy can be ensured, and the screening efficiency can be improved.
As shown in fig. 1, according to an embodiment of the present disclosure, a method of solving a slope occurrence includes:
arbitrarily selecting three reference points from the point selection, and preliminarily drawing up a slope surface equation S0：
A0X+B0Y+C0Z+1＝0 (1)
Wherein: a. the0、B0、C0Is the coefficient of the plane equation;
x, Y, Z is a spatial rectangular three-dimensional coordinate;
assuming that the final slope surface equation fitted by all the selected datum points is S:
AX+BY+CZ+1＝0 (2)
wherein: A. b, C are coefficients of a plane equation;
x, Y, Z is a spatial rectangular three-dimensional coordinate;
then S0And S satisfies the following relationship:
A＝A0+A*,B＝B0+B*,C＝C0+C* (3)
wherein: a. the*、B*、C*Respectively representing A, B, C corresponding error coefficients;
bringing formula (3) into formula (2) to obtain a reference point (X)i,Yi,Zi) The fitted slope surface equation S is:
S＝A*Xi+B*Yi+C*Zi+A0Xi+B0Yi+C0Zi+1＝A*Xi+B*Yi+C*Zi+Ni (4)
wherein: n is a radical ofiThe error caused by the equation introduced by the ith slope point is represented as a constant term;
reference point is not less than 3, using least square method, using point (X)i,Yi,Zi) I-0-n-1, the slope equation should be fitted, so that equation (5) is minimized,
to make it
The system of equations is obtained as:
solving and calculating a linear equation set (6), solving to obtain an error coefficient, and finally obtaining a slope equation S:
AfX+BfY+CfZ+1＝0 (7)
wherein: a. thef、Bf、CfRepresents the final coefficients of slope equation X, Y, Z;
the trend of the slope is obtained, the intersection line of the slope S and the plane Z is 0, and the calculation formula is as follows:
(1) when B is presentfWhen the value is 0, the trend of the slope is as follows:
θ1＝90 (8)
(2) when B is presentfWhen not equal to 0, the trend calculation formula of the side slope is as follows:
θ1＝arctan(-Af/Bf) (9)
and (3) calculating the inclination angle of the side slope, wherein the calculation formula is as follows:
and (3) solving the slope inclination, considering the quadrant of the projection of the slope equation normal vector on the X-Y plane, and specifically calculating the following steps:
(1) when C is presentf＞0：
(2) When C is presentf＜0：
The slope equation is fitted through the datum points by adopting a least square method, and then the attitude of the side slope is solved through the slope equation, the whole solving process is strict in logic, the accuracy of the numerical value of the attitude of the side slope is improved, and accurate basic data are provided for the later analysis of geological prospecting personnel.
As shown in fig. 1, according to an embodiment of the present disclosure, a method of determining a triangular surface closest to a slope includes:
combining the reserved slope surface points and the reserved reference points, and encrypting by using a Delaunay triangulation network;
solving the occurrence of each triangular surface through cyclic iteration;
based on all the triangular surface attitude and the side slope attitude, any attribute of the trend, the inclination and the dip angle is selected, and all the triangular surfaces and the side slopes are visualized in a thermodynamic diagram mode, so that the triangular surface closest to the side slope attitude elements is determined.
According to an embodiment of the present disclosure, a method of visualizing using thermodynamic diagrams includes: by calculating:
wherein, XfRepresenting any occurrence attribute on any triangle surface;
Xarearepresenting the corresponding occurrence attribute on the final slope surface equation;
Utilize the form of thermodynamic diagram, visual to all triangle faces and side slopes, can make things convenient for geological survey personnel to discern and confirm and the analysis to the occurrence of side slope, improved the efficiency and the accuracy that the occurrence of side slope was measured, had fine practicality.
The embodiments of the present disclosure have been described above, but the embodiments are only for illustrative purposes and are not intended to limit the scope of the present disclosure. Although the embodiments are described separately above, this does not mean that the measures in the embodiments cannot be used in advantageous combination. The scope of the disclosure is defined by the appended claims and equivalents thereof. Various alternatives and modifications can be devised by those skilled in the art without departing from the scope of the disclosure, and these alternatives and modifications are intended to be within the scope of the disclosure.
Claims (10)
1. A high-precision side slope occurrence measuring method comprises the following steps:
carrying out aerial photography on the area to be measured by using an unmanned aerial vehicle aerial photography technology, and acquiring geological attribute information;
converting the acquired geological attribute information into a tilted photography model by using a graphic workstation, and then importing the tilted photography model into mobile equipment;
determining a natural side slope to be measured in the oblique photography model, randomly selecting at least three reference points on the natural side slope, and determining a side slope range;
randomly generating n slope point data within the slope range, and then screening and filtering the outlier slope points to obtain reserved slope points, wherein n is an integer greater than or equal to 1;
the datum points adopt a least square method to fit a slope surface equation, and the occurrence of the side slope is solved through the slope surface equation;
encrypting the reserved slope points by using a Delou internal triangular network, circularly traversing all the triangular surfaces, and solving the occurrence of all the triangular surfaces;
and comparing the side slope attitude with the attitude of all the triangular faces, and determining the triangular face closest to the attitude elements of the side slope.
2. The surveying method according to claim 1, wherein the oblique photography model comprises a three-dimensional real scene model formed by superposition processing of images acquired by aerial photography by the unmanned aerial vehicle and elevation data and capable of representing a field terrain environment.
3. The measurement method according to claim 1, the determination method of the side slope range comprising: and clicking a screen through the mobile equipment to select the reference points, clicking three or more representative reference points, and determining the approximate range of the slope.
4. The method of measurement according to claim 3, the fiducial points being uniformly distributed.
5. The method of measurement according to claim 1, wherein the method of screening and filtering outlier slope points comprises:
determining the range of a horizontal coordinate X, Y and a vertical coordinate X, Y according to the clicked reference points, randomly combining n groups of (X, Y) data point pairs in the range, and extracting elevation information on the oblique photography model according to the oblique photography model and the n groups of (X, Y) data point pairs to obtain a first set of three-dimensional coordinates of the n groups of (X, Y) data point pairs;
and then automatically screening the outlier slope points by using an isolated forest algorithm, and filtering to obtain the retained second set of three-dimensional coordinates of the slope points.
6. The method of measuring of claim 1, the method of solving for the slope attitude comprising:
selecting three reference points from the point selection, and preliminarily drawing up a slope surface equation S0：
A0X+B0Y+C0Z+1＝0 (1)
Wherein: a. the0、B0、C0Is the coefficient of the plane equation;
x, Y, Z is a spatial rectangular three-dimensional coordinate;
assuming that the final slope surface equation fitted by all the selected datum points is S:
AX+BY+CZ+1＝0 (2)
wherein: A. b, C are coefficients of a plane equation;
x, Y, Z is a spatial rectangular three-dimensional coordinate;
then S0And S satisfies the following relationship:
A＝A0+A*,B＝B0+B*,C＝C0+C* (3)
wherein: a. the*、B*、C*Respectively representing A, B, C corresponding error coefficients;
bringing formula (3) into formula (2) to obtain the reference point (X)i,Yi,Zi) The fitted slope surface equation S is:
S＝A*Xi+B*Yi+C*Zi+A0Xi+B0Yi+C0Zi+1＝A*Xi+B*Yi+C*Zi+Ni (4)
wherein: n is a radical ofiThe error caused by the equation introduced by the ith slope point is represented as a constant term;
reference point is not less than 3, using least square method, using point (X)i,Yi,Zi) I-0-n-1, the slope equation should be fitted, so that equation (5) is minimized,
to make it
The system of equations is obtained as:
solving and calculating a linear equation set (6), solving to obtain an error coefficient, and finally obtaining a slope equation S:
AfX+BfY+CfZ+1＝0 (7)
wherein: a. thef、Bf、CfFinal coefficients representing slope equation X, Y, Z;
the trend of the slope is obtained, the intersection line of the slope S and the plane Z is 0, and the calculation formula is as follows:
(1) when B is presentfWhen the value is 0, the trend of the slope is as follows:
θ1＝90 (8)
(2) when B is presentfWhen not equal to 0, the trend calculation formula of the side slope is as follows:
θ1＝arctan(-Af/Bf) (9)
and (3) calculating the inclination angle of the side slope, wherein the calculation formula is as follows:
and (3) solving the slope inclination, considering the quadrant of the projection of the slope equation normal vector on the X-Y plane, and specifically calculating the following steps:
(1) when C is presentf＞0：
(2) When C is presentf＜0：
7. The method of measuring of claim 1, the method of determining the triangular face closest to the slope comprising:
combining the reserved slope surface points and the reference points, and encrypting by using a Delaunay triangulation network;
solving the occurrence of each triangular surface through cyclic iteration;
based on all the triangular surface attitude and the side slope attitude, any attribute of the trend, the inclination and the dip angle is selected, and all the triangular surfaces and the side slopes are visualized in a thermodynamic diagram mode, so that the triangular surface closest to the side slope attitude elements is determined.
8. The measurement method of claim 7, the method of visualizing using thermodynamic diagrams comprising: by calculating:
wherein, XfRepresenting any occurrence attribute on any triangle surface;
Xarearepresenting the corresponding occurrence attribute on the final slope surface equation;
9. The measurement method of claim 1, the mobile device comprising a mobile tablet.
10. The measurement method of claim 1, wherein the system employed by the mobile device comprises a Windows operating system.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110783751.5A CN113375644B (en) | 2021-07-12 | 2021-07-12 | High-precision measuring method for side slope occurrence |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110783751.5A CN113375644B (en) | 2021-07-12 | 2021-07-12 | High-precision measuring method for side slope occurrence |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113375644A true CN113375644A (en) | 2021-09-10 |
| CN113375644B CN113375644B (en) | 2023-06-02 |
Family
ID=77581734
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110783751.5A Active CN113375644B (en) | 2021-07-12 | 2021-07-12 | High-precision measuring method for side slope occurrence |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113375644B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113821850A (en) * | 2021-09-29 | 2021-12-21 | 中水北方勘测设计研究有限责任公司 | Steep slope oblique photography model geological boundary optimization method using point migration technology |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106846476A (en) * | 2017-02-10 | 2017-06-13 | 中国电建集团成都勘测设计研究院有限公司 | Rock block stability fast appraisement method based on three-dimensional live and stereographic projection |
| CN107146283A (en) * | 2017-05-10 | 2017-09-08 | 南京师范大学 | A kind of Automated Partition Method of Stratified Rock Slope type |
| CN108489403A (en) * | 2018-06-08 | 2018-09-04 | 绍兴文理学院 | The quick fine obtaining value method of open mine side slope ROCK MASS JOINT occurrence based on 3 D laser scanning |
| CN108801221A (en) * | 2018-06-08 | 2018-11-13 | 绍兴文理学院 | The quick fine obtaining value method of open mine side slope ROCK MASS JOINT scale based on digital photogrammetry |
| CN108830317A (en) * | 2018-06-08 | 2018-11-16 | 绍兴文理学院 | The quick fine obtaining value method of open mine side slope ROCK MASS JOINT occurrence based on digital photogrammetry |
-
2021
- 2021-07-12 CN CN202110783751.5A patent/CN113375644B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106846476A (en) * | 2017-02-10 | 2017-06-13 | 中国电建集团成都勘测设计研究院有限公司 | Rock block stability fast appraisement method based on three-dimensional live and stereographic projection |
| CN107146283A (en) * | 2017-05-10 | 2017-09-08 | 南京师范大学 | A kind of Automated Partition Method of Stratified Rock Slope type |
| CN108489403A (en) * | 2018-06-08 | 2018-09-04 | 绍兴文理学院 | The quick fine obtaining value method of open mine side slope ROCK MASS JOINT occurrence based on 3 D laser scanning |
| CN108801221A (en) * | 2018-06-08 | 2018-11-13 | 绍兴文理学院 | The quick fine obtaining value method of open mine side slope ROCK MASS JOINT scale based on digital photogrammetry |
| CN108830317A (en) * | 2018-06-08 | 2018-11-16 | 绍兴文理学院 | The quick fine obtaining value method of open mine side slope ROCK MASS JOINT occurrence based on digital photogrammetry |
Non-Patent Citations (2)
| Title |
|---|
| 贾曙光: "无人机摄影测量在高陡边坡地质调查中的应用", 《岩土力学》 * |
| 邓一: "基于非接触测量技术的隧道开挖对杨家山边坡稳定性影响分析", 《中国优秀博硕学位论文全文数据库（硕士）基础科学辑，》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113821850A (en) * | 2021-09-29 | 2021-12-21 | 中水北方勘测设计研究有限责任公司 | Steep slope oblique photography model geological boundary optimization method using point migration technology |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113375644B (en) | 2023-06-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Carvajal-Ramírez et al. | Virtual reconstruction of damaged archaeological sites based on Unmanned Aerial Vehicle Photogrammetry and 3D modelling. Study case of a southeastern Iberia production area in the Bronze Age | |
| CN106846476B (en) | Rock block stability rapid evaluation method based on three-dimensional live-action and red-horizontal projection | |
| Yakar et al. | Performance of photogrammetric and terrestrial laser scanning methods in volume computing of excavtion and filling areas | |
| CN114998536A (en) | Model generation method and device based on novel basic mapping and storage medium | |
| CN102800127A (en) | Light stream optimization based three-dimensional reconstruction method and device | |
| CN109961510B (en) | High-cut-slope geological rapid recording method based on three-dimensional point cloud reconstruction technology | |
| Gold et al. | A terrestrial lidar-based workflow for determining three–dimensional slip vectors and associated uncertainties | |
| CN107709926A (en) | The mobile ground mapping of automation | |
| Viana et al. | Algorithms for extraction of structural attitudes from 3D outcrop models | |
| Moemen et al. | 3-D reconstruction and measurement system based on multimobile robot machine vision | |
| CN111006645A (en) | Unmanned aerial vehicle surveying and mapping method based on motion and structure reconstruction | |
| CN113375644A (en) | High-precision side slope occurrence measuring method | |
| Xiao et al. | Monitoring excavation slope stability using drones | |
| WO2022104251A1 (en) | Image analysis for aerial images | |
| Yilmaz | Close range photogrammetry in volume computing | |
| Feng | Novel methods for 3-D semi-automatic mapping of fracture geometry at exposed rock faces | |
| TW202001615A (en) | Method for using point cloud to determine rock discontinuity attitude which includes acquiring proper point cloud; and determining an attitude of rock discontinuity by using the point cloud according to geomagnetism and geographic directions | |
| Zahm et al. | Use of unmanned aerial vehicles (UAVs) to create digital outcrop models: An example from the Cretaceous Cow Creek Formation, Central Texas | |
| Stylianidis et al. | A digital close-range photogrammetric technique for monitoring slope displacements | |
| CN115393537A (en) | Precision evaluation system and method for three-dimensional visual modeling of power transmission channel | |
| Smith et al. | Absolute and exterior orientation using linear features | |
| CN115018973A (en) | Low-altitude unmanned-machine point cloud modeling precision target-free evaluation method | |
| JP2813767B2 (en) | Topographic map creation method and device | |
| CN113670266A (en) | Technology for measuring real estate title by utilizing unmanned aerial vehicle oblique photography | |
| Huang et al. | An Innovative Approach of Evaluating the Accuracy of Point Cloud Generated by Photogrammetry-Based 3D Reconstruction |
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 |