BJR
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Avilés Lucas, P
Right arrow Articles by Vañó Carruana, E
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Avilés Lucas, P
Right arrow Articles by Vañó Carruana, E
British Journal of Radiology 74 (2001),1128-1136 © 2001 The British Institute of Radiology

Full paper

Analysis of surface dose variation in CT procedures

P Avilés Lucas, MSc 1 I A Castellano, MSc, MIPEM 1 D R Dance, PhD, FIPEM 1 and E Vañó Carruana, PhD 2

1Physics Department, The Royal Marsden NHS Trust, Fulham Road, London SW3 6JJ, UK and 2Departamento de Radiología, Facultad de Medicina, Universidad Complutense de Madrid, Madrid 28040, Spain


    Abstract
 Top
 Abstract
 Introduction
 Method
 Results and discussion
 Conclusions
 References
 
An experimental and theoretical study has been made of the variations in air kerma–length product (AKLP) at the surface of a phantom exposed in a CT scanner, using clinical parameters. For the theoretical part, a computer simulation was developed, based on a simple analytical model, requiring information generic to the scanner model. The effects of patient size, position within the gantry plane and beam-shaping filter type were studied using three different elliptical phantoms. A dose reduction technique based on a sinusoidal tube current modulation system was also investigated. The surface AKLP was shown to be independent of phantom size (within experimental error) but decreases as the surface moves vertically away from the isocentre. The major contributor to this variation is the beam-shaping filter. A maximum difference of 19% between the values of surface AKLP was calculated for the two beam-shaping filters available. When the tube current modulation system was used, the maximum surface dose reduction was 18%. A maximum difference of 2% was found between measurements and computer simulations. It is therefore possible to predict the behaviour of AKLP using the analytical model. As the dose measured on the surface of a patient is simply related to AKLP, the model can be applied to data obtained from patient surface dose surveys and can be helpful in interpreting the sources of variation in the latter.


    Introduction
 Top
 Abstract
 Introduction
 Method
 Results and discussion
 Conclusions
 References
 
CT is a high dose diagnostic technique, estimated to be responsible for 40% of the total collective dose from all diagnostic procedures in the UK [13]. State-of-the-art CT scanners (multislice CT scanners) now allow a wider choice of slice thickness and a decrease in exposure time. Such scanners have the potential to improve visualization of small structures and so examination frequency is expected to further increase [4]. Consequently, there is a need to evaluate patient doses and the risk associated with these procedures. It is very important that CT examination exposure parameters are optimized and patient dose is kept as low as reasonably achievable consistent with the need to obtain adequate diagnostic information. This is in accordance with a recent European Directive [5]. Therefore, measurement of patient dose forms an essential part of the CT quality assurance programme.

Patient dose can be estimated from the CT dose index (CTDI) [6], which for this purpose is defined as Go


where D(z) is the dose from a single slice expressed as a function of position along an axis parallel to the axis of rotation of the scanner (z-axis) and T is the nominal slice thickness. The integral is evaluated over 10 cm, which is the active length of the pencil ionization chambers used routinely for this measurement.

In practice, CTDI is evaluated free-in-air or in a standard American Association of Physicists in Medicine (AAPM) head or body phantom (cylindrical phantoms of 16 cm or 32 cm diameter, and 15 cm thickness made from polymethylmethacrylate (PMMA)) [7]. In the case of measurements in AAPM phantoms, CTDI can be averaged over the phantom radius [8]. The resulting weighted CTDI (CTDIw) may then be combined, as follows, with the known scan parameters for a complete CT examination to produce a dose–length product (DLP), which is more representative of the total risk than the CTDI itself: Go


where N is the number of slices and C is the product of tube current and exposure time used for each slice.

Alternatively, CTDI may be combined with conversion factors (calculated using Monte Carlo techniques) that relate it to organ doses or effective doses for single slices. The conversion factors currently available [9, 10] are scanner-specific and require the measurement of CTDI free-in-air on the axis of the scanner, a quantity that does not account for the quality and geometry of the X-ray beam (including the effect of the bow tie filter) in a given scanner. The EU CT protocol [8] proposes a simplified approach that applies examination-specific conversion factors, for any given CT scanner, to the DLP values.

In both the above approaches, the measured dosimetric quantity is not dependent upon patient size or their position within the gantry plane. However, CT procedures and technical parameters are (or should be) tailored to individual patient sizes and clinical indications [11]. It has been recognized [1214] that patient size needs to be taken into account when designing imaging protocols for patients undergoing CT examinations. Thus, it can be considered desirable that dose measurements for a particular examination and scanner reflect the variability that should be found in practice.

A complementary approach to the in-air or in-phantom measurements is the direct measurement of patient surface dose using thermoluminescent dosemeters (TLDs). We are developing such an approach [15, 16], which has the advantage of taking into account individual customized CT parameters and patient size and position. This method also has the potential to be used in routine patient surveys since it is inexpensive (the dosemeters could simply be sent to the clinic with a protocol) and may be more convenient for use in centres with a high workload as it does not require measurement of CTDI. As part of the development of this technique, it was considered important to study the variation of surface dose with patient size and gantry position. This paper reports a study of their effects for single CT slices, based on both phantom measurements and computer simulation. The influence of tube current modulation techniques, which are available on some CT scanners for saving dose, is also investigated.


    Method
 Top
 Abstract
 Introduction
 Method
 Results and discussion
 Conclusions
 References
 
The variation of surface dose with patient size and position within the gantry plane, the choice of bow tie filter and the application of tube current modulation techniques have been studied with measurements and computer simulations of exposures in a GE HiSpeed CT/i scanner (GE Medical Systems, Milwaukee, WI). This is a third generation, helical CT scanner with a single-row, solid state ceramic detector array.

The scanner has two beam-shaping filters available for use in different clinical scanning protocols. The standard filter is referred to as the "small" bow tie filter whilst the other is referred to as the "large" or SmartBeamTM, bow tie filter. These filters cause a reduction in radiation dose at the phantom or patient edges to compensate for the shorter path lengths of the fan beam through the periphery of the patient. The small bow tie filter is designed for general scanning but is most suited to small cross-sections through the patient (such as head examinations). The large bow tie filter is suited to scanning larger cross-sections through the patient and is selected for most examinations in our hospital.

The GE HiSpeed CT/i scanner also provides an optional tube current modulation system called SmartScanTM, which allows dose reduction in projections where attenuation of the X-ray beam decreases. In this mode, tube current varies sinusoidally throughout the entire rotation of the tube, decreasing tube current (mA) for anteroposterior (AP) projections (assuming these to be the least attenuating regions) and increasing it for lateral projections (assuming these to be the most attenuating regions), up to a maximum value selected by the operator. The CT scanner calculates tube current reduction using the ratio of patient attenuations as seen by the central part of the detector array from AP and lateral views. The maximum reduction possible is 50%.

In this work, measurements and calculations were made with three different elliptical phantoms. Their dimensions are given in Table 1Go and were chosen from our database, which contains measurements of 82 patients examined using standard CT procedures. Phantom A is more representative of the human dimensions encountered in the shoulders or upper thorax examinations, whereas phantoms B and C reproduce the mean size of the lower thorax and the abdomen. The phantoms used for the measurements were specially constructed from PMMA and, for computational purposes, were modelled as such.


View this table:
[in this window]
[in a new window]
 
Table 1. Cross-sectional dimensions of the three elliptical phantoms used for measurements and calculations. The phantoms had a 10 cm depth

 
As CTDI measurements are always based on the line-integral of a dose-related quantity along the length of a pencil chamber or set of TLDs, all the results and analysis presented are expressed in terms of air kerma–length product (AKLP), defined as follows: Go


where Ds(z) is the surface dose from a single slice expressed as a function of position along the z-axis. The integral is over the active length of the chamber. AKLP is simply the surface CTDI multiplied by the nominal slice thickness.

Measurements of the variations in surface AKLP
Every experiment was initially set up by positioning the selected phantom symmetrically in the centre of the gantry with a 10 cm long pencil ionization chamber (model No 20X5-CT with a MDH 2025 electrometer; Radcal Corporation, Monrovia, CA) on its anterior surface, positioned so that its axis was parallel to the rotation axis of the scanner. The chamber was aligned with the help of lasers. The anterior surface was chosen, as this is the part of the patient most easily accessible for dose surveys. The position of the phantom and chamber within the gantry could be verified from the CT image obtained. Therefore to evaluate the accuracy in positioning, we could make use of the grid that is overlapped onto the phantom image.

