CN108875755A - A kind of horizontal stratum extraction method based on oriented parallel feature - Google Patents

A kind of horizontal stratum extraction method based on oriented parallel feature Download PDF

Info

Publication number
CN108875755A
CN108875755A CN201810521612.3A CN201810521612A CN108875755A CN 108875755 A CN108875755 A CN 108875755A CN 201810521612 A CN201810521612 A CN 201810521612A CN 108875755 A CN108875755 A CN 108875755A
Authority
CN
China
Prior art keywords
feature
boundary
formula
stratigraphic
horizontal stratum
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810521612.3A
Other languages
Chinese (zh)
Other versions
CN108875755B (en
Inventor
郑嘉康
李安波
李安营
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Normal University
Original Assignee
Nanjing Normal University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing Normal University filed Critical Nanjing Normal University
Priority to CN201810521612.3A priority Critical patent/CN108875755B/en
Publication of CN108875755A publication Critical patent/CN108875755A/en
Application granted granted Critical
Publication of CN108875755B publication Critical patent/CN108875755B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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/40—Extraction of image or video features
    • G06V10/42—Global feature extraction by analysis of the whole pattern, e.g. using frequency domain transformations or autocorrelation
    • G06V10/422—Global feature extraction by analysis of the whole pattern, e.g. using frequency domain transformations or autocorrelation for representing the structure of the pattern or shape of an object therefor
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/24—Classification techniques

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Multimedia (AREA)
  • Image Analysis (AREA)
  • Processing Or Creating Images (AREA)

Abstract

The invention discloses a kind of horizontal stratum extraction methods based on oriented parallel feature, specifically include step：Linearly stratum boundary line number evidence is generated based on planar geology diagram data；The corresponding relationship of stratigraphic boundary element and contour line feature is obtained based on geology face component attributes；The characteristic coefficient of the two is calculated based on the corresponding relationship of stratigraphic boundary element and contour line feature；Based on characteristic coefficient, judgment curves it is parallel；Based on parallel curves, horizontal stratum boundary line is extracted；Horizontal stratum is generated based on horizontal stratum boundary line.Compared with prior art, the present invention extraction process is converted to Computer Automatic Extraction from artificial extraction, extraction efficiency and quality are improved.

Description

A kind of horizontal stratum extraction method based on oriented parallel feature
Technical field
The present invention relates to geographical information technology application fields, and in particular to a kind of horizontal stratum based on oriented parallel feature Extraction method.
Background technique
Horizontal stratum (horizontal stratum) refer to occurrence in horizontal or close horizontal rock stratum (<5°).Horizontal rock Layer is only integral elevating and the stratum without pitching with the original level occurrence of fold after deposition diagenesis, also includes by construction It changes, but still there is nearly flat lying stratum, such as the rock stratum of Huge Recumbent Foldss both wings.
Horizontal stratum shows as stratigraphic boundary on geologic map and is roughly parallel to topographic contour, so the knowledge of horizontal stratum It Zhu Yaoshi not be realized by the parallel judgement of curve.
Since the nineties in last century, just having many people to similarity of curves both at home and abroad, this project expands research, bent One of the direction of the parallel and similarity of curves research of line.With the development of computer technology, it is similar to there are many judgements The algorithm of property.Mainly function method of identification and method of characteristic more commonly used at present.Function method of identification mainly passes through definition one A similarity definition, then carry out the calculating between distance, according to definition come the similitude of judgment curves, due to need by curve into Line function, thus it is sufficiently complex, but accuracy rate is higher.And characteristic value rule is by carrying out research comparison to characteristic parameter Carry out the similitude between judgment curves, although accuracy rate is relatively simple practical not as good as function method of identification, but existing spy Value indicative method is all using manual identified method, so that the extraction efficiency of horizontal stratum is lower.
Summary of the invention
Goal of the invention：In order to overcome the deficiencies in the prior art, a kind of level based on oriented parallel feature is provided Rock stratum extraction method, by existing manual identified method promotion arrived Computer Automatic Recognition, using Hausdorff away from From improving horizontal stratum extraction efficiency using the method for judging characteristic coefficient.
Technical solution：To achieve the above object, it is automatic to provide a kind of horizontal stratum based on oriented parallel feature by the present invention Extracting method includes the following steps：
1) it generates stratigraphic boundary data Boundary and is handled；
2) corresponding relationship of stratigraphic boundary data Boundary and contour line data Elevation are obtained；
3) primary dcreening operation is carried out based on length relation coefficient Log；
4) primary dcreening operation is carried out based on angular relationship coefficient Ang；
5) postsearch screening is carried out based on distance relation coefficient Dis；
6) horizontal stratum boundary line HorizontalLine is extracted based on parallel coefficient Par；
7) horizontal stratum figure layer HorizontalPolygonal is generated based on horizontal stratum boundary line.
Further, the specific method step of the step 1 is：
Geology face data 1-1) is read in Polygonal={ pk| k=0,1,2 ..., h-1 }, h is that prime number is wanted in geology face Amount, is converted to Linear element data to Boundary={ l for area pattern data Polygonali| i=0,1,2 ..., m-1 }, Wherein, m is stratigraphic boundary element quantity, and Boundary is stratigraphic boundary data；
Field P_ID 1-2) is added in the attribute list of Boundary data, field P_ID is for recording stratigraphic boundary element liCorresponding geology face element pkId information, generate liAnd pkCorrespondence set T1={ (li, pk) | i ∈ [0, m-1], k∈[0,h-1]}。
Further, the specific method step of the step 2 is：
Contour line data 2-1) is read in Elevation={ ej| j=0,1,2 ..., n-1 }, wherein n is contour line feature Quantity does overlay analysis to contour line data Elevation and geology face data Polygonal.Work as ejWith pkWhen intersection, according to Intersection point interrupts contour line feature ejAnd it is generated as new segmental arc fj, and assign new segmental arc ejAnd pkAttribute value.Work as ejIn pkWhen interior, Directly by contour line feature ejIt is generated as new segmental arc fjAnd assign new segmental arc ejWith pkAttribute value；
2-2) the new segmental arc f in integration step 2-1j, generate overlay analysis figure layer NewElevation, NewElevation Attribute list in include contour line feature ejID attribute and geology face element pkID attribute (as shown in the table), generate ejWith pkCorrespondence set T2={ (ej, pk)|j∈[0,n-1],k∈[0,h-1]}；
FID Shape Id_Ele Id_Pol
0 esriGemotryPolygon 0 5
1 esriGemotryPolygon 1 13
... ... ... ...
2-3) it is based on correspondence set T1And T2, according to formula (1), generate liAnd ejCorresponding relationship matrix R={ rij| I=0,1,2 ..., m-1；J=0,1,2 ..., n-1 }, wherein m is stratigraphic boundary element quantity, and n is contour line feature quantity：
2-4), it is based on corresponding relationship matrix R, generates distance relation matrix D={ di according to formula (2)j| i=0,1, 2 ..., m-1；J=0,1,2 ..., n-1 }, wherein m is stratigraphic boundary element quantity, and n is contour line feature quantity：
Wherein, max is the variable that user specifies, in rijIt is d in the case where=0ijAssignment, H (li,ej) it is according to public affairs Formula (3), the stratigraphic boundary element l of calculatingiWith corresponding contour line feature ejHausdorff distance, calculation method is as follows：
H (A, B)=max (h (A, B), h (B, A)) (3)
Wherein the ‖ ‖ in formula (4) and formula (5) is the distance between point set A and point set B normal form, the h in formula (3) (A, B) and h (B, A) be respectively maximum value in the minimum range of the point in the point-to-point collection B in point set A and point set B to A most Maximum value in small distance, H (A, B) are the maximum value in h (A, B) and h (B, A) the two, i.e. Hausdorff distance；
It 2-5) is based on distance relation matrix D, finds stratigraphic boundary element liHausdorff apart from minimum value dijIt is corresponding Contour line feature ej, generate liWith ejCorrespondence set T3={ (li, ej)|i∈[0,m-1],j∈[0,n-1]}。
Further, the specific method step of the step 3 is：
3-1) it is based on correspondence set T3, according to formula (6), calculate stratigraphic boundary element liLength DlWith it is corresponding etc. High line feature ejLength De：
Wherein, aiIt include x coordinate and y-coordinate for the point for constituting curve；
3-2) according to formula (7), computational length coefficient of relationship Log (Log ∈ [0,1]) sets length relation coefficient threshold TLog(TLog∈[0,0.6]).As Log≤TLogWhen, execute step 4；As Log > TLogWhen, execute step 3；
Further, the specific method step of the step 4 is：
4-1) it is based on correspondence set T3, according to formula (8), calculate stratigraphic boundary element liAngle [alpha]lWith it is corresponding etc. High line feature ejAngle [alpha]e：
Wherein, an-1For the terminal for constituting curve, a0It include x coordinate and y-coordinate for the starting point for constituting curve；
4-2) according to formula (9), calculate angular relationship coefficient Ang (Ang ∈ [0,1]), set angle coefficient of relationship threshold value TAng(TAng∈[0,0.5]).As Ang≤TAngWhen, execute step 5；As Ang > TAngWhen, execute step 3；
Ang=| sin (αl-αe)| (9)
Further, the specific method step of the step 5 is：
5-1) it is based on correspondence set T3, according to formula (10), calculate distance relation coefficient Dis (Dis ∈ [0,1]). Set distance coefficient of relationship threshold value TDis(TDis∈[0.2,0.6]).As Dis≤TDisWhen, execute step 6；As Dis > TDisWhen, Execute step 3；
Wherein, diIt is stratigraphic boundary element liWith corresponding contour line feature ejHausdorff distance, dmaxAnd dminIt is institute There are stratigraphic boundary and the maximum value and minimum value in the Hausdorff of corresponding contour distance.
Further, the specific method step of the step 6 is：
6-1) it is based on set of relationship T3, according to formula (11), calculate parallel degree P；
Wherein, diIt is stratigraphic boundary element liIn each point to corresponding contour line feature ejThe shortest distance, davgIt is these The average value of the shortest distance；
6-2) parallel coefficient Par (Par ∈ [0,1]) is calculated according to formula (12).Set parallel coefficient threshold TPar(TPar∈ [0.2,0.6]).As Par≤TParWhen, extract liIt is put into new horizontal stratum boundary line figure layer HorizontalLine, and holds Row step 3 is until be fully completed；As Par > TParWhen, execute step 3；
Wherein, PiIt is stratigraphic boundary element liWith corresponding contour line feature ejParallel degree, PmaxAnd PminIt is all stratum Boundary line and the maximum value and minimum value in the parallel degree of corresponding contour.
Further, the specific method step of the step 7 is：
It 7-1) is based on horizontal stratum boundary line HorizontalLine and map boundary line line, generates the horizontal stratum figure of planar Layer HorizontalPolygonal.
Design principle of the invention is：The present invention uses method of characteristic, on the basis of topographic and geologic diagram data, by obtaining The corresponding relationship of stratigraphic boundary and contour is taken, the characteristic coefficient of stratigraphic boundary and corresponding contour is calculated, to judge stratum Whether boundary line is parallel with contour, last to generate horizontal stratum according to the stratigraphic boundary for being parallel to contour, is mentioned automatically with reaching Take the purpose of horizontal stratum.
Beneficial effect：Compared with prior art, the present invention can be on the basis of topographic and geologic diagram data, by by geology Face data is converted to stratigraphic boundary data, regenerates the corresponding relationship of stratigraphic boundary and contour, calculates stratigraphic boundary and correspondence The characteristic coefficient of contour judges whether stratigraphic boundary is parallel with contour, last according to the stratigraphic boundary for being parallel to contour Horizontal stratum boundary line is generated, horizontal stratum is generated based on horizontal stratum boundary line.It realizes existing manual identified method Computer Automatic Recognition is promoted, so that the efficiency for extracting horizontal stratum greatly improved, and the extraction matter of horizontal stratum Amount is also obviously improved.
Detailed description of the invention
Fig. 1 is linghe terrain and geologic map；
Fig. 2 is the flow diagram of the method for the present invention；
Fig. 3 is the distribution situation of linghe area stratigraphic boundary and contour；
Fig. 4 is the horizontal stratum boundary line schematic diagram extracted；
Fig. 5 is the horizontal stratum schematic diagram generated.
Specific embodiment
In the following with reference to the drawings and specific embodiments, the present invention is furture elucidated.
As shown in Figure 1, the experimental data of this example is using 1:20000 linghe topographic and geologic diagram data goes out in area The figure layer of dew is roughly divided into 4 ages, is well-developed land in gully.
As shown in Fig. 2, the present embodiment provides a kind of horizontal stratum extraction methods based on oriented parallel feature, including Following steps：
1) it generates stratigraphic boundary data Boundary and is handled；
Geology face data 1-1) is read in Polygonal={ pk| k=0,1,2 ..., h-1 }, in the present embodiment, h=21. Area pattern data Polygonal is converted into Linear element data to Boundary={ li| i=0,1,2 ..., m-1 }, this reality It applies in example, m=42.Boundary is stratigraphic boundary data, specific as shown in Figure 3；
Field P_ID 1-2) is added in the attribute list of Boundary data, field P_ID is for recording stratigraphic boundary element liCorresponding geology face element pkID attribute, generate liWith pkCorrespondence set T1={ (l0, p13), (l1, p20), (l2, p20) ..., (l41, p6)}；
2) corresponding relationship of stratigraphic boundary data Boundary and contour line data Elevation are obtained：
Contour line data 2-1) is read in Elevation={ ej| j=0,1,2 ..., n-1 }, in the present embodiment, n=18. Overlay analysis is done to contour line data Elevation and geology face data Polygonal.Work as ejWith pkWhen intersection, beaten according to intersection point Break contour line feature ejAnd it is generated as new segmental arc fj, and assign new segmental arc ejAnd pkAttribute value.Work as ejIn pkWhen interior, directly will Contour line feature ejIt is generated as new segmental arc fjAnd assign new segmental arc ejWith pkAttribute value；
2-2) the new segmental arc f in integration step 2-1j, generate figure layer NewElevation, the attribute list of NewElevation In include contour line feature ejID attribute and geology face element pkID attribute (as shown in the table), generate ejWith pkCorrespondence Set of relationship T2={ (e0, p13), (e0, p17) ..., (e17, p0)}；
FID Shape Id_Ele Id_Pol
0 esriGemotryPolygon 0 13
1 esriGemotryPolygon 0 17
... ... ... ...
52 esriGemotryPolygon 17 0
2-3) it is based on correspondence set T1And T2, according to formula (1), generate liAnd ejCorresponding relationship matrix R={ rij| I=0,1,2 ..., m-1；J=0,1,2 ..., n-1 }, in the present embodiment, m=42, n=18；
It 2-4) is based on corresponding relationship matrix R, generates distance relation matrix D={ d according to formula (2)ij| i=0,1,2 ..., m-1；J=0,1,2 ..., n-1 }, in the present embodiment, m=42, n=18：
Wherein, max is the variable that user specifies, in rijIt is d in the case where=0ijAssignment, H (li,ej) it is according to public affairs Formula (3), the stratigraphic boundary element l of calculatingiWith corresponding contour line feature ejHausdorff distance, calculation method is as follows：
H (A, B)=max (h (A, B), h (B, A)) (3)
Wherein the ‖ ‖ in formula (4) and formula (5) is the distance between point set A and point set B normal form, the h in formula (3) (A, B) and h (B, A) be respectively maximum value in the minimum range of the point in the point-to-point collection B in point set A and point set B to A most Maximum value in small distance, H (A, B) are the maximum value in h (A, B) and h (B, A) the two, i.e. Hausdorff distance；
It 2-5) is based on distance relation matrix D, finds stratigraphic boundary element liHausdorff apart from minimum value dijIt is corresponding Contour line feature ej, generate liWith ejCorrespondence set T3={ (l0, e1), (l1, e9) ..., (l41, e13)}；
3) primary dcreening operation is carried out based on length relation coefficient Log：
3-1) it is based on correspondence set T3, according to formula (6), calculate stratigraphic boundary element liLength DlWith it is corresponding etc. High line feature eiLength De；
Wherein, aiIt include x coordinate and y-coordinate for the point for constituting curve.
3-2) according to formula (7), computational length coefficient of relationship Log (Log ∈ [0,1]).Set length relation coefficient threshold TLog=0.6, as Log≤TLogWhen；Execute step 4；As Log > TLogWhen, execute step 3；
4) primary dcreening operation is carried out based on angular relationship coefficient Ang：
4-1) it is based on correspondence set T3, according to formula (8), calculate stratigraphic boundary element liAngle [alpha]lWith it is corresponding etc. High line feature ejAngle [alpha]e；
Wherein, an-1For the terminal for constituting curve, a0It include x coordinate and y-coordinate for the starting point for constituting curve；
4-2) according to formula (9), calculate angular relationship coefficient Ang (Ang ∈ [0,1]).Set angle coefficient of relationship threshold value TAng=0.2, as Ang≤TAngWhen, execute step 5；As Ang > TAngWhen, execute step 3；
Ang=| sin (αl-αe)| (9)
5) postsearch screening is carried out based on distance relation coefficient Dis：
5-1) it is based on correspondence set T3, according to formula (10), calculate distance relation coefficient Dis (Dis ∈ [0,1]). Set distance coefficient of relationship threshold value TDis=0.2, as Dis≤TDisWhen, execute step 6；As Dis > TDisWhen, execute step 3；
Wherein, diIt is stratigraphic boundary element liWith corresponding contour line feature ejHausdorff distance, dmaxAnd dminIt is institute There are stratigraphic boundary and the maximum value and minimum value in the Hausdorff of corresponding contour distance；
6) horizontal stratum boundary line HorizontalLine is extracted based on parallel coefficient Par：
6-1) it is based on set of relationship T3, according to formula (11), calculate parallel degree P；
Wherein, diIt is stratigraphic boundary element liIn each point to corresponding contour line feature ejThe shortest distance, davgIt is these The average value of the shortest distance；
6-2) according to formula (12), parallel coefficient Par (Par ∈ [0,1]) is calculated.Set parallel coefficient threshold Tpar= 0.2, as Par≤TparWhen, extract liIt is put into new horizontal stratum boundary line figure layer HorizontalLine, and executes step 3 Until being fully completed；As Par > TpWhen ar, step 3 is executed；
Wherein, PiIt is stratigraphic boundary element liWith corresponding contour line feature ejParallel degree, PmaxAnd PminIt is all stratum Boundary line and the maximum value and minimum value in the parallel degree of corresponding contour.
7) horizontal stratum HorizontalPolygonal is generated based on horizontal stratum boundary line；
7-1) based on new stratigraphic boundary HorizontalLine (specific as shown in Figure 4) and map boundary line line, face is generated The horizontal stratum HorizontalPolygonal (specific as shown in Figure 5) of shape.
Interpretation of result：In present case, it is related to 21 geology faces, 42 stratigraphic boundaries, correct result of extracting should be Fig. 1 Shown in upper Cretaceous series sandstone and Lower Cretaceous Series conglomerate region, totally 8 geology faces, 7 stratigraphic boundaries.Actual extracting result is 7 geology faces, 6 stratigraphic boundaries.Compared to correct result, it is extracted 1 piece of horizontal stratum region less, is to be extracted shown in Fig. 5 There is no the case where error extraction or extra extraction in the blank parts in horizontal stratum region.Source of error is mainly lost The stratigraphic boundary length factor for leaking horizontal stratum is excessive, is error caused by special circumstances.On the whole, experimental result accuracy rate It is higher, it can satisfy practical application request.

