CN106709439A - Monocline structure landform automatic identification method - Google Patents
Monocline structure landform automatic identification method Download PDFInfo
- Publication number
- CN106709439A CN106709439A CN201611164730.0A CN201611164730A CN106709439A CN 106709439 A CN106709439 A CN 106709439A CN 201611164730 A CN201611164730 A CN 201611164730A CN 106709439 A CN106709439 A CN 106709439A
- Authority
- CN
- China
- Prior art keywords
- model
- key element
- place
- summit
- center
- 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
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000012216 screening Methods 0.000 claims abstract description 19
- 239000011435 rock Substances 0.000 claims description 9
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 9
- 238000003325 tomography Methods 0.000 description 3
- 230000019771 cognition Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 241000282898 Sus scrofa Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V30/00—Character recognition; Recognising digital ink; Document-oriented image-based pattern recognition
- G06V30/40—Document-oriented image-based pattern recognition
- G06V30/42—Document-oriented image-based pattern recognition based on the type of document
- G06V30/422—Technical drawings; Geographical maps
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/74—Image or video pattern matching; Proximity measures in feature spaces
- G06V10/75—Organisation of the matching processes, e.g. simultaneous or sequential comparisons of image or video features; Coarse-fine approaches, e.g. multi-scale approaches; using context analysis; Selection of dictionaries
- G06V10/751—Comparing pixel values or logical combinations thereof, or feature values having positional relevance, e.g. template matching
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Multimedia (AREA)
- Artificial Intelligence (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Databases & Information Systems (AREA)
- Health & Medical Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Image Analysis (AREA)
Abstract
The invention discloses a monocline structure landform automatic identification method. The method comprises the following steps: 1, extracting a bedrock stratum object, according to the rule of adjacent, inclination and tendency consistency, establishing a relation for an extracted adjacent object, and drafting an adjacent ARG graph; 2, based on the adjacent ARG graph, according to the adjacent relation, constructing a scene model, and simplifying the scene model; and 3, according to a fold screening rule, screening the simplified scene model, and reserving a model of an incomplete fold structure, i.e., a monocline structure. According to the method, automatic identification of a monocline structure landform is realized.
Description
Technical field
Digital geological map is based on the present invention relates to geographical information technology application field, more particularly to one kind, is tied by space
Structure pattern match, the method for automatic identification monoclinal strate tectonic landform.
Background technology
The rock stratum of level is influenceed run-off the straight by diastrophism, at a certain angle with horizontal plane shape, and rock stratum is to same side
To inclination, referred to as monoclinal structure.Common monoclinal strate tectonic landform includes swine back, cuesta etc..Monocline is likely to occur
In destroyed flank of anticline, or the horizontal stratum for appearing in the surrounding of disrupted dome structure, the periphery in basin, tilting
Or tomography is tilted etc..Monoclinal strate is mainly shown as three kinds of forms：One is a wing of fold tectonic；Two is fault tectonic
One disk；What uneven lifting or descending motion in three regions were formed.
Traditional conformation identification is mainly artificial operation, and its workload is big, lasts length, and be limited to the technology of making drawing personnel
And the level of understanding.The identification of monoclinal strate tectonic landform depends on the result of geologic information or ground observation, in geologic map
Upper lookup monocline is mainly to look for the serial tilted stratum in fold or tomography side.Old principal columns of a hall etc. (the old principal columns of a hall, Li Anbo, Yao Mengmeng,
Fold geomorphic type automatic identification [J] Earth Information Science journals of spatial structure model matching, 2016,18 are based on Deng
(11)) the method matched by spatial structure model, is realized based on digital geological dignity data to anticline, synclinal structure landforms
The method for carrying out automatic identification.But the tactic pattern only accounts for structure symmetrically, repeatedly, for monoclinal strate tectonic landform
It is not identified well then.
The formation of monoclinal strate tectonic landform has close ties with the changing of physiognomy of surrounding and structural evolution, studies monocline rock
Layer tectonic landform can provide effectively support for geologic map cognition, structural evolution reasoning, be tectonic knot, Parametric three-dimensional modeling etc.
The premise of work.
The content of the invention
Goal of the invention：The present invention for prior art exist problem, there is provided a kind of monoclinal strate tectonic landform it is automatic
Recognition methods.The method uses digital geological map, and consistent rule is inclined to for monoclinal stratum, with the side that tactic pattern is matched
Method, realizes in basement rock mountain region area, to the automatic identification of monoclinal strate tectonic landform substantially.Compared to traditional manual identified
Method, realizes the automatic identification based on geologic map, and also knowing figure for the geologic map interpreting blueprints of layman provides help.
Technical scheme：The automatic identifying method of monoclinal strate tectonic landform of the present invention is comprised the following steps：
(1) bedrock formation object (not including plutone) is extracted, according to adjacent, inclination and the rule of tendency uniformity,
The contiguous object opening relationships obtained to extraction, draws adjacent ARG figures (attributed relational graph, attributed relational
graph)；
(2) based on adjacent ARG figures, model of place is built according to syntople, and model of place is simplified；
(3) model of place after simplification is screened according to fold screening rule, retains the mould of incomplete pleated structure
Type, as monoclinal strate are constructed.
Wherein, step (1) is specifically included：
Whether (1-1) sets angle threshold α as judging two orientation angles or angle of strike difference in tolerance interval
Interior foundation, general value [22.5,90)；Tilt threshold β is set used as judging whether rock stratum reaches the foundation of oblique standard, one
As value [15,90), specifically can voluntarily be determined according to user's request；
(1-2) loads the stratal surface key element figure layer of vector format, obtains all face elements combination Stra={ si| i=1,2,
3,…,n}；Wherein, siI-th face key element is represented, face key element includes numbering attribute Id, stratum age attribute Age, occurrence tendency
Attribute OccT and occurrence dip angle attribute OccA, n are the quantity of face key element；
(1-3) calculates the center of mass point of face key element respectively, obtains stratal surface key element center of mass point elements combination OriGrav={ ogi
(xogi,yogi) | i=1,2,3 ..., n }；Wherein, ogiRepresent i-th center of mass point key element of face key element, (xogi,yogi) it is ogi
Coordinate, center of mass point key element inherits the numbering attribute Id of face key element, stratum age attribute Age, occurrence tendency attribute OccT and product
Shape dip angle attribute OccA；
(1-4) center of mass point key element is screened：To set OriGrav, according to stratum age attribute Age, non-Quaternary Period, non-is chosen
The center of mass point key element on the stratum of plutone, creates center of mass point elements combination Grav={ gi(xgi,ygi) | i=1,2,3 ...,
c}；Wherein, giBe retain i-th center of mass point key element as (xgi,ygi) it is giCoordinate, c be retain center of mass point key element
Quantity；
(1-5) creates stratum object relationship：For center of mass point elements combination Grav, center of mass point key element g is takeniAnd gj, j=
1,2,3 ..., c, j>I, judges whether to meet following condition：
A) center of mass point key element giAnd gjOccurrence dip angle attribute OccA be more than or equal to tilt threshold β；
B) center of mass point key element giAnd gjStratum age attribute Age differ；
C) center of mass point key element giAnd gjCorresponding stratal surface key element siAnd sjIt is adjacent；
D) center of mass point key element giAnd gjThe difference of occurrence tendency attribute OccT be no more than angle threshold α；
E) center of mass point key element giAnd gjThe trend of place straight line and center of mass point key element giOr gjOccurrence tendency attribute OccT
Difference is no more than angle threshold α；
For the center of mass point key element g for meeting all conditions aboveiAnd gjEstablishment relation rk<gi,gj>, and store relation rkIn
It is left
Object giWith right object gjNumbering attribute Id；
The relation that (1-6) is completed between all center of mass point key elements is created, and obtains set of relationship Rel={ rk<gi,gj>| k=1,
2,3 ..., m }, m is the quantity of the relation for creating；
(1-7) draws side according to center of mass point Element Drawing vertex set according to the relation between key element, then will based on center of mass point
Element set Grav and set of relationship Rel completes adjacent ARG figures.
Wherein, step (2) is specifically included：
(2-1) is based on adjacent ARG figures, and model of place is created to all connection sides by condition of syntople, obtains scene
Model set OriM={ OMi| i=1,2,3 ..., p }, OMiI-th model of place is represented, p is the quantity of model of place, scene
Model OMiIncluding vertex set and line set.Side r for example to there is syntopleu<gi, gj>With rv<gi,gk>(u ≠ v) builds
Model of place OM, model of place OM include vertex set OMGrav={ gi, g, gk, and line set OMRel={ ru,rv, and
Calculate and recorded summit gi, g, gkThe quantity on side be respectively 2,1,1；
(2-2) searches end points：In the vertex set of model of place, there is simplified end points two kinds of with popular endpoints
End points, searches simplify end points and popular endpoints according to the following steps successively：
1) summit for meeting following all conditions is searched as simplified end points：
A) quantity for crossing the side on the summit is more than or equal to 2；
B) there is stratum age attribute identical situation in the Extreme points set being connected with the summit；
C) there is side of at least one pair of angle less than 90 ° in the line set being connected with the summit；
2) summit of condition of the quantity equal to 1 on the side for meeting connection is searched as popular endpoints；
(2-3) recursive reconstruction model：All simplified end points and unlabelled popular endpoints are pressed as starting endpoint successively
Carry out recurrence according to following steps, in the absence of simplify end points model of place directly using popular endpoints as starting endpoint according to following
Step carries out recurrence；
1) by starting endpoint giSummit chain set Link is added, by connected other summits gijIt is considered as giChild, j
=1,2,3 ..., ti, tiIt is and end points giThere is the summit quantity of connection：
If 2) Link length<2, then by summit gijAdd Link, step 4 is carried out to it) judgement；If Link length>=
2, then step 3 is carried out to it) judge；
3) opposite vertexes gijCarry out rule judgment, j=1,2,3 ..., ti：If summit gij1. condition is met simultaneously to simplify end
Point is not labeled；2. do not exist in Link；3. in Link in the orientation angles and Link on the side of most latter two summit composition
Point to g in last summitijOrientation angles difference be less than predetermined angle threshold value α, then by summit gijAdd Link, mark top
Point gij, continue step 4) judge；Conversely, then continuing to next child summit gi(j+1)The judgement of this step is carried out, until giIt is all
Child judges to finish, and jumps to step 5)；
4) opposite vertexes gijCarry out end points judgement：If summit gijTo simplify end points or popular endpoints, then the Link terminates, root
Gather reconstruction model M according to current Link, by current vertex gijRemoved from Link；Conversely, then by summit gijIt is considered as gi, to its child
Sub- gijRepeat step 2) judge；
5) current parent summit giAll children traversal finish, by summit giRemoved from Link, continued to next top
Point g(i+1) or next starting endpoint repeat step 2) judge, to the last an end points judges to finish, and recurrence terminates；
(2-4) is screened by following rule to the model of place that recurrence is obtained：
1) two summits and the model of place on side are only included in deleting madel；
2) it is multiple models of place that beginning or end is rebuild for same simplified end points, is rejected by following principle：
If a) the vertex set number of elements of model is identical, the corresponding stratum age attribute Age difference phase of set inner vertex
Together, then the big model of place of corresponding stratal surface area is retained；
If b) vertex set of model is incomplete same, retain the more model of place of set element；
(2-5) enters line renumbering to the model of place after screening, obtains model set AllM={ M1,M2,M3,…,
Mq, q is model of place quantity.
Wherein, step (3) is specifically included：
(3-1) model Vonoroi element subdivisions：Voronoi element subdivisions are carried out according to the model of place that step (2) is obtained
(Wang is newborn, and Liu Jiyuan, the village is generous, and Approximation Algorithm for Constructing Voronoi Diagrams with General Generators [J] the geographical sciences for waiting to be based on GIS enter
Exhibition, 2004,23 (4)), the Voronoi polygonal elements after subdivision inherit pattern number attribute；
(3-2) carries out fold screening to model：For the model of place M in the Voronoi polygonal elements that adjoin each otheri
With MjIt is compared, judges that two models of place constitute complete fold if following condition is met simultaneously, rejects the two models：
A) model of place MiWith MjVertex set it is consistent；
B) according to consistent summit sequentially, the trend of two models of place is opposite；
The model of place that (3-3) completes to retain after fold screening is considered as monocline model, and monocline model is obtained after full figure screening
Set FinM={ FMh| h=1,2,3 ..., w }, w is monocline model quantity, FMhRepresent h-th monocline model, FMhIncluding summit
Set FMGravhWith line set FMRelh；
(3-4) carries out stratal surface key element merging to each model of place in the model set FinM after screening, obtains right
The monocline face key element answered, forms monocline face elements combination Poly={ pr|r=1,2,3 ..., w }, prIt is model of place FMrList
Inclined-plane key element.
Beneficial effect：Compared with prior art, its remarkable advantage is the present invention：The present invention uses digital geological map, for
Consistent rule is inclined in monoclinal stratum, in the method that tactic pattern is matched, realizes in basement rock mountain region area, to monoclinal strate structure
Make the automatic identification of landforms.Compared to traditional manual identified method, the automatic identification based on geologic map is realized, also for non-
The geologic map interpreting blueprints of professional knows figure and provides help.The formation of monoclinal strate tectonic landform and the changing of physiognomy of surrounding and structure
Making evolution has close ties, and research monoclinal strate tectonic landform can provide effectively support for geologic map cognition, structural evolution reasoning,
It is the premise of the work such as tectonic knot, Parametric three-dimensional modeling.
Brief description of the drawings
Fig. 1 is the flow chart of the inventive method；
Fig. 2 is embodiment geology dignity figure layer schematic diagram data；
Fig. 3 is embodiment center of mass point key element figure layer schematic diagram data；
Fig. 4 is that embodiment screening center of mass point key element builds relation schematic diagram；
Fig. 5 is example (a) model of place figure layer (b) vertex set attribute (c) line set attribute schematic diagram；
Fig. 6 is the model of place exemplary plot under complex situations
Fig. 7 is model schematic diagram data after simplified screening；
Fig. 8 is model data Voronoi subdivision schematic diagrames；
Fig. 9 is complete fold model schematic；
Figure 10 is the monocline face figure layer schematic diagram data of identification.
Specific embodiment
Further illustrate below in conjunction with the accompanying drawings and by describing an example for automatic identification monoclinal strate tectonic landform
Effect of the invention.Flow chart of the method for the present invention is as shown in figure 1, comprise the following steps：
(1) adjacent ARG figures generation
Step 11：It is 45 ° to set angle threshold α, and it is 15 ° to set tilt threshold β；
Step 12：The shp formatted datas of the geology dignity figure layer of the present embodiment selection Pattern of Zijinshan are example, such as Fig. 2
It is shown, the figure stratal surface key element figure layer of vector format is loaded, obtain all face elements combination Stra={ si| i=1,2,
3,…,78}.Attribute OccT and occurrence inclination angle category are inclined in the numbering attribute Id of storage face key element, stratum age attribute Age, occurrence
Property OccA.Data represent that occurrence is inclined to using 16 orientation models, convert thereof into numerical value, here, with the interval center in orientation
Orientation angles are substituted；
Step 13：The center of mass point of face key element is calculated respectively, obtains stratal surface key element center of mass point elements combination OriGrav=
{ogi(xogi,yogi) | i=1,2,3 ..., 78 }.Center of mass point key element inherits numbering attribute Id, the stratum age attribute of face key element
Age, occurrence tendency attribute OccT and occurrence dip angle attribute OccA；
Step 14：Screening center of mass point, according to stratum age attribute Age, chooses the non-Quaternary Period, the stratum of non-intruding rock mass
Center of mass point key element, obtains center of mass point object set Grav={ gi(xgi,ygi) | i=1,2,3 ..., 62 }, the barycenter that screening is obtained
Point object figure layer is as shown in Figure 3；
Step 15：Create stratum object relationship.For center of mass point object set Grav, barycenter point object g is takeniAnd gj(j
=1,2,3 ..., 62, j>I) judge whether to meet following condition：
A) center of mass point key element giAnd gjOccurrence dip angle attribute OccA be more than or equal to tilt threshold β；
B) center of mass point key element giAnd gjStratum age attribute Age differ；
C) center of mass point key element giAnd gjCorresponding stratal surface key element siAnd sjIt is adjacent；
D) center of mass point key element giAnd gjThe difference of occurrence tendency attribute OccT be no more than angle threshold α；
E) center of mass point key element giAnd gjThe trend of place straight line and center of mass point key element giOr gjOccurrence tendency attribute OccT
Difference is no more than angle threshold α；
For the center of mass point key element g for meeting all conditions aboveiAnd gjEstablishment relation rk<gi,gj>, and store relation rkIn
It is left
Object giWith right object gjNumbering attribute Id, obtain set Rel={ rk<gi,gj>| k=1,2,3 ..., 21 }, root
According to matter
Heart point object and set of relationship draw summit and side, complete ARG figures, as shown in Figure 4.
(2) model of place builds and simplifies
Step 21：Based on ARG figures, model of place is created by condition of neighbouring relations.With model OM5As a example by (such as Fig. 5 (a)
It is shown), its vertex set OMGrav5Including 8 object OMGrav5={ g23, g24, g27, g28, g30, g31, g32, g33, attribute is such as
Shown in Fig. 5 (b), its line set OMRel5Including 7 side OMRel5={ r9<g23, g28>, r10<g23, g30>, r11<g23, g31>, r12
<g23, g33>, r13<g27, g28>, r14<g31, g32>, r15<g24, g27>, shown in attribute such as Fig. 5 (c).Summit connection in a model
The quantity on side be respectively 4,1,2,2,1,2,1,1；
Step 22：Full figure carries out model of place structure, obtains model of place set OriM={ OMi| i=1,2,3 ...,
8}；
Step 23：Model of place is simplified.There may be is influenceed by later stage tomography and causes model of place structure to answer
Miscellaneous situation (as shown in Figure 6), it is necessary to simplify to model of place, with model of place OM5As a example by, model of place simplified
Journey is illustrated as shown in step 24~26.
Step 24：Search end points.
1) in model of place OM5Vertex set OMGrav5In, summit g23Meet following condition, then summit g23It is simplification
End points：
A) summit g is crossed23Side quantity be 4, more than 2；
B) with summit g23In the vertex set of connection, g31Withg33Stratum age attribute it is consistent
C) summit g is included23All line sets in there is angle less than 90 ° of side；
2) secondly, condition of the quantity according to the side passed through equal to 1 finds popular endpoints, then summit g24, g30, g32, g33For
Popular endpoints；
Step 25：Recursive reconstruction model.A summit chain set Link is created, for the simplified end points g of model of place23Enter
Row recurrence.End points g will be simplified first23Link is added, its other summits g connected by side23j(j=1,2,3 ..., 4) respectively
It is g28、g30、g31、g33, then recurrence is carried out as follows：
1) g is now only included in Link23, length is less than 2, then by summit g28Link is added, g is marked28；
2) opposite vertexes g28Carry out end points judgement, g28Neither simplify end points nor popular endpoints, then find it and pass through side
Other summits g of connection28j(j=1,2) it is respectively g23, g27；
3) g is now included in Link23、g28, length is equal to 2, then opposite vertexes g23Carry out rule judgment.If summit gijSimultaneously
Condition is met 1. to simplify end points or not being labeled；2. do not exist in Link；3. the side that most latter two summit constitutes in Link
Orientation angles and Link in last summit point to gijOrientation angles difference be less than predetermined angle threshold value α, then by summit
gijLink is added, summit g is markedij.For summit g23, it is to simplify end points, is present in Link, therefore is unsatisfactory for condition 2.,
The point is then skipped, and then to g27Carry out rule judgment；
4) for summit g27, its is unmarked, and does not exist in Link, calculates g23Point to g28Angle be 341.3 °, g28
Point to g27Angle be 90.3 °, two angles differ by more than 45 °, are unsatisfactory for condition 3., then skip the point, now g28Child
Judgement is finished；
5) by summit g28Removed from Link, g is now only included in Link23, length is less than 2, then by summit g30Add
Link, marks g30；
6) opposite vertexes g30Carry out end points judgement, g30It is popular endpoints, then the Link terminates, newly-built model M, MGrav=
{g23, g30, MRel={ r10<g23, g30>}.By g30Removed from Link, g is now only included in Link23, length be less than 2, then
By summit g31Link is added, g is marked31；
7) opposite vertexes g31Carry out end points judgement, g31Neither simplify end points, nor popular endpoints, then find it and pass through
Other summits g of side connection31j(j=1,2) it is respectively g23, g32；
8) g is now included in Link23、g31, length is equal to 2, then opposite vertexes g23Carry out rule judgment.For summit g23, its
To simplify end points, exist in Link and be unsatisfactory for condition 2., then skip the point, and then to g32Carry out rule judgment；
9) for summit g32, its is unmarked, and does not exist in Link, calculates g23Point to g31Angle be 251.4 °, g31
Point to g32Angle be 256.0 °, two angles difference is less than 45 ° of threshold value, then by summit g32Link is added, g is marked32；
10) opposite vertexes g32Carry out end points judgement, g32It is popular endpoints, then the Link terminates, newly-built model M, MGrav=
{g23, g31, g32, MRel={ r11<g23, g31>, r14<g31, g32>}.By g32Removed from Link, now, g31Child judge
Terminate；
11) by summit g31Removed from Link, g is now only included in Link23, length is less than 2, then by summit g33Add
Link, marks g33；
12) opposite vertexes g33Carry out end points judgement, g33It is popular endpoints, then the Link terminates, newly-built model M, MGrav=
{g23, g33, MRel={ r12<g23, g33>}.By g33Removed from Link, now, g23Child judge terminate；
13) unlabelled popular endpoints are traveled through.According to judgment rule and flow, to not labeled popular endpoints
It is g24Child judged, model M is obtained, including vertex set MGrav={ g24, g27、g28, line set MRel=
{r13<g27, g28>, r15<g24, g27>}；
14) popular endpoints g30、g32、g33All it is labeled, no longer travels through, current scene model recursion terminates.Exemplary scene mould
Type OriM5It is as follows that simplified completion obtains 4 models：
M1={ MGrav1, MRel1, MGrav1={ g23, g30, MRel1={ r10<g23, g30>}
M2={ MGrav2, MRel2, MGrav2={ g23, g31, g32, MRel2={ r11<g23, g31>, r14<g31, g32>}
M3={ MGrav3, MRel3, MGrav3={ g23, g33, MRel3={ r12<g23, g33>}
M4={ MGrav4, MRel4, MGrav4={ g24, g27, g28, MRel4={ r13<g27, g28>, r15<g24, g27>}
Step 26：By rule to model of place OM54 simplified models are screened, wherein, model M1With M3Only include
Two summits and a line, then deleting madel M1With M3, finally retain M2With M4Two models；
Step 27：Full figure obtains 5 models after simplifying screening, is marked after renumbeing as shown in Figure 7.
(3) monocline pattern match and screening
Step 31：Voronoi element subdivisions are carried out according to the model after combing, the polygon after subdivision inherits pattern number
Attribute, the Voronoi polygon key elements of generation are as shown in Figure 8；
Step 32：Model to being adjoined each other in Voronoi diagram is screened.Voronoi for adjoining each other is polygon
Model of place M in shape unitiWith MjIt is compared, judges that two model of place compositions are complete if following condition is met simultaneously
Fold, rejects the two models：
A) model of place MiWith MjVertex set it is consistent；
B) according to consistent summit sequentially, two trends of model of place are opposite (as shown in Figure 9)；
Do not exist complete fold in example, without the monocline object for needing to reject, therefore, obtain final model set FinM
={ FMh| h=1,2,3 ..., 5 }, the quantity with model M is consistent；
Step 33：Model FM to meeting monoclineh, the corresponding stratal surface feature object of vertex set is merged, newly
Build a monocline face key element ph.Stratal surface key element merging is carried out to the model set FinM after screening, monocline face key element collection is obtained
Close Poly={ ph| h=1,2,3 ..., 5 }, as shown in Figure 10.
Above disclosed is only a kind of preferred embodiment of the invention, it is impossible to the right model of the present invention is limited with this
Enclose, therefore the equivalent variations made according to the claims in the present invention, still belong to the scope that the present invention is covered.
Claims (4)
1. a kind of automatic identifying method of monoclinal strate tectonic landform, it is characterised in that the method is comprised the following steps：
(1) bedrock formation object is extracted, according to adjacent, inclination and the rule of tendency uniformity, the contiguous object obtained to extraction
Opening relationships, draws adjacent ARG figures；
(2) based on adjacent ARG figures, model of place is built according to syntople, and model of place is simplified；
(3) model of place after simplification is screened according to fold screening rule, retains the model of incomplete pleated structure, i.e.,
For monoclinal strate is constructed.
2. the automatic identifying method of monoclinal strate tectonic landform according to claim 1, it is characterised in that：Step (1) has
Body includes：
Whether (1-1) sets angle threshold α as judging two orientation angles or angle of strike difference in tolerance interval
Foundation, sets tilt threshold β as judging whether rock stratum reaches the foundation of oblique standard；
(1-2) loads the stratal surface key element figure layer of vector format, obtains all face elements combination Stra={ si| i=1,2,3 ...,
n}；Wherein, siI-th face key element is represented, face key element includes numbering attribute Id, stratum age attribute Age, occurrence tendency attribute
OccT and occurrence dip angle attribute OccA, n are the quantity of face key element；
(1-3) calculates the center of mass point of face key element respectively, obtains stratal surface key element center of mass point elements combination OriGrav={ ogi
(xogi,yogi) | i=1,2,3 ..., n }；Wherein, ogiRepresent i-th center of mass point key element of face key element, (xogi,yogi) it is ogi
Coordinate, center of mass point key element inherits the numbering attribute Id of face key element, stratum age attribute Age, occurrence tendency attribute OccT and product
Shape dip angle attribute OccA；
(1-4) center of mass point key element is screened：To set OriGrav, according to stratum age attribute Age, the non-Quaternary Period, non-intruding are chosen
The center of mass point key element on the stratum of rock mass, creates center of mass point elements combination Grav={ gi(xgi,ygi) | i=1,2,3 ..., c }；Its
In, giBe retain i-th center of mass point key element as (xgi,ygi) it is giCoordinate, c be retain center of mass point key element quantity；
(1-5) creates stratum object relationship：For center of mass point elements combination Grav, center of mass point key element g is takeniAnd gj, j=1,2,
3 ..., c, j>I, judges whether to meet following condition：
A) center of mass point key element giAnd gjOccurrence dip angle attribute OccA be more than or equal to tilt threshold β；
B) center of mass point key element giAnd gjStratum age attribute Age differ；
C) center of mass point key element giAnd gjCorresponding stratal surface key element siAnd sjIt is adjacent；
D) center of mass point key element giAnd gjThe difference of occurrence tendency attribute OccT be no more than angle threshold α；
E) center of mass point key element giAnd gjThe trend of place straight line and center of mass point key element giOr gjOccurrence tendency attribute OccT difference
No more than angle threshold α；
For the center of mass point key element g for meeting all conditions aboveiAnd gjEstablishment relation rk<gi,gj>, and store relation rkA middle left side is right
As giWith right object gjNumbering attribute Id；
The relation that (1-6) is completed between all center of mass point key elements is created, and obtains set of relationship Rel={ rk<gi,gj>| k=1,2,
3 ..., m }, m is the quantity of the relation for creating；
(1-7) draws side, then based on center of mass point key element collection according to center of mass point Element Drawing vertex set according to the relation between key element
Close Grav and set of relationship Rel and complete adjacent ARG figures.
3. the automatic identifying method of monoclinal strate tectonic landform according to claim 1, it is characterised in that：Step (2) has
Body includes：
(2-1) is based on adjacent ARG figures, and model of place is created to all connection sides by condition of syntople, obtains model of place
Set OriM={ OMi| i=1,2,3 ..., p }, OMiI-th model of place is represented, p is the quantity of model of place, model of place
OMiIncluding vertex set and line set；
(2-2) searches end points：In the vertex set of model of place, there is simplified end points and the two kinds of end of popular endpoints
Point, searches simplify end points and popular endpoints according to the following steps successively：
1) summit for meeting following all conditions is searched as simplified end points：
A) quantity for crossing the side on the summit is more than or equal to 2；
B) there is stratum age attribute identical situation in the Extreme points set being connected with the summit；
C) there is side of at least one pair of angle less than 90 ° in the line set being connected with the summit；
2) summit of condition of the quantity equal to 1 on the side for meeting connection is searched as popular endpoints；
(2-3) recursive reconstruction model：Successively using all simplified end points and unlabelled popular endpoints as starting endpoint according to
Lower step carries out recurrence, in the absence of simplify end points model of place directly using popular endpoints as starting endpoint according to following steps
Carry out recurrence；
1) by starting endpoint giSummit chain set Link is added, by connected other summits gijIt is considered as giChild, j=1,
2,3,…,ti, tiIt is and end points giThere is the summit quantity of connection：
If 2) Link length<2, then by summit gijAdd Link, step 4 is carried out to it) judgement；If Link length>=2, then
Step 3 is carried out to it) judge；
3) opposite vertexes gijCarry out rule judgment, j=1,2,3 ..., ti：If summit gijSimultaneously meet condition 1. for simplify end points or
It is not labeled；2. do not exist in Link；3. the orientation angles on the side of most latter two summit composition are last with Link in Link
Point to g in one summitijOrientation angles difference be less than predetermined angle threshold value α, then by summit gijLink is added, summit g is markedij,
Continue step 4) judge；Conversely, then continuing to next child summit gi(j+1)The judgement of this step is carried out, until giAll children
Judgement is finished, and jumps to step 5)；
4) opposite vertexes gijCarry out end points judgement：If summit gijTo simplify end points or popular endpoints, then the Link terminates, according to current
Link gathers reconstruction model M, by current vertex gijRemoved from Link；Conversely, then by summit gijIt is considered as gi, to its child gij
Repeat step 2) judge；
5) current parent summit giAll children traversal finish, by summit giRemoved from Link, continued to next summit
g(i+1) or next starting endpoint repeat step 2) judge, to the last an end points judges to finish, and recurrence terminates；
(2-4) is screened by following rule to the model of place that recurrence is obtained：
1) two summits and the model of place on side are only included in deleting madel；
2) it is multiple models of place that beginning or end is rebuild for same simplified end points, is rejected by following principle：
If a) the vertex set number of elements of model is identical, the corresponding stratum age attribute Age difference of set inner vertex is identical, then
Retain the big model of place of corresponding stratal surface area；
If b) vertex set of model is incomplete same, retain the more model of place of set element；
(2-5) enters line renumbering to the model of place after screening, obtains model set AllM={ M1,M2,M3,…,Mq, q is
Model of place quantity.
4. the automatic identifying method of monoclinal strate tectonic landform according to claim 1, it is characterised in that：Step (3) has
Body includes：
(3-1) model Vonoroi element subdivisions：Voronoi element subdivisions are carried out according to the model of place that step (2) is obtained, is cutd open
Voronoi polygonal elements after point inherit pattern number attribute；
(3-2) carries out fold screening to model：For the model of place M in the Voronoi polygonal elements that adjoin each otheriWith Mj
It is compared, judges that two models of place constitute complete fold if following condition is met simultaneously, rejects the two models：
A) model of place MiWith MjVertex set it is consistent；
B) according to consistent summit sequentially, the trend of two models of place is opposite；
The model of place that (3-3) completes to retain after fold screening is considered as monocline model, and monocline model set is obtained after full figure screening
FinM={ FMh| h=1,2,3 ..., w }, w is monocline model quantity, FMhRepresent h-th monocline model, FMhIncluding vertex set
FMGravhWith line set FMRelh；
(3-4) carries out stratal surface key element merging to each model of place in the model set FinM after screening, obtains corresponding
Monocline face key element, forms monocline face elements combination Poly={ ph|h=1,2,3 ..., w }, phIt is model of place FMhMonocline face will
Element.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611164730.0A CN106709439B (en) | 2016-12-16 | 2016-12-16 | Automatic identification method for monoclinic rock stratum structure landform |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611164730.0A CN106709439B (en) | 2016-12-16 | 2016-12-16 | Automatic identification method for monoclinic rock stratum structure landform |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106709439A true CN106709439A (en) | 2017-05-24 |
| CN106709439B CN106709439B (en) | 2020-04-03 |
Family
ID=58938876
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611164730.0A Active CN106709439B (en) | 2016-12-16 | 2016-12-16 | Automatic identification method for monoclinic rock stratum structure landform |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106709439B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111583776A (en) * | 2020-04-28 | 2020-08-25 | 南京师范大学 | Method for acquiring development time sequence of invaded rock mass |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101770026A (en) * | 2009-01-07 | 2010-07-07 | 中国科学院电子学研究所 | Method for estimating terrain by polarization interference of data of synthetic aperture radar and software thereof |
| CN102288944A (en) * | 2011-05-12 | 2011-12-21 | 西安电子科技大学 | Super-resolution height measuring method based on topographic matching for digital array meter wave radar |
| CN105956066A (en) * | 2016-04-28 | 2016-09-21 | 南京师范大学 | Automated identification method for fold landform type |
| CN106023197A (en) * | 2016-05-18 | 2016-10-12 | 南京师范大学 | Automated identification and extraction method of vertical stratum |
-
2016
- 2016-12-16 CN CN201611164730.0A patent/CN106709439B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101770026A (en) * | 2009-01-07 | 2010-07-07 | 中国科学院电子学研究所 | Method for estimating terrain by polarization interference of data of synthetic aperture radar and software thereof |
| CN102288944A (en) * | 2011-05-12 | 2011-12-21 | 西安电子科技大学 | Super-resolution height measuring method based on topographic matching for digital array meter wave radar |
| CN105956066A (en) * | 2016-04-28 | 2016-09-21 | 南京师范大学 | Automated identification method for fold landform type |
| CN106023197A (en) * | 2016-05-18 | 2016-10-12 | 南京师范大学 | Automated identification and extraction method of vertical stratum |
Non-Patent Citations (2)
| Title |
|---|
| 杨树文 等: ""基于SPOT 5图像的岩溶地貌单元自动提取方法"", 《国土资源遥感》 * |
| 陈楹 等: "基于空间结构模式匹配的褶皱地貌类型自动识别", 《地球信息科学学报》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111583776A (en) * | 2020-04-28 | 2020-08-25 | 南京师范大学 | Method for acquiring development time sequence of invaded rock mass |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106709439B (en) | 2020-04-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105574259B (en) | A kind of Urban cognition ground drawing generating method based on internet word frequency | |
| CN112598796B (en) | Method for constructing and automatically updating three-dimensional building information model based on generalized point cloud | |
| CN109147045A (en) | City section geology three-dimensional method for automatic modeling based on topological subregion | |
| CN103699654B (en) | A kind of across engineer's scale map vector network of rivers data target matching method of the same name | |
| CN105258704A (en) | Multi-scale space-time hot point path detection method based on rapid road network modeling | |
| CN116152461B (en) | Geological modeling method, device, computer equipment and computer readable storage medium | |
| CN106484943A (en) | A kind of online generation method of the log sheet based on nonlinear fitting control | |
| CN111858810B (en) | Modeling elevation point screening method for road DEM construction | |
| CN107330030A (en) | A kind of hierarchical network construction method that storage is compressed towards magnanimity road net data | |
| CN107092978A (en) | A kind of shortest path hierarchical reconfiguration planning method of the Virtual earth | |
| CN107180452A (en) | The stratified geological objects modeling method of the multi-source data coupling on the stratum containing standard | |
| CN113360989B (en) | City green roof planning method | |
| CN101924647A (en) | Local area topology rebuilding method for updating navigation road network increment | |
| Déderix | Communication networks, interactions, and social negotiation in prepalatial south-central Crete | |
| CN115757604A (en) | GDP (generalized projection) space-time evolution analysis method based on noctilucent image data | |
| CN106777117A (en) | A kind of automatic identifying method of horizontal stratum tectonic landform | |
| CN106649776B (en) | A kind of method of semi-automation synthetic vector polygon | |
| CN115169842A (en) | Urban historical street promotion strategy generation method and system | |
| CN112581558B (en) | Model construction method and system for intrusion structure in cut geological section | |
| CN106709439A (en) | Monocline structure landform automatic identification method | |
| KR101063827B1 (en) | Semi-automated Conjugated Point Pair Extraction Method for Transformation of Geometric Maps between Continuous and Digital Topographic Maps | |
| CN112116709B (en) | Terrain characteristic line processing method for improving terrain expression precision | |
| CN106023161B (en) | A kind of automatic correlation method of broken ridge line | |
| CN109190255B (en) | Three-dimensional reconstruction method for urban three-dimensional property space | |
| CN116258612A (en) | Urban land utilization distribution evaluation method based on knowledge embedding model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |