""" 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()