An evaluation of an energy independent CT reconstruction algorithm for use in radiotherapy treatment planning

OBJECTIVE

Radiotherapy treatment planning relies upon density information provided by CT for accurate dose calculations. Hounsfield units (HUs) are converted to electron/physical density via an energy dependant calibration curve. Multiple curves are required to make full use of the available accelerating potentials (kVp). The curves are bi-linear with a discontinuity occurring at soft-tissue densities. The commercial algorithm, DirectDensityTM (Siemens Healthcare GmbH), constructs a single calibration curve covering all available kVp. This enables the optimisation of the CT image quality, e.g. in terms of contrast, or the reduction of the imaging dose, whilst rendering the radiotherapy treatment dose calculation robust to the energy used to acquire the CT image. We report our investigations on the clinical utilisation of the DirectDensityTM algorithm for radiotherapy treatments, by using all accelerating potentials, i.e. from 70 kVp up to 140 kVp, available at our CT treatment simulator, in contrast to previous studies that were limited to accelerating potentials spanning a subset of the available kVp.

METHODS

The DirectDensityTM (DD) reconstruction algorithm available on a SOMATOM go.Open Pro CT scanner (Siemens Healthineers) was evaluated using the RayStation v. 9 treatment planning system (RaySearch Laboratories, Stockholm, Sweden) and a CIRS Model 002LFC IMRT Thorax Phantom (SunNuclear, Melbourne, FL), which was imaged at all available kVp with clinical protocols corresponding to various anatomical sites. The DD images were compared to those with the standard reconstruction algorithm acquired only at 120 kVp, as per our routine clinical practice. The effect of increasing kVp on HU is investigated for relevant tissue substitutes. In addition, a dosimetric comparison is performed for a VMAT plan technique with 6 MV X-rays using retrospective patient CT data sets representing four anatomical sites (pelvis, thorax, brain and "head and neck") with five patients for each site. The original dose distributions were calculated on images acquired at 120 kVp using the standard clinical iterative reconstruction (Qr40) and compared with dose distributions recalculated on images reconstructed with the new DD (Sm40) algorithm.

RESULTS

The maximum difference for radiotherapy doses calculated using images of the phantom reconstructed with Qr40 (120 kVp) or DD (all available kVp) was 0.73%. The patient plans on the anatomically representative sites studied here showed a mean PTV dose difference of -0.2% (s.d. 0.7) for D99%, -0.4% (s.d. 0.4%) for D50% and -0.3% (s.d. 0.4%) for D2%. Incidentally, we found a previously unreported decrease in HU, mostly notable for bone type inserts (~34 HU (cortical bone)), at 110 kVp for the DD reconstructed images. The effect was not noted for the standard Qr40 reconstructions.

CONCLUSION

DD has a minimal dosimetric impact in the dose calculations for radiotherapy treatments and could be implemented with existing clinical workflows. Attention should be paid to the HU values for images acquired at 110 kVp (DD algorithm), which warrants further investigation.

ADVANCES IN KNOWLEDGE

This is the first paper where DD was evaluated at all available kVp, leading to the incidental discovery of abnormal HU values at 110 kVp for this algorithm.

Evaluation of the RadCalc collapsed cone dose calculation algorithm against measured data

This work describes the experimental validation of the RadCalc (Lifeline software Inc, Tyler) collapsed cone dose calculation algorithm against measured data for a range of scenarios. 6 MV photon beam data were measured in a large water tank on a Varian TrueBeam linear accelerator. These were input into the RadCalc software, in conjunction with head geometry and output calibration information, then used to create a collapsed cone beam model. The model performance was assessed by comparison against measurement, using a selection of homogeneous and inhomogeneous geometries not incorporated into the original beam model. Dose calculations generated using the collapsed cone algorithm are generally in good agreement with measurement. However, the primary collimating of the linac is not accounted for in the RadCalc model and hence dose in the corners of large fields is significantly overestimated. Percentage depth doses were within 0.5% beyond a depth of 2 cm. The dose in the build-up region was underestimated by RadCalc Version 7.1.4.1, with (Distance To Agreement) discrepancies of up to 3 mm which were corrected in Version 7.2.2.0. Beam profiles for homogeneous phantom comparisons were within 2% in the central 80% of the field with out of field dose underestimated by no more than 3%. Dose comparisons in heterogeneous geometries were acceptable and generally within 3.5%. The largest observed differences were found at density interfaces and a result of the RadCalc dose resolution of 2 mm against 1 mm measured. Absolute dose comparisons demonstrated that RadCalc agreed with measurement to within 1.2% under homogeneous media irradiation geometries. For static beam IMRT deliveries agreement was within 2% or 2 mm of measured data, and for complex VMAT deliveries within 3% or 2 mm. The implementation of the (model-based) photon collapsed cone algorithm in RadCalc shows generally good agreement with measured data over a range of simple and complex scenarios considered.