The height of the selected phantom could be adjusted by raising or lowering the couch along the vertical axis of the tomographic plane. The error in positioning of phantoms at the isocentre of the CT scanner was ±3 mm. The table has a load capacity of 180 kg, with an error of ±0.25 mm in the position along the vertical axis of the tomographic plane. For each measurement of AKLP, a single exposure consisting of three 360° rotations with no table motion was performed. This strategy was chosen to reduce any error arising from changes between exposures in the starting point of the tube and, therefore, the discontinuity in X-ray emission that occurs at this location. Results were expressed for a single 360° rotation and the coefficient of variation of AKLP from a single rotation was estimated to be ±0.4%. Positional and AKLP uncertainties were combined and the total error in AKLP was estimated to be 1%. Exposure settings were chosen to be consistent with clinical scan protocols, but for practical purposes were kept constant for all phantoms regardless of their size. Tube current was 200 mA, tube voltage 120 kVp, slice thickness 10 mm and rotation speed 360° s-1.

To study the influence of phantom size and phantom position within the tomographic plane on AKLP, a series of exposures was carried out with the elliptical phantoms A and B. The large bow tie filter was used for all measurements. AKLP was measured at positions along the vertical axis of the tomographic plane (y-axis). This experimental configuration is illustrated in Figure 1Go.



View larger version (17K):
[in this window]
[in a new window]
 
Figure 1. Configuration used to measure or calculate the air kerma–length product on the surface of an elliptical phantom (for clarity, the figure is not drawn to scale). The chamber is placed on the upper phantom surface. The horizontal and vertical axes of the tomographic plane are given by x and y, whereas z refers to the rotation axis. R is the distance between the focus and the isocentre, {theta} is the projection angle, {psi} is the angle between the chamber–focus and the focus–isocentre vectors, r is the distance between the focus and the centre of the chamber and h is distance of the ionisation chamber from the isocentre.

 
The same configuration was used to study the effect of the bow tie filter on AKLP. The AKLP integral was measured on the top of phantom A, with and without the bow tie filter, at a series of positions along the y-axis.

To study the effect of dose reduction caused by the tube current modulation system, AKLPs were also measured on the surface of the three elliptical phantoms with and without the tube current modulation feature. The maximum or prescribed mA and the average mA (Smart mATM) were recorded for measurements with the tube current modulation technique. For this experiment, the phantoms were set up symmetrically in the centre of the tomographic plane since the system will then carry out the reduction on the basis of difference in attenuation between the minor and major axes of the ellipses.

Calculation of the variations in surface AKLP
To understand and further investigate sources of variation in surface AKLP along the vertical axis of the tomographic plane, a computer simulation based on a straightforward model has been developed. In this model the AKLP measured with the CT chamber is calculated as an integral along a line on the surface of the patient (z-direction) parallel to the axis of rotation of the scanner. Let Ko(z,{theta},h) X({psi})d{theta} be the kerma at position z along this line, in the absence of the phantom and beam shaping filter, for a tube rotation through angle d{theta} from CT projection angle {theta} and at a given chamber position along the y-axis h. The function X({psi}) takes account of variation of the attenuation through the flat filter with angle of the X-ray beam ({psi} is shown in Figure 1Go). The flat filter consists of inherent and additional filtration. When the radiation is normal to the filter, X({psi}) is unity. The AKLP, P, evaluated along this line with the phantom and beam-shaping filter present is given by Go


where L is the length of the ionization chamber, which is positioned symmetrically about z=0 and at distance h from the isocentre (Figure 1Go). Attenuation through the beam-shaping filter is modelled by the weighting factor w({psi}). This factor is defined as the ratio of air kerma free-in-air at a given value of {psi} with the bow tie filter present to that with the bow tie filter removed. The exponential term in Equation (4)Go accounts for the attenuation through the phantom. Beam hardening is neglected and a constant linear attenuation coefficient µ is assumed. The length q({theta},h) is the distance the ray traverses through the phantom (z-dependence is neglected). Attenuation through the couch is not modelled. The term S(z,{theta},h) allows for the contribution of scattered radiation to the AKLP. In practice this function was not known and it was necessary to assume that it could be replaced by some average value S for those rays making a significant contribution to the integral. The factor C({theta}) introduces sinusoidal modulation of the tube current.

