Is there a way to create a watertight mesh from scalar field？

[qing-1998]

Hello,

First of all, this package is really great, and is very useful for me to create a model from drilling data .And , i was wandering is there a way to create a watertight mesh from the data: model.solutions.raw_arrays.lith_block, i have already try MC(Marching Cubes) method, but not smooth enough.

watting ur reply.

Jo

[javoha]

Hello @qing-1998,

thanks for your question. Maybe before answering we have to agree on what kind of mesh you are talking about. The surfaces meshes currently generated by gempy (basically the ones you see in the 3D plot) are created with a Dual Contouring approach. You can create similar meshes using using MC and they can be wuite smooth when you use the scalar fields from gempy to generate them. These meshes are surface meshes, meaning 2d objects.

Many people are interested in 3D volume meshes (e.g. for process simulation). These are hard to generate automatically and gempy does not have the functionaility to do that.

Which of these do you need?

[qing-1998]

Yes, i'm interested in 3D volume meshs, and also it's really hard to generate automatically. i trid to use Dual Contouring to generate the mesh , but the model shape is wrong, maybe i did a wrong step during write the method of Dual Contouring .
Could u teach me how u use Dual Contouring to generate a surface mesh? may be we can find a way to generate 3D meshs.

[javoha]

The surface meshes that gempy creates using Dual Contoruing can be accessed an used using (0 indicating the first element and so on):

geo_model.solutions.dc_meshes[0].vertices
geo_model.solutions.dc_meshes[0].edges

# vertices need to be back-transformed to original coordinate system
back_transformed_vertices = geo_model.input_transform.apply_inverse(geo_model.solutions.dc_meshes[0].vertices)

Regarding 3D volume meshes: We are working on routines to provide this option in a different project but this will not release in the super near future. If you need this urgently you should probably look into meshing software like gmsh or meshio.

Cheers,
jan

[NilsChudalla]

As long as you are not having a very faulted model, you could consider checking out GeoMeshPy. Please note: We are not maintaining or developing this code, so you will have to ask the author in the case of questions.

[qing-1998]

Thanks for your advice @javoha and @NilsChudalla
I will try the way as u say and will reply latter, thanks again!
Have a good day!
Jo