Objective. To allow the estimation of secondary cancer risks from radiation therapy treatment plans in a comprehensive and user-friendly Monte Carlo (MC) framework.Method. Patient planning computed tomography scans were extended superior-inferior using the International Commission on Radiological Protection\'s Publication 145 computational mesh phantoms and skeletal matching. Dose distributions were calculated with the TOPAS MC system using novel mesh capabilities and the digital imaging and communications in medicine radiotherapy extension interface. Finally, in-field and out-of-field cancer risk was calculated using both sarcoma and carcinoma risk models with two alternative parameter sets.Result. The TOPAS MC framework was extended to facilitate epidemiological studies on radiation-induced cancer risk. The framework is efficient and allows automated analysis of large datasets. Out-of-field organ dose was small compared to in-field dose, but the risk estimates indicate a non-negligible contribution to the total radiation induced cancer risk.Significance. This work equips the TOPAS MC system with anatomical extension, mesh geometry, and cancer risk model capabilities that make state-of-the-art out-of-field dose calculation and risk estimation accessible to a large pool of users. Furthermore, these capabilities will facilitate further refinement of risk models and sensitivity analysis of patient specific treatment options.

Objective.Virtual imaging trials enable efficient assessment and optimization of medical image devices and techniques via simulation rather than physical studies. These studies require realistic, detailed ground-truth models or phantoms of the relevant anatomy or physiology. Anatomical structures within computational phantoms are typically based on medical imaging data; however, for small and intricate structures (e.g. trabecular bone), it is not reasonable to use existing clinical data as the spatial resolution of the scans is insufficient. In this study, we develop a mathematical method to generate arbitrary-resolution bone structures within virtual patient models (XCAT phantoms) to model the appearance of CT-imaged trabecular bone.Approach. Given surface definitions of a bone, an algorithm was implemented to generate stochastic bicontinuous microstructures to form a network to define the trabecular bone structure with geometric and topological properties indicative of the bone. For an example adult male XCAT phantom (50th percentile in height and weight), the method was used to generate the trabecular structure of 46 chest bones. The produced models were validated in comparison with published properties of bones. The utility of the method was demonstrated with pilot CT and photon-counting CT simulations performed using the accurate DukeSim CT simulator on the XCAT phantom containing the detailed bone models.Main results. The method successfully generated the inner trabecular structure for the different bones of the chest, having quantiative measures similar to published values. The pilot simulations showed the ability of photon-counting CT to better resolve the trabecular detail emphasizing the necessity for high-resolution bone models.Significance.As demonstrated, the developed tools have great potential to provide ground truth simulations to access the ability of existing and emerging CT imaging technology to provide quantitative information about bone structures.

Preclinical dosimetry is essential for guiding the design of animal radiopharmaceutical biodistribution, imaging, and therapy experiments, evaluating efficacy and/or toxicities in such experiments, ensuring compliance with ethical standards for animal research, and, perhaps most importantly, providing reasonable initial estimates of normal-organ doses in humans, required for clinical translation of new radiopharmaceuticals. This MIB Guide provides a basic protocol for obtaining preclinical dosimetry estimates with organ-level dosimetry software.

OBJECTIVE: Steadily increasing use of computational/virtual phantoms in medical physics has motivated expanding development of new simulation methods and data representations for modelling human anatomy. This has emphasized the need for increased realism, user control, and availability. In breast cancer research, virtual phantoms have gained an important role in evaluating and optimizing imaging systems. For this paper, we have developed an algorithm to model breast abnormalities based on fractal Perlin noise. We demonstrate and characterize the extension of this approach to simulate breast lesions of various sizes, shapes, and complexity.
METHODS: Recently, we developed an algorithm for simulating the 3D arrangement of breast anatomy based on Perlin noise. In this paper, we have expanded the method to also model soft tissue breast lesions. We simulated lesions within the size range of clinically representative breast lesions (masses, 5-20 mm in size). Simulated lesions were blended into simulated breast tissue backgrounds and visualized as virtual digital mammography images. The lesions were evaluated by observers following the BI-RADS assessment criteria.
RESULTS: Observers categorized the lesions as round, oval or irregular, with circumscribed, microlobulated, indistinct or obscured margins. The majority of the simulated lesions were considered by the observers to have a realism score of moderate to well. The simulation method provides almost real-time lesion generation (average time and standard deviation: 1.4 ± 1.0 s).
CONCLUSIONS: We presented a novel algorithm for computer simulation of breast lesions using Perlin noise. The algorithm enables efficient simulation of lesions, with different sizes and appearances.

