Evaluation of accuracy of B-spline transformation-based deformable image registration with different parameter settings for thoracic images

Deformable image registration (DIR) is fundamental technique for adaptive radiotherapy and image-guided radiotherapy. However, further improvement of DIR is still needed. We evaluated the accuracy of B-spline transformation-based DIR implemented in elastix. This registration package is largely based on the Insight Segmentation and Registration Toolkit (ITK), and several new functions were implemented to achieve high DIR accuracy. The purpose of this study was to clarify whether new functions implemented in elastix are useful for improving DIR accuracy. Thoracic 4D computed tomography images of ten patients with esophageal or lung cancer were studied. Datasets for these patients were provided by DIR-lab (dir-lab.com) and included a coordinate list of anatomical landmarks that had been manually identified. DIR between peak-inhale and peak-exhale images was performed with four types of parameter settings. The first one represents original ITK (Parameter 1). The second employs the new function of elastix (Parameter 2), and the third was created to verify whether new functions improve DIR accuracy while keeping computational time (Parameter 3). The last one partially employs a new function (Parameter 4). Registration errors for these parameter settings were calculated using the manually determined landmark pairs. 3D registration errors with standard deviation over all cases were 1.78 (1.57), 1.28 (1.10), 1.44 (1.09) and 1.36 (1.35) mm for Parameter 1, 2, 3 and 4, respectively, indicating that the new functions are useful for improving DIR accuracy, even while maintaining the computational time, and this B-spline-based DIR could be used clinically to achieve high-accuracy adaptive radiotherapy.     

Development of fast patient position verification software using 2D-3D image registration and its clinical experience

To improve treatment workflow, we developed a graphic processing unit (GPU)-based patient positional verification software application and integrated it into carbon-ion scanning beam treatment. Here, we evaluated the basic performance of the software. The algorithm provides 2D/3D registration matching using CT and orthogonal X-ray flat panel detector (FPD) images. The participants were 53 patients with tumors of the head and neck, prostate or lung receiving carbon-ion beam treatment. 2D/3D-ITchi-Gime (ITG) calculation accuracy was evaluated in terms of computation time and registration accuracy. Registration calculation was determined using the similarity measurement metrics gradient difference (GD), normalized mutual information (NMI), zero-mean normalized cross-correlation (ZNCC), and their combination. Registration accuracy was dependent on the particular metric used. Representative examples were determined to have target registration error (TRE) = 0.45 ± 0.23 mm and angular error (AE) = 0.35 ± 0.18° with ZNCC + GD for a head and neck tumor; TRE = 0.12 ± 0.07 mm and AE = 0.16 ± 0.07° with ZNCC for a pelvic tumor; and TRE = 1.19 ± 0.78 mm and AE = 0.83 ± 0.61° with ZNCC for lung tumor. Calculation time was less than 7.26 s.The new registration software has been successfully installed and implemented in our treatment process. We expect that it will improve both treatment workflow and treatment accuracy.     

A slice classification model-facilitated 3D encoder–decoder network for segmenting organs at risk in head and neck cancer

For deep learning networks used to segment organs at risk (OARs) in head and neck (H&N) cancers, the class-imbalance problem between small volume OARs and whole computed tomography (CT) images results in delineation with serious false-positives on irrelevant slices and unnecessary time-consuming calculations. To alleviate this problem, a slice classification model-facilitated 3D encoder–decoder network was developed and validated. In the developed two-step segmentation model, a slice classification model was firstly utilized to classify CT slices into six categories in the craniocaudal direction. Then the target categories for different OARs were pushed to the different 3D encoder–decoder segmentation networks, respectively. All the patients were divided into training (n = 120), validation (n = 30) and testing (n = 20) datasets. The average accuracy of the slice classification model was 95.99%. The Dice similarity coefficient and 95% Hausdorff distance, respectively, for each OAR were as follows: right eye (0.88 ± 0.03 and 1.57 ± 0.92 mm), left eye (0.89 ± 0.03 and 1.35 ± 0.43 mm), right optic nerve (0.72 ± 0.09 and 1.79 ± 1.01 mm), left optic nerve (0.73 ± 0.09 and 1.60 ± 0.71 mm), brainstem (0.87 ± 0.04 and 2.28 ± 0.99 mm), right temporal lobe (0.81 ± 0.12 and 3.28 ± 2.27 mm), left temporal lobe (0.82 ± 0.09 and 3.73 ± 2.08 mm), right temporomandibular joint (0.70 ± 0.13 and 1.79 ± 0.79 mm), left temporomandibular joint (0.70 ± 0.16 and 1.98 ± 1.48 mm), mandible (0.89 ± 0.02 and 1.66 ± 0.51 mm), right parotid (0.77 ± 0.07 and 7.30 ± 4.19 mm) and left parotid (0.71 ± 0.12 and 8.41 ± 4.84 mm). The total segmentation time was 40.13 s. The 3D encoder–decoder network facilitated by the slice classification model demonstrated superior performance in accuracy and efficiency in segmenting OARs in H&N CT images. This may significantly reduce the workload for radiation oncologists.     

