Biomedical subjects

C S Chui

Publications and source records attributed to C S Chui.

At least 19 recordsLinked to original sources

PubMed1998

A gradient inverse planning algorithm with dose-volume constraints.

An inverse planning algorithm for determining the intensity-modulated beams that will most closely generate a desired dose distribution is presented. The algorithm is three-dimensional and does not explicitly depend on beam energies and modalities. It allows a single prescription dose or a window of acceptable doses to be specified for the target, with additional constraints to account for under- or over-dosing. For the protection of organs at risk, it provides maximum-dose and dose-volume constraints. The latter apply to the entire volume of the organ exposed to the corresponding dose levels. Several levels of each type of constraint, with varying penalty weights, may be specified for each organ. The objective function that serves as the measure of the goodness of the solution is of the least-squares type and is minimized using conjugate gradient methods. Typical clinical cases involving 40,000 points and 4000 rays to be determined require about 10 min of CPU time on a DEC AlphaStation. Results are presented for two clinical sites, prostate and lung. The optimization algorithm yielded plans that featured higher target dose homogeneity, compared with the human planner's plan, while selectively sparing more of the normal organs at the desired dose regions.

Aged↗

PubMed1998

A patient-specific Monte Carlo dose-calculation method for photon beams.

A patient-specific, CT-based, Monte Carlo dose-calculation method for photon beams has been developed to correctly account for inhomogeneity in the patient. The method employs the EGS4 system to sample the interaction of radiation in the medium. CT images are used to describe the patient geometry and to determine the density and atomic number in each voxel. The user code (MCPAT) provides the data describing the incident beams, and performs geometry checking and energy scoring in patient CT images. Several variance reduction techniques have been implemented to improve the computation efficiency. The method was verified with measured data and other calculations, both in homogeneous and inhomogeneous media. The method was also applied to a lung treatment, where significant differences in dose distributions, especially in the low-density region, were observed when compared with the results using an equivalent pathlength method. Comparison of the DVHs showed that the Monte Carlo calculated plan predicted an underdose of nearly 20% to the target, while the maximum doses to the cord and the heart were increased by 25% and 33%, respectively. These results suggested that the Monte Carlo method may have an impact on treatment designs, and also that it can be used as a benchmark to assess the accuracy of other dose calculation algorithms. The computation time for the lung case employing five 15-MV wedged beams, with an approximate field size of 13 X 13 cm and the dose grid size of 0.375 cm, was less than 14 h on a 175-MHz computer with a standard deviation of 1.5% in the high-dose region.

Algorithms↗

PubMed1998

Physical and dosimetric aspects of a multileaf collimation system used in the dynamic mode for implementing intensity modulated radiotherapy.

The use of a multileaf collimator in the dynamic mode to perform intensity modulated radiotherapy became a reality at our institution in 1995. Unlike treatment with static fields using a multileaf collimator, there are significant dosimetric issues which must be assessed before dynamic therapy can be implemented. We have performed a series of calculations and measurements to quantify head scatter for small fields, collimator transmission, and the transmission through rounded leaf ends. If not accounted for, these factors affect the delivered dose to the prostate by 5%-20% for a typical plan. Data obtained with ion chambers and radiographic film are presented for both 6 and 15 MV x-ray beams. The impact on the delivered dose of the mechanical accuracy of the multileaf collimator, achieved during leaf position calibration and maintained during dose delivery, is also discussed.

Biophysical Phenomena↗

PubMed1997

Total body irradiation with an arc and a gravity-oriented compensator.

