Real-time Dose Calculation Research Articles

Time-resolved clinical dose volume metrics, calculations and predictions based on source tracking measurements and uncertainties to aid treatment verification and error detection for HDR brachytherapy—a proof-of-principle study

Objective. High-dose-rate (HDR) brachytherapy lacks routinely available treatment verification methods. Real-time tracking of the radiation source during HDR brachytherapy can enhance treatment verification capabilities. Recent developments in source tracking allow for measurement of dwell times and source positions with high accuracy. However, more clinically relevant information, such as dose discrepancies, is still needed. To address this, a real-time dose calculation implementation was developed to provide more relevant information from source tracking data. A proof-of-principle of the developed tool was shown using source tracking data obtained from a 3D-printed anthropomorphic phantom. Approach. Software was developed to calculate dose-volume-histograms (DVH) and clinical dose metrics from experimental HDR prostate treatment source tracking data, measured in a realistic pelvic phantom. Uncertainty estimation was performed using repeat measurements to assess the inherent dose measuring uncertainty of the in vivo dosimetry (IVD) system. Using a novel approach, the measurement uncertainty can be incorporated in the dose calculation, and used for evaluation of cumulative dose and clinical dose-volume metrics after every dwell position, enabling real-time treatment verification. Main results. The dose calculated from source tracking measurements aligned with the generated uncertainty bands, validating the approach. Simulated shifts of 3 mm in 5/17 needles in a single plan caused DVH deviations beyond the uncertainty bands, indicating errors occurred during treatment. Clinical dose-volume metrics could be monitored in a time-resolved approach, enabling early detection of treatment plan deviations and prediction of their impact on the final dose that will be delivered in real-time. Significance. Integrating dose calculation with source tracking enhances the clinical relevance of IVD methods. Phantom measurements show that the developed tool aids in tracking treatment progress, detecting errors in real-time and post-treatment evaluation. In addition, it could be used to define patient-specific action limits and error thresholds, while taking the uncertainty of the measurement system into consideration.

Real-time motion-including dose estimation of simulated multi-leaf collimator-tracked magnetic resonance-guided radiotherapy.

Real-time dose estimation is a key-prerequisite to enable online intra-fraction treatment adaptation in magnetic resonance (MR)-guided radiotherapy (MRgRT). It is an essential component for the assessment of the dosimetric benefits and risks of online adaptive treatments, such as multi-leaf collimator (MLC)-tracking. We present a proof-of-concept for a software workflow for real-time dose estimation of MR-guided adaptive radiotherapy based on real-time data-streams of the linac delivery parameters and target positions. A software workflow, combining our in-house motion management software DynaTrack, a real-time dose calculation engine that connects to a research version of the treatment planning software (TPS) Monaco (v.6.09.00, Elekta AB, Stockholm, Sweden) was developed and evaluated. MR-guided treatment delivery on the Elekta Unity MR-linac was simulated with and without MLC-tracking for three prostate patients, previously treated on the Elekta Unity MR-linac (36.25Gy/five fractions). Three motion scenarios were used: no motion, regular motion, and erratic prostate motion. Accumulated monitor units (MUs), centre of mass target position and MLC-leaf positions, were forwarded from DynaTrack at a rate of 25Hz to a Monte Carlo (MC) based dose calculation engine which utilises the research GPUMCD-library (Elekta AB, Stockholm, Sweden). A rigid isocentre shift derived from the selected motion scenarios was applied to a bulk density-assigned session MR-image. The respective electron density used for treatment planning was accessed through the research Monaco TPS. The software workflow including the online dose reconstruction was validated against offline dose reconstructions. Our investigation showed that MC-based real-time dose calculations that account for all linac states (including MUs, MLC positions and target position) were infeasible, hence states were randomly sampled and used for calculation as follows; Once a new linac state was received, a dose calculation with 106 photons was started. Linac states that arrived during the time of the ongoing calculation were put into a queue. After completion of the ongoing calculation, one new linac state was randomly picked from the queue and assigned the MU accumulated from the previous state until the last sample in the queue. The queue was emptied, and the process repeated throughout treatment simulation. On average 27% (23%-30%) of received samples were used in the real-time calculation, corresponding to a calculation time for one linac state of 148ms. Median gamma pass rate (2%/3mm local) was 100.0% (99.9%-100%) within the PTV volume and 99.1% (90.1%-99.4.0%) with a 15% dose cut off. Differences in PTVDmean , CTVDmean , RectumD2% , and BladderD2% (offline-online, % of prescribed dose) were below 0.64%. Beam-by-beam comparisons showed deviations below 0.07Gy. Repeated simulations resulted in standard deviations below 0.31% and 0.12Gy for the investigated volume and dose criteria respectively. Real-time dose estimation was successfully performed using the developed software workflow for different prostate motion traces with and without MLC-tracking. Negligible dosimetric differences were seen when comparing online and offline reconstructed dose, enabling online intra-fraction treatment decisions based on estimates of the delivered dose.

