Search PubMed⌕ Search

Biomedical subjects

C S Chui

Publications and source records attributed to C S Chui.

51 records · Page 3Linked to original sources

Stereotactic treatment of brain tumors with radioactive implants or external photon beams: radiobiophysical aspects.

We perform calculations, based on the linear-quadratic model, to assess the biologically effective doses (BED) of tumor and normal tissue in the stereotactic irradiation of brain tumors with either radioactive implants or radiosurgery techniques. Treatment protocols for radiosurgery and radioactive implants, as obtained from the literature, are reviewed and compared. A figure of merit is defined to be the ratio of tumor to normal tissue BED, expressed in units of Gy10/Gy3. These comparisons indicate a clear radiobiological advantage for brachytherapy, unless the radiosurgery is to be delivered in a large number of fractions. The differences in dose uniformity, and in the volume of normal tissue encompassed by the high dose regions, are factors that may also influence clinical results.

Brachytherapy↗

Perspectives of multidimensional conformal radiation treatment.

We consider the present technological advancement that underlies the implementation of computer-controlled conformal radiotherapy. We also consider the developments in modern biology that may provide input to therapy planning. The concept of multidimensional conformal radiotherapy is advanced, which integrates geometrical precision and biological conformality, to optimize the treatment planning for individual patients, with a view to improve the overall success of radiotherapy.

Humans↗

Interinstitutional experience in verification of external photon dose calculations.

Under the auspices of NCI contracts, four institutions have collaborated to assess the accuracy of the pixel-based dose calculation methods they employ for external photon treatment planning. The approach relied on comparing calculations using each group's algorithm with measurements in phantoms of increasing complexity. The first set of measurements consisted of ionization chamber measurements in water phantoms in normally incident square fields, an elongated field, a wedged field, a blocked field, and an obliquely incident beam. The second group of measurements was carried out using thermoluminescent dosimeters in phantoms designed to investigate the effects of surface curvature, high density heterogeneities, and low density heterogeneities. The final study tested the entire treatment planning system, including CT data conversion, in an anthropomorphic phantom. Overall, good agreement between calculation and measurements was found for all algorithms. Regions in which discrepancies were observed are pointed out, areas for algorithm improvement are identified and the clinical import of algorithm accuracy is discussed.

Humans↗

Computer graphics tools for radiation treatment planning.

The objective of radiation therapy treatment is to eradicate a cancerous tumor while keeping the damage to nearby healthy organs to a minimum. A variety of tools employing computer graphics exist to aid in the planning and verification of treatments. Three-dimensional (3D) image information available from sources such as computerized tomography (CT) scanners is used to define the sizes, shapes, and spatial locations of the tumor and normal structures in the form of transverse contours. These object definitions are displayed in 3D perspective to enable the determination of the best possible directions from which to aim radiation beams at the tumor. The beams may be shaped to match the outline of the tumor, and their intensities may be modified using compensating devices. The results of calculations done to predict the distribution of radiation dose throughout the body due to a given set-up of beams can be displayed to the user in many ways. Dose may be shown in the form of isodose contours overlaid on transverse CT images, or on reconstructed image planes of arbitrary orientation in space. There are also a number of methods of 3D display; dose can be shown on the surface of objects, or in the form of isodose surfaces relative to anatomical structures. Computer-generated beam film images may be used to verify patient set-up and tumor coverage.

Computer Graphics↗

A comprehensive three-dimensional radiation treatment planning system.

A comprehensive software system has been developed to allow 3-dimensional planning of radiation therapy treatments using the extensive anatomical information made available by imaging modalities such as CT and MR. Biological structures of interest and tumor volumes are defined by outlines drawn on a sequence of CT slices. Beam set-ups may then be determined in three dimensions by displaying the structure contours in a beam's eye view, or in two dimensions using a single CT cut. Each beam defined may be shaped by the specification of block aperture contours, and its intensity may be modified with the use of planar compensators. 3D dose calculation algorithms are discussed. To evaluate the calculation results, dose volume histograms are provided, as well as various types of displays in two and three dimensions, including dose on arbitrarily oriented planes, dose on the surface of anatomical objects, and isodose surfaces. Computer generated beam films are also available as an aid in patient set-up verification. These tools, and others, provide the basis for a comprehensive 3D system that can be used throughout the treatment planning process.

Humans↗

