Optimizing MRI sequences and images for MRI-based stereotactic radiosurgery treatment planning

AimDevelopment of MRI sequences and processing methods for the production of images appropriate for direct use in stereotactic radiosurgery (SRS) treatment planning.BackgroundMRI is useful in SRS treatment planning, especially for patients with brain lesions or anatomical targets that are poorly distinguished by CT, but its use requires further refinement. This methodology seeks to optimize MRI sequences to generate distortion-free and clinically relevant MR images for MRI-only SRS treatment planning.Materials and methodsWe used commercially available SRS MRI-guided radiotherapy phantoms and eight patients to optimize sequences for patient imaging. Workflow involved the choice of correct MRI sequence(s), optimization of the sequence parameters, evaluation of image quality (artifact free and clinically relevant), measurement of geometrical distortion, and evaluation of the accuracy of our offline correction algorithm.ResultsCT images showed a maximum deviation of 1.3 mm and minimum deviation of 0.4 mm from true fiducial position for SRS coordinate definition. Interestingly, uncorrected MR images showed maximum deviation of 1.2 mm and minimum of 0.4 mm, comparable to CT images used for SRS coordinate definition. After geometrical correction, we observed a maximum deviation of 1.1 mm and minimum deviation of only 0.3 mm.ConclusionOur optimized MRI pulse sequences and image correction technique show promising results; MR images produced under these conditions are appropriate for direct use in SRS treatment planning.     

Technical note: Unexpected external markers artifact in 3D k-space based parallel imaging turbo spin-echo magnetic resonance imaging

PurposeMRI for radiotherapy planning requires spatial referencing using immobilization devices and markers. Clinical images of a difficult-to-interpret artifact are presented, resembling a metastasis, which occurs when combining CAIPIRINHA k-space-based parallel imaging (PI), 3D distortion correction, and external markers.MethodsA 3D variable flip angle Turbo Spin Echo sequence was used on a 1.5 T and 3 T MRI using flexible and head and neck coils. Two types of markers were tested: Liquimark LM1 and Spee-D-Mark. A silicone oil phantom was used that represents low signal intensity, such as gray matter. 3D Fourier transforms were also used to show the issue’s origin.ResultsThe markers can appear in an unexpected region of a patient, not in the same original or reconstructed slice nor in a rectilinear direction in a slice, especially when using CAIPIRINHA acceleration with 3D distortion correction. The probability of occurrence was respectively 13% and 80% for distances of <=2 mm and >2 mm between marker and patient, for example when using thermoplastic masks. Clinical cases are shown where this semi-randomly occurring artifact appears post contrast only, and thus can be interpreted as metastases. The artifact did not appear when using compressed sensing acceleration.ConclusionMarkers used for radiotherapy MRI application can introduce additional artifacts that can be interpreted as metastases. However, other high signal intensity structures on the surface of a patient, such as the ear, can lead to an equivalent error.     

Geometric inaccuracy and co-registration errors for CT,DynaCT and MRI images used in robotic stereotactic radiosurgery treatment planning

PurposeTo measure the combined errors due to geometric inaccuracy and image co-registration on secondary images (dynamic CT angiography (dCTA), 3D DynaCT angiography (DynaCTA), and magnetic resonance images (MRI)) that are routinely used to aid in target delineation and planning for stereotactic radiosurgery (SRS).MethodsThree phantoms (one commercial and two in-house built) and two different analysis approaches (commercial and MATLAB based) were used to quantify the magnitude of geometric image distortion and co-registration errors for different imaging modalities within CyberKnife’s MultiPlan treatment planning software. For each phantom, the combined errors were reported as a mean target registration error (TRE). The mean TRE’s for different intramodality imaging parameters (e.g., mAs, kVp, and phantom set-ups) and for dCTA, DynaCTA, and MRI systems were measured.ResultsOnly X-ray based imaging can be performed with the commercial phantom, and the mean TRE ± standard deviation values were large compared to the in-house analysis using MATLAB. With the 3D printed phantom, even drastic changes in treatment planning CT imaging protocols did not greatly influence the mean TRE (<0.5 mm for a 1 mm slice thickness CT). For all imaging modalities, the largest mean TRE was found on DynaCT, followed by T2-weighted MR images (albeit all <1 mm).ConclusionsThe user may overestimate the mean TRE if the commercial phantom and MultiPlan were used solely. The 3D printed phantom design is a sensitive and suitable quality assurance tool for measuring 3D geometric inaccuracy and co-registration errors across all imaging modalities.     

