A set of tools for analyzing 3D images of porous materials
Project description
PoreSpy
What is PoreSpy?
PoreSpy is a collection of image analysis tool used to extract information from 3D images of porous materials (typically obtained from X-ray tomography). There are many packages that offer generalized image analysis tools (i.e Skimage and Scipy.NDimage in the Python environment, ImageJ, MatLab's Image Processing Toolbox), but the all require building up complex scripts or macros to accomplish tasks of specific use to porous media. The aim of PoreSpy is to provide a set of pre-written tools for all the common porous media measurements.
Capabilities
PoreSpy consists of the following modules:
generators
: Routines for generating artificial images of porous materials useful for testing and illustrationfilters
: Functions that accept an image and return an altered imagemetrics
: Tools for quantifying properties of imagessimulations
: More complex calculations based on physical processesnetwork_extraction
: Tools for obtaining pore network representations of imagesvisualization
: Helper functions for creating useful views of the imageio
: Functions for output image data in various formats for use in common softwaretools
: Various useful tools for working with images
Installation
PoreSpy depends heavily on the Scipy Stack. The best way to get a fully functional environment is the Anaconda distribution. Be sure to get the Python 3.6+ version.
PoreSpy is available on the Python Package Index and can be installed using PIP as follows:
C:\> pip install porespy
If you think you may be interested in contributing to PoreSpy and wish to both use and edit the source code, then you should clone the repository to your local machine, and install it using the following PIP command:
C:\> pip install -e "C:\path\to\the\local\files\"
Examples
A Github repository of examples is available here. The following code snippets illustrate generating a 2D image, applying several filters, and calculating some common metrics.
Generating an image
>>> import porespy as ps
>>> import matplotlib.pyplot as plt
>>> im = ps.generators.blobs(shape=[200, 200], porosity=0.5, blobiness=2)
>>> plt.imshow(im)
Applying filters
>>> lt = ps.filters.local_thickness(im)
>>> plt.imshow(lt)
>>> cr = ps.filters.apply_chords(im)
>>> cr = ps.filters.flood(cr, mode='size')
>>> plt.imshow(cr)
Calculating metrics
>>> data = ps.metrics.two_point_correlation_fft(im)
>>> plt.plot(*data, 'b.-')
>>> mip = ps.filters.porosimetry(im)
>>> data = ps.metrics.pore_size_distribution(mip)
>>> plt.imshow(mip)
>>> plt.plot(*data, 'b.-') # Note: small image results in noisy curve
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.