The kerma function Ko(z,{theta},h) can be rewritten to make its inverse-square law dependence explicit: Go


For slice widths of 10 mm or less the z2 term can be neglected. The length r is the distance between the focus and the centre of the chamber (Figure 1Go). The function F(z,{theta},h) can be regarded as a beam profile along the z-direction, and for the free-in-air situation will scale according to the value of rcos{psi}. Therefore, Go


where k is given by: Go


The length R is the distance between the focus and the isocentre. The assumption of free-in-air scaling of profile shape is valid provided a rectangular collimator is used to limit the extent of the beam in the z-direction. Figure 2Go shows schematically how the beam profile shape scales with increasing distance from the focus. The scaling factor k is the ratio of the two distances indicated in the figure.



View larger version (12K):
[in this window]
[in a new window]
 
Figure 2. Two schematic free-in-air kerma profiles parallel to the axis of rotation of the CT scanner (z-axis), shown at distance rcos{psi} from the focus, at the isocentre. When rectangular collimation is used, the shapes of the two profiles are simply related using similar triangles once the inverse-square law dependence is explicitly stated in accordance with Equation (5)Go. The distances rcos{psi} and R are shown in Figure 1Go.

 
Equations (4)Go to (6)Go can be combined using the above approximations to give: Go


Substituting Z=kz gives: Go


Provided the limits of integration are sufficiently far from the edge of the in-air beam profile (which is always the case), P becomes Go


where Equation (7)Go has been used to substitute for k and the constant f is given by Go


It will be seen that in this approximation, AKLP has an inverse distance dependence. In addition, it is not necessary to know S or f to compare calculations and experiments since only relative values are needed.

Finally, the factor C({theta}) in Equation (4)Go is given by the following expression [17]: Go


where mAmin is the minimum tube current selected by the system and mAmax is the maximum tube current prescribed by the user of the CT scanner. In the case of no tube current modulation, minimum and maximum tube current are equal and C({theta}) is then unity. The average tube current (mAmean) is: Go


Since the value of the mean tube current is provided by the system, mAmin can be calculated using this equation.

Equation (10)Go has been evaluated numerically for various configurations. The product of X({psi}) and w({psi}) was determined experimentally for the two bow tie filters available in the HiSpeed CT/i scanner. For this experiment the tube was stationary and AKLP free-in-air was recorded with the large and small bow tie filter for a series of positions, regularly spaced in intervals of 1 cm, along an axis normal to the central beam and containing the isocentre. The experimental set-up was based on that used by Jansen et al for the same purpose [18]. For the calculation of the product of X({psi}) and w({psi}), account was taken of the variation of the focus-to-chamber distance.

Equation (4)Go was formulated on the assumption that spectral variations can be accounted for by the choice of a suitable value of the linear attenuation coefficient µ. The value used for the majority of calculations (0.21 cm-1 for PMMA) was determined on the basis of previous work [17, 19] and was chosen to represent a typical 120 kVp spectrum. However, to investigate the sensitivity of the results to the value of the attenuation coefficient, some calculations were repeated with values of 0.154 cm-1 and 0.314 cm-1, as an error in the choice of the attenuation coefficient of 30% or greater was deemed to be very unlikely.

To validate the model, the variation of AKLP with position was calculated for phantoms A and B and the results were compared with the corresponding measurements. The validated model was then used to further investigate sources of variation in the surface AKLP, including the presence or absence of the bow tie filter, the dependence on size (for all phantoms) and the effect of the GE HiSpeed CT/i dose reduction technique.


    Results and discussion
 Top
 Abstract
 Introduction
 Method
 Results and discussion
 Conclusions
 References
 
