Note
Click here to download the full example code
Example 03: Generate parallel-beam projectionsΒΆ
import foam_ct_phantom
import numpy as np
import h5py
import pylab as pl
pl.gray()
phantom = foam_ct_phantom.FoamPhantom('test_phantom.h5')
geom = foam_ct_phantom.ParallelGeometry(256,256,np.linspace(0,np.pi,128,False),3/256)
phantom.generate_projections('test_projs_par.h5',geom)
projs = foam_ct_phantom.load_projections('test_projs_par.h5')
pl.imshow(projs[0])
pl.show()