Derek W. Jokisch

An image-based skeletal dosimetry model for the ICRP reference adult male—internal electron sources

In this study, a comprehensive electron dosimetry model of the adult male skeletal tissues is presented. The model is constructed using the University of Florida adult male hybrid phantom of Lee et al (2010 Phys. Med. Biol. 55 339-63) and the EGSnrc-based Paired Image Radiation Transport code of Shah et al (2005 J. Nucl. Med. 46 344-53). Target tissues include the active bone marrow, associated with radiogenic leukemia, and total shallow marrow, associated with radiogenic bone cancer. Monoenergetic electron emissions are considered over the energy range 1 keV to 10 MeV for the following sources: bone marrow (active and inactive), trabecular bone (surfaces and volumes), and cortical bone (surfaces and volumes). Specific absorbed fractions are computed according to the MIRD schema, and are given as skeletal-averaged values in the paper with site-specific values reported in both tabular and graphical format in an electronic annex available from http://stacks.iop.org/0031-9155/56/2309/mmedia. The distribution of cortical bone and spongiosa at the macroscopic dimensions of the phantom, as well as the distribution of trabecular bone and marrow tissues at the microscopic dimensions of the phantom, is imposed through detailed analyses of whole-body ex vivo CT images (1 mm resolution) and spongiosa-specific ex vivo microCT images (30 µm resolution), respectively, taken from a 40 year male cadaver. The method utilized in this work includes: (1) explicit accounting for changes in marrow self-dose with variations in marrow cellularity, (2) explicit accounting for electron escape from spongiosa, (3) explicit consideration of spongiosa cross-fire from cortical bone, and (4) explicit consideration of the ICRPs change in the surrogate tissue region defining the location of the osteoprogenitor cells (from a 10 µm endosteal layer covering the trabecular and cortical surfaces to a 50 µm shallow marrow layer covering trabecular and medullary cavity surfaces). Skeletal-averaged values of absorbed fraction in the present model are noted to be very compatible with those weighted by the skeletal tissue distributions found in the ICRP Publication 110 adult male and female voxel phantoms, but are in many cases incompatible with values used in current and widely implemented internal dosimetry software.

NMR microscopy of trabecular bone and its role in skeletal dosimetry

One of the more intractable problems in internal dosimetry is the assessment of energy deposition by alpha and beta particles within trabecular, or cancellous, bone. In the past few years, new technologies have emerged that allow for the direct and nondestructive 3D imaging of trabecular bone with sufficient spatial resolution to characterize trabecular bone structure in a manner needed for radiation dosimetry models. High-field proton nuclear magnetic resonance (NMR) imaging is one such technology. NMR is an ideal modality for imaging trabecular bone due to the sharp contrast in proton density between the bone matrix and bone marrow regions. In this study, images of the trabecular regions within the bodies of a human thoracic vertebra have been obtained at a field strength of 14.1 T. These images were digitally processed to measure chord length distribution data through both the bone trabeculae and marrow cavities. These distributions, which were found to be qualitatively consistent with those measured by F. W. Spiers and colleagues at the University of Leeds using physical sectioning and automated light microscopy, yielded a mean trabecular thickness of 201 microm and a mean marrow cavity thickness of 998 microm. The NMR techniques developed here for vertebral imaging may be extended to other skeletal sites, allowing for improved site-specific skeletal dosimetry.

Voxel effects within digital images of trabecular bone and their consequences on chord-length distribution measurements

Chord-length distributions through the trabecular regions of the skeleton have been investigated since the early 1960s. These distributions have become important features for bone marrow dosimetry; as such, current models rely on the accuracy of their measurements. Recent techniques utilize nuclear magnetic resonance (NMR) microscopy to acquire 3D images of trabecular bone that are then used to measure 3D chord-length distributions by Monte Carlo methods. Previous studies have shown that two voxel effects largely affect the acquisition of these distributions within digital images. One is particularly pertinent as it dramatically changes the shape of the distribution and reduces its mean. An attempt was made to reduce this undesirable effect and good results were obtained for a single-sphere model using minimum acceptable chord (MAC) methods (Jokisch et al 2001 Med. Phys. 28 1493-504). The goal of the present work is to extend the study of these methods to more general models in order to better quantify their consequences. First, a mathematical model of a trabecular bone sample was used to test the usefulness of the MAC methods. The results showed that these methods were not efficient for this simulated bone model. These methods were further tested on a single voxelized sphere over a large range of voxel sizes. The results showed that the MAC methods are voxel-size dependent and overestimate the mean chord length for typical resolutions used with NMR microscopy. The study further suggests that bone and marrow chord-length distributions currently utilized in skeletal dosimetry models are most likely affected by voxel effects that yield values of mean chord length lower than their true values.

