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']])