Hydroxyurea Pharmacokinetics and Fetal Hemoglobin Response in Ugandan Children with Sickle Cell Anemia

Introduction Over 300,000 infants are born with sickle cell anemia (SCA) in Africa every year, and the under-5 mortality of this population is estimated to be 50-90% without treatment. Hydroxyurea is a safe, oral medication that decreases morbidity and mortality in SCA by increasing protective fetal hemoglobin (HbF) production and reducing inflammation. Hydroxyurea dosing to maximum tolerated dose (MTD) offers the optimal treatment with greatest laboratory responses and clinical benefits. Patients have variable hydroxyurea absorption and elimination rates, so the traditional approach to determining MTD is through incremental dose escalation to achieve mild myelosuppression and reduction in absolute neutrophil count (ANC). In the United States, pharmacokinetic (PK)-guided dosing of hydroxyurea for children with SCA has been shown to decrease time to MTD and to augment treatment responses, with most patients achieving over 30% HbF. Maximizing clinical benefits and minimizing logistical burdens of dosing hydroxyurea would be advantageous for areas with high SCA prevalence, but PK parameters of African children with SCA who might have higher rates of malnutrition, infections, and organ damage that affect drug absorption, distribution, and clearance are yet to be determined. Here we describe hydroxyurea PK exposure-response in Ugandan children with SCA from the Optimizing Hydroxyurea Therapy in Children with Sickle Cell Anemia In Malaria Endemic Areas: The NOHARM MTD Study (NOHARM MTD, NCT03128515). Methods NOHARM MTD was a randomized controlled trial that compared dosing strategies between hydroxyurea at fixed, weight-based dosing (~20 mg/kg/day) and dose escalation to MTD in children with SCA living in Uganda. NOHARM MTD participants underwent PK testing after 18 months of treatment. Children were administered their prescribed hydroxyurea dose in the morning at clinic, and blood was obtained by fingerstick using Mitra® microsampler devices after 20, 60, and 180 minutes. The microsampler devices were frozen at -80C in Uganda and then transported to Cincinnati Children's Hospital for analysis of hydroxyurea concentration via liquid chromatography with tandem mass spectrometry. Using our previously published population PK model, individual hydroxyurea exposures expressed as area-under-the-curve (AUC) were generated by Bayesian estimation using the PK modeling software NONMEM (Version 7.5, ICON Development Solutions, Ellicott City, USA). A hydroxyurea exposure-response analysis was completed with HbF% as the primary pharmacodynamic outcome. Results A total of 90 NOHARM MTD participants completed PK testing; samples from 12 were excluded from full analysis due to an elevated baseline concentration suggesting the participant already received the daily hydroxyurea dose before PK testing (9), supraphysiologic concentration (1), or mislabeling (2). PK profiles were analyzed for the remaining 78 participants (mean age 6.5 years, SD 0.9); of these, 41 were on fixed dose hydroxyurea and 37 were at MTD with mean doses of 18.7 mg/kg/day (SD 1.6) and 28.3 mg/kg/day (SD 4.3), respectively (Table). Those at MTD had a higher AUC with a significantly higher %HbF (mean 19% HbF on fixed dose vs 29% at MTD, p=0.0002) and lower ANC (mean 4.7 x109/L at fixed dose vs 3.2 at MTD, p=0.0045). Eighteen participants (23%) had malnutrition with a body mass index Z-score less than -1.0, and 12 (15%) participants had renal hyperfiltration with an estimated glomerular filtration rate over 150 mL/min/1.73m2; neither condition was associated with differences in AUC. Conclusion This study demonstrates the feasibility of collecting hydroxyurea PK data with timed, sparse sampling in Uganda, and the exposure-response analysis confirms that higher hydroxyurea exposure in Ugandan children with SCA is associated with greater induction of protective HbF and reduction in ANC. Malnutrition and glomerular hyperfiltration did not affect hydroxyurea exposure in this cohort. Further studies of PK-guided dosing to optimize response and time to response in this population are needed. Future directions include building capacity for local analysis of hydroxyurea concentrations using HPLC technology, real-time dose calculation via a PK model-informed algorithm, and prospective evaluation of PK-guided dosing. Figure 1View largeDownload PPTFigure 1View largeDownload PPT Close modal

