is2D = False data = TXRMDataReader(file_name=filename).read() if is2D: data = data.get_slice(vertical='centre') else: binner = Binner(roi={ 'horizontal': (None, None, 4), 'vertical': (None, None, 4) }) data = binner(data) # show_geometry(data.geometry) data = TransmissionAbsorptionConverter()(data) data.reorder(order='tigre') ig = data.geometry.get_ImageGeometry() fbp = FBP(ig, data.geometry) recon = fbp(data) from cil.io import NEXUSDataWriter writer = NEXUSDataWriter() writer.set_up(data=recon, file_name=os.path.abspath('noce.nxs')) writer.write() # show2D([data, recon], fix_range=(-0.01,0.06)) from ccpi.viewer import viewer2D, viewer3D from ccpi.viewer.utils.conversion import Converter import vtk from functools import partial
#%% # neg log data = TransmissionAbsorptionConverter()(data) #%% from cil.utilities.jupyter import islicer islicer(data, direction='angle', cmap='inferno', origin='upper-right') #%% ig = data.geometry.get_ImageGeometry() #%% data.reorder(order='astra') #%% fbp = aFBP(ig, data.geometry) # fbp.set_input(data2d) # recon = fbp.get_output() recon = fbp(data) #%% from cil.utilities.jupyter import islicer islicer(recon, direction='vertical', cmap='inferno', origin='lower') #%% show2D([data, recon], title=['sinogram', 'ASTRA FBP'],\ slice_list=[['vertical','centre'],['vertical','410']])