Search PubMed⌕ Search

Biomedical subjects

M van Herk

Publications and source records attributed to M van Herk.

At least 19 recordsLinked to original sources

Comparison of ghosting effects for three commercial a-Si EPIDs.

Many studies have reported dosimetric characteristics of amorphous silicon electronic portal imaging devices (EPIDs). Some studies ascribed a non-linear signal to gain ghosting and image lag. Other reports, however, state the effect is negligible. This study compares the signal-to-monitor unit (MU) ratio for three different brands of EPID systems. The signal was measured for a wide range of monitor units (5-1000), dose-rates, and beam energies. All EPIDs exhibited a relative under-response for beams of few MUs; giving 4 to 10% lower signal-to-MU ratios relative to that of 1000 MUs. This under-response is consistent with ghosting effects due to charge trapping.

Calibration↗

Clinical experience with EPID dosimetry for prostate IMRT pre-treatment dose verification.

The aim of this study was to demonstrate how dosimetry with an amorphous silicon electronic portal imaging device (a-Si EPID) replaced film and ionization chamber measurements for routine pre-treatment dosimetry in our clinic. Furthermore, we described how EPID dosimetry was used to solve a clinical problem. IMRT prostate plans were delivered to a homogeneous slab phantom. EPID transit images were acquired for each segment. A previously developed in-house back-projection algorithm was used to reconstruct the dose distribution in the phantom mid-plane (intersecting the isocenter). Segment dose images were summed to obtain an EPID mid-plane dose image for each field. Fields were compared using profiles and in two dimensions with the y evaluation (criteria: 3%/3 mm). To quantify results, the average gamma (gamma avg), maximum gamma (gamma max), and the percentage of points with gamma < 1(P gamma < 1) were calculated within the 20% isodose line of each field. For 10 patient plans, all fields were measured with EPID and film at gantry set to 0 degrees. The film was located in the phantom coronal mid-plane (10 cm depth), and compared with the back-projected EPID mid-plane absolute dose. EPID and film measurements agreed well for all 50 fields, with (gamma avg) =0.16, (gamma max)=1.00, and (P gamma < 1)= 100%. Based on these results, film measurements were discontinued for verification of prostate IMRT plans. For 20 patient plans, the dose distribution was re-calculated with the phantom CT scan and delivered to the phantom with the original gantry angles. The planned isocenter dose (plan(iso)) was verified with the EPID (EPID(iso)) and an ionization chamber (IC(iso)). The average ratio, (EPID(iso)/IC(iso)), was 1.00 (0.01 SD). Both measurements were systematically lower than planned, with (EPID(iso)/plan(iso)) and (IC(iso)/plan(iso))=0.99 (0.01 SD). EPID mid-plane dose images for each field were also compared with the corresponding plane derived from the three dimensional (3D) dose grid calculated with the phantom CT scan. Comparisons of 100 fields yielded (gamma avg)=0.39, gamma max=2.52, and (P gamma < 1)=98.7%. Seven plans revealed under-dosage in individual fields ranging from 5% to 16%, occurring at small regions of overlapping segments or along the junction of abutting segments (tongue-and-groove side). Test fields were designed to simulate errors and gave similar results. The agreement was improved after adjusting an incorrectly set tongue-and-groove width parameter in the treatment planning system (TPS), reducing (gamma max) from 2.19 to 0.80 for the test field. Mid-plane dose distributions determined with the EPID were consistent with film measurements in a slab phantom for all IMRT fields. Isocenter doses of the total plan measured with an EPID and an ionization chamber also agreed. The EPID can therefore replace these dosimetry devices for field-by-field and isocenter IMRT pre-treatment verification. Systematic errors were detected using EPID dosimetry, resulting in the adjustment of a TPS parameter and alteration of two clinical patient plans. One set of EPID measurements (i.e., one open and transit image acquired for each segment of the plan) is sufficient to check each IMRT plan field-by-field and at the isocenter, making it a useful, efficient, and accurate dosimetric tool.

Algorithms↗

