Ruijiang Li

Journal Articles

• First study of on-treatment volumetric imaging during respiratory gated VMAT. Medical physics Choi, K., Xing, L., Koong, A., Li, R. 2013; 40 (4): 040701-?

  Abstract

  To obtain on-treatment volumetric patient anatomy during respiratory gated volumetric modulated arc therapy (VMAT).On-board imaging device integrated with Linacs offers a viable tool for obtaining patient anatomy during radiation treatment delivery. In this study, the authors acquired beam-level kV images during gated VMAT treatments using a Varian TrueBeam™STx Linac. These kV projection images are triggered by a respiratory gating signal and can be acquired immediately before treatment MV beam on at every breathing cycle during delivery. Because the kV images are acquired with an on-board imaging device during a rotational arc therapy, they provide the patient anatomical information from many different angles or projection views (typically 20-40). To reconstruct the volumetric image representing patient anatomy during the VMAT treatment, the authors used a compressed sensing method with a fast first-order optimization algorithm. The conventional FDK reconstruction was also used for comparison purposes. The method was tested on a dynamic anthropomorphic physical phantom as well as a lung patient.The reconstructed volumetric images for a dynamic anthropomorphic physical phantom and a lung patient showed clearly visible soft-tissue target as well as other anatomical structures, with the proposed compressed sensing-based image reconstruction method. Compared with FDK, the compressed sensing method leads to a ? two and threefold increase in contrast-to-noise ratio around the target area in the phantom and patient case, respectively.The proposed technique provides on-treatment volumetric patient anatomy, with only a fraction (<10%) of the imaging dose used in conventional CBCT procedures. This anatomical information may be valuable for geometric verification and treatment guidance, and useful for verification of treatment dose delivery, accumulation, and adaptation in the future.

  View details for DOI 10.1118/1.4794925

  View details for PubMedID 23556870

• Real-time volumetric image reconstruction and 3D tumor localization based on a single x-ray projection image for lung cancer radiotherapy MEDICAL PHYSICS Li, R., Jia, X., Lewis, J. H., Gu, X., Folkerts, M., Men, C., Jiang, S. B. 2010; 37 (6): 2822-2826

  Abstract

  To develop an algorithm for real-time volumetric image reconstruction and 3D tumor localization based on a single x-ray projection image for lung cancer radiotherapy.Given a set of volumetric images of a patient at N breathing phases as the training data, deformable image registration was performed between a reference phase and the other N-1 phases, resulting in N-1 deformation vector fields (DVFs). These DVFs can be represented efficiently by a few eigenvectors and coefficients obtained from principal component analysis (PCA). By varying the PCA coefficients, new DVFs can be generated, which, when applied on the reference image, lead to new volumetric images. A volumetric image can then be reconstructed from a single projection image by optimizing the PCA coefficients such that its computed projection matches the measured one. The 3D location of the tumor can be derived by applying the inverted DVF on its position in the reference image. The algorithm was implemented on graphics processing units (GPUs) to achieve real-time efficiency. The training data were generated using a realistic and dynamic mathematical phantom with ten breathing phases. The testing data were 360 cone beam projections corresponding to one gantry rotation, simulated using the same phantom with a 50% increase in breathing amplitude.The average relative image intensity error of the reconstructed volumetric images is 6.9% +/- 2.4%. The average 3D tumor localization error is 0.8 +/- 0.5 mm. On an NVIDIA Tesla C1060 GPU card, the average computation time for reconstructing a volumetric image from each projection is 0.24 s (range: 0.17 and 0.35 s).The authors have shown the feasibility of reconstructing volumetric images and localizing tumor positions in 3D in near real-time from a single x-ray image.

  View details for DOI 10.1118/1.3426002

  View details for Web of Science ID 000278573100045

  View details for PubMedID 20632593

• Cone beam CT imaging with limited angle of projections and prior knowledge for volumetric verification of non-coplanar beam radiation therapy: a proof of concept study. Physics in medicine and biology Meng, B., Xing, L., Han, B., Koong, A., Chang, D., Cheng, J., Li, R. 2013; 58 (21): 7777-7789

  Abstract

  Non-coplanar beams are important for treatment of both cranial and noncranial tumors. Treatment verification of such beams with couch rotation/kicks, however, is challenging, particularly for the application of cone beam CT (CBCT). In this situation, only limited and unconventional imaging angles are feasible to avoid collision between the gantry, couch, patient, and on-board imaging system. The purpose of this work is to develop a CBCT verification strategy for patients undergoing non-coplanar radiation therapy. We propose an image reconstruction scheme that integrates a prior image constrained compressed sensing (PICCS) technique with image registration. Planning CT or CBCT acquired at the neutral position is rotated and translated according to the nominal couch rotation/translation to serve as the initial prior image. Here, the nominal couch movement is chosen to have a rotational error of 5° and translational error of 8 mm from the ground truth in one or more axes or directions. The proposed reconstruction scheme alternates between two major steps. First, an image is reconstructed using the PICCS technique implemented with total-variation minimization and simultaneous algebraic reconstruction. Second, the rotational/translational setup errors are corrected and the prior image is updated by applying rigid image registration between the reconstructed image and the previous prior image. The PICCS algorithm and rigid image registration are alternated iteratively until the registration results fall below a predetermined threshold. The proposed reconstruction algorithm is evaluated with an anthropomorphic digital phantom and physical head phantom. The proposed algorithm provides useful volumetric images for patient setup using projections with an angular range as small as 60°. It reduced the translational setup errors from 8 mm to generally <1 mm and the rotational setup errors from 5° to <1°. Compared with the PICCS algorithm alone, the integration of rigid registration significantly improved the reconstructed image quality, with a reduction of mostly 2-3 folds (up to 100) in root mean square image error. The proposed algorithm provides a remedy for solving the problem of non-coplanar CBCT reconstruction from limited angle of projections by combining the PICCS technique and rigid image registration in an iterative framework. In this proof of concept study, non-coplanar beams with couch rotations of 45° can be effectively verified with the CBCT technique.

  View details for DOI 10.1088/0031-9155/58/21/7777

  View details for PubMedID 24140954

• Automatic prostate tracking and motion assessment in volumetric modulated arc therapy with an electronic portal imaging device. International journal of radiation oncology, biology, physics Azcona, J. D., Li, R., Mok, E., Hancock, S., Xing, L. 2013; 86 (4): 762-768

  Abstract

  PURPOSE: To assess the prostate intrafraction motion in volumetric modulated arc therapy treatments using cine megavoltage (MV) images acquired with an electronic portal imaging device (EPID). METHODS AND MATERIALS: Ten prostate cancer patients were treated with volumetric modulated arc therapy using a Varian TrueBeam linear accelerator equipped with an EPID for acquiring cine MV images during treatment. Cine MV images acquisition was scheduled for single or multiple treatment fractions (between 1 and 8). A novel automatic fiducial detection algorithm that can handle irregular multileaf collimator apertures, field edges, fast leaf and gantry movement, and MV image noise and artifacts in patient anatomy was used. All sets of images (approximately 25,000 images in total) were analyzed to measure the positioning accuracy of implanted fiducial markers and assess the prostate movement. RESULTS: Prostate motion can vary greatly in magnitude among different patients. Different motion patterns were identified, showing its unpredictability. The mean displacement and standard deviation of the intrafraction motion was generally less than 2.0 ± 2.0 mm in each of the spatial directions. In certain patients, however, the percentage of the treatment time in which the prostate is displaced more than 5 mm from its planned position in at least 1 spatial direction was 10% or more. The maximum prostate displacement observed was 13.3 mm. CONCLUSION: Prostate tracking and motion assessment was performed with MV imaging and an EPID. The amount of prostate motion observed suggests that patients will benefit from its real-time monitoring. Megavoltage imaging can provide the basis for real-time prostate tracking using conventional linear accelerators.

  View details for DOI 10.1016/j.ijrobp.2013.03.007

  View details for PubMedID 23608236

