How to specify label to mesh in 3D array example? about pygalmesh HOT 2 CLOSED

@thoth291:

The domain to be meshed consists of only the voxels with label strictly bigger than 0.

That is, let arr be the 3D image with labels: pore space (0), material space (1), fuild(2).

Then if you want to mesh only the material region simply define:

arr_material = arr
arr_material[arr==2] = 0

Now arr_material will have non-zero labels only to the material region, so you'll get a mesh only for that region.

Note:
If you want to mesh all three regions, please make sure to increment each label by 1.

arr_new = arr + 1

Hope this help.

from pygalmesh.

The domain to be meshed consists of only the voxels with label strictly bigger than 0.

This is exactly what I was looking for !
Thank you for this confirmation!

from pygalmesh.