Example 05: Generate a 4D expanding phantom

import foam_ct_phantom
import h5py
import numpy as np
import pylab as pl
pl.gray()

random_seed = 9876

foam_ct_phantom.ExpandingFoamPhantom.generate('test_expand.h5','test_phantom.h5',random_seed)

expanding_phantom = foam_ct_phantom.ExpandingFoamPhantom('test_expand.h5')

vol_geom = foam_ct_phantom.VolumeGeometry(256,256,1,3/256)

expanding_phantom.generate_volume('test_expand_vol_start.h5', vol_geom, time=0)
expanding_phantom.generate_volume('test_expand_vol_end.h5', vol_geom, time=1)

vol0 = foam_ct_phantom.load_volume('test_expand_vol_start.h5')
vol1 = foam_ct_phantom.load_volume('test_expand_vol_end.h5')
pl.imshow(vol0[0])
pl.figure()
pl.imshow(vol1[0])
pl.show()

proj_geom = foam_ct_phantom.ParallelGeometry(256,128,np.linspace(0,4*np.pi,128,False),3/256)
expanding_phantom.generate_projections('test_expand_projs.h5', proj_geom)

projs = foam_ct_phantom.load_projections('test_expand_projs.h5')

pl.imshow(projs[0])
pl.figure()
pl.imshow(projs[-1])
pl.show()