Fusion of respiration-correlated PET and CT scans: correlated lung tumour motion in anatomical and functional scans.

Lower lobe lung tumours in particular can move up to 2 cm in the cranio-caudal direction during the respiration cycle. This breathing motion causes image artefacts in conventional free-breathing computed tomography (CT) and positron emission tomography (PET) scanning, rendering delineation of structures for radiotherapy inaccurate. The purpose of this study was to develop a method for four-dimensional (4D) respiration-correlated (RC) acquisition of both CT and PET scans and to develop a framework to fuse these modalities. The breathing signal was acquired using a thermometer in the breathing airflow of the patient. Using this breathing signal, the acquired CT and PET data were grouped to the corresponding respiratory phases, thereby obtaining 4D CT and PET scans. Tumour motion curves were assessed in both image modalities. From these tumour motion curves, the deviation with respect to the mean tumour position was calculated for each phase. The absolute position of the centre of the tumour, relative to the bony anatomy, in the RCCT and gated PET scans was determined. This 4D acquisition and 4D fusion methodology was performed for five patients with lower lobe tumours. The peak-to-peak amplitude range in this sample group was 1-2 cm. The 3D tumour motion curve differed less than 1 mm between PET and CT for all phases. The mean difference in amplitude was less than 1 mm. The position of the centre of the tumour (relative to the bony anatomy) in the RCCT and gated PET scan was similar (difference <1 mm) when no atelectasis was present. Based on these results, we conclude that the method described in this study allows for accurate quantification of tumour motion in CT and PET scans and yields accurate respiration-correlated 4D anatomical and functional information on the tumour region.

Algorithms↗

Three-dimensional dose reconstruction of breast cancer treatment using portal imaging.

In this study, we present an algorithm for three-dimensional (3-D) dose reconstruction using portal images obtained with an electronic portal imaging device (EPID). For this purpose an algorithm for 2-D dose reconstruction, which was previously developed in our institution, was adapted. The external contour of the patient was used to correct for absorption of primary photons, but the presence of inhomogeneities was not taken into account. The accuracy of the algorithm was determined by irradiating two anthropomorphic breast phantoms with 6 MV photons. The dose values derived from portal images were compared with results from 3-D dose calculations, which, in turn, were verified with data obtained with an ionization chamber and film dosimetry. It was found that the application of contour information significantly improves the accuracy of 2-D dose reconstruction. If the total dose at the isocenter plane resulting from all treatment beams is reconstructed, the average deviation from the planned dose is 0.1%+/-1.7% (1 SD). If contour information is not available, the differences increase up to +/-20% for the individual beams. In that case, the dose can only be reconstructed with reasonable accuracy when (nearly) opposing beams are used. The average deviation of the 3-D reconstructed dose from the planned dose in the irradiated volume is 1.4%+/-5.4% (1 SD). If the irradiated volume is enclosed by planes less than 5 cm distant from the isocenter plane, then the average deviation is only 0.5%+/-3.4% (1 SD). It can be concluded that the proposed algorithm for a 3-D dose reconstruction allows a determination of the dose at the isocenter plane and the dose-volume histogram with an accuracy acceptable for an independent verification of the treatment.

Algorithms↗

Field size reduction enables iso-NTCP escalation of tumor control probability for irradiation of lung tumors.

