A Reconstruction Toolkit for Multichannel CT

Lead Research Organisation: University of Manchester
Department Name: Materials

Abstract

Currently, conventional Computed Tomographic (CT) imaging is still in a black and white (1 channel) era, just as it was with the first image Rontgen captured in 1895!

Conventional X-ray imaging entails a polychromatic X-ray source (i.e. with a full spectrum of energies) but with energy-indiscriminate detectors (registering a single grey-scale channel). However, technological breakthroughs in energy-sensitive detectors enable a new era of tomographic imaging in 'colours' (multiple channels). Each pixel of the energy-selective detector records a spectrum consisting of hundreds or thousands of energy channels. Currently available software only allows us to reconstruct each (noisy) channel independently in turn, which is a significant limitation. We need to unlock the power of next-generation correlative reconstruction methods for multi-channel tomography. Notably, the registered energy channels are mutually correlated, just like the red-green-blue (RGB) channels of the color image. Therefore, noise and other inaccuracies in spectral measurements can be treated holistically across the channels, leading to massive improvements in imaging quality (higher signal-to-noise ratio and resolution) in addition to fundamentally new opportunities such as spectroscopic imaging, i.e., direct decomposition into fundamental elements.

The overall goal of this CCP Software Flagship project is to expand upon existing single-channel image reconstruction software (already developed by the CCPi project) to enable sophisticated multi-channel correlative reconstruction methods. A novel Reconstruction Toolkit for Multichannel CT (RT-MCT) will be developed and become a part of the end-user data pipeline Savu (a modular Python-based platform for tomographic data processing developed at Diamond Light Source (DLS) at Harwell, UK). Three major imaging facilities are key collaborators and committed initial users of the RT-MCT: 1) Manchester X-ray Imaging Facility (MXIF) is a leader of laboratory-based X-ray CT imaging and has developed the unique multi-channel instrument "The Colour Bay" (cone-beam geometry scanner which uses HEXITEC hyper-spectral detectors); 2) A new national Neutron Imaging and Diffraction Facility (IMAT) at the ISIS pulsed neutron spallation source (Harwell). IMAT will take advantage of the neutron time-of-flight (TOF) measurement technique for effective energy discrimination into thousands of channels making this unique technique hyper-spectral; 3) Diamond Light Source (DLS), the national synchrotron facility at Harwell, has a number of imaging beamlines including I18 and I14, dedicated to X-ray fluorescence, X-ray spectroscopy and diffraction, all of which entail multi-channel data sets.

The main aim is to deliver the RT-MCT to these facilities to provide much more efficient data reconstruction and analysis. Several work packages are identified which constitute the RT-MCT, namely a) accurate mathematical modelling of multi-channel imaging; b) formulation of optimal reconstruction problems; c) efficient algorithm implementation and integration in existing software framework; d) deployment to facilities and use in proof-of-concept case studies.

Planned Impact

Societal: This research will focus on improving spectral imaging through reconstruction software development which will revolutionise materials science analysis. The proposed project will ensure closer collaborations with other CCPs; especially the CCP PET/MR (hybrid imaging) and CCP-EM (electron microscopy) groups. At present X-ray Computed Tomography is used for innovative non-destructive inspection, analysis of new materials, and threat detection in airport security screening. Multichannel imaging will provide a sea-change in research capability providing faster and more accurate image analysis with chemical mapping overlay generated in parallel. The full potential of multispectral imaging can only be realised through the development of a software toolkit for end users, the RT-MCT. With accurate image analysis provided by the RT-MCT industrial connections have already led to new funded projects with coders from the University of Manchester (MXIF) accessing raw and metadata outputs and reconstruction expertise, for example writing code for Nikon X-ray scanner. Industrial impact will be enhanced through continual collaborations with members of the CCP network and Working Group: in collaboration with industrial network members. The following industrial companies and national organisations are active members of the CCP network and on the Working Group; Nikon Metrology, FEI, LaVision, Zeiss, Bruker, Deben, Rapiscan, Simpleware, Tessella and the BSI Group.

Economic: Industrial X-ray Inspection is worth in excess of $450M worldwide and emerging applications are driving new demand for CT inspection. We will maximise the impact to industry by: 1) extending our current collaborations with the equipment manufacturers Nikon Metrology, FEI, LaVision, Zeiss, Bruker, Deben, Rapiscan, Simpleware, Tessella and the BSI Group. We will focus on developing both the software, ensuring we are at the forefront of new technology. 2) The University will utilise its network of commercial partners as pathways to industrial impact, specifically in the oil & gas, nuclear and aerospace sectors through dedicated research institutes: BP International Centre for Advanced Materials (BP-ICAM), Dalton Nuclear Institute, Aerospace Research Institute. 3) Collaborate with the SME Spectral-X, established in early 2014 to commercialise colour CT imaging, to capture value from truly "game-changing'' developments.

