porespy
porespy copied to clipboard
PMEAL
A set of tools for characterizing and analying 3D images of porous materials
We have quite a few functions in `tools` which are more like `utils`. The difference being that tools are used to work with images, while utils are used to help...
I find myself constantly filling blind pores before computing porosity, so realized it would helpful if the pre-built function for this actually did it for me...on demand of course.
At the moment the image 'top left' corner of the foam is returned...but it would be handy to return the centered version if possible. Not sure how much extra coding...
This would make it much easier to scale the color map. Setting vmin=04 and cmap.set_under('k') makes it look quite nice with a good depth of field, but the `sem` result...
It's time to renovate the local thickness and porosimetry functions. The changes I have in mind are backwards incompatible so this is probably more like "add to the wish list...
The [documentation](https://porespy.org/modules/generated/generated/porespy.generators.pseudo_gravity_packing.html#porespy.generators.pseudo_gravity_packing) and the [function header](https://porespy.org/_modules/porespy/generators/_pseudo_packings.html#pseudo_gravity_packing) include the argument 'edges'. Changing this argument from default value does not change anything in the output. Checking the source code revealed: it is...
The function points_to_spheres is not correctly adding spheres for a 3D image. ```python radius=10 im=np.full((41,41,41),0) im[20,20,20]=radius res_ps=ps.tools.points_to_spheres(im=im) res_fixed = points_to_spheres_fixed(im=im) print(f'Number voxels porespy: {np.sum(res_ps)}') # Number voxels porespy: 1 print(f'Number...
Use numba to scan the xy plane, and drill down in the z direction until a zero is found. This could also be parallelized.