The first patients treated with MR-CBCT soft-tissue matching in a MR-only prostate radiotherapy pathway

INTRODUCTION

Magnetic Resonance (MR)-only radiotherapy for prostate cancer has previously been reported using fiducial markers for on-treatment verification. MR-Cone Beam Computed Tomography (CBCT) soft-tissue matching does not require invasive fiducial markers and enables MR-only treatments to other pelvic cancers. This study evaluated the first clinical implementation of MR-only prostate radiotherapy using MR-CBCT soft-tissue matching.

METHODS

Twenty prostate patients were treated with MR-only radiotherapy using a synthetic (s)CT-optimised plan with MR-CBCT soft-tissue matching. Two MR sequences were acquired: small Field Of View (FOV) for target delineation and large FOV for organs at risk delineation, sCT generation and on-treatment verification. Patients also received a CT for validation. The prostate was independently contoured on the small FOV MR, copied to the registered CT and modified if there were MR-CT soft-tissue alignment differences (MR-CT volume). This was compared to the MR-only volume with a paired t-test. The treatment plan was recalculated on CT and the doses compared. Independent offline CT-CBCT matches for 5/20 fractions were performed by three therapeutic radiographers using the MR-only contours and compared to the online MR-CBCT matches using two one-sided paired t-tests for equivalence within ±1 mm.

RESULTS

The MR-only volumes were significantly smaller than MR-CT (p = 0.003), with a volume ratio 0.92 ± 0.02 (mean ± standard error). The sCT isocentre dose difference to CT was 0.2 ± 0.1%. MR-CBCT soft-tissue matching was equivalent to CT-CBCT (p < 0.001), with differences of 0.1 ± 0.2 mm (vertical), -0.1 ± 0.2 mm (longitudinal) and 0.0 ± 0.1 mm (lateral).

CONCLUSIONS

MR-only radiotherapy with soft-tissue matching has been successfully clinically implemented. It produced significantly smaller target volumes with high dosimetric and on-treatment matching accuracy.

IMPLICATIONS FOR PRACTICE

MR-only prostate radiotherapy can be safely delivered without using invasive fiducial markers. This enables MR-only radiotherapy to be extended to other pelvic cancers where fiducial markers cannot be used.

Cardiac stereotactic ablative radiotherapy for control of refractory ventricular tachycardia: initial UK multicentre experience

Background

Options for patients with ventricular tachycardia (VT) refractory to antiarrhythmic drugs and/or catheter ablation remain limited. Stereotactic radiotherapy has been described as a novel treatment option.

Methods

Seven patients with recurrent refractory VT, deemed high risk for either first time or redo invasive catheter ablation, were treated across three UK centres with non-invasive cardiac stereotactic ablative radiotherapy (SABR). Prior catheter ablation data and non-invasive mapping were combined with cross-sectional imaging to generate radiotherapy plans with aim to deliver a single 25 Gy treatment. Shared planning and treatment guidelines and prospective peer review were used.

Results

Acute suppression of VT was seen in all seven patients. For five patients with at least 6 months follow-up, overall reduction in VT burden was 85%. No high-grade radiotherapy treatment-related side effects were documented. Three deaths (two early, one late) occurred due to heart failure.

Conclusions

Cardiac SABR showed reasonable VT suppression in a high-risk population where conventional treatment had failed.

Comparison of the RayStation photon Monte Carlo dose calculation algorithm against measured data under homogeneous and heterogeneous irradiation geometries

PURPOSE

This work compares Monte Carlo dose calculations performed using the RayStation treatment planning system against data measured on a Varian Truebeam linear accelerator with 6 MV and 10 MV FFF photon beams.

METHODS

The dosimetric performance of the RayStation Monte Carlo calculations was evaluated in a variety of irradiation geometries employing homogeneous and heterogeneous phantoms. Profile and depth dose comparisons against measurement were carried out in relative mode using the gamma index as a quantitative measure of similarity within the central high dose regions.

RESULTS

