Skip to main content

Cone beam CT-based dose accumulation and analysis of delivered dose to the dominant intraprostatic lesion in primary radiotherapy of prostate cancer



Evaluation of delivered dose to the dominant intraprostatic lesion (DIL) for moderately hypofractionated radiotherapy of prostate cancer by cone beam computed tomography (CBCT)-based dose accumulation and target coverage analysis.


Twenty-three patients with localized prostate cancer treated with moderately hypofractionated prostate radiotherapy with simultaneous integrated boost (SIB) between December 2016 and February 2020 were retrospectively analyzed. Included patients were required to have an identifiable DIL on bi-parametric planning magnetic resonance imaging (MRI). After import into the RayStation treatment planning system and application of a step-wise density override, the fractional doses were computed on each CBCT and were consecutively mapped onto the planning CT via a deformation vector field derived from deformable image registration. Fractional doses were accumulated for all CBCTs and interpolated for missing CBCTs, resulting in the delivered dose for PTVDIL, PTVBoost, PTV, and the organs at risk. The location of the index lesions was recorded according to the sector map of the Prostate Imaging Reporting and Data System (PIRADS) Version 2.1. Target coverage of the index lesions was evaluated and stratified for location.


In total, 338 CBCTs were available for analysis. Dose accumulation target coverage of PTVDIL, PTVBoost, and PTV was excellent and no cases of underdosage in DMean, D95%, D02%, and D98% could be detected. Delivered rectum DMean did not significantly differ from the planned dose. Bladder mean DMean was higher than planned with 19.4 ± 7.4 Gy versus 18.8 ± 7.5 Gy, p < 0.001. The penile bulb showed a decreased delivered mean DMean with 29.1 ± 14.0 Gy versus 29.8 ± 14.4 Gy, p < 0.001. Dorsal DILs, defined as DILs in the posterior medial peripheral zone of the prostate, showed a significantly lower delivered dose with a mean DMean difference of 2.2 Gy (95% CI 1.3–3.1 Gy, p < 0.001) compared to ventral lesions.


CBCT-based dose accumulation showed an adequate delivered dose to the dominant intraprostatic lesion and organs at risk within planning limits. Cautious evaluation of the target coverage for index lesions adjacent to the rectum is warranted to avoid underdosage.


Radiotherapy is an established treatment option for one of the most common cancer types in men: prostate cancer [1]. While whole-gland radiotherapy is the standard treatment modality for prostate cancer, focal therapy may result in decreased toxicity while maintaining high tumor control. For the prostate, the relevant target structure for focal therapy is the dominant intraprostatic lesion (DIL) or index lesion as local recurrence often occurs at the site of the initial primary tumor which influences prognosis and disease progression [2,3,4,5]. The implementation of focal ablative micro boosts to the DIL, as investigated in the FLAME trial, depends on accurate planning and treatment delivery [6]. To achieve the goal of precise and highly conformal dose distributions, MRI- and CBCT-based positioning monitoring methods are commonly implemented [7]. Nevertheless, planned and delivered doses may differ due to motion and daily variations in rectum and bladder filling states. Especially for focal targeted therapy of the index lesion with tight margins, adequate dose delivery is of utmost importance to avoid target miss and underdosage. Dose accumulation strategies may therefore guide the clinical decision to adapt the current radiotherapy plan during treatment. In this study, a CBCT-based dose accumulation for the dominant intraprostatic lesion in moderately hypofractionated volumetric arc radiotherapy (VMAT) with SIB was performed as a step towards dose-guided adaptive prostate radiotherapy.


Patient characteristics and treatment planning

