W Mao

Combined kV and MV imaging for real-time tracking of implanted fiducial markers

In the presence of intrafraction organ motion, target localization uncertainty can greatly hamper the advantage of highly conformal dose techniques such as intensity modulated radiation therapy (IMRT). To minimize the adverse dosimetric effect caused by tumor motion, a real-time knowledge of the tumor position is required throughout the beam delivery process. The recent integration of onboard kV diagnostic imaging together with MV electronic portal imaging devices on linear accelerators can allow for real-time three-dimensional (3D) tumor position monitoring during a treatment delivery. The aim of this study is to demonstrate a near real-time 3D internal fiducial tracking system based on the combined use of kV and MV imaging. A commercially available radiotherapy system equipped with both kV and MV imaging systems was used in this work. A hardware video frame grabber was used to capture both kV and MV video streams simultaneously through independent video channels at 30 frames per second. The fiducial locations were extracted from the kV and MV images using a software tool. The geometric tracking capabilities of the system were evaluated using a pelvic phantom with embedded fiducials placed on a moveable stage. The maximum tracking speed of the kV/MV system is approximately 9 Hz, which is primarily limited by the frame rate of the MV imager. The geometric accuracy of the system is found to be on the order of less than 1 mm in all three spatial dimensions. The technique requires minimal hardware modification and is potentially useful for image-guided radiation therapy systems.

Fast internal marker tracking algorithm for onboard MV and kV imaging systems

Intrafraction organ motion can limit the advantage of highly conformal dose techniques such as intensity modulated radiation therapy (IMRT) due to target position uncertainty. To ensure high accuracy in beam targeting, real-time knowledge of the target location is highly desired throughout the beam delivery process. This knowledge can be gained through imaging of internally implanted radio-opaque markers with fluoroscopic or electronic portal imaging devices (EPID). In the case of MV based images, marker detection can be problematic due to the significantly lower contrast between different materials in comparison to their kV-based counterparts. This work presents a fully automated algorithm capable of detecting implanted metallic markers in both kV and MV images with high consistency. Using prior CT information, the algorithm predefines the volumetric search space without manual region-of-interest (ROI) selection by the user. Depending on the template selected, both spherical and cylindrical markers can be detected. Multiple markers can be simultaneously tracked without indexing confusion. Phantom studies show detection success rates of 100% for both kV and MV image data. In addition, application of the algorithm to real patient image data results in successful detection of all implanted markers for MV images. Near real-time operational speeds of approximately 10 frames/sec for the detection of five markers in a 1024 x 768 image are accomplished using an ordinary PC workstation.

Predicting respiratory tumor motion with multi-dimensional adaptive filters and support vector regression

Intra-fraction tumor tracking methods can improve radiation delivery during radiotherapy sessions. Image acquisition for tumor tracking and subsequent adjustment of the treatment beam with gating or beam tracking introduces time latency and necessitates predicting the future position of the tumor. This study evaluates the use of multi-dimensional linear adaptive filters and support vector regression to predict the motion of lung tumors tracked at 30 Hz. We expand on the prior work of other groups who have looked at adaptive filters by using a general framework of a multiple-input single-output (MISO) adaptive system that uses multiple correlated signals to predict the motion of a tumor. We compare the performance of these two novel methods to conventional methods like linear regression and single-input, single-output adaptive filters. At 400 ms latency the average root-mean-square-errors (RMSEs) for the 14 treatment sessions studied using no prediction, linear regression, single-output adaptive filter, MISO and support vector regression are 2.58, 1.60, 1.58, 1.71 and 1.26 mm, respectively. At 1 s, the RMSEs are 4.40, 2.61, 3.34, 2.66 and 1.93 mm, respectively. We find that support vector regression most accurately predicts the future tumor position of the methods studied and can provide a RMSE of less than 2 mm at 1 s latency. Also, a multi-dimensional adaptive filter framework provides improved performance over single-dimension adaptive filters. Work is underway to combine these two frameworks to improve performance.

Real-time 3D internal marker tracking during arc radiotherapy by the use of combined MV?kV imaging