Claims (8)

1. a kind of horizontal stratum extraction method based on oriented parallel feature, it is characterised in that：Include the following steps：
1) it generates stratigraphic boundary data Boundary and is handled；
2) corresponding relationship of stratigraphic boundary data Boundary and contour line data Elevation are obtained；
3) primary dcreening operation is carried out based on length relation coefficient Log；
4) primary dcreening operation is carried out based on angular relationship coefficient Ang；
5) postsearch screening is carried out based on distance relation coefficient Dis；
6) horizontal stratum boundary line HorizontalLine is extracted based on parallel coefficient Par；
7) horizontal stratum figure layer HorizontalPolygonal is generated based on horizontal stratum boundary line.
2. a kind of horizontal stratum extraction method based on oriented parallel feature according to claim 1, feature exist In：The specific method step of the step 1 is：
Geology face data 1-1) is read in Polygonal={ pk| k=0,1,2 ..., h-1 }, h is geology face element quantity, by face Shape factor data Polygonal is converted to Linear element data to Boundary={ li| i=0,1,2 ..., m-1 }, wherein m is Stratigraphic boundary element quantity, Boundary are stratigraphic boundary data；
Field P_ID 1-2) is added in the attribute list of Boundary data, field P_ID is for recording stratigraphic boundary element liInstitute Corresponding geology face element pkId information, generate liAnd pkCorrespondence set T1={ (li, pk)|i∈[0,m-1],k∈ [0,h-1]}。
3. a kind of horizontal stratum extraction method based on oriented parallel feature according to claim 1 or 2, feature It is：The specific method step of the step 2 is：
Contour line data 2-1) is read in Elevation={ ej| j=0,1,2 ..., n-1 }, wherein n is contour line feature quantity, Overlay analysis is done to contour line data Elevation and geology face data Polygonal.Work as ejWith pkWhen intersection, beaten according to intersection point Break contour line feature ejAnd it is generated as new segmental arc fj, and assign new segmental arc ejAnd pkAttribute value.Work as ejIn pkWhen interior, directly will Contour line feature ejIt is generated as new segmental arc fjAnd assign new segmental arc ejWith pkAttribute value；
2-2) the new segmental arc f in integration step 2-1j, generate the category of overlay analysis figure layer NewElevation, NewElevation Property table in include contour line feature ejID attribute and geology face element pkID attribute, generate ejWith pkCorrespondence set T2 ={ (ej, pk)|j∈[0,n-1],k∈[0,h-1]}；
2-3) it is based on correspondence set T1And T2, according to formula (1), generate liAnd ejCorresponding relationship matrix R={ rij| i= 0,1,2 ..., m-1；J=0,1,2 ..., n-1 }, wherein m is stratigraphic boundary element quantity, and n is contour line feature quantity：
2-4), it is based on corresponding relationship matrix R, generates distance relation matrix D={ d according to formula (2)ij| i=0,1,2 ..., m-1； J=0,1,2 ..., n-1 }, wherein m is stratigraphic boundary element quantity, and n is contour line feature quantity：
Wherein, max is the variable that user specifies, in rijIt is d in the case where=0ijAssignment, H (li,ej) it is according to formula (3), the stratigraphic boundary element l of calculatingiWith corresponding contour line feature ejHausdorff distance, calculation method is as follows：
H (A, B)=max (h (A, B), h (B, A)) (3)
Wherein the ‖ ‖ in formula (4) and formula (5) is the distance between point set A and point set B normal form, h in formula (3) (A, B) and h (B, A) be respectively maximum value and point set B to A in the minimum range of point in point-to-point collection B in point set A minimum Maximum value in distance, H (A, B) are the maximum value in h (A, B) and h (B, A) the two, i.e. Hausdorff distance；
It 2-5) is based on distance relation matrix D, finds stratigraphic boundary element liHausdorff apart from minimum value dijCorresponding etc. High line feature ej, generate liWith ejCorrespondence set T3={ (li, ej)|i∈[0,m-1[,j∈[0,n-1]}。
4. a kind of horizontal stratum extraction method based on oriented parallel feature according to claim 3, feature exist In：The specific method step of the step 3 is：
3-1) it is based on correspondence set T3, according to formula (6), calculate stratigraphic boundary element liLength DlWith corresponding contour Element ejLength De：
Wherein, aiIt include x coordinate and y-coordinate for the point for constituting curve；
3-2) according to formula (7), computational length coefficient of relationship Log (Log ∈ [0,1]) sets length relation coefficient threshold TLog (TLog∈[0,0.6]).As Log≤TLogWhen, execute step 4；As Log > TLogWhen, execute step 3；
5. a kind of horizontal stratum extraction method based on oriented parallel feature according to claim 3, feature exist In：The specific method step of the step 4 is：
4-1) it is based on correspondence set T3, according to formula (8), calculate stratigraphic boundary element liAngle [alpha]lWith corresponding contour Element ejAngle [alpha]e：
Wherein, an-1For the terminal for constituting curve, a0It include x coordinate and y-coordinate for the starting point for constituting curve；
4-2) according to formula (9), calculate angular relationship coefficient Ang (Ang ∈ [0,1]), set angle coefficient of relationship threshold value TAng (TAng∈[0,0.5]).As Ang≤TAngWhen, execute step 5；As Ang > TAngWhen, execute step 3；
Ang=| sin (αl-αe)|
(9)
6. a kind of horizontal stratum extraction method based on oriented parallel feature according to claim 3, feature exist In：The specific method step of the step 5 is：
5-1) it is based on correspondence set T3, according to formula (10), calculate distance relation coefficient Dis (Dis ∈ [0,1]).Setting away from From coefficient of relationship threshold value TDis(TDis∈[0.2,0.6]).As Dis≤TDisWhen, execute step 6；As Dis > TDisWhen, execute step Rapid 3；
Wherein, diIt is stratigraphic boundary element liWith corresponding contour line feature ejHausdorff distance, dmaxAnd dminIt is allly Maximum value and minimum value of the stratum boundary line with the Hausdorff of corresponding contour in.
7. a kind of horizontal stratum extraction method based on oriented parallel feature according to claim 3, feature exist In：The specific method step of the step 6 is：
6-1) it is based on set of relationship T3, according to formula (11), calculate parallel degree P；
Wherein, diIt is stratigraphic boundary element liIn each point to corresponding contour line feature ejThe shortest distance, davgIt is that these are most short The average value of distance；
6-2) parallel coefficient Par (Par ∈ [0,1]) is calculated according to formula (12).Set parallel coefficient threshold TPar(TPar∈ [0.2,0.6]).As Par≤TParWhen, extract liIt is put into new horizontal stratum boundary line figure layer HorizontalLine, and holds Row step 3 is until be fully completed；As Par > TParWhen, execute step 3；
Wherein, PiIt is stratigraphic boundary element liWith corresponding contour line feature ejParallel degree, PmaxAnd PminIt is all stratigraphic boundaries With the maximum value and minimum value in the parallel degree of corresponding contour.
8. a kind of horizontal stratum extraction method based on oriented parallel feature according to claim 1 or claim 7, feature It is：The specific method step of the step 7 is：
It 7-1) is based on horizontal stratum boundary line HorizontalLine and map boundary line line, generates the horizontal stratum figure layer of planar HorizontalPolygonal。
CN201810521612.3A 2018-05-28 2018-05-28 Automatic horizontal rock stratum extraction method based on curve parallel characteristics Active CN108875755B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810521612.3A CN108875755B (en) 2018-05-28 2018-05-28 Automatic horizontal rock stratum extraction method based on curve parallel characteristics

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810521612.3A CN108875755B (en) 2018-05-28 2018-05-28 Automatic horizontal rock stratum extraction method based on curve parallel characteristics