Scientific: The outputs from this Flagship will assist some of the major national facilities (ISIS and DLS) and lab-based systems ensuring that the UK is at the forefront of developments in this new analytical capability. It will also utilise related infrastructures including compute clusters within STFC and in the regions (N8, Archer and Hartree Centre). Diamond Light Source (DLS) scientists and the hundreds of academic and industrial beamline users annually will greatly benefit from the RT-MCT methods and software applied to their experimental data produced on energy-dispersive X-ray spectroscopy and diffraction (I18, I14) beamlines. IMAT ISIS neutron facility will accommodate the Savu-packaged RT-MCT for all future users. MXIF where laboratory based colour imaging was pioneered and the RT-MCT will substantially reduce the amount of time required to process data, increasing the research capability of the instrument.

People: The retention and development skilled people able to work in multidisciplinary teams is not only essential to the future of the CCPi but also to develop the leaders of tomorrow. We will provide the postdoctoral researchers with opportunities to acquire and improve upon a wide range of technical and transferable skills, which will be of value both during this project, and in their subsequent careers. The PDRAs will benefit from working in a multidisciplinary research environment and we place a high priority on training our researchers in communication and engagement with end users of our research.

Publications

10 25 50
 
Description Until recently, X-ray CT has been undertaken in black and white, providing simply a black and white 3-D image of the attenuation through the object. Hyper-spectral or colour imaging records a full profile at each pixel. In most cases, these were reconstructed individually, which because of the poor signal gave many low quality 3-D images. In this work we have developed reconstruction codes capable of reconstructing colour 3-D images.

We have expanded the Core Imaging Library to handle different tasks such as such as de-noising, de-blurring, in-painting (i.e. filling in blank spaces) and tomography reconstruction. In addition, tools have been developed to process extremely large amounts of data received from a scanner; they speed up reconstruction of the final image by slicing the data into very small chunks. There are also pre-processing tools capable of preventing unwanted artefacts appearing in the final image.
Exploitation Route Currently the software has been installed on the IMAT beamline at the ISIS national Neutron facility for colour time of flight imaging and will be used to enable colour X-ray CT imaging within the National Research Facility for lab X-ray CT, which became operational in November 2020.

The development team have also enabled certain features from the CIL to be incorporated in the Synergistic Image Reconstruction Framework (SIRF) from the Collaborative Computational Project in Synergistic Reconstruction for Biomedical Imaging (CCP SyneRBI), suitable for PET and MRI reconstruction.

The work funded by this grant has been further developed and exploited by a number of other researchers. Dr. Evangelos Papoutsellis was a finalist in the CoSeC Impact Award in 2020 for his work on the Core Imaging Library, while Ryan Warr, a PhD student at the University of Manchester, won the award the following year for his work developing Spectral X-ray Computed Tomography imaging techniques.

Dr Tim Burnett has been awarded funding with ARTC and SIMTech in Singapore to investigate the potential of machine learning to exploit the use of hyperspectral X-ray CT.

Finally, we are in the process of preparing an ERC advanced grant to exploit colour imaging.
Sectors Aerospace

Defence and Marine

Energy

Manufacturing

including Industrial Biotechology

Pharmaceuticals and Medical Biotechnology

Transport

URL http://www.ccpi.ac.uk/cil
 
Description This has contributed to a £2m innovate/industry project to develop colour imaging for medical applications. It has also led to software that is now being used by IMAT neutron beamline at ISIS to enable their users to interpret Bragg edge signals. Furthermore, it is going to drive software for the new colour x-ray instrument which is being established with NXCT for the UK user community. This has led to a collaboration with Adaptix (SME) resulting in a £1.6m innovate grant application to develop colour tomosynthesis.
First Year Of Impact 2017
Sector Aerospace, Defence and Marine,Digital/Communication/Information Technologies (including Software),Energy,Healthcare,Manufacturing, including Industrial Biotechology,Culture, Heritage, Museums and Collections,Pharmaceuticals and Medical Biotechnology
Impact Types Economic

 
Description EPSRC XCT Equipment Roadmap Panel
Geographic Reach National 
Policy Influence Type Membership of a guideline committee
 
Description Reader for BSI/ISO standard XCT TDF/4/4/1
Geographic Reach National 
Policy Influence Type Membership of a guideline committee
 
Description Alan Turing Institute: Summer School Application (Large Data Image Reconstruction and Visualisation)
Amount £4,500 (GBP)
Organisation Alan Turing Institute 
Sector Academic/University
Country United Kingdom
Start 05/2020 
End 06/2020
 
Description EPSRC: Multiscale and In Situ Laboratory X-ray Computed Tomography National Research Facility
Amount £10,000,000 (GBP)
Organisation Engineering and Physical Sciences Research Council (EPSRC) 
Sector Public
Country United Kingdom
Start 07/2020 
End 07/2025
 
