Skip to main content

<!– PROJECT SHIELDS –>

[![CircleCI](https://circleci.com/gh/PolymerGuy/AXITOM.svg?style=svg)](https://circleci.com/gh/PolymerGuy/AXITOM) [![codecov](https://codecov.io/gh/PolymerGuy/AXITOM/branch/master/graph/badge.svg)](https://codecov.io/gh/PolymerGuy/AXITOM) [![MIT License][license-shield]][license-url] [![Documentation Status](https://readthedocs.org/projects/axitom/badge/?version=latest)](https://axitom.readthedocs.io/en/latest/?badge=latest) [![PyPI version](https://badge.fury.io/py/axitom.svg)](https://badge.fury.io/py/axitom)

<!– PROJECT LOGO –> <br /> <p align=”center”> <img src=”./docs/logo.png” alt=”Logo” width=”500” height=”500”> </p>

<h3 align=”center”>AXITOM</h3>

<p align=”center”> Tomographic reconstruction of axisymmetric fields acquired by a cone beam <br /> <a href=”https://axitom.readthedocs.io/en/latest/”><strong>Explore the docs </strong></a> <br /> </p>

<!– ABOUT THE PROJECT –> About The Project —————– This python package provides tools for axis-symmetric cone-beam computed tomography. A Feldkamp David Kress algorithm performs the reconstruction which have been adapted such that is exploits the axis-symmetric nature of the tomogram.

This toolkit has a highly specialised usage, and there are plenty of more general and excellent frameworks for tomographic reconstruction, such as:

This project aims at providing a simple, accessible toolkit for forward-projection and reconstruction of axis-symmetric tomograms based on a conical beam geometry.

### Built With This project is heavily based on the following packages: * [Numpy](https://numpy.org/) * [Scipy](https://www.scipy.org/) * [Scikit-image](https://scikit-image.org/) * [Matplotlib](https://matplotlib.org/)

<!– GETTING STARTED –> Getting Started ————— To get a local copy up and running follow these steps.

### Install via package manager:

The toolkit is available via PIP, and the instructions below shows how a virtual environment can be created and the toolkit installed.

Prerequisites:

This toolkit is tested on Python 3.6 We recommend the use of virtualenv

Installing:

$ virtualenv -p python3.6 env $ source ./env/bin/activate #On Linux and Mac OS $ envScriptsactivate.bat #On Windows $ pip install axitom

Now the toolkit is installed and ready for use.

Run the tests:

$ nosetests axitom

If you want to check out the examples, then download the files in the examples folder and run the examples.

### Clone the repo:

These instructions will get you a copy of the project up and running on your local machine for development and testing purposes.

Prerequisites:

This toolkit is tested on Python 3.6 We recommend the use of virtualenv

Clone this repo to your preferred location

$ git init $ git clone https://github.com/PolymerGuy/axitom.git

We recommend that you always use virtual environments, either by virtualenv or by Conda env

$ virtualenv -p python3.6 env $ source ./env/bin/activate #On Linux and Mac OS $ envScriptsactivate.bat #On Windows $ pip install -r requirements.txt

You can now run an example:

$ python <path_to_axitom>/examples/comparison_to_Nikon.py

### Run the tests The tests should always be launched to check your installation. These tests are integration and unit tests.

If you cloned the repo, you have to call pytest from within the folder

$ pytest

Example

Let us now go through the necessary steps for doing a reconstruction of a tomogram based on a single image. First, we need to import the tools

import axitom as tom from scipy.ndimage.filters import median_filter

The example data can be downloaded from the AXITOM/tests/example_data/ folder. The dataset was collected during tensile testing of a polymer specimen. Assuming that the example data from the repo is located in root folder, we can make a config object from the .xtekct file

config = tom.config_from_xtekct(“radiogram.xtekct”)

We now import the projection

projection = tom.read_image(r”radiogram.tif”, flat_corrected=True)

As we will use a single projection only in this reconstruction, we will reduce the noise content of the projection by employing a median filter. Using such a filter works fine since the density gradients within the specimen are relatively small. You may here choose any filter of your liking.

projection = median_filter(projection, size=21)

Now, the axis of rotation has to be determined. The axis of rotation is found by first binarizing of the image into object and background, and subsequently determining the centre of gravity of the object

_, center_offset = tom.object_center_of_rotation(projection, background_internsity=0.9)

The config object has to be updated with the correct values

config = config.with_param(center_of_rot=center_offset)

We are now ready to initiate the reconstruction

tomo = tom.fdk(projection, config)

The results can then be visualized

import matplotlib.pyplot as plt plt.title(“Radial slice”) plt.imshow(tomo.transpose(), cmap=plt.cm.magma)

<img src=”./docs/results.png” alt=”Results” width=”300”/>

<!– CONTRIBUTING –> Contributing ————

Contributions are what makes the open-source community such a fantastic place to learn, inspire, and create. Any contributions you make are greatly appreciated.

  1. Fork the Project

  2. Create your Feature Branch (git checkout -b feature/AmazingFeature)

  3. Commit your Changes (git commit -m ‘Add some AmazingFeature)

  4. Push to the Branch (git push origin feature/AmazingFeature)

  5. Open a Pull Request

<!– LICENSE –> License ——-

Distributed under the MIT License. See LICENSE for more information.

<!– CONTACT –> Contact ——-

Sindre Nordmark Olufsen (PolymerGuy) - sindre.n.olufsen@ntnu.no

<!– ACKNOWLEDGEMENTS –> Acknowledgements —————- We are in great debt to the open-source community and all the contributors the projects on which this toolkit is based.

<!– MARKDOWN LINKS & IMAGES –> [license-shield]: https://img.shields.io/badge/license-MIT-blue.svg?style=flat-square [license-url]: https://choosealicense.com/licenses/mit

Release files for axitom 0.1.3

For a detailed explanation of source distributions (sdists) and built distributions (wheels), please see the package formats documentation.

Source distribution (sdist)

Source distribution for axitom 0.1.3
File Size Uploaded
axitom-0.1.3.tar.gz 19.0 MB Details

Built distribution (wheel)

Table of built distributions (wheels) for axitom 0.1.3
File Interpreter ABI Platform
axitom-0.1.3-py2-none-any.whl Python 2 none any Details

Total release size: 38.0 MB

Release files / axitom-0.1.3.tar.gz

Download URL axitom-0.1.3.tar.gz
Size 19.0 MB
Tags Source
SHA-256 checksum
How to use checksums
1bd9549d1e720f68fe414f2dd9105ad244589d3980f0f6da22b1f08f9cce9a6b
BLAKE2b-256 checksum
How to use checksums
8d828ab261d228489880fd4f8a35deeba4377af7e8e6498f997459f184ec48ed
Upload date
Uploaded using Trusted Publishing?
What is trusted publishing?
No
Uploaded via twine/2.0.0 pkginfo/1.5.0.1 requests/2.22.0 setuptools/41.2.0 requests-toolbelt/0.9.1 tqdm/4.36.1 CPython/3.6.8

Release files / axitom-0.1.3-py2-none-any.whl

Download URL axitom-0.1.3-py2-none-any.whl
Size 19.0 MB
Tags Python 2
SHA-256 checksum
How to use checksums
cbdfc9208896cc3a2d0c33b504240c8a9fae65363d69a962c21fbca257eb5f28
BLAKE2b-256 checksum
How to use checksums
916c6ea3d51b0d200e9fc3a0fa257e3af3c10f522314db69cc7c8e73cbdf88e4
Upload date
Uploaded using Trusted Publishing?
What is trusted publishing?
No
Uploaded via twine/2.0.0 pkginfo/1.5.0.1 requests/2.22.0 setuptools/41.2.0 requests-toolbelt/0.9.1 tqdm/4.36.1 CPython/3.6.8

Release history Release notifications | RSS feed

This release

0.1.3 This release

2 release files

0.1.2

2 release files

0.1.1

2 release files

0.1.0

1 release file

Anthropic, PBC Visionary sponsor Bloomberg Visionary sponsor Hudson River Trading Visionary sponsor Meta Visionary sponsor NVIDIA Visionary sponsor Microsoft Sustainability sponsor Depot Continuous Integration AWS Cloud computing and Security Sponsor Datadog Monitoring Fastly CDN Google Download Analytics Sentry Error logging StatusPage Status page