PURPOSE: With the mean lung dose (MLD) as an estimator for the normal tissue complication probability (NTCP) of the lung, we assessed whether the probability of tumor control of lung tumors might be increased by dose escalation in combination with a reduction of field sizes, thus increasing target dose inhomogeneity while maintaining a constant MLD. METHODS AND MATERIALS: An 8-MV AP-PA irradiation of a lung tumor, located in a cylindrically symmetric lung-equivalent phantom, was modeled using numerical simulation. Movement of the clinical target volume (CTV) due to patient breathing and setup errors was simulated. The probability of tumor control, expressed as the equivalent uniform dose (EUD) of the CTV, was assessed as a function of field size, under the constraint of a constant MLD. The approach was tested for a treatment of a non-small cell lung cancer (NSCLC) patient using the beam directions of the clinically applied treatment plan. RESULTS: In the phantom simulation it was shown that by choosing field sizes that ensured a minimum dose of 95% in the CTV ("conventional" plan) taking into account setup errors and tumor motion, an EUD of the CTV of 43.8 Gy can be obtained for a prescribed dose of 44.2 Gy. By reducing the field size and thus shifting the 95% isodose surface inwards, the EUD increases to a maximum of 68.3 Gy with a minimum dose in the CTV of 55.2 Gy. This increase in EUD is caused by the fact that field size reduction enables escalation of the prescribed dose while maintaining a constant MLD. Further reduction of the field size results in decrease of the EUD because the minimum dose in the CTV becomes so low that it has a predominant effect on the EUD, despite further escalation of the prescribed dose. For the NSCLC patient, the EUD could be increased from an initial 62.2 Gy for the conventional plan, to 83.2 Gy at maximum. In this maximum, the prescribed dose is 88.1 Gy, and the minimum dose in the CTV is 67.4 Gy. In this case, the 95% isodose surface is conformed closely to the "static" CTV during treatment planning. CONCLUSIONS: Iso-NTCP escalation of the probability of tumor control is possible for lung tumors by reducing field sizes and allowing a larger dose inhomogeneity in the CTV. Optimum field sizes can be derived, having the highest EUD and highest minimum dose in the CTV under condition of a constant NTCP of the lungs. We conclude that the concept of homogeneous dose in the target volume is not the best approach to reach the highest probability of tumor control for lung tumors.

Humans↗

Accurate measurement of the dynamic response of a scanning electronic portal imaging device.

An important condition for the safe introduction of dynamic intensity modulated radiotherapy (IMRT) using a multileaf collimator (MLC) is the ability to verify the leaf trajectories. In order to verify IMRT using an electronic portal imaging device (EPID), the EPID response should be accurate and fast. Noninstantaneous dynamic response causes motion blurring. The aim of this study is to develop a measurement method to determine the magnitude of the geometrical error as a result of motion blurring for imagers with scanning readout. The response of a liquid-filled ionization chamber EPID, as an example of a scanning imager, on a moving beam is compared with the response of a diode placed at the surface of the EPID. The signals are compared under the assumption that all EPID rows measure the same dose rate when a straight moving field edge is imaged. The measurements are performed at several levels of attenuation to investigate the influence of dose rate on the response of the detector. The accuracy of the measurement method is better than 0.25 mm. We found that the liquid-filled ionization chamber EPID does not suffer from significant motion blurring under clinical circumstances. Using a maximum gradient edge detector to determine the field edge in an image obtained by a liquid-filled ionization chamber EPID, errors smaller than 1 mm are found at a dose rate of 105 MU/min and a field edge speed of 1.1 cm/s. The errors reduce at higher dose rates. The presented method is capable of quantifying the geometrical errors in determining the position of the edge of a moving field with subpixel accuracy. The errors in field edge position determined by a liquid-filled ionization chamber EPID are negligible in clinical practice. Consequently, these EPIDs are suitable for geometric IMRT verification, as far as dynamic response is concerned.

Image Processing, Computer-Assisted↗

Physical aspects of a real-time tumor-tracking system for gated radiotherapy.