Image segmentation of trabecular spongiosa by visual inspection of the gradient magnitude

Recent advances in physical models of skeletal dosimetry utilize high-resolution 3-dimensional microscopic computed tomography images of trabecular spongiosa. These images are coupled to radiation transport codes to assess energy deposition within active bone marrow and trabecular endosteum. These transport codes rely primarily on the segmentation of the spongiosa images into bone and marrow voxels. Image thresholding has been the segmentation of choice for bone sample images because of its extreme simplicity. However, the ability of the segmentation to reproduce the physical boundary between bone and marrow depends on the selection of the threshold value. Statistical models, as well as visual inspection of the image, have been employed extensively to determine the correct threshold. Both techniques are affected by partial volume effect and can provide unexpected results if performed without care. In this study, we propose a new technique to threshold trabecular spongiosa images based on visual inspection of the image gradient magnitude. We first show that the gradient magnitude of the image reaches a maximum along a surface that remains almost independent of partial volume effect and that is a good representation of the physical boundary between bone and marrow. A computer program was then developed to allow a user to compare the position of the iso-surface produced by a threshold with the gradient magnitude. The threshold that produces the iso-surface that best coincides with the maximum gradient is chosen. The technique was finally tested with a set of images of a true bone sample with different resolutions, as well as with three images of a cube of Duocell aluminium foam of known mass and density. Both tests demonstrate the ability of the gradient magnitude technique to retrieve sample volumes or media volume fractions with 1% accuracy at 30 microm voxel size.

Skeletal dosimetry via NMR microscopy: Investigations of sample reproducibility and signal source

Nuclear magnetic resonance microscopy has been used for several years as a means of quantifying the 3D microarchitecture of the cancellous regions of the skeleton. These studies were originally undertaken for the purpose of developing non-invasive techniques for the early detection of osteoporosis and other bone structural changes. Recently, nuclear magnetic resonance microscopy has also been used to acquire this same 3D data for the purpose of both (1) generating chord length data across bone trabeculae and marrow cavities and (2) generating 3D images for direct coupling to Monte Carlo radiation transport codes. In both cases, one is interested in the reproducibility of the dosimetric data obtained from nuclear magnetic resonance microscopy. In the first of two studies, a trabecular bone sample from the femoral head of a 51-y-old male cadaver was subjected to repeated image acquisition, image processing, image coupling, and radiation transport simulations. The resulting absorbed fractions at high electron energies (4 MeV) were shown to vary less than 4% among four different imaging sessions of the same sample. In a separate study, two femoral head samples were imaged under differing conditions of the NMR signal source. In the first case, the samples were imaged with intact marrow. These samples were then subjected to marrow digestion and immersed in Gd-doped water, which then filled the marrow cavities. Energy-dependent absorbed fraction profiles for both the marrow-intact and marrow-free samples showed essentially equivalent results. These studies thus provide encouragement that skeletal dosimetry models of improved patient specificity can be achieved via NMR microscopy in vivo.

Site-specific variability in trabecular bone dosimetry: Considerations of energy loss to cortical bone

With continual advances in radionuclide therapies, increasing emphasis is being placed on improving the patient specificity of dose estimates to marrow tissues. While much work has been focused on determining patient-specific assessments of radionuclide uptake in the skeleton, few studies have been initiated to explore the individual variability of absorbed fraction data for electron and beta-particle sources in various skeletal sites. The most recent values of radionuclide S values used in clinical medicine continue to utilize a formalism in which electrons are transported under a trabecular bone geometry of infinite extent. No provisions are thus made for the fraction of energy lost to the cortical bone cortex of the skeletal site and its surrounding tissues. In the present study, NMR microscopy was performed on trabecular bone samples taken from the femoral head and humeral proximal epiphysis of three subjects: a 51-year male, an 82-year female, and an 86-year female. Following image segmentation and coupling to EGS4, electrons were transported within macrostructural models of the various skeletal sites that explicitly include the spatial extent of the spongiosa, as well as the thickness of the surrounding cortical bone. These energy-dependent profiles of absorbed fractions to marrow tissues were then compared to transport simulations made within an infinite region of spongiosa. Ratios of mean absorbed fraction, as weighted by the beta energy spectra, under both transport methodologies were then assembled for the radionuclides 32P and 90Y. These ratios indicate that corrections to existing radionuclide S values for 32P can vary by as much as 5% for the male, 6% for the 82-year female, and 8% for the 86-year female. For the higher-energy beta spectrum of 90Y, these same corrections can reach 8%, 10%, and 11%, respectively.