Description EPSRC: Tomographic Imaging: UK Collaborative Computational Projects
Amount £296,693 (GBP)
Funding ID EP/T026677/1 
Organisation Engineering and Physical Sciences Research Council (EPSRC) 
Sector Public
Country United Kingdom
Start 08/2020 
End 03/2025
 
Description Rich Nonlinear Tomography for advanced materials (LEAD)
Amount £635,422 (GBP)
Funding ID EP/V007742/1 
Organisation Engineering and Physical Sciences Research Council (EPSRC) 
Sector Public
Country United Kingdom
Start 05/2021 
End 05/2024
 
Title Hyperspectral X-ray CT Voxelized TV reconstruction of a single, iodine-stained lizard head sample 
Description Dataset description These datasets are voxel based reconstructions of hyperspectral CT data using the Core Imaging Library (CIL). They are stored as NeXus files (derived from hdf5) which can be read in, visualised and manipulated using CIL. - PDHG_TV_1000_Sp_alpha_0.004.nxs Is the solution after 1000 iterations of PDHG with TV applied in the spatial domain. - PDHG_TV_1000_SpCh_alpha_0.003_beta_0.5.nxs Is the solution after 1000 iterations of PDHG with TV applied both in the spatial domain, and in the energy (channel) domain. Dataset intended use These datasets are used in the CIL training notebook: https://github.com/TomographicImaging/CIL-Demos/blob/main/examples/3_Multichannel/03_Hyperspectral_reconstruction.ipynb They can be imported using CIL, with the following code snippet:
from cil.io import NEXUSDataReader reader = NEXUSDataReader(file_name='path/to/data/PDHG_TV_1000_Sp_alpha_0.004.nxs') data = reader.read()
 
Type Of Material Database/Collection of data 
Year Produced 2019 
Provided To Others? Yes  
URL https://zenodo.org/record/7016573
 
Title Hyperspectral X-ray CT data set of mineralised ore sample with Au and Pb deposits 
Description General data description: This is a hyperspectral (energy-resolved) X-ray CT projection data set of a mineralised ore sample with small gold and galena deposits. It was acquired in a laboratory micro-CT scanner with an energy-sensitive HEXITEC detector in the Henry Moseley X-ray Imaging Facility at The University of Manchester. The data included contains all the relevant files required for reconstruction, following a hyperspectral scan of a mineralised ore sample. The sample contains a number of mineral phases, of varying concentration, distributed throughout. Some phases (including gold, and lead-based Galena) produce unique absorption edges, which act as spectral identifiers that can be measured by an energy-sensitive detector. File descriptions: The data set consists of one .txt file and three .mat (MATLAB) data files. Au_rock_scan_geometry.txt gives a breakdown of the full sample and detector geometry used when acquiring the raw projections. The number of horizontal detector pixels accounts for the fact that a set of 5 tiled scans of the sample were collected and later stitched together. Au_rock_sinogram_full.mat contains the full 4D sinogram constructed following flat-field normalisation of the raw projection data. The data matrix contains the total number of energy channels acquired during scanning, as well as the conventional elements of vertical/horizontal detector pixel number and total projection angles. commonX.mat provides a direct conversion between the energy channels, and the energies (in keV) that they correspond to, following a calibration procedure prior to scanning. FF.mat contains the 4D flatfield data acquired when no sample was present. This data was used to normalise the projection datasets, as the sinogram was constructed. 
Type Of Material Database/Collection of data 
Year Produced 2020 
Provided To Others? Yes  
URL https://zenodo.org/record/4157614
 
Title Hyperspectral X-ray CT dataset of a single, iodine-stained lizard head sample 
Description General Data description: This is a hyperspectral (energy-resolved) X-ray CT projection dataset of a lizard head sample, stained with an iodine contrast agent. It was acquired in a custom-built, laboratory micro-CT scanner with an energy-sensitive HEXITEC detector in the Henry Moseley X-ray Imaging Facility at The University of Manchester. The following data contains all the files necessary for reconstruction, after a hyperspectral scan was taken of a single, iodine-stained lizard head sample. The iodine contrast agent provided a spectral marker, measured by an energy-sensitive detector, which may be used for spatial mapping and segmentation of stained soft tissue regions. File descriptions: Contained are four MATLAB (.mat) data files, as well as a single text (.txt) file. Lizard_head_scan_parameters.txt provides the full sample and detector geometry of the scan acquisition. lizard_180Proj_noSupp_1_180.mat contains the full 4D sinogram constructed following flatfield normalisation of the raw projection data. The 4D array contains the total number of energy channels acquired during scanning, vertical and horizontal pixel number, and total projections angles acquired. The data provided is prior to application of any post-processing filters. The first 180 energy channels are included. lizard_180Proj_Supp_1_180.mat contains the full 4D sinogram constructed following flatfield normalisation of the raw projection data. This dataset is identical to the .mat file above, however here we have also applied a ring-reduction filter, using a wavelet-based Fourier filter which suppresses the presence of ring artefacts in every energy channel. The first 180 energy channels are included. Energy_axis.mat provides a direct conversion between the energy channels, and the energies (in keV) that they correspond to, following a calibration procedure prior to scanning. FF.mat contains the 4D flatfield data acquired when no sample was present. This data was used to normalise the projection datasets, as the sinogram was constructed. The first 180 energy channels are included. 
Type Of Material Database/Collection of data 
Year Produced 2022 
Provided To Others? Yes  
Impact Not yet realised 
URL https://zenodo.org/record/4352943
 
