Monte Carlo package for simulating radiographic images of realistic anthropomorphic phantoms described by triangle meshes

X-ray imaging system optimization increases the benefit-to-cost ratio by reducing the radiation dose to the patient while maximizing image quality. We present a new simulation tool for the generation of realistic medical x-ray images for assessment and optimization of complete imaging systems. The Monte Carlo code simulates radiation transport physics using the subroutine package PENELOPE, which accurately simulates the transport of electrons and photons within the typical medical imaging energy range. The new code implements a novel object-oriented geometry package that allows simulations with homogeneous objects of arbitrary shapes described by triangle meshes. The flexibility of this code, which uses the industry standard PLY input-file format, allows the use of detailed anatomical models developed using computer-aided design tools applied to segmented CT and MRI data. The use of triangle meshes highly simplifies the ray-tracing algorithm without reducing the generality of the code, since most surface models can be tessellated into triangles while retaining their geometric details. Our algorithm incorporates an octree spatial data structure to sort the triangles and accelerate the simulation, reaching execution speeds comparable to the original quadric geometry model of PENELOPE. Coronary angiograms were simulated using a tessellated version of the NURBS-based Cardiac-Torso (NCAT) phantom. The phantom models 330 objects, comprised in total of 5 million triangles. The dose received by each organ and the contribution of the different scattering processes to the final image were studied in detail.

[1]  J. Sempau,et al.  Experimental benchmarks of the Monte Carlo code penelope , 2003 .

[2]  Stephen Rudin,et al.  Study of the generalized MTF and DQE for a new microangiographic system , 2004, SPIE Medical Imaging.

[3]  Erik Fosse,et al.  «The Visible Human Project» – et etisk dilemma , 2009 .

[4]  J. Sempau,et al.  PENELOPE-2006: A Code System for Monte Carlo Simulation of Electron and Photon Transport , 2009 .

[5]  Donald Meagher,et al.  Geometric modeling using octree encoding , 1982, Comput. Graph. Image Process..

[6]  Peter A. Wieringa,et al.  Light transport in trabecular bone: Monte Carlo simulation based on 3D triangle meshes , 2006, SPIE Medical Imaging.

[7]  Kyle J. Myers,et al.  Gender-specific statistical models of pathological coronary arteries for generating simulated angiograms , 2006, SPIE Medical Imaging.

[8]  S. Evans Catalogue of Diagnostic X-Ray Spectra and Other Data , 1998 .

[9]  Sung-Joon Ye,et al.  Benchmark of PENELOPE code for low-energy photon transport: dose comparisons with MCNP4 and EGS4. , 2004, Physics in medicine and biology.

[10]  Kyle J. Myers,et al.  Monte Carlo simulated coronary angiograms of realistic anatomy and pathology models , 2007, SPIE Medical Imaging.

[11]  D. Badouel An efficient ray-polygon intersection , 1990 .

[12]  J. Baró,et al.  An algorithm for Monte Carlo simulation of coupled electron-photon transport , 1997 .

[13]  J. Tabary,et al.  Coupling Photon Monte Carlo Simulation and CAD Software. Application to X-ray Nondestructive Evaluation , 2001 .

[14]  Donald Meagher,et al.  Geometric modeling using octree encoding , 1982, Computer Graphics and Image Processing.

[15]  D. R. White,et al.  The composition of body tissues. , 1986, The British journal of radiology.

[16]  Allen Y. Chang A Survey of Geometric Data Structures for Ray Tracing , 2001 .

[17]  J. Sulkimo,et al.  Particle tracking in sophisticated CAD models for simulation purposes , 1996 .

[18]  Andrew S. Glassner,et al.  Space subdivision for fast ray tracing , 1984, IEEE Computer Graphics and Applications.

[19]  Carlos Ureña,et al.  An Efficient Parametric Algorithm for Octree Traversal , 2000, WSCG.

[20]  Tomas Akenine-Möller,et al.  Fast, minimum storage ray/triangle intersection , 1997, J. Graphics, GPU, & Game Tools.

[21]  Richard H. Moore,et al.  Monte Carlo simulation of x-ray scatter based on patient model from digital breast tomosynthesis , 2006, SPIE Medical Imaging.

[22]  C Orsholm,et al.  Geometry package for Monte Carlo simulations on CAD files , 2004 .

[23]  P J Scanlon,et al.  ACC/AHA guidelines for coronary angiography: executive summary and recommendations. A report of the American College of Cardiology/American Heart Association Task Force on Practice Guidelines (Committee on Coronary Angiography) developed in collaboration with the Society for Cardiac Angiography and , 1999, Circulation.

[24]  Andreu Badal,et al.  A package of Linux scripts for the parallelization of Monte Carlo simulations , 2006, Comput. Phys. Commun..

[25]  John Amanatides,et al.  A Fast Voxel Traversal Algorithm for Ray Tracing , 1987, Eurographics.

[26]  William Paul Segars,et al.  Development of a new dynamic NURBS-based cardiac-torso (NCAT) phantom , 2001 .

[27]  B.M.W. Tsui,et al.  Extension of the 4D NCAT phantom to dynamic X-ray CT simulation , 2003, 2003 IEEE Nuclear Science Symposium. Conference Record (IEEE Cat. No.03CH37515).

[28]  G. G. Stokes "J." , 1890, The New Yale Book of Quotations.