Dose delivery accuracy on helical tomotherapy for 4-dimensional tumor motion — a phantom study

BackgroundThe advances in image guidance and capability of highly conformal dose deliveries made possible the use of helical tomotherapy (HT) for lung cancer treatment. To determine the effect of respiratory motion on the delivered dose in HT, film dosimetry using a dynamic phantom was performed. This was a phantom study to determine the effect of motion on the delivered dose in HT.Materials and methods4D computed tomography (4DCT) was acquired for various target motions of CIRS dynamic phantom (CIRS Inc., Norfolk, USA) with 2.5cm diameter spherical target of volume 8.2 cc moving in the COS⁴ motion pattern. AveIP images and treatment plans were generated in the HT planning system. Target excursions during treatment delivery were changed in the superior-inferior, anteroposterior and lateral directions. The breathing cycle time was varied from 4 to 5 sec. and also the delivery interruptions were introduced. A film was exposed for each delivery and gamma analysis was performed.ResultsThe gamma pass rate (GPR) with 3%, 2 mm criteria for the target motion in the S-I direction showed a significant reduction from 97.5% to 54.4% as the motion increased from 3 mm to 8 mm (p = 0.03). For the target motion in S-I = 8 mm, L-R = A-P = 3 mm, the percentage decrease in the GPR was 74% (p = 0.001) for three interruptions.ConclusionThe ITV based approach in HT is ideal for a shallow breathing situation when the tumor excursions were confined to 5 mm in the S-I and 3 mm in L-R and A-P directions.     

A practical methodology to improve the dosimetric accuracy of MR-based radiotherapy simulation for brain tumors

PurposeTo investigate the dosimetric accuracy of synthetic computed tomography (sCT) images generated by a clinically-ready voxel-based MRI simulation package, and to develop a simple and feasible method to improve the accuracy.Methods20 patients with brain tumor were selected to undergo CT and MRI simulation. sCT images were generated by a clinical MRI simulation package. The discrepancy between planning CT and sCT in CT number and body contour were evaluated. To resolve the discrepancies, an sCT specific CT-relative electron density (RED) calibration curve was used, and a layer of pseudo-skin was created on the sCT. The dosimetric impact of these discrepancies, and the improvement brought about by the modifications, were evaluated by a planning study. Volumetric modulated arc therapy (VMAT) treatment plans for each patient were created and optimized on the planning CT, which were then transferred to the original sCT and the modified-sCT for dose re-calculation. Dosimetric comparisons and gamma analysis between the calculated doses in different images were performed.ResultsThe average gamma passing rate with 1%/1 mm criteria was only 70.8% for the comparison of dose distribution between planning CT and original sCT. The mean dose difference between the planning CT and the original sCT were −1.2% for PTV D₉₅ and −1.7% for PTV D_max, while the mean dose difference was within 0.7 Gy for all relevant OARs. After applying the modifications on the sCT, the average gamma passing rate was increased to 92.2%. Mean dose difference in PTV D₉₅ and D_max were reduced to −0.1% and −0.3% respectively. The mean dose difference was within 0.2 Gy for all OAR structures and no statistically significant difference were found.ConclusionsThe modified-sCT demonstrated improved dosimetric agreement with the planning CT. These results indicated the overall dosimetric accuracy and practicality of this improved MR-based treatment planning method.     

Definition of the margin of major coronary artery bifurcations during radiotherapy with electrocardiograph-gated 4D-CT