Title Hyperspectral X-ray CT datasets for a set of multiply-stained mouse limb specimens 
Description General Data description: The following are hyperspectral (energy-resolved) X-ray CT datasets for a set of mouse limb specimens, each stained with multiple contrast agents. All scans were acquired with an energy-sensitive HEXITEC detector in the Henry Moseley X-ray Imaging Facility at The University of Manchester. The following data contains all the files necessary for reconstruction of each dataset. The biological specimens were produced as they each contain multiple contrast agents, with distinct spectral markers. When measured by an energy-sensitive detector, each contrast agent may be identified and segmented individually following spectral analysis. A mouse hindlimb was double-stained with elemental iodine and BaSO4. A mouse forelimb was triple-stained with I2KI, BaSO4 and PTA. File descriptions: Contained are two HDF5 (.h5) data files, as well as two (.txt) metadata files and a MATLAB (.mat) file. Hindlimb_scan_parameters.txt provides the full sample and detector geometry of the scan acquisition for the double-stained hindlimb. Forelimb_scan_parameters.txt provides the full sample and detector geometry of the scan acquisition for the triple-stained forelimb. DS_Mouse_hindlimb_sinogram.h5 contains the full 4D sinogram constructed following flatfield normalisation of the raw projection data for the double-stained hindlimb specimen. The 4D array contains the total number of energy channels acquired during scanning, followed by vertical and horizontal pixel number, and finally total projections angles acquired. In addition, a ring artefact reduction filter was applied. TS_Mouse_forelimb_sinogram.h5 contains the full 4D sinogram constructed following flatfield normalisation of the raw projection data for the triple-stained forelimb specimen. The 4D array contains the total number of energy channels acquired during scanning, followed by vertical and horizontal pixel number, and finally total projections angles acquired. In addition, a ring artefact reduction filter was applied. Energy_axis.mat provides a direct conversion between the energy channels, and the energies (in keV) that they correspond to, following a calibration procedure prior to scanning. 
Type Of Material Database/Collection of data 
Year Produced 2022 
Provided To Others? Yes  
URL https://zenodo.org/record/6787593
 
Title Hyperspectral X-ray CT datasets of an aluminium phantom containing three metal-based powders 
Description General Data description: This is a set of two hyperspectral (energy-resolved) X-ray CT projection datasets of a multi-phase phantom. It was acquired in a custom-built, laboratory micro-CT scanner with an energy-sensitive HEXITEC detector in the Henry Moseley X-ray Imaging Facility at The University of Manchester. The following data contains all the files necessary for reconstruction, following two hyperspectral scans of a metal, multi-phase phantom. The phantom consists of an external aluminium cylinder, with three holes, each filled with a different metal-based powder (CeO2, ZnO, Fe). Each powder provides a unique attenuation signal, with CeO2 in particular producing a distinct spectral marker which can be measured by an energy-sensitive detector. Two identical scans were acquired, with only the exposure time per projection changed. Note: Zenodo Version 2 of this dataset contains the incorrect version of the 180s, 180 projection phantom dataset, if wishing to analyse the dataset used in the associated hyperspectral paper. This version (Version 3) contains the correct dataset from the paper. File descriptions: Contained is an image (.jpg) of the sample, along with five MATLAB (.mat) data files, as well as a single text (.txt) file. Where necessary, the files have been named to match the dataset they belong to, based on the different exposure times used for each dataset. Phantom_design_measurements.jpg contains a photograph of the physical phantom, combined with a diagram showing full sample measurements. Powder_phantom_scan_geometry.txt gives a breakdown of the full sample and detector geometry used when acquiring the raw projections for both scans. Powder_phantom_30s_30Proj_sinogram.mat contains the 4D sinogram constructed following flatfield normalisation of the raw projection data, where an exposure time of 30 s was used for each projection. The 4D array contains the total number of energy channels acquired during scanning, followed by vertical and horizontal pixel number, and finally total projections angles acquired during scanning. The total number of channels in the file is 200. Powder_phantom_180s_180Proj_sinogram.mat is the 4D sinogram for the dataset, when exposure times of 180 s were used for each projection, following flatfield normalisation. A discontinuity occurs at projection 137 due to an interruption in the scan procedure. The total number of channels in the file is 200. Energy_axis.mat provides a direct conversion between the energy channels, and the energies (in keV) that they correspond to, following a calibration procedure prior to scanning. This is the same for both datasets. FF_30s.mat contains the 4D flatfield data acquired when no sample was present, in the case of 30 s exposure times. This data was used to normalise the projection datasets, as the sinogram was constructed. The first 200 channels are included. FF_180s.mat contains the 4D flatfield data for the dataset where 180 s exposure times were used. The first 200 channels are included. 
Type Of Material Database/Collection of data 
Year Produced 2022 
Provided To Others? Yes  
Impact N/A 
URL https://zenodo.org/record/4354815
 