Dosimetric outcomes of preoperative treatment planning with intraoperative optimization using stranded seeds in prostate brachytherapy.

This study aimed to evaluate the quality of low-dose-rate (LDR) prostate brachytherapy (BT) based on treatment-related dosimetric outcomes. Data of 100 patients treated using LDR BT with stranded seeds from November 2012 to November 2017 were collected. The prescription dose for the prostate was 145 Gy. The dose constraints for the preoperative plan were: V100% ≥ 95%, V150% ≤ 60%, V200% ≤ 20% for the prostate; V100% for rectum, ≤ 1 cc; and V200 Gy for urethra, 0.0 cc. Intraoperative real-time dose calculation and postoperative dose distribution analysis on days 0 and 30 were performed. Median dosimetric outcomes on days 0 and 30 respective were: V100% 92.28% and 92.23%, V200% 18.63% and 25.02%, and D90% 150.88 Gy and 151.46 Gy for the prostate; V100% for the rectum, 0.11 cc and 0.22 cc; and V200 Gy for the urethra, 0.00 cc and 0.00 cc, respectively. Twenty patients underwent additional seed implantation to compensate for insufficient dose coverage of the prostate. No loss or substantial migration of seeds or severe toxicity was reported. With stranded seed implantation and intraoperative optimization, appropriate dose delivery to the prostate without excessive dose to the organs at risk could be achieved.

XIORT-MC: A real-time MC-based dose computation tool for low- energy X-rays intraoperative radiation therapy.

The INTRABEAM system is a miniature accelerator for low-energy X-ray Intra-Operative Radiation Therapy (IORT), and it could benefit from a fast and accurate dose computation tool. With regards to accuracy, dose computed with Monte Carlo (MC) simulations are the gold standard, however, they require a large computational effort and consequently they are not suitable for real-time dose planning. This work presents a comparison of the implementation on Graphics Processing Unit (GPU) of two different dose calculation algorithms based on MC phase-space (PHSP) information to compute dose distributions for the INTRABEAM device within seconds and with the accuracy of realistic MC simulations. The MC-based algorithms we present incorporate photoelectric, Compton and Rayleigh effects for the interaction of low-energy X-rays. XIORT-MC (X-ray Intra-Operative Radiation Therapy Monte Carlo) includes two dose calculation algorithms; a Woodcock-based MC algorithm (WC-MC) and a Hybrid MC algorithm (HMC), and it is implemented in CPU and in GPU. Detailed MC simulations have been generated to validate our tool in homogeneous and heterogeneous conditions with all INTRABEAM applicators, including three clinically realistic CT-based simulations. A performance study has been done to determine the acceleration reached with the code, in both CPU and GPU implementations. Dose distributions were obtained with the HMC and the WC-MC and compared to standard reference MC simulations with more than 95% voxels fulfilling a 7%-0.5mm gamma evaluation in all the cases considered. The CPU-HMC is 100 times more efficient than the reference MC, and the CPU-WC-MC is about 50 times more efficient. With the GPU implementation, the particle tracking of the WC-MC is faster than the HMC, with the extraction of the particle's information from the PHSP file taking a major part of the time. However, thanks to the variance reduction techniques implemented in the HMC, up to 400 times less particles are needed in the HMC to reach the same level of noise than the WC-MC. Therefore, in our implementation for INTRABEAM energies, the HMC is about 1.3 times more efficient than the WC-MC in an NVIDIA GeForce GTX 1080 Ti card and about 5.5 times more efficient in an NVIDIA GeForce RTX 3090. Dose with noise below 5% has been obtained in realistic situations in less than 5 s with the WC-MC and in less than 0.5 s with the HMC. The XIORT-MC is a dose computation tool designed to take full advantage of modern GPUs, making possible to obtain MC-grade accurate dose distributions within seconds. Its high speed allows a real-time dose calculation that includes the realistic effects of the beam in voxelized geometries of patients. It can be used as a dose-planning tool in the operating room during a XIORT treatment with any INTRABEAM device.