Twenty-three patients with histologically proven low- to high-risk prostate cancer and identifiable DIL in the pre-radiotherapy planning MRI underwent image-guided radiotherapy (IGRT) with CBCT-guidance at an Elekta Synergy Agility linear accelerator (Elekta AB, Stockholm, Sweden) from December 2016 to February 2020. All patients received computed tomography (CT) and bi-parametric MRI scans for radiotherapy planning in a supine position. For the CT scan and the radiotherapy appliance, the patients were instructed to have a full bladder and an empty rectum. Planning MRI protocol consisted of a T2-weighted turbo spin-echo sequence and diffusion-weighted imaging [8]. Radiotherapy was delivered with image-guided VMAT in 33 fractions with simultaneous integrated boost. Regarding the patient set-up workflow, a bone to bone match restricted to the target area was applied first. Afterwards, a consecutive manual soft-tissue matching to the target volumes was performed. Three degrees of freedom shifts were applied without rotational compensation. CBCT scans were performed on each of the first five fractions and once every three fractions afterwards. Contouring, dose prescription, and treatment planning have been reported in detail in earlier publications [8, 9]. In short, prostate radiotherapy was planned with two dose levels of 1.82 Gy for the low dose PTV and 2.31 Gy for the high dose PTVBoost per fraction, resulting in a prescribed PTV dose of 60.06 Gy (D95%) and a PTVBoost mean dose of 76.23 Gy. A clinical target volume (CTVP−SV) was generated consisting of the prostate and only the base of the seminal vesicles, whereas the CTVP+SV included the prostate and the whole seminal vesicles. PTVBoost was defined by a 5 mm margin around CTVP−SV with avoidance of the rectum. The PTV was created by a 10 mm margin around CTVP+SV in all but the dorsal direction, where a 7 mm margin was used. The planning MRI was rigidly co-registered to the planning CT with soft-tissue matching in the target area. Contouring of the clinical target volumes took place on the fused MRI and planning CT without additional contour transferring. PTVDIL was defined as the DIL without margin and contoured by an experienced radiation oncologist according to the co-registered planning MRI sequences. The contouring of PTVDIL is illustrated in Fig. 1. Pinnacle3 version 16.2 (Philips Radiation Oncology Systems, Fitchburg, WI, USA) was used as treatment planning system (TPS). RayStation version 6.99 (RaySearch Laboratories AB, Stockholm, Sweden) was used for analysis.

Fig. 1

Dominant intraprostatic lesion contouring. Shown is the index lesion in the peripheral zone on the a planning CT, b native T2-weighted image, c apparent diffusion coefficient map, d planned dose distribution. Abbreviations: pCT planning CT, MRI Magnetic resonance imaging, T2w T2-weighted, ADC map apparent diffusion coefficient map

Dose accumulation

The workflow for generating a CBCT-based dose calculation and dose accumulation is described in the following section and illustrated in Fig. 2:

  1. 1.

    Import into the RayStation TPS: The original patient’s plan, which was generated with the Pinnacle3 TPS, as well as all CBCTs were transferred to the RayStation TPS. The Raystation TPS was chosen as primary TPS for this study because of RayStation`s deformable image registration algorithm, which was proven to be accurate and reliable in the literature [10, 11]. The target volumes, the organ at risk (OAR) contours and couch shift data from the kilovoltage X-ray volume imaging system (Elekta AB, Stockholm, Sweden) were imported into RayStation. The extent of the couch shift in the three main axes (left–right, superior–inferior, anterior–posterior) was analyzed.

  2. 2.

    Density override of CBCT/planning CT: A density override for each CBCT and planning CT (pCT) was performed within the RayStation TPS to reduce the influence of Houndsfield unit (HU) variability on the dose calculation. In the RayStation TPS, HU ranges of all CBCTs and planning CTs were used to select areas with the respective HU range before density override. Thereby, a density table with five density levels (air, lung, adipose tissue, cartilage/bone, and higher density for prosthesis) was created for all CBCTs and pCTs, resulting in a stepwise HU-to-density conversion curve [12]. The corrected pCT served as reference CT (CTref). The density override was applied to both the CBCTs as well as the pCTs. Therefore, the dose calculation error by density overriding is found in both calculated delivered and calculated planned dose in the RayStation TPS, which allows the comparison of both. The density overridden CBCTs and pCTs were utilized in the deformable image registration (DIR) process as described in the following sections.

  3. 3.

    Dose calculation on CTref: The patient`s plan was recalculated on the CTref with overrridden step-wise density in the RayStation TPS and served as the reference plan.

  4. 4.

    Rigid registration and dose calculation on CBCT: Next, each CBCT was rigidly co-registered to CTref according to the couch shift information from the integrated XVI system. This ensured that the actual treatment position of the patient with reference to the isocenter was taken into account. Afterwards, the dose distribution of the patient’s plan was calculated on each CBCT.

  5. 5.

    Deformable image registration (DIR): A DIR of the CBCT on CTref was performed which generated a deformation vector field by the appliance of the RayStation ANAtomically CONstrained Deformation Algorithm (ANACONDA)(10, 11). The ANACONDA algorithm measures image similarity by the appliance of a correlation coefficient. The DIR algorithm as described here offers a hybrid registration solution, combining both geometric and intensity information [10]. Because density overriding results in areas with homogeneous density, the algorithm mostly acts on organ boundaries. Target volumes did not change during treatment and DIR was not challenged by changing target volumes during treatment.

  6. 6.

    Dose mapping: The CBCT dose distribution was remapped back on CTref by utilizing the obtained deformation vector field. This resulted in the CBCT dose distribution on CTref, calculated out of the original plan on the RayStation TPS.

  7. 7.

    Dose accumulation: The CBCT dose distribution on the CTref, calculated out of the original plan on the RayStation TPS as described above, was used to calculate the fractional dose for each CBCT. In case that no CBCT was performed on a given day, the previous CBCT was selected for fractional dose calculation instead, and it was assumed that organ deformation stayed similar. By summation of all fractional doses a dose accumulation for PTVDIL, PTVBoost, PTV, and the organs at risk rectum, bladder, femoral heads, and penile bulb was achieved.

