"""
FBP reconstruction example for simulated Skull CT data
"""
import numpy as np
import adutils
# Discretization
reco_space = adutils.get_discretization()
# Forward operator (in the form of a broadcast operator)
A = adutils.get_ray_trafo(reco_space, use_subset=True)
# Define fbp
fbp = adutils.get_fbp(A)
# Data
rhs = adutils.get_data(A, use_subset=True)
# Reconstruct
x = fbp(rhs)
# Show result
x.show(coords=[None, 0, None])
x.show(coords=[0, None, None])
x.show(coords=[None, None, 90])
# Save
saveReco = False
if saveReco:
saveName = '/home/user/Simulated/120kV/reco/Reco_HelicalSkullCT_70100644Phantom_no_bed_Dose150mGy_FBP.npy'
FBP reconstruction example for simulated Skull CT data - 2D Note: So far only works for phantom_number = '70100644' """ import numpy as np import adutils # Discretization reco_space = adutils.get_discretization(use_2D=True) # Forward operator (in the form of a broadcast operator) A = adutils.get_ray_trafo(reco_space, use_2D=True) # Define fbp fbp = adutils.get_fbp(A, use_2D=True) # Data rhs = adutils.get_data(A, use_2D=True) # Reconstruct x = fbp(rhs) # Show result x.show() # Compare to phantom phantom = reco_space.element(adutils.get_phantom(use_2D=True)) phantom.show()