• Improving IMRT delivery efficiency with reweighted L1-minimization for inverse planning MEDICAL PHYSICS Kim, H., Becker, S., Lee, R., Lee, S., Shin, S., Candes, E., Xing, L., Li, R. 2013; 40 (7)

  Abstract

  Purpose: This study presents an improved technique to further simplify the fluence-map in intensity modulated radiation therapy (IMRT) inverse planning, thereby reducing plan complexity and improving delivery efficiency, while maintaining the plan quality.Methods: First-order total-variation (TV) minimization (min.) based on L1-norm has been proposed to reduce the complexity of fluence-map in IMRT by generating sparse fluence-map variations. However, with stronger dose sparing to the critical structures, the inevitable increase in the fluence-map complexity can lead to inefficient dose delivery. Theoretically, L0-min. is the ideal solution for the sparse signal recovery problem, yet practically intractable due to its nonconvexity of the objective function. As an alternative, the authors use the iteratively reweighted L1-min. technique to incorporate the benefits of the L0-norm into the tractability of L1-min. The weight multiplied to each element is inversely related to the magnitude of the corresponding element, which is iteratively updated by the reweighting process. The proposed penalizing process combined with TV min. further improves sparsity in the fluence-map variations, hence ultimately enhancing the delivery efficiency. To validate the proposed method, this work compares three treatment plans obtained from quadratic min. (generally used in clinic IMRT), conventional TV min., and our proposed reweighted TV min. techniques, implemented by a large-scale L1-solver (template for first-order conic solver), for five patient clinical data. Criteria such as conformation number (CN), modulation index (MI), and estimated treatment time are employed to assess the relationship between the plan quality and delivery efficiency.Results: The proposed method yields simpler fluence-maps than the quadratic and conventional TV based techniques. To attain a given CN and dose sparing to the critical organs for 5 clinical cases, the proposed method reduces the number of segments by 10-15 and 30-35, relative to TV min. and quadratic min. based plans, while MIs decreases by about 20%-30% and 40%-60% over the plans by two existing techniques, respectively. With such conditions, the total treatment time of the plans obtained from our proposed method can be reduced by 12-30 s and 30-80 s mainly due to greatly shorter multileaf collimator (MLC) traveling time in IMRT step-and-shoot delivery.Conclusions: The reweighted L1-minimization technique provides a promising solution to simplify the fluence-map variations in IMRT inverse planning. It improves the delivery efficiency by reducing the entire segments and treatment time, while maintaining the plan quality in terms of target conformity and critical structure sparing.

  View details for DOI 10.1118/1.4811100

  View details for Web of Science ID 000321272200023

  View details for PubMedID 23822423

• An adaptive planning strategy for station parameter optimized radiation therapy (SPORT): Segmentally boosted VMAT. Medical physics Li, R., Xing, L. 2013; 40 (5): 050701-?

  Abstract

  Conventional volumetric modulated arc therapy (VMAT) discretizes the angular space into equally spaced control points during planning and then optimizes the apertures and weights of the control points. The aperture at an angle in between two control points is obtained through interpolation. This approach tacitly ignores the differential need for intensity modulation of different angles. As such, multiple arcs are often required, which may oversample some angle(s) and undersample others. The purpose of this work is to develop a segmentally boosted VMAT scheme to eliminate the need for multiple arcs in VMAT treatment with improved dose distribution and?or delivery efficiency.The essence of the new treatment scheme is how to identify the need of individual angles for intensity modulation and to provide the necessary beam intensity modulation for those beam angles that need it. We introduce a "demand metric" at each control point to decide which station or control points need intensity modulation. To boost the modulation at selected stations, additional segments are added in the vicinity of the selected stations. The added segments are then optimized together with the original set of station or control points as a whole. The authors apply the segmentally boosted planning technique to four previously treated clinical cases: two head and neck (HN) cases, one prostate case, and one liver case. The proposed planning technique is compared with conventional one-arc and two-arc VMAT.The proposed segmentally boosted VMAT technique achieves better critical structure sparing than one-arc VMAT with similar or better target coverage in all four clinical cases. The segmentally boosted VMAT also outperforms two-arc VMAT for the two complicated HN cases, yet with ?30% reduction in the machine monitor units (MUs) relative to two-arc VMAT, which leads to less leakage?scatter dose to the patient and can potentially translate into faster dose delivery. For the less challenging prostate and liver cases, similar critical structure sparing as the two-arc VMAT plans was obtained using the segmentally boosted VMAT. The benefit for the two simpler cases is the reduction of MUs and improvement of treatment delivery efficiency.Segmentally boosted VMAT achieves better dose conformality and?or reduced MUs through effective consideration of the need of individual beam angles for intensity modulation. Elimination of the need for multiple arcs in rotational arc therapy while improving the dose distribution should lead to improved workflow and treatment efficacy, thus may have significant implication to radiation oncology practice.

  View details for DOI 10.1118/1.4802748

  View details for PubMedID 23635247

• Development and clinical evaluation of automatic fiducial detection for tumor tracking in cine megavoltage images during volumetric modulated arc therapy MEDICAL PHYSICS Azcona, J. D., Li, R., Mok, E., Hancock, S., Xing, L. 2013; 40 (3)

  Abstract

  Real-time tracking of implanted fiducials in cine megavoltage (MV) imaging during volumetric modulated arc therapy (VMAT) delivery is complicated due to the inherent low contrast of MV images and potential blockage of dynamic leaves configurations. The purpose of this work is to develop a clinically practical autodetection algorithm for motion management during VMAT.The expected field-specific segments and the planned fiducial position from the Eclipse (Varian Medical Systems, Palo Alto, CA) treatment planning system were projected onto the MV images. The fiducials were enhanced by applying a Laplacian of Gaussian filter in the spatial domain for each image, with a blob-shaped object as the impulse response. The search of implanted fiducials was then performed on a region of interest centered on the projection of the fiducial when it was within an open field including the case when it was close to the field edge or partially occluded by the leaves. A universal template formula was proposed for template matching and normalized cross correlation was employed for its simplicity and computational efficiency. The search region for every image was adaptively updated through a prediction model that employed the 3D position of the fiducial estimated from the localized positions in previous images. This prediction model allowed the actual fiducial position to be tracked dynamically and was used to initialize the search region. The artifacts caused by electronic interference during the acquisition were effectively removed. A score map was computed by combining both morphological information and image intensity. The pixel location with the highest score was selected as the detected fiducial position. The sets of cine MV images taken during treatment were analyzed with in-house developed software written in MATLAB (The Mathworks, Inc., Natick, MA). Five prostate patients were analyzed to assess the algorithm performance by measuring their positioning accuracy during treatment.The algorithm was able to accurately localize the fiducial position on MV images with success rates of more than 90% per case. The percentage of images in which each fiducial was localized in the studied cases varied between 23% and 65%, with at least one fiducial having been localized between 40% and 95% of the images. This depended mainly on the modulation of the plan and fiducial blockage. The prostate movement in the presented cases varied between 0.8 and 3.5 mm (mean values). The maximum displacement detected among all patients was of 5.7 mm.An algorithm for automatic detection of fiducial markers in cine MV images has been developed and tested with five clinical cases. Despite the challenges posed by complex beam aperture shapes, fiducial localization close to the field edge, partial occlusion of fiducials, fast leaf and gantry movement, and inherently low MV image quality, good localization results were achieved in patient images. This work provides a technique for enabling real-time accurate fiducial detection and tumor tracking during VMAT treatments without the use of extra imaging dose.

  View details for DOI 10.1118/1.4791646

  View details for Web of Science ID 000316369400011

  View details for PubMedID 23464303