PurposeThe aim was to measure the cardiac motion-induced displacements of major coronary artery bifurcations utilizing electrocardiography (ECG)-gated four-dimensional computed tomography (4D-CT) and to determine the margin of coronary artery bifurcations.MethodsThirty-seven female patients who underwent retrospective ECG-gated 4D-CT in inspiratory breath hold (IBH) were enrolled. The left main coronary artery bifurcation (LM), the obtuse marginal branch bifurcation (OM), the first diagonal branch bifurcation (D₁), the second diagonal branch bifurcation (D₂), the caudal portion of the left anterior descending branch (APX), the first right ventricular artery bifurcation (V) and the acute marginal branch bifurcation (AM) were contoured. The center of the contour of the coronary arterial bifurcations at end systole was defined as the standard, and the margin were then calculated.ResultsThe margin in the left–right (LR), cranio-caudal (CC), and anterior-posterior (AP) coordinates were as follows: LM 3, 3, and 3 mm; D₁ 6, 3, and 3 mm; D₂ 3, 3, and 3 mm; APX 4, 4, and 4 mm; OM 4, 6, and 5 mm; V 6, 8, and 7 mm; and AM 6, 8, and 7 mm, respectively.ConclusionCoronary artery bifurcations should be considered a separate organ at risk (OAR), and different margin should be provided due to the differences resulting from motion displacement. The maximum margin in the LR, CC, and AP coordinates of left coronary artery bifurcations were 6, 6, and 5 mm, and those of the right coronary artery bifurcations were 6, 8, and 7 mm, respectively.     

Monte Carlo simulation of the dose distribution of ICRP adult reference computational phantoms for acquisitions with a 320 detector-row cone-beam CT scanner

PurposeThe purpose of this study was to develop and validate a Monte Carlo (MC) simulation tool for patient dose assessment for a 320 detector-row CT scanner, based on the recommendations of International Commission on Radiological Protection (ICRP). Additionally, the simulation was applied on four clinical acquisition protocols, with and without automatic tube current modulation (TCM).MethodsThe MC simulation was based on EGS4 code and was developed specifically for a 320 detector-row cone-beam CT scanner. The ICRP adult reference phantoms were used as patient models. Dose measurements were performed free-in-air and also in four CTDI phantoms: 150 mm and 350 mm long CT head and CT body phantoms. The MC program was validated by comparing simulations results with these actual measurements acquired under the same conditions. The measurements agreed with the simulations across all conditions within 5%. Patient dose assessment was performed for four clinical axial acquisitions using the ICRP adult reference phantoms, one of them using TCM.ResultsThe results were nearly always lower than those obtained from other dose calculator tools or published in other studies, which were obtained using mathematical phantoms in different CT systems. For the protocol with TCM organ doses were reduced by between 28 and 36%, compared to the results obtained using a fixed mA value.ConclusionsThe developed simulation program provides a useful tool for assessing doses in a 320 detector-row cone-beam CT scanner using ICRP adult reference computational phantoms and is ready to be applied to more complex protocols.     

Characterization and Correction of Geometric Distortions in 814 Diffusion Weighted Images

Introduction

Diffusion Weighted Imaging (DWI), which is based on Echo Planar Imaging (EPI) protocols, is becoming increasingly important for neurosurgical applications. However, its use in this context is limited in part by significant spatial distortion inherent to EPI.

Method

We evaluated an efficient algorithm for EPI distortion correction (EPIC) across 814 DWI scans from 250 brain tumor patients and quantified the magnitude of geometric distortion for whole brain and multiple brain regions.

Results

Evaluation of the algorithm’s performance revealed significantly higher mutual information between T1-weighted pre-contrast images and corrected b = 0 images than the uncorrected b = 0 images (p < 0.001). The distortion magnitude across all voxels revealed a median EPI distortion effect of 2.1 mm, ranging from 1.2 mm to 5.9 mm, the 5^th and 95^th percentile, respectively. Regions adjacent to bone-air interfaces, such as the orbitofrontal cortex, temporal poles, and brain stem, were the regions most severely affected by DWI distortion.

Conclusion

Using EPIC to estimate the degree of distortion in 814 DWI brain tumor images enabled the creation of a topographic atlas of DWI distortion across the brain. The degree of displacement of tumors boundaries in uncorrected images is severe but can be corrected for using EPIC. Our results support the use of distortion correction to ensure accurate and careful application of DWI to neurosurgical practice.     

Tridimensional dose evaluation of the respiratory motion influence on breast radiotherapy treatments using conformal radiotherapy,forward IMRT,and inverse IMRT planning techniques