The results demonstrate that the treatment planning system dose calculation engine agrees with measurement to within 2%/1 mm for more than 95% of the data points in the high dose regions for all test cases. A systematic underestimation was observed at the tail of the profile penumbra and out of field, with mean differences generally <0.5 mm or 1% of curve dose maximum respectively. Out of field agreement varied between evaluated beam models.

CONCLUSIONS

The RayStation implementation of photon Monte Carlo dose calculations show good agreement with measured data for the range of scenarios considered in this work and is deemed sufficiently accurate for introduction into clinical use.

A systematic review and practical considerations of stereotactic body radiotherapy in the treatment of head and neck cancer

OBJECTIVES

Stereotactic radiotherapy (SBRT) is gaining popularity although its use in head and neck cancer (HNC) is not well defined. The primary objective was to review the published evidence regarding the use of stereotactic radiotherapy in HNC.

METHODS

A literature search was performed by using MEDLINE and EMBASE databases for eligible studies from 2000 to 2019 and 26 relevant studies were identified.

RESULTS

Literature demonstrates a heterogeneous use of this technique with regards to patient population, primary or salvage treatment, dose fractionation regimens, outcomes and follow-up protocols. Carotid blow out syndrome is a risk as with other forms of reirradiation but alternative treatment regimens may reduce this risk.

CONCLUSION

At present there is a lack of evidence regarding SBRT as a primary treatment option for HNC and definitive answers regarding efficacy and tolerability cannot be provided but there is growing evidence that SBRT reirradiation regimens are safe and effective. In lieu of evidence from large Phase III trials, we define appropriate organ at risk constraints and prescription doses, with accurate plan summation approaches. Prospective randomised trials are warranted to validate improved treatment outcomes and acceptable treatment morbidity.

ADVANCES IN KNOWLEDGE

This article provides a comprehensive review of evidence of use of stereotactic radiotherapy in HNC site (either as a primary treatment or as reirradiation). We also provide an evidence-based approach to the implementation and practical consideration of stereotactic radiotherapy in HNC.

Technical Note: Efficient and accurate MRI-only based treatment planning of the prostate using bulk density assignment through atlas-based segmentation

PURPOSE

This study investigates the dosimetric accuracy as well as the robustness of a bulk density assignment approach to magnetic resonance imaging (MRI)-only based treatment planning of the prostate, with bulk density regions automatically identified using atlas-based segmentation (ABS).

METHODS

Twenty prostate radiotherapy patients received planning computed tomography (CT) and MRI scans and were treated with volumetric modulated arc therapy (VMAT). Two bulk densities were set, one for bone and one for soft tissue. The bone contours were created by using ABS followed by manual modification if considered necessary. A range of soft tissue and bone density pairs, between 0.95 and 1.03 g/cm³ with increments of 0.01 for soft tissue, and between 1.15 and 1.65 g/cm³ with increments of 0.05 for bone, were evaluated. Using the density pair giving the lowest dose difference compared to the CT-based dose, dose differences were calculated using both the manually modified bone contours and the bone contours from ABS. Contour overlap measurements between the ABS contours and the manually modified contours were calculated.

RESULTS

The dose comparison shows a very good agreement with the CT when using 0.98 g/cm³ for soft tissue and 1.20 g/cm³ for bone, with a dose difference less than 1 % in average dose in all regions of interest. The mean Dice similarity coefficient for bone was 0.94 and the Mean Distance to Agreement was <1 mm in most cases.

CONCLUSIONS

Using bulk density assignment on MR images with suitable densities for bone and soft tissue results in clinically acceptable dose differences compared to dose calculated on the CT, for both atlas-based and manual bone contours. This demonstrates that an integrated MRI-only pathway utilizing a bulk density assignment for two tissue types is a feasible and robust approach for patients with prostate cancer treated with VMAT.

A UK wide study of current prostate planning practice

Objectives:

The objective of this work was to undertake a non-judgemental study of prostate planning practice across the UK by inviting all departments to undertake the same case.

Methods:

An invitation to take part in the study was sent to the Heads of all UK radiotherapy departments and posted on the UK Medical Physics mailbase. Individuals interested in participating were able to access a single anonymised CT dataset for download with the prostate gland, seminal vesicles, bladder, rectum, bowel, femoral heads, and penile bulb outlined. A brief patient history was also supplied. Participants were asked to create planning target volumes (PTVs) according to their local clinical protocol and plan to give 60 Gy in 20 fractions to the PTV receiving the highest dose. No guidance was given for acceptable organ at risk doses. Dicom plan and dose information was loaded back into ProKnow for analysis by contributors.

Results:

There were 102 plan submissions made to the study representing 48 different UK radiotherapy departments. Seventeen distinct methodologies for creating the prescription PTV from the prostate and seminal vesicles were identified with the ethos of the CHHIP trial protocol for margin growing followed in nearly two-thirds of cases. Positive correlations were found when assessing the doses received by the bladder and rectum against the volume of the PTV to which 60 Gy was prescribed.

Conclusions:

A national planning study whereby staff from a multitude of radiotherapy departments create plans based solely on a single dataset is feasible. The cohort of data was made available to all participants following the study to enable self-assessment and benchmarking against that of their peers.

Advances in knowledge:

This is the first UK wide treatment planning study to investigate local clinical prostate planning practice. This has given UK departments the opportunity to evaluate their planning practices against those of their peers.

Evaluation of the RayStation electron Monte Carlo dose calculation algorithm

The aim of this work was to evaluate the accuracy of the RayStation treatment planning system electron Monte Carlo algorithm against measured data for a range of clinically relevant scenarios. This was done by comparing measured percentage depth dose data (PDD) in water, profiles at oblique incidence and with heterogeneities in the beam path, and output factor data and that generated using the RayStation treatment planning system Monte Carlo VMC++ based calculation algorithm. While electron treatments are widely employed in the radiotherapy setting accurate modelling is challenging (TPS) in the presence of patient being both heterogeneous and nonrectangular. Watertank-based measurements were made on a Varian TrueBeam linear accelerator covering electron beam energies 6 to 18 MeV. These included both normal and oblique incidence, heterogeneous geometries, and irregular shaped cut-outs. The measured geometries were replicated in RayStation and the Monte Carlo dose calculation engine used to generate dosimetric data for comparison against measurement in what were considered clinically relevant settings. Water-based PDDs and profile comparisons showed excellent agreement for all electron beam energies. Profiles measured with oblique beam incidence demonstrated acceptable agreement to the treatment planning system calculations although the correspondence worsened as the angle increased with the planning system overestimating the dose in the shoulder region. Profile measurements under inhomogeneities were generally good. The planning system had a tendency to overestimate dose under the heterogeneity and also demonstrated a broader penumbra than measurement. Of the 170 different output factors calculated in RayStation over the range of electron energies commissioned, 141 were within ± 3% of measured values and 164 within ± 5%. Four of the 6 comparisons beyond 5% were at 18 MeV and all had a cut-out edge within 3 cm of the beam central axis/measurement point. The RayStation implementation of a VMC++ electron Monte Carlo dose calculation algorithm shows good agreement with measured data for a range of scenarios studied and represented sufficient accuracy for clinical use.

Dosimetric evaluation of VMAT for palliative radiotherapy for non-small cell lung carcinoma

OBJECTIVE:

To compare the dosimetric consequences of volumetric modulated arc therapy (VMAT) for high-dose palliative thoracic radiotherapy through comparison with conventionally used isocentric parallel opposed pair (POP) of fields.

METHODS:

20 consecutive patients with non-small cell lung cancer who received 36 Gy in 12 fractions using a POP technique were re-planned using a single VMAT arc. Salient dosimetric parameters were compared between the plans using a paired t-test.

RESULTS:

VMAT demonstrated dosimetric superiority; all PTV dose parameters were significantly improved and importantly the volume of normal lung receiving a high dose was also significantly reduced (mean volume of normal lung receiving 36 Gy was 12.9% in POP vs 1.8% in VMAT, p < 0.005).

CONCLUSION:

The standard POP technique does not take into account tissue densities which results in higher doses to the normal tissue outside the target volume and reduced conformity to the PTV.

ADVANCES IN KNOWLEDGE:

With the help of modern VMAT techniques, it is possible to effectively achieve highly conformal dose delivery which may provide an opportunity to escalate the dose to the tumour in this group of patients.

The accuracy of treatment planning system dose modelling in the presence of brass mesh bolus

AIM

This work assesses the dosimetric accuracy of three commercial treatment planning system (TPS) photon dose calculation algorithms in the presence of brass mesh used as a bolus.

BACKGROUND

Bolus material is used in radiotherapy to provide dose build-up where superficial tissues require irradiation. They are generally water equivalent but high density materials can also be used.

MATERIALS AND METHODS

Dose calculations were performed on Monaco and Masterplan TPS (Elekta AB, Sweden) using phantoms defined by the three DICOM CT image sets of water equivalent blocks (no bolus, 1 layer and 2 layers of brass mesh) exported from the CT scanner. The effect of the mesh on monitor units, build-up dose, phantom exit dose and beam penumbra were compared to measured data.