• Evaluation of 3D fluoroscopic image generation from a single planar treatment image on patient data with a modified XCAT phantom PHYSICS IN MEDICINE AND BIOLOGY Mishra, P., Li, R., St James, S., Mak, R. H., Williams, C. L., Yue, Y., Berbeco, R. I., Lewis, J. H. 2013; 58 (4): 841-858

  Abstract

  Accurate understanding and modeling of respiration-induced uncertainties is essential in image-guided radiotherapy. Explicit modeling of the overall lung motion and interaction among different organs promises to be a useful approach. Recently, preliminary studies on 3D fluoroscopic treatment imaging and tumor localization based on principal component analysis motion models and cost function optimization have shown encouraging results. However, the performance of this technique for varying breathing parameters and under realistic conditions remains unclear and thus warrants further investigation. In this work, we present a systematic evaluation of a 3D fluoroscopic image generation algorithm via two different approaches. In the first approach, the model's accuracy is tested for changing parameters for sinusoidal breathing. These parameters include changing respiratory motion amplitude, period and baseline shift. The effects of setup error, imaging noise and different tumor sizes are also examined. In the second approach, we test the model for anthropomorphic images obtained from a modified XCAT phantom. This set of experiments is important as all the underlying breathing parameters are simultaneously tested, as in realistic clinical conditions. Based on our simulation results for more than 250 s of breathing data for eight different lung patients, the overall tumor localization accuracies of the model in left-right, anterior-posterior and superior-inferior directions are 0.1 ± 0.1, 0.5 ± 0.5 and 0.8 ± 0.8 mm, respectively. 3D tumor centroid localization accuracy is 1.0 ± 0.9 mm.

  View details for DOI 10.1088/0031-9155/58/4/841

  View details for Web of Science ID 000314396800008

• Accurate Respiration Measurement Using DC-Coupled Continuous-Wave Radar Sensor for Motion-Adaptive Cancer Radiotherapy IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING Gu, C., Li, R., Zhang, H., Fung, A. Y., Torres, C., Jiang, S. B., Li, C. 2012; 59 (11): 3117-3123

  Abstract

  Accurate respiration measurement is crucial in motion-adaptive cancer radiotherapy. Conventional methods for respiration measurement are undesirable because they are either invasive to the patient or do not have sufficient accuracy. In addition, measurement of external respiration signal based on conventional approaches requires close patient contact to the physical device which often causes patient discomfort and undesirable motion during radiation dose delivery. In this paper, a dc-coupled continuous-wave radar sensor was presented to provide a noncontact and noninvasive approach for respiration measurement. The radar sensor was designed with dc-coupled adaptive tuning architectures that include RF coarse-tuning and baseband fine-tuning, which allows the radar sensor to precisely measure movement with stationary moment and always work with the maximum dynamic range. The accuracy of respiration measurement with the proposed radar sensor was experimentally evaluated using a physical phantom, human subject, and moving plate in a radiotherapy environment. It was shown that respiration measurement with radar sensor while the radiation beam is on is feasible and the measurement has a submillimeter accuracy when compared with a commercial respiration monitoring system which requires patient contact. The proposed radar sensor provides accurate, noninvasive, and noncontact respiration measurement and therefore has a great potential in motion-adaptive radiotherapy.

  View details for DOI 10.1109/TBME.2012.2206591

  View details for Web of Science ID 000310154700016

  View details for PubMedID 22759434

• 4D cone beam CT via spatiotemporal tensor framelet MEDICAL PHYSICS Gao, H., Li, R., Lin, Y., Xing, L. 2012; 39 (11): 6943-6946

  Abstract

  On-board 4D cone beam CT (4DCBCT) offers respiratory phase-resolved volumetric imaging, and improves the accuracy of target localization in image guided radiation therapy. However, the clinical utility of this technique has been greatly impeded by its degraded image quality, prolonged imaging time, and increased imaging dose. The purpose of this letter is to develop a novel iterative 4DCBCT reconstruction method for improved image quality, increased imaging speed, and reduced imaging dose.The essence of this work is to introduce the spatiotemporal tensor framelet (STF), a high-dimensional tensor generalization of the 1D framelet for 4DCBCT, to effectively take into account of highly correlated and redundant features of the patient anatomy during respiration, in a multilevel fashion with multibasis sparsifying transform. The STF-based algorithm is implemented on a GPU platform for improved computational efficiency. To evaluate the method, 4DCBCT full-fan scans were acquired within 30 s, with a gantry rotation of 200°; STF is also compared with a state-of-art reconstruction method via spatiotemporal total variation regularization.Both the simulation and experimental results demonstrate that STF-based reconstruction achieved superior image quality. The reconstruction of 20 respiratory phases took less than 10 min on an NVIDIA Tesla C2070 GPU card. The STF codes are available at https://sites.google.com/site/spatiotemporaltensorframelet.By effectively utilizing the spatiotemporal coherence of the patient anatomy among different respiratory phases in a multilevel fashion with multibasis sparsifying transform, the proposed STF method potentially enables fast and low-dose 4DCBCT with improved image quality.

  View details for DOI 10.1118/1.4762288

  View details for Web of Science ID 000310726300042

  View details for PubMedID 23127087

• Real-time tumor motion estimation using respiratory surrogate via memory-based learning PHYSICS IN MEDICINE AND BIOLOGY Li, R., Lewis, J. H., Berbeco, R. I., Xing, L. 2012; 57 (15): 4771-4786

  Abstract

  Respiratory tumor motion is a major challenge in radiation therapy for thoracic and abdominal cancers. Effective motion management requires an accurate knowledge of the real-time tumor motion. External respiration monitoring devices (optical, etc) provide a noninvasive, non-ionizing, low-cost and practical approach to obtain the respiratory signal. Due to the highly complex and nonlinear relations between tumor and surrogate motion, its ultimate success hinges on the ability to accurately infer the tumor motion from respiratory surrogates. Given their widespread use in the clinic, such a method is critically needed. We propose to use a powerful memory-based learning method to find the complex relations between tumor motion and respiratory surrogates. The method first stores the training data in memory and then finds relevant data to answer a particular query. Nearby data points are assigned high relevance (or weights) and conversely distant data are assigned low relevance. By fitting relatively simple models to local patches instead of fitting one single global model, it is able to capture highly nonlinear and complex relations between the internal tumor motion and external surrogates accurately. Due to the local nature of weighting functions, the method is inherently robust to outliers in the training data. Moreover, both training and adapting to new data are performed almost instantaneously with memory-based learning, making it suitable for dynamically following variable internal/external relations. We evaluated the method using respiratory motion data from 11 patients. The data set consists of simultaneous measurement of 3D tumor motion and 1D abdominal surface (used as the surrogate signal in this study). There are a total of 171 respiratory traces, with an average peak-to-peak amplitude of ?15 mm and average duration of ?115 s per trace. Given only 5 s (roughly one breath) pretreatment training data, the method achieved an average 3D error of 1.5 mm and 95th percentile error of 3.4 mm on unseen test data. The average 3D error was further reduced to 1.4 mm when the model was tuned to its optimal setting for each respiratory trace. In one trace where a few outliers are present in the training data, the proposed method achieved an error reduction of as much as ?50% compared with the best linear model (1.0 mm versus 2.1 mm). The memory-based learning technique is able to accurately capture the highly complex and nonlinear relations between tumor and surrogate motion in an efficient manner (a few milliseconds per estimate). Furthermore, the algorithm is particularly suitable to handle situations where the training data are contaminated by large errors or outliers. These desirable properties make it an ideal candidate for accurate and robust tumor gating/tracking using respiratory surrogates.

  View details for DOI 10.1088/0031-9155/57/15/4771

  View details for Web of Science ID 000306521900007

  View details for PubMedID 22772042