Technical Note: A feasibility study on deep learning-based radiotherapy dose calculation.

Various dose calculation algorithms are available for radiation therapy for cancer patients. However, these algorithms are faced with the tradeoff between efficiency and accuracy. The fast algorithms are generally less accurate, while the accurate dose engines are often time consuming. In this work, we try to resolve this dilemma by exploring deep learning (DL) for dose calculation. We developed a new radiotherapy dose calculation engine based on a modified Hierarchically Densely Connected U-net (HD U-net) model and tested its feasibility with prostate intensity-modulated radiation therapy (IMRT) cases. Mapping from an IMRT fluence map domain to a three-dimensional (3D) dose domain requires a deep neural network of complicated architecture and a huge training dataset. To solve this problem, we first project the fluence maps to the dose domain using a broad beam ray-tracing (RT) algorithm, and then we use the HD U-net to map the RT dose distribution into an accurate dose distribution calculated using a collapsed cone convolution/superposition (CS) algorithm. The model is trained on 70 patients with fivefold cross validation, and tested on a separate 8 patients. It takes about 1s to compute a 3D dose distribution for a typical 7-field prostate IMRT plan, which can be further reduced to achieve real-time dose calculation by optimizing the network. The average Gamma passing rate between DL and CS dose distributions for the 8 test patients are 98.5% (±1.6%) at 1mm/1% and 99.9% (±0.1%) at 2mm/2%. For comparison of various clinical evaluation criteria (dose-volume points) for IMRT plans between two dose distributions, the average difference for dose criteria is less than 0.25Gy while for volume criteria is <0.16%, showing that the DL dose distributions are clinically identical to the CS dose distributions. We have shown the feasibility of using DL for calculating radiotherapy dose distribution with high accuracy and efficiency.

33 Entropic model for real-time dose calculation

Introduction This work proposes a completely new Grid-Based Boltzmann Solver (GBBS) developed for the transport and energy deposition by energetic particles. Its entropic closure and mathematical formulation allow our entropic model to calculate the delivered dose with an accuracy comparable to Monte Carlo (MC) codes with a computational time that is reduced to the order of few seconds without any special processing power requirement. Methods In contrast to discrete ordinates angular discretization methods, such as Acuros®, our method is based on a reduced number of moment equations closed with Boltzmann’s H-theorem. Keeping a good accuracy of calculations, the algorithm can simulate different treatment techniques such as the external radiotherapy even in presence of magnetic field [1] , [2] (e.g., MRI-guided radiotherapy), brachytherapy or intra-operative radiation therapy. The first validation step consists in simulating dose distributions in complex numerical phantoms including a large number of heterogeneity shapes such as bone, lung and air. For both brachytherapy and external beam radiotherapy, simulations based on CT scan, using the real phase-space of the source, have been performed. The entropic model is validated by a direct comparison with the reference MC code PENELOPE. Results The code is capable of calculating 3D dose distributions with 1 mm3voxels without statistical uncertainties in few seconds instead of several minutes like PENELOPE. In brachytherapy applications the dose distributions significantly differ from those calculated with the TG-43 approximations, thanks to its capability to account for inhomogeneities and strong density gradients. Moreover, for both applications the code shows an excellent agreement with PENELOPE within the 1%/1 mm gamma-index criterion. Conclusions In the comparison with the MC results the excellent accuracy of the model is demonstrated. Thanks to its reduced computational time and its accuracy, this model is a promising candidate to become a real-time dose calculation algorithm.