Fig. 2

Illustration of the dose computation workflow. After import into the RayStation TPS and application of a density override, the dose is computed on the density overridden pCT (CTref) (a) and on the density overridden CBCT (b). The dose computed on each CBCT is consecutively mapped via a deformation vector field (c), derived from deformable image registration, onto CTref (d). The vector field is color-coded from blue (0.0 cm deformation) to red (0.5 cm deformation). Abbreviation: CTref reference CT

Dominant intraprostatic lesion location and accumulated dose

The location of all index lesions was reviewed and classified according to the sector map of the American College of Radiology Prostate Imaging Reporting and Data System (PIRADS) Version 2.1, which is illustrated in Fig. 3 [13]. The accumulated dose for PTVDIL was analyzed for differences between ventral and dorsal index lesions. Dorsal location was defined as location inside the posterior medial peripheral zone (PZpm) of the prostate. Ventral location was defined as all other locations inside the prostate.

Fig. 3

PIRADS 2.1 sector map [14]. Permission for publication obtained from American College of Radiology Committee on PI-RADS®, licensed under Creative Commons Attribution-NoDerivatives 4.0 International (CC BY-ND 4.0). Abbreviations: AFS anterior fibromuscular stroma, CZ central zone, TZa/p transition zone, anterior/posterior, PZa/pm/pl peripheral zone, anterior/posterior medial/posterior lateral, US urethra


Matched-pairs t-tests were applied for the comparison of planned and accumulated doses for normally distributed parameters according to the Shapiro–Wilk test. In the case of non-normally distributed parameters, Wilcoxon matched-pairs signed-rank tests were applied instead. An unpaired t-test was used for the evaluation of differences in the accumulated dose according to index lesion location. Statistical significance was declared in the case of a two-sided p < 0.05. Statistical analysis was conducted using IBM SPSS v.26.0 (IBM Corp., Armonk, NY, USA). Descriptive quantitative values are expressed as mean ± standard deviation or as median with the corresponding range as appropriate.


Patient characteristics