Title Hyperspectral X-ray CT datasets of three chemical phantoms 
Description General Data description: The following are hyperspectral (energy-resolved) X-ray CT datasets for a set of chemical phantom samples, each containing multiple phases of an aqueous contrast agent at different concentrations. All scans were acquired with an energy-sensitive HEXITEC detector in the Henry Moseley X-ray Imaging Facility at The University of Manchester. The following data contains all the files necessary for reconstruction of each dataset. The phantom samples were produced as they each offer a distinct spectral marker which, when measured by an energy-sensitive detector, may be used as a form of calibration for spectral analysis. The phantoms were for the common contrast agents of I2KI, BaSO4 and PTA. File descriptions: Contained are four MATLAB (.mat) data files, as well as three text (.txt) metadata files. Iodine_Phantom_scan_parameters.txt provides the full sample and detector geometry of the scan acquisition for the I2KI phantom. The concentrations for the iodine phases were 25, 50, 76 and 101 mg/ml of aqueous I3- ions respectively. Barium_Phantom_scan_parameters.txt provides the full sample and detector geometry of the scan acquisition for the BaSO4 phantom. The concentrations for the BaSO4 phases were 100, 200 and 400 mg/ml of BaSO4 respectively. Tungsten_Phantom_scan_parameters.txt provides the full sample and detector geometry of the scan acquisition for the PTA phantom. The concentrations for the PTA phases were 50, 100 and 200 mg/ml of PTA respectively. Iodine_phantom_sinogram.mat contains the full 4D sinogram constructed following flatfield normalisation of the raw projection data for the I2KI phantom. The 4D array contains the total number of energy channels acquired during scanning, followed by vertical and horizontal pixel number, and finally total projections angles acquired. In addition, a ring artefact reduction filter was applied, as well as a centre-of-rotation correction. Barium_phantom_sinogram.mat contains the full 4D sinogram constructed following flatfield normalisation of the raw projection data for the BaSO4 phantom. The 4D array contains the total number of energy channels acquired during scanning, followed by vertical and horizontal pixel number, and finally total projections angles acquired. In addition, a ring artefact reduction filter was applied, as well as a centre-of-rotation correction. Tungsten_phantom_sinogram.mat contains the full 4D sinogram constructed following flatfield normalisation of the raw projection data for the PTA phantom. The 4D array contains the total number of energy channels acquired during scanning, followed by vertical and horizontal pixel number, and finally total projections angles acquired. In addition, a ring artefact reduction filter was applied, as well as a centre-of-rotation correction. Energy_axis.mat provides a direct conversion between the energy channels, and the energies (in keV) that they correspond to, following a calibration procedure prior to scanning. 
Type Of Material Database/Collection of data 
Year Produced 2022 
Provided To Others? Yes  
URL https://zenodo.org/record/6787488
 
Title Neutron tomography data of high-purity metal rods using golden-ratio angular acquisition (IMAT, ISIS) 
Description General description: This is a neutron tomography projection data set acquired at the IMAT beamline of the ISIS Neutron and Muon Source, Harwell, UK. It was acquired using a 512-by-512 pixel MCP time-of-flight detector (pixel size 0.055 mm) binned into 2332 energy intervals between 1.7Å and 5.5Å using a parallel-beam scan geometry. 186 projections, each with 30 min exposure time (20 micro A), in the golden-ratio angular acquisition-mode. That is, a constant angular increment of (sqrt(5)-1)/2*180 degrees = 111.2461 degrees was used, and if an angle was outside 0 to 180 it was wrapped to that interval. Data was preprocessed to compensate for distortions caused by detector electronics and then summed over energy bins to simulate a white-beam data set. The data has not been corrected for centre-of-rotation offset. This can be achieved for example by cropping 52 pixels off all images from the left. The sample was an Al cylinder (Ø 22mm) with cylindrical inset holes drilled from the top: A central 10mm hole (left empty) and 5 holes of 3mm in pentagon corner positions and 5 holes of 1mm in pentagon corner positions offset from the 3mm set. High-purity elemental Cu, Fe, Ni, Ti, and Zn were inserted in the 3mm and 1mm holes with a single 1mm hole left empty. File descriptions: The data consists of a single ZIP file with the following contents: 186 TIFF files - these are the individual projection images. golden_ratio_angles.txt - this text file lists the 186 angles in degrees at which the projections were acquired. scan_parameters.txt - this text file contains geometric metadata about the experimental scan configuration. 
Type Of Material Database/Collection of data 
Year Produced 2021 
Provided To Others? Yes  
Impact Impact not yet realised 
URL https://zenodo.org/record/4273968
 