ECG-based 4D-dose reconstruction of cardiac arrhythmia ablation with carbon ion beams: application in a porcine model

Noninvasive ablation of cardiac arrhythmia by scanned particle radiotherapy is highly promising, but especially challenging due to cardiac and respiratory motion. Irradiations for catheter-free ablation in intact pigs were carried out at the GSI Helmholtz Center in Darmstadt using scanned carbon ions. Here, we present real-time electrocardiogram (ECG) data to estimate time-resolved (4D) delivered dose. For 11 animals, surface ECGs and temporal structure of beam delivery were acquired during irradiation. R waves were automatically detected from surface ECGs. Pre-treatment ECG-triggered 4D-CT phases were synchronized to the R–R interval. 4D-dose calculation was performed using GSI’s in-house 4D treatment planning system. Resulting dose distributions were assessed with respect to coverage (D95 and V95), heterogeneity (HI = D5–D95) and normal tissue exposure. Final results shown here were performed offline, but first calculations were started shortly after irradiation The D95 for TV and PTV was above 95% for 10 and 8 out of 11 animals, respectively. HI was reduced for PTV versus TV volumes, especially for some of the animals targeted at the atrioventricular junction, indicating residual interplay effects due to cardiac motion. Risk structure exposure was comparable to static and 4D treatment planning simulations. ECG-based 4D-dose reconstruction is technically feasible in a patient treatment-like setting. Further development of the presented approach, such as real-time dose calculation, may contribute to safe, successful treatments using scanned ion beams for cardiac arrhythmia ablation.

PO-0946: Entropic model for real-time dose calculation: I-125 prostate brachytherapy application

PO-0946: Entropic model for real-time dose calculation: I-125 prostate brachytherapy application

SU-F-303-17: Real Time Dose Calculation of MRI Guided Co-60 Radiotherapy Treatments On Free Breathing Patients, Using a Motion Model and Fast Monte Carlo Dose Calculation

Purpose: To demonstrate real-time dose calculation of free-breathing MRI guided Co−60 treatments, using a motion model and Monte-Carlo dose calculation to accurately account for the interplay between irregular breathing motion and an IMRT delivery. Methods: ViewRay Co-60 dose distributions were optimized on ITVs contoured from free-breathing CT images of lung cancer patients. Each treatment plan was separated into 0.25s segments, accounting for the MLC positions and beam angles at each time point. A voxel-specific motion model derived from multiple fast-helical free-breathing CTs and deformable registration was calculated for each patient. 3D images for every 0.25s of a simulated treatment were generated in real time, here using a bellows signal as a surrogate to accurately account for breathing irregularities. Monte-Carlo dose calculation was performed every 0.25s of the treatment, with the number of histories in each calculation scaled to give an overall 1% statistical uncertainty. Each dose calculation was deformed back to the reference image using the motion model and accumulated. The static and real-time dose calculations were compared. Results: Image generation was performed in real time at 4 frames per second (GPU). Monte-Carlo dose calculation was performed at approximately 1frame per second (CPU), giving a total calculation time of approximately 30 minutes per treatment. Results show both cold- and hot-spots in and around the ITV, and increased dose to contralateral lung as the tumor moves in and out of the beam during treatment. Conclusion: An accurate motion model combined with a fast Monte-Carlo dose calculation allows almost real-time dose calculation of a free-breathing treatment. When combined with sagittal 2D-cine-mode MRI during treatment to update the motion model in real time, this will allow the true delivered dose of a treatment to be calculated, providing a useful tool for adaptive planning and assessing the effectiveness of gated treatments.