• Intrafraction Verification of Gated RapidArc by Using Beam-Level Kilovoltage X-Ray Images INTERNATIONAL JOURNAL OF RADIATION ONCOLOGY BIOLOGY PHYSICS Li, R., Mok, E., Chang, D. T., Daly, M., Loo, B. W., Diehn, M., Quynh-Thu Le, Q. T., Koong, A., Xing, L. 2012; 83 (5): E709-E715

  Abstract

  To verify the geometric accuracy of gated RapidArc treatment using kV images acquired during dose delivery.Twenty patients were treated using the gated RapidArc technique with a Varian TrueBeam STx linear accelerator. One to 7 metallic fiducial markers were implanted inside or near the tumor target before treatment simulation. For patient setup and treatment verification purposes, the internal target volume (ITV) was created, corresponding to each implanted marker. The gating signal was generated from the Real-time Position Management (RPM) system. At the beginning of each fraction, individualized respiratory gating amplitude thresholds were set based on fluoroscopic image guidance. During the treatment, we acquired kV images immediately before MV beam-on at every breathing cycle, using the on-board imaging system. After the treatment, all implanted markers were detected, and their 3-dimensional (3D) positions in the patient were estimated using software developed in-house. The distance from the marker to the corresponding ITV was calculated for each patient by averaging over all markers and all fractions.The average 3D distance between the markers and their ITVs was 0.8 ± 0.5 mm (range, 0-1.7 mm) and was 2.1 ± 1.2 mm at the 95th percentile (range, 0-3.8 mm). On average, a left-right margin of 0.6 mm, an anterior-posterior margin of 0.8 mm, and a superior-inferior margin of 1.5 mm is required to account for 95% of the intrafraction uncertainty in RPM-based RapidArc gating.To our knowledge, this is the first clinical report of intrafraction verification of respiration-gated RapidArc treatment in stereotactic ablative radiation therapy. For some patients, the markers deviated significantly from the ITV by more than 2 mm at the beginning of the MV beam-on. This emphasizes the need for gating techniques with beam-on/-off controlled directly by the actual position of the tumor target instead of external surrogates such as RPM.

  View details for DOI 10.1016/j.ijrobp.2012.03.006

  View details for Web of Science ID 000306128100022

  View details for PubMedID 22554582

• Efficient IMRT inverse planning with a new L1-solver: template for first-order conic solver PHYSICS IN MEDICINE AND BIOLOGY Kim, H., Suh, T., Lee, R., Xing, L., Li, R. 2012; 57 (13): 4139-4153

  Abstract

  Intensity modulated radiation therapy (IMRT) inverse planning using total-variation (TV) regularization has been proposed to reduce the complexity of fluence maps and facilitate dose delivery. Conventionally, the optimization problem with L-1 norm is solved with quadratic programming (QP), which is time consuming and memory expensive due to the second-order Newton update. This study proposes to use a new algorithm, template for first-order conic solver (TFOCS), for fast and memory-efficient optimization in IMRT inverse planning. The TFOCS utilizes dual-variable updates and first-order approaches for TV minimization without the need to compute and store the enlarged Hessian matrix required for Newton update in the QP technique. To evaluate the effectiveness and efficiency of the proposed method, two clinical cases were used for IMRT inverse planning: a head and neck case and a prostate case. For comparison, the conventional QP-based method for the TV form was adopted to solve the fluence map optimization problem in the above two cases. The convergence criteria and algorithm parameters were selected to achieve similar dose conformity for a fair comparison between the two methods. Compared with conventional QP-based approach, the proposed TFOCS-based method shows a remarkable improvement in computational efficiency for fluence map optimization, while maintaining the conformal dose distribution. Compared with QP-based algorithms, the computational speed using TFOCS for fluence optimization is increased by a factor of 4 to 6, and at the same time the memory requirement is reduced by a factor of 3 to 4. Therefore, TFOCS provides an effective, fast and memory-efficient method for IMRT inverse planning. The unique features of the approach should be particularly important in inverse planning involving a large number of beams, such as in VMAT and dense angularly sampled and sparse intensity modulated radiation therapy (DASSIM-RT).

  View details for DOI 10.1088/0031-9155/57/13/4139

  View details for Web of Science ID 000305803600006

  View details for PubMedID 22683930

• Dose optimization with first-order total-variation minimization for dense angularly sampled and sparse intensity modulated radiation therapy (DASSIM-RT) MEDICAL PHYSICS Kim, H., Li, R., Lee, R., Goldstein, T., Boyd, S., Candes, E., Xing, L. 2012; 39 (7): 4316-4327

  Abstract

  A new treatment scheme coined as dense angularly sampled and sparse intensity modulated radiation therapy (DASSIM-RT) has recently been proposed to bridge the gap between IMRT and VMAT. By increasing the angular sampling of radiation beams while eliminating dispensable segments of the incident fields, DASSIM-RT is capable of providing improved conformity in dose distributions while maintaining high delivery efficiency. The fact that DASSIM-RT utilizes a large number of incident beams represents a major computational challenge for the clinical applications of this powerful treatment scheme. The purpose of this work is to provide a practical solution to the DASSIM-RT inverse planning problem.The inverse planning problem is formulated as a fluence-map optimization problem with total-variation (TV) minimization. A newly released L1-solver, template for first-order conic solver (TFOCS), was adopted in this work. TFOCS achieves faster convergence with less memory usage as compared with conventional quadratic programming (QP) for the TV form through the effective use of conic forms, dual-variable updates, and optimal first-order approaches. As such, it is tailored to specifically address the computational challenges of large-scale optimization in DASSIM-RT inverse planning. Two clinical cases (a prostate and a head and neck case) are used to evaluate the effectiveness and efficiency of the proposed planning technique. DASSIM-RT plans with 15 and 30 beams are compared with conventional IMRT plans with 7 beams in terms of plan quality and delivery efficiency, which are quantified by conformation number (CN), the total number of segments and modulation index, respectively. For optimization efficiency, the QP-based approach was compared with the proposed algorithm for the DASSIM-RT plans with 15 beams for both cases.Plan quality improves with an increasing number of incident beams, while the total number of segments is maintained to be about the same in both cases. For the prostate patient, the conformation number to the target was 0.7509, 0.7565, and 0.7611 with 80 segments for IMRT with 7 beams, and DASSIM-RT with 15 and 30 beams, respectively. For the head and neck (HN) patient with a complicated target shape, conformation numbers of the three treatment plans were 0.7554, 0.7758, and 0.7819 with 75 segments for all beam configurations. With respect to the dose sparing to the critical structures, the organs such as the femoral heads in the prostate case and the brainstem and spinal cord in the HN case were better protected with DASSIM-RT. For both cases, the delivery efficiency has been greatly improved as the beam angular sampling increases with the similar or better conformal dose distribution. Compared with conventional quadratic programming approaches, first-order TFOCS-based optimization achieves far faster convergence and smaller memory requirements in DASSIM-RT.The new optimization algorithm TFOCS provides a practical and timely solution to the DASSIM-RT or other inverse planning problem requiring large memory space. The new treatment scheme is shown to outperform conventional IMRT in terms of dose conformity to both the targetand the critical structures, while maintaining high delivery efficiency.

  View details for DOI 10.1118/1.4729717

  View details for Web of Science ID 000306893000029

  View details for PubMedID 22830765

