Note
Click here to download the full example code
Example 08: Apply Poission noise to projection dataΒΆ
import foam_ct_phantom
import numpy as np
import h5py
import pylab as pl
pl.gray()
# Run 03_generate_projections_parallel.py first
fac = foam_ct_phantom.estimate_absorption_factor('test_projs_par.h5',0.5)
foam_ct_phantom.apply_poisson_noise('test_projs_par.h5','test_projs_noisy.h5',1234,100,fac)
projs = foam_ct_phantom.load_projections('test_projs_noisy.h5')
pl.imshow(projs[0])
pl.show()