Engineering topics

Verhaegen, Frank

Publications and source records attributed to Verhaegen, Frank.

15 records

DOE OSTI2016

Simulation of pseudo-CT images based on deformable image registration of ultrasound images: A proof of concept for transabdominal ultrasound imaging of the prostate during radiotherapy

Purpose: Imaging of patient anatomy during treatment is a necessity for position verification and for adaptive radiotherapy based on daily dose recalculation. Ultrasound (US) image guided radiotherapy systems are currently available to collect US images at the simulation stage (US{sub sim}), coregistered with the simulation computed tomography (CT), and during all treatment fractions. The authors hypothesize that a deformation field derived from US-based deformable image registration can be used to create a daily pseudo-CT (CT{sub ps}) image that is more representative of the patients’ geometry during treatment than the CT acquired at simulation stage (CT{sub sim}). Methods: The three prostate patients, considered to evaluate this hypothesis, had coregistered CT and US scans on various days. In particular, two patients had two US–CT datasets each and the third one had five US–CT datasets. Deformation fields were computed between pairs of US images of the same patient and then applied to the corresponding US{sub sim} scan to yield a new deformed CT{sub ps} scan. The original treatment plans were used to recalculate dose distributions in the simulation, deformed and ground truth CT (CT{sub gt}) images to compare dice similarity coefficients, maximum absolute distance, and mean absolute distance on CT delineations and gamma index (γ) evaluations on both the Hounsfield units (HUs) and the dose. Results: In the majority, deformation did improve the results for all three evaluation methods. The change in gamma failure for dose (γ{sub Dose}, 3%, 3 mm) ranged from an improvement of 11.2% in the prostate volume to a deterioration of 1.3% in the prostate and bladder. The change in gamma failure for the CT images (γ{sub CT}, 50 HU, 3 mm) ranged from an improvement of 20.5% in the anus and rectum to a deterioration of 3.2% in the prostate. Conclusions: This new technique may generate CT{sub ps} images that are more representative of the actual patient anatomy than the CT{sub sim} scan.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

A validated tumor control probability model based on a meta-analysis of low, intermediate, and high-risk prostate cancer patients treated by photon, proton, or carbon-ion radiotherapy

Purpose: A fully heterogeneous population averaged mechanistic tumor control probability (TCP) model is appropriate for the analysis of external beam radiotherapy (EBRT). This has been accomplished for EBRT photon treatment of intermediate-risk prostate cancer. Extending the TCP model for low and high-risk patients would be beneficial in terms of overall decision making. Furthermore, different radiation treatment modalities such as protons and carbon-ions are becoming increasingly available. Consequently, there is a need for a complete TCP model. Methods: A TCP model was fitted and validated to a primary endpoint of 5-year biological no evidence of disease clinical outcome data obtained from a review of the literature for low, intermediate, and high-risk prostate cancer patients (5218 patients fitted, 1088 patients validated), treated by photons, protons, or carbon-ions. The review followed the preferred reporting item for systematic reviews and meta-analyses statement. Treatment regimens include standard fractionation and hypofractionation treatments. Residual analysis and goodness of fit statistics were applied. Results: The TCP model achieves a good level of fit overall, linear regression results in a p-value of <0.000 01 with an adjusted-weighted-R{sup 2} value of 0.77 and a weighted root mean squared error (wRMSE) of 1.2%, to the fitted clinical outcome data. Validation of the model utilizing three independent datasets obtained from the literature resulted in an adjusted-weighted-R{sup 2} value of 0.78 and a wRMSE of less than 1.8%, to the validation clinical outcome data. The weighted mean absolute residual across the entire dataset is found to be 5.4%. Conclusions: This TCP model fitted and validated to clinical outcome data, appears to be an appropriate model for the inclusion of all clinical prostate cancer risk categories, and allows evaluation of current EBRT modalities with regard to tumor control prediction.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2016

Comparison of proton therapy treatment planning for head tumors with a pencil beam algorithm on dual and single energy CT images