Fully automated rigid image registration versus human registration in postoperative spine stereotactic body radiation therapy: a multicenter non-inferiority study

To confirm the fully automated rigid image registration (A-RIR) accuracy in postoperative spine stereotactic body radiation therapy (SBRT), we conducted a multicenter non-inferiority study compared to the human rigid image registration (H-RIR). Twenty-eight metastatic cancer patients who underwent postoperative spine SBRT are enrolled—image registration (IR) of planning computed tomography (CT) and CT-myelogram for delineating the spinal cord. The adopted A-RIR workflow is a contour-focused algorithm performing a rigid registration by maximizing normalized mutual information (NMI) restricted to the data contained within the automatically extracted contour. Three radiation oncologists (ROs) from multicenters were prompted to review two blinded registrations and choose one for clinical use. Indistinguishable cases were allowed to vote equivalent, counted A-RIR side. A-RIR is considered non-inferior to H-RIR if the lower limit of the 95% confidence interval (CI) of A-RIR preferable/equivalent is greater than 0.45. We also evaluated the NMI improvement from the baseline and the translational/rotational errors between A-RIR and H-RIR. The A-RIR preferable/equivalent was selected in 21 patients (0.75, 95% CI: 0.55–0.89), demonstrating non-inferiority to H-RIR. The A-RIR’s NMI improvement was greater than H-RIR in 24 patients: the mean value ± SD was 0.225 ± 0.115 in A-RIR and 0.196 ± 0.114 in H-RIR (P < 0.001). The absolute translational error was 0.38 ± 0.31 mm. The rotational error was −0.03 ± 0.20, 0.05 ± 0.19, −0.04 ± 0.20 degrees in axial, coronal, and sagittal planes (range: −0.66–0.52). In conclusion, A-RIR shows non-inferior to H-RIR in CT and CT-myelogram registration for postoperative spine SBRT planning.     

Dosimetric and geometric evaluation of the use of deformable image registration in adaptive intensity-modulated radiotherapy for head-and-neck cancer

The aim of this study was to carry out geometric and dosimetric evaluation of the usefulness of a deformable image registration algorithm utilized for adaptive head-and-neck intensity-modulated radiotherapy. Data consisted of seven patients, each with a planning CT (pCT), a rescanning CT (ReCT) and a cone beam CT (CBCT). The CBCT was acquired on the same day (±1 d) as the ReCT (i.e. at Fraction 17, 18, 23, 24 or 29). The ReCT served as ground truth. A deformed CT (dCT) with structures was created by deforming the pCT to the CBCT. The geometrical comparison was based on the volumes of the deformed, and the manually delineated structures on the ReCT. Likewise, the center of mass shift (CMS) and the Dice similarity coefficient were determined. The dosimetric comparison was performed by recalculating the initial treatment plan on the dCT and the ReCT. Dose–volume histogram (DVH) points and a range of conformity measures were used for the evaluation. We found a significant difference in the median volume of the dCT relative to that of the ReCT. Median CMS values were ∼2–5 mm, except for the spinal cord, where the median CMS was 8 mm. Dosimetric evaluation of target structures revealed small differences, while larger differences were observed for organs at risk. The deformed structures cannot fully replace manually delineated structures. Based on both geometrical and dosimetrical measures, there is a tendency for the dCT to overestimate the need for replanning, compared with the ReCT.     

Estimation of the total rectal dose of radical external beam and intracavitary radiotherapy for uterine cervical cancer using the deformable image registration method