Title Neutron tomography data of high-purity metal rods using golden-ratio angular acquisition (IMAT, ISIS) 
Description This is a neutron tomography projection data set acquired at the IMAT beamline of the ISIS Neutron and Muon Source, Harwell, UK. 
Type Of Material Database/Collection of data 
Year Produced 2021 
Provided To Others? Yes  
Impact N/A 
URL https://zenodo.org/record/4273969#.ZA8jsnbP0uV
 
Title SIRF: Synergistic Image Reconstruction Framework 
Description The combination of positron emission tomography (PET) with magnetic resonance (MR) imaging opens the way to more accurate diagnosis and improved patient management. At present, the data acquired by PET-MR scanners are essentially processed separately, but the opportunity to improve accuracy of the tomographic reconstruction via synergy of the two imaging techniques is an active area of research. In this paper, we present Release 2.1.0 of the CCP-PETMR Synergistic Image Reconstruction Framework (SIRF) software suite, providing an open-source software platform for efficient implementation and validation of novel reconstruction algorithms. SIRF provides user-friendly Python and MATLAB interfaces built on top of C++ libraries. SIRF uses advanced PET and MR reconstruction software packages and tools. Currently, for PET this is Software for Tomographic Image Reconstruction (STIR); for MR, Gadgetron and ISMRMRD; and for image registration tools, NiftyReg. The software aims to be capable of reconstructing images from acquired scanner data, whilst being simple enough to be used for educational purposes. The most recent version of the software can be downloaded from http://www.ccppetmr.ac.uk/downloads and https://github.com/CCPPETMR/. 
Type Of Material Database/Collection of data 
Year Produced 2020 
Provided To Others? Yes  
URL https://data.mendeley.com/datasets/s45f5jh55j
 
Title SIRF: Synergistic Image Reconstruction Framework 
Description The combination of positron emission tomography (PET) with magnetic resonance (MR) imaging opens the way to more accurate diagnosis and improved patient management. At present, the data acquired by PET-MR scanners are essentially processed separately, but the opportunity to improve accuracy of the tomographic reconstruction via synergy of the two imaging techniques is an active area of research. In this paper, we present Release 2.1.0 of the CCP-PETMR Synergistic Image Reconstruction Framework (SIRF) software suite, providing an open-source software platform for efficient implementation and validation of novel reconstruction algorithms. SIRF provides user-friendly Python and MATLAB interfaces built on top of C++ libraries. SIRF uses advanced PET and MR reconstruction software packages and tools. Currently, for PET this is Software for Tomographic Image Reconstruction (STIR); for MR, Gadgetron and ISMRMRD; and for image registration tools, NiftyReg. The software aims to be capable of reconstructing images from acquired scanner data, whilst being simple enough to be used for educational purposes. The most recent version of the software can be downloaded from http://www.ccppetmr.ac.uk/downloads and https://github.com/CCPPETMR/. 
Type Of Material Database/Collection of data 
Year Produced 2020 
Provided To Others? Yes  
URL https://data.mendeley.com/datasets/s45f5jh55j/1
 
Title Anaconda of the Python software CIL 
Description Good software development practices have been implemented, including software code project management, version control, issue tracking, and systematic code testing and builds. 
Type Of Technology Software 
Year Produced 2018 
Impact We have made public releases through Anaconda of the Python software CIL. 
URL http://cil.readthedocs.io/en/latest/
 
Title Code to reproduce results of "Core Imaging Library Part I: a versatile python framework for tomographic imaging" 
Description This code reproduces all the results presented in the article Core Imaging Library Part I: a versatile python framework for tomographic imaging by Jakob S. Jørgensen, Evelina Ametova, Genoveva Burca, Gemma Fardell, Evangelos Papoutsellis, Edoardo Pasca, Kris Thielemans, Martin Turner, Ryan Warr, William R. B. Lionheart, and Philip J. Withers which will be available from 5 July 2021 at https://doi.org/10.1098/rsta.2020.0192 A preprint is available from arXiv: https://arxiv.org/abs/2102.04560 Instructions are available in the file README.md as well as at the source GitHub repository https://github.com/TomographicImaging/Paper-2021-RSTA-CIL-Part-I 
Type Of Technology Software 
Year Produced 2021 
Impact This has grown to be used in many institutions across the UK, Europe and the US, many collaborations formed including with the ISIS neutron facility, EPAC laser facility. 
URL https://zenodo.org/record/4744394
 