RESULTS

Dose calculations for 6 and 15 MV photon beams on plain water equivalent phantoms were seen to agree well with measurement validating the basic planning system algorithms and models. Dose in the build-up region, phantom exit dose and beam penumbra were poorly modelled in the presence of the brass mesh. The beam attenuation created by the bolus material was overestimated by all three calculation algorithms, at both photon energies, e.g. 1.6% for one layer and up to 3.1% for two layers at 6 MV. The poor modelling of the physical situation in the build-up region is in part a consequence of the high HU artefact caused by the mesh in the CT image.

CONCLUSIONS

CT imaging is not recommended with the brass mesh bolus in situ due to the poor accuracy of the subsequent TPS modelling.

A practical implementation of the 2010 IPEM high dose rate brachytherapy code of practice for the calibration of 192Ir sources

This paper details a practical method for deriving the reference air kerma rate calibration coefficient for Farmer NE2571 chambers using the U.K. Institute of Physics and Engineering in Medicine (IPEM) code of practice for the determination of the reference air kerma rate for HDR (192)Ir brachytherapy sources based on the National Physical Laboratory (NPL) air kerma standard. The reference air kerma rate calibration coefficient was derived using pressure, temperature and source decay corrected ionization chamber response measurements over three successive (192)Ir source clinical cycles. A secondary standard instrument (a Standard Imaging 1000 Plus well chamber) and four tertiary standard instruments (one additional Standard Imaging 1000 Plus well chamber and three Farmer NE2571 chambers housed in a perspex phantom) were used to provide traceability to the NPL primary standard and enable comparison of performance between the chambers. Conservative and optimized estimates on the expanded uncertainties (k = 2) associated with chamber response, ion recombination and reference air kerma rate calibration coefficient were determined. This was seen to be 2.3% and 0.4% respectively for chamber response, 0.2% and 0.08% respectively for ion recombination and 2.6% and 1.2% respectively for the calibration coefficient. No significant change in ion recombination with source decay was observed over the duration of clinical use of the respective 192Ir sources.

Explaining racial disparities in incidence of and survival from out-of-hospital cardiac arrest

A prospective observational study of 4,653 consecutive cases of out-of-hospital cardiac arrest (OOHCA) occurring in New York City from April 1, 2002, to March 31, 2003, was used to assess racial/ethnic differences in the incidence of OOHCA and 30-day survival after hospital discharge among OOHCA patients. The age-adjusted incidence of OOHCA per 10,000 adults was higher among Blacks than among persons in other racial/ethnic groups, and age-adjusted survival from OOHCA was higher among Whites compared with other groups. In analyses restricted to 3,891 patients for whom complete data on all variables were available, the age-adjusted relative odds of survival from OOHCA among Blacks were 0.4 (95% confidence interval: 0.2, 0.7) as compared with Whites. A full multivariable model accounting for demographic factors, prior functional status, initial cardiac rhythm, and characteristics of the OOHCA event explained approximately 41 percent of the lower age-adjusted survival among Blacks. The lower prevalence of ventricular fibrillation as the initial cardiac rhythm among Blacks relative to Whites was the primary contributor. A combination of factors probably accounts for racial/ethnic disparities in OOHCA survival. Previously hypothesized factors such as delays in emergency medical service response or differences in the likelihood of receipt of cardiopulmonary resuscitation did not appear to be substantial contributors to these racial/ethnic disparities.

Modelling approaches for coastal simulation based on cellular automata: the need and potential

The paper summarizes the theoretical and practical needs for cellular automata (CA)-type models in coastal simulation, and describes early steps in the development of a CA-based model for estuarine sedimentation. It describes the key approaches and formulae used for tidal, wave and sediment processes in a prototype integrated cellular model for coastal simulation designed to simulate estuary sedimentary responses during the tidal cycle in the short-term and climate driven changes in sea-level in the long-term. Results of simple model testing for both one-dimensional and two-dimensional models, and a preliminary parameterization for the Blackwater Estuary, UK, are shown. These reveal a good degree of success in using a CA-type model for water and sediment transport as a function of water level and wave height, but tidal current vectors are not effectively simulated in the approach used. The research confirms that a CA-type model for the estuarine sediment system is feasible, with a real prospect for coupling to existing catchment and nearshore beach/cliff models to produce integrated coastal simulators of sediment response to climate, sea-level change and human actions.