Converting corrupted 3d geomodels into voxel mode.

B. Wieczoreck and Georg Semmler and Jan Gietzel and Helmut Schaeben. ( 2015 )
in: 35th Gocad Meeting - 2015 RING Meeting, ASGA

Abstract

Visualization, statistical analysis and numerical process simulation are just a few amongst a wide range of applications of 3d digital geomodels. Depending on the specified purpose and application, one of the various representations of a geomodel may be more appropriate than any other. The obvious representations are (i) triangulated geological surfaces, horizons or faults (see Figure 1), (ii) geological bodies in terms of their boundaries (BREPs), (iii) a solid partitioned in tetrahedrons, and (iv) a solid partitioned in regular hexahedrons, which we usually refer to as voxels. The first three representations may be the result of the computer aided construction of the model itself, while the voxel representation is usually the result of an additional process of “voxelization” applied to the initial geomodel. With respect to geographic information systems, vector and raster data models, respectively, are classics. Their representation as well as their storage in database systems is standardized. If the geomodel is “watertight”, ray casting, a generalization of the “point-in-polygon” algorithm, applies as a corollary of the Jordan curve theorem. It does not apply if the representation of the initial geomodel is corrupted and not watertight. For the case that the geomodel is not watertight, this contribution introduces a new algorithm. In the first step, triangulated surfaces are converted into sets of voxels. Then voxels along columns of the geomodel are assigned to bodies of the geomodel whenever immediately possible. Voxels which are not assignable are tagged. Assigning initially tagged voxels requires additional provisions. It may not be possible to resolve all ambiguities due to the corruption of the initial geomodel, and the affected voxels remain tagged. Independently of the way the representation of a geomodel has been converted into voxel mode, the raster data manager “rasdaman” [Misev and Baumann [2014]] is used to store the geomodel in voxel mode. “rasdaman” features OGC standard compliant output formats like WCS or WMS. Then the voxel mode of 3d subsurface geomodels allows for statistical analyses and numerical processing through OGC interfaces.

Download / Links

BibTeX Reference

@inproceedings{Wieczoreck_GM2015,
 abstract = { Visualization, statistical analysis and numerical process simulation are just a few amongst a wide range of applications of 3d digital geomodels. Depending on the specified purpose and application, one of the various representations of a geomodel may be more appropriate than any other. The obvious representations are (i) triangulated geological surfaces, horizons or faults (see Figure 1), (ii) geological bodies in terms of their boundaries (BREPs), (iii) a solid partitioned in tetrahedrons, and (iv) a solid partitioned in regular hexahedrons, which we usually refer to as voxels. The first three representations may be the result of the computer aided construction of the model itself, while the voxel representation is usually the result of an additional process of “voxelization” applied to the initial geomodel. With respect to geographic information systems, vector and raster data models, respectively, are classics. Their representation as well as their storage in database systems is standardized. If the geomodel is “watertight”, ray casting, a generalization of the “point-in-polygon” algorithm, applies as a corollary of the Jordan curve theorem. It does not apply if the representation of the initial geomodel is corrupted and not watertight. For the case that the geomodel is not watertight, this contribution introduces a new algorithm. In the first step, triangulated surfaces are converted into sets of voxels. Then voxels along columns of the geomodel are assigned to bodies of the geomodel whenever immediately possible. Voxels which are not assignable are tagged. Assigning initially tagged voxels requires additional provisions. It may not be possible to resolve all ambiguities due to the corruption of the initial geomodel, and the affected voxels remain tagged. Independently of the way the representation of a geomodel has been converted into voxel mode, the raster data manager “rasdaman” [Misev and Baumann [2014]] is used to store the geomodel in voxel mode. “rasdaman” features OGC standard compliant output formats like WCS or WMS. Then the voxel mode of 3d subsurface geomodels allows for statistical analyses and numerical processing through OGC interfaces. },
 author = { Wieczoreck, B. AND Semmler, Georg AND Gietzel, Jan AND Schaeben, Helmut },
 booktitle = { 35th Gocad Meeting - 2015 RING Meeting },
 publisher = { ASGA },
 title = { Converting corrupted 3d geomodels into voxel mode. },
 year = { 2015 }
}