CN109325987A - A kind of tomography polygon rendering figure auto-correction method and device - Google Patents
A kind of tomography polygon rendering figure auto-correction method and device Download PDFInfo
- Publication number
- CN109325987A CN109325987A CN201811069856.9A CN201811069856A CN109325987A CN 109325987 A CN109325987 A CN 109325987A CN 201811069856 A CN201811069856 A CN 201811069856A CN 109325987 A CN109325987 A CN 109325987A
- Authority
- CN
- China
- Prior art keywords
- coordinate
- point
- spacing
- sampling point
- symmetrical
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/203—Drawing of straight lines or curves
Abstract
The invention discloses a kind of tomography polygon rendering figure auto-correction method and devices, comprising: obtains the first intersecting point coordinate between tomography polygon and main profile, and the second intersecting point coordinate between cross-track；The corresponding plane coordinates of two endpoints is determined respectively；Obtain the normal of preset quantity corresponding with the first straight line after two endpoints connections；Determine the third intersection point between normal that the second straight line between adjacent first intersection point and the second intersection point intersects with second straight line；When determining that third intersection point is between adjacent the first intersection point and the second intersection point, third intersection point is saved；Count the coordinate of all third intersection points, and be divided on symmetrical make an inventory coordinate and under make an inventory coordinate.Determine make an inventory on symmetrical coordinate and under the spacing made an inventory between coordinate；Using maximum spacing as benchmark spacing, after reducing benchmark spacing according to preset rules, symmetrical sampling point is adjusted to ordinate；It is smoothed after being encrypted using method for resampling, obtaining final tomography polygon can chart.
Description
Technical field
The present invention relates to technical field of data processing, and in particular to a kind of tomography polygon rendering figure auto-correction method and
Device.
Background technique
It is looked in miner's work in the oil prospecting of current geology industry, the establishment of structural map is an important action, mainly
It is the kinematics informations such as reflection interval, same to phasic property, the velocity of wave become using seismic wave, studies the distribution and fluctuating of bed boundary
Form, mature fault situation, and be transformed into the depth of bed boundary the hourage in seismic events section, draw ground texture
Figure is made, provides basic document to find structural deposit.Wherein, the plane drafting of tomography polygons (polygon) is that its is heavy
The component part wanted draws the mapping precision that precision directly affects structural map, and the precision of structural map is directly affected to trap
It is extremely important finely accurately therefore to draw tomography polygons with accurately identifying for oil reservoir.Current method for drafting mainly has
Following two: one is completing by manually, manual drawing is carried out according to tomography and the aperture position of layer position, this mode exists
It is more universal in practical application；Second is attached automatically using fracture of the layer position at each tomography.First way
It is drawn using artificial, it is time-consuming and laborious although accuracy is higher, especially area is big, the exploration acreage more than tomography, if people
Work editor's tomography polygons, every figure require a large amount of mechanicalness, repetitive operation, and time-consuming effort greatly affected
Job schedule；Although improving the speed of drafting for second by machine automatic identification, there are the following problems:
1, Polygons heave is changed greatly along tomography extending direction, and it is strong not meet general tomography intermediate stress effect
Heave is big, the small geologic rule of the weak turn-off of both sides stress；Mainly there are two aspect reasons: first is that human factor causes,
During practice interpretation, since layer position cannot snap on tomography completely, there is a fixed gap；Second is that current seism processing
Technology can't carry out accurate imaging to subsurface fault, cause in interpretation process that there are unbalanced factors for heave；
2, Polygons is not round and smooth, and catastrophe point is more；
3, Polygons wherein do not close by one end, and the trend of closed one end not instead of pinching, rectangle.
Summary of the invention
The embodiment of the present invention is designed to provide a kind of tomography polygon rendering figure auto-correction method and device, to
The low problem of the precision that geologic structure diagram is drawn automatically to existing machine, to promote ground while promoting working efficiency
The precision of matter structural map.
To achieve the above object, the technical solution of the embodiment of the present invention provides a kind of automatic school of tomography polygon rendering figure
Correction method, this method comprises:
It is more that each main profile and tomography in a plurality of main profile pre-seted are obtained from original tomography polygon rendering figure
The first intersecting point coordinate between the shape of side；And in the cross-track pre-seted between each cross-track and tomography polygon
Second intersecting point coordinate, original tomography polygon can chart as the tomography polygon rendering figure drawn by machine automatic identification；
The corresponding plane coordinates of two endpoints in original tomography polygon rendering figure, two of them endpoint pair are determined respectively
The plane coordinates answered is the first intersecting point coordinate and/or the second intersecting point coordinate；
According to the corresponding plane coordinates of two endpoints, present count corresponding with the first straight line after two endpoints connections is obtained
The normal of amount；
It determines between normal that the second straight line between adjacent first intersection point and the second intersection point intersects with second straight line
Third intersection point；
When determining that third intersection point is between adjacent the first intersection point and the second intersection point, third intersection point is saved, is otherwise lost
Abandon third intersection point；
All third intersection points are counted, and obtain the coordinate of third intersection point；The coordinate of third intersection point is divided into symmetrically
On make an inventory coordinate and under make an inventory coordinate；
Determine respectively make an inventory on symmetrical coordinate and under the spacing made an inventory between coordinate, and obtain maximum spacing；
By maximum spacing it is corresponding on make an inventory coordinate and under make an inventory symmetrical sampling point pair centered on coordinate, using maximum spacing as base
Quasi- spacing, after reducing benchmark spacing according to preset rules, as between the symmetrical sampling point pair in addition to central symmetry sampling point pair
Spacing adjusts the ordinate of the symmetrical sampling point pair in addition to central symmetry sampling point pair, obtains the first tomography polygon rendering figure,
In, it is gradually reduced far from the spacing of central symmetry sampling point pair, until be zero, symmetrical sampling point is to being to make an inventory on symmetrical and lower wall
The sampling point pair that point collectively constitutes；
Using the sampling point in method for resampling the two the first tomography polygon rendering figure of encryption；
Encrypted sampling point is smoothed, obtaining final tomography polygon can chart.
The embodiment of the present invention has the advantages that the first intersecting point coordinate obtained between main profile and tomography polygon, and
After the second intersecting point coordinate between cross-track and tomography polygon, it is first determined two extreme coordinates on tomography shape changeable,
Then first straight line is connected into, and obtains a plurality of normal corresponding with first straight line.And the friendship of a plurality of normal and tomography polygon
Point trend more levels off to the trend of tomography, therefore draws tomography polygon rendering using the intersection point of a plurality of normal and tomography polygon
Figure will be more accurate.Before anything is drawn, it is also necessary to using the spaced relationship between symmetrical sampling point pair, the different symmetrical sampling points of adjustment
Ordinate data between, change the horizontal space of tomography polygon more uniformly, and both ends can be realized nature and annihilate,
Meet geologic rule.Finally, using method for resampling encrypt the second tomography polygon can chart in sampling point and smoothly located
Reason, makes tomography that can finally move towards smooth nature.
To achieve the above object, the technical solution of the embodiment of the present invention provides a kind of automatic school of tomography polygon rendering figure
Equipment, the device include:
Acquiring unit, for obtaining each master in a plurality of main profile pre-seted from original tomography polygon rendering figure
The first intersecting point coordinate between survey line and tomography polygon；And each cross-track and tomography in the cross-track pre-seted
The second intersecting point coordinate between polygon, the tomography that original tomography polygon can chart to be drawn by machine automatic identification are polygon
Shape draws figure；
The corresponding plane coordinates of two endpoints in original tomography polygon rendering figure, two of them endpoint pair are determined respectively
The plane coordinates answered is the first intersecting point coordinate and/or the second intersecting point coordinate；
Processing unit, for according to the corresponding plane coordinates of two endpoints, first after acquisition is connect with two endpoints to be straight
The normal of the corresponding preset quantity of line；
It determines between normal that the second straight line between adjacent first intersection point and the second intersection point intersects with second straight line
Third intersection point；
When determining that third intersection point is between adjacent the first intersection point and the second intersection point, third intersection point is saved, is otherwise lost
Abandon third intersection point；
Statistic unit for counting all third intersection points, and obtains the coordinate of third intersection point；By the coordinate of third intersection point
Be divided on symmetrical make an inventory coordinate and under make an inventory coordinate；
Processing unit is also used to, determine respectively make an inventory on symmetrical coordinate and under the spacing made an inventory between coordinate, and obtain
Maximum spacing；
By maximum spacing it is corresponding on make an inventory coordinate and under make an inventory symmetrical sampling point pair centered on coordinate, using maximum spacing as base
Quasi- spacing, after reducing benchmark spacing according to preset rules, as between the symmetrical sampling point pair in addition to central symmetry sampling point pair
Spacing adjusts the ordinate of the symmetrical sampling point pair in addition to central symmetry sampling point pair, obtains the first tomography polygon rendering figure,
In, it is gradually reduced far from the spacing of central symmetry sampling point pair, until be zero, symmetrical sampling point is to being to make an inventory on symmetrical and lower wall
The sampling point pair that point collectively constitutes；
Encryption unit, for encrypting the sampling point in the first tomography polygon rendering figure using method for resampling；
Processing unit is also used to, and is smoothed to encrypted sampling point, and obtaining final tomography polygon can chart.
The embodiment of the present invention has the advantages that the first intersecting point coordinate obtained between main profile and tomography polygon, and
After the second intersecting point coordinate between cross-track and tomography polygon, it is first determined two extreme coordinates on tomography shape changeable,
Then first straight line is connected into, and obtains a plurality of normal corresponding with first straight line.And the friendship of a plurality of normal and tomography polygon
Point trend more levels off to the trend of tomography, therefore draws tomography polygon rendering using the intersection point of a plurality of normal and tomography polygon
Figure will be more accurate.Before anything is drawn, it is also necessary to using the spaced relationship between symmetrical sampling point pair, the different symmetrical sampling points of adjustment
Ordinate data between, change the horizontal space of tomography polygon more uniformly, and both ends can be realized nature and annihilate,
Meet geologic rule.Finally, using method for resampling encrypt the second tomography polygon can chart in sampling point and smoothly located
Reason, makes tomography that can finally move towards smooth nature.
Detailed description of the invention
Fig. 1 is a kind of tomography polygon rendering figure auto-correction method flow diagram provided in an embodiment of the present invention；
Fig. 2 is the original tomography polygon rendering figure provided by the invention comprising main profile and contact measuring point；
Fig. 3 is the drafting first straight line provided by the invention on the basis of Fig. 2 and the schematic diagram after a plurality of normal；
Fig. 4 is the schematic diagram provided by the invention on the basis of Fig. 3 after drafting second straight line；
Fig. 5 is the contrast schematic diagram before and after the heave progress equilibrium treatment provided by the invention to tomography polygon；
Fig. 6 is the front and back contrast schematic diagram provided by the invention being smoothed to tomography；
Fig. 7 is a kind of tomography polygon rendering figure automatic correction device structural schematic diagram provided in an embodiment of the present invention.
Specific embodiment
The following examples are used to illustrate the present invention, but are not intended to limit the scope of the present invention..
Embodiment 1
The embodiment of the present invention 1 provides a kind of tomography polygon rendering figure auto-correction method flow diagram.This method
Mainly it is corrected to having already passed through the original tomography polygon drawn after traditional machine automatic identification and can chart.So as to
Heave is not balanced enough in the tomography polygon rendering figure for overcoming conventional method to draw, and does not meet geologic rule, polygon
It is sufficiently oily, the problems such as catastrophe point is more and polygon both ends are not the trend annihilated.
This method flow diagram is specific as shown in Figure 1, this method may include:
Step 110, each main profile in a plurality of main profile pre-seted is obtained from original tomography polygon rendering figure
The first intersecting point coordinate between tomography polygon；And each cross-track and tomography are polygon in the cross-track pre-seted
The second intersecting point coordinate between shape.
Specifically, since different types of tomography is not quite similar upwards walking, some fault strike azimuths and main profile
Angle is small, and some is then bigger, in general lesser with main profile angle, and main profile tomography fracture obtained is relatively
More, cross-track fracture is then relatively fewer, and vice versa, therefore will be from main profile and cross-track two when calculating heave
A direction is respectively calculated, ensuring the planar distribution form of tomography polygon polygons.
Therefore, the original tomography polygon that machine automatic identification first is drawn can chart, and it includes pre- for obtaining coordinate data
The first intersecting point coordinate in a plurality of registration lines being arranged between each main profile and tomography polygon, and the connection pre-seted
The second intersecting point coordinate in network survey line between each cross-track and tomography polygon.It is specific as shown in Figure 2.It is shown in Fig. 2
Be during original tomography polygon can chart only have main profile and cross-track schematic diagram.It is laterally main profile, it is longitudinal for connection
Network survey line.
Step 120, the corresponding plane coordinates of two endpoints in original tomography polygon rendering figure is determined respectively.
Shown in Fig. 2, the corresponding plane coordinates of two of them endpoint is the first intersecting point coordinate and/or the second intersecting point coordinate.
Step 130, according to the corresponding plane coordinates of two endpoints, the first straight line obtained after connecting with two endpoints is corresponding
Preset quantity normal.
Specifically, as shown in figure 3, two o'clock determines straight line.And in present specification, then it is more by original tomography
The corresponding plane coordinates of two endpoints in the shape drafting figure of side determines first straight line y=k1x+b.And specific determination process is then
It is that can be obtained by data algorithm, it is known that the coordinate (x1, y1) of two o'clock and (x2, y2) can obtain the parameter in first straight line
K and b, and then can determine first straight line.In turn, it is known that the coordinate data of the fooled any point of this straight line.
Then, according to the corresponding plane coordinates of two endpoints, the first straight line obtained after connecting with two endpoints is corresponding
The normal of preset quantity.That is, points presetting in such as first straight line include 10 points, the transverse and longitudinal coordinate of each point is ok
Know, and normal corresponding with first straight line is the straight line vertical with first straight line, and, intersection point is a point in 10 points,
On the basis of each point in the default points in first straight line, the normal of preset quantity corresponding with first straight line is obtained.And
The coordinate of each of preset point can be known by mathematical algorithm, then can be according to the expression for determining each normal
Formula yMethod i=kMethod ix+bMethod i, i is i-th normal.
Optionally, in selected normal, the spacing that can be set between two adjacent normals is identical.It is demonstrated in Figure 3
Be that spacing is identical between normal.
Step 140, the method that the second straight line between adjacent first intersection point and the second intersection point intersects with second straight line is determined
Third intersection point between line.
Step 150, it when determining that third intersection point is between adjacent the first intersection point and the second intersection point, saves third and hands over
Otherwise point abandons third intersection point.
Specifically, then being needed to determine the intersecting point coordinate in normal and original tomography polygon rendering figure between curve
Determine the third intersection point between normal that the second straight line between adjacent first intersection point and the second intersection point intersects with second straight line.
Same reason, two o'clock determine straight line.Specifically as shown in figure 4, in Fig. 4, with one group of first adjacent intersection point
For second straight line between coordinate and the second intersecting point coordinate, intersect with wherein one or more normal, and after intersecting, it is inevitable
There are an intersection points.At this point it is possible to determine intersection point whether on the curve that original tomography polygon rendering figure is constituted.If
Then retain the intersection point, otherwise rejects the intersection point.And specific judgment method is then: judging whether third intersection point is located at first
Between intersection point and the second intersection point.And whether the foundation judged is then determining third intersection point to the distance between first straight line less than the
One intersection point arrives the distance between first straight line.If being less than, it is determined that third intersection point is located at original tomography polygon rendering figure institute structure
At curve on, otherwise, be not on the curve that original tomography polygon rendering figure is constituted.And determining third intersection point position
When between adjacent the first intersection point and the second intersection point, third intersection point is saved, otherwise abandons third intersection point.
Step 160, all third intersection points are counted, and obtain the coordinate of third intersection point.
Second straight line and the expression formula of normal are known, then the coordinate of its third intersection point can necessarily be calculated by mathematics
Method is found out.Therefore, the coordinate of all third intersection point and third intersection point can be counted.
Step 170, the coordinate of third intersection point is divided on symmetrical make an inventory coordinate and under make an inventory coordinate.
Specifically, the coordinate data in original tomography polygon rendering figure is number of coordinates of making an inventory in first output in output
According to then exporting down the coordinate data made an inventory again.Make an inventory on and and under make an inventory between there are symmetric relations, therefore be easily discriminated
For output, first half is coordinate data of above making an inventory, and latter half is lower wall point coordinate data.And coordinate of making an inventory on first
Data and first lower wall point coordinate data are corresponding.And so on, so that it may easily know that those points are made an inventory on being, which
A little points are made an inventory under being, and then obtain corresponding coordinate data.
Step 180, determine respectively make an inventory on symmetrical coordinate and under the spacing made an inventory between coordinate, and obtain between maximum
Away from.
Step 190, by maximum spacing it is corresponding on make an inventory coordinate and under make an inventory symmetrical sampling point pair centered on coordinate, with maximum
Spacing is benchmark spacing, after reducing benchmark spacing according to preset rules, as the symmetrical sampling point in addition to central symmetry sampling point pair
Spacing between adjusts the ordinate of the symmetrical sampling point pair in addition to central symmetry sampling point pair, obtains the first tomography polygon
Draw figure.
Specifically, based on many detailed observations (including to appear, the observation of mineral deposit and seismic reflection data etc.),
The turn-off idealized model of " independent " single tomography is, tomography is to connect paracentral position turn-off big, to tomography edge (end
Line) all directions reduce, in end, line position turn-off is zero.Based on this feature, it can determine that original tomography is polygon first
Shape draws the spacing in figure between every a pair of symmetrical sampling point pair, then obtains maximum spacing.And it is spacing is maximum a pair of symmetrical
Center of the sampling point to the position at place as tomography, this is to symmetrical sampling point pair it can be appreciated that central symmetry sampling point pair.
For sake of convenience, it will make an inventory in central symmetry sampling point centering and be customized for the first sampling point, the lower wall of central symmetry sampling point centering
Point is customized for the second sampling point.With the spacing between central symmetry sampling point pair for basic spacing.Symmetrical sampling point pair mentioned here,
Really be exactly make an inventory on symmetrical and under make an inventory the sampling point pair collectively constituted.
Because with the spacing between central symmetry sampling point pair to be maximum in tomography, between the symmetrical sampling point pair of the right and left
Spacing gradually decrease.Therefore, after basic spacing being reduced according to preset rules, as other than central symmetry sampling point pair
Symmetrical sampling point pair between spacing, adjust the coordinate of symmetrical sampling point pair in addition to central symmetry sampling point pair.To obtain the
One tomography polygon rendering figure.Wherein, the spacing far from central symmetry sampling point pair gradually decreases, until being zero.
Specific preset rules are as follows: using maximum spacing as benchmark spacing, reduce n multiplied by the numerical value of i after, centered on symmetrical sample
For point to the spacing of the n-th pair of symmetrical sampling point pair in the left and right sides, n is more than or equal to 1 and to be less than or equal to all symmetrical sampling points
To the positive integer of the half of total number, the value of n is incremented by successively, and being incremented by spacing is 1；I is between maximum spacing and n
Ratio.
But in order to guarantee that the longitudinal center of every a pair of symmetrical sampling point pair does not change, thus it is right adjusting every a pair
When claiming the ordinate of sampling point pair, need to adjust separately up and down on the basis of the center of original two symmetrical sampling points pair.For example, most
The value that big spacing is 10, n is 1-5, i 2.So, the ordinate of every a pair of symmetrical sampling point pair is in adjustment, vertical sit of above making an inventory
Numerical value after being designated as initial ordinate+(benchmark spacing-n × i)/2；The ordinate made an inventory under and is initial ordinate-(between benchmark
Away from-n × i)/2 after numerical value.Certainly, can be in adjustment successively reduces symmetrical sampling point pair at equal intervals in the manner described above
Between spacing, may not be equidistant adjustment, the objective followed in a word is pair of the central symmetry sampling point to the right and left
The spacing between sampling point pair is claimed to be gradually reduced this principle.Until spacing is between the last symmetrical sampling point pair of the right and left
Zero.
By the above-mentioned means, can guarantee that the heave of tomography polygon is more balanced, tomography is by between center
Away from maximum, the spacing to the right and left is gradually decreased, until annihilating, to meet geologic rule.It is polygon that Fig. 5 then gives tomography
The heave of shape carries out the comparative examples before and after equilibrium treatment, specific as shown in Figure 5.
Step 200, the sampling point in the second tomography polygon rendering figure is encrypted using method for resampling.
Step 210, encrypted sampling point is smoothed, obtaining final tomography polygon can chart.
Specifically, after obtaining the first tomography polygon rendering figure, in order to keep tomography polygon rendering figure more round
It is sliding, sampling number first can be encrypted using method for resampling.Increase symmetrical sampling point pair.Then to all symmetrical sampling points into
Row smoothing processing can be optionally smoothed using median filtering or mean filter method.To improve polygon
Smoothness.To enable the tomography polygon finally obtained can chart in curve it is more smooth natural.Specifically such as Fig. 6 institute
Show, Fig. 6 is the front and back comparative examples being smoothed to tomography.Curve in Fig. 6 is obviously more smooth.
A kind of tomography polygon rendering figure auto-correction method provided in an embodiment of the present invention, obtains main profile and tomography is more
After the second intersecting point coordinate between the first intersecting point coordinate between the shape of side and cross-track and tomography polygon, it is first determined disconnected
Then two extreme coordinates on layer shape changeable connect into first straight line, and obtain a plurality of normal corresponding with first straight line.And
The intersection point trend of a plurality of normal and tomography polygon more levels off to the trend of tomography, therefore utilizes a plurality of normal and tomography polygon
Intersection point draw tomography polygon rendering figure will be more accurate.Before anything is drawn, it is also necessary to using between symmetrical sampling point pair
Spaced relationship adjusts the ordinate data between different symmetrical sampling points pair, becomes the horizontal space of tomography polygon more uniformly
Change, and both ends can be realized nature and annihilate, and meet geologic rule.Finally, the second tomography polygon is encrypted using method for resampling
Sampling point in charting simultaneously is smoothed, and makes tomography that can finally move towards smooth nature.
Embodiment 2
Corresponding with embodiment 1, the embodiment of the invention also provides a kind of tomography polygon rendering figures to automatically correct dress
It sets, specifically as shown in fig. 7, the device includes: acquiring unit 701, processing unit 702, statistic unit 703 and encryption unit
704。
Acquiring unit 701, it is each in a plurality of main profile pre-seted for being obtained from original tomography polygon rendering figure
The first intersecting point coordinate between main profile and tomography polygon；And in the cross-track pre-seted each cross-track with
The second intersecting point coordinate between tomography polygon, original tomography polygon can chart as the tomography drawn by machine automatic identification
Polygon rendering figure；
The corresponding plane coordinates of two endpoints in original tomography polygon rendering figure, two of them endpoint pair are determined respectively
The plane coordinates answered is the first intersecting point coordinate and/or the second intersecting point coordinate；
Processing unit 702, for obtaining first after connecting with two endpoints according to the corresponding plane coordinates of two endpoints
The normal of the preset quantity of line correspondences；
It determines between normal that the second straight line between adjacent first intersection point and the second intersection point intersects with second straight line
Third intersection point；
When determining that third intersection point is between adjacent the first intersection point and the second intersection point, third intersection point is saved, is otherwise lost
Abandon third intersection point；
Statistic unit 703 for counting all third intersection points, and obtains the coordinate of third intersection point；By third intersection point
Coordinate be divided on symmetrical make an inventory coordinate and under make an inventory coordinate；
Processing unit 702 is also used to, determine respectively make an inventory on symmetrical coordinate and under the spacing made an inventory between coordinate, and obtain
Take maximum spacing；
By maximum spacing it is corresponding on make an inventory coordinate and under make an inventory symmetrical sampling point pair centered on coordinate, using maximum spacing as base
Quasi- spacing, after reducing benchmark spacing according to preset rules, as between the symmetrical sampling point pair in addition to central symmetry sampling point pair
Spacing adjusts the ordinate of the symmetrical sampling point pair in addition to central symmetry sampling point pair, obtains the first tomography polygon rendering figure,
In, it is gradually reduced far from the spacing of central symmetry sampling point pair, until be zero, symmetrical sampling point is to being to make an inventory on symmetrical and lower wall
The sampling point pair that point collectively constitutes；
Encryption unit 704, for encrypting the sampling point in the first tomography polygon rendering figure using method for resampling；
Processing unit 702 is also used to, and is smoothed to encrypted sampling point, and obtaining final tomography polygon can make
Figure.
Optionally, acquiring unit 701 are specifically used for, and according to the corresponding plane coordinates of two endpoints, determine that two endpoints connect
First straight line after connecing；
On the basis of each point in the default points in first straight line, preset quantity corresponding with first straight line is obtained
Normal.
Optionally, the spacing in the normal of preset quantity between adjacent two normals is identical.
Optionally, preset rules are as follows:
Using maximum spacing as benchmark spacing, after reducing n multiplied by the numerical value of i, centered on symmetrical sampling point to the left and right sides n-th
To the spacing of symmetrical sampling point pair, n be more than or equal to 1 and be less than or equal to all symmetrical sampling points to the two of total number/
The value of one positive integer, n is incremented by successively, and being incremented by spacing is 1；I is the ratio between maximum spacing and n, and is protected when adjustment
It is constant to longitudinal center to demonstrate,prove every a pair of symmetrical sampling point.
Optionally, processing unit 702 is specifically used for: smoothly being located using median filtering or mean filter to sampling point
Reason, obtains final tomography polygon rendering figure.
Function performed by each component is equal in tomography polygon rendering figure automatic correction device provided in an embodiment of the present invention
It has been described in detail in above-described embodiment 1, therefore is not described in detail here.
A kind of tomography polygon rendering figure automatic correction device provided in an embodiment of the present invention, obtains main profile and tomography is more
After the second intersecting point coordinate between the first intersecting point coordinate between the shape of side and cross-track and tomography polygon, it is first determined disconnected
Then two extreme coordinates on layer shape changeable connect into first straight line, and obtain a plurality of normal corresponding with first straight line.And
The intersection point trend of a plurality of normal and tomography polygon more levels off to the trend of tomography, therefore utilizes a plurality of normal and tomography polygon
Intersection point draw tomography polygon rendering figure will be more accurate.Before anything is drawn, it is also necessary to using between symmetrical sampling point pair
Spaced relationship adjusts the ordinate data between different symmetrical sampling points pair, becomes the horizontal space of tomography polygon more uniformly
Change, and both ends can be realized nature and annihilate, and meet geologic rule.Finally, the second tomography polygon is encrypted using method for resampling
Sampling point in charting simultaneously is smoothed, and makes tomography that can finally move towards smooth nature.
Embodiment 3
In addition, being stored thereon with computer program the embodiment of the invention also provides a kind of computer readable storage medium
Instruction, the program instruction realize the method and step of above-described embodiment 1 when being executed by processor.
Although above having used general explanation and specific embodiment, the present invention is described in detail, at this
On the basis of invention, it can be made some modifications or improvements, this will be apparent to those skilled in the art.Therefore,
These modifications or improvements without departing from theon the basis of the spirit of the present invention are fallen within the scope of the claimed invention.
Claims (10)
1. a kind of tomography polygon rendering figure auto-correction method, which is characterized in that the described method includes:
Each main profile and tomography polygon in a plurality of main profile pre-seted are obtained from original tomography polygon rendering figure
Between the first intersecting point coordinate；And in the cross-track pre-seted between each cross-track and tomography polygon second
Intersecting point coordinate, the original tomography polygon can chart as the tomography polygon rendering figure drawn by machine automatic identification；
The corresponding plane coordinates of two endpoints in the original tomography polygon rendering figure is determined respectively, wherein described two ends
The corresponding plane coordinates of point is the first intersecting point coordinate and/or the second intersecting point coordinate；
According to the corresponding plane coordinates of described two endpoints, obtain corresponding pre- with the first straight line after the connection of described two endpoints
If the normal of quantity；
It determines between normal that the second straight line between adjacent first intersection point and the second intersection point intersects with the second straight line
Third intersection point；
When determining that the third intersection point is between the first adjacent intersection point and second intersection point, the third is saved
Otherwise intersection point abandons the third intersection point；
All third intersection points are counted, and obtain the coordinate of the third intersection point；
The coordinate of the third intersection point is divided on symmetrical make an inventory coordinate and under make an inventory coordinate；
Determine respectively it is described it is symmetrical on make an inventory coordinate and under the spacing made an inventory between coordinate, and obtain maximum spacing；
By the maximum spacing it is corresponding on make an inventory coordinate and under make an inventory symmetrical sampling point pair centered on coordinate, with the maximum spacing
For benchmark spacing, after reducing the benchmark spacing according to preset rules, as symmetrical in addition to the central symmetry sampling point pair
Spacing between sampling point pair, the ordinate of the adjustment symmetrical sampling point pair in addition to the central symmetry sampling point pair obtain the
One tomography polygon rendering figure, wherein the spacing far from the central symmetry sampling point pair is gradually reduced, until be zero, it is described right
Claim sampling point to be make an inventory on symmetrical and under make an inventory the sampling point pair collectively constituted；
Sampling point in the first tomography polygon rendering figure is encrypted using method for resampling；
Encrypted sampling point is smoothed, obtaining final tomography polygon can chart.
2. the method according to claim 1, wherein described according to the corresponding plane coordinates of described two endpoints,
The normal for obtaining preset quantity corresponding with the first straight line after the connection of described two endpoints, specifically includes:
First straight line according to the corresponding plane coordinates of described two endpoints, after determining described two endpoint connections；
On the basis of each point in the default points in the first straight line, present count corresponding with the first straight line is obtained
The normal of amount.
3. the method according to claim 1, wherein in the normal of the preset quantity between adjacent two normals
Spacing it is identical.
4. method according to claim 1-3, which is characterized in that the preset rules are as follows:
Using the maximum spacing as benchmark spacing, after reducing n multiplied by the numerical value of i, as the central symmetry sampling point to left and right two
To the spacing of symmetrical sampling point pair, n is more than or equal to 1 and to be less than or equal to all symmetrical sampling points to total number for side n-th
The value of the positive integer of half, n is incremented by successively, and being incremented by spacing is 1；I is the ratio between maximum spacing and n, and is adjusted
Guarantee that every a pair of symmetrical sampling point is constant to longitudinal center when whole.
5. method according to claim 1-3, which is characterized in that described to be carried out to the encrypted sampled point
Smoothing processing, obtaining final tomography polygon can chart, and specifically include:
The sampling point is smoothed using median filtering or mean filter, obtains final tomography polygon rendering
Figure.
6. a kind of tomography polygon rendering figure automatic correction device, which is characterized in that described device includes:
Acquiring unit, for obtaining each main profile in a plurality of main profile pre-seted from original tomography polygon rendering figure
The first intersecting point coordinate between tomography polygon；And each cross-track and tomography are polygon in the cross-track pre-seted
The second intersecting point coordinate between shape, the tomography that the original tomography polygon can chart to be drawn by machine automatic identification are polygon
Shape draws figure；
The corresponding plane coordinates of two endpoints in the original tomography polygon rendering figure is determined respectively, wherein described two ends
The corresponding plane coordinates of point is the first intersecting point coordinate and/or the second intersecting point coordinate；
Processing unit, for according to the corresponding plane coordinates of described two endpoints, obtaining the after connecting with described two endpoints
The normal of the corresponding preset quantity of one straight line；
It determines between normal that the second straight line between adjacent first intersection point and the second intersection point intersects with the second straight line
Third intersection point；
When determining that the third intersection point is between the first adjacent intersection point and second intersection point, the third is saved
Otherwise intersection point abandons the third intersection point；
Statistic unit for counting all third intersection points, and obtains the coordinate of the third intersection point；By the third intersection point
Coordinate be divided on symmetrical make an inventory coordinate and under make an inventory coordinate；
The processing unit is also used to, determine respectively it is described it is symmetrical on make an inventory coordinate and under the spacing made an inventory between coordinate, and
Obtain maximum spacing；
By the maximum spacing it is corresponding on make an inventory coordinate and under make an inventory symmetrical sampling point pair centered on coordinate, with the maximum spacing
For benchmark spacing, after reducing the benchmark spacing according to preset rules, as symmetrical in addition to the central symmetry sampling point pair
Spacing between sampling point pair, the ordinate of the adjustment symmetrical sampling point pair in addition to the central symmetry sampling point pair obtain the
One tomography polygon rendering figure, wherein the spacing far from the central symmetry sampling point pair is gradually reduced, until be zero, it is described right
Claim sampling point to be make an inventory on symmetrical and under make an inventory the sampling point pair collectively constituted；
Encryption unit, for encrypting the sampling point in the first tomography polygon rendering figure using method for resampling；
The processing unit is also used to, and is smoothed to encrypted sampling point, and obtaining final tomography polygon can chart.
7. device according to claim 6, which is characterized in that the acquiring unit is specifically used for, according to described two ends
The corresponding plane coordinates of point, the first straight line after determining described two endpoint connections；
On the basis of each point in the default points in the first straight line, present count corresponding with the first straight line is obtained
The normal of amount.
8. device according to claim 6, which is characterized in that in the normal of the preset quantity between adjacent two normals
Spacing it is identical.
9. according to the described in any item devices of claim 6-8, which is characterized in that preset rules are as follows:
Using the maximum spacing as benchmark spacing, after reducing n multiplied by the numerical value of i, as the central symmetry sampling point to left and right two
To the spacing of symmetrical sampling point pair, n is more than or equal to 1 and to be less than or equal to all symmetrical sampling points to total number for side n-th
The value of the positive integer of half, n is incremented by successively, and being incremented by spacing is 1；I is the ratio between maximum spacing and n, and is adjusted
Guarantee that every a pair of symmetrical sampling point is constant to longitudinal center when whole.
10. according to the described in any item devices of claim 6-8, which is characterized in that the processing unit is specifically used for: in utilization
Value filtering or mean filter are smoothed the sampling point, obtain final tomography polygon rendering figure.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811069856.9A CN109325987B (en) | 2018-09-13 | 2018-09-13 | Automatic correction method and device for fault polygon drawing graph |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811069856.9A CN109325987B (en) | 2018-09-13 | 2018-09-13 | Automatic correction method and device for fault polygon drawing graph |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109325987A true CN109325987A (en) | 2019-02-12 |
| CN109325987B CN109325987B (en) | 2023-02-24 |
Family
ID=65265725
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811069856.9A Active CN109325987B (en) | 2018-09-13 | 2018-09-13 | Automatic correction method and device for fault polygon drawing graph |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109325987B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10969508B2 (en) * | 2016-11-18 | 2021-04-06 | Landmark Graphics Corporation | Automated MisTie analysis and correction across two-dimensional (“2D”) seismic surveys |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103500468A (en) * | 2013-09-09 | 2014-01-08 | 电子科技大学 | Space curved surface reconstruction method based on boundary constraint |
| CN104200423A (en) * | 2014-08-27 | 2014-12-10 | 大连海事大学 | Multi-control-point based automatic adjusting method for upper and lower border line shapes of vector graphics |
| CN104408755A (en) * | 2014-10-30 | 2015-03-11 | 中国电建集团成都勘测设计研究院有限公司 | Geologic plane graph drawing method based on structural plane attitude and geologic observation spot coordinate |
| CN105225273A (en) * | 2015-09-22 | 2016-01-06 | 中国石油天然气股份有限公司 | Method and the device of FAULT MODEL is set up based on space tomography polygon |
| US20160180582A1 (en) * | 2013-07-31 | 2016-06-23 | Total Sa | Method for adapting a mesh model of a geological subsurface |
-
2018
- 2018-09-13 CN CN201811069856.9A patent/CN109325987B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160180582A1 (en) * | 2013-07-31 | 2016-06-23 | Total Sa | Method for adapting a mesh model of a geological subsurface |
| CN103500468A (en) * | 2013-09-09 | 2014-01-08 | 电子科技大学 | Space curved surface reconstruction method based on boundary constraint |
| CN104200423A (en) * | 2014-08-27 | 2014-12-10 | 大连海事大学 | Multi-control-point based automatic adjusting method for upper and lower border line shapes of vector graphics |
| CN104408755A (en) * | 2014-10-30 | 2015-03-11 | 中国电建集团成都勘测设计研究院有限公司 | Geologic plane graph drawing method based on structural plane attitude and geologic observation spot coordinate |
| CN105225273A (en) * | 2015-09-22 | 2016-01-06 | 中国石油天然气股份有限公司 | Method and the device of FAULT MODEL is set up based on space tomography polygon |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10969508B2 (en) * | 2016-11-18 | 2021-04-06 | Landmark Graphics Corporation | Automated MisTie analysis and correction across two-dimensional (“2D”) seismic surveys |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109325987B (en) | 2023-02-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109800863B (en) | Logging phase identification method based on fuzzy theory and neural network | |
| Garnero et al. | Comparisons between different interpolation techniques | |
| CN106482700B (en) | A kind of digitlization House Property Area Surveying measurement method of sketch direct mapping | |
| CN105631935A (en) | BIM-based three-dimensional geological automatic modeling method | |
| CN105279790A (en) | Fracture network 3D digital core modeling method | |
| CN102759748A (en) | Seismic reservoir prediction method by cross analysis based on target analysis | |
| CN103884312B (en) | A kind of rock structural face appearance and size effect analysis method | |
| CN105654483A (en) | Three-dimensional point cloud full-automatic registration method | |
| CN112150582B (en) | Multi-modal data-oriented geological profile approximate expression method | |
| CN105093280B (en) | Surface-level model is to the low frequency of earthquake data influence and the decomposition method of radio-frequency component | |
| CN109325987A (en) | A kind of tomography polygon rendering figure auto-correction method and device | |
| CN108595749A (en) | A kind of resource reserve appraisal procedure using variation function single direction structural analysis | |
| CN103315770A (en) | Target parameter processing method and system based on three-dimensional ultrasonic image | |
| CN102783936A (en) | Method and system for determining brain state | |
| CN102353993B (en) | Method for checking closing error of velocity field during two-dimensional work area speed analysis in real time | |
| CN110610539A (en) | Stratum curved surface construction method, device, equipment and storage medium | |
| CN107507179B (en) | Rock-soil mass quantitative analysis method based on GOCAD | |
| CN104331389A (en) | Eight-point method-based contour-line tracing algorithm | |
| CN105931297A (en) | Data processing method applied to three-dimensional geological surface model | |
| CN109633781B (en) | Geological property acquisition method and device, electronic equipment and storage medium | |
| CN103853882B (en) | Drawing method and system for engineering geological profile | |
| CN108491482B (en) | Geological map dynamic synthesis method and system considering proximity relation | |
| Salah | Geostatistical analysis of groundwater levels in the south Al Jabal Al Akhdar area using GIS | |
| CN106526667B (en) | Prestack time migration method and device | |
| CN107015289B (en) | Trap evaluation stock number determines method and device |
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 |