PURPOSE: To deliver uniform dose distributions for total-body irradiation (TBI) with an arc field and a gravity-oriented compensator. This technique allows the patient to be treated lying on the floor in a small treatment room. METHODS AND MATERIALS: Through the sweeping motion of the gantry, a continuous arc field can deliver a large field to a patient lying on the floor. The dose profile, however, would not be uniform if no compensator were used, due to the effects of inverse square variation of beam intensity with distance as well as the slanted depth in patient. To solve this problem, a gravity-oriented compensator made of cerrobend alloy was designed. This compensator has a cross-section of an inverted isosceles triangle, with the apex always pointing downward, due to gravity. By properly selecting the thickness of the compensator, the width of the base, and the distance between the pivots to the base, the difference in the path length through the compensator can be made just right to compensate the effects of inverse-square and slanted depth, thus producing a uniform dose profile. RESULTS: Arc fields with a gravity-oriented compensator were used for 6, 10, 15, and 18 MV photon beams. The arc field can cover a patient with a height up to 180 cm. The field width was chosen from 32 to 40 cm at the machine isocenter. The optimal thickness of the compensator was found to be 2.5 cm, and its base was 25 cm wide. The distance from the pivot points to the flat surface of the compensator proximal to the beam ranges from 13 to 14 cm for different beam energies. The dose uniformity at a depth of 10 cm is within +/-5% for all beam energies used in this study. CONCLUSIONS: Highly uniform dose profiles for TBI treatments can be delivered with an arc and a gravity-oriented compensator. The proposed technique is simple and versatile. A single compensator can be used for all energies, because the amount of compensation can be adjusted by changing the distance to the pivot and/or the field size.

Equipment Design↗

PubMed1997

Planning, delivery, and quality assurance of intensity-modulated radiotherapy using dynamic multileaf collimator: a strategy for large-scale implementation for the treatment of carcinoma of the prostate.

PURPOSE: To improve the local control of patients with adenocarcinoma of the prostate we have implemented intensity modulated radiation therapy (IMRT) to deliver a prescribed dose of 81 Gy. This method is based on inverse planning and the use of dynamic multileaf collimators (DMLC). Because IMRT is a new modality, a major emphasis was on the quality assurance of each component of the process and on patient safety. In this article we describe in detail our procedures and quality assurance program. METHODS AND MATERIALS: Using an inverse algorithm, we have developed a treatment plan consisting five intensity-modulated (IM) photon fields that are delivered with DMLC. In the planning stage, the planner specifies the number of beams and their directions, and the desired doses for the target, the normal organs and the "overlap" regions. Then, the inverse algorithm designs intensity profiles that best meet the specified criteria. A second algorithm determines the leaf motion that would produce the designed intensity pattern and produces a DMLC file as input to the MLC control computer. Our quality assurance program for the planning and treatment delivery process includes the following components: 1) verification of the DMLC field boundary on localization port film, 2) verification that the leaf motion of the DMLC file produces the planned dose distribution (with an independent calculation), 3) comparison of dose distribution produced by DMLC in a flat phantom with that calculated by the treatment planning computer for the same experimental condition, 4) comparison of the planned leaf motions with that implemented for the treatment (as recorded on the MLC log files), 5) confirmation of the initial and final positions of the MLC for each field by a record-and-verify system, and 6) in vivo dose measurements. RESULTS: Using a five-field IMRT plan we have customized dose distribution to conform to and deliver 81 Gy to the PTV. In addition, in the overlap regions between the PTV and the rectum, and between the PTV and the bladder, the dose is kept within the tolerance of the respective organs. Our QA checks show acceptable agreement between the planned and the implemented leaf motions. Correspondingly, film and TLD dosimetry indicates that doses delivered agrees with the planned dose to within 2%. As of September 15, 1996, we have treated eight patients to 81 Gy with IMRT. CONCLUSION: For complex planning problems where the surrounding normal tissues place severe constraints on the prescription dose, IMRT provides a powerful and efficient solution. Given a comprehensive and rigorous quality-assurance program, the intensity-modulated fields can be efficaciously and accurately delivered using DMLC. IMRT treatment is now ready for routine implementation on a large scale in our clinic.

Adenocarcinoma↗

PubMed1997

Lying-on position of total skin electron therapy.