Chronic obstructive pulmonary disease (COPD) is one of the top three causes of death worldwide, characterized by emphysema and bronchitis. Airway measurements reflect the severity of bronchitis and other airway-related diseases. Airway structures can be objectively evaluated with quantitative computed tomography (CT). The accuracy of such quantifications is limited by the spatial resolution and image noise characteristics of the imaging system and can be potentially improved with the emerging photon-counting CT (PCCT) technology. This study evaluated the quantitative performance of PCCT against energy-integrating CT (EICT) systems for airway measurements, and further identified optimum CT imaging parameters for such quantifications. The study was performed using a novel virtual imaging framework by developing the first library of virtual patients with bronchitis. These virtual patients were developed based on CT images of confirmed COPD patients with varied bronchitis severity. The human models were virtually imaged at 6.3 and 12.6 mGy dose levels using a scanner-specific simulator (DukeSim), synthesizing clinical PCCT and EICT scanners (NAEOTOM Alpha, FLASH, Siemens). The projections were reconstructed with two algorithms and kernels at different matrix sizes and slice thicknesses. The CT images were used to quantify clinically relevant airway measurements (\"Pi10\" and \"WA%\") and compared against their ground truth values. Compared to EICT, PCCT provided more accurate Pi10 and WA% measurements by 63.1% and 68.2%, respectively. For both technologies, sharper kernels and larger matrix sizes led to more reliable bronchitis quantifications. This study highlights the potential advantages of PCCT against EICT in characterizing bronchitis utilizing a virtual imaging platform.

Studies of radiobiological effects in murine rodents exposed to internal radiation in the wild or in laboratory experiments require dosimetric support. The main problem of bone marrow (BM) dosimetry for bone-seeking β-emitters is dosimetric modeling, because the bone is a heterogeneous structure with complex microarchitecture. To date, there are several approaches to calculating the absorbed dose in BM, which mostly use rough geometric approximations. Recently, in the framework of studies of people exposed to 90Sr in the Urals, a new approach (SPSD) has been developed. The aim of the current study was to test for the first time the possibility of extension of the SPSD approach elaborated for humans to mice. For this, computational phantoms of femur bones of laboratory animals (C57BL/6, C57BL/6 J, BALB/c, BALB/cJ) aged 5-8 weeks (growing) and > 8 weeks (adults) were created. The dose factors DFSr-90(BM ← TBV + CBV) to convert the Sr isotope activity concentration in a bone tissue into units of dose rate absorbed in the bone marrow were 1.75 ± 0.42 and 2.57 ± 0.93 μGy day-1 per Bq g-1 for growing and adult animals, respectively, while corresponding values for DFSr-89(BM ← TBV + CBV) were 1.08 ± 0.27 and 1.66 ± 0.67 μGy day-1 per Bq g-1, respectively. These results are about 2.5 times lower than skeleton-average DFs calculated assuming homogenous bone, where source and target coincide. The results of the present study demonstrate the possibility of application of the SPSD approach elaborated for humans to non-human mammals. It is concluded that the study demonstrates the feasibility and appropriateness of application of the SPSD approach elaborated for humans to non-human mammals. This approach opens up new prospects for studying the radiobiological consequences of red bone marrow exposure for both laboratory and wildlife mammals.

Phantoms of different sizes, as indicated by several studies, have a significant impact on the accuracy of dose calculations. Therefore, it is necessary to establish a body-size-dependent series of Chinese standing adult phantoms to improve the accuracy of radiation dosimetry. In this study, the Chinese reference polygon-mesh phantomsCRAM_S/CRAF_Shave been refined and a method for automatically constructing lymph nodes in a mesh phantom has been proposed. Then, based on the refined phantoms, this study has developed 42 anthropometric standing adult computational phantoms, 21 models for each gender, with a height range of 145-185 cm and weight as a function of body mass index corresponding to healthy, overweight and obese. The parameters were extracted from the National Occupational Health Standards (GBZ) document of the People\'s Republic of China, which covers more than 90% of the Chinese population. For a given body height and mass, phantoms are scaled in proportion to a factor reflecting the change of adipose tissue and the internal organs. The remainder is adjusted manually to match the target parameters. In addition, the constructed body-size-specific phantoms have been implemented in the in-house THUDose Monte Carlo code to calculate the dose coefficients (DCs) for external photon exposures in the antero-posterior, postero-anterior and right lateral geometries. The results showed that organ DCs varied significantly with body size at low energies (<2MeV) and high energies (>8MeV) due to the differences in anatomy. Organ DC differences between a phantom of a given size and a reference phantom vary by up to 40% for the same height and up to 400% for the whole phantom. The influence of body size differences on the DCs demonstrates that the body-size-dependent Chinese adult phantoms hold great promise for a wide range of applications in radiation dosimetry.

Objective.CT-mesh hybrid phantoms (or \'hybrid(s)\') made from integrated patient CT data and mesh-type reference computational phantoms (MRCPs) can be beneficial for patient-specific whole-body dose evaluation, but this benefit has yet to be evaluated for second cancer risk prediction. The purpose of this study is to compare the hybrid\'s ability to predict risk throughout the body with a patient-scaled MRCP against ground truth whole-body CTs (WBCTs).Approach.Head and neck active scanning proton treatment plans were created for and simulated on seven hybrids and the corresponding scaled MRCPs and WBCTs. Equivalent dose throughout the body was calculated and input into five second cancer risk models for both excess absolute and excess relative risk (EAR and ERR). The hybrid phantom was evaluated by comparing equivalent dose and risk predictions against the WBCT.Main results.The hybrid most frequently provides whole-body second cancer risk predictions which are closer to the ground truth when compared to a scaled MRCP alone. The performance of the hybrid relative to the scaled MRCP was consistent across ERR, EAR, and all risk models. For all in-field organs, where the hybrid shares the WBCT anatomy, the hybrid was better than or equal to the scaled MRCP for both equivalent dose and risk prediction. For out-of-field organs across all patients, the hybrid\'s equivalent dose prediction was superior than the scaled MRCP in 48% of all comparisons, equivalent for 34%, and inferior for 18%. For risk assessment in the same organs, the hybrid\'s prediction was superior than the scaled MRCP in 51.8% of all comparisons, equivalent in 28.6%, and inferior in 19.6%.Significance.Whole-body risk predictions from the CT-mesh hybrid have shown to be more accurate than those from a reference phantom alone. These hybrids could aid in risk-optimized treatment planning and individual risk assessment to minimize second cancer incidence.

Objective. As cancer survivorship increases, there is growing interest in minimizing the late effects of radiation therapy such as radiogenic second cancer, which may occur anywhere in the body. Assessing the risk of late effects requires knowledge of the dose distribution throughout the whole body, including regions far from the treatment field, beyond the typical anatomical extent of clinical computed tomography (CT) scans.Approach. A hybrid phantom was developed which consists of in-field patient CT images extracted from ground truth whole-body CT scans, out-of-field mesh phantoms scaled to basic patient measurements, and a blended transition region. Four of these hybrid phantoms were created, representing male and female patients receiving proton therapy treatment in pelvic and cranial sites. To assess the performance of the hybrid approach, we simulated treatments using the hybrid phantoms, the scaled and unscaled mesh phantoms, and the ground truth whole-body CTs. We calculated absorbed dose and equivalent dose in and outside of the treatment field, with a focus on neutrons induced in the patient by proton therapy. Proton and neutron dose was calculated using a general purpose Monte Carlo code.Main results. The hybrid phantom provided equal or superior accuracy in calculated organ dose and equivalent dose values relative to those obtained using the mesh phantoms in 78% in all selected organs and calculated dose quantities. Comparatively the default mesh and scaled mesh were equal or superior to the other phantoms in 21% and 28% of cases respectively.Significance. The proposed methodology for hybrid synthesis provides a tool for whole-body organ dose estimation for individual patients without requiring CT scans of their entire body. Such a capability would be useful for personalized assessment of late effects and risk-optimization of treatment plans.

Purpose.Radiation epidemiology studies of childhood cancer survivors treated in the pre-computed tomography (CT) era reconstruct the patients\' treatment fields on computational phantoms. For such studies, the phantoms are commonly scaled to age at the time of radiotherapy treatment because age is the generally available anthropometric parameter. Several reference size phantoms are used in such studies, but reference size phantoms are only available at discrete ages (e.g.: newborn, 1, 5, 10, 15, and Adult). When such phantoms are used for RT dose reconstructions, the nearest discrete-aged phantom is selected to represent a survivor of a specific age. In this work, we (1) conducted a feasibility study to scale reference size phantoms at discrete ages to various other ages, and (2) evaluated the dosimetric impact of using exact age-scaled phantoms as opposed to nearest age-matched phantoms at discrete ages.Methods.We have adopted the University of Florida/National Cancer Institute (UF/NCI) computational phantom library for our studies. For the feasibility study, eight male and female reference size UF/NCI phantoms (5, 10, 15, and 35 years) were downscaled to fourteen different ages which included next nearest available lower discrete ages (1, 5, 10 and 15 years) and the median ages at the time of RT for Wilms\' tumor (3.9 years), craniospinal (8.0 years), and all survivors (9.1 years old) in the Childhood Cancer Survivor Study (CCSS) expansion cohort treated with RT. The downscaling was performed using our in-house age scaling functions (ASFs). To geometrically validate the scaling, Dice similarity coefficient (DSC), mean distance to agreement (MDA), and Euclidean distance (ED) were calculated between the scaled and ground-truth discrete-aged phantom (unscaled UF/NCI) for whole-body, brain, heart, liver, pancreas, and kidneys. Additionally, heights of the scaled phantoms were compared with ground-truth phantoms\' height, and the Centers for Disease Control and Prevention (CDC) reported 50th percentile height. Scaled organ masses were compared with ground-truth organ masses. For the dosimetric assessment, one reference size phantom and seventeen body-size dependent 5-year-old phantoms (9 male and 8 female) of varying body mass indices (BMI) were downscaled to 3.9-year-old dimensions for two different radiation dose studies. For the first study, we simulated a 6 MV photon right-sided flank field RT plan on a reference size 5-year-old and 3.9-year-old (both of healthy BMI), keeping the field size the same in both cases. Percent of volume receiving dose ≥15 Gy (V15) and the mean dose were calculated for the pancreas, liver, and stomach. For the second study, the same treatment plan, but with patient anatomy-dependent field sizes, was simulated on seventeen body-size dependent 5- and 3.9-year-old phantoms with varying BMIs. V15, mean dose, and minimum dose received by 1% of the volume (D1), and by 95% of the volume (D95) were calculated for pancreas, liver, stomach, left kidney (contralateral), right kidney, right and left colons, gallbladder, thoracic vertebrae, and lumbar vertebrae. A non-parametric Wilcoxon rank-sum test was performed to determine if the dose to organs of exact age-scaled and nearest age-matched phantoms were significantly different (p < 0.05).Results.In the feasibility study, the best DSCs were obtained for the brain (median: 0.86) and whole-body (median: 0.91) while kidneys (median: 0.58) and pancreas (median: 0.32) showed poorer agreement. In the case of MDA and ED, whole-body, brain, and kidneys showed tighter distribution and lower median values as compared to other organs. For height comparison, the overall agreement was within 2.8% (3.9 cm) and 3.0% (3.2 cm) of ground-truth UF/NCI and CDC reported 50th percentile heights, respectively. For mass comparison, the maximum percent and absolute differences between the scaled and ground-truth organ masses were within 31.3% (29.8 g) and 211.8 g (16.4%), respectively (across all ages). In the first dosimetric study, absolute difference up to 6% and 1.3 Gy was found for V15and mean dose, respectively. In the second dosimetric study, V15and mean dose were significantly different (p < 0.05) for all studied organs except the fully in-beam organs. D1and D95were not significantly different for most organs (p > 0.05).Conclusion.We have successfully evaluated our ASFs by scaling UF/NCI computational phantoms from one age to another age, which demonstrates the feasibility of scaling any CT-based anatomy. We have found that dose to organs of exact age-scaled and nearest aged-matched phantoms are significantly different (p < 0.05) which indicates that using the exact age-scaled phantoms for retrospective dosimetric studies is a better approach.