CN107146283A - A kind of Automated Partition Method of Stratified Rock Slope type - Google Patents
A kind of Automated Partition Method of Stratified Rock Slope type Download PDFInfo
- Publication number
- CN107146283A CN107146283A CN201710323875.9A CN201710323875A CN107146283A CN 107146283 A CN107146283 A CN 107146283A CN 201710323875 A CN201710323875 A CN 201710323875A CN 107146283 A CN107146283 A CN 107146283A
- Authority
- CN
- China
- Prior art keywords
- point
- slope
- key element
- side slope
- type
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 67
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000005192 partition Methods 0.000 title claims abstract description 13
- 230000011218 segmentation Effects 0.000 claims abstract description 5
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000012216 screening Methods 0.000 claims description 3
- 238000013459 approach Methods 0.000 claims description 2
- 230000000052 comparative effect Effects 0.000 claims description 2
- 125000004122 cyclic group Chemical group 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 239000000284 extract Substances 0.000 description 2
- 238000003780 insertion Methods 0.000 description 2
- 230000037431 insertion Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
Classifications
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Remote Sensing (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
- Processing Or Creating Images (AREA)
Abstract
The present invention discloses a kind of Automated Partition Method of Stratified Rock Slope type, comprises the following steps：(1) ridge line is pre-processed, then based on ridge line segmentation massif border, to form slope unit；(2) side slope occurrence and attitude of rocks information are calculated；(3) according to the Stratified Rock Slope criteria for classifying, side slope is classified.The present invention can divide the side slope type of Stratified Rock automatically, and automaticity is high, and objectivity is strong, disclosure satisfy that the Type division needs of a wide range of Stratified Rock Slope.
Description
Technical field
The invention belongs to geographical information technology application technology, and in particular to a kind of automatic division of Stratified Rock Slope type
Method.
Background technology
Side slope is nature or the slope that is artificially formed, is not only one of geological environment most basic in Human dried bloodstains,
Or most common project in engineering construction.According to different criteria for classifications, side slope can be divided into different types：Press
Formation lithology, can be divided into soil-slope and rock side slope；By strata structure, can be divided into monolithic, bulk, stratiform, cataclastic shape and
Big type of granular media shape structure 5 etc..And in various types of rock side slopes, Stratified Rock Slope is considered as that distribution is most wide,
The simplest class side slope of structure type, therefore the correlative study carried out to it is also more.
The edge slope structure parameter larger to Stratified Rock Slope stability influence, mainly includes：Side slope and formation strike it
Between angle, angle and formation dip etc. between tendency., first, can be by bedded rock according to side slope and the relation of formation strike
Matter side slope is divided into vertical slope (trend is vertical), tangential slope (moving towards oblique), nearly horizontal slope (moving towards parallel)；Further, according to side
The relation that slope and rock stratum are inclined to, nearly horizontal side slope can be divided into horizontal slope (about 0 ° of formation dip), dip slope (tendency is identical),
Anti- slope of inclining (tendency is opposite).In summary, Stratified Rock Slope can substantially be divided into following 5 type：1. horizontal slope；2. forward
Slope；3. tangential slope；4. vertical slope；5. anti-slope of inclining.
Automation divides the premise of Stratified Rock Slope type, is to obtain side slope data boundary and attitude of rocks data.Side
Slope data boundary can be based on dem data source, and by extracting the generation of massif border, (Li Mengyuan, Li Anbo, Xie Xianli wait a kind of mountains of
Body border extraction method [P] China:201710049100.7)；And attitude of rocks data can be based on terrain and geologic map, from
(Li Anbo, Chen Ying, Yao Mengmeng wait a kind of adaptive judgements of the earth's surface exposure attitude of rocks of for the calculating of the adaptation completion attitude of rocks
Method [P] China:201610912428.2).
The Type division of Stratified Rock Slope, is Correct Analysis its Mechanism of Deformation And Failure, carries out estimation of stability and use
The precondition of rational support design method.
However, the correlative study of existing side slope Type division, time-consuming generally using artificial or semi-automatic division methods
Arduously, and automaticity is low, it is difficult to the need for meeting a wide range of side slope Type division.
The content of the invention
Goal of the invention：It is an object of the invention to solve the deficiencies in the prior art, there is provided a kind of Stratified Rock side
The Automated Partition Method of slope type, the present invention is based on massif boundary face data, ridge line data and occurrence data, and then realizes layer
The automatic division of shape rock side slope type.
Technical scheme：A kind of Automated Partition Method of Stratified Rock Slope type of the present invention, specifically includes following successively
Step：
(1) ridge line is pre-processed, then based on ridge line segmentation massif border, to form slope unit；
(2) side slope occurrence and attitude of rocks information are calculated；
(3) according to the Stratified Rock Slope criteria for classifying, side slope is classified.
Further, the step (1) specifically includes procedure below：
(1-1) loading massif boundary face figure layer, ridge line figure layer and occurrence point figure layer, respectively obtain massif boundary face key element
Set Mount={ mi| i=1,2 ..., n, ridge line line feature set RidgeL={ rli| i=1,2 ..., a } and occurrence point
Elements combination Occur={ oi| i=1,2 ..., b }；
Wherein, miI-th of face key element is represented, n is the number of face key element, rliI-th of ridge line line feature is represented, a is mountain
The number of crestal line line feature, oiI-th of occurrence point key element is represented, its attribute includes tendency [Tendency] and inclination angle
[dipAngle], b is the number of occurrence point key element；
(1-2) calculates ridge line trend：For each ridge line line feature rli, read in point set Pi={ pj| j=
1,2 ..., t }, t is the number of a key element, and adjacent 2 points p is taken successivelyjAnd pj+1, calculate pjObject-point points to pj+1The side of object-point
Parallactic angle, obtains azimuth set Azimuth={ azi| i=1,2 ..., t-1 }, azi∈[0°,180°)；Again by maximum frequency method
Screen azimuth：
(1-3) two ends extend ridge line；
(1-4) splits massif into side slope：For each face key element mi, search itself and line feature rliUnder intersection point former point
Mark e1And e2：The intersection point of face key element and line feature is to refer to rl hereiniTwo point p from beginning to end1And pt；By face key element miRead in point set Ci=
{cj| j=1,2 ..., s }, adjacent 2 points c is taken successivelyj、cj+1If, first point p1In line segmentOn, then e1=j；If tail point ptOnline
SectionOn, then e2=j；Compare e again1With e2Size, s for point key element number.
Further, the azimuthal method of screening is in the step (1-2)：
(1-2-1) is by parameter T and 360 ° come computer azimuth angular spacing N, N=360 °/T；Wherein, T is 360 approximate number (T
Conventional value is 36 or 72)；
(1-2-2) by the azimuthal distribution in set Azimuth it is each it is uniform it is interval [0 °, N), [N, 2N), [2N,
3N)、…、[(T-1)×N,360°)；
(1-2-3) calculates azimuthal average in the maximum interval of frequency, and regard the average as rliSubstantially trend
Φi, Φi∈[0°,180°)。
Further, the step (1-3) includes following 2 links：
(1-3-1) extends ridge line tail end：Line taking key element rl successivelyiThree point point p of tail endt-2、pt-1And pt, calculate adjacent
Distance between two pointsWithAverage d；Inverted order travels through point set P againi, two point point p are taken successivelyjAnd pj-1, wherein j is circulation
Variable, initial value is t, makes following judgement：
If a) this 2 points in massif boundary face key element miIt is interior, then remove point pj；
If b) only point pj-1In face key element miIt is interior, then remove point pj, and by line segment at j+1 positionsWith face key element mi's
Intersection point insertion point set Pi, and terminate circulation；
If c) this 2 points in face key element miIt is interior, then the coordinate p'(X, Y of newly-increased point are calculated by formula (1)), if point p' exists
Face key element miIt is interior, point p' is added into point set Pi, j=j+1；Conversely, by line segmentWith face key element miIntersection point add point set Pi, and
Terminate circulation；
Wherein, pjAnd p .xj.y it is respectively point pjAbscissa and ordinate, ΦiFor line feature rliSubstantially trend；
(1-3-2) extends ridge line head end：Step is step (1-3-1) with the method for extension ridge line tail end.
Further, e in the step (1-4)1With e2Size comparative approach be：
(1-4-1) is if e1<e2, traversal point set Ci, take point cj, carry out following judgement and operate：
If I, j<e1Or j>e2, by point cjAdd point set RSlope；
If II, j=e1, by point cjAnd point set PiIn have an addition point collection RSlope；
If III, e1<j<e2, by point cjAdd point set LSlope；
If IV, j=e2, by point cjAnd point set PiIn backward point add point set LSlope；
(1-4-2) is if e1>e2, process step is with (1-4-1)；
Point set CiTraversal is finished, and obtains left side slope point set LSlope={ lsi| i=1,2 ..., q } and right side slope point set
RSlope={ rsi| i=1,2 ..., r }, wherein, q is the number of point key element, and r is the number of point key element.
Further, the detailed process of the step (2) is：
(2-1) calculates side slope occurrence：The trend of ridge line is used as the substantially trend of left and right side slope, i.e. ρl'=ρr'=Φi；
And the tendency of left and right side slope can be calculated by formula (2)：
Wherein, ρl' and θl' be left side slope trend and tendency, ρr' and θr' refer to right side slope trend and tendency；
(2-2) calculates the attitude of rocks：
The attitude of rocks calculation procedure on slope is as follows on the left of (2-2-1)：
A, traversal occurrence point elements combination Occur, screen all occurrence points in the point set LSlope of left side slope；
B, the tendency for reading occurrence point in the slope of left side, the tendency in the maximum interval of frequency is searched by maximum frequency method, then
Corresponding occurrence point is further searched for, the dominant orientation point set LOccP={ lp on left side slope are obtainedi| i=1,2 ..., h }, h is excellent
The number of gesture occurrence point；
C, traversal point set LOccP, calculate tendency average a little, inclination angle average respectively as left side slope rock stratum tendency
θl, inclination angle δl；ρ is moved towards by formula (3) calculating left side slope rock stratuml, ρl∈[0°,180°)；
(2-2-2) calculates the occurrence (ρ of right side slope rock stratumr,θr,δr), step is with (2-2-1), ρr、θrAnd δrIt is right side respectively
Trend, tendency and the inclination angle of slope rock stratum.
Further, the step (3) specifically includes：
(3-1) divides Stratified Rock Slope：
(3-1-1) sets the type attribute on left side slope：Left side slope and the absolute value of left side slope formation strike difference are calculated respectively
α, tendency difference absolute value β, if β>180 °, then β=360 °-β；If moving towards uniformity tolerance for u, identical or opposite tolerance is inclined to
For w, left side slope type is set by following Rule of judgment：
If i, 60 °<α<120 °, then it is " vertical slope " to remember type attribute Type；
If ii, u≤α≤60 ° or 120 °≤α≤180 °-u, note type attribute Type is " tangential slope "；
If iii, 0 °≤α<U or 180 ° of-u<α<180 °, 3 kinds of situations can be divided into：When left side slope inclination angle meets 0 °≤δl≤5°
When, then it is " horizontal slope " to remember type attribute Type；When the poor absolute value β of tendency meets 0 °≤β<During w, then remember that type attribute Type is
" dip slope "；When β meets 180 ° of-w<During β≤180 °, then it is " anti-slope of inclining " to remember type attribute Type；
The type attribute setting steps on slope are with (3-1-1) on the right side of (3-1-2), until the Type division of all left and right side slopes
Complete；
Side slope face key element after division is saved as new figure layer by (3-2).
Beneficial effect：The present invention is based on massif boundary face data, ridge line data and attitude of rocks data, is drawn with reference to side slope
Minute mark is accurate, obtains side slope occurrence and attitude of rocks information, realizes the automatic division of Stratified Rock Slope type；Compared to existing
Manually-semi-automatic division methods, automaticity is high, and objectivity is strong, disclosure satisfy that the type of a wide range of Stratified Rock Slope is drawn
Dividing needs.
Brief description of the drawings
Fig. 1 is the flow chart of the inventive method；
Fig. 2 is massif boundary face figure layer and ridge line figure layer schematic diagram in embodiment；
Fig. 3 is stratum line chart layer and occurrence point figure layer schematic diagram in embodiment；
Fig. 4 is the ridge line figure layer schematic diagram after extension processing in embodiment；
Fig. 5 is the side slope face figure layer schematic diagram formed behind segmentation massif border in embodiment；
Fig. 6 is the side slope face figure layer schematic diagram of classified types in embodiment.
Embodiment
Technical solution of the present invention is described in detail below, but protection scope of the present invention is not limited to the implementation
Example.
Embodiment：
The experimental data of this example comes from Mount Lushan Mountainous Area of North dem data, and grid size is 2162 × 2260.Based on DEM
Data source, extracts ridge line (such as Fig. 2), valley route and contour line data；Basin boundary face data are extracted with anti-landform DEM again,
It is then based on ridge line and valley route automatically extracts massif boundary face data (such as Fig. 2)；It is final to combine contour and geology dignity
Data, calculate attitude of rocks information (such as Fig. 3).Below in conjunction with the accompanying drawings, and by describing an automatic division Stratified Rock Slope
The example of type, to further illustrate the effect of the present invention.
As shown in figure 1, in the present embodiment, the Automated Partition Method of Stratified Rock Slope type comprises the following steps：
Step (1), is specifically included：
(1-1) as shown in Figures 2 and 3, loading massif boundary face figure layer, ridge line figure layer and occurrence point figure layer are obtained respectively
To massif boundary face elements combination Mount={ mi| i=1,2 ..., 8, ridge line line feature set RidgeL={ rli| i=1,
2 ..., 8 } and occurrence point elements combination Occur={ oi| i=1,2 ..., 507 }；
(1-2) calculates ridge line trend：With line feature rl1Exemplified by, read in point set P1={ pj| j=1,2 ...,
924 }, adjacent 2 points p is taken successively1、p2, calculate p1Object-point points to p2The azimuth of object-point is 68 °；Adjacent 2 points p is taken again2、
p3, calculate p2Object-point points to p3The azimuth of object-point is 73 °；……；Treat point set P1Last 2 points p923、p924Azimuth
Calculating is finished, and obtains azimuth set Azimuth={ 68 °, 73 ° ..., 44 ° } (azimuth number is 923)；Again by maximum frequency
Number method screening azimuth：
1) by parameter T=36 and 360 ° come computer azimuth angle at intervals of：360 °/36=10 °；
2) by the azimuthal distribution in set Azimuth it is each it is uniform it is interval [0 °, 10 °), [10 °, 20 °), [20 °,
30°)、…、[350°,360°)；
3) the maximum interval of statistics frequency for [50 °, 60 °) (frequency is 113), take 54 ° of the average at the interval inner orientation angle
It is used as rl1Substantially trend；
(1-3) two ends extend ridge line：The step includes following 2 links：
1. ridge line tail end is extended：Line taking key element rl1The point of tail end 3 p922、p923And p924, calculate distance between two pointsAverage be 7.71205216370225；Inverted order travels through point set P again1, take two point p924、p923, 2 points exist
Face key element m1It is interior, then the coordinate p' of newly-increased point is calculated by formula (1)：
p'(405511.09571363218,3278760.43547032245)；
Wherein, p' is in face key element m1It is interior, p' is added into point set P1, j=925；Point set P is taken again1The point of tail end two p925、p924, two
Point is in face key element m1, then the coordinate of newly-increased point is calculated, the point is not in face key element m1It is interior, by the point and point p925The line segment of composition with
Face key element m1Intersection point add point set P1, terminate circulation；
2. ridge line head end is extended：Line taking key element rl1The point of head end 3 p1、p2And p3, calculate distance between two points 's
Average is 5.9827602150087351；Travel through point set P1, take two point p1、p2, 2 points in face key element m1It is interior, then by formula (1)
Calculate the coordinate p' of newly-increased point：
p'(401633.04610488308,3273824.3561243751)；
Wherein, p' is in face key element m1It is interior, p' is inserted into point set P1Head end, j is constant；Point set P is taken again1The point of head end two p1、p2, two
Point is in face key element m1It is interior, then the coordinate of newly-increased point is calculated, the point is in face key element m1It is interior, the point is inserted into point set P1Head end, j is not
Become；Continue aforesaid operations, until the 7th time takes point set P1During head end 2,2 points in face key element m1It is interior, but newly-increased point is not or not face
Key element m1It is interior, by newly-increased point and point p1The line segment of composition and face key element m1Intersection point insertion P1Head end, terminates circulation；
After the completion of extension operation, line feature rl1Include point set P1Point key element number be 933.Its in RidgeL to be gathered
After remaining line feature is finished by step (1-3), the ridge line line feature figure layer after generation extension processing, as shown in Figure 4；
(1-4) splits massif into side slope：For face key element m1, search itself and line feature rl1Intersection point (i.e. rl1Head and the tail point p1、
p933) former point subscript e1、e2：By face key element m1Read in point set C1={ cj| j=1,2 ..., 2078 }, take successively at adjacent 2 points,
It is c when adjacent 2 points29、c30When, tail point p933In line segmentOn, then e2=29；It is c when adjacent 2 points880、c881When, first point
p1In line segmentOn, then e1=880；Compare e1With e2Size understand：
e1>e2, traversal point set C1, by c1-c28These points add point set LSlope；Then will point c29Point set LSlope is added,
Simultaneously by point set P1In point backward add point set LSlope；Again by c30-c879These points add point set RSlope；Again by c880
Point set RSlope is added, while by point set P1In point sequentially add point set RSlope；Finally by c881-c2078These addition points
Collect LSlope；
Point set C1Traversal is finished, and obtains left side slope point set LSlope (point key element number is 2193) and right side slope point set
RSlope (point key element number is 1751)；
Its lap key element in the elements combination Mount of face is carried out the segmentation of massif by step (1-4), treats whole divisions
After left and right side side slope, generation side slope face key element figure layer, as shown in Figure 5.
Step (2), is specifically included：
(2-1) calculates side slope occurrence：Ridge line trend can as left and right side slope substantially trend.For
Massif boundary face key element m1And ridge line line feature rl1, the trend of its left and right side slope of correspondence is ρl'=ρr'=Φ1
=
54 °, and the tendency of left and right side slope can be calculated θ by formula (2)l'=144 °, θr'=324 °；Remaining massif
The occurrence of the corresponding left and right side slope of boundary face key element is calculated as stated above；
(2-2) calculates the attitude of rocks：The step includes following 2 links：
1. the attitude of rocks calculation procedure on left side slope is as follows：
A, traversal occurrence point elements combination Occur, screen all occurrence points (totally 167 in the point set LSlope of left side slope
Individual occurrence point)；
B, the tendency for reading occurrence point in the slope of left side, the tendency in the maximum interval of frequency is searched by maximum frequency method, then
Corresponding occurrence point is further searched for, the dominant orientation point set LOccP={ lp on left side slope are obtained1,lp2,…,lp32}；
C, traversal point set LOccP, calculate tendency average a little, inclination angle average respectively as left side slope rock stratum tendency
θl=313 °, inclination angle δl=75 °；ρ is moved towards by formula (3) calculating left side slope rock stratuml=43 °；
2. the attitude of rocks calculation procedure on right side slope is as follows：
A, traversal occurrence point elements combination Occur, screen all occurrence points (totally 69 in the point set RSlope of right side slope
Occurrence point)；
B, the tendency for reading occurrence point in the slope of right side, the tendency in the maximum interval of frequency is searched by maximum frequency method, then
Corresponding occurrence point is further searched for, the dominant orientation point set ROccP={ rp on right side slope are obtained1,rp2,…,rp9}；
C, traversal point set ROccP, calculate tendency average a little, inclination angle average respectively as right side slope rock stratum tendency
θr=114 °, inclination angle δr=80 °；ρ is moved towards by formula (3) calculating left side slope rock stratumr=24 °；
Continue to calculate the attitude of rocks of remaining left and right side slope by step (2-2), until the attitude of rocks of all side slopes is calculated
Finish.
Step (3), is specifically included：
(3-1) divides Stratified Rock Slope：The step includes following 2 links：
1. the type attribute on left side slope is set：Calculate respectively the absolute value α of left side slope and left side slope formation strike difference=
11 °, absolute value β=169 ° of tendency difference；If it is 10 ° to move towards uniformity tolerance, be inclined to identical or opposite tolerance, by judgement bar
Part is known：10°<α<60 °, meet condition ii, therefore the type attribute Type=on left side slope " tangential slope "；
2. the type attribute on right side slope is set：Calculate respectively the absolute value α of right side slope and right side slope formation strike difference=
30 °, absolute value β=210 ° of tendency difference, β>180 °, then -210 ° of β=360 °=150 °, are known by Rule of judgment：10°<α<60 °,
Meet condition ii, therefore the type attribute Type=on right side slope " tangential slope "；
Continue executing with aforesaid operations to set the type of remaining left and right side slope, until completing the type attribute of all side slopes
Assignment；
(3-2) by the side slope face key element after division as shown in fig. 6, save as new figure layer.
It can be seen from above-described embodiment that this patent combination side slope criteria for classifying, and using edge slope structure parameter as incision
Point, from geoscience applications angle, forms a kind of automation division methods of Stratified Rock Slope type, efficiently precisely and automatic
Division degree is high.
Claims (7)
1. a kind of Automated Partition Method of Stratified Rock Slope type, it is characterised in that：Specifically comprise the following steps successively：
(1) ridge line is pre-processed, then based on ridge line segmentation massif border, to form slope unit；
(2) side slope occurrence and attitude of rocks information are calculated；
(3) according to the Stratified Rock Slope criteria for classifying, side slope is classified.
2. the Automated Partition Method of Stratified Rock Slope type according to claim 1, it is characterised in that：The step
(1) procedure below is specifically included：
(1-1) loading massif boundary face figure layer, ridge line figure layer and occurrence point figure layer, respectively obtain massif boundary face elements combination
Mount={ mi| i=1,2 ..., n, ridge line line feature set RidgeL={ rli| i=1,2 ..., a } and occurrence point key element
Set Occur={ oi| i=1,2 ..., b }；
Wherein, miI-th of face key element is represented, n is the number of face key element, rliI-th of ridge line line feature is represented, a is ridge line
The number of line feature, oiI-th of occurrence point key element is represented, its attribute includes tendency [Tendency] and inclination angle [dipAngle], b
For the number of occurrence point key element；
(1-2) calculates ridge line trend：For each ridge line line feature rli, read in point set Pi={ pj| j=1,
2 ..., t }, t is the number of a key element, and adjacent 2 points p is taken successivelyjAnd pj+1, calculate pjObject-point points to pj+1The orientation of object-point
Angle, obtains azimuth set Azimuth={ azi| i=1,2 ..., t-1 }, azi∈[0°,180°)；Again by maximum frequency method sieve
Select azimuth：
(1-3) two ends extend ridge line；
(1-4) splits massif into side slope：For each face key element mi, search itself and line feature rliThe subscript e of intersection point former point1
And e2：The intersection point of face key element and line feature is to refer to rl hereiniTwo point p from beginning to end1And pt；By face key element miRead in point set Ci={ cj|
J=1,2 ..., s }, adjacent 2 points c is taken successivelyj、cj+1If, first point p1In line segmentOn, then e1=j；If tail point ptIn line segmentOn, then e2=j；Compare e again1With e2Size, s for point key element number.
3. the Automated Partition Method of Stratified Rock Slope type according to claim 2, it is characterised in that：The step
The azimuthal method of screening is in (1-2)：
(1-2-1) is by parameter T and 360 ° come computer azimuth angular spacing N, N=360 °/T；Wherein, T is 360 approximate number；
(1-2-2) by the azimuthal distribution in set Azimuth it is each it is uniform it is interval [0 °, N), [N, 2N), [2N, 3N) ...,
[(T-1)×N,360°)；
(1-2-3) calculates azimuthal average in the maximum interval of frequency, and regard the average as rliSubstantially move towards Φi,
Φi∈[0°,180°)。
4. the Automated Partition Method of Stratified Rock Slope type according to claim 2, it is characterised in that：The step
(1-3) includes procedure below：
(1-3-1) extends ridge line tail end：Line taking key element rl successivelyiThree point point p of tail endt-2、pt-1And pt, calculate at adjacent 2 points
Between distanceWithAverage d；Inverted order travels through point set P againi, two point point p are taken successivelyjAnd pj-1, wherein j is cyclic variable,
Initial value is t, makes following judgement：
If a) this 2 points in massif boundary face key element miIt is interior, then remove point pj；
If b) only point pj-1In face key element miIt is interior, then remove point pj, and by line segment at j+1 positionsWith face key element miIntersection point
Insert point set Pi, and terminate circulation；
If c) this 2 points in face key element miIt is interior, then the coordinate p'(X, Y of newly-increased point are calculated by formula (1)), if point p' is in face key element
miIt is interior, point p' is added into point set Pi, j=j+1；Conversely, by line segmentWith face key element miIntersection point add point set Pi, and terminate follow
Ring；
Wherein, pjAnd p .xj.y it is respectively point pjAbscissa and ordinate, ΦiFor line feature rliSubstantially trend；
(1-3-2) extends ridge line head end：Step is step (1-3-1) with the method for extension ridge line tail end.
5. the Automated Partition Method of Stratified Rock Slope type according to claim 2, it is characterised in that：The step
E in (1-4)1With e2Size comparative approach be：
(1-4-1) is if e1<e2, traversal point set Ci, take point cj, carry out following judgement and operate：
If I, j<e1Or j>e2, by point cjAdd point set RSlope；
If II, j=e1, by point cjAnd point set PiIn have an addition point collection RSlope；
If III, e1<j<e2, by point cjAdd point set LSlope；
If IV, j=e2, by point cjAnd point set PiIn backward point add point set LSlope；
(1-4-2) is if e1>e2, process step is with (1-4-1)；
Point set CiTraversal is finished, and obtains left side slope point set LSlope={ lsi| i=1,2 ..., q } and right side slope point set RSlope=
{rsi| i=1,2 ..., r }, wherein, q is the number of point key element, and r is the number of point key element.
6. the Automated Partition Method of Stratified Rock Slope type according to claim 1, it is characterised in that：The step
(2) detailed process is：
(2-1) calculates side slope occurrence：The trend of ridge line is used as the substantially trend of left and right side slope, i.e. ρl'=ρr'=Φi；And
The tendency of left and right side slope can be calculated by formula (2)：
Wherein, ρl' and θl' be left side slope trend and tendency, ρr' and θr' refer to right side slope trend and tendency；
(2-2) calculates the attitude of rocks：
The attitude of rocks calculation procedure on slope is as follows on the left of (2-2-1)：
A, traversal occurrence point elements combination Occur, screen all occurrence points in the point set LSlope of left side slope；
B, the tendency for reading occurrence point in the slope of left side, search the tendency in the maximum interval of frequency, then enter one by maximum frequency method
Step searches corresponding occurrence point, obtains the dominant orientation point set LOccP={ lp on left side slopei| i=1,2 ..., h }, h produces for advantage
The number of shape point；
C, traversal point set LOccP, calculate tendency average a little, inclination angle average respectively as left side slope rock stratum tendency θl, incline
Angle δl；ρ is moved towards by formula (3) calculating left side slope rock stratuml, ρl∈[0°,180°)；
(2-2-2) calculates the occurrence (ρ of right side slope rock stratumr,θr,δr), step is with (2-2-1), ρr、θrAnd δrIt is right side slope rock respectively
Trend, tendency and the inclination angle of layer.
7. the Automated Partition Method of Stratified Rock Slope type according to claim 1, it is characterised in that：The step
(3) specifically include：
(3-1) divides Stratified Rock Slope：
(3-1-1) sets the type attribute on left side slope：The left side slope absolute value α poor with left side slope formation strike is calculated respectively, inclined
To the absolute value β of difference, if β>180 °, then β=360 °-β；If moving towards uniformity tolerance for u, it is w to be inclined to identical or opposite tolerance,
By following Rule of judgment, left side slope type is set：
If i, 60 °<α<120 °, then it is " vertical slope " to remember type attribute Type；
If ii, u≤α≤60 ° or 120 °≤α≤180 °-u, note type attribute Type is " tangential slope "；
If iii, 0 °≤α<U or 180 ° of-u<α<180 °, 3 kinds of situations can be divided into：When left side slope inclination angle meets 0 °≤δlAt≤5 °, then
It is " horizontal slope " to remember type attribute Type；When the poor absolute value β of tendency meets 0 °≤β<During w, then remember type attribute Type for " forward
Slope "；When β meets 180 ° of-w<During β≤180 °, then it is " anti-slope of inclining " to remember type attribute Type；
The type attribute setting steps on slope are with (3-1-1) on the right side of (3-1-2), until the Type division of all left and right side slopes is completed；
Side slope face key element after division is saved as new figure layer by (3-2).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710323875.9A CN107146283B (en) | 2017-05-10 | 2017-05-10 | Automatic dividing method for layered rock slope types |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710323875.9A CN107146283B (en) | 2017-05-10 | 2017-05-10 | Automatic dividing method for layered rock slope types |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107146283A true CN107146283A (en) | 2017-09-08 |
| CN107146283B CN107146283B (en) | 2020-05-19 |
Family
ID=59777935
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710323875.9A Active CN107146283B (en) | 2017-05-10 | 2017-05-10 | Automatic dividing method for layered rock slope types |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107146283B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108829742A (en) * | 2018-05-24 | 2018-11-16 | 南京师范大学 | A kind of slope pattern Type division method |
| CN113375644A (en) * | 2021-07-12 | 2021-09-10 | 天津大学 | High-precision side slope occurrence measuring method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101538861A (en) * | 2009-04-21 | 2009-09-23 | 中国科学院武汉岩土力学研究所 | Highway slope stability grading evaluation method |
| CN105701848A (en) * | 2016-01-14 | 2016-06-22 | 南京师范大学 | Automatic generation method of stratum boundary map layer |
| CN105894587A (en) * | 2016-04-01 | 2016-08-24 | 南京师范大学 | Ridge line and valley line filtering method based on rule constraints |
| CN105956066A (en) * | 2016-04-28 | 2016-09-21 | 南京师范大学 | Automated identification method for fold landform type |
-
2017
- 2017-05-10 CN CN201710323875.9A patent/CN107146283B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101538861A (en) * | 2009-04-21 | 2009-09-23 | 中国科学院武汉岩土力学研究所 | Highway slope stability grading evaluation method |
| CN105701848A (en) * | 2016-01-14 | 2016-06-22 | 南京师范大学 | Automatic generation method of stratum boundary map layer |
| CN105894587A (en) * | 2016-04-01 | 2016-08-24 | 南京师范大学 | Ridge line and valley line filtering method based on rule constraints |
| CN105956066A (en) * | 2016-04-28 | 2016-09-21 | 南京师范大学 | Automated identification method for fold landform type |
Non-Patent Citations (1)
| Title |
|---|
| 姚蒙蒙等: "地质界线图层的自动化生成方法", 《地球信息科学学报》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108829742A (en) * | 2018-05-24 | 2018-11-16 | 南京师范大学 | A kind of slope pattern Type division method |
| CN108829742B (en) * | 2018-05-24 | 2021-12-17 | 南京师范大学 | Slope form type dividing method |
| CN113375644A (en) * | 2021-07-12 | 2021-09-10 | 天津大学 | High-precision side slope occurrence measuring method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107146283B (en) | 2020-05-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103592681B (en) | A kind of seismic image tracing of horizons method based on signal classification | |
| CN107807387B (en) | Acquisition methods when seismic first break neural network based is walked | |
| CN106194181A (en) | Intelligent work surface coal-rock interface identification method based on geologic data | |
| CN104809266A (en) | Spline based accurate predicating method for face coal seam occurrence condition | |
| CN104459782A (en) | Horizon chromatography inversion earthquake speed modeling method with thin layers used and modeling unit | |
| CN105701848A (en) | Automatic generation method of stratum boundary map layer | |
| CN107146283A (en) | A kind of Automated Partition Method of Stratified Rock Slope type | |
| CN107766978A (en) | The intelligent optimization method of skewed pattern | |
| CN105403913A (en) | Pre-stack depth migration method and device | |
| CN108957530B (en) | A kind of crack automatic testing method based on Acceleration Algorithm in Seismic Coherence Cube slice | |
| CN105607119B (en) | Near-surface model construction method and static correction value acquiring method | |
| CN103901467A (en) | Method for tracking positions of three-dimensional seismic data | |
| CN105116447A (en) | A geological river direction discrimination method based on curvature-abnormal stripes | |
| CN109655890A (en) | A kind of shallow mid-deep strata joint chromatography inversion speed modeling method of Depth Domain and system | |
| CN112329099B (en) | Intelligent matching method for main tunnel of railway mountain tunnel based on custom database | |
| CN106777117A (en) | A kind of automatic identifying method of horizontal stratum tectonic landform | |
| CN105089659B (en) | A kind of Conglomerate Reservoir permeable unit recognition methods | |
| CN116152512A (en) | Height measuring and calculating method based on building shadow restoration | |
| CN115294290A (en) | Method for extracting data in geological mineral exploration | |
| CN109117534A (en) | A kind of fracture and crack prediction method based on geological map | |
| CN114608476A (en) | Intelligent analysis and extraction method for three-dimensional point cloud structural surface of complex rock mass | |
| CN111754421B (en) | Improved guiding filtering three-dimensional scattered point cloud rapid fairing method | |
| CN101980055B (en) | Digital modeling method for seismic exploration | |
| CN102096102B (en) | Digital modeling method for seismic exploration | |
| CN111737867A (en) | Method for rapidly updating three-dimensional entity model of ore body |
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 |