Measurements of the variations in surface AKLP
Figure 3Go shows the variation of AKLP with offset from the isocentre of the centre of phantoms A and B. It might be expected that differences between the two sets of results could be attributed to the effect of phantom size. However, a different conclusion can be drawn from Figure 4Go, where AKLP is plotted against chamber position along the y-axis of the scanner. Thus, Figure 4Go shows that AKLP for both phantoms varies in the same way. Indeed, results for the two phantoms fall on the same curve. For three pairs of results, phantoms were positioned so that the chamber ended up at almost the same point on the y-axis (within 0.15 cm). In these cases the difference between AKLP values is within 2%. Consequently, AKLP depends upon chamber position and not phantom size. These data also indicate that AKLP decreases with increasing distance from the isocentre.



View larger version (11K):
[in this window]
[in a new window]
 
Figure 3. Surface air kerma–length product (AKLP) as a function of the offset of the phantom centre along the vertical y-axis of the tomographic plane. {blacktriangleup}, measurements with phantom A; •, measurements with phantom B.

 


View larger version (10K):
[in this window]
[in a new window]
 
Figure 4. Normalized surface air kerma–length product (AKLP) as a function of the chamber position along the y-axis. Results are normalized to the maximum surface AKLP measured for phantom A at the point closest to the isocentre (the latter being the origin). {blacktriangleup}, measurements with phantom A; •, measurements with phantom B.

 
The source of variation with position shown in Figure 4Go can be understood by reference to Figure 5Go, which shows the results of measurements with and without the large bow tie filter. It can be seen that AKLP on the surface of phantom A obtained without the bow tie filter is, within experimental error, independent of chamber position along the y-axis. Consequently the major contributor to AKLP variation along the y-axis of the CT scanner is the beam-shaping filter. For results with the bow tie filter, maximum surface AKLP is obtained at the isocentre. In this position, attenuation of the beam-shaping filter is least and is constant throughout the entire rotation of the tube. The symmetry in shape of the bow tie filter will cause surface AKLP to decrease symmetrically below and above this point.



View larger version (10K):
[in this window]
[in a new window]
 
Figure 5. Normalized air kerma–length product (AKLP) on the surface of phantom A as a function of the chamber position along the y-axis. Results are normalized to the value measured at the isocentre (origin of the coordinate system). {blacktriangleup}, measurements with the large bow tie filter; •, measurements without the bow tie filter.

 
AKLP is shown in Figure 6Go for the three elliptical phantoms exposed with and without the tube current modulation technique. Results are normalized to the maximum AKLP value obtained with phantom A when the tube current modulation technique is not selected. As expected, surface AKLPs obtained with and without the tube current technique agree for phantoms with nearly circular cross-sections, since in that case there is very little variation between the attenuation in lateral and anterior projections. Differences between normalized values obtained with and without the current modulation system increase with the increase in eccentricity and decrease in thickness of the phantom, reaching a value of 18% for phantom A.



View larger version (13K):
[in this window]
[in a new window]
 
Figure 6. Normalized air kerma–length product (AKLP) on the surface of elliptical phantoms A, B and C. {blacktriangleup},measurements with the tube current modulation technique; •, measurements without the tube current modulation technique.

 
Calculation of the variations in surface AKLP
A validation of the simple theoretical model is given in Table 2Go, which shows the variation along the y-axis of both calculated and experimental AKLPs on the surface of phantoms A and B. The maximum discrepancy between measurements and simulations is 2%. The error associated with experimental measurements is 1%. It is more difficult to estimate the systematic error associated with the theoretical model, especially the assumptions about the scatter function S({theta},h,z). The limitations arising from the use of a single value of the linear attenuation coefficient µ could, however, be analysed by performing simulations with several µ values. It was found that a change in µ of at least 30% caused a maximum variation in the absolute value of AKLP of 5.3%, thus demonstrating that AKLP is insensitive to the value of µ. Another possible source of discrepancy is the absence of the couch in the model. Caon et al [20] considered the effect of the couch on radiation dose in phantoms by simulations and measurements using a fixed geometry. They found that when attenuation of the table was not included in their simulations, the dose at the phantom's centre was about 5% greater, whereas doses at the surface were almost unchanged. A 2–6% difference in the effective dose equivalent was found in a study by Jones and Shrimpton [9] and so the couch should be considered when organ and effective doses are calculated.


View this table:
[in this window]
[in a new window]
 