PURPOSE: A technique for whole-body electron therapy with the patient in a lying position has been developed. This technique allows Total Skin Electron Therapy (TSET) for those patients who were previously unable to be treated in a conventional standing position. METHODS AND MATERIALS: This study was carried out on a Varian 2100C linear accelerator with a 6 MeV high dose rate electron beam. The collimator was open to a width of 36 x 36 cm. There were two main procedures, with six dual-field techniques: 1) two static AP/PA vertical dual fields (VDF): the patient laid on the floor transversely under the collimator when the gantry was in a vertical position. A 0.6 cm acrylic board was placed 15 cm away from the patient, then the gantry was rotated 25 degrees clockwise and counterclockwise to treat the patient in the supine and prone positions, respectively. 2) Four oblique junction fields (OJF): the patient laid on the floor in a prone and supine position parallel to the wave guide at (227 - body thickness x tan 60 degrees) cm away from the vertical axis of the gantry, then the gantry was rotated 60 degrees toward the patient. A 0.6 cm acrylic board was placed 15 cm away from the patient perpendicular to the beam. The patient was move along the field central axis. It allowed the patient's body to be within the 160 cm effective treatment profile. When the patient's body axis move 5 degrees toward the lateral side of the field central axis, we could obtain a better dose distribution in the vertex of the scalp and the soles of the feet. The angle of the VDF was measured by chamber detectors to obtain the effective treatment profile. Likewise, the optimal profile for the OJF was determined by the same procedures. The Rando phantom was used to measure the superficial dose of the body. RESULTS: The dimension of effective treatment profile for the VDF was 188 x 72 cm at 87% dose level For the OJF, we had to move the patient along the field central axis to obtain the effective treatment profile in a 180 x 85 cm dimension at a 87% dose level. The vertex and sole dose measured in this setup was in the range of 80-88%. CONCLUSIONS: The empirical data showed that the lying-on position for TSET was technically feasible. The dose distribution in the body surface was also compatible with the Stanford standing technique. The nonambulatory skin malignancy patient can be treated in a comfortable and reproducible position.

Electrons↗

PubMed1997

Implementation of a Monte Carlo dosimetry method for patient-specific internal emitter therapy.

In internal emitter therapy, an accurate description of the absorbed dose distribution is necessary to establish an administered dose-response relationship, as well as to avoid critical organ toxicity. This work describes the implementation of a dosimetry method that accounts for the radionuclide decay spectrum, and patient-specific activity and density distributions. The dosimetry algorithm is based on a Monte Carlo procedure that simulates photon and electron transport and scores energy depositions within the patient. The necessary input information may be obtained from a registered set of CT and SPECT or PET images. The algorithm provides the absorbed dose rate for the radioactivity distribution provided by the SPECT or PET image. The algorithm was benchmarked by reproducing dosimetric quantities using the Medical Internal Radionuclide Dose (MIRD) Committee's Standard Man phantom and was used to calculate absorbed dose distributions for representative case studies.

Humans↗

PubMed1996

Conformal radiation treatment of prostate cancer using inversely-planned intensity-modulated photon beams produced with dynamic multileaf collimation.

PURPOSE: To implement radiotherapy with intensity-modulated beams, based on the inverse method of treatment design and using a multileaf collimation system operating in the dynamic mode. METHODS AND MATERIALS: An algorithm, based on the inverse technique, has been integrated into the radiotherapy treatment-planning computer system in our Center. This method of computer-assisted treatment design was used to derive intensity-modulated beams to optimize the boost portion of the treatment plan for a patient with a T1c cancer of the prostate. A dose of 72 Gy (in 40 fractions) was given with a six-field plan, and an additional 9 Gy (in five fractions) with six intensity-modulated beams. The intensity-modulated fields were delivered using dynamic multileaf collimation, that is, individual leaves were in motion during radiation delivery, with the treatment machine operating in the clinical mode. Exhaustive quality assurance measurement and monitoring were carried out to ensure safe and accurate implementation. RESULTS: Dose distribution and dose-volume histogram of the "inverse method" boost plan and of the composite (72 Gy primary + 9 Gy boost) plan were judged clinically acceptable. Compared to a manually designed boost plan, the inverse treatment design gave improved conformality and increased dose homogeneity in the planning target volume. Film and ion chamber dosimetry, performed prior to the first treatment, indicated that each of the six intensity-modulated fields was accurately produced. Thermoluminescent dosimeter (TLD) measurements performed on the patient confirmed that the intended dose was delivered in the treatment. In addition, computer-aided treatment-monitoring programs assured that the multileaf collimator (MLC) position file was executed to the specified precision. In terms of the overall radiation treatment process, there will likely be labor savings in the planning and the treatment phases. CONCLUSIONS: We have placed into clinical use an integrated system of conformal radiation treatment that incorporated the inverse method of treatment design and the use of dynamic multileaf collimation to deliver intensity-modulated beams. The system can provide better treatment design, which can be implemented reliably and safely. We are hopeful that improved treatment efficacy will result.