PurposeTo evaluate the respiratory motion influence on the tridimensional (3D) dose delivery to breast-shaped phantoms using conformal radiotherapy (3D-RT), Field-in Field (FiF), and IMRT planning techniques.MethodsThis study used breast-shaped phantoms filled with MAGIC-f gel dosimeter to simulate the breast, and an oscillation platform to simulate the respiratory motion. The platform allowed motion in the anterior-posterior direction with oscillation amplitudes of 0.34 cm, 0.88 cm, and 1.22 cm. CT images of the static phantom were used for the 3D-RT, FiF, and IMRT treatment planning. Five phantoms were prepared and irradiated for each planning technique evaluated. Phantom 1 was irradiated static, phantoms 2–4 were irradiated moving with the three different motion amplitudes, and phantom 5 was used as a reference. The 3D dose distributions were obtained by relaxometry of magnetic resonance imaging, and the respiratory motion influence in the doses distribution was accessed by gamma evaluations (3%/3mm/15% threshold) comparing the measurements of the phantoms irradiated under movement with the static ones.ResultsThe mean gamma approvals for three oscillatory amplitudes were 96.44%, 93.23%, and 91.65%; 98.42%, 95.66%, and 94.31%; and 94.49%, 93.51%, and 86.62% respectively for 3D-RT, FiF and IMRT treatments. A gamma results profile per slice along the phantom showed that for FiF and IMRT irradiations, most of the failures occurred in the central region of the phantom.ConclusionsBy increasing the respiratory motion movement, the dose distribution variations for the three planning techniques were more pronounced, being the FiF technique variations the smallest one.     

Translational and rotational localization errors in cone-beam CT based image-guided lung stereotactic radiotherapy

PurposeAccurate localization is crucial in delivering safe and effective stereotactic body radiation therapy (SBRT). The aim of this study was to analyse the accuracy of image-guidance using the cone-beam computed tomography (CBCT) of the VERO system in 57 patients treated for lung SBRT and to calculate the treatment margins.Materials and methodsThe internal target volume (ITV) was obtained by contouring the tumor on maximum and mean intensity projection CT images reconstructed from a respiration correlated 4D-CT. Translational and rotational tumor localization errors were identified by comparing the manual registration of the ITV to the motion-blurred tumor on the CBCT and they were corrected by means of the robotic couch and the ring rotation. A verification CBCT was acquired after correction in order to evaluate residual errors.ResultsThe mean 3D vector at initial set-up was 6.6 ± 2.3 mm, which was significantly reduced to 1.6 ± 0.8 mm after 6D automatic correction. 94% of the rotational errors were within 3°. The PTV margins used to compensate for residual tumor localization errors were 3.1, 3.5 and 3.3 mm in the LR, SI and AP directions, respectively.ConclusionsOn-line image guidance with the ITV–CBCT matching technique and automatic 6D correction of the VERO system allowed a very accurate tumor localization in lung SBRT.     

The effect of scatter correction and radius of rotation on semiquantitative measurements in SPECT 123I-FP-CIT imaging. A phantom study

PurposeThe high energy emissions of ¹²³I and the suboptimal radius of rotation affect the semiquantitative measurements performed during ¹²³I-FP-CIT tomographic imaging. An in-house extra low cost striatum phantom with brain and striatum compartments was constructed and was used to study the effects of Triple Energy Window scatter correction (TEW-SC) and radius of rotation on the Specific Binding Ratio (SBR) measurements.Materials and methodsThe phantom compartments were filled with radioactive ¹²³I solutions with varying concentrations, in a series of experiments. Tomographic images were acquired at six different radii of rotation, with and without TEW-SC and the SBRs were calculated using appropriate regions of interest, as in clinical imaging.ResultsSBRs decreased with increasing radius of rotation in both non-SC and TEW-SC images, the decrease being more pronounced in the latter. The application of TEW-SC increases SBR values by 40% on average. A maximum %Recovery of 42.7% of the true SBR value was achieved in the non-SC images, which increased to 64.6% after TEW-SC. Appropriate correction factors (CF) were calculated in order to make the SBR values independent on the radius of rotation, which could be used to correct SBR values obtained from tomographic acquisitions with suboptimal radius of rotation.ConclusionThe use of appropriate CF can provide more consistent SBR values and a more meaningful comparison between SBRs calculated from images acquired at different radii of rotation.     

In phantom evaluation of targeting accuracy in MRI-based brain radiosurgery