Arbitrary oblique image sections for 3-D radiation treatment planning.

Methods for selecting and computing arbitrary image sections for displaying anatomic and isodose information for three-dimensional treatment planning are investigated. Selection of the desired plane may be made by defining a plane that is perpendicular to an existing image section (called the base image) and passing through a line on the base image. Alternatively, the anatomic structures displayed perspectively in three dimensions as a series of contours that can be rotated and translated may be used to define an arbitrary plane for image reconstruction. The viewing screen is considered to be the plane of interest. As a typical three-dimensional image of 30 to 60 sections requires considerable computer storage (on the order of 25 megabytes), a reconstruction algorithm may need extensive memory space or CPU and disk I/O time. Of the schemes examined, we believe the following is the most efficient. One pair of images is read from the disk at a time in sequence and intersections of the rows of the cutting plane with the box formed by the consecutive images are computed. Pixel values of all points between the given images are computed by interpolation. Special cases, such as the cutting plane being parallel to or coincident with an existing image, must be considered separately.

Algorithms↗

Validity of the concept of separating primary and scatter dose.

The technique of separating dose into primary and scatter components for calculating photon dose distributions is widely used. The primary and scatter dose models ignore ignore the fact that electrons have a finite range. This may be a good approximation for 60Co photons but not for higher energies. For the latter, the range of electrons may be several centimeters in soft tissue and even longer in lungs and will lead to errors in computed dose in regions where electronic equilibrium does not exist. Ignoring the finite range of electrons will affect dose at points such as those near the beam boundaries, near inhomogeneity boundaries, and at bone-soft-tissue interfaces. Other possible problems associated with the definition and use of "primary" and "scatter" dose in dose distribution calculations result from extrapolation of measured data to obtain data for zero and very large field sizes and from the use of these quantities, which are defined for central axis, for points at large distances from the central axis. This paper examines the limits of the validity of these assumptions.

Cobalt Radioisotopes↗

Off-center ratios for three-dimensional dose calculations.

A new method is proposed for computing the off-center ratios (OCR's) in three-dimensional dose calculations. For an open field, the OCR at a point is computed as the product of the primary OCR (POCR) and the boundary factors (BF's). The POCR describes the beam profile for an infinite field, that is, without the effect of the collimators. It is defined as the ratio of the dose at a point off the central ray to the dose at the point on the central ray at the same depth for an infinite field. The POCR is a function of radial distance from the beam central ray and depth. The BF describes the shape of the beam in the neighborhood of the field boundary defined by the collimators. It is defined as the ratio of the OCR at a point for a finite field to the OCR at the same point for an infinite field. The BF is a function of distance from the field boundary, depth, and field size. For a wedged field, we assume that the boundary factors remain the same as for open fields but the POCR's are altered. The changes in beam profiles are described by a factor called the wedge profile factor (WPF), defined as the ratio of the dose at a point for the largest wedged field to the dose at the same point for an open field of the same field size. The WPF is a function of lateral distance from the beam central plane and depth. Calculated OCR's using this new method are in agreement with the measured data along both the transverse and the diagonal directions of the field.

Biometry↗

Use of fast Fourier transforms in calculating dose distributions for irregularly shaped fields for three-dimensional treatment planning.

In three-dimensional radiation treatment planning, essentially all fields are irregular and compensated. Consequently, it is important to predict accurately dose for such fields to ensure adequate coverage of the target region and sparing of healthy tissues. Traditional approaches, namely, those involving scatter integration and extended source and those utilizing negatively weighted fields, are inaccurate, especially near the boundaries defined by blocks and collimators. In the method presented in this paper, dose distributions for arbitrarily shaped beams are calculated by two-dimensional convolution of the relative primary photon fluence distributions and kernels representing the cross-sectional profiles of a pencil beam at a series of depths. The pencil beam dose distributions are computed, once and for all, with the Monte Carlo method for photon energy spectrum for each treatment machine. The finite size of the source, which is important for cobalt machines, is also taken into account using convolution of the source with the relative primary fluence distribution. Convolutions are performed using fast Fourier transforms on an array processor. Results of calculations are in excellent agreement with measured data. While no data are presented for fields modified by compensators, the method of calculation should apply at least as well for such fields since the variations in fluence distribution for compensated fields are not as sharp as for points near the block boundaries.

Biophysical Phenomena↗

Dose computations for asymmetric fields defined by independent jaws.