To minimize the adverse dosimetric effect caused by tumor motion, it is desirable to have real-time knowledge of the tumor position throughout the beam delivery process. A promising technique to realize the real-time image guided scheme in external beam radiation therapy is through the combined use of MV and onboard kV beam imaging. The success of this MV-kV triangulation approach for fixed-gantry radiation therapy has been demonstrated. With the increasing acceptance of modern arc radiotherapy in the clinics, a timely and clinically important question is whether the image guidance strategy can be extended to arc therapy to provide the urgently needed real-time tumor motion information. While conceptually feasible, there are a number of theoretical and practical issues specific to the arc delivery that need to be resolved before clinical implementation. The purpose of this work is to establish a robust procedure of system calibration for combined MV and kV imaging for internal marker tracking during arc delivery and to demonstrate the feasibility and accuracy of the technique. A commercially available LINAC equipped with an onboard kV imager and electronic portal imaging device (EPID) was used for the study. A custom built phantom with multiple ball bearings was used to calibrate the stereoscopic MV-kV imaging system to provide the transformation parameters from imaging pixels to 3D world coordinates. The accuracy of the fiducial tracking system was examined using a 4D motion phantom capable of moving in accordance with a pre-programmed trajectory. Overall, spatial accuracy of MV-kV fiducial tracking during the arc delivery process for normal adult breathing amplitude and period was found to be better than 1 mm. For fast motion, the results depended on the imaging frame rates. The RMS error ranged from approximately 0.5 mm for the normal adult breathing pattern to approximately 1.5 mm for more extreme cases with a low imaging frame rate of 3.4 Hz. In general, highly accurate real-time tracking of implanted markers using hybrid MV-kV imaging is achievable and the technique should be useful to improve the beam targeting accuracy of arc therapy.

Development of a QA phantom and automated analysis tool for geometric quality assurance of on-board MV and kV x-ray imaging systems.

The medical linear accelerator (linac) integrated with a kilovoltage (kV) flat-panel imager has been emerging as an important piece of equipment for image-guided radiation therapy. Due to the sagging of the linac head and the flexing of the robotic arms that mount the x-ray tube and flat-panel detector, geometric nonidealities generally exist in the imaging geometry no matter whether it is for the two-dimensional projection image or three-dimensional cone-beam computed tomography. Normally, the geometric parameters are established during the commissioning and incorporated in correction software in respective image formation or reconstruction. A prudent use of an on-board imaging system necessitates a routine surveillance of the geometric accuracy of the system like the position of the x-ray source, imager position and orientation, isocenter, rotation trajectory, and source-to-imager distance. Here we describe a purposely built phantom and a data analysis software for monitoring these important parameters of the system in an efficient and automated way. The developed tool works equally well for the megavoltage (MV) electronic portal imaging device and hence allows us to measure the coincidence of the isocenters of the MV and kV beams of the linac. This QA tool can detect an angular uncertainty of 0.1 degrees of the x-ray source. For spatial uncertainties, such as the source position, the imager position, or the kV/MV isocenter misalignment, the demonstrated accuracy of this tool was better than 1.6 mm. The developed tool provides us with a simple, robust, and objective way to probe and monitor the geometric status of an imaging system in a fully automatic process and facilitate routine QA workflow in a clinic.

A fiducial detection algorithm for real-time image guided IMRT based on simultaneous MV and kV imaging

