Fredrik Nordström

Technical Note: MRI only prostate radiotherapy planning using the statistical decomposition algorithm.

PURPOSE In order to enable a magnetic resonance imaging (MRI) only workflow in radiotherapy treatment planning, methods are required for generating Hounsfield unit (HU) maps (i.e., synthetic computed tomography, sCT) for dose calculations, directly from MRI. The Statistical Decomposition Algorithm (SDA) is a method for automatically generating sCT images from a single MR image volume, based on automatic tissue classification in combination with a model trained using a multimodal template material. This study compares dose calculations between sCT generated by the SDA and conventional CT in the male pelvic region. METHODS The study comprised ten prostate cancer patients, for whom a 3D T2 weighted MRI and a conventional planning CT were acquired. For each patient, sCT images were generated from the acquired MRI using the SDA. In order to decouple the effect of variations in patient geometry between imaging modalities from the effect of uncertainties in the SDA, the conventional CT was nonrigidly registered to the MRI to assure that their geometries were well aligned. For each patient, a volumetric modulated arc therapy plan was created for the registered CT (rCT) and recalculated for both the sCT and the conventional CT. The results were evaluated using several methods, including mean average error (MAE), a set of dose-volume histogram parameters, and a restrictive gamma criterion (2% local dose/1 mm). RESULTS The MAE within the body contour was 36.5 ± 4.1 (1 s.d.) HU between sCT and rCT. Average mean absorbed dose difference to target was 0.0% ± 0.2% (1 s.d.) between sCT and rCT, whereas it was -0.3% ± 0.3% (1 s.d.) between CT and rCT. The average gamma pass rate was 99.9% for sCT vs rCT, whereas it was 90.3% for CT vs rCT. CONCLUSIONS The SDA enables a highly accurate MRI only workflow in prostate radiotherapy planning. The dosimetric uncertainties originating from the SDA appear negligible and are notably lower than the uncertainties introduced by variations in patient geometry between imaging sessions.

Control chart analysis of data from a multicenter monitor unit verification study.

BACKGROUND AND PURPOSE This study aims to investigate the process of monitor unit verification using control charts. Control charts is a key tool within statistical process control (SPC), through which process characteristics can be visualized, usually chronologically with statistically determined limits. MATERIAL AND METHODS Our group has developed a monitor unit verification software that has been adopted at several Swedish institutions for pre-treatment verification of radiotherapy treatments. Deviations between point dose calculations using the treatment planning systems and using the independent monitor unit verification software from 9219 treatment plans and five different institutions were included in this multicenter study. The process of monitor unit verification was divided into subprocesses. Each subprocess was analyzed using probability plots and control charts. RESULTS Differences in control chart parameters for the investigated subprocesses were found between different treatment sites and different institutions, as well as between different treatment techniques. 19 of 37 subprocesses met the clinical specification (± 5%), i.e. process capability index was equal to or above one. CONCLUSIONS Control charts were found to be a useful tool for continuous analysis of data from the monitor unit verification software for patient specific quality control, as well as for comparisons between different institutions and treatment sites. The derived control chart limits were in agreement with AAPM TG114 guidelines on action levels.

Assessment of dosimetric impact of system specific geometric distortion in an MRI only based radiotherapy workflow for prostate

Dosimetric errors in a magnetic resonance imaging (MRI) only radiotherapy workflow may be caused by system specific geometric distortion from MRI. The aim of this study was to evaluate the impact on planned dose distribution and delineated structures for prostate patients, originating from this distortion. A method was developed, in which computer tomography (CT) images were distorted using the MRI distortion field. The displacement map for an optimized MRI treatment planning sequence was measured using a dedicated phantom in a 3 T MRI system. To simulate the distortion aspects of a synthetic CT (electron density derived from MR images), the displacement map was applied to CT images, referred to as distorted CT images. A volumetric modulated arc prostate treatment plan was applied to the original CT and the distorted CT, creating a reference and a distorted CT dose distribution. By applying the inverse of the displacement map to the distorted CT dose distribution, a dose distribution in the same geometry as the original CT images was created. For 10 prostate cancer patients, the dose difference between the reference dose distribution and inverse distorted CT dose distribution was analyzed in isodose level bins. The mean magnitude of the geometric distortion was 1.97 mm for the radial distance of 200-250 mm from isocenter. The mean percentage dose differences for all isodose level bins, were  ⩽0.02% and the radiotherapy structure mean volume deviations were  <0.2%. The method developed can quantify the dosimetric effects of MRI system specific distortion in a prostate MRI only radiotherapy workflow, separated from dosimetric effects originating from synthetic CT generation. No clinically relevant dose difference or structure deformation was found when 3D distortion correction and high acquisition bandwidth was used. The method could be used for any MRI sequence together with any anatomy of interest.