Purpose: Dual energy CT (DECT) has recently been proposed as an improvement over single energy CT (SECT) for stopping power ratio (SPR) estimation for proton therapy treatment planning (TP), thereby potentially reducing range uncertainties. Published literature investigated phantoms. This study aims at performing proton therapy TP on SECT and DECT head images of the same patients and at evaluating whether the reported improved DECT SPR accuracy translates into clinically relevant range shifts in clinical head treatment scenarios. Methods: Two phantoms were scanned at a last generation dual source DECT scanner at 90 and 150 kVp with Sn filtration. The first phantom (Gammex phantom) was used to calibrate the scanner in terms of SPR while the second served as evaluation (CIRS phantom). DECT images of five head trauma patients were used as surrogate cancer patient images for TP of proton therapy. Pencil beam algorithm based TP was performed on SECT and DECT images and the dose distributions corresponding to the optimized proton plans were calculated using a Monte Carlo (MC) simulation platform using the same patient geometry for both plans obtained from conversion of the 150 kVp images. Range shifts between the MC dose distributions from SECT and DECT plans were assessed using 2D range maps. Results: SPR root mean square errors (RMSEs) for the inserts of the Gammex phantom were 1.9%, 1.8%, and 1.2% for SECT phantom calibration (SECT{sub phantom}), SECT stoichiometric calibration (SECT{sub stoichiometric}), and DECT calibration, respectively. For the CIRS phantom, these were 3.6%, 1.6%, and 1.0%. When investigating patient anatomy, group median range differences of up to −1.4% were observed for head cases when comparing SECT{sub stoichiometric} with DECT. For this calibration the 25th and 75th percentiles varied from −2% to 0% across the five patients. The group median was found to be limited to 0.5% when using SECT{sub phantom} and the 25th and 75th percentiles varied from −1% to 2%. Conclusions: Proton therapy TP using a pencil beam algorithm and DECT images was performed for the first time. Given that the DECT accuracy as evaluated by two phantoms was 1.2% and 1.0% RMSE, it is questionable whether the range differences reported here are significant.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2015

HDR {sup 192}Ir source speed measurements using a high speed video camera

Purpose: The dose delivered with a HDR {sup 192}Ir afterloader can be separated into a dwell component, and a transit component resulting from the source movement. The transit component is directly dependent on the source speed profile and it is the goal of this study to measure accurate source speed profiles. Methods: A high speed video camera was used to record the movement of a {sup 192}Ir source (Nucletron, an Elekta company, Stockholm, Sweden) for interdwell distances of 0.25–5 cm with dwell times of 0.1, 1, and 2 s. Transit dose distributions were calculated using a Monte Carlo code simulating the source movement. Results: The source stops at each dwell position oscillating around the desired position for a duration up to (0.026 ± 0.005) s. The source speed profile shows variations between 0 and 81 cm/s with average speed of ∼33 cm/s for most of the interdwell distances. The source stops for up to (0.005 ± 0.001) s at nonprogrammed positions in between two programmed dwell positions. The dwell time correction applied by the manufacturer compensates the transit dose between the dwell positions leading to a maximum overdose of 41 mGy for the considered cases and assuming an air-kerma strength of 48 000 U. The transit dose component is not uniformly distributed leading to over and underdoses, which is within 1.4% for commonly prescribed doses (3–10 Gy). Conclusions: The source maintains its speed even for the short interdwell distances. Dose variations due to the transit dose component are much lower than the prescribed treatment doses for brachytherapy, although transit dose component should be evaluated individually for clinical cases.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2014

A Longitudinal Evaluation of Partial Lung Irradiation in Mice by Using a Dedicated Image-Guided Small Animal Irradiator

Purpose: In lung cancer radiation therapy, the dose constraints are determined mostly by healthy lung toxicity. Preclinical microirradiators are a new tool to evaluate treatment strategies closer to clinical irradiation devices. In this study, we quantified local changes in lung density symptomatic of radiation-induced lung fibrosis (RILF) after partial lung irradiation in mice by using a precision image-guided small animal irradiator integrated with micro-computed tomography (CT) imaging. Methods and Materials: C57BL/6 adult male mice (n=76) were divided into 6 groups: a control group (0 Gy) and groups irradiated with a single fraction of 4, 8, 12, 16, or 20 Gy using 5-mm circular parallel-opposed fields targeting the upper right lung. A Monte Carlo model of the small animal irradiator was used for dose calculations. Following irradiation, all mice were imaged at regular intervals over 39 weeks (10 time points total). Nonrigid deformation was used to register the initial micro-CT scan to all subsequent scans. Results: Significant differences could be observed between the 3 highest (>10 Gy) and 3 lowest irradiation (<10 Gy) dose levels. A mean difference of 120 ± 10 HU between the 0- and 20-Gy groups was observed at week 39. RILF was found to be spatially limited to the irradiated portion of the lung. Conclusions: The data suggest that the severity of RILF in partial lung irradiation compared to large field irradiation in mice for the same dose is reduced, and therefore higher doses can be tolerated.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2014