The advantage of highly conformal dose techniques such as 3DCRT and IMRT is limited by intrafraction organ motion. A new approach to gain near real-time 3D positions of internally implanted fiducial markers is to analyze simultaneous onboard kV beam and treatment MV beam images (from fluoroscopic or electronic portal image devices). Before we can use this real-time image guidance for clinical 3DCRT and IMRT treatments, four outstanding issues need to be addressed. (1) How will fiducial motion blur the image and hinder tracking fiducials? kV and MV images are acquired while the tumor is moving at various speeds. We find that a fiducial can be successfully detected at a maximum linear speed of 1.6 cm/s. (2) How does MV beam scattering affect kV imaging? We investigate this by varying MV field size and kV source to imager distance, and find that common treatment MV beams do not hinder fiducial detection in simultaneous kV images. (3) How can one detect fiducials on images from 3DCRT and IMRT treatment beams when the MV fields are modified by a multileaf collimator (MLC)? The presented analysis is capable of segmenting a MV field from the blocking MLC and detecting visible fiducials. This enables the calculation of nearly real-time 3D positions of markers during a real treatment. (4) Is the analysis fast enough to track fiducials in nearly real time? Multiple methods are adopted to predict marker positions and reduce search regions. The average detection time per frame for three markers in a 1024 x 768 image was reduced to 0.1 s or less. Solving these four issues paves the way to tracking moving fiducial markers throughout a 3DCRT or IMRT treatment. Altogether, these four studies demonstrate that our algorithm can track fiducials in real time, on degraded kV images (MV scatter), in rapidly moving tumors (fiducial blurring), and even provide useful information in the case when some fiducials are blocked from view by the MLC. This technique can provide a gating signal or be used for intra-fractional tumor tracking on a Linac equipped with a kV imaging system. Any motion exceeding a preset threshold can warn the therapist to suspend a treatment session and reposition the patient.

Scatter correction for cone-beam computed tomography using moving blocker strips: A preliminary study

Purpose: One well-recognized challenge of cone-beam computed tomography(CBCT) is the presence of scatter contamination within the projection images.Scatter degrades the CBCTimage quality by decreasing the contrast, introducing shading artifacts, and leading to inaccuracies in the reconstructedCT number. The authors propose a blocker-based approach to simultaneously estimate scatter signal and reconstruct the complete volume within the field of view (FOV) from a single CBCT scan. Methods: A physical strip attenuator (i.e., “blocker”), consisting of lead strips, is inserted between the x-ray source and the patient. The blocker moves back and forth along the z -axis during the gantry rotation. With such a design, the data required for the filtering step of the Feldkamp–Davis–Kress (FDK) algorithm are complete in the unblocked region and the entire volume within the FOV has the measurements at different projection views. The two-dimensional scatter fluence is estimated by interpolating the signal from the blocked regions. A modified FDK algorithm and an iterative reconstruction based on the constraint optimization are used to reconstructCBCTimages from unblocked projection data after the scatter signal is subtracted. A simulation study and an experimental study are performed to evaluate the performance of the proposed scatter correction scheme. Results: The scatter-induced shading/cupping artifacts are substantially reduced in CBCT using the proposed strategy. In the simulation study, the mean relative error is reduced from 25% to 3% and 2% in the imagesreconstructed by the modified FDK and constraint optimization, respectively. In the experimental study using a CatPhan® 600 phantom, CT number errors in the selected regions of interest are reduced from 256 to less than 20. Conclusions: An effective scatter correction scheme is proposed for CBCT. A moving blocker consisting of lead strips is inserted between the x-ray source and the patient during CBCT acquisition. The proposed method allows the authors to simultaneously estimate the scatter signal in projection data, reduce the imaging dose, and obtain complete volumetric information within the FOV.

A novel markerless technique to evaluate daily lung tumor motion based on conventional cone-beam CT projection data.

