CIL-HTC2022-Algo5

Authors:

• Gemma Fardell (STFC), United Kingdom
• Jakob Sauer Jørgensen (DTU), Denmark
• Laura Murgatroyd (STFC), United Kingdom
• Evangelos Papoutsellis (STFC, Finden), United Kingdom
• Edoardo Pasca (STFC), United Kingdom

Addresses

STFC: United Kingdom Research and Innovation Scientific Computing Department of Science Technology Facilities Council Rutherford Appleton Laboratory Harwell Campus Didcot OX11 0QX

DTU: Technical University of Denmark, Department of Applied Mathematics and Computer Science Richard Petersens Plads Building 324 2800 Kgs. Lyngby Denmark

Finden: Building R71, Rutherford Appleton Laboratory, Harwell, Oxford, OX11 0QX

Description of the algorithm

This is an entry for the HTC2022 competition. The algorithm in algo.py is developed using CIL, a toolkit for tomographic imaging and optimisation. The main steps of the algorithm is:

1. Pre-processing: renormalisation, single material beam hardening correction, zero padding
2. Generation of pixelwise lower and upper bound circular masks (upper bound mask has a diameter of 97% of the image height)
3. Regularised iterative reconstruction algorithm using tools from CIL: Least-squares data fidelity with TV regularisation
4. Post-processing: segmentation of the reconstruction with multi-Otsu threshold

Installation instructions

Installation instructions, including any requirements.

conda env create --file environment.yml

Usage instructions.

conda activate htc-22-cil-algo5
python main.py path/to/input/files path/to/output/files 3

Examples

Examples of reconstructing the example datasets, where we have limited the angles to 90, 60 and 30 degree ranges. The 'ref' column is the given segmented result from the full dataset.

	Ref	90	60	30
ta
tb
tc
td

Scores for each sample and angle:

	90	60	30
ta	0.951	0.879	0.806
tb	0.888	0.823	0.737
tc	0.941	0.883	0.769
td	0.940	0.924	0.889

Repository content

• utils.py
• algo.py
• main.py
• environment.yml
• README.md
• recalc_score.py
• test_data
  • htc2022_ta_full_recon_fbp_seg.png
  • htc2022_ta_sparse_example.mat
• results
  • htc2022_ta_sparse_example.png
  • AR30
    • htc2022_ta_full.png
    • htc2022_tb_full.png
    • htc2022_tc_full.png
    • htc2022_td_full.png
  • AR60
    • htc2022_ta_full.png
    • htc2022_tb_full.png
    • htc2022_tc_full.png
    • htc2022_td_full.png
  • AR90
    • htc2022_ta_full.png
    • htc2022_tb_full.png
    • htc2022_tc_full.png
    • htc2022_td_full.png
  • AR360
    • htc2022_ta_full.png
    • htc2022_tb_full.png
    • htc2022_tc_full.png
    • htc2022_td_full.png

License

All files in the repository come with the Apache-v2.0 license unless differently specified.