Construction of a Voxel Model from CT Images with Density Derived from CT Numbers

Abstract

The voxel models representing human anatomy have been developed to calculate dose distribution in human body, while the density and elemental composition are the most important physical properties of voxel model. Usually, when creating the Monte Carlo input files, the average tissue densities recommended in ICRP Publication were used to assign each voxel in the existing voxel models. As each tissue consists of many voxels with different densities, the conventional method of average tissue densities failed to take account of the voxel's discrepancy, and therefore could not represent human anatomy faithfully. To represent human anatomy more faithfully, a method was implemented to assign each voxel, the densities of which were derived from CT number. In order to compare with the traditional method, we constructed two models from the cadaver specimen dataset. A CT-based pelvic voxel model called Pelvis-CT model was constructed, the densities of which were derived from the CT numbers. A color photograph-based pelvic voxel model called Pelvis-Photo model was also constructed, the densities of which were taken from ICRP Publication. The CT images and the color photographs were obtained from the same female cadaver specimen. The Pelvis-CT and Pelvis-Photo models were both ported into Monte Carlo code MCNP to calculate the conversion coefficients from kerma free-in-air to absorbed dose for external monoenergetic photon beams with energies of 0.1, 1 and 10 MeV under anterior-posterior (AP) geometry. The results were compared with those of given in ICRP Publication 74. Differences of up to 50% were observed between conversion coefficients of Pelvis-CT and Pelvis-Photo models, moreover the discrepancies decreased for the photon beams with higher energies. The overall trend of conversion coefficients of the Pelvis-CT model agreed well with that of ICRP Publication 74 data.

CAD-based Monte Carlo program for integrated simulation of nuclear system SuperMC

“…''Invertor'' can reconstruct CAD models or facet models from simulation models for observing and modifying models by 3D visualization Li et al, 2007;Zeng et al, 2006;Hu et al, 2007). Besides CAD geometry, CT, MRI, segmented images and other scan data can be converted for human body dose assessment in radiation shielding or medical physics (Cheng et al, 2011).…”