Search

The purpose of this work is to demonstrate an ultra-fast reconstruction technique for digital tomosynthesis (DTS) imaging based on the algorithm proposed by Feldkamp, Davis, and Kress (FDK) using standard general-purpose graphics processing unit (GPGPU) programming interface. To this end, the FDK-based DTS algorithm was programmed “in-house” with C language with utilization of 1) GPU and 2) central processing unit (CPU) cards. The GPU card consisted of 480 processing cores (2 × 240 dual chip) with 1,242 MHz processing clock speed and 1,792 MB memory space. In terms of CPU hardware, we used 2.68 GHz clock speed, 12.0 GB DDR3 RAM, on a 64-bit OS. The performance of proposed algorithm was tested on twenty-five patient cases (5 lung, 5 liver, 10 prostate, and 5 head-and-neck) scanned either with a full-fanor half-fanmode on our cone-beam computed tomography (CBCT) system. For the full-fanscans, the projections from 157.5°–202.5° (45°-scan) were used to reconstruct coronal DTS slices, whereas for the half-fanscans, the projections from both 157.5°–202.5° and 337.5°–22.5° (2 × 45°-scan) were used to reconstruct larger FOV coronal DTS slices. For this study, we chose 45°-scan angle that contained ~80 projections for the full-fanand ~160 projections with 2 × 45°-scan angle for the half-fanmode, each with 1024 × 768 pixels with 32-bit precision. Absolute pixel value differences, profiles, and contrast-to-noise ratio (CNR) calculations were performed to compare and evaluate the images reconstructed using GPU- and CPU-based implementations. The time dependence on the reconstruction volume was also tested with (512 × 512) × 16, 32, 64, 128, and 256 slices. In the end, the GPU-based implementation achieved, at most, 1.3 and 2.5 seconds to complete full reconstruction of 512 × 512 × 256 volume, for the full-fan and half-fan modes, respectively. In turn, this meant that our implementation can process > 13 projections-per-second (pps) and > 18 pps for the full-fanand half-fanmodes, respectively. Since commercial CBCT system nominally acquires 11 pps (with 1 gantry-revolution-per-minute), our GPU-based implementation is sufficient to handle the incoming projections data as they are acquired and reconstruct the entire volume immediately after completing the scan. In addition, on increasing the number of slices (hence volume) to be reconstructed from 16 to 256, only minimal increases in reconstruction time were observed for the GPU-based implementation where from 0.73 to 1.27 seconds and 1.42 to 2.47 seconds increase were observed for the full-fanand half-fanmodes, respectively. This resulted in speed improvement of up to 87 times compared with the CPU-based implementation (for 256 slices case), with visually identical images and small pixel-value discrepancies (< 6.3%), and CNR differences (< 2.3%). With this achievement, we have shown that time allocation for DTS image reconstruction is virtually eliminated and that clinical implementation of this approach has become quite appealing. In addition, with the speed achievement, further image processing and real-time applications that was prohibited prior due to time restrictions can now be tempered with.

The purpose of this study is to accurately estimate the difference between the planned and the delivered dose due to respiratory motion and free breathing helical CT artefacts for lung IMRT treatments, and to estimate the impact of this difference on clinical outcome. Six patients with representative tumour motion, size and position were selected for this retrospective study. For each patient, we had acquired both a free breathing helical CT and a ten-phase 4D-CT scan. A commercial treatment planning system was used to create four IMRT plans for each patient. The first two plans were based on the GTV as contoured on the free breathing helical CT set, with a GTV to PTV expansion of 1.5 cm and 2.0 cm, respectively. The third plan was based on the ITV, a composite volume formed by the union of the CTV volumes contoured on free breathing helical CT, end-of-inhale (EOI) and end-of-exhale (EOE) 4D-CT. The fourth plan was based on GTV contoured on the EOE 4D-CT. The prescribed dose was 60 Gy for all four plans. Fluence maps and beam setup parameters of the IMRT plans were used by the Monte Carlo dose calculation engine MCSIM for absolute dose calculation on both the free breathing CT and 4D-CT data. CT deformable registration between the breathing phases was performed to estimate the motion trajectory for both the tumour and healthy tissue. Then, a composite dose distribution over the whole breathing cycle was calculated as a final estimate of the delivered dose. EUD values were computed on the basis of the composite dose for all four plans. For the patient with the largest motion effect, the difference in the EUD of CTV between the planed and the delivered doses was 33, 11, 1 and 0 Gy for the first, second, third and fourth plan, respectively. The number of breathing phases required for accurate dose prediction was also investigated. With the advent of 4D-CT, deformable registration and Monte Carlo simulations, it is feasible to perform an accurate calculation of the delivered dose, and compare our delivered dose with doses estimated using prior techniques.