Asymmetric fields defined by independent jaws can be used to split a beam or to match adjacent fields. We have extended a method originally developed for symmetric fields to calculate the dose for asymmetric fields. The dose to a point is computed as the product of the tissue maximum ratio (TMR), the off center ratio (OCR), and the inverse square factor. The TMR is computed from the measured central axis depth doses for symmetric fields. The OCR is obtained by multiplying the primary OCR (POCR) and the boundary factors (BF's) for the four jaws. The POCR's and BF's were derived from measured beam profiles, which include the effect of off-axis beam quality variations. Using this method, the beam profiles and isodose distributions for asymmetric fields of a 6-MV accelerator were calculated and compared with the measured data. The agreement is within experimental errors both in the penumbra region and along the central ray of the asymmetric field.

Humans↗

The effect of angular spread on the intensity distribution of arbitrarily shaped electron beams.

Knowledge of the relative intensity distribution at the patient's surface is essential for pencil beam calculations of three-dimensional dose distributions for arbitrarily shaped electron beams. To calculate the relative intensity distribution, the spatial spread resulting from angular spread is convolved with a two-dimensional step function whose shape corresponds to the applicator aperture. Two different approaches to obtain angular spread or the equivalent spatial spread are investigated. In the first method, the pencil beam angular spread is assumed to be Gaussian in shape. The angular spread constants (sigma theta) are then obtained from the slopes of measured intensity profiles. In the second method, the angular spread, in the form of an array of numerical values, is obtained by the deconvolution of measured intensity profiles. After obtaining the angular spread, the calculation for convolution is done in a number of parallel planes normal to the central axis at various distances from the electron collimator. Intensity at any arbitrary point in space is computed by interpolating between intensity distributions in adjacent planes on either side of the point. The effects of variations in angular spread as a function of field size for two treatment machines, one with a scanned electron beam and the other with a scattering foil, have been studied. The consequences of assuming angular spread to be of Gaussian shape are also examined. The electron intensity calculation techniques described in this paper apply primarily to methods of dose calculations that employ pencil beams generated using Monte Carlo simulations.

Electrons↗

Extraction of pencil beam kernels by the deconvolution method.

A method has been developed to extract pencil beam kernels from measured broad beam profiles. In theory, the convolution of a symmetric kernel with a step function will yield a function that is symmetric about the inflection point. Conversely, by deconvolution, the kernel may be extracted from a measured distribution. In practice, however, due to the uncertainties and errors associated with the measurements and due to the singularities produced in the fast Fourier transforms employed in the deconvolution process, the kernels thus obtained and the dose distributions calculated therefrom, often exhibit erratic fluctuations. We propose a method that transforms measured profiles to new, modified distributions so that they satisfy the theoretical symmetry condition. The resultant kernel from the deconvolution is then free of fluctuations. We applied this method to compute photon and electron dose distributions at various depths in water and electron fluence distributions in air. The agreement between measured and computed profiles is within 1% in dose or 1 mm in distance in high dose gradient regions.

Fourier Analysis↗

Beam characteristics of a new model of 6-MV linear accelerator.

This paper describes the beam characteristics and dosimetry measurements performed on the 6-MV photon beam of a new model of linear accelerator, three of which were recently introduced and installed in our institution. Percent depth dose and tissue maximum ratio tables for a variety of field sizes and depths, as well as other parameters used for treatment planning are presented. These accelerators are the first of their kind using both hardware and software tools to control interlocks. Checking procedures for these interlocks are available from the authors upon request. Comparison of characteristic parameters between these three new 6-MV linear accelerators and with the 6-MV beams of two other accelerators is also made.

Calibration↗

Diode dosimetry of models 6711 and 6712 125I seeds in a water phantom.

Two-dimensional relative dose distributions have been measured around 125I brachytherapy seeds. The two seed models studied, models 6711 and 6712, were manufactured by the 3M Company. Silicon detectors immersed in water phantoms were used to measure the dose. A computerized data acquisition system that controlled the radial position of the diode and the angular rotation of the seed, as well as a manually controlled system were used to collect and store the data. Our results show that the two seed models have relative dose distributions which are quite similar; however, the absolute dose distributions are sufficiently different to warrant separate look-up tables for the two seed models. Additionally, our results are compared with dose distribution data previously obtained for the model 6711 seed.

Brachytherapy↗