PubMed Health⌕ Search

Biomedical subjects

Ben J M Heijmen

Publications and source records attributed to Ben J M Heijmen.

16 recordsLinked to original sources

Computer optimization of noncoplanar beam setups improves stereotactic treatment of liver tumors.

PURPOSE: To investigate whether computer-optimized fully noncoplanar beam setups may improve treatment plans for the stereotactic treatment of liver tumors. METHODS: An algorithm for automated beam orientation and weight selection (Cycle) was extended for noncoplanar stereotactic treatments. For 8 liver patients previously treated in our clinic using a prescription isodose of 65%, Cycle was used to generate noncoplanar and coplanar plans with the highest achievable minimum planning target volume (PTV) dose for the clinically delivered isocenter and mean liver doses, while not violating the clinically applied hard planning constraints. The clinical and the optimized coplanar and noncoplanar plans were compared, with respect to D(PTV,99%), the dose received by 99% of the PTV, the PTV generalized equivalent uniform dose (gEUD), and the compliance with the clinical constraints. RESULTS: For each patient, the ratio between D(PTV,99%) and D(isoc), and the gEUD(-5) and gEUD(-20) values of the optimized noncoplanar plan were higher than for the clinical plan with an average increase of respectively 18.8% (range, 7.8-24.0%), 6.4 Gy (range, 3.4-11.8 Gy), and 10.3 Gy (range, 6.7-12.5). D(PTV,99%)/D(isoc), gEUD(-5), and gEUD(-20) of the optimized noncoplanar plan was always higher than for the optimized coplanar plan with an average increase of, respectively, 4.5% (range, 0.2-9.7%), 2.7 Gy (range, 0.6-9.7 Gy), and 3.4 Gy (range, 0.6-9.9 Gy). All plans were within the imposed hard constraints. On average, the organs at risk were better spared with the optimized noncoplanar plan than with the optimized coplanar plan and the clinical plan. CONCLUSIONS: The use of automatically generated, fully noncoplanar beam setups results in plans that are favorable compared with coplanar techniques. Because of the automation, we found that the planning workload can be decreased from 1 to 2 days to 1 to 2 h.

Algorithms↗

Dosimetric pre-treatment verification of IMRT using an EPID; clinical experience.

BACKGROUND AND PURPOSE: In our clinic a QA program for IMRT verification, fully based on dosimetric measurements with electronic portal imaging devices (EPID), has been running for over 3 years. The program includes a pre-treatment dosimetric check of all IMRT fields. During a complete treatment simulation at the linac, a portal dose image (PDI) is acquired with the EPID for each patient field and compared with a predicted PDI. In this paper, the results of this pre-treatment procedure are analysed, and intercepted errors are reported. An automated image analysis procedure is proposed to limit the number of fields that need human intervention in PDI comparison. MATERIALS AND METHODS: Most of our analyses are performed using the gamma index with 3% local dose difference and 3mm distance to agreement as reference values. Scalar parameters are derived from the gamma values to summarize the agreement between measured and predicted 2D PDIs. Areas with all pixels having gamma values larger than one are evaluated, making decisions based on clinically relevant criteria more straightforward. RESULTS: In 270 patients, the pre-treatment checks revealed four clinically relevant errors. Calculation of statistics for a group of 75 patients showed that the patient-averaged mean gamma value inside the field was 0.43 +/- 0.13 (1SD) and only 6.1 +/- 6.8% of pixels had a gamma value larger than one. With the proposed automated image analysis scheme, visual inspection of images can be avoided in 2/3 of the cases. CONCLUSION: EPIDs may be used for high accuracy and high resolution routine verification of IMRT fields to intercept clinically relevant dosimetric errors prior to the start of treatment. For the majority of fields, PDI comparison can fully rely on an automated procedure, avoiding excessive workload.

Electronic Data Processing↗

Limitations of the planning organ at risk volume (PRV) concept.

PURPOSE: Previously, we determined a planning target volume (PTV) margin recipe for geometrical errors in radiotherapy equal to M(T) = 2 Sigma + 0.7 sigma, with Sigma and sigma standard deviations describing systematic and random errors, respectively. In this paper, we investigated margins for organs at risk (OAR), yielding the so-called planning organ at risk volume (PRV). METHODS AND MATERIALS: For critical organs with a maximum dose (D(max)) constraint, we calculated margins such that D(max) in the PRV is equal to the motion averaged D(max) in the (moving) clinical target volume (CTV). We studied margins for the spinal cord in 10 head-and-neck cases and 10 lung cases, each with two different clinical plans. For critical organs with a dose-volume constraint, we also investigated whether a margin recipe was feasible. RESULTS: For the 20 spinal cords considered, the average margin recipe found was: M(R) = 1.6 Sigma + 0.2 sigma with variations for systematic and random errors of 1.2 Sigma to 1.8 Sigma and -0.2 sigma to 0.6 sigma, respectively. The variations were due to differences in shape and position of the dose distributions with respect to the cords. The recipe also depended significantly on the volume definition of D(max). For critical organs with a dose-volume constraint, the PRV concept appears even less useful because a margin around, e.g., the rectum changes the volume in such a manner that dose-volume constraints stop making sense. CONCLUSION: The concept of PRV for planning of radiotherapy is of limited use. Therefore, alternative ways should be developed to include geometric uncertainties of OARs in radiotherapy planning.

Algorithms↗

Fast, daily linac verification for segmented IMRT using electronic portal imaging.

PURPOSE: Intensity modulated radiotherapy (IMRT) requires dedicated quality assurance (QA). Recently, we have published a method for fast (1-2 min) and accurate linac quality control for dynamic multileaf collimation, using a portal imaging device. This method is in routine use for daily leaf motion verification. The purpose of the present study was to develop an equivalent procedure for QA of IMRT with segmented (static) multileaf collimation (SMLC). MATERIALS AND METHODS: The QA procedure is based on measurements performed during 3- to 8-month periods at Elekta, Siemens and Varian accelerators. On each measurement day, images were acquired for a field consisting of five 3 x 22 cm(2) segments. These 10 monitor unit (MU) segments were delivered in SMLC mode, moving the leaves from left to right. Deviations of realized leaf gap widths from the prescribed width were analysed to study the leaf positioning accuracy. To assess hysteresis in leaf positioning, the sequential delivery of the SMLC segments was also inverted. A static 20 x 20 cm(2) field was delivered with exposures between 1 and 50 MU to study the beam output and beam profile at low exposures. Comparisons with an ionisation chamber were made to verify the EPID dose measurements at low MU. Dedicated software was developed to improve the signal-to-noise ratio and to correct for image distortion. RESULTS AND CONCLUSIONS: The observed long-term leaf gap reproducibility (1 standard deviation) was 0.1 mm for the Varian, and 0.2 mm for the Siemens and the Elekta accelerators. In all cases the hysteresis was negligible. Down to the lowest MU, beam output measurements performed with the EPID agreed within 1+/-1% (1SD) with ionisation chamber measurements. These findings led to a fast (3-4 min) procedure for accurate, daily linac quality control for SMLC.

Calibration↗

Fast, multiple optimizations of quadratic dose objective functions in IMRT.

Inverse treatment planning for intensity-modulated radiotherapy may include time consuming, multiple minimizations of an objective function. In this paper, methods are presented to speed up the process of (repeated) minimization of the well-known quadratic dose objective function, extended with a smoothing term that ensures generation of clinically acceptable beam profiles. In between two subsequent optimizations, the voxel-dependent importance factors of the quadratic terms will generally be adjusted, based on an intermediate plan evaluation. The objective function has been written in matrix-vector format, facilitating the use of a recently published, fast quadratic minimization algorithm, instead of commonly applied gradient-based methods. This format also reduces the calculation time in between subsequent minimizations, related to adjustment of the voxel-dependent importance factors. Sparse matrices are used to limit the required amount of computer memory. For three patients, comparisons have been made with a gradient method. Mean speed improvements of up to a factor of 37 have been achieved.

Algorithms↗

Stereotactic body radiation therapy for primary and metastatic liver tumors: A single institution phase i-ii study.

The feasibility, toxicity and tumor response of stereotactic body radiation therapy (SBRT) for treatment of primary and metastastic liver tumors was investigated. From October 2002 until June 2006, 25 patients not suitable for other local treatments were entered in the study. In total 45 lesions were treated, 34 metastases and 11 hepatocellular carcinoma (HCC). Median follow-up was 12.9 months (range 0.5-31). Median lesion size was 3.2 cm (range 0.5-7.2) and median volume 22.2 cm3 (range 1.1-322). Patients with metastases, HCC without cirrhosis, and HCC < 4 cm with cirrhosis were mostly treated with 3 x 12.5 Gy. Patients with HCC > or =4 cm and cirrhosis received 5 x 5 Gy or 3 x 10 Gy. The prescription isodose was 65%. Acute toxicity was scored following the Common Toxicity Criteria and late toxicity with the SOMA/LENT classification. Local failures were observed in two HCC and two metastases. Local control rates at 1 and 2 years for the whole group were 94% and 82%. Acute toxicity grade > or =3 was seen in four patients; one HCC patient with Child B developed a liver failure together with an infection and died (grade 5), two metastases patients presented elevation of gamma glutamyl transferase (grade 3) and another asthenia (grade 3). Late toxicity was observed in one metastases patient who developed a portal hypertension syndrome with melena (grade 3). SBRT was feasible, with acceptable toxicity and encouraging local control. Optimal dose-fractionation schemes for HCC with cirrhosis have to be found. Extreme caution should be used for patients with Child B because of a high toxicity risk.

Adult↗

Automated selection of beam orientations and segmented intensity-modulated radiotherapy (IMRT) for treatment of oesophagus tumors.

BACKGROUND AND PURPOSE: For some treatment sites, there is evidence in the literature that five to nine equi-angular input beam directions are enough for generating IMRT plans. For oesophagus cancer, there is a report showing that going from four to nine beams may even result in lower quality plans. In this paper, our previously published algorithm for automated beam angle selection (Cycle) has been extended to include segmented IMRT. For oesophagus cancer patients, we have investigated whether automated orientation selection from a large number of equi-angular input beam directions (up to thirty-six) for IMRT optimisation can result in improved lung sparing. MATERIALS AND METHODS: CT-data from five oesophagus patients treated recently in our institute were used for this study. For a prescribed mean PTV dose of 55 Gy, Cycle was used in an iterative procedure to minimise the mean lung dose under the following hard constraints: standard deviation for PTV dose inhomogeneity 2% (1,1 Gy), maximum spinal cord dose 45 Gy. Conformal radiotherapy (CFRT) and IMRT plans for a standard four field oesophagus beam configuration were compared with IMRT plans generated by automated selection from nine or thirty-six equi-angular input beam orientations. Comparisons were also made with dose distributions generated with our commercial treatment planning system (TPS), and with observations in the literature. RESULTS: Using Cycle, automated orientation selection from nine or thirty-six input beam directions resulted in improved lung sparing compared to the four field set-ups. Compared to selection from nine input orientations, selection from thirty-six directions did always result in lower mean lung doses, sometimes with even fewer non-zero weight beams. On average only seven beams with a non-zero weight were enough for obtaining the lowest mean lung dose, yielding clinically feasible plans even in case of thirty-six input directions for the optimisation process. With our commercial TPS we observed the same contra-intuitive, unfavourable results as reported in the literature; nine field equi-angular IMRT plans had substantially higher mean lung doses than plans for the conventional four field set-ups. For all cases, the Cycle plans generated from nine equi-angular input directions were superior compared to similar plans generated with our commercial TPS. CONCLUSIONS: For the studied oesophagus cancer patients the best plans for IMRT were obtained with Cycle, using automated beam orientation selection from thirty-six input beam directions. The lowest mean lung doses could be obtained with, on average, a selection of only seven beams with non-zero weight.

Algorithms↗

Application of the No Action Level (NAL) protocol to correct for prostate motion based on electronic portal imaging of implanted markers.

PURPOSE: To evaluate the efficacy of the No Action Level (NAL) off-line correction protocol in the reduction of systematic prostate displacements as determined from electronic portal images (EPI) using implanted markers. METHODS AND MATERIALS: Four platinum markers, two near the apex and two near the base of the prostate, were implanted for localization purposes in patients who received fractionated high dose rate brachytherapy. During the following course of 25 fractions of external beam radiotherapy, the position of each marker relative to the corresponding position in digitally reconstructed radiographs (DRRs) was measured in EPI in 15 patients for on average 17 fractions per patient. These marker positions yield the composite displacements due to both setup error and internal prostate motion, relative to the planning computed tomography scan. As the NAL protocol is highly effective in reducing systematic errors (recurring each fraction) due to setup inaccuracy alone, we investigated its efficacy in reducing systematic composite displacements. The analysis was performed for the center of mass (COM) of the four markers, as well as for the cranial and caudal markers separately. Furthermore, the impact of prostate rotation on the achieved positioning accuracy was determined. RESULTS: In case of no setup corrections, the standard deviations of the systematic composite displacements of the COM were 3-4 mm in the craniocaudal and anterior-posterior directions, and 2 mm in the left-right direction. The corresponding SDs of the random displacements (interfraction fluctuations) were 2-3 mm in each direction. When applying a NAL protocol based on three initial treatment fractions, the SDs of the systematic COM displacements were reduced to 1-2 mm. Displacements at the cranial end of the prostate were slightly larger than at the caudal end, and quantitative analysis showed this originates from left-right axis rotations about the prostate apex. Further analysis revealed that significant time trends are present in these prostate rotations. No significant trends were observed for the prostate translations. CONCLUSIONS: The NAL protocol based on marker positions in EPI halved the composite systematic displacements using only three imaged fractions per patient, and thus allowed for a significant reduction of planning margins. Although large rotations of the prostate, and time trends therein, were observed, the net impact on the measured displacements and on the accuracy obtained with NAL was small.

Algorithms↗

SIFT: a method to verify the IMRT fluence delivered during patient treatment using an electronic portal imaging device.

PURPOSE: Radiotherapy patients are increasingly treated with intensity-modulated radiotherapy (IMRT) and high tumor doses. As part of our quality control program to ensure accurate dose delivery, a new method was investigated that enables the verification of the IMRT fluence delivered during patient treatment using an electronic portal imaging device (EPID), irrespective of changes in patient geometry. METHODS AND MATERIALS: Each IMRT treatment field is split into a static field and a modulated field, which are delivered in sequence. Images are acquired for both fields using an EPID. The portal dose image obtained for the static field is used to determine changes in patient geometry between the planning CT scan and the time of treatment delivery. With knowledge of these changes, the delivered IMRT fluence can be verified using the portal dose image of the modulated field. This method, called split IMRT field technique (SIFT), was validated first for several phantom geometries, followed by clinical implementation for a number of patients treated with IMRT. RESULTS: The split IMRT field technique allows for an accurate verification of the delivered IMRT fluence (generally within 1% [standard deviation]), even if large interfraction changes in patient geometry occur. For interfraction radiological path length changes of 10 cm, deliberately introduced errors in the delivered fluence could still be detected to within 1% accuracy. Application of SIFT requires only a minor increase in treatment time relative to the standard IMRT delivery. CONCLUSIONS: A new technique to verify the delivered IMRT fluence from EPID images, which is independent of changes in the patient geometry, has been developed. SIFT has been clinically implemented for daily verification of IMRT treatment delivery.

Head and Neck Neoplasms↗

Automated beam angle and weight selection in radiotherapy treatment planning applied to pancreas tumors.

PURPOSE: To extend and investigate the clinical value of a recently developed algorithm for automatic beam angle and beam weight selection, for irradiation of pancreas tumors. METHODS AND MATERIALS: The algorithm aims at generation of acceptable treatment plans, i.e., delivering the prescribed tumor dose while strictly obeying the imposed hard constraints for organs at risk and target. Extensions were made to minimize the beam number and/or to escalate the tumor dose. For 5 pancreas patients, the clinical value and the potential for beam number reduction and dose escalation were investigated. Comparisons were made with clinical plans and equiangular plans. RESULTS: Compared to clinical plans, the generated plans with the same number of beams yielded a substantial reduction in the dose to critical tissues. Using the algorithm, an escalated tumor dose of 58 Gy could be achieved for two cases. Maximum dose escalations required a minimum of 3 to 4 beam orientations. For 13 CT slices and an in-slice resolution of 0.5 cm, the total calculation times were 23-55 min, including precalculation of 180 input dose distributions (15 min). CONCLUSIONS: The algorithm yielded acceptable treatment plans with clinically feasible numbers of beams, even for escalated tumor doses. Generated plans were superior to the clinically applied plans and to equiangular setups. Calculation times were clinically acceptable. The algorithm is now increasingly used in clinical routine.

Algorithms↗

Procedures for high precision setup verification and correction of lung cancer patients using CT-simulation and digitally reconstructed radiographs (DRR).

PURPOSE: In a recent study, large systematic setup errors were detected in patients with lung cancer when a conventional simulation procedure was used to define and mark the treatment isocenter. In the present study, we describe a procedure to omit the session at a conventional simulator to remove simulation errors entirely. Isocenter definition and verification was performed at a computed tomography (CT) simulator, and digitally reconstructed radiographs (DRRs) were used for setup verification and correction at the treatment unit. METHODS AND MATERIALS: A CT simulation protocol was developed, in which radiopaque markers were used to verify the coincidence of the isocenter marked on the patients' skin with the isocenter defined in the planning CT scan. This protocol was evaluated for 20 patients. Subsequently, electronic portal images were acquired at the treatment unit. The three-dimensional setup error was established from a template match of the appropriate anatomy visible in two orthogonal beams with the corresponding anatomy in DRRs. An offline setup correction protocol was applied to reduce systematic setup errors. RESULTS: For all patients, the skin marks defined the planning CT scan isocenter to within +/- 1.5 mm in each of the three main directions. Random setup errors at the treatment unit were 1.8, 2.0, and 1.9 mm (1 SD) for the lateral (x), the superior-inferior (y), and the anterior-posterior (z) directions, respectively. With the use of the correction protocol, the systematic errors for x, y, and z were 1.5, 1.5, and 1.3 mm (1 SD). CONCLUSIONS: Because the distributions of treatment setup errors measured against DRRs obtained in our CT simulation were equal to previously obtained distributions measured against simulator films, conventional simulation can be omitted and DRRs are well-suited for setup verification. By adopting our CT simulation procedure, the large systematic simulation setup errors, which remain hidden if a conventional simulation is performed, can be avoided.

Carcinoma, Non-Small-Cell Lung↗

Two-dimensional measurement of photon beam attenuation by the treatment couch and immobilization devices using an electronic portal imaging device.

In our institution, an individualized dosimetric quality assurance protocol for intensity modulated radiotherapy (IMRT) is being implemented. This protocol includes dosimetric measurements with a fluoroscopic electronic portal imaging device (EPID) for all IMRT fields while the patient is being irradiated. For some of the first patients enrolled in this protocol, significant beam attenuation by (carbon fiber) components of the treatment couch was observed. To study this beam attenuation in two-dimensional, EPID images were also acquired in absence of the patient, both with and without treatment couch and immobilization devices, as positioned during treatment. For treatments of head and neck cancer patients with a 6 MV photon beam, attenuation of up to 15% was detected. These findings led to the development of new tools and procedures for planning and treatment delivery to avoid underdosages in the tumor.

Artifacts↗

Detection of organ movement in cervix cancer patients using a fluoroscopic electronic portal imaging device and radiopaque markers.

PURPOSE: To investigate the use of a fluoroscopic electronic portal imaging device (EPID) and radiopaque markers to detect internal cervix movement. METHODS AND MATERIALS: For 10 patients with radiopaque markers clamped to the cervix, electronic portal images were made during external beam irradiation. Bony structures and markers in the portal images were registered with the same structures in the corresponding digitally reconstructed radiographs of the planning computed tomogram. RESULTS: The visibility of the markers in the portal images was good, but their fixation should be improved. Generally, the correlation between bony structure displacements and marker movement was poor, the latter being substantially larger. The standard deviations describing the systematic and random bony anatomy displacements were 1.2 and 2.6 mm, 1.7 and 2.9 mm, and 1.6 and 2.7 mm in the lateral, cranial-caudal, and dorsal-ventral directions, respectively. For the marker movement those values were 3.4 and 3.4 mm, 4.3 and 5.2 mm, 3.2 and 5.2 mm, respectively. Estimated clinical target volume to planning target volume (CTV-PTV) planning margins (approximately 11 mm) based on the observed overall marker displacements (bony anatomy + internal cervix movement) are only marginally larger than the margins required to account for internal marker movement alone. CONCLUSIONS: With our current patient setup techniques and methods of setup verification and correction, the required CTV-PTV margins are almost fully determined by internal organ motion. Setup verification and correction using radiopaque markers might allow decreasing those margins, but technical improvements are needed.

Cervix Uteri↗

Geometrical uncertainties, radiotherapy planning margins, and the ICRU-62 report.

In this paper, we elaborate on the proposals in the ICRU-62 report concerning planning target volume (PTV) margins for geometrical uncertainties during radiotherapy, such as variations in patient set-up and internal organ motion. According to the ICRU, these margins should be such that the planned dose in the PTV is representative of the real dose in the 'moving' clinical target volume (CTV). We demonstrate that the dosimetrical consequences of systematic and random geometrical uncertainties are fundamentally different, which should be reflected in margin calculations. The recommendation in the ICRU-62 report, to quadratically add standard deviations for systematic (Sigma(tot)) and random (sigma(tot)) errors to determine an overall standard deviation for margin calculations, is therefore generally not valid. Instead, a previously published recipe for PTV margin calculation, M = 2Sigma(tot) + 0.7sigma(tot), does indeed account for the different impact of systematic and random errors on the dose in the CTV. If, for both random and systematic uncertainties, the internal and external errors are uncorrelated and quantified by the standard deviations sigma(int), sigma(ext), Sigma(int), Sigma(ext), then Sigma(tot) = square root (Sigma(int)(2) + Sigma(ext)(2)) and sigma(tot) = square root (sigma(int)(2) + sigma(ext)(2)). If the PTV margin thus acquired is deliberately reduced to spare normal tissues, the planned PTV dose is not representative of the CTV anymore. Therefore, we recommend to also report the minimum dose in the volume originally defined by the recipe (designated RTV, i.e. representative target volume).

Humans↗

Portal dose image prediction for dosimetric treatment verification in radiotherapy. II. An algorithm for wedged beams.

A method is presented for calculation of a two-dimensional function, T(wedge)(x,y), describing the transmission of a wedged photon beam through a patient. This in an extension of the method that we have published for open (nonwedged) fields [Med. Phys. 25, 830-840 (1998)]. Transmission functions for open fields are being used in our clinic for prediction of portal dose images (PDI, i.e., a dose distribution behind the patient in a plane normal to the beam axis), which are compared with PDIs measured with an electronic portal imaging device (EPID). The calculations are based on the planning CT scan of the patient and on the irradiation geometry as determined in the treatment planning process. Input data for the developed algorithm for wedged beams are derived from (the already available) measured input data set for transmission prediction in open beams, which is extended with only a limited set of measurements in the wedged beam. The method has been tested for a PDI plane at 160 cm from the focus, in agreement with the applied focus-to-detector distance of our fluoroscopic EPIDs. For low and high energy photon beams (6 and 23 MV) good agreement (approximately 1%) has been found between calculated and measured transmissions for a slab and a thorax phantom.

Algorithms↗

Fast and accurate leaf verification for dynamic multileaf collimation using an electronic portal imaging device.

A prerequisite for accurate dose delivery of IMRT profiles produced with dynamic multileaf collimation (DMLC) is highly accurate leaf positioning. In our institution, leaf verification for DMLC was initially done with film and ionization chamber. To overcome the limitations of these methods, a fast, accurate and two-dimensional method for daily leaf verification, using our CCD-camera based electronic portal imaging device (EPID), has been developed. This method is based on a flat field produced with a 0.5 cm wide sliding gap for each leaf pair. Deviations in gap widths are detected as deviations in gray scale value profiles derived from the EPID images, and not by directly assessing leaf positions in the images. Dedicated software was developed to reduce the noise level in the low signal images produced with the narrow gaps. The accuracy of this quality assurance procedure was tested by introducing known leaf position errors. It was shown that errors in leaf gap as small as 0.01-0.02 cm could be detected, which is certainly adequate to guarantee accurate dose delivery of DMLC treatments, even for strongly modulated beam profiles. Using this method, it was demonstrated that both short and long term reproducibility in leaf positioning were within 0.01 cm (1sigma) for all gantry angles, and that the effect of gravity was negligible.

Calibration↗