We adapted the deformable image registration (DIR) technique to accurately calculate the cumulative intracavitary brachytherapy (ICBT) and external beam radiotherapy (EBRT) rectal dose for treating uterine cervical cancer. A total of 14 patients with primary cervical cancer radically treated with ICRT and EBRT were analysed using the Velocity AI^TM software. Computed tomography (CT) images were registered, and EBRT and ICBT dose distributions were determined. Cumulative D_2cm³, D_1cm³ and D_0.1cm³ were calculated by simple addition of fractional values or by DIR. The accuracy of DIR was evaluated by means of a virtual phantom mimicking the rectum. The dice similarity coefficient (DSC) was calculated to evaluate rectal contour concordance between CT images before and after DIR. Virtual phantom analysis revealed that the average difference between the DIR-based phantom D_mean and the simple phantom D_mean was 1.9 ± 2.5 Gy (EQD₂), and the DIR method included an uncertainty of ∼8.0%. The mean DSC between reference CT and CT was significantly improved after DIR (EBRT: 0.43 vs 0.85, P < 0.005; ICBT: 0.60 vs 0.87, P < 0.005). The average simple rectal D_2cm³, D_1cm³ and D_0.1cm³ values were 77.6, 81.6 and 91.1 Gy (EQD₂), respectively; the DIR-based values were 76.2, 79.5 and 87.6 Gy, respectively. The simple addition values were overestimated, on average, by 3.1, 3.7 and 5.5 Gy, respectively, relative to the DIR-based values. In conclusion, the difference between the simple rectal dose–volume histogram (DVH) parameter addition and DIR-based cumulative rectal doses increased with decreasing DVH parameters.     

Comparison of predictive performance for toxicity by accumulative dose of DVH parameter addition and DIR addition for cervical cancer patients

We compared predictive performance between dose volume histogram (DVH) parameter addition and deformable image registration (DIR) addition for gastrointestinal (GI) toxicity in cervical cancer patients. A total of 59 patients receiving brachytherapy and external beam radiotherapy were analyzed retrospectively. The accumulative dose was calculated by three methods: conventional DVH parameter addition, full DIR addition and partial DIR addition. , and (minimum doses to the most exposed 2 cm³, 1cm³ and 0.1 cm³ of tissue, respectively) of the rectum and sigmoid were calculated by each method. V50, V60 and V70 Gy (volume irradiated over 50, 60 and 70 Gy, respectively) were calculated in full DIR addition. The DVH parameters were compared between toxicity (≥grade1) and non-toxicity groups. The area under the curve (AUC) of the receiver operating characteristic (ROC) curves were compared to evaluate the predictive performance of each method. The differences between toxicity and non-toxicity groups in were 0.2, 5.7 and 3.1 Gy for the DVH parameter addition, full DIR addition and partial DIR addition, respectively. The AUCs of were 0.51, 0.67 and 0.57 for DVH parameter addition, full DIR addition and partial DIR addition, respectively. In full DIR addition, the difference in dose between toxicity and non-toxicity was the largest and AUC was the highest. AUCs of V50, V60 and V70 Gy were 0.51, 0.63 and 0.62, respectively, and V60 and V70 were high values close to the value of of the full DIR addition. Our results suggested that the full DIR addition may have the potential to predict toxicity more accurately than the conventional DVH parameter addition, and that it could be more effective to accumulate to all pelvic irradiation by DIR.     

Assessment of anatomical and dosimetric changes by a deformable registration method during the course of intensity-modulated radiotherapy for nasopharyngeal carcinoma

The aim of this study was to quantify the anatomic variations and the dosimetric effects accessed by a deformable registration method throughout the entire course of radiotherapy, and to evaluate the necessity of replanning for patients with nasopharyngeal carcinoma (NPC). Plan1(CT1) was based on the original CT, and Plan2(CT2) was generated from the midtreatment CT scan acquired after 25 fractions of IMRT of Plan1. Both sets of CTs, RT structures and RT doses for the two group plans were transferred to a workstation, and then a hybrid IMRT plan, Plan1(CT2), was generated by deforming doses of Plan1 to CT2. Subsequently, the accumulated plan, Plan1 + 2(CT2), was generated to quantify the actual dosimetric effects during the course. The transverse diameter of the neck at the center of the odontoid process was (15.4 ± 1.0) cm and (14.4 ± 1.1) cm in CT1 and CT2, respectively (P < 0.05). Compared with CT1, the mean volumes of the right and left parotid glands were significantly decreased by (24.6 ± 11.9)% and (35.1 ± 20.1)%, respectively. Comparison of Plan1 (CT1) with Plan1 (CT2) indicated that the doses to targets decreased without replanning. With repeated CT and replanning after 25 fractions, the doses to targets would be improved. The doses to normal tissue were increased without replanning. For eight patients out of 12, the doses to the spinal cord and brainstem exceeded the constraints without replanning, while the corresponding values decreased with replanning. During the entire course of IMRT, the volumes of the targets and the parotid glands would be reduced significantly. Midtreatment CT scanning and replanning are recommended to ensure adaptive doses to the targets and critical normal tissues.     