Comparison of TG-43 and TG-186 in breast irradiation using a low energy electronic brachytherapy source

Purpose: The recently updated guidelines for dosimetry in brachytherapy in TG-186 have recommended the use of model-based dosimetry calculations as a replacement for TG-43. TG-186 highlights shortcomings in the water-based approach in TG-43, particularly for low energy brachytherapy sources. The Xoft Axxent is a low energy (<50 kV) brachytherapy system used in accelerated partial breast irradiation (APBI). Breast tissue is a heterogeneous tissue in terms of density and composition. Dosimetric calculations of seven APBI patients treated with Axxent were made using a model-based Monte Carlo platform for a number of tissue models and dose reporting methods and compared to TG-43 based plans. Methods: A model of the Axxent source, the S700, was created and validated against experimental data. CT scans of the patients were used to create realistic multi-tissue/heterogeneous models with breast tissue segmented using a published technique. Alternative water models were used to isolate the influence of tissue heterogeneity and backscatter on the dose distribution. Dose calculations were performed using Geant4 according to the original treatment parameters. The effect of the Axxent balloon applicator used in APBI which could not be modeled in the CT-based model, was modeled using a novel technique that utilizes CAD-based geometries. These techniques were validated experimentally. Results were calculated using two dose reporting methods, dose to water (D{sub w,m}) and dose to medium (D{sub m,m}), for the heterogeneous simulations. All results were compared against TG-43-based dose distributions and evaluated using dose ratio maps and DVH metrics. Changes in skin and PTV dose were highlighted. Results: All simulated heterogeneous models showed a reduced dose to the DVH metrics that is dependent on the method of dose reporting and patient geometry. Based on a prescription dose of 34 Gy, the average D{sub 90} to PTV was reduced by between ∼4% and ∼40%, depending on the scoring method, compared to the TG-43 result. Peak skin dose is also reduced by 10%–15% due to the absence of backscatter not accounted for in TG-43. The balloon applicator also contributed to the reduced dose. Other ROIs showed a difference depending on the method of dose reporting. Conclusions: TG-186-based calculations produce results that are different from TG-43 for the Axxent source. The differences depend strongly on the method of dose reporting. This study highlights the importance of backscatter to peak skin dose. Tissue heterogeneities, applicator, and patient geometries demonstrate the need for a more robust dose calculation method for low energy brachytherapy sources.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2014

On the Sensitivity of α/β Prediction to Dose Calculation Methodology in Prostate Brachytherapy

Purpose: To study the relationship between the accuracy of the dose calculation in brachytherapy and the estimations of the radiosensitivity parameter, α/β, for prostate cancer. Methods and Materials: In this study, Monte Carlo methods and more specifically the code ALGEBRA was used to produce accurate dose calculations in the case of prostate brachytherapy. Equivalent uniform biologically effective dose was calculated for these dose distributions and was used in an iso-effectiveness relationship with external beam radiation therapy. Results: By considering different levels of detail in the calculations, the estimation for the α/β parameter varied from 1.9 to 6.3 Gy, compared with a value of 3.0 Gy suggested by the American Association of Physicists in Medicine Task Group 137. Conclusions: Large variations of the α/β show the sensitivity of this parameter to dose calculation modality. The use of accurate dose calculation engines is critical for better evaluating the biological outcomes of treatments.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2013

Critical assessment of intramodality 3D ultrasound imaging for prostate IGRT compared to fiducial markers

Purpose: A quantitative 3D intramodality ultrasound (US) imaging system was verified for daily in-room prostate localization, and compared to prostate localization based on implanted fiducial markers (FMs).Methods: Thirteen prostate patients underwent multiple US scans during treatment. A total of 376 US-scans and 817 matches were used to determine the intra- and interoperator variability. Additionally, eight other patients underwent daily prostate localization using both US and electronic portal imaging (EPI) with FMs resulting in 244 combined US-EPI scans. Scanning was performed with minimal probe pressure and a correction for the speed of sound aberration was performed. Uncertainties of both US and FM methods were assessed. User variability of the US method was assessed.Results: The overall US user variability is 2.6 mm. The mean differences between US and FM are: 2.5 {+-} 4.0 mm (LR), 0.6 {+-} 4.9 mm (SI), and -2.3 {+-} 3.6 mm (AP). The intramodality character of this US system mitigates potential errors due to transducer pressure and speed of sound aberrations.Conclusions: The overall accuracy of US (3.0 mm) is comparable to our FM workflow (2.2 mm). Since neither US nor FM can be considered a gold standard no conclusions can be drawn on the superiority of either method. Because US imaging captures the prostate itself instead of surrogates no invasive procedure is required. It requires more effort to standardize US imaging than FM detection. Since US imaging does not involve a radiation burden, US prostate imaging offers an alternative for FM EPI positioning.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2013