TH-F-211-04: A Fast Finite Size Pencil Beam Algorithm for Dose Calculation Using GPUs

Purpose: To facilitate real-time dose calculation for online adaptive radiation therapy, we have implemented the Monte Carlo-based finite size pencil beam (FSPB) dose calculation algorithm on Graphics Processing Units (GPUs). Methods: The FSPB dose calculation algorithm is based on the analytical kernel calculated from Monte Carlo simulations using the EGSnrc Monte Carlo code. This provides the benefit of Monte Carlo accuracy while the GPU implementation and the exponential kernel significantly reduces execution time. The GPU-based FSPB implementation was executed on the NVidia GTX480 GPU card. Both CPU and GPU versions of the FSPB algorithm were implemented operationally identically and benchmarked for accuracy and execution time. The dose calculation was performed in heterogeneous phantoms consisting of water, lung and bone slabs, and in clinical cases. The GPU-based FSPB implementation was compared with the CPU version and with a commercial treatment planning system for different situations including IMRT plans. Results: The FSPB algorithm with heterogeneity correction and anisotropic analytical algorithm (AAA) in heterogeneous phantom and clinical treatment planning cases agree within 2%. The GPU implementation was found to be up to 477 times faster than the single threaded CPU while achieving identical result in terms of dosimetric accuracy. The GPU-based FSPB was close to six times faster than the AAA, with single beam calculations completed within 0.45 s for typical scenarios. Conclusions: This work describes the implementation of a GPU-based FSPB algorithm for dose calculation in radiation therapy. The results demonstrated that the GPU-based FSPB algorithm can provide accurate 3D dose distributions in heterogeneous media in less than a second per beam for clinically relevant cases. This research is supported by CPRIT Individual Investigator Award RP110329.

Verification and Validation of GPU-based TomoTherapy Dose Calculation Engine

Dose calculation is essential in radiation therapy treatment planning. TomoTherapy uses accurate Collapsed Cone Convolution/superposition (CCCS) algorithm as its dose calculation engine. As the CCCS calculation is computationally demanding, computer cluster with 28 cpu cores or 56 cpu cores is provided for the optimization and dose calculation. The recent generation of graphic processing units (GPUs) possesses great computing power due to its massively parallel architecture. We have developed an ultrafast GPU-based Tomotherapy® dose calculation engine capable of performing near real-time dose calculation with a single PC.

Accelerated ray tracing for radiotherapy dose calculations on a GPU

The graphical processing unit (GPU) on modern graphics cards offers the possibility of accelerating arithmetically intensive tasks. By splitting the work into a large number of independent jobs, order-of-magnitude speedups are reported. In this article, the possible speedup of PLATO's ray tracing algorithm for dose calculations using a GPU is investigated. A GPU version of the ray tracing algorithm was implemented using NVIDIA's CUDA, which extends the standard C language with functionality to program graphics cards. The developed algorithm was compared based on the accuracy and speed to a multithreaded version of the PLATO ray tracing algorithm. This comparison was performed for three test geometries, a phantom and two radiotherapy planning CT datasets (a pelvic and a head-and-neck case). For each geometry, four different source positions were evaluated. In addition to this, for the head-and-neck case also a vertex field was evaluated. The GPU algorithm was proven to be more accurate than the PLATO algorithm by elimination of the look-up table for z indices that introduces discretization errors in the reference algorithm. Speedups for ray tracing were found to be in the range of 2.1-10.1, relative to the multithreaded PLATO algorithm running four threads. For dose calculations the speedup measured was in the range of 1.5-6.2. For the speedup of both the ray tracing and the dose calculation, a strong dependency on the tested geometry was found. This dependency is related to the fraction of air within the patient's bounding box resulting in idle threads. With the use of a GPU, ray tracing for dose calculations can be performed accurately in considerably less time. Ray tracing was accelerated, on average, with a factor of 6 for the evaluated cases. Dose calculation for a single beam can typically be carried out in 0.6-0.9 s for clinically realistic datasets. These findings can be used in conventional planning to enable (nearly) real-time dose calculations. Also the importance for treatment optimization techniques is evident.