PurposeTo determine the targeting accuracy of brain radiosurgery when planning procedures employing different MRI and MRI + CT combinations are adopted.Materials and methodA new phantom, the BrainTool, has been designed and realized to test image co-registration and targeting accuracy in a realistic anatomical situation. The phantom was created with a 3D printer and materials that mimic realistic brain MRI and CT contrast using a model extracted from a synthetic MRI study of a human brain. Eight markers distributed within the BrainTool provide for assessment of the accuracy of image registrations while two cavities that host an ionization chamber are used to perform targeting accuracy measurements with an iterative cross-scan method. Two procedures employing 1.5 T MRI-only or a combination of MRI (taken with 1.5 T or 3 T scanners) and CT to carry out Gamma Knife treatments were investigated. As distortions can impact targeting accuracy, MR images were preliminary evaluated to assess image deformation extent using GammaTool phantom.ResultsMR images taken with both scanners showed average and maximum distortion of 0.3 mm and 1 mm respectively. The marker distances in co-registered images resulted below 0.5 mm for both MRI scans. The targeting mismatches obtained were 0.8 mm, 1.0 mm and 1.2 mm for MRI-only and MRI + CT (1,5T and 3 T), respectively.ConclusionsProcedures using a combination of MR and CT images provide targeting accuracies comparable to those of MRI-only procedures. The BrainTool proved to be a suitable tool for carrying out co-registration and targeting accuracy of Gamma Knife brain radiosurgery treatments.     

High resolution MR eye protocol optimization: Comparison between 3D-CISS, 3D-PSIF and 3D-VIBE sequences

PurposeThe purpose of this study was to compare selected MRI pulse sequences and to evaluate their utility for depicting specific anatomic regions in the eye.MethodsA High-Resolution (HR) 0.08 × 0.08 × 0.60 mm³ MRI protocol was developed on a 1.5-T clinical system and applied in the left eye of an albino rabbit, utilizing a small field of view surface coil. The comprehensive MRI protocol consisted of two 3D (T2/T1)w sequences (3D-PSIF and 3D-CISS), and one 3D T1w sequence (3D-VIBE). The T1w 3D-VIBE sequence was acquired, before and after intravenous injection of 0.2 mmol/kgr gadolinium-DTPA. Signal-to-Noise Ratios (SNR) and Contrast-to-Noise Ratios (CNR) amongst specific eye anatomical areas were calculated for each sequence. The presence of artifacts was rated subjectively utilizing a 5 point scale.Results3D-PSIF and 3D-CISS provide better delineation and visualization of the eye as compared with 3D-VIBE sequences. 3D-CISS images present the highest SNR and revealed better discrimination of the ocular surrounding tissues; its basic drawback though is related to ghost artifacts appearing in the anterior chamber and resulting in the lowest image quality. In post-contrast imaging, the T1w 3D-VIBE sequence provided the best overall image quality. Moreover, 3D (T2/T1)w sequences can provide good anatomic depiction of the eye segments. Agreement between the two independent readers was good.ConclusionsOptimization of a comprehensive MR eye imaging protocol is achieved. A higher SNR, a better spatial resolution and a reduction of the total scan time were obtained, thus making clinical MRI systems more reliable in eye imaging.     

Impact of lung density on isolated lung tumor dose in VMAT using inline MR-Linac

PurposeThis study investigated the impact of lung density on the isolated lung tumor dose for volumetric modulated arc therapy (VMAT) in an inline magnetic resonance linear accelerator (MR-Linac) using the Monte Carlo (MC) simulation.MethodsCT images of the thorax phantoms with lung tumors of 1, 2, and 3 cm diameters were converted into voxel-base phantoms with lung densities of 0.1, 0.2, and 0.3 g/cm³, respectively. The dose distributions were calculated for partial-arc VMAT. The dose distributions were compared using dose differences, dose volume histograms, and dose volume indices.ResultsIn all cases, the inline magnetic field significantly enhanced the lung tumor dose compared to that at 0 T. For the 1 cm lung tumor, the inline magnetic field of 1 T increased the minimum dose of 95% of the Planning target volume (PTV D₉₅) by 14.0% in 0.1 g/cm³ lung density as compared to that in 0.3 g/cm³ at 0 T. In contrast, at 0 and 0.5 T, the PTV D₉₅ in 0.3 g/cm³ lung density was larger than that in lung density of 0.1 g/cm³. For the 2 cm lung tumor, a similar tendency to 1 cm was observed, whereas the dose impact of lung density was smaller than that for 1 cm. For the 3 cm lung tumor, the lung tumor dose was independent of lung density at 0.5 T and 1.0 T.ConclusionThe inline MR-Linac with the magnetic field over 1 T can enhance the PTV D₉₅ for VMAT regardless of the lung density.     