Active Breathing Control in Combination With Ultrasound Imaging: A Feasibility Study of Image Guidance in Stereotactic Body Radiation Therapy of Liver Lesions

Purpose: Accurate tumor positioning in stereotactic body radiation therapy (SBRT) of liver lesions is often hampered by motion and setup errors. We combined 3-dimensional ultrasound imaging (3DUS) and active breathing control (ABC) as an image guidance tool. Methods and Materials: We tested 3DUS image guidance in the SBRT treatment of liver lesions for 11 patients with 88 treatment fractions. In 5 patients, 3DUS imaging was combined with ABC. The uncertainties of US scanning and US image segmentation in liver lesions were determined with and without ABC. Results: In free breathing, the intraobserver variations were 1.4 mm in left-right (L-R), 1.6 mm in superior-inferior (S-I), and 1.3 mm anterior-posterior (A-P). and the interobserver variations were 1.6 mm (L-R), 2.8 mm (S-I), and 1.2 mm (A-P). The combined uncertainty of US scanning and matching (inter- and intraobserver) was 4 mm (1 SD). The combined uncertainty when ABC was used reduced by 1.7 mm in the S-I direction. For the L-R and A-P directions, no significant difference was observed. Conclusion: 3DUS imaging for IGRT of liver lesions is feasible, although using anatomic surrogates in the close vicinity of the lesion may be needed. ABC-based breath-hold in midventilation during 3DUS imaging can reduce the uncertainty of US-based 3D table shift correction.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2012

On the significance of density-induced speed of sound variations on US-guided radiotherapy