Publications (2)

Publication Number Publication Date
CN108875755A true CN108875755A (en) 2018-11-23
CN108875755B CN108875755B (en) 2021-11-12

Family

ID=64335218

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810521612.3A Active CN108875755B (en) 2018-05-28 2018-05-28 Automatic horizontal rock stratum extraction method based on curve parallel characteristics

Country Status (1)

Country Link
CN (1) CN108875755B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113539051A (en) * 2021-06-24 2021-10-22 南京师范大学 Geological map-based stratum boundary point-by-point rock stratum attitude obtaining method and device

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6298155B1 (en) * 1997-07-31 2001-10-02 Sony Corporation Apparatus and method for generating time-series data, apparatus and method for editing curves and recording medium having program recorded thereon
CN106023197A (en) * 2016-05-18 2016-10-12 南京师范大学 Automated identification and extraction method of vertical stratum
CN106934357A (en) * 2017-02-28 2017-07-07 南京师范大学 A kind of automatic identifying method of parallel fault

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6298155B1 (en) * 1997-07-31 2001-10-02 Sony Corporation Apparatus and method for generating time-series data, apparatus and method for editing curves and recording medium having program recorded thereon
CN106023197A (en) * 2016-05-18 2016-10-12 南京师范大学 Automated identification and extraction method of vertical stratum
CN106934357A (en) * 2017-02-28 2017-07-07 南京师范大学 A kind of automatic identifying method of parallel fault

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113539051A (en) * 2021-06-24 2021-10-22 南京师范大学 Geological map-based stratum boundary point-by-point rock stratum attitude obtaining method and device