Table 2. Comparison between the experimental and calculated air kerma–length product (AKLP) on the surface of phantoms A and B at positions along the y-axis. Results are normalized to the maximum value obtained when the chamber is closest to the isocentre

 
In Figure 7Go, AKLP calculated for phantom B, in several positions along the y-axis, is shown for the two bow tie filters. Results are normalized to the maximum value obtained for each filter. The figure shows that normalized AKLP values decrease faster when the small bow tie filter is used. A maximum difference of 19% was obtained between the calculated AKLPs for the large and the small bow tie filters. This occurs at 12 cm from the isocentre. Differences between values for the two filters start to decrease when the chamber is far away from the isocentre (at 20 cm or more from the isocentre) where the filters offer similar attenuation. This exercise confirms the experimental finding that the behaviour of AKLP with distance from the isocentre is determined by the beam-shaping filter, and furthermore that variation is specific to filter shape.



View larger version (12K):
[in this window]
[in a new window]
 
Figure 7. Normalized air kerma–length product (AKLP) calculated on the surface of phantom B as a function of the chamber position along the y-axis. For each filter, results are normalized to the maximum AKLP. {blacktriangleup}, calculations with large bow tie filter; {square}, calculations with the small bow tie filter.

 
Table 3Go shows the normalized AKLPs calculated for the three elliptical phantoms. Data are normalized to the maximum value found for phantom C at the closest position to the isocentre. Differences between AKLPs calculated for the three phantoms for a given chamber position depend in principle upon attenuation through the phantom in each projection, which is related to its geometrical shape, the chamber position along the y-axis and the beam-shaping filter. The maximum difference calculated for a given chamber position between AKLP on phantoms A and C was found to be below 2%. This is in very good agreement with the experimental result found when comparing AKLP for phantoms A and B, which also showed a maximum difference of 2%. Therefore, for a given chamber position, phantom size has no effect on surface AKLP (within experimental errors of 1% and within the limitations of the model). However, because the patient dimension determines the position of the upper surface region of measurement in patient examinations (the patient is approximately centred within the gantry plane), patient size should be considered when surface air kerma is evaluated. For instance, in the case of two elliptical phantoms of equivalent horizontal thickness and centred at the isocentre of the CT scanner, a 33% change in vertical thickness (from 16 cm to 24 cm) leads to a 27% change in surface AKLP.


View this table:
[in this window]
[in a new window]
 
Table 3. Calculated values of the air kerma–length product (AKLP) on the surface of phantoms A, B and C. The positions chosen match those selected in the experiment carried out with phantom A and B. Results are normalized to the maximum value, which was obtained for phantom C

 
Finally, AKLP calculated with and without tube current modulation for the three elliptical phantoms has been compared with experimental values. The agreement between simulated and experimental results was found to be within 2%. Consequently, the model developed in this work is also able to simulate AKLP variation along the y-axis when dose reduction techniques are used.

Application to patient dose surveys
The extension of the above results to patient measurements with TLDs is straightforward. In clinical procedures, the practical dosimetric quantity to be measured is the surface multiple scan average dose (MSAD) [6]. Surface MSAD is related to AKLP as follows: Go


where I is the interval between slices. As this is a simple linear relationship, the results worked out in this paper are equally applicable to surface MSAD measured on patients, as long as sufficient slices are taken [6]; this should not pose a problem if the measurement site is chosen carefully.

The experiments and calculations above indicate that surface MSAD will depend on the offset of the dosemeter from the isocentre (as determined by patient dimensions and positioning), as well as any tailoring of the technical parameters to the size of that patient. The simple model developed in this work can be readily used to derive a normalized function of surface MSAD with offset from the isocentre. The only data required are the focus-to-isocentre distance and the X-ray attenuation distal from the central beam axis caused by the flat and beam-shaping filters. These features are generic to a scanner model and do not need to be measured directly on every scanner. If the offset of the dosemeter from the isocentre is obtained at the time of the survey, e.g. from the clinical images, then variation in surface MSAD owing to offset of the dosemeter from the isocentre can be corrected for, thus leaving a distribution that reflects variations in scanning practices rather than patient positioning.


    Conclusions
 Top
 Abstract
 Introduction
 Method
 Results and discussion
 Conclusions
 References
 