High Dynamic Range Processing for Magnetic Resonance Imaging

Purpose

To minimize feature loss in T₁- and T₂-weighted MRI by merging multiple MR images acquired at different T_R and T_E to generate an image with increased dynamic range.

Materials and Methods

High Dynamic Range (HDR) processing techniques from the field of photography were applied to a series of acquired MR images. Specifically, a method to parameterize the algorithm for MRI data was developed and tested. T₁- and T₂-weighted images of a number of contrast agent phantoms and a live mouse were acquired with varying T_R and T_E parameters. The images were computationally merged to produce HDR-MR images. All acquisitions were performed on a 7.05 T Bruker PharmaScan with a multi-echo spin echo pulse sequence.

Results

HDR-MRI delineated bright and dark features that were either saturated or indistinguishable from background in standard T₁- and T₂-weighted MRI. The increased dynamic range preserved intensity gradation over a larger range of T₁ and T₂ in phantoms and revealed more anatomical features in vivo.

Conclusions

We have developed and tested a method to apply HDR processing to MR images. The increased dynamic range of HDR-MR images as compared to standard T₁- and T₂-weighted images minimizes feature loss caused by magnetization recovery or low SNR.     

Implementation and evaluation of a transit dosimetry system for treatment verification

PurposeTo evaluate a formalism for transit dosimetry using a phantom study and prospectively evaluate the protocol on a patient population undergoing 3D conformal radiotherapy.MethodsAmorphous silicon EPIDs were calibrated for dose and used to acquire images of delivered fields. The measured EPID dose map was back-projected using the planning CT images to calculate dose at pre-specified points within the patient using commercially available software, EPIgray (DOSIsoft, France). This software compared computed back-projected dose with treatment planning system dose. A series of tests were performed on solid water phantoms (linearity, field size effects, off-axis effects). 37 patients were enrolled in the prospective study.ResultsThe EPID dose response was stable and linear with dose. For all tested field sizes the agreement was good between EPID-derived and treatment planning system dose in the central axis, with performance stability up to a measured depth of 18 cm (agreement within −0.5% at 10 cm depth on the central axis and within −1.4% at 2 cm off-axis). 126 transit images were analysed of 37 3D-conformal patients. Patient results demonstrated the potential of EPIgray with 91% of all delivered fields achieved the initial set tolerance level of Δ_D of 0 ± 5-cGy or %Δ_D of 0 ± 5%.ConclusionsThe in vivo dose verification method was simple to implement, with very few commissioning measurements needed. The system required no extra dose to the patient, and importantly was able to detect patient position errors that impacted on dose delivery in two of cases.     

Comparison of three different phantoms used for Winston-Lutz test with Artiscan software

BackgroundOne of the most important test in every quality assurances process of medical linear accelerators is the Winston-Lutz test, allowing an evaluation of the treatment isocentre in the light of uncertainty of the position of the collimator, the gantry and the couch.AimThe purpose of this work was analysis of the results of the Winston-Lutz test performed with three different phantoms for two different accelerators.Materials and methodsMeasurements were performed on two Varian machines: TrueBeam equipped with aS1200 EPID and TrueBeam equipped with aS1000 EPID. During the study three different phantoms dedicated for verification of the radiation isocentre were used: PTW Isoball, AQUILAB Isocentre Phantom and Varian Isocentre Cube. Analysis of the DICOM images was performed in Artiscan software.ResultsFor TrueBeam with as1200 EPID, gantry MV isocentre was about 0.18 mm larger for Varian Isocentre Cube than for two other phantoms used in this study. The largest variability of this parameter was observed for the couch. The results differed to 1.16 mm. For TrueBeam with as1000 EPID, results for collimator isocentre with PTW Isoball phantom were about 0.10 mm larger than for two other phantoms. For the gantry, results obtained with Varian Isocentre Cube were 0.21 mm larger.ConclusionThe obtained results for all three phantoms are within the accepted tolerance range. The largest differences were observed for treatment couch, which may be related to the phantom mobility during couch movement.     