• Evaluation of the geometric accuracy of surrogate-based gated VMAT using intrafraction kilovoltage x-ray images MEDICAL PHYSICS Li, R., Mok, E., Han, B., Koong, A., Xing, L. 2012; 39 (5): 2686-2693

  Abstract

  To evaluate the geometric accuracy of beam targeting in external surrogate-based gated volumetric modulated arc therapy (VMAT) using kilovoltage (kV) x-ray images acquired during dose delivery.Gated VMAT treatments were delivered using a Varian TrueBeam STx Linac for both physical phantoms and patients. Multiple gold fiducial markers were implanted near the target. The reference position was created for each implanted marker, representing its correct position at the gating threshold. The gating signal was generated from the RPM system. During the treatment, kV images were acquired immediately before MV beam-on at every breathing cycle, using the on-board imaging system. All implanted markers were detected and their 3D positions were estimated using in-house developed software. The positioning error of a marker is defined as the distance of the marker from its reference position for each frame of the images. The overall error of the system is defined as the average over all markers. For the phantom study, both sinusoidal motion (1D and 3D) and real human respiratory motion was simulated for the target and surrogate. In the baseline case, the two motions were synchronized for the first treatment fraction. To assess the effects of surrogate-target correlation on the geometric accuracy, a phase shift of 5% and 10% between the two motions was introduced. For the patient study, intrafraction kV images of five stereotactic body radiotherapy (SBRT) patients were acquired for one or two fractions.For the phantom study, a high geometric accuracy was achieved in the baseline case (average error: 0.8 mm in the superior-inferior or SI direction). However, the treatment delivery is prone to geometric errors if changes in the target-surrogate relation occur during the treatment: the average error was increased to 2.3 and 4.7 mm for the phase shift of 5% and 10%, respectively. Results obtained with real human respiratory curves show a similar trend. For a target with 3D motion, the technique is able to detect geometric errors in the left-right (LR) and anterior-posterior (AP) directions. For the patient study, the average intrafraction positioning errors are 0.8, 0.9, and 1.4 mm and 95th percentile errors are 1.7, 2.1, and 2.7 mm in the LR, AP, and SI directions, respectively.The correlation between external surrogate and internal target motion is crucial to ensure the geometric accuracy of surrogate-based gating. Real-time guidance based on kV x-ray images overcomes the potential issues in surrogate-based gating and can achieve accurate beam targeting in gated VMAT.

  View details for DOI 10.1118/1.4704729

  View details for Web of Science ID 000303604300039

  View details for PubMedID 22559639

• Improved compressed sensing-based cone-beam CT reconstruction using adaptive prior image constraints PHYSICS IN MEDICINE AND BIOLOGY Lee, H., Xing, L., Davidi, R., Li, R., Qian, J., Lee, R. 2012; 57 (8): 2287-2307

  Abstract

  Volumetric cone-beam CT (CBCT) images are acquired repeatedly during a course of radiation therapy and a natural question to ask is whether CBCT images obtained earlier in the process can be utilized as prior knowledge to reduce patient imaging dose in subsequent scans. The purpose of this work is to develop an adaptive prior image constrained compressed sensing (APICCS) method to solve this problem. Reconstructed images using full projections are taken on the first day of radiation therapy treatment and are used as prior images. The subsequent scans are acquired using a protocol of sparse projections. In the proposed APICCS algorithm, the prior images are utilized as an initial guess and are incorporated into the objective function in the compressed sensing (CS)-based iterative reconstruction process. Furthermore, the prior information is employed to detect any possible mismatched regions between the prior and current images for improved reconstruction. For this purpose, the prior images and the reconstructed images are classified into three anatomical regions: air, soft tissue and bone. Mismatched regions are identified by local differences of the corresponding groups in the two classified sets of images. A distance transformation is then introduced to convert the information into an adaptive voxel-dependent relaxation map. In constructing the relaxation map, the matched regions (unchanged anatomy) between the prior and current images are assigned with smaller weight values, which are translated into less influence on the CS iterative reconstruction process. On the other hand, the mismatched regions (changed anatomy) are associated with larger values and the regions are updated more by the new projection data, thus avoiding any possible adverse effects of prior images. The APICCS approach was systematically assessed by using patient data acquired under standard and low-dose protocols for qualitative and quantitative comparisons. The APICCS method provides an effective way for us to enhance the image quality at the matched regions between the prior and current images compared to the existing PICCS algorithm. Compared to the current CBCT imaging protocols, the APICCS algorithm allows an imaging dose reduction of 10-40 times due to the greatly reduced number of projections and lower x-ray tube current level coming from the low-dose protocol.

  View details for DOI 10.1088/0031-9155/57/8/2287

  View details for Web of Science ID 000302567100013

  View details for PubMedID 22460008

• Response to "Comment on 'Bridging the gap between IMRT and VMAT: Dense angularly sampled and sparse intensity modulated radiation therapy'" [Med. Phys. 38, 4912-4919 (2011)] MEDICAL PHYSICS Li, R., Xing, L. 2012; 39 (3): 1676-1676

  View details for DOI 10.1118/1.3687906

  View details for Web of Science ID 000301503400053

  View details for PubMedID 22380399

• On a PCA-based lung motion model PHYSICS IN MEDICINE AND BIOLOGY Li, R., Lewis, J. H., Jia, X., Zhao, T., Liu, W., Wuenschel, S., Lamb, J., Yang, D., Low, D. A., Jiang, S. B. 2011; 56 (18): 6009-6030

  Abstract

  Respiration-induced organ motion is one of the major uncertainties in lung cancer radiotherapy and is crucial to be able to accurately model the lung motion. Most work so far has focused on the study of the motion of a single point (usually the tumor center of mass), and much less work has been done to model the motion of the entire lung. Inspired by the work of Zhang et al (2007 Med. Phys. 34 4772-81), we believe that the spatiotemporal relationship of the entire lung motion can be accurately modeled based on principle component analysis (PCA) and then a sparse subset of the entire lung, such as an implanted marker, can be used to drive the motion of the entire lung (including the tumor). The goal of this work is twofold. First, we aim to understand the underlying reason why PCA is effective for modeling lung motion and find the optimal number of PCA coefficients for accurate lung motion modeling. We attempt to address the above important problems both in a theoretical framework and in the context of real clinical data. Second, we propose a new method to derive the entire lung motion using a single internal marker based on the PCA model. The main results of this work are as follows. We derived an important property which reveals the implicit regularization imposed by the PCA model. We then studied the model using two mathematical respiratory phantoms and 11 clinical 4DCT scans for eight lung cancer patients. For the mathematical phantoms with cosine and an even power (2n) of cosine motion, we proved that 2 and 2n PCA coefficients and eigenvectors will completely represent the lung motion, respectively. Moreover, for the cosine phantom, we derived the equivalence conditions for the PCA motion model and the physiological 5D lung motion model (Low et al 2005 Int. J. Radiat. Oncol. Biol. Phys. 63 921-9). For the clinical 4DCT data, we demonstrated the modeling power and generalization performance of the PCA model. The average 3D modeling error using PCA was within 1 mm (0.7 ± 0.1 mm). When a single artificial internal marker was used to derive the lung motion, the average 3D error was found to be within 2 mm (1.8 ± 0.3 mm) through comprehensive statistical analysis. The optimal number of PCA coefficients needs to be determined on a patient-by-patient basis and two PCA coefficients seem to be sufficient for accurate modeling of the lung motion for most patients. In conclusion, we have presented thorough theoretical analysis and clinical validation of the PCA lung motion model. The feasibility of deriving the entire lung motion using a single marker has also been demonstrated on clinical data using a simulation approach.

  View details for DOI 10.1088/0031-9155/56/18/015

  View details for Web of Science ID 000294787300017

  View details for PubMedID 21865624