Technical evaluation of a laser-based optical surface scanning system for prospective and retrospective breathing adapted computed tomography

Abstract Background. For breathing adapted radiotherapy, the same motion monitoring system can be used for imaging and triggering of the accelerator. Purpose. To evaluate a new technique for prospective gated computed tomography (CT) and four-dimensional CT (4DCT) using a laser based surface scanning system (Sentinel™, C-RAD, Uppsala, Sweden). The system was compared to the AZ-733V respiratory gating system (Anzai Medical, Tokyo, Japan) and the Real-Time Position Management System (RPM™) (Varian Medical Systems, Palo Alto, CA, USA). Material and methods. Temporal accuracy was evaluated using a moving phantom programmed to move a platform along trajectories following a sin6(ωt) function with amplitudes from 6 to 20 mm and periods from 2 to 5 s during 120 s while the motion was recorded. The recorded data was Fourier transformed and the peak area at the fundamental and harmonic frequencies compared to data generated using the same sinusoidal function. For verification of the 4DCT reconstruction process, the phantom was programmed to move along a sinusoidal trajectory. Ten phase series were reconstructed. The distance from the couch to the platform was measured in each image. By fitting the function sin(ωt-φ) to the values measured in the images corresponding to each slice, the phase of each image was verified. Results and conclusion. In the recorded data, the peak area at the fundamental frequency covered on average 104 ± 4%, 102 ± 4% and 91 ± 27% of the peak area in the generated data for the Sentinel™, RPM™ and AZ-733V systems, respectively. All systems managed to resolve both harmonic frequencies. The second experiment showed that all images were sorted into the correct series using breathing data recorded by each system. The systems generated very similar results, however, it is preferable to use the same system both for imaging and treatment.

4D dosimetry and its applications to pre-treatment quality control and real-time in vivo dosimetry of VMAT treatments

In this study, a 4D dosimetry concept was developed. This concept included a method for calculation of 3D reference absorbed dose matrices at every control point of the delivery using a clinical treatment planning system (TPS). Further, the gamma evaluation method was extended to incorporate the 4th dimension of the TPS calculated dose distributions. The applications of the 4D dosimetry concept on pre-treatment quality control and real-time in vivo dosimetry were investigated.

Ensuring the integrity of treatment parameters throughout the radiotherapy process.

BACKGROUND AND PURPOSE Ensuring data integrity in radiotherapy is of major importance and a complex task. The aim of this study was to compare three different combinations of treatment planning and record and verify systems with respect to data integrity. MATERIALS AND METHODS A software for comparison of treatment parameters in DICOM-RT files was developed using the MATLAB R2010a (MathWorks Inc.) environment. One hundred treatment plans were analyzed for each system combination. In the first step of the analysis, all parameters were compared and a normal condition for each system combination was identified. The second step focused on the discovery of potential special cause deviations, e.g. by applying tolerance levels. RESULTS In total, 15% and 0.37% of all comparisons failed to meet the defined integrity demands in step 1 and step 2 of the analysis, respectively. Differences in the data integrity level between the systems were observed, ranging on average from 3.1 to 11.9 discrepancies per beam for the different RV-TPS combinations. CONCLUSIONS The proposed method can be used to increase the safety for individual patients by ensuring that the intended treatment is delivered. The system combination with the highest level of data integrity was found to be the one which shares a single database.

Tumor-tracking radiotherapy of moving targets; verification using 3D polymer gel, 2D ion-chamber array and biplanar diode array