In this study, measurements and simple model-based calculations of the variations in surface AKLPs were made for different elliptical phantoms, beam-shaping filters and a dose reduction technique. The maximum difference between calculations and measurements was 2%. The diversity of experiments simulated proves the versatility of our simple model. Limitations in the model were also investigated. A 30% change in the value of the linear attenuation coefficient caused a 5% variation in surface AKLP. Therefore, AKLP is determined mainly by the primary radiation arriving at the phantom surface and the contribution of scattered radiation.

For a given chamber position, surface AKLP is essentially independent of phantom size for the phantoms investigated. Surface AKLP decreases as the surface moves away from the isocentre along the vertical axis of the tomographic plane. The decrease is owing to the effect of the beam-shaping filter and, for the large filter for the GE CT/i HiSpeed scanner, amounts to 19% at 12 cm from the isocentre. When a dose reduction technique based on a sinusoidal modulation of tube current was applied, the maximum dose saving measured on the surface of the elliptical phantoms was 18%. AKLP is simply related to surface CTDI, hence conclusions about the behaviour of the former apply equally to the latter.

The simple model has not only enhanced our understanding of the factors that determine surface AKLP or CTDI, but we believe that its predictions can be simply extended to patient dose measurements. AKLP is linearly related to surface MSAD (measurable with TLDs), hence the above conclusions about the behaviour of AKLP will be equally valid for MSAD. It is envisaged that the model can be easily implemented for any CT scanner, requiring only information generic to the scanner model. TLD data from a patient surface dose survey can then be interpreted; variations arising from changes in practice, as opposed to variations in patient size and position, can be isolated.


    Acknowledgments
 
P Avilés Lucas has been funded by a EU Marie Curie grant (Proposal number ERB4001GT), which is gratefully acknowledged. We thank the Medical Physics Group (Complutense University of Madrid) and the Medical Physics Service at the Clínico Hospital of Madrid for discussions and Prof. José Manuel Calleja Pardo (Autonomous University of Madrid) for his support of the project. We are indebted to Mr Bruce Bembrick for construction of the elliptical phantoms.