Analysis on the accuracy of CT-guided radioactive I-125 seed implantation with 3D printing template assistance in the treatment of thoracic malignant tumors

This article analyzes the accuracy of needle track and dose of a 3-dimensional printing template (3DPT) in the treatment of thoracic tumor with radioactive I-125 seed implantation (RISI). A total of 28 patients were included. The technical process included: (i) preoperative CT positioning, (ii) preoperative planning design, (iii) 3DPT design and printing, (iv) 3DPT alignment, (v) puncture and seed implantation. The errors of needle position and dosimetric parameters were analyzed. A total of 318 needles were used. The mean errors in needle depth, needle insertion point, needle tip and needle angle were 0.52 ± 0.48 cm, 3.4 ± 1.7 mm, 4.4 ± 2.9 mm and 2.8 ± 1.7°, respectively. The differences between actual needle insertion angle and needle depth and those designed in the preoperative were statistically significant (p < 0.05). The mean values of all the errors of the chest wall cases were smaller than those of the lungs, and the differences were statistically significant (p < 0.05). There was no significant difference between the D90 calculated in the postoperative plan and those designed in the preoperative and intraoperative plans (p > 0.05). Some dosimetric parameters of preoperative plans such as V100, V200, CI and HI were not consistent with that of preoperative plans, and the difference was statistically significant (p < 0.05). However, there were no statistical difference in the dosimetric parameters between the postoperative plans and intraoperative plans (p > 0.05). We conclude that for thoracic tumors, even under the guidance of 3DPT, there will be errors. The plan should be optimized in real time during the operation.     

Development of a novel airbag system of abdominal compression for reducing respiratory motion: preliminary results in healthy volunteers

This study used cine-magnetic resonance imaging (cine-MRI) to evaluate the safety and efficacy of a novel airbag system combined with a shell-type body fixation system in reducing respiratory motion in normal volunteers. The airbag system consists of a six-sided polygon inflatable airbag, a same shape plate, a stiff air supply tube, an air-supply pump and a digital pressure load cell monitor. Piezoelectric sensors were installed in the plate to detect compression pressure load changes; pressure load data were transferred to the digital pressure load cell monitor through Bluetooth. Five volunteers underwent cine-MRI with and without airbag compression to detect differences in the respiratory motion of the organs. The volunteers’ physiologic signs were stable during the experiment. The maximum inspiration pressure load was 4.48 ± 0.86 kgf (range, 4.00–6.00 kgf), while the minimum expiration pressure load was 3.69 ± 0.95 kgf (range, 2.8–5.3 kgf). Under airbag compression, the right diaphragm movement was reduced from 19.50 ± 6.43 mm to 9.60 ± 3.61 mm (P < 0.05) in the coronal plane and 23.12 ± 6.30 mm to 11.00 ± 3.69 mm (P < 0.05) in the sagittal plane. The left diaphragm, pancreas and liver in the coronal plane and the right kidney and liver in the sagittal plane also showed significant movement reduction. This novel airbag abdominal compression system was found to be safe during the experiment and successful in the reduction of internal organ respiratory motion and promises to be a convenient and efficient tool for clinical radiotherapy.     

First clinical experience in carbon ion scanning beam therapy: retrospective analysis of patient positional accuracy