Title Code to reproduce results of "Core Imaging Library Part II: multichannel reconstruction for dynamic and spectral tomography" 
Description This code reproduces all the results presented in the article Core Imaging Library Part II: multichannel reconstruction for dynamic and spectral tomography by Evangelos Papoutsellis, Evelina Ametova, Claire Delplancke, Gemma Fardell, Jakob S. Jørgensen, Edoardo Pasca, Martin Turner, Ryan Warr, William R. B. Lionheart, and Philip J. Withers which will be available from 5 July 2021 at https://doi.org/10.1098/rsta.2020.0193 A preprint is available from arXiv: https://arxiv.org/abs/2102.06126 Instructions are available in the file README.md as well as at the source GitHub repository https://github.com/TomographicImaging/Paper-2021-RSTA-CIL-Part-II 
Type Of Technology Software 
Year Produced 2021 
URL https://zenodo.org/record/4744745
 
Title The Core Imaging Library 
Description CIL is an open-source mainly Python framework for tomographic imaging for cone and parallel beam geometries. It comes with tools for loading, preprocessing, reconstructing and visualising tomographic data. 
Type Of Technology Webtool/Application 
Year Produced 2021 
Open Source License? Yes  
Impact This has grown to be used in many institutions across the UK, Europe and the US, many collaborations formed including with the ISIS neutron facility, EPAC laser facility. 
 
Description 12th annual Tomography for Scientific Advancement (ToScA) symposium 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact Team member gave presentation at virtual conference
Year(s) Of Engagement Activity 2021
 
Description Analysis Workshop 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Professional Practitioners
Results and Impact Workshop on digital volume correlation (DVC) analysis, outlining state of the art and applications in materials science.
Year(s) Of Engagement Activity 2021
 
Description CCPi Virtual Seminar Series 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Postgraduate students
Results and Impact Seminar series over zoom in which researchers and industry discussed how they used the Core Imaging Library. There were over 25 events covering a total of 369 attendees with 45 speakers.
Year(s) Of Engagement Activity 2020,2021
 
Description CCPi Working Group event 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Professional Practitioners
Results and Impact Discussion with SPIERS looking software product range
Year(s) Of Engagement Activity 2021
 
Description CIL @ ToScA UK & Europe 2022, Queen Mary University London, 7 - 9 September 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact We present the Core Imaging Library (CIL) - a versatile open-source Python package for tomographic imaging reconstruction and analysis. CIL can be used for the entire tomography data workflow from loading common lab and synchrotron X-ray CT data formats, pre-processing for various imperfections, reconstruction by filtered back-projection or a variety of iterative methods as well as data visualisation and exploration. CIL offers tools for conventional as well as challenging setups including in-situ experiments with imaging rigs, fast dynamic CT, hyper-spectral CT, and non-standard scan geometries such as laminography. In this talk, we give an overview of the functionality provided by CIL and provide multiple application examples, including comparisons of filtered back-projection reconstruction with emerging iterative reconstruction methods.
Year(s) Of Engagement Activity 2022
URL https://www.ccpi.ac.uk/CIL
 
Description CIL Training @ ToScA UK 2021 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Professional Practitioners
Results and Impact Hands-on training for the Core Imaging Library (CIL), an open-source reconstruction platform for challenging and novel data.
Year(s) Of Engagement Activity 2021
 
Description CIL Training Course @ IBSim 18 October 2022 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Professional Practitioners
Results and Impact Session 1: Introduction to XCT and the CIL framework for reconstruction of commercial lab-based XCT data with FBP/FDK
Session 2: Pre-processing of XCT data, e.g. centre of rotation correction
Session 3: Iterative reconstruction methods for standard XCT data
Session 4 (choose one of):
Reconstruction of data simulated from gVXR
Advanced iterative reconstruction with spectral XCT data
Advanced iterative reconstruction on laminography XCT data
Year(s) Of Engagement Activity 2022
URL https://ibsim.co.uk/events/ibsim-4i/registration/
 
Description CT image based modelling of composite materials' failure 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact An invited lecture to an international conference.
Year(s) Of Engagement Activity 2022
 
Description Dataset training 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Professional Practitioners
Results and Impact A training course was held over two days at Harwell Campus to help users with the Visualisation & quantification of tomographic datasets
Year(s) Of Engagement Activity 2021
 
Description Digital Volume Correlation. Training on iDVC 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Professional Practitioners
Results and Impact Digital volume correlation (DVC) is used in material science to quantify the internal displacements and strain fields from in situ tomography experiments. This workshop will cover the theory of in situ tomography experiments and DVC analysis, followed by a practical session using the Collaborative Computational Project in Tomographic Imaging (CCPi) iDVC app.
Year(s) Of Engagement Activity 2022
 
