Metallic artifact mitigation and organ-constrained tissue assignment for Monte Carlo calculations of permanent implant lung brachytherapy.

PURPOSE To investigate methods of generating accurate patient-specific computational phantoms for the Monte Carlo calculation of lung brachytherapy patient dose distributions. METHODS Four metallic artifact mitigation methods are applied to six lung brachytherapy patient computed tomography (CT) images: simple threshold replacement (STR) identifies high CT values in the vicinity of the seeds and replaces them with estimated true values; fan beam virtual sinogram replaces artifact-affected values in a virtual sinogram and performs a filtered back-projection to generate a corrected image; 3D median filter replaces voxel values that differ from the median value in a region of interest surrounding the voxel and then applies a second filter to reduce noise; and a combination of fan beam virtual sinogram and STR. Computational phantoms are generated from artifact-corrected and uncorrected images using several tissue assignment schemes: both lung-contour constrained and unconstrained global schemes are considered. Voxel mass densities are assigned based on voxel CT number or using the nominal tissue mass densities. Dose distributions are calculated using the EGSnrc user-code BrachyDose for (125)I, (103)Pd, and (131)Cs seeds and are compared directly as well as through dose volume histograms and dose metrics for target volumes surrounding surgical sutures. RESULTS Metallic artifact mitigation techniques vary in ability to reduce artifacts while preserving tissue detail. Notably, images corrected with the fan beam virtual sinogram have reduced artifacts but residual artifacts near sources remain requiring additional use of STR; the 3D median filter removes artifacts but simultaneously removes detail in lung and bone. Doses vary considerably between computational phantoms with the largest differences arising from artifact-affected voxels assigned to bone in the vicinity of the seeds. Consequently, when metallic artifact reduction and constrained tissue assignment within lung contours are employed in generated phantoms, this erroneous assignment is reduced, generally resulting in higher doses. Lung-constrained tissue assignment also results in increased doses in regions of interest due to a reduction in the erroneous assignment of adipose to voxels within lung contours. Differences in dose metrics calculated for different computational phantoms are sensitive to radionuclide photon spectra with the largest differences for (103)Pd seeds and smallest but still considerable differences for (131)Cs seeds. CONCLUSIONS Despite producing differences in CT images, dose metrics calculated using the STR, fan beam + STR, and 3D median filter techniques produce similar dose metrics. Results suggest that the accuracy of dose distributions for permanent implant lung brachytherapy is improved by applying lung-constrained tissue assignment schemes to metallic artifact corrected images.

[1]  S. Kalnicki,et al.  Intraoperative 125I brachytherapy for high-risk stage I non-small cell lung carcinoma. , 1999, International journal of radiation oncology, biology, physics.

[2]  J G H Sutherland,et al.  A Monte Carlo investigation of lung brachytherapy treatment planning , 2013, Physics in medicine and biology.

[3]  R. E. P. Taylor,et al.  Sci—Sat AM(2): Brachy — 05: Fast Monte Carlo Dose Calculations for Brachytherapy with BrachyDose , 2010 .

[4]  R. Landreneau,et al.  Intraoperative brachytherapy following thoracoscopic wedge resection of stage I lung cancer. , 1998, Chest.

[5]  L. Anderson,et al.  Dosimetry of interstitial brachytherapy sources: Recommendations of the AAPM Radiation Therapy Committee Task Group No. 43 , 1995 .

[6]  R. Mohan,et al.  Converting absorbed dose to medium to absorbed dose to water for Monte Carlo based photon beam dose calculations. , 2000, Physics in medicine and biology.

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

[8]  Harald Paganetti,et al.  Dose to water versus dose to medium in proton beam therapy , 2009, Physics in medicine and biology.

[9]  R. Keenan,et al.  Mature follow-up for high-risk stage I non-small-cell lung carcinoma treated with sublobar resection and intraoperative iodine-125 brachytherapy. , 2011, International journal of radiation oncology, biology, physics.

[10]  J. Williamson,et al.  Update of AAPM Task Group No. 43 Report: A revised AAPM protocol for brachytherapy dose calculations. , 2004 .

[11]  J. Williamson,et al.  Report of the Task Group 186 on model-based dose calculation methods in brachytherapy beyond the TG-43 formalism: current status and recommendations for clinical implementation. , 2012, Medical physics.

[12]  R. Cormack,et al.  Permanent iodine-125 interstitial planar seed brachytherapy for close or positive margins for thoracic malignancies. , 2010, International journal of radiation oncology, biology, physics.

[13]  Y. Garces,et al.  Model-based dose calculations for 125 I lung brachytherapy. , 2012, Medical physics.

[14]  A. Ahmadian,et al.  Reduction of dental filling metallic artifacts in CT-based attenuation correction of PET data using weighted virtual sinograms optimized by a genetic algorithm. , 2010, Medical physics.