• Mitigation of motion artifacts in CBCT of lung tumors based on tracked tumor motion during CBCT acquisition PHYSICS IN MEDICINE AND BIOLOGY Lewis, J. H., Li, R., Jia, X., Watkins, W. T., Lou, Y., Song, W. Y., Jiang, S. B. 2011; 56 (17): 5485-5502

  Abstract

  An algorithm capable of mitigating respiratory motion blurring artifacts in cone-beam computed tomography (CBCT) lung tumor images based on the motion of the tumor during the CBCT scan is developed. The tumor motion trajectory and probability density function (PDF) are reconstructed from the acquired CBCT projection images using a recently developed algorithm Lewis et al (2010 Phys. Med. Biol. 55 2505-22). Assuming that the effects of motion blurring can be represented by convolution of the static lung (or tumor) anatomy with the motion PDF, a cost function is defined, consisting of a data fidelity term and a total variation regularization term. Deconvolution is performed through iterative minimization of this cost function. The algorithm was tested on digital respiratory phantom, physical respiratory phantom and patient data. A clear qualitative improvement is evident in the deblurred images as compared to the motion-blurred images for all cases. Line profiles show that the tumor boundaries are more accurately and clearly represented in the deblurred images. The normalized root-mean-squared error between the images used as ground truth and the motion-blurred images are 0.29, 0.12 and 0.30 in the digital phantom, physical phantom and patient data, respectively. Deblurring reduces the corresponding values to 0.13, 0.07 and 0.19. Application of a -700 HU threshold to the digital phantom results in tumor dimension measurements along the superior-inferior axis of 2.8, 1.8 and 1.9 cm in the motion-blurred, ground truth and deblurred images, respectively. Corresponding values for the physical phantom are 3.4, 2.7 and 2.7 cm. A threshold of -500 HU applied to the patient case gives measurements of 3.1, 1.6 and 1.7 cm along the SI axis in the CBCT, 4DCT and deblurred images, respectively. This technique could provide more accurate information about a lung tumor's size and shape on the day of treatment.

  View details for DOI 10.1088/0031-9155/56/17/003

  View details for Web of Science ID 000294786400006

  View details for PubMedID 21813959

• Bridging the gap between IMRT and VMAT: Dense angularly sampled and sparse intensity modulated radiation therapy MEDICAL PHYSICS Li, R., Xing, L. 2011; 38 (9): 4912-4919

  Abstract

  To propose an alternative radiation therapy (RT) planning and delivery scheme with optimal angular beam sampling and intrabeam modulation for improved dose distribution while maintaining high delivery efficiency.In the proposed approach, coined as dense angularly sampled and sparse intensity modulated RT (DASSIM-RT), a large number of beam angles are used to increase the angular sampling, leading to potentially more conformal dose distributions as compared to conventional IMRT. At the same time, intensity modulation of the incident beams is simplified to eliminate the dispensable segments, compensating the increase in delivery time caused by the increased number of beams and facilitating the plan delivery. In a sense, the proposed approach shifts and transforms, in an optimal fashion, some of the beam segments in conventional IMRT to the added beams. For newly available digital accelerators, the DASSIM-RT delivery can be made very efficient by concatenating the beams so that they can be delivered sequentially without operator's intervention. Different from VMAT, the level of intensity modulation in DASSIS-RT is field specific and optimized to meet the need of each beam direction. Three clinical cases (a head and neck (HN) case, a pancreas case, and a lung case) are used to evaluate the proposed RT scheme. DASSIM-RT, VMAT, and conventional IMRT plans are compared quantitatively in terms of the conformality index (CI) and delivery efficiency.Plan quality improves generally with the number and intensity modulation of the incident beams. For a fixed number of beams or fixed level of intensity modulation, the improvement saturates after the intensity modulation or number of beams reaches to a certain level. An interplay between the two variables is observed and the saturation point depends on the values of both variables. For all the cases studied here, the CI of DASSIM-RT with 15 beams and 5 intensity levels (0.90, 0.79, and 0.84 for the HN, pancreas, and lung cases, respectively) is similar with that of conventional IMRT with seven beams and ten intensity levels (0.88, 0.79, and 0.83) and is higher than that of single-arc VMAT (0.75, 0.75, and 0.82). It is also found that the DASSIM-RT plans generally have better sparing of organs-at-risk than IMRT plans. It is estimated that the dose delivery time of DASSIM-RT with 15 beams and 5 intensity levels is about 4.5, 4.4, and 4.2 min for the HN, pancreas, and lung case, respectively, similar to that of IMRT plans with 7 beams and 10 intensity levels.DASSIS-RT bridges the gap between IMRT and VMAT and allows optimal sampling of angular space and intrabeam modulation, thus it provides improved conformity in dose distributions while maintaining high delivery efficiency.

  View details for DOI 10.1118/1.3618736

  View details for Web of Science ID 000294482900002

  View details for PubMedID 21978036

• A Bayesian approach to real-time 3D tumor localization via monoscopic x-ray imaging during treatment delivery MEDICAL PHYSICS Li, R., Fahimian, B. P., Xing, L. 2011; 38 (7): 4205-4214

  Abstract

  Monoscopic x-ray imaging with on-board kV devices is an attractive approach for real-time image guidance in modern radiation therapy such as VMAT or IMRT, but it falls short in providing reliable information along the direction of imaging x-ray. By effectively taking consideration of projection data at prior times and/or angles through a Bayesian formalism, the authors develop an algorithm for real-time and full 3D tumor localization with a single x-ray imager during treatment delivery.First, a prior probability density function is constructed using the 2D tumor locations on the projection images acquired during patient setup. Whenever an x-ray image is acquired during the treatment delivery, the corresponding 2D tumor location on the imager is used to update the likelihood function. The unresolved third dimension is obtained by maximizing the posterior probability distribution. The algorithm can also be used in a retrospective fashion when all the projection images during the treatment delivery are used for 3D localization purposes. The algorithm does not involve complex optimization of any model parameter and therefore can be used in a "plug-and-play" fashion. The authors validated the algorithm using (1) simulated 3D linear and elliptic motion and (2) 3D tumor motion trajectories of a lung and a pancreas patient reproduced by a physical phantom. Continuous kV images were acquired over a full gantry rotation with the Varian TrueBeam on-board imaging system. Three scenarios were considered: fluoroscopic setup, cone beam CT setup, and retrospective analysis.For the simulation study, the RMS 3D localization error is 1.2 and 2.4 mm for the linear and elliptic motions, respectively. For the phantom experiments, the 3D localization error is < 1 mm on average and < 1.5 mm at 95th percentile in the lung and pancreas cases for all three scenarios. The difference in 3D localization error for different scenarios is small and is not statistically significant.The proposed algorithm eliminates the need for any population based model parameters in monoscopic image guided radiotherapy and allows accurate and real-time 3D tumor localization on current standard LINACs with a single x-ray imager.

  View details for DOI 10.1118/1.3598435

  View details for Web of Science ID 000292521100037

  View details for PubMedID 21859022

• GPU-based fast low-dose cone beam CT reconstruction via total variation JOURNAL OF X-RAY SCIENCE AND TECHNOLOGY Jia, X., Lou, Y., Lewis, J., Li, R., Gu, X., Men, C., Song, W. Y., Jiang, S. B. 2011; 19 (2): 139-154

  Abstract

  X-ray imaging dose from serial Cone-beam CT (CBCT) scans raises a clinical concern in most image guided radiation therapy procedures. The goal of this paper is to develop a fast GPU-based algorithm to reconstruct high quality CBCT images from undersampled and noisy projection data so as to lower the imaging dose. The CBCT is reconstructed by minimizing an energy functional consisting of a data fidelity term and a total variation regularization term. We develop a GPU-friendly version of a forward-backward splitting algorithm to solve this problem. A multi-grid technique is also employed. We test our CBCT reconstruction algorithm on a digital phantom and a head-and-neck patient case. The performance under low mAs is also validated using physical phantoms. It is found that 40 x-ray projections are sufficient to reconstruct CBCT images with satisfactory quality for clinical purposes. Phantom experiments indicate that CBCT images can be successfully reconstructed under 0.1 mAs/projection. Comparing with the widely used head-and-neck scanning protocol of about 360 projections with 0.4 mAs/projection, an overall 36 times dose reduction has been achieved. The reconstruction time is about 130 sec on an NVIDIA Tesla C1060 GPU card, which is estimated ? 100 times faster than similar regularized iterative reconstruction approaches.

  View details for DOI 10.3233/XST-2011-0283

  View details for Web of Science ID 000292735700001

  View details for PubMedID 21606579

