This is a tool that fits Gaussian Mixture Models to grey value distributions from X-ray micro-computed tomography images. This tool was developed as part of a PhD project at the University of Southampton in the 3-D X-ray histology team.
GaussQuality fits a 1-D Gaussian Mixture Model (GMM) to the grey value distribution of a 2D X-ray micro-computed tomography (μCT) image. The GMM fits a specified number of Gaussian components to the grey value distribution, where each component represents the distribution of grey values for a single material in the specimen. The mean (μ), standard deviation (σ), and weight (φ) of each Gaussian component are estimated.
- μ is the location of the Gaussian along the grey value axis, dependent on X-ray interactions with the material
- σ is its spread which is dependent on heterogeneity of the material and image noise
- φ is the proportion of the image taken up by this material
μ and σ can be used to calculate signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR).
μ, σ, and φ can be estimated for several 2D images in a 3D stack, showing how the image varies in 3D. Only one 2D image is held in memory at once, so images larger than memory can be processed. 3D images must be saved as sequences of 2D images.
Download the .exe for Windows or .app for MacOS.
Simply double click the executable and the GUI should start. If you have any trouble, try opening the executable in a terminal to debug any errors.
GaussQuality is available on conda. In the terminal:
conda install -c elainehoml gaussquality
Alternatively, you can install GaussQuality with pip. In the terminal:
pip install gaussquality
It's best to install GaussQuality into a fresh virtual environment to avoid conflicts with other projects you may have.
See GaussQuality_Example.ipynb for a demo, and explanation of what each parameter controls.