Accounting for beta-particle energy loss to cortical bone via paired-image radiation transport (PIRT)

Current methods of skeletal dose assessment in both medical physics (radionuclide therapy) and health physics (dose reconstruction and risk assessment) rely heavily on a single set of bone and marrow cavity chord-length distributions in which particle energy deposition is tracked within an infinite extent of trabecular spongiosa, with no allowance for particle escape to cortical bone. In the present study, we introduce a paired-image radiation transport (PIRT) model which provides a more realistic three-dimensional (3D) geometry for particle transport in the skeletal site at both microscopic and macroscopic levels of its histology. Ex vivo CT scans were acquired of the pelvis, cranial cap, and individual ribs excised from a 66-year male cadaver (BMI of 22.7 kg m(-2)). For the three skeletal sites, regions of trabecular spongiosa and cortical bone were identified and segmented. Physical sections of interior spongiosa were taken and subjected to microCT imaging. Voxels within the resulting microCT images were then segmented and labeled as regions of bone trabeculae, endosteum, active marrow, and inactive marrow through application of image processing algorithms. The PIRT methodology was then implemented within the EGSNRC radiation transport code whereby electrons of various initial energies are simultaneously tracked within both the ex vivo CT macroimage and the CT microimage of the skeletal site. At initial electron energies greater than 50-200 keV, a divergence in absorbed fractions to active marrow are noted between PIRT model simulations and those estimated under existing techniques of infinite spongiosa transport. Calculations of radionuclide S values under both methodologies imply that current chord-based models may overestimate the absorbed dose to active bone marrow in these skeletal sites by 0% to 27% for low-energy beta emitters (33P, 169Er, and 177Lu), by approximately 4% to 49% for intermediate-energy beta emitters (153Sm, 186Re, and 89Sr), and by approximately 14% to 76% for high-energy beta emitters (32p, 188Re, and 90Y). The PIRT methodology allows for detailed modeling of the 3D macrostructure of individual marrow-containing bones within the skeleton thus permitting improved estimates of absorbed fractions and radionuclide S values for intermediate-to-high energy beta emitters.

Considerations of anthropometric, tissue volume, and tissue mass scaling for improved patient specificity of skeletal S values

It is generally acknowledged that reference man (70 kg in mass and 170 cm in height) does not adequately represent the stature and physical dimensions of many patients undergoing radionuclide therapy, and thus scaling of radionuclide S values is required for patient specificity. For electron and beta sources uniformly distributed within internal organs, the mean dose from self-irradiation is noted to scale inversely with organ mass, provided no escape of electron energy occurs at the organ boundaries. In the skeleton, this same scaling approach is further assumed to be correct for marrow dosimetry; nevertheless, difficulties in quantitative assessments of marrow mass in specific skeletal regions of the patient make this approach difficult to implement clinically. Instead, scaling of marrow dose is achieved using various anthropometric parameters that presumably scale in the same proportion. In this study, recently developed three-dimensional macrostructural transport models of the femoral head and humeral epiphysis in three individuals (51-year male, 82-year female, and 86-year female) are used to test the abilities of different anthropometric parameters (total body mass, body surface area, etc.) to properly scale radionuclide S values from reference man models. The radionuclides considered are 33P, 177Lu, 153Sm, 186Re, 89Sr, 166Ho, 32P, 188Re, and 90Y localized in either the active marrow or endosteal tissues of the bone trabeculae. S value scaling is additionally conducted in which the 51-year male subject is assigned as the reference individual; scaling parameters are then expanded to include tissue volumes and masses for both active marrow and skeletal spongiosa. The study concludes that, while no single anthropometric parameter emerges as a consistent scaler of reference man S values, lean body mass is indicated as an optimal scaler when the reference S values are based on 3D transport techniques. Furthermore, very exact patient-specific scaling of radionuclide S values can be achieved if measurements of spongiosa volume and marrow volume fraction (high-resolution CT with image segmentation) are known in both the patient and the reference individual at skeletal sites for which dose estimates are sought. However, the study indicates that measurements of the spongiosa volume alone may be sufficient for reasonable patient-specific scaling of S values for the majority of radionuclides of interest in internal-emitter therapy.