Needle position during 125I seed implantation: accurately recognized by sagittal transrectal ultrasonography

The aim of this study was to assess the variation of probe rotation angles for detecting a single needle using sagittal images of transrectal ultrasonography (TRUS). A phantom study was performed. One needle was inserted through each of 10 holes of the template, and variations in the probe rotation angles for detecting the needle were measured. The mean variation of probe rotation for detecting a single needle was 17.0 degrees (range 4 degrees -25 degrees ). Slightly broader variation was seen for the needle in holes farther away from the probe. Probe rotation angles for detecting a single needle displayed considerable variation. Seed locations recognized on sagittal imaging by TRUS are thus indeterminate, and real-time dose calculations using TRUS for (125)I seed implantation should be used with care.

Real-time dose calculation and visualization for the proton therapyof ocular tumours

A new real-time dose calculation and visualization was developed as part ofthe new 3D treatment planning tool OCTOPUS for proton therapy of oculartumours within a national research project together with the Hahn-MeitnerInstitut Berlin. The implementation resolves the common separation betweenparameter definition, dose calculation and evaluation and allows a directexamination of the expected dose distribution while adjusting the treatmentparameters. The new tool allows the therapist to move the desired dosedistribution under visual control in 3D to the appropriate place. Thevisualization of the resulting dose distribution as a 3D surface model, on any2D slice or on the surface of specified ocular structures is doneautomatically when adapting parameters during the planning process. Inaddition, approximate dose volume histograms may be calculated with littleextra time. The dose distribution is calculated and visualized in 200 ms withan accuracy of 6% for the 3D isodose surfaces and 8% for other objects. Thispaper discusses the advantages and limitations of this new approach.

Real-time 3D dose calculation and display: A tool for plan optimization

Purpose : Both human and computer optimization of treatment plans have advantages; humans are much better at global pattern recognition, and computers are much better at detailed calculations. A major impediment to human optimization of treatment plans by manipulation of beam parameters is the long time required for feedback to the operator on the effectiveness of a change in beam parameters. Our goal was to create a real-time dose calculation and display system that provides the planner with immediate (fraction of a second) feedback with displays of three-dimensional (3D) isodose surfaces, digitally reconstructed radiographs (DRRs), dose-volume histograms, and/or a figure of merit (FOM) (i.e., a single value plan score function). This will allow the experienced treatment planner to optimize a plan by adjusting beam parameters based on a direct indication of plan effectiveness, the FOM value, and to use 3D display of target, critical organs, DRRs, and isodose contours to guide changes aimed at improving the FOM value. Methods and Materials : We use computer platforms that contain easily utilized parallel processors and very tight coupling between calculation and display. We ported code running on a network of two workstations and an array of transputers to a single multiprocessor workstation. Our current high-performance graphics workstation contains four 150-MHz processors that can be readily used in a shared-memory multithreaded calculation. Results : When a 10 × 10-cm beam is moved, using an 8-mm dose grid, the full 3D dose matrix is recalculated using a Bentley—Milan-type dose calculation algorithm, and the 3D dose surface display is then updated, all in < 0.1 s. A 64 × 64-pixel DRR calculation can be performed in < 0.1 s. Other features, such as automated aperture calculation, are still required to make real-time feedback practical for clinical use. Conclusion : We demonstrate that real-time plan optimization using general pyrpose multiprocessor workstations is a practical goal. Parallel processing technology provides this capability for 3D planning systems, and when combined with objective plan ranking algorithms should prove effective for optimizing 3D conformal radiation therapy. Compared to our earlier transputer work, multiprocessor workstations are more easily programmed, making software development costs more reasonable compared with uniprocessor development costs. How the dose calculation is partitioned into parallel tasks on a multiprocessor work station can make a significant difference in performance. Shared-memory multiprocessor workstations are our first choice for future work, because they require minimum programming effort and continue to be driven to higher performance by competition in the workstation arena.