Description Expanding 3D Nondestructive X-ray Microscopy Through Laboratory Diffraction Contrast Tomography (LabDCT) 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact An introductory webinar to LabDCT is and how it works. It was aimed at materials scientists, engineers, and researchers working in either the academic or industrial environments interested in understanding the microstructural and crystallographic information of (single/poly)-crystalline samples.
Year(s) Of Engagement Activity 2020
 
Description Fully3D Training School for the Synergistic Image Reconstruction Framework (SIRF) and Core Imaging Library (CIL) 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Professional Practitioners
Results and Impact This was an online training course for the image reconstruction, optimisation and regularisation software packages: Synergistic Image Reconstruction Framework (SIRF) and the Core Imaging Library (CIL).
Year(s) Of Engagement Activity 2021
 
Description Hackathon 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Professional Practitioners
Results and Impact A hackathon took place at STFC to develop an algorithms benchmark for tomographic image reconstruction
Year(s) Of Engagement Activity 2021
 
Description IBSim-4i (Image-Based Simulation for Industry) 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact Team member gave presentation at virtual conference
Year(s) Of Engagement Activity 2020
URL https://ibfem.co.uk/events/
 
Description Lunch-and-learn virtual seminar series 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact A series of weekly virtual seminars with invited speakers.
Year(s) Of Engagement Activity 2020,2021
 
Description Network t-conf meetings 
Form Of Engagement Activity A formal working group, expert panel or dialogue
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Other audiences
Results and Impact Network t-conf meetings occur continually, with updated website, and regular newsletters, comprises talks and training events. The size of our community has risen from ~250 in 2013 to over 400 by the end of 2018; over 60% growth.
Year(s) Of Engagement Activity 2013,2014,2015,2016,2017,2018
 
Description Online Training - Introduction to the Core Imaging Library (CIL) 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact This short training course will introduce you to using CIL for pre-processing and reconstruction of tomographic data. You will be introduced to iterative reconstruction algorithms and compare these with the typical FBP/FDK reconstructions.
The course will be held remotely. The training material will be provided for you through the STFC cloud which you will access through a browser.
Year(s) Of Engagement Activity 2022
 
Description Online Training - Introduction to the Core Imaging Library (CIL) 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact This short training course will introduce you to using CIL for pre-processing and reconstruction of tomographic data. You will be introduced to iterative reconstruction algorithms and compare these with the typical FBP/FDK reconstructions.
The course will be held remotely. The training material will be provided for you through the STFC cloud which you will access through a browser.
Year(s) Of Engagement Activity 2021
 
Description Poster Presentation at STEM for Britain Event (Westminster) 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Policymakers/politicians
Results and Impact STEM for BRITAIN exists to raise the profile of Britain's early-stage researchers at Westminster by engaging Members of both Houses of Parliament with current science, engineering and mathematics research being undertaken in the UK, especially that by their local constituents and in their local University.
Year(s) Of Engagement Activity 2020
 
Description Software Training on CCPi Core Imaging Library @ Synergistic Symposium, Chester UK, Nov 2019 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Other audiences
Results and Impact 45 international researchers from postgraduate to professor level attended a 2-day software training school in a software output from the project, which received stellar feedback and increased uptake of software in the imaging community.
Year(s) Of Engagement Activity 2019
 
Description Talk at IBSim-4i 2021 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact Gave a talk on the CIL at IBSim-4i 2021 to promote the resource to a wider audience
Year(s) Of Engagement Activity 2021
 
Description Talk at the CoSeC Annual Conference 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Professional Practitioners
Results and Impact A talk was given at the CoSeC annual conference entitled "The CCPi core imaging library, a versatile software for tomographic imaging", which helped to raise the profile of the resource.
Year(s) Of Engagement Activity 2021
URL https://www.scd.stfc.ac.uk/Pages/CIUK-2021-Breakout-Sessions.aspx
 
Description Training Course 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Postgraduate students
Results and Impact Avizo Course at University of Manchester designed to guide beginners on how to use the programme.
Year(s) Of Engagement Activity 2021
 
Description Training Course at Conference 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact Three week training course prior to the 16th International Meeting on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine. The course covered image reconstruction for PET, MRI and x-ray CT, emphasising commonalities but allowing participants to concentrate on topics that fit their interest.
Year(s) Of Engagement Activity 2021
 
Description Training at conference 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Postgraduate students
Results and Impact Training day on reconstruction toolbox which increased understanding of software.
Year(s) Of Engagement Activity 2021
URL https://www.npl.co.uk/research/dimensional/dxct-conference/dxct-conference-2021
 
Description Training school at PSMR2022 - Total body PET workshop 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact Workshop and training session for CIL
Year(s) Of Engagement Activity 2022
 
Description Training session at conference 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact Provided two 3hr hands-on CIL training sessions at the 2021 ToScA Europe conference, which led to an increased understanding of how to use the software
Year(s) Of Engagement Activity 2021
URL https://www.toscainternational.org