jewettaij
commented
3 years ago Here's a workaround to address the slow computation time:
The computation time can be reduced by using a mesh with a smaller number of polygons. You can reduce the number of polygons in the mesh by:
- opening the PLY file in meshlab,
- using the "Filters"->"Remeshing, Simplification, and Reconstruction"->"Quadratic Edge Collapse Decimation" menu option,
- reducing the number of polygons "faces" (eg, to less than 20000), and
- exporting the mesh as a new PLY file.
This might result in a mesh with topological problems (such as holes or cavities), so it does not always work. To reduce the chance of these kinds of problems, you might want to smooth the mesh beforehand (for example by opening the PLY file in meshlab, and selecting the "Smoothing, Fairing and Deformation"->"HC Laplacian Smoothing" menu option).
Reducing the number of polygons in the mesh makes the "voxelize_mesh.py" program run faster, but it does not seem to reduce memory consumption. I still don't know how to solve that problem.
The RAM required by this program is 25-100 times larger than the size of the original 3D image (tomogram) that we used to extract the surface mesh. This occurs because I am using 3rd party tools (pyvista, vtk) to handle the computation, instead of writing a new program from scratch. The computation is also slow. Unfortunately fixing this is not a priority for me yet -A 2020-12-15