PURPOSE In this study, we present a novel markerless technique, based on cone beam computed tomography (CBCT) raw projection data, to evaluate lung tumor daily motion. METHOD AND MATERIALS The markerless technique, which uses raw CBCT projection data and locates tumors directly on every projection, consists of three steps. First, the tumor contour on the planning CT is used to create digitally reconstructed radiographs (DRRs) at every projection angle. Two sets of DRRs are created: one showing only the tumor, and another with the complete anatomy without the tumor. Second, a rigid two-dimensional image registration is performed to register the DRR set without the tumor to the CBCT projections. After the registration, the projections are subtracted from the DRRs, resulting in a projection dataset containing primarily tumor. Finally, a second registration is performed between the subtracted projection and tumor-only DRR. The methodology was evaluated using a chest phantom containing a moving tumor, and retrospectively in 4 lung cancer patients treated by stereotactic body radiation therapy. Tumors detected on projection images were compared with those from three-dimensional (3D) and four-dimensional (4D) CBCT reconstruction results. RESULTS Results in both static and moving phantoms demonstrate that the accuracy is within 1 mm. The subsequent application to 22 sets of CBCT scan raw projection data of 4 lung cancer patients includes about 11,000 projections, with the detected tumor locations consistent with 3D and 4D CBCT reconstruction results. This technique reveals detailed lung tumor motion and provides additional information than conventional 4D images. CONCLUSION This technique is capable of accurately characterizing lung tumor motion on a daily basis based on a conventional CBCT scan. It provides daily verification of the tumor motion to ensure that these motions are within prior estimation and covered by the treatment planning volume.

The use of EPID‐measured leaf sequence files for IMRT dose reconstruction in adaptive radiation therapy

For intensity modulated radiation treatment (IMRT) dose reconstruction, multileaf collimator (MLC) log files have been shown applicable for deriving delivered fluence maps. However, MLC log files are dependent on the accuracy of leaf calibration and only available from one linear accelerator manufacturer. This paper presents a proof of feasibility and principles in (1) using an amorphous silicon electronic portal imaging device (aSi-EPID) to capture the MLC segments during an IMRT delivery and (2) reconstituting a leaf sequence (LS) file based on the leaf end positions calculated from the MLC segments and their associated fractional monitor units. These EPID-measured LS files are then used to derive delivered fluence maps for dose reconstruction. The developed approach was tested on a pelvic phantom treated with a typical prostate IMRT plan. The delivered fluence maps, which were derived from the EPID-measured LS files, showed slight differences in the intensity levels compared with the corresponding planned ones. The dose distribution calculated with the delivered fluence maps showed a discernible difference in the high dose region when compared to that calculated with the planned fluence maps. The maximum dose in the former distribution was also 2.5% less than that in the latter one. The EPID-measured LS file can serve the same purpose as a MLC log files does for the derivation of the delivered fluence map and yet is independent of the leaf calibration. The approach also allows users who do not have access to MLC log files to probe the actual IMRT delivery and translate the information gained for dose reconstruction in adaptive radiation therapy.

IMAGE-GUIDED RADIOTHERAPY IN NEAR REAL TIME WITH INTENSITY- MODULATED RADIOTHERAPY MEGAVOLTAGE TREATMENT BEAM IMAGING

PURPOSE To utilize image-guided radiotherapy (IGRT) in near real time by obtaining and evaluating the online positions of implanted fiducials from continuous electronic portal imaging device (EPID) imaging of prostate intensity-modulated radiotherapy (IMRT) delivery. METHODS AND MATERIALS Upon initial setup using two orthogonal images, the three-dimensional (3D) positions of all implanted fiducial markers are obtained, and their expected two-dimensional (2D) locations in the beams-eye-view (BEV) projection are calculated for each treatment field. During IMRT beam delivery, EPID images of the megavoltage treatment beam are acquired in cine mode and subsequently analyzed to locate 2D locations of fiducials in the BEV. Simultaneously, 3D positions are estimated according to the current EPID image, information from the setup portal images, and images acquired at other gantry angles (the completed treatment fields). The measured 2D and 3D positions of each fiducial are compared with their expected 2D and 3D setup positions, respectively. Any displacements larger than a predefined tolerance may cause the treatment system to suspend the beam delivery and direct the therapists to reposition the patient. RESULTS Phantom studies indicate that the accuracy of 2D BEV and 3D tracking are better than 1 mm and 1.4 mm, respectively. A total of 7330 images from prostate treatments were acquired and analyzed, showing a maximum 2D displacement of 6.7 mm and a maximum 3D displacement of 6.9 mm over 34 fractions. CONCLUSIONS This EPID-based, real-time IGRT method can be implemented on any external beam machine with portal imaging capabilities without purchasing any additional equipment, and there is no extra dose delivered to the patient.