FDK and ProjOp for RTK

[epapoutsellis]

Description

CIL Wrapper for RTK: The Reconstruction ToolKit

• FDK (3D-CPU)
• FDK (3D-GPU)
• ProjectionOperator (Cone)

@SimonRit, I would really appreciate your help for the RTK part and syntax.

Example Usage

Need to run bash scripts/create_local_env_for_cil_development.sh -t without

• ipp'>=2021.10'
• ipp-devel'>=2021.10'
• ipp-include'>=2021.10'
  and

pip install --no-deps .
python -m pip install itk-rtk

import os
from cil.io import ZEISSDataReader
from cil.plugins.rtk import FDK
from cil.processors import TransmissionAbsorptionConverter, Slicer, Binner

path = 'valnut/valnut_2014-03-21_643_28/tomo-A'
filename = os.path.join(path, "valnut_tomo-A.txrm")
data3D = ZEISSDataReader(file_name=filename).read()

binned_data = Slicer(roi={'angle':(0,1600,2)})(data3D)
binned_data = Binner(roi={'vertical':(None,None,2), 'horizontal':(None,None,2)}, accelerated=False)(binned_data)
absorption_data = TransmissionAbsorptionConverter()(binned_data) 

ag3D = absorption_data.geometry
ig3D = ag3D.get_ImageGeometry()

fdk_recon = FDK(ig3D, ag3D)(absorption_data)

See also CIL-RTK demo

Contribution Notes

• The content of this Pull Request (the Contribution) is intentionally submitted for inclusion in CIL (the Work) under the terms and conditions of the Apache-2.0 License
• I confirm that the contribution does not violate any intellectual property rights of third parties

❤️ Thanks for your contribution!

[lauramurgatroyd]

Hi @epapoutsellis, thanks for your PR. Can we convert this to a draft for now whilst you're still working on FDK (3D-GPU) and
ProjectionOperator (Cone)?