PURPOSE: To reduce uncertainty due to setup error and organ motion during radiotherapy of tumors in or near the lung, by means of real-time tumor tracking and gating of a linear accelerator. METHODS AND MATERIALS: The real-time tumor-tracking system consists of four sets of diagnostic X-ray television systems (two of which offer an unobstructed view of the patient at any time), an image processor unit, a gating control unit, and an image display unit. The system recognizes the position of a 2.0-mm gold marker in the human body 30 times per second using two X-ray television systems. The marker is inserted in or near the tumor using image guided implantation. The linear accelerator is gated to irradiate the tumor only when the marker is within a given tolerance from its planned coordinates relative to the isocenter. The accuracy of the system and the additional dose due to the diagnostic X-ray were examined in a phantom, and the geometric performance of the system was evaluated in 4 patients. RESULTS: The phantom experiment demonstrated that the geometric accuracy of the tumor-tracking system is better than 1.5 mm for moving targets up to a speed of 40 mm/s. The dose due to the diagnostic X-ray monitoring ranged from 0.01% to 1% of the target dose for a 2.0-Gy irradiation of a chest phantom. In 4 patients with lung cancer, the range of the coordinates of the tumor marker during irradiation was 2.5-5.3 mm, which would have been 9.6-38.4 mm without tracking. CONCLUSION: We successfully implemented and applied a tumor-tracking and gating system. The system significantly improves the accuracy of irradiation of targets in motion at the expense of an acceptable amount of diagnostic X-ray exposure.

Algorithms↗

The probability of correct target dosage: dose-population histograms for deriving treatment margins in radiotherapy.