Humans↗

PubMed1996

The use of diode dosimetry in quality improvement of patient care in radiation therapy.

The purpose of this work is to improve the quality of patient care in radiation therapy by implementing a comprehensive quality assurance (QA) program aiming to enhance patient in vivo dosimetry on a routine basis. The characteristics of two commercially available semi-conductor diode dosimetry systems were evaluated. The diodes were calibrated relative to an ionization chamber-electrometer system with calibrations traceable to the National Institute of Standards and Technology (NIST). Correction factors of clinical relevance were quantified to convert the diode readings into patient dose. The results of dose measurements on 6 patients undergoing external beam radiation therapy for carcinoma of the prostate on three different therapy units are presented. Field shaping during treatments was accomplished either by multileaf collimation or by cerrobend blocking. A deviation of less than +/-4% between the measured and prescribed patient doses was observed. The results indicate that the diodes exhibit excellent linearity, dose reproducibility, minimal anisotropy, and can be used with confidence for patient dose verification. Furthermore, diodes render real time verification of dose delivered to patients.

Calibration↗

PubMed1996

Dosimetric verification of intensity-modulated fields.

The optimization of intensity distributions and the delivery of intensity-modulated treatments with dynamic multi-leaf collimators (MLC) offer important improvements to three-dimensional conformal radiotherapy. In this study, a nine-beam intensity-modulated prostate plan was generated using the inverse radiotherapy technique. The resulting fluence profiles were converted into dynamic MLC leaf motions as functions of monitor units. The leaf motion pattern data were then transferred to the MLC control computer and were used to guide the motions of the leaves during irradiation. To verify that the dose distribution predicted by the optimization and planning systems was actually delivered, a homogeneous polystyrene phantom was irradiated with each of the nine intensity-modulated beams incident normally on the phantom. For each exposure, a radiographic film was placed normal to the beam in the phantom to record the deposited dose. The films were calibrated and scanned to generate 2-D isodose distributions. The dose was also calculated by convolving the incident fluence pattern with pencil beams. The measured and calculated dose distributions were compared and found to have discrepancies in excess of 5% of the central axis dose. The source of discrepancies was suspected to be the rounded edges of the leaves and the scattered radiation from the various components of the collimation system. After approximate corrections were made for these effects, the agreement between the two dose distributions was within 2%. We also studied the impact of the "tongue-and-groove" effect on dynamic MLC treatments and showed that it is possible to render this effect inconsequential by appropriately synchronizing leaf motions. This study also demonstrated that accurate and rapid delivery of realistic intensity-modulated plans is feasible using a dynamic multi-leaf collimator.

Biophysical Phenomena↗

PubMed1996

Radionuclide photon dose kernels for internal emitter dosimetry.