Also Published As

Publication number Publication date
CN108875755B (en) 2021-11-12

Similar Documents

Publication Publication Date Title
CN109800863B (en) Logging phase identification method based on fuzzy theory and neural network
CN105139379B (en) Based on the progressive extracting method of classified and layered airborne Lidar points cloud building top surface
CN111583407B (en) Efficient three-dimensional geological modeling intelligent processing method based on paper drilling
CN107977664B (en) Road vanishing point detection method based on single image
CN113344019A (en) K-means algorithm for improving decision value selection initial clustering center
CN101916373A (en) Road semiautomatic extraction method based on wavelet detection and ridge line tracking
CN105740521B (en) Small grid elimination method and device during reservoir numerical simulation system solution
CN104881852A (en) Image segmentation method based on immune clone algorithm and fuzzy kernel-clustering algorithm
CN112945196A (en) Strip mine step line extraction and slope monitoring method based on point cloud data
CN109102538A (en) Method, the synoptic analysis method of piston ring land characteristic point and land features line are extracted using isopleth data
CN112487537A (en) Building surface multistage optimization extraction method based on full convolution neural network
CN112241676A (en) Method for automatically identifying terrain sundries
CN115546116A (en) Method and system for extracting and calculating spacing of discontinuous surface of fully-covered rock mass
CN108875755A (en) A kind of horizontal stratum extraction method based on oriented parallel feature
CN105740878B (en) The Symbol Recognition and device of GIS-Geographic Information System
CN113689453A (en) Method, device and equipment for automatically identifying well logging image cracks and storage medium
CN107066997B (en) A kind of electrical component price quoting method based on image recognition
CN103208003B (en) Geometric graphic feature point-based method for establishing shape descriptor
CN104123382B (en) A kind of image set abstraction generating method under Social Media
CN104102718A (en) Method for unbalanced data classification by unbalanced data set oriented cluster boundary sampling and dynamic pruning
CN109543846B (en) MVO (mechanical vapor deposition) improvement based DBSCAN (direct species analysis controller area network) mine water inrush spectrum identification method
CN111340000A (en) Method and system for extracting and optimizing PDF document table
CN112269951B (en) Linear shape space retrieval method for vector line data
CN109033746A (en) A kind of protein complex recognizing method based on knot vector
Liu et al. The classification of airborne LiDAR building point clouds based on multi‐scale and multi‐level cloth simulation

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