Our institute has constructed a new treatment facility for carbon ion scanning beam therapy. The first clinical trials were successfully completed at the end of November 2011. To evaluate patient setup accuracy, positional errors between the reference Computed Tomography (CT) scan and final patient setup images were calculated using 2D-3D registration software. Eleven patients with tumors of the head and neck, prostate and pelvis receiving carbon ion scanning beam treatment participated. The patient setup process takes orthogonal X-ray flat panel detector (FPD) images and the therapists adjust the patient table position in six degrees of freedom to register the reference position by manual or auto- (or both) registration functions. We calculated residual positional errors with the 2D-3D auto-registration function using the final patient setup orthogonal FPD images and treatment planning CT data. Residual error averaged over all patients in each fraction decreased from the initial to the last treatment fraction [1.09 mm/0.76° (averaged in the 1st and 2nd fractions) to 0.77 mm/0.61° (averaged in the 15th and 16th fractions)]. 2D-3D registration calculation time was 8.0 s on average throughout the treatment course. Residual errors in translation and rotation averaged over all patients as a function of date decreased with the passage of time (1.6 mm/1.2° in May 2011 to 0.4 mm/0.2° in December 2011). This retrospective residual positional error analysis shows that the accuracy of patient setup during the first clinical trials of carbon ion beam scanning therapy was good and improved with increasing therapist experience.     

Intrafractional setup errors in patients undergoing non-invasive fixation using an immobilization system during hypofractionated stereotactic radiotherapy for lung tumors

Intrafractional setup errors during hypofractionated stereotactic radiotherapy (SRT) were investigated on the patient under voluntary breath-holding conditions with non-invasive immobilization on the CT-linac treatment table. A total of 30 patients with primary and metastatic lung tumors were treated with the hypofractionated SRT with a total dose of 48–60 Gy with four treatment fractions. The patient was placed supine and stabilized on the table with non-invasive patient fixation. Intrafractional setup errors in Right/Left (R.L.), Posterior/Anterior (P.A.), and Inferior/Superior (I.S.) dimensions were analyzed with pre- and post-irradiation CT images. The means and one standard deviation of the intrafractional errors were 0.9 ± 0.7mm (R.L.), 0.9 ± 0.7mm (P.A.) and 0.5 ± 1.0 mm (I.S.). Setup errors in each session of the treatment demonstrated no statistically significant difference in the mean value between any two sessions. The frequency within 3mm displacement was 98% in R.L., 98% in P.A. and 97% in I.S. directions. SRT under the non-invasive patient fixation immobilization system with a comparatively loose vacuum pillow demonstrated satisfactory reproducibility of minimal setup errors with voluntary breath-holding conditions that required a small internal margin.     

High-dose-rate interstitial brachytherapy for gynecologic malignancies—dosimetric changes during treatment period

To overcome cranio-caudal needle displacement in pelvic high-dose-rate interstitial brachytherapy (HDRIB), we have been utilizing a fullystretched elastic tape to thrust the template into the perineum. The purpose of the current study was to evaluate dosimetric changes during the treatment period using this thrusting method, and to explore reproducible planning methods based on the results of the dosimetric changes. Twenty-nine patients with gynecologic malignancies were treated with HDRIB at the Cancer Institute Hospital. Pre-treatment and post-treatment computed tomography (CT) scans were acquired and a virtual plan for post-treatment CT was produced by applying the dwell positions/times of the original plan. For the post-treatment plan, D90 for the clinical target volume (CTV) and D2cc for the rectum and bladder were assessed and compared with that for the original plan. Cranio-caudal needle displacement relative to CTV during treatment period was only 0.7 ± 1.9 mm. The mean D90 values for the CTV in the pre- and post-treatment plans were stable (6.8 Gy vs. 6.8 Gy) and the post-treatment/pre-treatment D90 ratio was 1.00 ± 0.08. The post-/pre-treatment D2cc ratio was 1.14 ± 0.22 and the mean D2cc for the rectum increased for the post-treatment plan (5.4 Gy vs. 6.1 Gy), especially when parametrial infiltration was present. The mean D2cc for the bladder was stable (6.3 Gy vs. 6.6 Gy) and the ratio was 1.06 ± 0.20. Our thrusting method achieved a stable D90 for the CTV, in contrast to previous prostate HDRIB reports displaying reductions of 35–40% for D90 during the treatment period.     

Use of lung treatment plans to evaluate DIR algorithms