The median patient age was 71 years (57–83 years) with a median Karnofsky Performance Status of 100% (70–100%). The mean initial prostate-specific antigen (PSA) value reached 9.2 ± 5.9 ng/mL and 5/13/5 patients had a Gleason score of 6/7/8–10, respectively. There were 5 cases of low-risk, 11 cases of intermediate-risk, and 7 cases of high-risk prostate cancer after the D`Amico risk classification [15]. The patient characteristics are summarized in Table 1.

Table 1 Patient characteristics

X-ray volume imaging couch shift analysis

Overall, 338 CBCTs were acquired with mean 14.7 ± 1.9 CBCTs per patient. The mean (± SD) couch shift in the three main axes were: 0.1 ± 2.9 mm (left–right), − 0.4 ± 2.7 mm (superior-inferior), − 0.0 ± 2.8 mm (anterior–posterior). Median shift for left–right was 0 cm (− 7 to 8 mm), median shift for superior-inferior was 0 cm (− 9 to 7 mm), median shift for anterior–posterior was 0 cm (− 9 to 8 mm). Regarding shift distribution, 95.8% of all shifts in the three axes were within 5 mm and 64% were within 2 mm. 0.8% of all shifts were greater than 7 mm.

Dose accumulation

Target coverage of PTVBoost and PTV was excellent and no cases of underdosage in the accumulated doses for PTVDIL, PTVBoost, and PTV in DMean, D95%, D02%, and D98% could be detected. For the organ at risk rectum, mean DMean was not significantly different between dose accumulation and initial treatment plans. Delivered bladder mean DMean was higher than planned with 19.4 ± 7.4 Gy versus 18.8 ± 7.5 Gy, p < 0.001. The penile bulb showed a slightly decreased delivered mean DMean with 29.1 ± 14.0 Gy versus 29.8 ± 14.4 Gy, p < 0.001. The target coverage and organ at risk parameters are summarized in Table 2.

Table 2 Target volume coverage and dose to organs at risk

Dominant intraprostatic lesion

The mean volume of the PTVDIL was 2.2 ± 2.1 cm3. 23 index lesions were identified. The sector distribution is summarized in Table 1. Eight out of 23 index lesions (34.8%) were located dorsally adjacent to the OAR rectum in the posterior medial peripheral zone (PZpm) according to the classification of PIRADS 2.1. Location of the index lesion in the PZpm had a statistically significant impact on the applied dose: DMean of the index lesions in the PZpm was significantly lower (76.5 ± 0.4 Gy, n = 8) than for index lesions in the rest of the prostate (78.7 ± 0.3 Gy, n = 15), p < 0.001, t-value = 4.921 (21 df). The mean DMean difference was 2.2 Gy (95% CI 1.3–3.1 Gy). The effect strength according to Cohen was 0.73, which represents a strong effect (Fig. 4).

Fig. 4

Effect of DIL location on the applied dose. Shown is the boxplot of the delivered mean dose in Gy for PTVDIL for ventral and dorsal index lesions. Dorsal index lesions, defined as index lesions in the posterior medial peripheral zone of the prostate, showed a significantly lower mean DMean (p < 0.001)


Dose-guided adaptive radiotherapy for external beam radiotherapy is increasingly used to adjust to anatomy changes and to avoid underdosage of the target volumes due to motion and organ deformation [16]. Dose accumulation and the CBCT-derived fractional daily dose may offer a feedback mechanism to trigger the adaption of radiotherapy. For focal prostate radiotherapy, the critical target structure is the index lesion and to our knowledge, this is the first study to investigate the dose coverage of the DIL in prostate cancer by CBCT-based dose accumulation.

Dose accumulation for prostate radiotherapy in general has been evaluated before: In a recent analysis, Ong et al. presented a workflow for dose accumulation of prostate and pelvic lymph node radiotherapy with two isocenters. Due to different coordinate systems in the utilized MIM v.6.9 (MIMVista Corp., Cleveland OH) and MosaiQ (Elekta AB, Stockholm, Sweden) software Ong et al. had to perform phantom measurements for couch shift commissioning which allowed the import of MosaiQ extracted, correctly aligned shift information into the contouring software [17]. In our current analysis, the RayStation TPS allowed the direct import and usage of couch shift data and therefore phantom measurements were omitted. The analysis of shift distribution demonstrated significantly, but non-excessive couch shifts. Over 95% of all shifts were adequately covered within 5 mm and therefore within the high dose PTV (PTVBoost). 4.2% and 0.8% of all shifts exceeded 5 mm, respective 7 mm. This underlines the importance of inter-fractional dose accumulation and also confirms the utilized institutional margin concept. The delivered dose demonstrated excellent target coverage of the target volumes and the dominant intraprostatic lesion. No cases of underdosage for the PTV, PTVBoost as well as PTVDIL could be detected.

The delivered mean dose was significantly lower for dorsal DILs which may be attributable to the dose fall-off at the edges of PTVBoost as a steep dose gradient for sparing the OAR rectum was utilized. This is in line with an earlier work that showed decreased planned dose for dorsal index lesions in prostate radiotherapy [8]. Despite no margin between the rectum and PTVBoost, a long-term analysis of the utilized regime of moderately hypofractionated prostate radiotherapy with the institutional target volume and margin concept showed excellent biochemical relapse-free survival and low rates of gastrointestinal toxicity [9]. In the current work, dose accumulation showed excellent agreement between planned and delivered rectum mean DMean. The delivered dose to the bladder was higher than planned which emphasizes the importance of strict adherence to bladder filling protocols to avoid genitourinary toxicity. In a future analysis, the presented workflow should be explored further for the correlation between delivered dose and clinical toxicity as delivered dose has been prognostic for unwanted sequelae in the literature [18]. In case of severe CBCT artifacts in the target area our approach may not be possible. CBCT artifacts in the target area which need to be mitigated for our approach are especially scatter artifacts, extinction artifacts, and beam hardening artifacts [19]. Scatter artifacts reduce soft-tissue contrast and will affect the density values of the tissues. For highly absorbing material in the target area, the signal recorded in the detector behind the material may be close to zero, resulting in extinction artifacts. The absorption of lower energetic rays, for example by titanium in hip total endoprosthesis, will cause massive beam hardening artifacts for pelvic CBCTs [19]. In our study cohort, we did not observe severe CBCT artifacts. A quality assurance procedure for the presented workflow would have to incorporate the verification of the density override and DIR algorithm accuracy. A possible solution would be to use another dedicated software or treatment planning system to compare the results with the RayStation generated ones. To quantify the differences, a metric would have to be developed which would be used for defining action thresholds. The development of a quality assurance procedure for the dose accumulation workflow was beyond the scope of this study.

Our study has limitations: The accuracy of CBCT derived dose is determined by relative electron density and therefore correlated with HU. CBCT quality has a major influence on density calculation and therefore may have introduced uncertainty in the dose accumulation analysis. As detailed in the methods section, although a density assignment as correction strategy was applied, the calculated absolute dose values suffer from uncertainties due to step-wise instead of continuous conversion curves [12]. Dunlop et al. showed an absolute mean dose error of 0.7% for pelvic CBCT-based step-wise dose calculation in case of an anterior–posterior distance < 25 cm [12]. In our study, we estimated the error of our approach by comparing the mean dose difference between pCT and CTref, which amounted to 0.7% for PTVBoost. DIR quality also influences the robustness of the dose accumulation, but the ANACONDA DIR algorithm was shown to have high accuracy in the literature [11]. Another limitation is the lack of daily CBCTs according to the institutional imaging protocol which foresees daily CBCT only in the beginning when the patient is not accustomed to the organ preparation protocols. Later on, when the patient is trained, organ filling states remain mostly stable in clinical practice. We recently published our data for definitive prostate radiotherapy using this IGRT protocol showing no detrimental effect on biochemical relapse-free survival [9]. As strict adherence to organ preparation protocols is still mandatory, selecting the closest temporal CBCT in case of missing CBCTs is expected to have no excessive impact on the dose accumulation accuracy. Nonetheless, daily instead of non-daily CBCT may further improve the dose accumulation workflow and will be explored in sequential studies [20, 21]. Our data is also limited by the absence of intra-fractional motion analysis and correction, which was not available for the analyzed patient cohort. In particular for stereotactic radiotherapy, but as well as moderately fractionated radiotherapy, target motion during delivery can be significant and influence delivered dose. Using VMAT instead of step & shoot IMRT, the treatment time per fraction is reduced, which translates to a decreased risk of intra-fractional organ changes [22]. As next step, incorporating intra-fractional prostate motion into the dose accumulation workflow will be investigated. Lastly, the retrospective character and small patient sample of the current analysis warrants confirmation in a bigger patient cohort.

To look further, artificial intelligence may improve CBCT correction and dose calculation as has been highlighted recently: Maspero et al. trained a single neural network to generate synthetic CTs (sCT) out of CBCT for head-and-neck, lung, and breast cancer, enabling sCT-based dose calculation with under 0.5% mean dose difference to re-planning CTs which is lower than the estimated error of 0.7% in our study, although a direct comparison was outside the scope of our analysis [23]. Deep learning-based dose calculation, together with automatic segmentation, may represent crucial steps on the way to fully automatic dose-guided adaptive radiotherapy in general and especially for focal boost radiotherapy with tight margins. Regarding dose-escalation to the index lesion, the FLAME trial recently showed that a high dose focal boost improves biochemical disease-free survival in intermediate- and high-risk localized prostate cancer without additional toxicity [6]. The correlation between delivered dose to the index lesion, derived by dose accumulation, and clinical outcome remains an interesting research topic.


A workflow for the calculation of the delivered dose to the dominant intraprostatic lesion using CBCT-based dose accumulation on the RayStation TPS was established. The delivered dose to the dominant intraprostatic lesion and organs at risk was adequate and well within planning limits. Dorsal index lesions demonstrated a lower applied mean dose than ventral lesions, but no cases of underdosage were observable.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.


ADC map:

Apparent diffusion coefficient map


Anterior fibromuscular stroma


Cone beam computed tomography


Confidence interval


Computed tomography

CTref :

Reference computed tomography


Clinical target volume


Central zone


Degrees of freedom

DMean :

Mean dose


Dominant intraprostatic lesion


Deformable image registration

DX% :

Dose to X% of the target volume


Hounsfield units


Image-guided radiotherapy


Karnofsky Performance Status


Magnetic resonance imaging


Organs at risk


Planning CT


Prostate Imaging Reporting and Data System


Prostate-specific antigen


Planning target volume


Peripheral zone, anterior/posterior medial/posterior lateral


Synthetic CT


Simultaneous integrated boost




Treatment planning system


Transition zone, anterior/posterior


Volumetric arc radiotherapy




X-ray volume imaging


  1. 1.

    Rawla P. Epidemiology of prostate cancer. World J Oncol. 2019;10(2):63–89.

    CAS  Article  Google Scholar 

  2. 2.

    Quivrin M, Loffroy R, Cormier L, Mazoyer F, Bertaut A, Chambade D, et al. Multiparametric MRI and post implant CT-based dosimetry after prostate brachytherapy with iodine seeds: the higher the dose to the dominant index lesion, the lower the PSA bounce. Radiother Oncol. 2015;117(2):258–61.

    Article  Google Scholar 

  3. 3.

    Cellini N, Morganti AG, Mattiucci GC, Valentini V, Leone M, Luzi S, et al. Analysis of intraprostatic failures in patients treated with hormonal therapy and radiotherapy: implications for conformal therapy planning. Int J Radiat Oncol Biol Phys. 2002;53(3):595–9.

    Article  Google Scholar 

  4. 4.

    Pucar D, Hricak H, Shukla-Dave A, Kuroiwa K, Drobnjak M, Eastham J, et al. Clinically significant prostate cancer local recurrence after radiation therapy occurs at the site of primary tumor: magnetic resonance imaging and step-section pathology evidence. Int J Radiat Oncol Biol Phys. 2007;69(1):62–9.

    Article  Google Scholar 

  5. 5.

    Arrayeh E, Westphalen AC, Kurhanewicz J, Roach M 3rd, Jung AJ, Carroll PR, et al. Does local recurrence of prostate cancer after radiation therapy occur at the site of primary tumor? Results of a longitudinal MRI and MRSI study. Int J Radiat Oncol Biol Phys. 2012;82(5):e787–93.

    Article  Google Scholar 

  6. 6.

    Kerkmeijer LGW, Groen VH, Pos FJ, Haustermans K, Monninkhof EM, Smeenk RJ, et al. Focal boost to the intraprostatic tumor in external beam radiotherapy for patients with localized prostate cancer: results from the FLAME randomized phase III trial. J Clin Oncol. 2021;39(7):787–96.

    Article  Google Scholar 

  7. 7.

    Tocco BR, Kishan AU, Ma TM, Kerkmeijer LGW, Tree AC. MR-guided radiotherapy for prostate cancer. Front Oncol. 2020;10:616291.

    Article  Google Scholar 

  8. 8.

    Tamihardja J, Zenk M, Flentje M. MRI-guided localization of the dominant intraprostatic lesion and dose analysis of volumetric modulated arc therapy planning for prostate cancer. Strahlenther Onkol. 2019;195(2):145–52.

    Article  Google Scholar 

  9. 9.

    Tamihardja J, Schortmann M, Lawrenz I, Weick S, Bratengeier K, Flentje M, et al. Moderately hypofractionated radiotherapy for localized prostate cancer: updated long-term outcome and toxicity analysis. Strahlenther Onkol. 2021;197(2):124–32.

    Article  Google Scholar 

  10. 10.

    Weistrand O, Svensson S. The ANACONDA algorithm for deformable image registration in radiotherapy. Med Phys. 2015;42(1):40–53.

    Article  Google Scholar 

  11. 11.

    Takayama Y, Kadoya N, Yamamoto T, Ito K, Chiba M, Fujiwara K, et al. Evaluation of the performance of deformable image registration between planning CT and CBCT images for the pelvic region: comparison between hybrid and intensity-based DIR. J Radiat Res. 2017;58(4):567–71.

    Article  Google Scholar 

  12. 12.

    Dunlop A, McQuaid D, Nill S, Murray J, Poludniowski G, Hansen VN, et al. Comparison of CT number calibration techniques for CBCT-based dose calculation. Strahlenther Onkol. 2015;191(12):970–8.

    Article  Google Scholar 

  13. 13.

    Turkbey B, Rosenkrantz AB, Haider MA, Padhani AR, Villeirs G, Macura KJ, et al. Prostate Imaging Reporting and Data System Version 2.1: 2019 update of Prostate Imaging Reporting and Data System Version 2. Eur Urol. 2019;76(3):340–51.

    Article  Google Scholar 

  14. 14.

    American College of Radiology Committee on PI-RADS®. PI-RADS 2.1 sector map. Accessed 22 June 2021.

  15. 15.

    D’Amico AV, Whittington R, Malkowicz SB, Schultz D, Blank K, Broderick GA, et al. Biochemical outcome after radical prostatectomy, external beam radiation therapy, or interstitial radiation therapy for clinically localized prostate cancer. JAMA. 1998;280(11):969–74.

    CAS  Article  Google Scholar 

  16. 16.

    McPartlin AJ, Li XA, Kershaw LE, Heide U, Kerkmeijer L, Lawton C, et al. MRI-guided prostate adaptive radiotherapy—a systematic review. Radiother Oncol. 2016;119(3):371–80.

    CAS  Article  Google Scholar 

  17. 17.

    Ong A, Knight K, Panettieri V, Dimmock M, Tuan JKL, Tan HQ, et al. Development of an automated radiotherapy dose accumulation workflow for locally advanced high-risk prostate cancer—a technical report. J Med Radiat Sci. 2020;6:66.

    Google Scholar 

  18. 18.

    Shelley LEA, Scaife JE, Romanchikova M, Harrison K, Forman JR, Bates AM, et al. Delivered dose can be a better predictor of rectal toxicity than planned dose in prostate radiotherapy. Radiother Oncol. 2017;123(3):466–71.

    CAS  Article  Google Scholar 

  19. 19.

    Schulze R, Heil U, Gross D, Bruellmann DD, Dranischnikow E, Schwanecke U, et al. Artefacts in CBCT: a review. Dentomaxillofac Radiol. 2011;40(5):265–73.

  20. 20.

    Splinter M, Sachpazidis I, Bostel T, Fechter T, Zamboglou C, Thieke C, et al. Dosimetric impact of the positional imaging frequency for hypofractionated prostate radiotherapy—a voxel-by-voxel analysis. Front Oncol. 2020;10:564068.

    Article  Google Scholar 

  21. 21.

    Bostel T, Sachpazidis I, Splinter M, Bougatf N, Fechter T, Zamboglou C, et al. Dosimetric impact of interfractional variations in prostate cancer radiotherapy-implications for imaging frequency and treatment adaptation. Front Oncol. 2019;9:940.

    Article  Google Scholar 

  22. 22.

    Ballhausen H, Li M, Ganswindt U, Belka C. Shorter treatment times reduce the impact of intra-fractional motion: A real-time 4DUS study comparing VMAT vs. step-and-shoot IMRT for prostate cancer. Strahlenther Onkol. 2018;194(7):664–74.

  23. 23.

    Maspero M, Houweling AC, Savenije MHF, van Heijst TCF, Verhoeff JJC, Kotte A, et al. A single neural network for cone-beam computed tomography-based radiotherapy of head-and-neck, lung and breast cancer. Phys Imaging Radiat Oncol. 2020;14:24–31.

    Article  Google Scholar 

Download references


Not applicable.


Open Access funding enabled and organized by Projekt DEAL.

Author information




Conception and design of the study: JT, MF. Data collection: JT, SC, PK, GR, FE, AR. Data analysis and interpretation: JT, SC, PK. Drafting the article: JT, PK. Critical revision of the article: JT, SC, PK, GR, FE, AR, BP, MF. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jörg Tamihardja.

Ethics declarations

Ethics approval and consent to participate

Ethical review and written informed consent were not required for participation in this retrospective analysis in accordance with the local legislation (BayKrG Art. 27 (4)) and institutional requirements. Written informed consent for treatment was provided by all patients.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Tamihardja, J., Cirsi, S., Kessler, P. et al. Cone beam CT-based dose accumulation and analysis of delivered dose to the dominant intraprostatic lesion in primary radiotherapy of prostate cancer. Radiat Oncol 16, 205 (2021).

Download citation


  • Adaptive radiotherapy
  • Deformable image registration
  • Dominant intraprostatic lesion
  • Dose accumulation
  • Prostate cancer
  • Prostate Imaging Reporting and Data System