Image guidance in radiotherapy and extracranial radiosurgery offers the potential for precise radiation dose delivery to a moving tumour. Recent work has demonstrated how to locate and track the position of a tumour in real-time using diagnostic x-ray imaging to find implanted radio-opaque markers. However, the delivery of a treatment plan through gating or beam tracking requires adequate consideration of treatment system latencies, including image acquisition, image processing, communication delays, control system processing, inductance within the motor, mechanical damping, etc. Furthermore, the imaging dose given over long radiosurgery procedures or multiple radiotherapy fractions may not be insignificant, which means that we must reduce the sampling rate of the imaging system. This study evaluates various predictive models for reducing tumour localization errors when a real-time tumour-tracking system targets a moving tumour at a slow imaging rate and with large system latencies. We consider 14 lung tumour cases where the peak-to-peak motion is greater than 8 mm, and compare the localization error using linear prediction, neural network prediction and Kalman filtering, against a system which uses no prediction. To evaluate prediction accuracy for use in beam tracking, we compute the root mean squared error between predicted and actual 3D motion. We found that by using prediction, root mean squared error is improved for all latencies and all imaging rates evaluated. To evaluate prediction accuracy for use in gated treatment, we present a new metric that compares a gating control signal based on predicted motion against the best possible gating control signal. We found that using prediction improves gated treatment accuracy for systems that have latencies of 200 ms or greater, and for systems that have imaging rates of 10 Hz or slower.

The design of an integrated radiotherapy imaging system (IRIS), consisting of gantry mounted diagnostic (kV) x-ray tubes and fast read-out flat-panel amorphous-silicon detectors, has been studied. The system is meant to be capable of three main functions: radiographs for three-dimensional (3D) patient set-up, cone-beam CT and real-time tumour/marker tracking. The goal of the current study is to determine whether one source/panel pair is sufficient for real-time tumour/marker tracking and, if two are needed, the optimal position of each relative to other components and the isocentre. A single gantry-mounted source/imager pair is certainly capable of the first two of the three functions listed above and may also be useful for the third, if combined with prior knowledge of the target's trajectory. This would be necessary because only motion in two dimensions is visible with a single imager/source system. However, with previously collected information about the trajectory, the third coordinate may be derived from the other two with sufficient accuracy to facilitate tracking. This deduction of the third coordinate can only be made if the 3D tumour/marker trajectory is consistent from fraction to fraction. The feasibility of tumour tracking with one source/imager pair has been theoretically examined here using measured lung marker trajectory data for seven patients from multiple treatment fractions. The patients' selection criteria include minimum mean amplitudes of the tumour motions greater than 1 cm peak-to-peak. The marker trajectory for each patient was modelled using the first fraction data. Then for the rest of the data, marker positions were derived from the imager projections at various gantry angles and compared with the measured tumour positions. Our results show that, due to the three dimensionality and irregular trajectory characteristics of tumour motion, on a fraction-to-fraction basis, a ‘monoscopic’ system (single source/imager) is inadequate for consistent real-time tumour tracking, even with prior knowledge. We found that, among the seven patients studied with peak-to-peak marker motion greater than 1 cm, five cases have mean localization errors greater than 2 mm and two have mean errors greater than 3 mm. Because of this uncertainty associated with a monoscopic system, two source/imager pairs are necessary for robust 3D target localization. Dual orthogonal x-ray source/imager pairs mounted on the linac gantry are chosen for the IRIS. We further studied the placement of the x-ray sources/panel based on the geometric specifications of the Varian 21EX Clinac. The best configuration minimizes the localization error while maintaining a large field of view and avoiding collisions with the floor/ceiling or couch.

Respiration-induced tumour motion can potentially compromise the use of intensity-modulated radiotherapy (IMRT) as a dose escalation tool for lung tumour treatment. We have experimentally investigated the intra-fractional organ motion effects in lung IMRT treatments delivered by multi-leaf collimator (MLC). An in-house made motor-driven platform, which moves sinusoidally with an amplitude of 1 cm and a period of 4 s, was used to mimic tumour motion. Tumour motion was simulated along cranial-caudal direction while MLC leaves moved across the patient from left to right, as in most clinical cases. The dose to a point near the centre of the tumour mass was measured according to geometric and dosimetric parameters from two five-field lung IMRT plans. For each field, measurement was done for two dose rates (300 and 500 MU min⁻¹), three MLC delivery modes (sliding window, step-and-shoot with 10 and 20 intensity levels) and eight equally spaced starting phases of tumour motion. The dose to the measurement point delivered from all five fields was derived for both a single fraction and 30 fractions by randomly sampling from measured dose values of each field at different initial phases. It was found that the mean dose to a moving tumour differs slightly (&lt;2–3&percent;) from that to a static tumour. The variation in breathing phase at the start of dose delivery results in a maximum variation around the mean dose of greater than 30&percent; for one field. The full width at half maximum for the probability distribution of the point dose is up to 8&percent; for all five fields in a single fraction, but less than 1–2&percent; after 30 fractions. In general, lower dose rate can reduce the motion-caused dose variation and therefore might be preferable for lung IMRT when no motion mitigation techniques are used. From the two IMRT cases we studied where tumour motion is perpendicular to MLC leaf motion, the dose variation was found to be insensitive to the MLC delivery mode.

How to speed up Monte Carlo (MC) simulation in dose calculation without losing its intrinsic accuracy is one of the key issues of making a clinical MC dose engine. In this study we intensively investigated a special parallel computation technique, the vectorization technique, to boost simulation efficiency on a personal computer (PC) without extra hardware investment. A MC code, dose planning method (DPM), was extensively modified into a vectorized code, V-DPM, using the streaming single-instruction–multiple-data extension (SSE) parallel computation model. Comparative simulations were conducted for typical simulation cases in both DPM and V-DPM codes. We found that in every case the V-DPM code runs 1.5 times faster than the DPM code with variance of 0.6&percent;.

Synchronized moving aperture radiation therapy (SMART) is a new technique for treating mobile tumours under development at Massachusetts General Hospital (MGH). The basic idea of SMART is to synchronize the moving radiation beam aperture formed by a dynamic multileaf collimator (DMLC) with the tumour motion induced by respiration. SMART is based on the concept of the average tumour trajectory (ATT) exhibited by a tumour during respiration. During the treatment simulation stage, tumour motion is measured and the ATT is derived. Then, the original IMRT MLC leaf sequence is modified using the ATT to compensate for tumour motion. During treatment, the tumour motion is monitored. The treatment starts when leaf motion and tumour motion are synchronized at a specific breathing phase. The treatment will halt when the tumour drifts away from the ATT and will resume when the synchronization between tumour motion and radiation beam is re-established. In this paper, we present a method to derive the ATT from measured tumour trajectory data. We also investigate the validity of the ATT concept for lung tumours during normal breathing. The lung tumour trajectory data were acquired during actual radiotherapy sessions using a real-time tumour-tracking system. SMART treatment is simulated by assuming that the radiation beam follows the derived ATT and the tumour follows the measured trajectory. In simulation, the treatment starts at exhale phase. The duty cycle of SMART delivery was calculated for various treatment times and gating thresholds, as well as for various exhale phases where the treatment begins. The simulation results show that in the case of free breathing, for 4 out of 11 lung datasets with tumour motion greater than 1 cm from peak to peak, the error in tumour tracking can be controlled to within a couple of millimetres while maintaining a reasonable delivery efficiency. That is to say, without any breath coaching/control, the ATT is a valid concept for some lung tumours. However, to make SMART an efficient technique in general, it is found that breath coaching techniques are required.

Total body irradiation is an important part of the conditioning regimens frequently used to prepare patients for allogeneic hematopoietic stem cell transplantation (SCT). Volumetric-modulated arc therapy enabled total body irradiation (VMAT-TBI), an alternative to conventional TBI (cTBI), is a novel radiotherapy treatment technique that has been implemented and investigated in our institution. The purpose of this study is to (1) report our six-year clinical experience in terms of treatment planning strategy and delivery time and (2) evaluate the clinical outcomes and toxicities in our cohort of patients treated with VMAT-TBI. This is a retrospective single center study. Forty-four patients at our institution received VMAT-TBI and chemotherapy conditioning followed by allogeneic SCT between 2014 and 2020. Thirty-two patients (73%) received standard-dose TBI (12-13.2 Gy in 6-8 fractions twice daily), whereas 12 (27%) received low-dose TBI (2-4 Gy in one fraction). Treatment planning, delivery, and treatment outcome data including overall survival (OS), relapse-free survival (RFS), and toxicities were analyzed. The developed VMAT-TBI planning strategy consistently generated plans satisfying our dose constraints, with planning target volume coverage >90%, mean lung dose ∼50% to 75% of prescription dose, and minimal hotspots in critical organs. Most of the treatment deliveries were <100 minutes (range 33-147, mean 72). The median follow-up was 26 months. At the last follow-up, 34 of 44 (77%) of patients were alive, with 1- and 2-year OS of 90% and 79% and RFS of 88% and 71%, respectively. The most common grade 3+ toxicities observed were mucositis (31 patients [71%]) and nephrotoxicity (6 patients [13%]), both of which were deemed multifactorial in cause. Four patients (9%) in standard-dose cohort developed grade 3+ pneumonitis, with 3 cases in the setting of documented respiratory infection and only 1 (2%) deemed likely related to radiation alone. VMAT-TBI provides a safe alternative to cTBI. The dose modulation capability of VMAT-TBI may lead to new treatment strategies, such as simultaneous boost and further critical organ sparing, for better malignant cell eradication, immune suppression, and lower toxicities.