The purpose of the study is to evaluate the accuracy of two deformable image registration algorithms by examining their influence on the dose summation results obtained using 4DCT (four dimensional computed tomography) dose distributions based on ‘4D’ planned and ‘4D optimal’ IMRT (intensity modulated radiation therapy) plans. For ten lung cancer patients, 4D step and shoot IMRT plans were produced. The breathing cycle was divided into ten parts and for each part a set of CT images was acquired. For each patient the treatment plan was copied to the CTs of each phase and subsequently recalculated. Each phase CT was then registered to the average intensity projection (AIP) CT using a deformable image registration (DIR) algorithm and the composite dose distribution was then calculated by summing up the deformed dose distributions from all the phases (‘4D’ treatment plan). The ‘4D optimal’ treatment plan was created by producing an optimal plan on the CTs of each phase of the respiratory cycle and summing up the deformed dose distributions from all the phases. The results indicate that it is possible to map the dose distributions of different breathing phases in lung using DIR, and that different DIR methods and target characteristics (motion amplitude, size, location) affect the differences between original plan, ‘4D’ and ‘4D optimal’ dose distributions. Although the ‘4D optimal’ plans were designed to achieve 95% target coverage, both of the used DIR methods failed to translate that coverage in some instances. The same variation between these methods was also observed in the ‘4D’ plan comparison. This study shows that it is feasible to perform an acceptably accurate calculation of the composite deformed dose. However, it is important to account for tumor motion and body deformation especially when the tumor volume is small and/or located in the lower lobe of the lung.     

Analysis of late toxicity associated with external beam radiation therapy for prostate cancer with uniform setting of classical 4-field 70 Gy in 35 fractions: a survey study by the Osaka Urological Tumor Radiotherapy Study Group

We aimed to analyse late toxicity associated with external beam radiation therapy (EBRT) for prostate cancer using uniform dose-fractionation and beam arrangement, with the focus on the effect of 3D (CT) simulation and portal field size. We collected data concerning patients with localized prostate adenocarcinoma who had been treated with EBRT at five institutions in Osaka, Japan, between 1998 and 2006. All had been treated with 70 Gy in 35 fractions, using the classical 4-field technique with gantry angles of 0°, 90°, 180° and 270°. Late toxicity was evaluated strictly in terms of the Common Terminology Criteria for Adverse Events Version 4.0. In total, 362 patients were analysed, with a median follow-up of 4.5 years (range 1.0–11.6). The 5-year overall and cause-specific survival rates were 93% and 96%, respectively. The mean ± SD portal field size in the right–left, superior–inferior, and anterior–posterior directions was, respectively, 10.8 ± 1.1, 10.2 ± 1.0 and 8.8 ± 0.9 cm for 2D simulation, and 8.4 ± 1.2, 8.2 ± 1.0 and 7.7 ± 1.0 cm for 3D simulation (P < 0.001). No Grade 4 or 5 late toxicity was observed. The actuarial 5-year Grade 2–3 genitourinary and gastrointestinal (GI) late toxicity rates were 6% and 14%, respectively, while the corresponding late rectal bleeding rate was 23% for 2D simulation and 7% for 3D simulation (P < 0.001). With a uniform setting of classical 4-field 70 Gy/35 fractions, the use of CT simulation and the resultant reduction in portal field size were significantly associated with reduced late GI toxicity, especially with less rectal bleeding.     

Usefulness of double dose contrast-enhanced magnetic resonance imaging for clear delineation of gross tumor volume in stereotactic radiotherapy treatment planning of metastatic brain tumors: a dose comparison study

The purpose of this study was to compare the size and clearness of gross tumor volumes (GTVs) of metastatic brain tumors on T1-weighted magnetic resonance images between a single dose contrast administration protocol and a double dose contrast administration protocol to determine the optimum dose of contrast-enhancement for clear delineation of GTV in stereotactic radiotherapy (SRT). A total of 28 small metastatic brain tumors were evaluated in 13 patients by intra-individual comparison of GTV measurements using single dose and double dose contrast-enhanced thin-slice (1-mm) magnetic resonance imaging (MRI). All patients had confirmed histological types of primary tumors and had undergone hypo-fractionated SRT for metastatic brain tumors. The mean tumor diameter with single dose and double dose contrast-enhancement was 12.0 ± 1.1 mm and 13.2 ± 1.1 mm respectively (P < 0.001). The mean incremental ratio (MIR) obtained by comparing mean tumor diameters was 11.2 ± 0.02 %. The mean volume of GTV-1 (single dose contrast-enhancement) and GTV-2 (double dose contrast-enhancement) was 1.38 ± 0.41 ml and 1.59 ± 0.45 ml respectively (P < 0.01). The MIR by comparing mean tumor volumes was 32.3 ± 0.4 %. The MIR of GTV-1 with < 1ml volume and GTV-1 with > 1ml volume was 41.8 ± 0.05 % and 12.4 ± 0.03 % respectively (P < 0.001). We conclude that double dose contrast-enhanced thin-slice MRI is a more useful technique than single dose contrast-enhanced thin-slice MRI, especially for clear delineation of GTVs of small metastatic brain tumors in treatment planning of highly precise SRT.     