[15]  T. Bortfeld,et al.  Correlation between CT numbers and tissue parameters needed for Monte Carlo simulations of clinical dose distributions. , 2000, Physics in medicine and biology.

[16]  Luc Beaulieu,et al.  Correction of CT artifacts and its influence on Monte Carlo dose calculations. , 2007, Medical physics.

[17]  Tae-Suk Suh,et al.  Overlapped seed localization in seed implant brachytherapy. , 2009, Medical engineering & physics.

[18]  D. Heron,et al.  The impact of adjuvant brachytherapy with sublobar resection on pulmonary function and dyspnea in high-risk patients with operable disease: preliminary results from the American College of Surgeons Oncology Group Z4032 trial. , 2011, The Journal of thoracic and cardiovascular surgery.

[19]  M. Rivard.  Brachytherapy dosimetry parameters calculated for a Cs131 source: Brachytherapy dosimetry parameters calculated for a Cs131 source , 2007 .

[20]  M. Rivard,et al.  Limited resection for non-small cell lung cancer: observed local control with implantation of I-125 brachytherapy seeds. , 2003, The Annals of thoracic surgery.

[21]  H. Fernando,et al.  Sublobar resection with brachytherapy mesh for stage I non-small cell lung cancer. , 2010, Seminars in thoracic and cardiovascular surgery.

[22]  Mehran Yazdi,et al.  An adaptive approach to metal artifact reduction in helical computed tomography for radiation therapy treatment planning: experimental and clinical studies. , 2005, International journal of radiation oncology, biology, physics.

[23]  R. Taylor,et al.  An EGSnrc Monte Carlo-calculated database of TG-43 parameters. , 2008, Medical physics.

[24]  F. Verhaegen,et al.  Sensitivity study for CT image use in Monte Carlo treatment planning , 2005, Physics in medicine and biology.

[25]  B. L. Kirk,et al.  Brachytherapy dosimetry parameters calculated for a new 103Pd source. , 2004, Medical physics.

[26]  R. Landreneau,et al.  Comparison between sublobar resection and 125Iodine brachytherapy after sublobar resection in high-risk patients with Stage I non-small-cell lung cancer. , 2003, Surgery.

[27]  Xiaochuan Pan,et al.  Reduction of computed tomography metal artifacts due to the Fletcher-Suit applicator in gynecology patients receiving intracavitary brachytherapy. , 2003, Brachytherapy.

[28]  I. Kawrakow.  Accurate condensed history Monte Carlo simulation of electron transport. I. EGSnrc, the new EGS4 version. , 2000, Medical physics.

[29]  R. Robb,et al.  Prostate brachytherapy seed localization by analysis of multiple projections: identifying and addressing the seed overlap problem. , 2004, Medical physics.

[30]  Ravinder Nath,et al.  AAPM recommendations on dose prescription and reporting methods for permanent interstitial brachytherapy for prostate cancer: Report of Task Group 137. , 2009, Medical physics.

[31]  P. Devlin,et al.  Intraoperative seed placement for thoracic malignancy-A review of technique, indications, and published literature. , 2009, Brachytherapy.

[32]  Frank Verhaegen,et al.  Sensitivity of low energy brachytherapy Monte Carlo dose calculations to uncertainties in human tissue composition. , 2010, Medical physics.

[33]  D. Heron,et al.  Intraoperative (125)I Vicryl mesh brachytherapy after sublobar resection for high-risk stage I non-small cell lung cancer. , 2004, Brachytherapy.

[34]  W. Kalender,et al.  Reduction of CT artifacts caused by metallic implants. , 1987 .

[35]  Claus F Behrens,et al.  CTC-ask: a new algorithm for conversion of CT numbers to tissue parameters for Monte Carlo dose calculations applying DICOM RS knowledge. , 2011, Physics in medicine and biology.

[36]  Moyed Miften,et al.  Dosimetric and technical aspects of intraoperative I-125 brachytherapy for stage I non-small cell lung cancer , 2007, Physics in medicine and biology.

[37]  W. Butler,et al.  Supplement to the 2004 update of the AAPM Task Group No. 43 Report. , 2007, Medical physics.

[38]  P. Basran,et al.  CT image artifacts from brachytherapy seed implants: a postprocessing 3D adaptive median filter. , 2011, Medical physics.

[39]  R. Taylor,et al.  Benchmarking brachydose: Voxel based EGSnrc Monte Carlo calculations of TG-43 dosimetry parameters. , 2007, Medical physics.

[40]  Habib Zaidi,et al.  A virtual sinogram method to reduce dental metallic implant artefacts in computed tomography-based attenuation correction for PET , 2010, Nuclear medicine communications.

[41]  Frank Verhaegen,et al.  Simulation study on potential accuracy gains from dual energy CT tissue segmentation for low-energy brachytherapy Monte Carlo dose calculations , 2011, Physics in medicine and biology.