Search PubMed⌕ Search

PubMed · 8994174

A source model for efficient brachytherapy computations with Monte Carlo.

Abstract

Monte Carlo techniques have the potential for producing accurate brachytherapy dose distributions in heterogeneous finite geometries. However, for routine clinical use, computational speed must be adequate. A fast, all-particle, CT-based Monte Carlo code called PEREGRINE is being developed at Lawrence Livermore National Laboratory for radiation treatment planning. As one feature, the code will produce accurate dose distributions from brachytherapy sources in heterogeneous geometries. For efficiency, brachytherapy sources in this model are treated as points or line segments. Radiation is emitted with the proper energy spectrum and (perhaps anisotropic) angular distribution. In particular, for anisotropic emission the polar angle is determined by a random-number driven empirical function constructed from a source's measured or precomputed fluence emission pattern. Source model parameters are presented for iodine and iridium sources. While designed for the PEREGRINE program, this source model can be used in any Monte Carlo code.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

K Weaver, C H Siantar, W Chandler, R M White. 1996. A source model for efficient brachytherapy computations with Monte Carlo.. https://doi.org/10.1118/1.597890

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

Searching for the best model: ambiguity of inverse solutions and application to fetal magnetoencephalography.

Fetal brain signals produce weak magnetic fields at the maternal abdominal surface. In the presence of much stronger interference these weak fetal fields are often nearly indistinguishable from noise. Our initial objective was to validate these weak fetal brain fields by demonstrating that they agree with the electromagnetic model of the fetal brain. The fetal brain model is often not known and we have attempted to fit the data to not only the brain source position, orientation and magnitude, but also to the brain model position. Simulation tests of this extended model search on fetal MEG recordings using dipole fit and beamformers revealed a region of ambiguity. The region of ambiguity consists of a family of models which are not distinguishable in the presence of noise, and which exhibit large and comparable SNR when beamformers are used. Unlike the uncertainty of a dipole fit with known model plus noise, this extended ambiguity region yields nearly identical forward solutions, and is only weakly dependent on noise. The ambiguity region is located in a plane defined by the source position, orientation, and the true model centre, and will have a diameter approximately 0.67 of the modelled fetal head diameter. Existence of the ambiguity region allows us to only state that the fetal brain fields do not contradict the electromagnetic model; we can associate them with a family of models belonging to the ambiguity region, but not with any specific model. In addition to providing a level of confidence in the fetal brain signals, the ambiguity region knowledge in combination with beamformers allows detection of undistorted temporal waveforms with improved signal-to-noise ratio, even though the source position cannot be uniquely determined.

Biophysical Phenomena↗

Recovery of total I-131 activity within focal volumes using SPECT and 3D OSEM.

We experimentally investigated the SPECT recovery of I-131 activity in multiple spheres located simultaneously at different locations within a cylindrical phantom that had an elliptical cross section. The sphere volumes ranged from 209 cc down to 4.2 cc. A Prism 3000 camera and two types of parallel-hexagonal-hole collimation were employed: high energy (HE) and ultra high energy (UHE). Using appropriately-different 3D models of the point source response function for the two types of collimation, approximately the same recovery of activity could be achieved with either collimation by 3D OSEM reconstruction. The recovery coefficient was greater with no background activity in the phantom by 0.10, on average, compared to that with background. In the HE collimation case, the activity recovery was considerably better for all volumes using 3D OSEM reconstruction than it had been in the past using 1D SAGE reconstruction. Recovery-coefficient-based correction in a simulated patient case involving spherical tumours moderately improved the activity estimates (average error reduced from 14% to 9% for UHE collimation, and from 15% to 11% for HE collimation). For a test case with HE collimation, increasing the projection-image sampling density while decreasing the image voxel size increased the recovery coefficient by 0.075 on average, and, if used in a full set of calibration measurements of recovery coefficient versus volume, might lead to further improvement in accuracy for the patient case.

Biophysical Phenomena↗

Improvement of spatial resolution in the longitudinal direction for isotropic imaging in helical CT.

Experiments were conducted to confirm the isotropic spatial resolution of multislice CT with a 0.5 mm slice thickness. Isotropic spatial resolution means that the spatial resolution in the transaxial plane (X-Y plane) and that in the longitudinal direction (Z direction) are equivalent. To obtain point spread function (PSF) values in the X-Y-Z directions, three-dimensional voxel data were obtained by helical scanning of a bead phantom. The modulation transfer function (MTF) values were then obtained by three-dimensional Fourier transform of the PSF. Evaluation of the spatial resolution in the X-Y-Z directions by the MTF values showed that the spatial resolution in the Z direction does not depend on the reconstruction kernel used. It was also found that the spatial resolution in the Z direction, as compared with that in the X-Y plane, is superior with the standard kernel for the abdomen and is inferior with the high-definition kernel for the ears/bones. By performing sharpening filter processing in the Z direction with a high-definition kernel, comparable spatial resolution could be obtained in the X-Y-Z directions. It was confirmed that adjusting the spatial resolution in the Z direction with the reconstruction kernel used is an effective method for isotropic imaging.

Biophysical Phenomena↗