Assessment of transposed ovarian movement: how much of a safety margin should be added during pelvic radiotherapy?

The purpose of this study was to analyze transposed ovarian movement. Data from 27 patients who underwent ovarian transposition after surgical treatment for uterine cancer were retrospectively analyzed. Computed tomography (CT) images including transposed ovaries were superimposed on other CT images acquired at different times, and were matched on bony structures. Differences in ovarian position between the CT images were measured. The planning organ at risk volume (PRV) margins were calculated from the formula of the 90% reference intervals (RIs) and the 95% RI, which were defined as mean ± 1.65 standard deviation (SD) and mean ± 1.96 SD, respectively. The 90% RI in the cranial, caudal, anterior, posterior, left and right directions were 1.5, 1.5, 1.4, 1.0, 1.7 and 0.9 cm, respectively. The 95% RI in the corresponding directions were 1.5, 2.0, 1.7, 1.2, 1.9 and 1.2 cm, respectively. These data suggest that bilateral ovaries need a PRV margin of ∼2 cm in all directions. The present study suggests that a transposed ovary needs the same PRV margin as a normal ovary (∼2 cm). Even after transposition, ovaries should be kept away from the radiation field to take into consideration the degree of ovarian movement.     

Prediction of the minimum spacer thickness required for definitive radiotherapy with carbon ions and photons for pelvic tumors: an in silico planning study using virtual spacers

We aimed to predict the minimum distance between a tumor and the gastrointestinal (GI) tract that can satisfy the dose constraint by creating simulation plans with carbon-ion (C-ion) radiotherapy (RT) and photon RT for each case assuming insertion of virtual spacers of various thicknesses. We enrolled 55 patients with a pelvic tumor adjacent to the GI tract. Virtual spacers were defined as the overlap volume between the GI tract and the volume expanded 7–17 mm from the gross tumor volume (GTV). Simulation plans (70 Gy in 35 fractions for at least 95% of the planning target volume [PTV]) were created with the lowest possible dose to the GI tract under conditions that meet the dose constraints of the PTV. We defined the minimum thickness of virtual spacers meeting D2 cc of the GI tract <50 Gy as ‘MTS’. Multiple regression was used with explanatory variables to develop a model to predict MTS. We discovered that MTSs were at most 9 mm and 13 mm for C-ion RT and photon RT plans, respectively. The volume of overlap between the GI tract and a virtual spacer of 14 mm in thickness (OV14)-PTV was found to be the most important explanatory variable in the MTS prediction equation for both C-ion and photon RT plans. Multiple R² values for the regression model were 0.571 and 0.347 for C-ion RT and photon RT plans, respectively. In conclusion, regression equations were developed to predict MTS in C-ion RT and photon RT.     

Development of raster scanning IMRT using a robotic radiosurgery system

Treatment time with the CyberKnife frameless radiosurgery system is prolonged due to the motion of the robotic arm. We have developed a novel scanning irradiation method to reduce treatment time. We generated treatment plans mimicking eight-field intensity-modulated radiotherapy (IMRT) plans generated for the Novalis radiosurgery system. 2D dose planes were generated with multiple static beam spots collimated by a fixed circular cone. The weights of the uniformly distributed beam spots in each dose plane were optimized using the attraction–repulsion model. The beam spots were converted to the scanning speed to generate the raster scanning plan. To shorten treatment time, we also developed a hybrid scanning method which combines static beams with larger cone sizes and the raster scanning method. Differences between the Novalis and the scanning plan’s dose planes were evaluated with the criterion of a 5% dose difference. The mean passing rates of three cases were > 85% for cone sizes ≤ 12.5 mm. Although the total monitor units (MU) increased for smaller cone sizes in an inverse-square manner, the hybrid scanning method greatly reduced the total MU, while maintaining dose distributions comparable to those with the Novalis plan. The estimated treatment time of the hybrid scanning with a 12.5 mm cone size was on average 22% shorter than that of the sequential plans. This technique will be useful in allowing the CyberKnife with conventional circular cones to achieve excellent dose distribution with a shortened treatment time.