Photon point dose kernels and absorbed fractions were generated in water for the full photon emission spectrum of each radionuclide of interest in nuclear medicine, by simulating the transport of particles using Monte Carlo. The kernels were then fitted to a mathematical expression. Absorbed fractions for point sources were obtained by integrating the kernels over spheres. Photon dose kernels and absorbed fractions were generated for the following radionuclides: I-123, I-124, I-125, I-131, In-111, Cu-64, Cu-67, Ga-67, Ga-68, Re-186, Re-188, Sm-153, Sn-117m, Tc-99m. The Monte Carlo simulation was verified by comparing the dose kernels to published monoenergetic photon kernels. Further validation was obtained by generating an I-125 brachytherapy seed kernel and comparing it with published data. Since Monte Carlo simulation was initialized by sampling from the complete photon spectra of these radionuclides, interpolation between monoenergetic kernels and absorbed fractions was not required. The absorbed-fraction due to uniform spherical distributions can be directly applied for use in internal dosimetry. In addition, the kernels can be used as input for three-dimensional internal dosimetry calculations.

Brachytherapy↗

PubMed1996

Generation of arbitrary intensity profiles by combining the scanning beam with dynamic multileaf collimation.

An algorithm, which combines the scanning beam with dynamic collimation to generate any arbitrary intensity profile, is presented. The desired intensity profile is assumed to be piecewise linear. The dynamic collimation method used is the "sliding window." The algorithm can be used either for a given scanning beam profile or to simultaneously determine the scanning beam profile and the leaf motions required to generate the desired intensity profile, which minimize the total treatment time. The limitations imposed by the physics of an elementary beam are taken into account. The algorithm is an iterative one, with typical calculation times being of the order of a few milliseconds.

Algorithms↗

PubMed1996

Testing of dynamic multileaf collimation.

It has been shown that intensity-modulated fields have the potential to deliver optimum dose distributions, i.e., high dose uniformity in the target and lower doses in the surrounding critical organs. One way to deliver such fields is by using dynamic multileaf collimation (DMLC). This capability is already available in research mode on some treatment machines. While much effort has been devoted to developing algorithms for DMLC, the mechanical reliability of this new treatment delivery mode has not been fully studied. In this work, we report a series of tests designed to investigate the mechanical aspects of DMLC and their implications on dosimetry. Specifically, these tests were designed to examine (1) the stability of leaf speed, (2) the effect of lateral disequilibrium on dose profiles between adjacent leaves, (3) the significance of acceleration and deceleration of leaf motion, (4) the effect of positional accuracy and rounded-end of the leaves, and (5) create a simple test pattern that may serve as a basis for routine quality assurance checks. Results of these tests are presented. The implications on dosimetry and consideration for the design of leaf motion are discussed.

Humans↗

PubMed1996

A Monte Carlo approach to patient-specific dosimetry.

In internal emitter therapy, an accurate description of the absorbed dose distribution is necessary to establish an administered dose-response relationship, as well as to avoid critical organ toxicity. Given a spatial distribution of cumulated activity, an absorbed dose distribution that accounts for the effects of attenuation and scatter can be obtained using a Monte Carlo method that simulates particle transport across the various densities and atomic numbers encountered in the human body. Patient-specific information can be obtained from CT and SPECT or PET imaging. Since the data from these imaging modalities is discrete, it is necessary to develop a technique to efficiently transport particles across discrete media. The Monte Carlo-based algorithm presented in this article produces accurate absorbed dose distributions due to patient-specific density and radionuclide activity distributions. The method was verified by creating CT and SPECT arrays for the Medical Internal Radionuclide Dose (MIRD) Committee's Standard Man phantom, and reproducing the spatially averaged specific absorbed fractions reported in MIRD Pamphlet 5. The algorithm was used to investigate the implications of replacing a mean absorbed dose with a distribution, and of neglecting atomic number and density variations for various patient geometries and energies. For example, the I-131 specific absorbed fraction for spleen to liver is the same as for liver to spleen, yet the distributions were different. Furthermore, neglecting atomic number variations across the vertebral bone led to an overestimation of I-125 absorbed dose by an order of magnitude, while no error was observed for I-131.

Algorithms↗

PubMed1995

Effect of ocular implants of different materials on the dosimetry of external beam radiation therapy.

