Import Gempy Model into Revit

[jc-barreto]

Hello,
I'd like to import the geological model created with GemPy into Revit. Is it possible to export it to IFC, or are there other formats and intermediate steps that could facilitate this process?

[NilsChudalla]

That is a very exciting question! Are you interested in having surfaces or volumes?

After my quick web search, I assume, Revit is not compatible with PyVista's file formats (which is how I would usually export 3d files). You can access the vertices and edges of your surfaces and write the file, if you are familiar with the IFC file structure. I also saw a python package to write dxf files, which could be useful. Can you maybe give me a list of all the different file formats, Revit is compatible with?

Another idea came to my mind: speckle provides different plugins to exchange files between different software systems. pyvista>blender>revit would be a bit tedious, but maybe worth the try!

Edit: If your geology is simple, you could also just create georeferenced depth maps as tif. This does not work for heavily folded/faulted areas, where you have doubling of the layers. I am no revit user, so please take this with a grain of salt.

Please keep me/us updated, if any of those work for you! I am really curious!

[NilsChudalla]

Found the question so exciting, I tried it myself. I wrote a small code snippet with some gpt help. The nuances are to back transform the vertices to real coordinates. I was struggling a bit with the visualisation of dxf on my machine right now, maybe you can give me a feedback if this worked for you.

import numpy as np
import ezdxf

def write_dxf_from_triangles(vertices:dict, triangles, filename="triangles_wireframe.dxf"):
    """
    Utility function to write files with ezdf, created with help of chat gpt
    """
    doc = ezdxf.new()
    msp = doc.modelspace()

    # Store edges as lines
    edges = set()  # Use a set to avoid duplicate lines

    for tri in triangles:
        v1, v2, v3 = tri
        edges.add(tuple(sorted((v1, v2))))
        edges.add(tuple(sorted((v2, v3))))
        edges.add(tuple(sorted((v3, v1))))
    
    # Add lines to DXF
    for edge in edges:
        print(edge)
        start, end = edge
        msp.add_line(vertices[start], vertices[end])

    doc.saveas(filename)

# return the solution from compute
solution = gp.compute_model(geo_model)

# access dual contouring result
meshes = solution.dc_meshes

# iterate over surfaces in meshes
for i in range(len(meshes)):
    # retransform coordinates of vertices
    vertices = geo_model.input_transform.apply_inverse(meshes[i].vertices)

    vertice_dict = {}

    # create dictionary with enurmated vertices
    for j, ivert in enumerate(vertices):
        # consider rounding, if you have large files
        vertice_dict[j] = np.round(ivert,2)
    edges = meshes[i].edges
    
    write_dxf_from_triangles(vertice_dict, edges, filename="output_%02d.dxf"%i)
 

For the volume, you can get the "lithblock" with ids of the different surfaces. You will have to reshape it and fit it to the extent of the model:

lithblock = geo_model.solutions.raw_arrays.lith_block

As for the volume data, I will have to leave it there (edit: for today). But this could maybe solve your surfaces issue.

[NilsChudalla]

I do not have experience with the packages, that you mentioned.

[NilsChudalla]

Hi @jc-barreto ,

did you have any success with the volume data? And did my solution for surfaces help you?

[jc-barreto]

Hello @NilsChudalla Sorry, I haven't had the opportunity to test it yet. Once I do, I'll share my solution here. Thank you for your advice!

[NilsChudalla]

Hi @jc-barreto ,
you mentioned in #1001 that you are getting only lines with my code suggestion. Can you share a screenshot?

If you only want a slim volume out of the surfaces and not a volume of each unit, pyvista indeed has a potential solution. If you are familiar with blender, you could also consider automating these things (select surface out of collection, extrude faces in z direction by value x) through a small script there.

[jc-barreto]

Hi @NilsChudalla,
Sure, below is an image of one of the five surfaces exported from this example using the code you provided above. The output consists of line segments rather than polylines or meshes (triangles).

In my usual workflow, I utilize Empirical Bayesian Kriging in ArcGIS Pro, exporting the results to a NetCDF4 format for further analysis in Python. However, I have yet to import these outputs into Revit or Civil 3D. Alternatively, I often work with Geo5 and export a point cloud that needs processing in Civil 3D. I'm considering whether it's feasible to perform all these tasks solely in Python using GemPy. However, I'm concerned that this approach might result in a more time-consuming workflow. I have experience with Python, but none with Blender.

What I would like to have is something like the image below, to import into Revit/Civil 3D: