| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Full paper |
1 Department of Radiation Oncology, University Hospital of Heidelberg and German Cancer Research Center, Heidelberg, 2 Department of Medical Physics, German Cancer Research Center, Heidelberg, Germany
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
There are different approaches to image guidance aimed at the reduction of uncertainties of target position. For acquisition of 3D anatomic information in the treatment position, megavoltage CT [3] and cone beam kilovoltage CT fixed at the gantry [4] or a separate in-room CT scanner [57] has been tested so far. When using an in-room CT, the system consists of a linear accelerator (linac) and a conventional CT scanner connected via a conventional treatment couch. The CT scanner is mounted on rails. The obvious advantage of this system is that all components are separately established for clinical application. The accuracy has been demonstrated by Cheng and coworkers [8].
We would like to point out three particularities of the system. (1) The patient is immobilized in the usual way on the treatment couch and CT scanning is performed for confirmation of the correct isocentre and patient (target) position. By using a gantry mounted on rails, scanning is performed without either couch or patient movement, eliminating reductions in accuracy. After performing the scan, the treatment couch is rotated to the linac side for irradiation. (2) The kilovoltage technique allows an image quality with high soft-tissue contrast. This is a requirement for detecting deviations of soft tissue targets from their position at the initial CT scan used for treatment planning, even when bony landmarks are repositioned correctly. (3) The CT data set can be used without any transformation for treatment planning and dose calculation. The basic input for radiotherapy planning systems considering tissue heterogeneities is the relationship between CT Hounsfield units and electron densities, which is determined using common CT-calibration methods. Thus the in-room CT is well suited for target point verification, correction of setup errors and interfraction target deviations due to organ motion, as well as for recalculation of the dose actually given. The simplest correction is to correct the target point without changing the treatment plan. Beyond this, systematic changes of the patient's anatomy under therapy (e.g. weight loss, tumour mass reduction) can be considered and, if necessary, the radiation treatment plan can be re-optimized until the next fraction. The highest level of adaptation achievable with a linac-CT combination is the adaptation/optimization of the treatment plan to the actual given situation, i.e. between the CT scan and the linac irradiation. Since a CT scan is not possible in treatment position, the described system cannot be used to detect intrafraction motion.
The purpose of this study is to demonstrate the potential of an integrated linac-CT scanner system for different tumour sites and to derive a workflow which is reasonable for clinical routine with target point correction and re-optimization of the dose distribution.
| Materials and methods |
|---|
|
|
|---|
|
Treatment
For all single-fractionated treatments, a control CT scan is performed directly prior to irradiation. By comparing the control CT scan to the CT scan made for treatment planning, any translational displacement is detected and the target point corrected accordingly.
For fractionated treatment of intracranial targets, the repositioning accuracy of the patient's head with the customized head mask was found to be better than 2 mm in all three dimensions [9]. We usually therefore include a 2 mm safety margin in the planning target volume (PTV), and apart from the verification of the target point on the first fraction, we do not perform further control CT scans.
Fractionated treatment of extracranial targets is accompanied by frequent CT control scans during the treatment course, since even with a rigid fixation significant interfractional discrepancies are observed (see results below). The scan frequency depends on the individual case and is between once weekly for unproblematic cases, up to daily CT scans for patients with large repositioning errors and critical proximity of target structures to organs at risk.
Analysis of CT scans
All CT scans generated with the Siemens Primatom are performed with stereotactic localizers attached to the fixation frame in exactly the same way as they were attached in the treatment planning CT scan. This defines a frame-based coordinate system independent of the patient's actual position, which is a prerequisite for any further analysis. For the first fraction, lead ball bearings (BB) were additionally put onto the laser adjustment lines for target point verification, so no further portal imaging was necessary.
When we started using the Siemens Primatom, the correct positioning of the patient could only be verified by manually highlighting representative anatomical landmarks on selected CT slices. For example, about three bony landmarks on the isocentre plane were pre-defined on the planning CT in the treatment planning system and then located on the control CT using the CT scanner console. The average coordinate difference between the scans then gave an estimation of the average displacement error. Although technically feasible, this procedure was quite time consuming and movements visible on only some CT slices could be missed.
We therefore developed a new workflow to both reduce the time requirements and improve the accuracy of the analysis, and integrated it into our in-house developed radiotherapy software environment. The control CT scan is transferred via a network to our treatment planning system VIRTUOS. Inside the treatment planning system, the localizers of the control CT and the planning CT are stereotactically correlated, making the coordinates in both cubes directly comparable. The region of the control CT cube containing the target volume is then automatically matched onto the planning CT. The transformation determined in this step immediately gives the current displacement error. At the moment, we use two rigid matching algorithms: Rigid correlation matching (RCM) based on bony anatomy considers only translational movements and provides the target point correction vector, and mutual information matching (MIM) for bony and soft structures considers translational and rotational movements. Both matching algorithms derive the transformation vector T between planning volume A and transformed control volume BT by searching for maximal similarity in the region of interest
, measured by the correlation coefficient CC and the mutual information MI, respectively.
The correlation coefficient between A and BT is calculated by:
|
|
is the mean value of A and A three-dimensional Fourier transformation of A and BT into the frequency space using a fast Fourier transformation (FFT) transforms the convolution in Equation (1) into a simple multiplication. After inverse transformation, we obtain the correlation volume that holds the CC values for all possible 3D translation vectors, and the global maximum can easily be found. Since we use discrete Fourier transformation, the precision is restricted to the voxel dimension. The RCM analysis of a typical CT dataset (512x512 points, 30 slices) by a standard PC takes less than 30 s. RCM matching was performed for the complete region around the target volume and separately for the upper und lower quarter. Significant differences between the three results indicate that the deviation is more complex than a pure translation.
The similarity between volume A and volume BT measured by mutual information is given by:
|
|
RCM is significantly faster than MIM, always finds the global optimum and gives a transformation that can be corrected for by a simple shift of the target point. We perform RCM matching of the upper and lower quarter of the region of interest and the MIM matching only to check for deviations that cannot be corrected for by translation only. In any case, the user has to visually check that the resulting match is correct. With standard tools, such as red-green overlay and checkerboard provided by VIRTUOS, this manual check usually takes less than 1 min.
Adaptive radiotherapy
In adaptive radiotherapy, deviations of patient repositioning and anatomy from the initial planning CT scan are detected and, if necessary, corrected for. This is in contrast to conventional static radiotherapy where the plan is based on one single CT scan and delivered to the patient throughout all fractions. Different levels of adaptation can be defined:
In level 0, the treatment plan is based on several CT scans rather than a single one, and the information about statistical movements of the target and organs at risk are integrated into the treatment plan. This can be accomplished by defining patient-specific or at least site-specific safety margins [13] or by including statistical methods into the inverse optimization process [14, 15].
Level 1 corrects for errors offline, i.e. the correction is done after several CT scans are gathered. Systematic interfractional errors can be detected and corrected. If the main deviation is a translation, usually it is sufficient to shift the target point and leave the plan unchanged (level 1A). If the error is more complex, the contours have to be adapted to the new geometry, and a new plan has to be generated by re-running the optimization of the inverse planning program (level 1B).
Level 2 uses the same methods as level 1, but here the correction is placed in between the CT scan and the directly-following linac irradiation. By analogy with level 1 definitions, level 2A stands for shifting the target point and 2B for generating a new plan. This way not only the systematic error, but also daily random interfractional errors can be corrected. Since the patient remains in the fixation device during the adaptation procedure, the time requirements for this step are more critical than for the lower adaptation levels.
Level 3 eventually also takes into account intrafractional errors. This task requires an imaging device operating during the irradiation itself and cannot be accomplished by a linac-CT scanner combination.
In the following discussion we will focus on our experiences with Primatom-based radiotherapy for prostate, paraspinal and head and neck cases. However, the general aspects and conclusion apply to other indications as well.
| Results and discussion |
|---|
|
|
|---|
Prostate cancer
Patients with prostate cancer treated at our institution are immobilized by a wrap-around body cast and a head mask. Usually they receive weekly control CT scans during the treatment course. Figure 2
shows the displacements of the bony structures around the prostate in the control scans relative to their position in the planning CT scan for 10 prostate patients (P1P10) as they were calculated by rigid correlation matching. The symbols (
, ,
) indicate the mean displacement values for each coordinate in x, y and z direction, with the range represented by the standard deviation of all values observed. Checks with RCM for the upper and lower quarter of the target volume separately showed differences below 1 mm, and checks with mutual information matching confirmed that the rotational error was negligibly small with only around 0.5° (maximum 1°) for all three rotational axes. As can be seen in Figure 2
, even without corrections of the target point the repositioning accuracy of the bony anatomy was very good, with displacements below 3 mm in almost all cases. The correct positioning of the bony anatomy alone already leads to good repositioning of the prostate itself, as is indicated by the results of Beard et al [16] who reported deviations of the prostate above 1 cm in 3% of all cases after bony matching, whereas Wong et al [7] reported deviations above 1 cm in 15% of all cases without bony matching. Manually checking the soft tissues also showed that in our data the displacements of the prostate itself were bigger than the bony variations. In a separate study, the positions of the prostate, bladder and rectum were statistically analysed and taken into account by adding a margin to the gross tumour volume (GTV) and clinical target volume (CTV) to obtain the planning target volume (PTV) (adaptation level 0) [13, 17]. To reduce these margins and further improve the accuracy of prostate treatments, we are currently working on matching algorithms that will enable the automatic correction of interfractional displacements of the prostate itself on adaptation level 2. Deformation of the shape of the prostate and seminal vesicles is reported to be small relative to organ motion [18], so simply shifting the target point (level 2A) will probably be sufficient.
|
) and craniocaudal (z:
) direction the typical displacement was around 3 mm. However, in these directions maximal errors above 10 mm were observed. Lateral displacements especially can become critical for these patients, as it is exemplarily shown in Figure 4
|
|
Head and neck tumours
The target volume for head and neck tumours regularly includes the base of skull and extends to the upper thoracic aperture. The patients are fixated with a head mask and a vacuum pillow. The cranial part inside the head mask is very accurately repositioned during the whole treatment course. In contrast, the lower extracranial part shows more variations. The result is a complex deformation of the target volume that cannot be described by a translation and cannot be corrected easily by shifting the target point without changing the treatment plan.
Here we show an exemplary case of a patient who was treated for teratocarcinosarcoma of the paranasal sinuses. During the treatment course frequent control CT scans were performed, and while the intracranial part was accurately positioned throughout all fractions, the lower, extracranial part of the body was systematically shifted approximately 1.5 cm along the y-axis from the middle of the treatment course onwards, see Figure 5
comparing the planning CT with the control CT of fraction 20. Note that the contours in Figure 5
refer to the planning CT and do not fit to the anatomical situation at fraction 20. We therefore re-drew the organ contours for fraction 20 and re-calculated the dose of the original plan on the CT of fraction 20. The results are shown in the form of dosevolume histograms (DVHs) in Figure 6
. Figure 6a
is the original plan based on the planning CT scan and Figure 6b
is the same plan applied to the situation at fraction 20 with updated organ contours. One can see that especially the coverage of the lower target volume has significantly deteriorated. The shift of the extracranial part of the body was systematically seen on three successive control CTs, so we decided to adapt the plan to the new geometry by re-running the optimization in our inverse planning program (KonRad by Siemens OCS) for the control CT scan with the new contours, but without changing the original setting of dose constraints and weighting factors. The new plan was calculated in approximately 2 min, and as shown in Figure 6c
the adapted plan resembles the original plan much better than the uncorrected one. Because the contours had to be re-drawn manually, at the moment this procedure takes too much time to fit in between the CT scan and the directly-following irradiation to correct for random interfractional setup errors. However, it allows for a very good adaptation to complex systematic variations occurring during the treatment course (adaptation level 1B).
|
|
Patients with paraspinal tumours can greatly benefit from the higher accuracy of the treatment. Target point correction seems sufficient; elastic deformation or plan re-optimization was not necessary for the patients we treated. For prostate patients, the bony anatomy is already quite precisely repositioned due to our rigid immobilization device. Here we expect further improvements by elastic matching algorithms that automatically detect the position of the prostate itself (based on Primatom CT scans) and adapt the plan either by target point correction or re-optimization. Also, for patients with head and neck tumours the quality of the radiotherapy could be significantly improved. For these cases, plan re-optimization appeared to be more important than for other tumour sites due to the complex nature of the interfractional deformations.
The linac-CT scanner combination already meets all hardware requirements to completely eliminate interfractional setup errors from the treatment course. However, a fast and robust workflow is necessary for widespread use in clinical practice, and to accomplish this further development of algorithms and software tools (e.g. automatic elastic matching) is needed. Real-time plan adaptation to elastic deformations in particular is under current investigation [20] and not in clinical practice yet.
Concerning the documentation of an adapted radiation treatment course, the most desirable final record would be a treatment plan where the doses to each volume element delivered throughout the course are superimposed, resulting in concise dose statistics for each structure and a single DVH for the complete treatment. However, this requires the tracking of each voxel throughout all CT scans. Simply averaging the DVHs of the single fractions to obtain a final DVH would lead to erroneous results since, in a DVH, the spatial information is lost. Techniques for tracking the voxels are under current development; at the moment, each modified (adapted) treatment plan is documented separately (as in Figure 6a,c
).
Our patient immobilization device with customized wrap-around body casts alone leads to quite high repositioning accuracy, but is also quite labour and time intensive to build. When using adaptive radiotherapy strategies, the fixation can probably be made less sophisticated, lowering the overall workload and further strengthening the role of a combined imaging/treatment device such as the linac-CT scanner combination.
| Conclusions |
|---|
|
|
|---|
Received for publication August 22, 2005. Revision received January 6, 2006. Accepted for publication January 25, 2006.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
R M HARRISON Doses to organs and tissues from concomitant imaging in radiotherapy: a suggested framework for clinical justification Br. J. Radiol., December 1, 2008; 81(972): 970 - 974. [Abstract] [Full Text] [PDF] |
||||
![]() |
R M Harrison, M Wilkinson, D J Rawlings, and M Moore Doses to critical organs following radiotherapy and concomitant imaging of the larynx and breast Br. J. Radiol., December 1, 2007; 80(960): 989 - 995. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| BJR | DMFR | IMAGING | ALL BIR JOURNALS |