Accuracy of dose calculation on iterative CBCT for head and neck radiotherapy

PurposeTo evaluate the feasibility of the use of iterative cone-beam computed tomography (CBCT) for dose calculation in the head and neck region.MethodsThis study includes phantom and clinical studies. All acquired CBCT images were reconstructed with Feldkamp–Davis–Kress algorithm-based CBCT (FDK-CBCT) and iterative CBCT (iCBCT) algorithm. The Hounsfield unit (HU) consistency between the head and body phantoms was determined in both reconstruction techniques. Volumetric modulated arc therapy (VMAT) plans were generated for 16 head and neck patients on a planning CT scan, and the doses were recalculated on FDK-CBCT and iCBCT with Anisotropic Analytical Algorithm (AAA) and Acuros XB (AXB). As a comparison of the accuracy of dose calculations, the absolute dosimetric difference and 1%/1 mm gamma passing rate analysis were analyzed.ResultsThe difference in the mean HU values between the head and body phantoms was larger for FDK-CBCT (max value: 449.1 HU) than iCBCT (260.0 HU). The median dosimetric difference from the planning CT were <1.0% for both FDK-CBCT and iCBCT but smaller differences were found with iCBCT (planning target volume D_50%: 0.38% (0.15–0.59%) for FDK-CBCT, 0.28% (0.13–0.49%) for iCBCT, AAA; 0.14% (0.04–0.19%) for FDK-CBCT, 0.07% (0.02–0.20%) for iCBCT). The mean gamma passing rate was significantly better in iCBCT than FDK-CBCT (AAA: 98.7% for FDK-CBCT, 99.4% for iCBCT; AXB: 96.8% for FDK_CBCT, 97.5% for iCBCT).ConclusionThe iCBCT-based dose calculation in VMAT for head and neck cancer was accurate compared to FDK-CBCT.     

Clinical radiotherapy application of N-vinylpyrrolidone-containing 3D polymer gel dosimeters with remote external MR-reading

PurposeAdvanced 3D dosimetry is required for verifications of complex dose distributions in modern radiotherapy. Two 3D polymer gel dosimeters, coupled with magnetic resonance (MR) imaging (3 T MRI) readout and data processing with polyGeVero® software, were tested for the verification of calculated 3D dose distributions by a treatment planning system (TPS) and ArcCHECK®–3DVH®, related to eradication of a lung tumour.MethodsN-vinylpyrrolidone-containing 3D polymer gel dosimeters were used: VIC (containing ascorbic acid and copper sulfate pentahydrate) and VIC-T (containing tetrakis(hydroxymethyl)phosphonium chloride). Three remote centers were involved in the dosimeters preparation and irradiation (Poland), and MRI (Austria). Cross beam calibration of the dosimeters and verification of a 3D dose distribution calculated with an Eclipse External Beam TPS and ArcCHECK®–3DVH® were performed. The 3D-to-3D comparisons of the VIC and VIC-T with TPS and ArcCHECK®–3DVH® along with ArcCHECK®–3DVH® versus TPS dose matrixes were performed with the aid of the polyGeVero® by analyzing dose profiles, isodoses lines, gamma index, gamma angle, dose difference, and related histograms.ResultsThe measured MR-relaxation rate (R₂ = 1/T₂) for the dosimeters relates to the dose, as follows: R₂ = 0.0928 ± 0.0008 [Gy⁻¹ s⁻¹] × D [Gy] + 2.985 ± 0.012 [s⁻¹] (VIC) and 0.1839 ± 0.0044 [Gy⁻¹ s⁻¹] × D [Gy] + 2.519 ± 0.053 [s⁻¹] (VIC-T). The 3D-to-3D comparisons revealed a good agreement between the measured and calculated 3D dose distributions.ConclusionsVIC and VIC-T with 3T MRI readout and polyGeVero® showed potential for verifications of calculated irradiation plans. The results obtained suggest the implementation of the irradiation plan for eradication of the lung tumour.