PURPOSE: To study the attenuation and scattering effects of ocular implants, made from different materials, on the dose distributions of a 6 MV photon beam, and 6, 9, and 12 MeV electron beams used in orbital radiotherapy. METHODS AND MATERIALS: Central axis depth-dose measurements were performed in a polystyrene phantom with embedded spherical ocular implants using film dosimetry of a 6 MV photon beam and electron beams of 6, 9, and 12 MeV energy. The isodose distributions were also calculated by a computerized treatment planning system using computerized tomography (CT) scans of a polystyrene phantom that had silicone, acrylic, and hydroxyapatite ocular implants placed into it. RESULTS: Electron beam dose distributions display distortions both on the measured and calculated data. This effect is most accentuated for the hydroxyapatite implants, for which the transmissions through ocular implants are on the order of 93% for the 6 MV photon beam, and range from 60% for 6 MeV electrons to 90% for 12 MeV electrons. CONCLUSION: We studied the effect of ocular implants of various materials, embedded in a polystyrene phantom, on the dose distributions of a 6 MV photon beam, and 6, 9, and 12 MeV electron beams. Our investigations show that while 6 MV photons experience only a few percent attenuation, lower energy electron beam with 60% transmission is not a suitable choice of treating tumors behind the ocular implants.

Colloids↗

PubMed1995

Mean mass energy absorption coefficient ratios for megavoltage x-ray beams.

Mean mass energy absorption coefficient ratios of acrylic, polystyrene, and water to air, were calculated using Monte Carlo generated energy spectra. The energy spectra were calculated for 4- to 50-MV x-ray beams, from machines using flattening filters and scanning beams. The validity of these spectra was verified by comparing the measured ionization ratios with the calculated values. The agreement was found to be within 1.9%. For beams of energy below 6 MV, our estimates of the mean mass energy absorption coefficient ratios agree well with those recommended by the TG-21 protocol. For higher energy beams, the discrepancy increases to about 3%. It was found that the discrepancy is attributable to the different spectra used in these calculations.

Humans↗

PubMed1995

Beam characteristics of a new generation 50 MeV racetrack microtron.

The first of a new generation of microtron accelerators has been installed and tested. It is currently in use for multisegment conformal radiotherapy at our institution. The unit produces x rays and electrons from 10 to 50 MeV in 5 MeV increments. It incorporates a 64 leaf, doubly focused multileaf collimator (MLC), which can be used to shape x-ray and electron beams. Both x-ray and electron beams are produced by magnetically scanning the electron beams from the accelerator. The new generation unit incorporates a purging magnet to sweep away any primary or secondary electrons that pass through the target(s). In this paper, the beam characteristics of the accelerator that were studied during acceptance testing are described. Representative examples of depth doses, beam profiles, output factors, and elementary beam distributions are presented and discussed, in comparison with the earlier generation of microtron accelerators and with other radiotherapy machines.

Humans↗

PubMed1995

A Monte Carlo model of photon beams used in radiation therapy.

A generic Monte Carlo model of a photon therapy machine is described. The model, known as McRad, is based on EGS4 and has been in use since 1991. Its primary function has been the characterization of the incident photon fluence for use by dose calculation algorithms. The accuracy of McRad is examined by comparing the dose distributions in a water phantom generated using only the Monte Carlo data with measured dose distributions for two machines in our clinic; a 6 MV Varian Clinac 600C and the 15 MV beam from a Clinac 2100C. The Monte Carlo generated dose distributions are computed using a dose calculation algorithm based on the use of differential pencil beam kernels. It was found that the match to measured data could be improved if the model is tuned by adjusting the energy of the electron beam incident on the target. The beam profiles were found to be more sensitive indicators of the electron beam energy than the depth dose curves. Beyond the depths reached by contaminant electrons, the computed and measured depth dose curves agree to better than 1%. The comparison of beam profiles indicate that in regions up to within 1 cm of the field edge, the measured and computed doses generally agree to within 2%-3%.

Algorithms↗