The aim of this study was to carry out a dosimetric verification of a dynamic multileaf collimator (DMLC)-based tumor-tracking delivery during respiratory-like motion. The advantage of tumor-tracking radiation delivery is the ability to allow a tighter margin around the target by continuously following and adapting the dose delivery to its motion. However, there are geometric and dosimetric uncertainties associated with beam delivery system constraints and output variations, and several investigations have to be accomplished before a clinical integration of this tracking technique. Two types of delivery were investigated in this study I) a single beam perpendicular to a target with a one dimensional motion parallel to the MLC moving direction, and II) an intensity modulated arc delivery (RapidArc®) with a target motion diagonal to the MLC moving direction. The feasibility study (I) was made using an 2D ionisation chamber array and a true 3D polymer gel. The arc delivery (II) was verified using polymer gel and a biplanar diode array. Good agreement in absorbed dose was found between delivery to a static target and to a moving target with DMLC tracking using all three detector systems. However, due to the limited spatial resolution of the 2D array a detailed comparison was not possible. The RapidArc® plan delivery was successfully verified using the biplanar diode array and true 3D polymer gel, and both detector systems could verify that the DMLC-based tumor-tracking delivery system has a very good ability to account for respiratory target motion.

Motion induced interplay effects for VMAT radiotherapy

The purpose of this study was to develop a method to simulate breathing motion induced interplay effects for volumetric modulated arc therapy (VMAT), to verify the proposed method with measurements, and to use the method to investigate how interplay effects vary with different patient- and machine specific parameters. VMAT treatment plans were created on a virtual phantom in a treatment planning system (TPS). Interplay effects were simulated by dividing each plan into smaller sub-arcs using an in-house developed software and shifting the isocenter for each sub-arc to simulate a sin6 breathing motion in the superior-inferior direction. The simulations were performed for both flattening-filter (FF) and flattening-filter free (FFF) plans and for different breathing amplitudes, period times, initial breathing phases, dose levels, plan complexities, CTV sizes, and collimator angles. The resulting sub-arcs were calculated in the TPS, generating a dose distribution including the effects of motion. The interplay effects were separated from dose blurring and the relative dose differences to 2% and 98% of the CTV volume (ΔD98% and ΔD2%) were calculated. To verify the simulation method, measurements were carried out, both static and during motion, using a quasi-3D phantom and a motion platform. The results of the verification measurements during motion were comparable to the results of the static measurements. Considerable interplay effects were observed for individual fractions, with the minimum ΔD98% and maximum ΔD2% being  -16.7% and 16.2%, respectively. The extent of interplay effects was larger for FFF compared to FF and generally increased for higher breathing amplitudes, larger period times, lower dose levels, and more complex treatment plans. Also, the interplay effects varied considerably with the initial breathing phase, and larger variations were observed for smaller CTV sizes. In conclusion, a method to simulate motion induced interplay effects was developed and verified with measurements, which allowed for a large number of treatment scenarios to be investigated. The simulations showed large interplay effects for individual fractions and that the extent of interplay effects varied with the breathing pattern, FFF/FF, dose level, CTV size, collimator angle, and the complexity of the treatment plan.

OC-0156: MRI only prostate radiotherapy using synthetic CT images

Purpose or Objective: Introducing an MRI-only workflow into the radiotherapy clinic, requires that MR-images can be used both for treatment planning calculations and for patient positioning. The two-fold aim of this study was to evaluate the use of MR-images with respect to 1) the accuracy of treatment planning dose calculations, and 2) the reliability of fiducial marker identification for patient positioning.

Verification of motion induced thread effect during tomotherapy using gel dosimetry

The purpose of the study was to evaluate how breathing motion during tomotherapy (Accuray, CA, USA) treatment affects the absorbed dose distribution. The experiments were carried out using gel dosimetry and a motion device simulating respiratory-like motion (HexaMotion, ScandiDos, Uppsala, Sweden). Normoxic polyacrylamide gels (nPAG) were irradiated, both during respiratory-like motion and in a static mode. To be able to investigate interplay effects the static absorbed dose distribution was convolved with the motion function and differences between the dynamic and convolved static absorbed dose distributions were interpreted as interplay effects. The expected dose blurring was present and the interplay effects formed a spiral pattern in the lower dose volume. This was expected since the motion induced affects the preset pitch and the theoretically predicted thread effect may emerge. In this study, the motion induced thread effect was experimentally verified for the first time.