A continuous penalty function method for inverse treatment planning.

D H Hristov; B G Fallone

doi:10.1118/1.598183

Home/Archive/PubMed

PubMed · 9507482

A continuous penalty function method for inverse treatment planning.

D H Hristov·B G Fallone

View original source ↗Export metadata ↓

Abstract

Conventional inverse treatment planning attempts to calculate dose distributions that may not be feasible given the specified dose levels to various anatomical structures. A technique for inverse treatment planning has been developed that uses only target dose levels which are easily selectable to be feasible. A nonlinear constrained minimization problem is formulated to reflect the goal of sparing critical organs as much as possible while delivering a certain target dose within specified uniformity. The objective function is the squared dose delivered to critical organs. The constraints require the delivery of certain target dose within specified uniformity and non-negative pencil beam weights. A continuous penalty function method is introduced as a method for solving the large-scale constrained minimization problem. The performance of the continuous penalty function method is optimized by numerical investigation of few numerical integration schemes and a pair of weighting functions which influence the utility of the method. Clinical examples are presented that illustrate several features of the technique. The properties of the continuous penalty function method suggest that it may be a viable alternative to conventional inverse treatment planning.

Explore related subjects

Biophysical Phenomena Biophysics Humans Models, Theoretical Neoplasms Phantoms, Imaging Photons Radiotherapy Planning, Computer-Assisted Radiotherapy, High-Energy Journal Article Research Support, Non-U.S. Gov't

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

D H Hristov, B G Fallone. 1998. A continuous penalty function method for inverse treatment planning.. https://doi.org/10.1118/1.598183

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

KEEP EXPLORING

Related citations

PubMed2007

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↗

PubMed2007

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↗

PubMed2007

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↗