PURPOSE: To provide an analytical description of the effect of random and systematic geometrical deviations on the target dose in radiotherapy and to derive margin rules. METHODS AND MATERIALS: The cumulative dose distribution delivered to the clinical target volume (CTV) is expressed analytically. Geometrical deviations are separated into treatment execution (random) and treatment preparation (systematic) variations. The analysis relates each possible preparation (systematic) error to the dose distribution over the CTV and allows computation of the probability distribution of, for instance, the minimum dose delivered to the CTV. RESULTS: The probability distributions of the cumulative dose over a population of patients are called dose-population histograms in short. Large execution (random) variations lead to CTV underdosage for a large number of patients, while the same level of preparation (systematic) errors leads to a much larger underdosage for some of the patients. A single point on the histogram gives a simple "margin recipe." For example, to ensure a minimum dose to the CTV of 95% for 90% of the patients, a margin between CTV and planning target volume (PTV) is required of 2.5 times the total standard deviation (SD) of preparation (systematic) errors (Sigma) plus 1.64 times the total SD of execution (random) errors (sigma') combined with the penumbra width, minus 1.64 times the SD describing the penumbra width (sigma(p)). For a sigma(p) of 3.2 mm, this recipe can be simplified to 2.5 Sigma + 0.7 sigma'. Because this margin excludes rotational errors and shape deviations, it must be considered as a lower limit for safe radiotherapy. CONCLUSION: Dose-population histograms provide insight into the effects of geometrical deviations on a population of patients. Using a dose-probability based approach, simple algorithms for choosing margins were derived.

Algorithms↗

3-D portal image analysis in clinical practice: an evaluation of 2-D and 3-D analysis techniques as applied to 30 prostate cancer patients.

PURPOSE: To investigate the clinical importance and feasibility of a 3-D portal image analysis method in comparison with a standard 2-D portal image analysis method for pelvic irradiation techniques. METHODS AND MATERIALS: In this study, images of 30 patients who were treated for prostate cancer were used. A total of 837 imaged fields were analyzed by a single technologist, using automatic 2-D and 3-D techniques independently. Standard deviations (SDs) of the random, systematic, and overall variations, and the overall mean were calculated for the resulting data sets (2-D and 3-D), in the three principal directions (left-right [L-R], cranial-caudal [C-C], anterior-posterior [A-P]). The 3-D analysis included rotations as well. For the translational differences between the three data sets, the overall SD and overall mean were computed. The influence of out-of-plane rotations on the 2-D registration accuracy was determined by analyzing the difference between the 2-D and 3-D translation data as function of rotations. To assess the reliability of the 2-D and 3-D methods, the number of times the automatic match was manually adjusted was counted. Finally, an estimate of the workload was made. RESULTS: The SDs of the random and systematic components of the rotations around the three orthogonal axes were 1. 1 (L-R), 0.6 (C-C), 0.5 (A-P) and 0.9 (L-R), 0.6 (C-C), 0.8 (A-P) degrees, respectively. The overall mean rotation around the L-R axis was 0.7 degrees, which deviated significantly from zero. Translational setup errors were comparable for 2-D and 3-D analysis (ranging from 1.4 to 2.2 mm SD and from 1.5 to 2.5 mm SD, respectively). The variation of the difference between the 2-D and 3-D translation data increased from 1.1 mm (SD) for zero rotations to 2.7 mm (SD) for out-of-plane rotations of 3 degrees, due to a reduced 2-D registration accuracy for large rotations. The number of times the analysis was not considered acceptable and was manually adjusted was 44% for the 2-D analysis, and 6% for the 3-D analysis. CONCLUSION: True 3-D analysis of setup errors for a group of 30 patients with prostate cancer demonstrated that setup rotations are rather small. The deformation of the projected anatomy in portal images caused by out-of-plane rotations leads to a reduced 2-D registration accuracy. For rotations larger than 3 degrees this effect can be quite pronounced, making 3-D registration the preferred method. Furthermore, the automatic 3-D registration has a higher success rate, most likely because this technique uses more information compared to the 2-D method.

Feasibility Studies↗

Target volumes in radiotherapy for high-grade malignant glioma of the brain.

Delineation of the clinical target volume (CTV) in radiation treatment planning of high-grade glioma is a controversial issue. The use of computerized tomography (CT) and magnetic resonance imaging (MRI) has greatly improved the accuracy of tumor localization in three-dimensional planning. This review aims at critically analyzing available literature data in which tumor extent of high-grade glioma has been assessed using CT and/or MRI and relating this to postmortem observations. Attention is given to the pattern of tumor spread at initial presentation and to tumor recurrence pattern after external beam irradiation. Special emphasis is given to the site of tumor regrowth after radiation treatment in relation to the boundaries of the CTV. Guidelines for delineating CTV will be inferred from this information, taking data on radiation effects on the normal brain into account.

Brain Neoplasms↗

The width of margins in radiotherapy treatment plans.

Publication of ICRU Reports 50 and 62 has highlighted the need to devise protocols for the process of drawing the planning target volume (PTV) around the clinical target volume (CTV). The margin surrounding the CTV should be wide enough to account for all geometric errors so that no part of the CTV accumulates a dose less than, for instance, 95% of that prescribed. One approach to the problem has been to draw a margin around the CTV delineated at the treatment preparation stage which is sufficiently wide that the mean position of the CTV will be encompassed in a specific percentage of cases, for example 90%. This accounts for the systematic errors. A further margin is then drawn to account for random set-up and organ-motion uncertainties during treatment. The width of this second margin has previously been shown to be 1.64(sigma - sigmap). Here sigma, a vector quantity, is the standard deviation which results from convolving the penumbra spread function of standard deviation sigmap with the Gaussian distributions of the daily positional uncertainties of organ motion and set-up error. However, it is shown in this paper that the calculation should take into account the beam configuration of the treatment plan. In a typical coplanar multibeam plan, usually in the transverse plane, any given edge of the target volume is normally defined by a single beam or two parallel and opposed beams. However, because of the presence of the other beams, the effect of the blurring of the edge-defining beam(s) is reduced, which changes the value of the required margin to beta (sigma - sigmap) where, for example, beta can be as low as 1.04 in the transverse plane of a three-beam plan. The width of the required margins is calculated for up to six beams and presented in a table. It is shown that, while the table was derived using an idealized plan of equally weighted plane beams irradiating a spherical target, it is also valid for non-uniform beam weightings, wedged-beam plans, target volumes of general shape and intensity-modulated radiotherapy (IMRT).

Humans↗

The accuracy of image registration for the brain and the nasopharynx using external anatomical landmarks.

We investigated the accuracy of 3D image registration using markers that are repeatedly applied to external anatomical landmarks on the head. The purpose of this study is to establish a lower limit of the errors that would occur in, for instance, MRI-SPECT matching, which in some situations can only be achieved using external landmarks. Marker matching was compared with (single-modality) volume matching for 20 MRI scans. The results were compared with a published expression for the target registration error (TRE) which gives the 3D distribution of the mismatch between both scans. It was found that the main error source is reapplying the external markers on the anatomical landmarks. The published expression describes the relative distribution of the TRE in space well, but tends to underestimate the actual registration error. This deviation is due to anisotropy in the error distribution of the marker position (errors in the direction perpendicular to the skin surface are in general much smaller than errors in other directions). A simulation of marker matching with anisotropy in the errors confirmed this finding. With four reapplied markers, the TRE is 6 mm or smaller in most regions of the head.

Brain↗

Leaf position verification during dynamic beam delivery: a comparison of three applications using electronic portal imaging.

The use of a dynamic multileaf collimator (MLC) to deliver intensity-modulated beams presents a problem for conventional verification techniques. The use of electronic portal imaging to track MLC leaves during beam delivery has been shown to provide a solution to this problem. An experimental comparison of three different verification systems, each using a different electronic portal imaging technology, is presented. Two of the systems presented are commercially available imagers with in-house modifications, with the third system being an in-house built experimental system. The random and systematic errors present in each of the verifications systems are measured and presented, together with the study of the effects of varying dose rate and leaf speed on verification system performance. The performance of the three systems is demonstrated to be very similar, with an overall accuracy in comparing measured and prescribed collimator trajectories of approximately +/-1.0 mm. Systematic errors in the percentage delivered dose signal provided by the accelerator are significant and must be corrected for good performance of the current systems. It is demonstrated that, with suitable modifications, commercially available portal imaging systems can be used to verify dynamic MLC beam delivery.

Algorithms↗

Comparison of prostate cancer treatment in two institutions: a quality control study.

PURPOSE: To minimize differences in the treatment planning procedure between two institutions within the context of a radiotherapy prostate cancer trial. PATIENTS AND METHODS: Twenty-two patients with N0 M0 prostate cancer underwent a computed tomography (CT) scan for radiotherapy treatment planning. For all patients, the tumor and organs at risk were delineated, and a treatment plan was generated for a three-field technique giving a dose of 78 Gy to the target volume. Ten of the 22 cases were delineated and planned in the other institution as well. The delineated volumes and dose distributions were compared. RESULTS: All treatments fulfilled the trial criteria. The mean volume ratio of the gross tumor volumes (GTVs) in both institutions was 1.01, while the mean volume ratio of the planning target volumes (PTVs) was 0.88. The three-dimensional (3D) PTV difference was 3 mm at the prostate apex and 6-8 mm at the seminal vesicles. This PTV difference was mainly caused by a difference in the method of 3D expansion, and disappeared when applying an improved algorithm in one institution. The treated volume (dose > or =95% of isocenter dose) reflects the size of the PTV and the conformity of the treatment technique. This volume was on average 66 cm3 smaller in institution A than in institution B; the effect of the PTV difference was 31 cm3 and the difference in technique accounted for 36 cm3. The mean delineated rectal volume including filling was 112 cm3 and 125 cm3 for institution A and B, respectively. This difference had a significant impact on the relative dose volume histogram (DVH) of the rectum. CONCLUSION: Differences in GTV delineation were small and comparable to earlier quantified differences between observers in one institution. Different expansion methods for generation of the PTV significantly influenced the amount of irradiated tissue. Strict definitions of target and normal structures are mandatory for reliable trial results.

Humans↗

Computerized design of target margins for treatment uncertainties in conformal radiotherapy.

PURPOSE: We describe a computerized method of determining target margins for beam aperture design in conformal radiotherapy plans. MATERIALS AND METHODS: The method uses previously measured data from a population of patients to simulate setup error and organ motion in the patient currently being planned. Starting with a clinical target volume (CTV) and nontarget organs from the patient's planning CT scan, the simulation is repeated many times to produce a spatial probability distribution for each organ in the treatment machine coordinate system. This is used to determine a prescribed dose volume (PDV), defined as the volume to receive the prescribed dose, which encompasses the CTV while restricting the volume of nontarget organs within it, according to planner-specified values. The PDV is used to design beam apertures using a conventional margin for beam penumbra. RESULTS: The method is applied to 6-field prostate conformal treatment plans, in which the PDV encloses the prostate and seminal vesicles while limiting the enclosed rectal wall volume. The effect of organ motion is assessed by applying the plans on subsequent CT scans of the same patients, calculating probabilities for tumor control (TCP) and normal tissue complication (NTCP), and comparing with plans designed from a physician-drawn planning target volume (PTV). Although prostate TCP and rectal wall NTCP are found to be similar in the two sets of plans, TCP for the seminal vesicles is significantly higher in the PDV-based plans. CONCLUSIONS: The method can improve the dose conformality of treatment plans by incorporating population-based measurements of treatment uncertainties and consideration of nontarget tissues in the design of nonuniform target margins.

Computer Simulation↗

Definition of the prostate in CT and MRI: a multi-observer study.

PURPOSE: To determine, in three-dimensions, the difference between prostate delineation in magnetic resonance (MR) and computer tomography (CT) images for radiotherapy treatment planning. PATIENTS AND METHODS: Three radiation oncologists, considered experts in the field, outlined the prostate without seminal vesicles both on CT, and axial, coronal, and sagittal MR images for 18 patients. To compare the resulting delineated prostates, the CT and MR scans were matched in three-dimensions using chamfer matching on bony structures. The volumes were measured and the interscan and interobserver variation was determined. The spatial difference between delineation in CT and MR (interscan variation) as well as the interobserver variation were quantified and mapped three-dimensionally (3D) using polar coordinates. A urethrogram was performed and the location of the tip of the dye column was compared with the apex delineated in CT and MR images. RESULTS: Interscan variation: CT volumes were larger than the axial MR volumes in 52 of 54 delineations. The average ratio between the CT and MR volumes was 1.4 (standard error of mean, SE: 0.04) which was significantly different from 1 (p < 0.005). Only small differences were observed between the volumes outlined in the various MR scans, although the coronal MR volumes were smallest. The CT derived prostate was 8 mm (standard deviation, SD: 6 mm) larger at the base of the seminal vesicles and 6 mm (SD 4 mm) larger at the apex of the prostate than the axial MRI. Similar figures were obtained for the CT and the other MRI scans. Interobserver variation: The average ratio between the volume derived by one observer for a particular scan and patient and the average volume was 0.95, 0.97, and 1.08 (SE 0.01) for the three observers, respectively. The 3D pattern of the overall observer variation (1 SD) for CT and axial MRI was similar and equal to 3.5 to 2.8 mm at the base of the seminal vesicles and 3 mm at the apex. CONCLUSION: CT-derived prostate volumes are larger than MR derived volumes, especially toward the seminal vesicles and the apex of the prostate. This interscan variation was found to be larger than the interobserver variation. Using MRI for delineation of the prostate reduces the amount of irradiated rectal wall, and could reduce rectal and urological complications.

Aged↗

Quantification and predictors of prostate position variability in 50 patients evaluated with multiple CT scans during conformal radiotherapy.

PURPOSE: To determine the extent and predictors for prostatic motion in a large number of patients evaluated with multiple CT scans during radiotherapy, and evaluate the implications of these data on the design of appropriate treatment margins for patients receiving high-dose three-dimensional conformal radiotherapy. MATERIALS AND METHODS: Fifty patients underwent four serial computerized tomography (CT) scans, consisting of an initial planning scan and subsequent scans at the beginning, middle, and end of the treatment course. Each scan was performed with the patient in the prone treatment position within an immobilization device used during therapy. Contours of the prostate and seminal vesicles were drawn on the axial CT slices of each scan, and the scans were matched by alignment of the pelvic bones with a chamfer matching algorithm. Using the contour information, distributions of the displacement of the organ center of mass and organ border from the planning position were determined separately for the prostate and seminal vesicles in each of the three principle directions: anterior-posterior (AP), superior-inferior (SI) and left-right (LR). Each distribution was fitted to a normal (Gaussian) distribution to determine confidence limits in the center of mass and border displacements and thereby evaluate for the optimal margins needed to contain target motion. RESULTS: The most common directions of displacement of the prostate center of mass (COM) were in the AP and SI directions and were significantly larger than any LR movement. The mean prostate COM displacement (+/- 1 standard deviation, SD) for the entire population was -1.2 +/- 2.9 mm, -0.5 +/- 3.3 mm and -0.6 +/- 0.8 mm in the, AP and SI and LR directions respectively (negative values indicate posterior, inferior or left displacement). The mean (+/- 1 SD) seminal vesicle COM displacement for the entire population was - 1.4 +/- 4.9 mm, 1.3 +/- 5.5 mm and -0.8 +/- 3.1 mm in the AP and SI and LR directions, respectively. The data indicate a tendency for the population towards posterior displacements of the prostate from the planning position and both posterior and superior displacements of the seminal vesicles. AP movement of both the prostate and seminal vesicles were correlated with changes in rectal volume (P = 0.0014 and < 0.0001, respectively) more than with changes in bladder volume (P = 0.030 for seminal vesicles and 0.19 for prostate). A logistic regression analysis identified the combination of rectal volume > 60 cm3 and bladder volumes > 40 cm3 as the only predictor of large ( > 3 mm) systematic deviations for the prostate and seminal vesicles (P = 0.05) defined for each patient as the difference between organ position in the planning scan and mean position as calculated from the three subsequent scans. CONCLUSIONS: Prostatic displacement during a course of radiotherapy is more pronounced among patients with initial planning scans with large rectal and bladder volumes. Such patients may require more generous margins around the CTV to assure its enclosure within the prescription dose region. Identification and correction of patients with large systematic errors will minimize the extent of the margin required and decrease the volume of normal tissue exposed to higher radiation doses.

Artifacts↗

A general methodology for three-dimensional analysis of variation in target volume delineation.

A generic method for three-dimensional (3-D) evaluation of target volume delineation in multiple imaging modalities is presented. The evaluation includes geometrical and statistical methods to estimate observer differences and variability in defining the Gross Tumor Volume (GTV) in relation to the diagnostic CT and MRI modalities. The geometrical method is based on mapping the 3-D shape of the target volume to a scalar representation, thus enabling a one-dimensional statistical analysis. The statistical method distinguishes observer and modality related uncertainties, which are expressed in terms of three error components: random observer deviations, systematic observer differences, and systematic modality differences. Monte Carlo simulations demonstrate that the standard errors of each of the three model parameters are inversely proportional to the square root of the product of the patient group size and the number of observers and proportional to the intraobserver variation. For 18 patients and 3 observers the standard errors of the estimated systematic modality and observer differences are 19% and 14% of the intraobserver standard deviation, respectively. A scalar representation of the shape of the prostate, delineated by 3 observers for 18 patients, was obtained by sampling the distance between the average center of gravity of the prostate in CT and the prostate surface for a large number of directions (2500), using polar coordinates. Observer variability and differences were obtained by applying the statistical method to the samples independently. The intraobserver variation for CT was largest in regions near the seminal vesicles (s.d: 3 mm) and the apex (s.d: 3 mm). The systematic observer variation in CT was largest in a region near the plexus Santorini, at the caudal-anterior side of the prostate (s.d.: 2 mm). The sensitivity for the choice of origin was tested by using the average center of gravity from axial MRI instead of CT. The results were almost identical. The polar map measures distances in the scanning directions. A correction procedure to get the variability in directions perpendicular to the surface of the prostate yielded variations that were a factor of 0.85 smaller for all directions. It is concluded that by separating the shape evaluation in a geometrical and a statistical part, the complexity of the analysis of 3-D shape differences can be significantly reduced. The method was successfully applied to a group of prostate patients, where we demonstrated that delineation variability is nonhomogeneous, with the largest variations occurring near the seminal vesicles and the apex.

Biophysical Phenomena↗