A python implementation of three-dimensional Shepp-Logan head phantom
Three variants, "Kak-Slaney" [1], "Yu-Ye-Wang"[2], "Toft-Schabel"[3] are implemented.
from matplotlib import pyplot as plt
nVoxelZ = 64
nVoxelY = 256
nVoxelX = 256
phantom0, e = shepp_logan_3d(size_out=[nVoxelZ,nVoxelY,nVoxelX], phantom_type="kak-slaney")
phantom1, e = shepp_logan_3d(size_out=[nVoxelZ,nVoxelY,nVoxelX], phantom_type="yu-ye-wang")
phantom2, e = shepp_logan_3d(size_out=[nVoxelZ,nVoxelY,nVoxelX], phantom_type="toft-schabel")
plt.imshow(np.concatenate(
[phantom0[3*phantom0.shape[0]//8],
phantom1[3*phantom1.shape[0]//8],
phantom2[ phantom2.shape[0]//2]],
axis=1))
plt.colorbar()
plt.title("KS z=-0.25 YYW z=-0.25 TS z=0.00")
plt.show()
3-clause BSD License
[1] Kak AC, Slaney M, Principles of Computerized Tomographic Imaging, 1988. p.102 http:https://www.slaney.org/pct/pct-errata.html
[2] Yu H, Ye Y, Wang G, Katsevich-Type Algorithms for Variable Radius Spiral Cone-Beam CT Proceedings of the SPIE, Volume 5535, p. 550-557 (2004)
[3] Matthias Schabel (2021). 3D Shepp-Logan phantom (https://www.mathworks.com/matlabcentral/fileexchange/9416-3d-shepp-logan-phantom), MATLAB Central File Exchange. Retrieved April 29, 2021.