CN100440258C - System and method for automatic creating data gridding of stratum and fault - Google Patents
System and method for automatic creating data gridding of stratum and fault Download PDFInfo
- Publication number
- CN100440258C CN100440258C CNB2006101521270A CN200610152127A CN100440258C CN 100440258 C CN100440258 C CN 100440258C CN B2006101521270 A CNB2006101521270 A CN B2006101521270A CN 200610152127 A CN200610152127 A CN 200610152127A CN 100440258 C CN100440258 C CN 100440258C
- Authority
- CN
- China
- Prior art keywords
- line
- rent
- stratum
- point
- grid
- 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.)
- Expired - Fee Related
Links
Images
Abstract
A method for automatically generating data lattice of stratum and fault includes fetching contur and fault line data, setting up uniform lattice height field as per surrounding box of input data, conforming fault region as per naming rule of fault line, confirming connection relation of lattice point as per regions passed by fault line, confirming height value of each lattice point in height field by utilizing secondary local growth manner, matching lattice points of multilayer fault line, plotting out section condition of fault in multilayer stratum and stratum fault on display.
Description
Technical field
The present invention relates to the system and method that a kind of stratum and layer data grid generate automatically.
Background technology
In recent years, in field of Computer Graphics, along with the introducing of the development of photo realism graphic generation technique and fractal geometry, wavelet analysis scheduling theory, people can utilize computer simulation to generate many natural scenes that are really true to life (comprising landform, plant, cloud, fire, wave, meadow, fountain etc.).In the virtual physical environment of flight simulator, military training simulator, ship maneuvering simulator etc., one of core technology is visual simulation, and dimensional topography is a virtual physical environment scenery the most basic in the what comes into a driver's, is the important component part of visual system.Because actual landform rises and falls and to vary, make that the data volume of three dimension realistic terrain simulation is very big, therefore aspect modeling and the demonstration two in real time all the graphics capability to computing machine very high requirement has been proposed.The development of computer graphical hardware accelerator has been played crucial effect to the real-time demonstration of dimensional topography, but in the actual landform modeling, required graph data amount is still than the hardware big one or more orders of magnitude of data presented amount in real time, the memory space that on the landform of 512*512 precision, just needs tens G as the bidirectional reflectance distribution function (BRDF) that proposes for the material that shows real world effectively and texture and bidirectional texture function (BTF), and the complexity of application model also usually surpasses the actual treatment ability of current graphics workstation.Therefore, in actual applications, attempt usually by algorithm and software engineering, on active computer performance and level of hardware, take all factors into consideration the relation between display speed and the generation graphical effect, to generate the three dimensional terrain figure that satisfies different geometric accuracies and fidelity requirement.A kind of as Lindstrom proposition in 1996 based on uniform grid, and the real-time landform rendering algorithm of the continuous level of detail relevant with viewpoint, they adopt quaternary tree to manage the landform scene, and adopt screen error decision condition to set up the tree node evaluation function.Hoppe was with its viewpoint associated change grid (View Dependent Progressive Meshes in 1998, VDPM) be applied in the landform, the real-time rendering and the roaming of large-scale terrain have been realized, they are divided into bulk with landform earlier, carry out pre-service with the thought that changes grid respectively to every, the information of record simplification process generates triangular mesh by these information in real time at different viewpoint positions.Losasso in 2004 etc. propose a kind of real-time landform rendering algorithm based on LOD, this algorithm adopts a kind of ring data structure, make move in viewpoint, when rotating, can upgrade terrain data efficiently, terrain data uses Vertex Buffer Object expansion to draw, this algorithm also compresses terrain data, makes that the terrain data of original 40G also can be as the real time roaming scene.
Level line is the more special and important line of a class in the terrain data, and it has expressed the projected footprint of point on surface level highly identical on the landform.Can be regarded as a special case of carrying out three-dimensionalreconstruction by the point in the three dimensions by level line reconstruct height field, and the existing a lot of achievements of this respect.Existent method has distance inverse ratio method, PDE method and method of conjugate gradient.And in many practical applications, because the relation of data volume and calculated amount, general layer model is simplifiedly all only handled level line and is not considered the zone of tomography.In existing method, the distance inverse ratio method is the reconstruct level line well, because the information of line of rent can not effectively be analyzed and handle to this method, and the result that PDE method and method of conjugate gradient can access, but the people can't be accepted for speed.Thereby the dimensional topography that these methods are set up is not owing to there being the information of tomography, so be not suitable for fields such as subsurface investigation and measurement.
Summary of the invention
Therefore, need a kind of automatic generation stratum to comprise the system and method for layer data grid in this technical field, this system and method can comprise tomography with true stratum and show on computer display screen apace.
For solving the problems of the technologies described above, the present invention at first uses for reference region growing method classical in the computer vision, and is aided with the surface tension threshold value, with the level line present position is initial value, line of rent is the border, once increases, and the part that increases gained is the height field part.Then, the present invention has overcome in the distance inverse ratio method because the existence of line of rent, the problem that point may pass through line of rent to isocontour vertical range, adopt " geodesic distance " to replace vertical range, with the distance inverse ratio method reconstruct height field after this improvement, and continue to use the region growing method, allow every level line increase as an initial entity, two nearest isocontour growth degree of depth under each some storage are carried out interpolation again.After second zone increases, obtained a height value more accurately on each point, the value of this value as final height field also compared accurately.But in order to obtain more level and smooth surface and height field more accurately, the height field value that the present invention will generate is used method of conjugate gradient optimization as initial value, because given initial value is more accurately restrained very fast.So just obtained the cross section situation of stratum interrupting layer, at last drawn in stratum and tomography on display
As another kind of improvement the of the present invention, because traditional level line generates the situation that the height field algorithm is not considered tomography, the present invention will put isocontour distance and redefine minimum fold line distance (being " geodesic distance " among the present invention) for not cut off by line of rent, can continue like this to use original algorithm to handle the height field that contains layer data, thereby keep former efficiency of algorithm height, fireballing characteristic.
The invention discloses the method that a kind of stratum and layer data grid generate automatically, this method comprises the steps:
The fault region determining step is set up uniform grid height field according to the input data, in conjunction with the line of rent Naming conventions, determines fault region;
The net point Connection Step, the grid that passes according to line of rent is determined the annexation of net point;
Grid point height value determining step, adopt the second zone growth method to determine each grid point height value, described second zone growth method is meant that employing " geodesic distance " replaces vertical range, with distance inverse ratio method reconstruct height field, the height field value that will generate is again used method of conjugate gradient optimization as initial value;
Plot step is mated the line of rent net point of multilayer, draws out the situation in stratum tomography and interrupting layer cross section, multilayer stratum at last on display.
In said method, can adopt the relief block of uniform grid to handle stratum and layer data.
In said method, can adopt from the classic method of level line generation height field and set up height field.
In said method, defining point is to the minimum fold line distance of isocontour geodesic distance for not cut off by line of rent.
In said method, preestablish the size of the region growing degree of depth.
In the plot step of said method, the normal direction value that can use the net point surface shows as the color of this point.
The invention also discloses the system that a kind of stratum and layer data grid generate automatically, this system comprises:
Fault region is determined device, and this device is set up uniform grid height field according to the input data, in conjunction with the line of rent Naming conventions, determines fault region;
The net point coupling arrangement, the grid that this device passes according to line of rent is determined the annexation of net point;
The grid point height value is determined device, this device adopts the second zone growth method to determine each grid point height value, described second zone growth method is meant that employing " geodesic distance " replaces vertical range, with distance inverse ratio method reconstruct height field, the height field value that will generate is again used method of conjugate gradient optimization as initial value;
Drawing apparatus, this device mates the line of rent net point of multilayer, draws out the situation in stratum tomography and interrupting layer cross section, multilayer stratum at last on display.
In said system, fault region determines that device has adopted the classic method that generates height field from level line to set up height field.
In said system, fault region determines that device can preestablish the size of the region growing degree of depth.
In said system, the normal direction value on drawing apparatus use net point surface shows as the color of this point.
Description of drawings
The present invention is further detailed explanation below in conjunction with the drawings and specific embodiments:
Fig. 1 is a structure flow chart of the present invention.
The uniform grid model synoptic diagram of Fig. 2 for adopting among the present invention.
Fig. 3 is for determining the synoptic diagram of fault region.
Fig. 4 is for redefining among the present invention a little to isocontour geodesic distance synoptic diagram.
Fig. 5 cuts off the synoptic diagram that neighbor mesh points connects for line of rent.
Fig. 6 determines the synoptic diagram of grid point height value for the second zone growth method.
Fig. 7 is the grid coupling synoptic diagram between the line of rent on multilayer stratum.
Fig. 8 imports data and output result's design sketch for the present invention.
Embodiment
Fig. 1 is a structure flow chart of the present invention.The input data are the formation data of multilayer, with Autocad R12/LT2dxf representation of file.The present invention can obtain contour line data and line of rent data by reading corresponding dxf file.Line of rent is divided into two kinds: there are not the line of rent and the line of rent that coupling is arranged (promptly in same stratum, two line of rent are corresponding mutually, and the zone between these two line of rent is fault region) of coupling, and can be referring to the hatched example areas among Fig. 3 (10).Be processing and the coupling that makes things convenient for line of rent, the present invention has done following standard to the storage that the dxf file interrupts layer line: every line of rent takies a layer in the dxf file separately, and the name of this layer must start with FAULTAGE.For the line of rent that does not have coupling, the called after FAULTAGE of line of rent (capitalization).As FAULTAGE_A, FAULTAGE_B.For the line of rent that coupling is arranged, the name of two line of rent is respectively FAULTAGE (capitalization) 1 and FAULTAGE (same capitalization) 2, and certain appearance in pairs.As FAULTAGE_C1 and FAULTAGE_C2.Between the Different Strata, the line of rent of same name mates automatically.As above the FAULTAGE_A in the FAULTAGE_A in the stratum and next stratum mates automatically, generates the tomography cross section.Bounding box according to the input data is set up uniform grid height field (11), in conjunction with our line of rent Naming conventions, determines fault region (12).The grid that passes according to line of rent is determined the annexation (13) of net point then.Adopt the second zone growth method to determine each grid point height value (14).Line of rent net point to multilayer mates (15), draws out the situation (16) in stratum tomography and interrupting layer cross section, multilayer stratum at last on display.
The uniform grid model synoptic diagram of Fig. 2 for adopting among the present invention.In the uniform grid model, data are represented by the height value of storing each net point on the square uniform grid.The height value of other positions should obtain by the method interpolation of nearest four the grid point height values in this position with bilinear interpolation.
Fig. 3 is for determining the synoptic diagram of fault region.Line of rent after the digitizing divides two kinds: a kind of not coupling, as line of rent (1); A kind of have a coupling, as line of rent (2 and 3), (4 and 5), and (6 and 7), (8 and 9) all have the corresponding of two couplings; Remaining lines is level line.The zone of coupling in the line of rent is fault region, and hatched example areas (10) is line of rent (8) and line of rent (9) the determined fault region afterwards that is complementary.
Fig. 4 is for redefining among the present invention a little to isocontour geodesic distance synoptic diagram.The present invention has adopted when the computed altitude field and a kind ofly easier has generated the height field algorithm from level line: for any 1 P, h1 and h2 are respectively from two nearest level lines of P point, d1 and d2 are respectively P o'clock to two isocontour distance, determine the height value that P is ordered with following formula so:
But, for the situation that line of rent is arranged, lines among Fig. 4 (43) and lines (44) are line of rent, line of rent separates the straight-line segment d1 of P point to level line h1, therefore, we redefine d1 and d2 for P o'clock to two isocontour geodesic distance (i.e. the minimum fold line distance of not cut off by line of rent), can still determine the height value that P is ordered with above-mentioned formula.
Fig. 5 cuts off the synoptic diagram that neighbor mesh points connects for line of rent.Line of rent (50) can cut off the connection of neighbor mesh points through part, as net point (51) and net point (52), net point (53) and net point (54), net point (53) and net point (55), net point (56) and net point (57), net point (58) and net point (59), net point (51 ') and net point (52 '), net point (52 ') and net point (53 ').Two grid end points are cut off by line of rent, make that two net points can not influence each other when the region growing method generates height.Can influence each other during all the other adjacent mesh dot generation height.
Fig. 6 determines the synoptic diagram of grid point height value for the second zone growth method.The region growing method is a kind of method that is used for color identification in the computer vision.We revise region growing method of the prior art among the present invention, and be used for computing grid and put isocontour geodesic distance, and then the height value of definite net point.Method step: with every level line is an entity, net point with this level line place is a starting point, the net point that is connecting is carried out region growing, and note nearest with it two level lines and, utilize the geodesic distance interpolation to obtain the height value of this net point to two isocontour geodesic distances (promptly increasing the degree of depth) at each net point.In addition, given one increases degree of depth limits value, and promptly every isocontour maximum region increases depth value, reaches this value if increase the degree of depth, promptly stops to increase.
Be purposes of simplicity of explanation, have only a line of rent and two level lines among Fig. 6 altogether, wherein: lines (60) are line of rent, and lines (61) are two different level lines with lines (62).From level line, increasing of a point of a point like this, can be noted the degree of depth of growth, i.e. geodesic distance at each net point place.Each the net point below that shows among Fig. 6 and the value of top are the geodesic distance that this is put level line (61) and arrives level line (62).By it can also be seen that among the figure, if predefined growth degree of depth limits value is 3, then the point outside the hachure (63) is no defined range.
Fig. 7 is the grid coupling synoptic diagram between the line of rent on multilayer stratum.Between the multilayer stratum, need between the tomography of coupling, draw the cross section.Owing between the stratum of different layers, specified the line of rent of coupling, got final product so only need between the line of rent of coupling, to draw triangular facet by the title of line of rent.Do suitable triangle division according to the length of line of rent, draw these triangular facets then and get final product.Line of rent among Fig. 7 (70) has 5 summits (70 ', 71 ', 72 ', 73 ', 74 '), line of rent (71) has 3 summits (75 ', 76 ', 77 '), we can access following a kind of triangle division: triangle (70 ', 71 ', 75 '), triangle (71 ', 75 ', 76 '), triangle (71 ', 72 ', 76 '), triangle (72 ', 76 ', 77 '), triangle (72 ', 73 ', 77 '), triangle (73 ', 74 ', 77 '), and the requirement that this triangle division satisfy to be drawn fully.
Fig. 8 imports data and output result's design sketch for the present invention.The first half among the figure (81) is the level line and the line of rent design sketch of input, and the latter half among the figure (82) is the 3 d effect graph that output is drawn by system.Not only show among the figure and defer to the isocontour visualization of 3 d landform of input, also clearly demonstrate the three-dimensional model of tomography, meet the line of rent of input fully, this is that reconstructing method institute was inaccessible in the past.Generate after the height field, for each net point, those skilled in the art can know: this net point is positioned at fault region, the height field zone still is the height value of undefined zone and this net point.For undefined zone, will not show; For fault region, a kind of color of unified use shows, is fault region to indicate this zone; For the height field zone, be the concave-convex sense on performance stratum, the normal direction value on system use net point surface shows as the color of this point.As shown in Figure 8, those skilled in the art can see that advantage of the present invention is the grid that can generate stratum and layer data automatically, and the effect in 3-D display tomography cross section apace.
Claims (10)
1. the method that generates automatically of stratum and layer data grid is characterized in that this method comprises the steps:
The fault region determining step is set up uniform grid height field (11) according to the input data, in conjunction with the line of rent Naming conventions, determines fault region (12);
The net point Connection Step, the grid that passes according to line of rent is determined the annexation (13) of net point;
Grid point height value determining step, adopt the second zone growth method to determine each grid point height value (14), described second zone growth method is meant that employing " geodesic distance " replaces vertical range, with distance inverse ratio method reconstruct height field, the height field value that will generate is again used method of conjugate gradient optimization as initial value;
Plot step is mated (15) to the line of rent net point of multilayer, draws out the situation (16) in stratum tomography and interrupting layer cross section, multilayer stratum at last on display.
2. according to the method for claim 1, it is characterized in that: adopted the relief block of uniform grid to handle stratum and layer data (21 and 22).
3. according to the method for claim 1, it is characterized in that: adopted from the classic method of level line generation height field and set up height field.
4. according to the method for claim 3, it is characterized in that: defining point is to the minimum fold line distance (41 and 42) of isocontour geodesic distance for not cut off by line of rent.
5. according to the method for claim 3, it is characterized in that: the size that preestablishes the region growing degree of depth.
6. according to the method for claim 1, it is characterized in that: the normal direction value on use net point surface shows as the color of this point in plot step.
7, the automatic system that generates of a kind of stratum and layer data grid is characterized in that this system comprises:
Fault region is determined device, and this device is set up uniform grid height field (11) according to the input data, in conjunction with the line of rent Naming conventions, determines fault region (12);
Net point coupling arrangement, the grid that this device passes according to line of rent are determined the annexation (13) of net point;
The grid point height value is determined device, this device adopts the second zone growth method to determine each grid point height value (14), described second zone growth method is meant that employing " geodesic distance " replaces vertical range, with distance inverse ratio method reconstruct height field, the height field value that will generate is again used method of conjugate gradient optimization as initial value;
Drawing apparatus, this device mates (15) to the line of rent net point of multilayer, draws out the situation (16) in stratum tomography and interrupting layer cross section, multilayer stratum at last on display.
8. system according to claim 7 is characterized in that: fault region determines that device has adopted the classic method that generates height field from level line to set up height field.
9, system according to claim 8 is characterized in that: fault region determines that device can preestablish the size of the region growing degree of depth.
10, system according to claim 7 is characterized in that: the normal direction value on drawing apparatus use net point surface shows as the color of this point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2006101521270A CN100440258C (en) | 2006-09-14 | 2006-09-14 | System and method for automatic creating data gridding of stratum and fault |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2006101521270A CN100440258C (en) | 2006-09-14 | 2006-09-14 | System and method for automatic creating data gridding of stratum and fault |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1916970A CN1916970A (en) | 2007-02-21 |
| CN100440258C true CN100440258C (en) | 2008-12-03 |
Family
ID=37737960
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB2006101521270A Expired - Fee Related CN100440258C (en) | 2006-09-14 | 2006-09-14 | System and method for automatic creating data gridding of stratum and fault |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN100440258C (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU2014306485B2 (en) * | 2013-08-16 | 2016-12-08 | Landmark Graphics Corporation | Generating a custom reservoir from multiple compartments representing one or more geological structures |
| CN106097445B (en) * | 2016-06-02 | 2019-10-29 | 广州市设计院 | A kind of method for drafting of three-dimensional formation curved surface |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0950538A (en) * | 1995-05-31 | 1997-02-18 | Nec Corp | Three-dimensional data preparing system from contour |
| JPH1076043A (en) * | 1996-09-04 | 1998-03-24 | P G P Internatl:Kk | Simulation method and device for golf course design |
| JP2000268195A (en) * | 1999-03-18 | 2000-09-29 | Ricoh Co Ltd | Method for compressing and generating curved surface of polygon data |
| JP2002048524A (en) * | 2000-08-04 | 2002-02-15 | Inst Of Physical & Chemical Res | Method and device for plane interpolation of measured data |
-
2006
- 2006-09-14 CN CNB2006101521270A patent/CN100440258C/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0950538A (en) * | 1995-05-31 | 1997-02-18 | Nec Corp | Three-dimensional data preparing system from contour |
| JPH1076043A (en) * | 1996-09-04 | 1998-03-24 | P G P Internatl:Kk | Simulation method and device for golf course design |
| JP2000268195A (en) * | 1999-03-18 | 2000-09-29 | Ricoh Co Ltd | Method for compressing and generating curved surface of polygon data |
| JP2002048524A (en) * | 2000-08-04 | 2002-02-15 | Inst Of Physical & Chemical Res | Method and device for plane interpolation of measured data |
Also Published As
| Publication number | Publication date |
|---|---|
| CN1916970A (en) | 2007-02-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Nishida et al. | Example‐driven procedural urban roads | |
| Takayama et al. | Lapped solid textures: filling a model with anisotropic textures | |
| CN109360260A (en) | A kind of cut-away restructing algorithm of triangle gridding threedimensional model | |
| Floriani et al. | Variant: A system for terrain modeling at variable resolution | |
| CN102044089A (en) | Method for carrying out self-adaption simplification, gradual transmission and rapid charting on three-dimensional model | |
| CN104778744B (en) | Extensive three-dimensional forest Visual Scene method for building up based on Lidar data | |
| Zhao | Application of 3D CAD in landscape architecture design and optimization of hierarchical details | |
| CN113593051A (en) | Live-action visualization method, dam visualization method and computer equipment | |
| Ribelles et al. | Multiresolution modeling of arbitrary polygonal surfaces: a characterization | |
| US20040201584A1 (en) | Spatial decomposition methods using bit manipulation | |
| Moran et al. | Visualization of AMR data with multi-level dual-mesh interpolation | |
| CN110349262A (en) | A kind of tunnel geological forecast three-dimensional modeling method and system based on BIM | |
| CN100440258C (en) | System and method for automatic creating data gridding of stratum and fault | |
| CN109684656B (en) | Assembly constraint inheritance method based on SolidWorks | |
| Amiri et al. | Connectivity maps for subdivision surfaces | |
| Chang et al. | Texture Tiling on 3D Models Using Automatic Polycube-maps and Wang Tiles. | |
| CN113096255A (en) | WebGL-based three-dimensional view synchronous fusion non-invasive integration technical method | |
| Fu et al. | Development and Application of Marine Environment Data Visualization Technology | |
| Shakaev et al. | View-Dependent Level of Detail for Real-Time Rendering of Large Isosurfaces | |
| Pan et al. | Level of detail and multi-resolution modeling for virtual prototyping | |
| Paris et al. | Modeling rocky scenery using implicit blocks | |
| Hossain et al. | Environment and object design for 3D simulation in context of commercial vehicles | |
| Westerteiger | Virtual Reality Methods for Research in the Geosciences | |
| Cruz et al. | High-level techniques for landscape creation | |
| Wei et al. | Research on a General Virtual Visualization System of Finite Element Analysis Result Data of Ship Structure |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20081203 Termination date: 20170914 |