An image-based skeletal dosimetry model for the ICRP reference newborn—internal electron sources

In this study, a comprehensive electron dosimetry model of newborn skeletal tissues is presented. The model is constructed using the University of Florida newborn hybrid phantom of Lee et al (2007 Phys. Med. Biol. 52 3309-33), the newborn skeletal tissue model of Pafundi et al (2009 Phys. Med. Biol. 54 4497-531) and the EGSnrc-based Paired Image Radiation Transport code of Shah et al (2005 J. Nucl. Med. 46 344-53). Target tissues include the active bone marrow (surrogate tissue for hematopoietic stem cells), shallow marrow (surrogate tissue for osteoprogenitor cells) and unossified cartilage (surrogate tissue for chondrocytes). Monoenergetic electron emissions are considered over the energy range 1 keV to 10 MeV for the following source tissues: active marrow, trabecular bone (surfaces and volumes), cortical bone (surfaces and volumes) and cartilage. Transport results are reported as specific absorbed fractions according to the MIRD schema and are given as skeletal-averaged values in the paper with bone-specific values reported in both tabular and graphic format as electronic annexes (supplementary data). The method utilized in this work uniquely includes (1) explicit accounting for the finite size and shape of newborn ossification centers (spongiosa regions), (2) explicit accounting for active and shallow marrow dose from electron emissions in cortical bone as well as sites of unossified cartilage, (3) proper accounting of the distribution of trabecular and cortical volumes and surfaces in the newborn skeleton when considering mineral bone sources and (4) explicit consideration of the marrow cellularity changes for active marrow self-irradiation as applicable to radionuclide therapy of diseased marrow in the newborn child.

Chord‐based versus voxel‐based methods of electron transport in the skeletal tissues

Anatomic models needed for internal dose assessment have traditionally been developed using mathematical surface equations to define organ boundaries, shapes, and their positions within the body. Many researchers, however, are now advocating the use of tomographic models created from segmented patient computed tomography (CT) or magnetic resonance (MR) scans. In the skeleton, however, the tissue structures of the bone trabeculae, marrow cavities, and endosteal layer are exceedingly small and of complex shape, and thus do not lend themselves easily to either stylistic representations or in-vivo CT imaging. Historically, the problem of modeling the skeletal tissues has been addressed through the development of chord-based methods of radiation particle transport, as given by studies at the University of Leeds (Leeds, U.K.) using a 44-year male subject. We have proposed an alternative approach to skeletal dosimetry in which excised sections of marrow-intact cadaver spongiosa are imaged directly via microCT scanning. The cadaver selected for initial investigation of this technique was a 66-year male subject of nominal body mass index (22.7 kg m(-2)). The objectives of the present study were to compare chord-based versus voxel-based methods of skeletal dosimetry using data from the UF 66-year male subject. Good agreement between chord-based and voxel-based transport was noted for marrow irradiation by either bone surface or bone volume sources up to 500-1000 keV (depending upon the skeletal site). In contrast, chord-based models of electron transport yielded consistently lower values of the self-absorbed fraction to marrow tissues than seen under voxel-based transport at energies above 100 keV, a feature directly attributed to the inability of chord-based models to account for nonlinear electron trajectories. Significant differences were also noted in the dosimetry of the endosteal layer (for all source tissues), with chord-based transport predicting a higher fraction of energy deposition than given by voxel-based transport (average factor of about 1.6). The study supports future use of voxel-based skeletal models which (1) permit nonlinear electron trajectories across the skeletal tissues, (2) do not rely on mathematical algorithms for treating the endosteal tissue layer, and (3) do not implicitly assume independence of marrow and bone trajectories as is the case for chord-based skeletal models.