• 3D tumor localization through real-time volumetric x-ray imaging for lung cancer radiotherapy Med. Phys. Li R, Lewis JH, Jia X, Gu X, Folkerts M, Men C, Song WY, Jiang SB 2011; 38 (5): 2783-2794
• Patient-specific motion artifacts in 4DCT MEDICAL PHYSICS Watkins, W. T., Li, R., Lewis, J., Park, J. C., Sandhu, A., Jiang, S. B., Song, W. Y. 2010; 37 (6): 2855-2861

  Abstract

  Four-dimensional computed tomography (4DCT) has enhanced images of the thorax and upper abdomen during respiration, but intraphase residual motion artifacts will persist in cine-mode scanning. In this study, the source and magnitude of projection artifacts due to intraphase target motion is investigated.A theoretical model of geometric uncertainty due to partial projection artifacts in cine-mode 4DCT was derived based on ideal periodic motion. Predicted artifacts were compared to measured errors with a rigid lung phantom attached to a programmable motion platform. Ideal periodic motion and actual patient breathing patterns were used as input for phantom motion. Reconstructed target dimensions were measured along the direction of motion and compared to the actual, known dimensions.Artifacts due to intraphase residual motion in cine-mode 4DCT range from a few mm up to a few cm on a given scanner, and can be predicted based on target motion and CT gantry rotation time. Errors in ITV and GTV dimensions were accurately characterized by the theoretical uncertainty at all phases when sinusoidal motion was considered, and in 96% of 300 measurements when patient breathing patterns were used as motion input. When peak-to-peak motion of 1.5 cm is combined with a breathing period of 4 s and gantry rotation time of 1 s, errors due to partial projection artifacts can be greater than 1 cm near midventilation and are a few mm in the inhale and exhale phases. Incorporation of such uncertainty into margin design should be considered in addition to other uncertainties.Artifacts due to intraphase residual motion exist in 4DCT, even for ideal breathing motions (e.g., sine waves). It was determined that these motion artifacts depend on patient-specific tumor motion and CT gantry rotation speed. Thus, if the patient-specific motion parameters are known (i.e., amplitude and period), a patient-specific margin can and should be designed to compensate for this uncertainty.

  View details for DOI 10.1118/1.3432615

  View details for Web of Science ID 000278573100049

  View details for PubMedID 20632597

• Markerless lung tumor tracking and trajectory reconstruction using rotational cone-beam projections: a feasibility study PHYSICS IN MEDICINE AND BIOLOGY Lewis, J. H., Li, R., Watkins, W. T., Lawson, J. D., Segars, W. P., Cervino, L. I., Song, W. Y., Jiang, S. B. 2010; 55 (9): 2505-2522

  Abstract

  Algorithms for direct tumor tracking in rotational cone-beam projections and for reconstruction of phase-binned 3D tumor trajectories were developed. The feasibility of the algorithm was demonstrated on a digital phantom, a physical phantom and two patients. Tracking results were obtained by comparing reference templates generated from 4DCT to rotational cone-beam projections. The 95th percentile absolute errors (e(95)) in phantom tracking results did not exceed 1.7 mm in either imager dimension, while e(95) in the patients was 3.3 mm or less. Accurate phase-binned trajectories were reconstructed in each case, with 3D maximum errors of no more than 1.0 mm in the phantoms and 2.0 mm in the patients. This work shows the feasibility of a direct tumor tracking technique for rotational images, and demonstrates that an accurate 3D tumor trajectory can be reconstructed from relatively less accurate tracking results. The ability to reconstruct the tumor's average trajectory from a 3D cone-beam CT scan on the day of treatment could allow for better patient setup and quality assurance, while direct tumor tracking in rotational projections could be clinically useful for rotational therapy such as volumetric modulated arc therapy (VMAT).

  View details for DOI 10.1088/0031-9155/55/9/006

  View details for Web of Science ID 000276816400006

  View details for PubMedID 20393232

• GPU-based fast cone beam CT reconstruction from undersampled and noisy projection data via total variation MEDICAL PHYSICS Jia, X., Lou, Y., Li, R., Song, W. Y., Jiang, S. B. 2010; 37 (4): 1757-1760

  Abstract

  Cone-beam CT (CBCT) plays an important role in image guided radiation therapy (IGRT). However, the large radiation dose from serial CBCT scans in most IGRT procedures raises a clinical concern, especially for pediatric patients who are essentially excluded from receiving IGRT for this reason. The goal of this work is to develop a fast GPU-based algorithm to reconstruct CBCT from undersampled and noisy projection data so as to lower the imaging dose.The CBCT is reconstructed by minimizing an energy functional consisting of a data fidelity term and a total variation regularization term. The authors developed a GPU-friendly version of the forward-backward splitting algorithm to solve this model. A multigrid technique is also employed.It is found that 20-40 x-ray projections are sufficient to reconstruct images with satisfactory quality for IGRT. The reconstruction time ranges from 77 to 130 s on an NVIDIA Tesla C1060 (NVIDIA, Santa Clara, CA) GPU card, depending on the number of projections used, which is estimated about 100 times faster than similar iterative reconstruction approaches. Moreover, phantom studies indicate that the algorithm enables the CBCT to be reconstructed under a scanning protocol with as low as 0.1 mA s/projection. Comparing with currently widely used full-fan head and neck scanning protocol of approximately 360 projections with 0.4 mA s/projection, it is estimated that an overall 36-72 times dose reduction has been achieved in our fast CBCT reconstruction algorithm.This work indicates that the developed GPU-based CBCT reconstruction algorithm is capable of lowering imaging dose considerably. The high computation efficiency in this algorithm makes the iterative CBCT reconstruction approach applicable in real clinical environments.

  View details for DOI 10.1118/1.3371691

  View details for Web of Science ID 000276211200040

  View details for PubMedID 20443497

• Single-projection based volumetric image reconstruction and 3D tumor localization in real time for lung cancer radiotherapy. Medical image computing and computer-assisted intervention : MICCAI ... International Conference on Medical Image Computing and Computer-Assisted Intervention Li, R., Jia, X., Lewis, J. H., Gu, X., Folkerts, M., Men, C., Jiang, S. B. 2010; 13: 449-456

  Abstract

  We have developed an algorithm for real-time volumetric image reconstruction and 3D tumor localization based on a single x-ray projection image. We first parameterize the deformation vector fields (DVF) of lung motion by principal component analysis (PCA). Then we optimize the DVF applied to a reference image by adapting the PCA coefficients such that the simulated projection of the reconstructed image matches the measured projection. The algorithm was tested on a digital phantom as well as patient data. The average relative image reconstruction error and 3D tumor localization error for the phantom is 7.5% and 0.9 mm, respectively. The tumor localization error for patient is approximately 2 mm. The computation time of reconstructing one volumetric image from each projection is around 0.2 and 0.3 seconds for phantom and patient, respectively, on an NVIDIA C1060 GPU. Clinical application can potentially lead to accurate 3D tumor tracking from a single imager.

  View details for PubMedID 20879431