Purpose: To show the effect of speed of sound (SOS) aberration on ultrasound guided radiotherapy (US-gRT) as a function of implemented workflow. US systems assume that SOS is constant in human soft tissues (at a value of 1540 m/s), while its actual nonuniform distribution produces small but systematic errors of up to a few millimeters in the positions of scanned structures. When a coregistered computerized tomography (CT) scan is available, the US image can be corrected for SOS aberration. Typically, image guided radiotherapy workflows implementing US systems only provide a CT scan at the simulation (SIM) stage. If changes occur in geometry or density distribution between SIM and treatment (TX) stage, SOS aberration can change accordingly, with a final impact on the measured position of structures which is dependent on the workflow adopted. Methods: Four basic scenarios were considered of possible changes between SIM and TX: (1) No changes, (2) only patient position changes (rigid rotation-translation), (3) only US transducer position changes (constrained on patient's surface), and (4) patient tissues thickness changes. Different SOS aberrations may arise from the different scenarios, according to the specific US-gRT workflow used: intermodality (INTER) where TX US scans are compared to SIM CT scans; intramodality (INTRA) where TX US scans are compared to SIM US scans; and INTERc and INTRAc where all US images are corrected for SOS aberration (using density information provided by SIM CT). For an experimental proof of principle, the effect of tissues thickness change was simulated in the different workflows: a dual layered phantom was filled with layers of sunflower oil (SOS 1478 m/s), water (SOS 1482 m/s), and 20% saline solution (SOS 1700 m/s). The phantom was US scanned, the layer thicknesses were increased and the US scans were repeated. The errors resulting from the different workflows were compared. Results: Theoretical considerations show that workflows implementing SOS correction based on SIM-CT scan (INTERc, INTRAc) give null errors in all scenarios except when tissues thickness changes, where an error proportional to the degree of change in SOS maps between SIM and TX ({Delta}SOS) occurs. An uncorrected workflow such as INTER produces in all scenarios a pure SOS error, while uncorrected INTRA produces a null error for rotation-translation of the patient, a {Delta}SOS error for changing tissues thickness and an error proportional to the degree of SOS distribution change along the different lines of view when shifting the transducer. The dual layered phantom demonstrated experimentally that the effect of SOS change between SIM and TX is clinically nonrelevant, being less than the intrinsic resolution of imaging systems, even when a substantial change in thicknesses is applied, provided that a SIM-CT-based SOS aberration correction is applied. Noncorrected workflows produce errors up to 4 mm for INTER and to 3 mm for INTRA in the phantom test. Conclusions: A SOS correction is advantageous for all US-gRT workflows and clinical cases, where the effect of SOS change can be considered a second order effect.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2012

Magnitude of speed of sound aberration corrections for ultrasound image guided radiotherapy for prostate and other anatomical sites

Purpose: The purpose of this work is to assess the magnitude of speed of sound (SOS) aberrations in three-dimensional ultrasound (US) imaging systems in image guided radiotherapy. The discrepancy between the fixed SOS value of 1540 m/s assumed by US systems in human soft tissues and its actual nonhomogeneous distribution in patients produces small but systematic errors of up to a few millimeters in the positions of scanned structures. Methods: A correction, provided by a previously published density-based algorithm, was applied to a set of five prostate, five liver, and five breast cancer patients. The shifts of the centroids of target structures and the change in shape were evaluated. Results: After the correction the prostate cases showed shifts up to 3.6 mm toward the US probe, which may explain largely the reported positioning discrepancies in the literature on US systems versus other imaging modalities. Liver cases showed the largest changes in volume of the organ, up to almost 9%, and shifts of the centroids up to more than 6 mm either away or toward the US probe. Breast images showed systematic small shifts of the centroids toward the US probe with a maximum magnitude of 1.3 mm. Conclusions: The applied correction in prostate and liver cancer patients shows positioning errors of several mm due to SOS aberration; the errors are smaller in breast cancer cases, but possibly becoming more important when breast tissue thickness increases.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2012

Dose reduction in LDR brachytherapy by implanted prostate gold fiducial markers

Purpose: The dosimetric impact of gold fiducial markers (FM) implanted prior to external beam radiotherapy of prostate cancer on low dose rate (LDR) brachytherapy seed implants performed in the context of combined therapy was investigated. Methods: A virtual water phantom was designed containing a single FM. Single and multi source scenarios were investigated by performing Monte Carlo dose calculations, along with the influence of varying orientation and distance of the FM with respect to the sources. Three prostate cancer patients treated with LDR brachytherapy for a recurrence following external beam radiotherapy with implanted FM were studied as surrogate cases to combined therapy. FM and brachytherapy seeds were identified on post implant CT scans and Monte Carlo dose calculations were performed with and without FM. The dosimetric impact of the FM was evaluated by quantifying the amplitude of dose shadows and the volume of cold spots. D{sub 90} was reported based on the post implant CT prostate contour. Results: Large shadows are observed in the single source-FM scenarios. As expected from geometric considerations, the shadows are dependent on source-FM distance and orientation. Large dose reductions are observed at the distal side of FM, while at the proximal side a dose enhancement is observed. In multisource scenarios, the importance of shadows appears mitigated, although FM at the periphery of the seed distribution caused underdosage (<prescription dose). In clinical cases, the FM reduced the dose to some voxels by up to 50% and generated shadows with extents of the order of 4 mm. Within the prostate contour, cold spots (<95% prescription dose) of the order of 20 mm{sup 3} were observed. D{sub 90} proved insensitive to the presence of FM for the cases selected. Conclusions: There is a major local impact of FM present in LDR brachytherapy seed implant dose distributions. Therefore, reduced tumor control could be expected from FM implanted in tumors, although our results are too limited to draw conclusions regarding clinical significance.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2011

Patient-Specific Monte Carlo-Based Dose-Kernel Approach for Inverse Planning in Afterloading Brachytherapy

Purpose: Brachytherapy planning software relies on the Task Group report 43 dosimetry formalism. This formalism, based on a water approximation, neglects various heterogeneous materials present during treatment. Various studies have suggested that these heterogeneities should be taken into account to improve the treatment quality. The present study sought to demonstrate the feasibility of incorporating Monte Carlo (MC) dosimetry within an inverse planning algorithm to improve the dose conformity and increase the treatment quality. Methods and Materials: The method was based on precalculated dose kernels in full patient geometries, representing the dose distribution of a brachytherapy source at a single dwell position using MC simulations and the Geant4 toolkit. These dose kernels are used by the inverse planning by simulated annealing tool to produce a fast MC-based plan. A test was performed for an interstitial brachytherapy breast treatment using two different high-dose-rate brachytherapy sources: the microSelectron iridium-192 source and the electronic brachytherapy source Axxent operating at 50 kVp. Results: A research version of the inverse planning by simulated annealing algorithm was combined with MC to provide a method to fully account for the heterogeneities in dose optimization, using the MC method. The effect of the water approximation was found to depend on photon energy, with greater dose attenuation for the lower energies of the Axxent source compared with iridium-192. For the latter, an underdosage of 5.1% for the dose received by 90% of the clinical target volume was found. Conclusion: A new method to optimize afterloading brachytherapy plans that uses MC dosimetric information was developed. Including computed tomography-based information in MC dosimetry in the inverse planning process was shown to take into account the full range of scatter and heterogeneity conditions. This led to significant dose differences compared with the Task Group report 43 approach for the Axxent source.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2011

The difference of scoring dose to water or tissues in Monte Carlo dose calculations for low energy brachytherapy photon sources

Purpose: The goal of this work is to compare D{sub m,m} (radiation transported in medium; dose scored in medium) and D{sub w,m} (radiation transported in medium; dose scored in water) obtained from Monte Carlo (MC) simulations for a subset of human tissues of interest in low energy photon brachytherapy. Using low dose rate seeds and an electronic brachytherapy source (EBS), the authors quantify the large cavity theory conversion factors required. The authors also assess whether applying large cavity theory utilizing the sources' initial photon spectra and average photon energy induces errors related to spatial spectral variations. First, ideal spherical geometries were investigated, followed by clinical brachytherapy LDR seed implants for breast and prostate cancer patients. Methods: Two types of dose calculations are performed with the GEANT4 MC code. (1) For several human tissues, dose profiles are obtained in spherical geometries centered on four types of low energy brachytherapy sources: {sup 125}I, {sup 103}Pd, and {sup 131}Cs seeds, as well as an EBS operating at 50 kV. Ratios of D{sub w,m} over D{sub m,m} are evaluated in the 0-6 cm range. In addition to mean tissue composition, compositions corresponding to one standard deviation from the mean are also studied. (2) Four clinical breast (using {sup 103}Pd) and prostate (using {sup 125}I) brachytherapy seed implants are considered. MC dose calculations are performed based on postimplant CT scans using prostate and breast tissue compositions. PTV D{sub 90} values are compared for D{sub w,m} and D{sub m,m}. Results: (1) Differences (D{sub w,m}/D{sub m,m}-1) of -3% to 70% are observed for the investigated tissues. For a given tissue, D{sub w,m}/D{sub m,m} is similar for all sources within 4% and does not vary more than 2% with distance due to very moderate spectral shifts. Variations of tissue composition about the assumed mean composition influence the conversion factors up to 38%. (2) The ratio of D{sub 90(w,m)} over D{sub 90(m,m)} for clinical implants matches D{sub w,m}/D{sub m,m} at 1 cm from the single point sources. Conclusions: Given the small variation with distance, using conversion factors based on the emitted photon spectrum (or its mean energy) of a given source introduces minimal error. The large differences observed between scoring schemes underline the need for guidelines on choice of media for dose reporting. Providing such guidelines is beyond the scope of this work.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2011

An algorithm for efficient metal artifact reductions in permanent seed implants

Purpose: In permanent seed implants, 60 to more than 100 small metal capsules are inserted in the prostate, creating artifacts in x-ray computed tomography (CT) imaging. The goal of this work is to develop an automatic method for metal artifact reduction (MAR) from small objects such as brachytherapy seeds for clinical applications. Methods: The approach for MAR is based on the interpolation of missing projections by directly using raw helical CT data (sinogram). First, an initial image is reconstructed from the raw CT data. Then, the metal objects segmented from the reconstructed image are reprojected back into the sinogram space to produce a metal-only sinogram. The Steger method is used to determine precisely the position and edges of the seed traces in the raw CT data. By combining the use of Steger detection and reprojections, the missing projections are detected and replaced by interpolation of non-missing neighboring projections. Results: In both phantom experiments and patient studies, the missing projections have been detected successfully and the artifacts caused by metallic objects have been substantially reduced. The performance of the algorithm has been quantified by comparing the uniformity between the uncorrected and the corrected phantom images. The results of the artifact reduction algorithm are indistinguishable from the true background value. Conclusions: An efficient algorithm for MAR in seed brachytherapy was developed. The test results obtained using raw helical CT data for both phantom and clinical cases have demonstrated that the proposed MAR method is capable of accurately detecting and correcting artifacts caused by a large number of very small metal objects (seeds) in sinogram space. This should enable a more accurate use of advanced brachytherapy dose calculations, such as Monte Carlo simulations.

60 APPLIED LIFE SCIENCES↗