Use of the fast Hartley transform for three-dimensional dose calculation in radionuclide therapy

Comparison of internal dose estimates obtained using organ-level, voxel S value, and Monte Carlo techniques

PURPOSE

The authors' objective was to compare internal dose estimates obtained using the Organ Level Dose Assessment with Exponential Modeling (OLINDA/EXM) software, the voxel S value technique, and Monte Carlo simulation. Monte Carlo dose estimates were used as the reference standard to assess the impact of patient-specific anatomy on the final dose estimate.

METHODS

Six patients injected with 99mTc-hydrazinonicotinamide-Tyr3-octreotide were included in this study. A hybrid planar/SPECT imaging protocol was used to estimate 99mTc time-integrated activity coefficients (TIACs) for kidneys, liver, spleen, and tumors. Additionally, TIACs were predicted for 131I, 177Lu, and 90Y assuming the same biological half-lives as the 99mTc labeled tracer. The TIACs were used as input for OLINDA/EXM for organ-level dose calculation and voxel level dosimetry was performed using the voxel S value method and Monte Carlo simulation. Dose estimates for 99mTc, 131I, 177Lu, and 90Y distributions were evaluated by comparing (i) organ-level S values corresponding to each method, (ii) total tumor and organ doses, (iii) differences in right and left kidney doses, and (iv) voxelized dose distributions calculated by Monte Carlo and the voxel S value technique.

RESULTS

The S values for all investigated radionuclides used by OLINDA/EXM and the corresponding patient-specific S values calculated by Monte Carlo agreed within 2.3% on average for self-irradiation, and differed by as much as 105% for cross-organ irradiation. Total organ doses calculated by OLINDA/EXM and the voxel S value technique agreed with Monte Carlo results within approximately ±7%. Differences between right and left kidney doses determined by Monte Carlo were as high as 73%. Comparison of the Monte Carlo and voxel S value dose distributions showed that each method produced similar dose volume histograms with a minimum dose covering 90% of the volume (D90) agreeing within ±3%, on average.

CONCLUSIONS

Several aspects of OLINDA/EXM dose calculation were compared with patient-specific dose estimates obtained using Monte Carlo. Differences in patient anatomy led to large differences in cross-organ doses. However, total organ doses were still in good agreement since most of the deposited dose is due to self-irradiation. Comparison of voxelized doses calculated by Monte Carlo and the voxel S value technique showed that the 3D dose distributions produced by the respective methods are nearly identical.

Overview of dosimetry for systemic targeted radionuclide therapy (STaRT)

The purposes of systemic targeted radionuclide therapy dosimetry include compiling a database of normal organ radiation-absorbed doses that are carrier- and radionuclide-specific, and assuring that the normal organ radiation doses are within a safe range before therapy. Also of importance is quantitation of radiation delivery to tumors vs. normal tissues to correlate absorbed dose with tumor control. For agents with significant and variable excretion, estimates of individual patient distribution/clearance may be needed to optimize the dose-response relationship.

Validation of the EGS usercode DOSE3D for internal beta dose calculation at the cellular and tissue levels

Internal radiotherapy is currently focusing on beta emitters such as 90Y or 131I because of their high-energy emissions. However, conventional dosimetric methods (MIRD) are known to be limited for such applications. They are unable to take into account microscopic radionuclide distribution because standardized anthropomorphic phantoms are used, and absorbed dose is calculated at the organ level. New tools are therefore required for dose assessment at cellular and tissue level (10-100 microm). The purpose of this study was to validate, at this scale, a Monte Carlo usercode (DOSE3D), based on the MORSE combinatorial geometry package and the EGS code system. Dose point-kernel calculations in water were compared to those published by Cross et al and Simpkin and Mackie. They confirm that DOSE3D is a reliable tool for cellular dosimetry in various geometric configurations.

A software tool for specifying voxel models for dosimetry estimation

Many dose estimation problems can be conveniently formulated in terms of finding the energy emitted and absorbed by a set of homogeneous volume elements (voxels) arranged in a rectilinear grid. The solution of these problems requires an accurate model of the source and target geometry to be established, whereupon conventional Monte Carlo simulation of radiation transport can be employed to determine energy deposition. A software application ("MrVoxel") has been developed to assist in the specification of the source and target models. This application includes tools for image segmentation and image registration (2D and 3D, intra- and inter-modality, interactive, and automatic). It employs a plug-in architecture to facilitate customization and future expansion: plug-ins can be written to perform image import and export as well as to implement specialized image processing routines. Using plug-ins, the package can, for example, import DICOM 3.10 files and export input files for a voxel-based Monte Carlo package. Standard dosimetric tools such as the geometric mean method, transmission based attenuation correction, and MIRD-style voxel dose kernel convolution are also implemented as plug-ins. MrVoxel was implemented on a Macintosh computer using a commercial software framework to produce a conventional document-centric application. Hence it includes useful features such as the ability to undo an operation or to save a processed image set at any point. This latter feature enables the production of a processing trail, to allow post-hoc auditing of the analysis process. This paper describes the MrVoxel application and its role in the analysis of a particular dosimetry problem.

Treatment planning for molecular targeted radionuclide therapy

Molecular targeted radionuclide therapy promises to expand the usefulness of radiation to successfully treat widespread cancer. The unique properties of radioactive tags make it possible to plan treatments by predicting the radiation absorbed dose to both tumors and normal organs, using a pre-treatment test dose of radiopharmaceutical. This requires a combination of quantitative, high-resolution, radiation-detection hardware and computerized dose-estimation software, and would ideally include biological dose-response data in order to translate radiation absorbed dose into biological effects. Data derived from conventional (external beam) radiation therapy suggests that accurate assessment of the radiation absorbed dose in dose-limiting normal organs could substantially improve the observed clinical response for current agents used in a myeloablative regimen, enabling higher levels of tumor control at lower tumor-to-normal tissue therapeutic indices. Treatment planning based on current radiation detection and simulations technology is sufficient to impact on clinical response. The incorporation of new imaging methods, combined with patient-specific radiation transport simulations, promises to provide unprecedented levels of resolution and quantitative accuracy, which are likely to increase the impact of treatment planning in targeted radionuclide therapy.

Dose-volume histograms computation comparisons using conventional methods and optimized fast Fourier transforms algorithms for brachytherapy

In anatomy based optimization procedures for large volume implants the calculation of dose-volume histograms (DVH) accounts for the major part of the time involved and can be as long as a few hours. This time is proportional to the number of seeds or source dwell positions required for the implant. A procedure for the calculation of brachytherapy seed dose distribution calculation employing fast Fourier transforms (FFT) and the convolution theorem has been described by others and was supposed to significantly improve the speed of the dose distribution computation. Using new significantly improved FFT algorithms and various other optimization techniques we have compared the calculated differential and integral DVHs in high dose rate (HDR) brachytherapy with a single stepping source using actual clinical implants. This is so that we could assess the efficiency and accuracy of the FFT method with that of conventional methods. Our results showed that the FFT based method of calculating DVHs in brachytherapy is comparable in speed with conventional dose calculation methods, but only for implants with more than 287 sources. It is therefore of limited practical use even for large implants. This result is in direct opposition to the claim by other authors.