Received for publication October 2, 2000. Revision received April 30, 2001. Accepted for publication June 14, 2001.


    References
 Top
 Abstract
 Introduction
 Method
 Results and discussion
 Conclusions
 References
 

  1. Shrimpton PC. Reference doses for computed tomography. Radiological Protection Bulletin 193. Chilton: NRPB, 1997.
  2. Shrimpton PC, Edyvean S. CT scanner dosimetry. Br J Radiol 1998;71:1–3.[Medline]
  3. Shrimpton PC, Wall BF. The increasing importance of X-ray computed tomography as a source of medical exposure. Radiat Prot Dosim 1995;57:413–5.[Abstract]
  4. Yamashita Y, Murashima S, Sakuma H, Matsumura K, Takeda K. Assessment of lung cancer using multidetector-row CT and multiplanar reformation. Eur Radiol 2000;10(Suppl. 1):162.[Medline]
  5. Council Directive of the European Communities 97/43/EURATOM. Directive of 30 June 1997 on health protection of individuals against the dangers of ionizing radiation in relation to medical exposures, and repealing Directive 84/466/EURATOM. Official Journal of the European Communities 1997;L180:22–7.
  6. Shope TB, Gagne RM, Johnson GC. A method for describing the doses delivered by transmission X-ray computed tomography. Med Phys 1981;8:488–95.[Medline]
  7. Shope TB, Morgan TS, Showalter CK, Pentlow KS, Rothenberg LN, White DR, et al. Radiation dosimetry survey of computed tomography systems from ten manufacturers. Br J Radiol 1982;55:60–9.[Abstract/Free Full Text]
  8. Menzel H-G, Schibilla H, Teunen DC, editors. European guidelines on quality criteria for computed tomography. EUR 16262 EN. Brussels: EU, 1998. http://www.drs.dk/guidelines/ct/quality
  9. Jones DG, Shrimpton PC. Normalised organ doses for X-ray computed tomography calculated using Monte Carlo techniques. NRPB-SR250. Chilton: NRPB, 1993.
  10. Zankl M, Panzer W, Drexler G. The calculation of dose from external photon exposures using reference human phantoms and Monte Carlo methods. Part VI: organ doses from computed tomographic examinations. GSF-Berich 30/91. Germany 1991.
  11. Wade JP, Weyman JC, Goldstone KE. CT standard protocols are of limited value in assessing actual patient dose. Br J Radiol 1997;70:1146–51.[Abstract]
  12. Huda W, Atherton JV, Ware DE, Cumming WA. An approach for the estimation of effective radiation dose at CT in paediatric patients. Radiology 1997;203:417–22.[Abstract/Free Full Text]
  13. Ware DE, Huda W, Mergo PJ, Litwiller AL. Radiation effective doses to patients undergoing abdominal CT examinations. Radiology 1999;210:645–50.[Abstract/Free Full Text]
  14. Huda W, Scalzetti EM, Roskopf M. Effective doses to patients undergoing thoracic computed tomography examinations. Med Phys 2000;27:838–44.[Medline]
  15. Avilés Lucas P, Dance DR, Castellano IA, Skinner CL, Vañó Carruana E. Dosis superficiales empleadas para estimar el riesgo radiológico en pacientes examinados en procedimientos de tomografia computarizada. XII Congreso Nacional de Física Médica, 1999. CD-ROM: Session A2.
  16. Avilés Lucas P, Dance DR, Castellano IA, Skinner CL, Vañó Carruana E. A method to audit individual patient doses in CT procedures using thermoluminescent dosimeters. Eur Radiol 2000;10 (Suppl. 1):151.
  17. Gies M, Kalender WA, Wolf H, Suess C. Dose reduction in CT by anatomically adapted tube current modulation. I. Simulation studies. Med Phys 1999;26:2235–47.[Medline]
  18. Jansen JThM, Geleijns J, Zweers D, Schultz FW, Zoetelief J. Calculation of computed tomography dose index to effective dose conversion factors based on measurements of the dose profile along the fan shaped beam. Br J Radiol 1996;69:33–41.[Abstract/Free Full Text]
  19. Harpen MD. A simple theorem relating noise and patient dose in computed tomography. Med Phys 1999;26:2231–4.[Medline]
  20. Caon M, Bibbo G, Pattison J. A comparison of radiation dose measured in CT dosimetry phantoms with calculations using EGS4 and voxel-based computational models. Phys Med Biol 1997;42:219–29.[Medline]



This article has been cited by other articles:


Home page
Am. J. Roentgenol.Home page
I. A. Tsalafoutas, V. Tsapaki, C. Triantopoulou, C. Pouli, V. Kouridou, I. Fagadaki, and J. Papailiou
Comparison of Measured and Calculated Skin Doses in CT-Guided Interventional Procedures
Am. J. Roentgenol., November 1, 2008; 191(5): 1601 - 1607.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Roentgenol.Home page
J. Li, U. K. Udayasankar, T. L. Toth, J. Seamans, W. C. Small, and M. K. Kalra
Automatic Patient Centering for MDCT: Effect on Radiation Dose
Am. J. Roentgenol., February 1, 2007; 188(2): 547 - 552.
[Abstract] [Full Text] [PDF]


Home page
Radiat Prot DosimetryHome page
P. Aviles Lucas, D. R. Dance, I. A. Castellano, and E. Vano
Estimation of the peak entrance surface air kerma for patients undergoing computed tomography-guided procedures
Radiat Prot Dosimetry, May 17, 2005; 114(1-3): 317 - 320.
[Abstract] [Full Text] [PDF]


Home page
RadiologyHome page
M. K. Kalra, C. Wittram, M. M. Maher, A. Sharma, G. B. Avinash, K. Karau, T. L. Toth, E. Halpern, S. Saini, and J.-A. Shepard
Can Noise Reduction Filters Improve Low-Radiation-Dose Chest CT Images? Pilot Study
Radiology, July 1, 2003; 228(1): 257 - 264.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Avilés Lucas, P
Right arrow Articles by Vañó Carruana, E
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Avilés Lucas, P
Right arrow Articles by Vañó Carruana, E


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
BJR DMFR IMAGING  ALL BIR JOURNALS