tomostream.solver
¶
Functions:
- class tomostream.solver.Solver(ntheta, n, nz, center, idx, idy, idz, rotx, roty, rotz, fbpfilter, dezinger, data_type)[source]¶
Bases:
object
Class for tomography reconstruction of ortho-slices through direct discreatization of circular integrals in the Radon transform.
- Parameters
ntheta (int) – The number of projections in the buffer (for simultaneous reconstruction)
n, nz (int) – The pixel width and height of the projection.
- recon_optimized(data, theta, ids, center, idx, idy, idz, rotx, roty, rotz, fbpfilter, dezinger, dbg=False)[source]¶
Optimized reconstruction of the object from the whole set of projections in the interval of size pi. Resulting reconstruction is obtained by replacing the reconstruction part corresponding to incoming projections, objnew = objold + recon(datanew) - recon(dataold) whenever the number of incoming projections is less than half of the buffer size. Reconstruction is done with using the whole buffer only when: the number of incoming projections is greater than half of the buffer size, idx/idy/idz, center, fbpfilter are changed, or new dark/flat fields are acquired.
- Parameters
data (np.array(nproj,nz,n)) – Projection data
theta (np.array(nproj)) – Angles corresponding to the projection data
ids (np.array(nproj)) – Ids of the data in the circular buffer array
center (float) – Rotation center for reconstruction
idx, idy, idz (int) – X-Y-Z ortho slices for reconstruction
rotx, roty, rotz (float) – Rotation angles for X-Y-Z slices
fbpfilter (str) – Reconstruction filter
dezinger (str) – None or radius for removing outliers
- Returns
obj (np.array(n,3*n)) – Concatenated reconstructions for X-Y-Z orthoslices