• A feasibility study of markerless fluoroscopic gating for lung cancer radiotherapy using 4DCT templates PHYSICS IN MEDICINE AND BIOLOGY Li, R., Lewis, J. H., Cervino, L. I., Jiang, S. B. 2009; 54 (20): N489-N500

  Abstract

  A major difficulty in conformal lung cancer radiotherapy is respiratory organ motion, which may cause clinically significant targeting errors. Respiratory-gated radiotherapy allows for more precise delivery of prescribed radiation dose to the tumor, while minimizing normal tissue complications. Gating based on external surrogates is limited by its lack of accuracy, while gating based on implanted fiducial markers is limited primarily by the risk of pneumothorax due to marker implantation. Techniques for fluoroscopic gating without implanted fiducial markers (markerless gating) have been developed. These techniques usually require a training fluoroscopic image dataset with marked tumor positions in the images, which limits their clinical implementation. To remove this requirement, this study presents a markerless fluoroscopic gating algorithm based on 4DCT templates. To generate gating signals, we explored the application of three similarity measures or scores between fluoroscopic images and the reference 4DCT template: un-normalized cross-correlation (CC), normalized cross-correlation (NCC) and normalized mutual information (NMI), as well as average intensity (AI) of the region of interest (ROI) in the fluoroscopic images. Performance was evaluated using fluoroscopic and 4DCT data from three lung cancer patients. On average, gating based on CC achieves the highest treatment accuracy given the same efficiency, with a high target coverage (average between 91.9% and 98.6%) for a wide range of nominal duty cycles (20-50%). AI works well for two patients out of three, but failed for the third patient due to interference from the heart. Gating based on NCC and NMI usually failed below 50% nominal duty cycle. Based on this preliminary study with three patients, we found that the proposed CC-based gating algorithm can generate accurate and robust gating signals when using 4DCT reference template. However, this observation is based on results obtained from a very limited dataset, and further investigation on a larger patient population has to be done before its clinical implementation.

  View details for DOI 10.1088/0031-9155/54/20/N03

  View details for Web of Science ID 000270563300028

  View details for PubMedID 19779221

• 4D CT sorting based on patient internal anatomy PHYSICS IN MEDICINE AND BIOLOGY Li, R., Lewis, J. H., Cervino, L. I., Jiang, S. B. 2009; 54 (15): 4821-4833

  Abstract

  Respiratory motion during free-breathing computed tomography (CT) scan may cause significant errors in target definition for tumors in the thorax and upper abdomen. A four-dimensional (4D) CT technique has been widely used for treatment simulation of thoracic and abdominal cancer radiotherapy. The current 4D CT techniques require retrospective sorting of the reconstructed CT slices oversampled at the same couch position. Most sorting methods depend on external surrogates of respiratory motion recorded by extra instruments. However, respiratory signals obtained from these external surrogates may not always accurately represent the internal target motion, especially when irregular breathing patterns occur. We have proposed a new sorting method based on multiple internal anatomical features for multi-slice CT scan acquired in the cine mode. Four features are analyzed in this study, including the air content, lung area, lung density and body area. We use a measure called spatial coherence to select the optimal internal feature at each couch position and to generate the respiratory signals for 4D CT sorting. The proposed method has been evaluated for ten cancer patients (eight with thoracic cancer and two with abdominal cancer). For nine patients, the respiratory signals generated from the combined internal features are well correlated to those from external surrogates recorded by the real-time position management (RPM) system (average correlation: 0.95+/-0.02), which is better than any individual internal measures at 95% confidence level. For these nine patients, the 4D CT images sorted by the combined internal features are almost identical to those sorted by the RPM signal. For one patient with an irregular breathing pattern, the respiratory signals given by the combined internal features do not correlate well with those from RPM (correlation: 0.68+/-0.42). In this case, the 4D CT image sorted by our method presents fewer artifacts than that from the RPM signal. Our 4D CT internal sorting method eliminates the need of externally recorded surrogates of respiratory motion. It is an automatic, accurate, robust, cost efficient and yet simple method and therefore can be readily implemented in clinical settings.

  View details for DOI 10.1088/0031-9155/54/15/012

  View details for Web of Science ID 000268191900012

  View details for PubMedID 19622855

• Single-trial P300 estimation with a spatiotemporal filtering method JOURNAL OF NEUROSCIENCE METHODS Li, R., Keil, A., Principe, J. C. 2009; 177 (2): 488-496

  Abstract

  A spatiotemporal filtering method for single-trial ERP component estimation is presented. Instead of modeling the entire ERP waveform, the method focuses on the ERP component local descriptors (amplitude and latency) thru the spatial diversity of multichannel recordings and thus it is tailored to extract signals in negative signal to noise ratio conditions. The model allows for both amplitude and latency variability in the ERP component under investigation. We applied the method to the estimation of the P300 component in an oddball target detection task and found that negative correlations exist between response time and single-trial P300 amplitude.

  View details for DOI 10.1016/j.jneumeth.2008.10.035

  View details for Web of Science ID 000263393300029

  View details for PubMedID 19041343

• Markerless gating for lung cancer radiotherapy based on machine learning techniques PHYSICS IN MEDICINE AND BIOLOGY Lin, T., Li, R., Tang, X., Dy, J. G., Jiang, S. B. 2009; 54 (6): 1555-1563

  Abstract

  In lung cancer radiotherapy, radiation to a mobile target can be delivered by respiratory gating, for which we need to know whether the target is inside or outside a predefined gating window at any time point during the treatment. This can be achieved by tracking one or more fiducial markers implanted inside or near the target, either fluoroscopically or electromagnetically. However, the clinical implementation of marker tracking is limited for lung cancer radiotherapy mainly due to the risk of pneumothorax. Therefore, gating without implanted fiducial markers is a promising clinical direction. We have developed several template-matching methods for fluoroscopic marker-less gating. Recently, we have modeled the gating problem as a binary pattern classification problem, in which principal component analysis (PCA) and support vector machine (SVM) are combined to perform the classification task. Following the same framework, we investigated different combinations of dimensionality reduction techniques (PCA and four nonlinear manifold learning methods) and two machine learning classification methods (artificial neural networks-ANN and SVM). Performance was evaluated on ten fluoroscopic image sequences of nine lung cancer patients. We found that among all combinations of dimensionality reduction techniques and classification methods, PCA combined with either ANN or SVM achieved a better performance than the other nonlinear manifold learning methods. ANN when combined with PCA achieves a better performance than SVM in terms of classification accuracy and recall rate, although the target coverage is similar for the two classification methods. Furthermore, the running time for both ANN and SVM with PCA is within tolerance for real-time applications. Overall, ANN combined with PCA is a better candidate than other combinations we investigated in this work for real-time gated radiotherapy.

  View details for DOI 10.1088/0031-9155/54/6/010

  View details for Web of Science ID 000263903600011

  View details for PubMedID 19229098

• A Spatiotemporal Filtering Methodology for Single-Trial ERP Component Estimation IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING Li, R., Principe, J. C., Bradley, M., Ferrari, V. 2009; 56 (1): 83-92

  Abstract

  A new spatiotemporal filtering method for single-trial event-related potential (ERP) estimation is proposed. Instead of attempting to model the entire ERP waveform, the method relies on modeling ERP component descriptors (amplitude and latency) thru the spatial diversity of multichannel recordings, and thus, it is tailored to extract signals in negative SNR conditions. The model allows for both amplitude and latency variability in the ERP component under investigation. The extracted ERP component is constrained through a spatial filter to have minimal distance (with respect to some metric) in the temporal domain from a user-designed template component. The spatial filter may be interpreted as a noise canceller in the spatial domain. Study with both simulated data and real cognitive ERP data shows the effectiveness of the proposed method.

  View details for DOI 10.1109/TBME.2008.2002153

  View details for Web of Science ID 000263640900011

  View details for PubMedID 19224722

• A unifying criterion for instantaneous blind source separation Signal Processing Li R, Liu W, Principe JC 2007; 87 (8): 1872-1881