Search

Objective. To demonstrate that complete cone beam CT (CBCT) scans from both MV-energy and kV-energy LINAC sources can reduce metal artifacts in radiotherapy guidance, while maintaining standard-of-care x-ray doses levels. Approach. MV-CBCT and kV-CBCT scans are acquired at half normal dose. The impact of lowered dose on MV-CBCT data quality is mitigated by the use of a 4-layer MV-imager prototype and reduced LINAC energy settings (2.5 MV) to improve photon capture. Additionally, the MV-CBCT is used to determine the 3D position and pose of metal implants, which in turn is used to guide model-based poly-energetic correction and interleaving of the kV-CBCT and MV-CBCT data. Certain edge-preserving regularization steps incorporated into the model-based correction algorithm further reduce MV data noise. Main results. The method was tested in digital phantoms and a real pelvis phantom with large 2.5″ spherical inserts, emulating hip replacements of different materials. The proposed method demonstrated an appealing compromise between the high contrast of kV-CBCT and low artifact content of MV-CBCT. Contrast-to-noise improved 3-fold compared to MV-CBCT with a clinical 1-layer architecture at matched dose (37 mGy) and edge blur levels. Visual delineation of the bladder and prostate improved noteably over kV- or MV-CBCT alone. Significance. The proposed method demonstrates that a full MV-CBCT scan can be combined with kV-CBCT to reduce metal artifacts without resorting to complicated beam collimation strategies to limit the MV-CBCT dose contribution. Additionally, significant improvements in CNR can be achieved as compared to metal artifact reduction through current clinical MV-CBCT practices. [ABSTRACT FROM AUTHOR]

Objective. CBCTs in image-guided radiotherapy provide crucial anatomy information for patient setup and plan evaluation. Longitudinal CBCT image registration could quantify the inter-fractional anatomic changes, e.g. tumor shrinkage, and daily OAR variation throughout the course of treatment. The purpose of this study is to propose an unsupervised deep learning-based CBCT-CBCT deformable image registration which enables quantitative anatomic variation analysis. Approach. The proposed deformable registration workflow consists of training and inference stages that share the same feed-forward path through a spatial transformation-based network (STN). The STN consists of a global generative adversarial network (GlobalGAN) and a local GAN (LocalGAN) to predict the coarse- and fine-scale motions, respectively. The network was trained by minimizing the image similarity loss and the deformable vector field (DVF) regularization loss without the supervision of ground truth DVFs. During the inference stage, patches of local DVF were predicted by the trained LocalGAN and fused to form a whole-image DVF. The local whole-image DVF was subsequently combined with the GlobalGAN generated DVF to obtain the final DVF. The proposed method was evaluated using 100 fractional CBCTs from 20 abdominal cancer patients in the experiments and 105 fractional CBCTs from a cohort of 21 different abdominal cancer patients in a holdout test. Main Results. Qualitatively, the registration results show good alignment between the deformed CBCT images and the target CBCT image. Quantitatively, the average target registration error calculated on the fiducial markers and manually identified landmarks was 1.91 ± 1.18 mm. The average mean absolute error, normalized cross correlation between the deformed CBCT and target CBCT were 33.42 ± 7.48 HU, 0.94 ± 0.04, respectively. Significance. In summary, an unsupervised deep learning-based CBCT-CBCT registration method is proposed and its feasibility and performance in fractionated image-guided radiotherapy is investigated. This promising registration method could provide fast and accurate longitudinal CBCT alignment to facilitate inter-fractional anatomic changes analysis and prediction. [ABSTRACT FROM AUTHOR]

Objective. In this study we developed an automatic method to predict tumour volume and shape in weeks 3 and 4 of radiotherapy (RT), using cone-beam computed tomography (CBCT) scans acquired up to week 2, allowing identification of large tumour changes. Approach. 240 non-small cell lung cancer (NSCLC) patients, treated with 55 Gy in 20 fractions, were collected. CBCTs were rigidly registered to the planning CT. Intensity values were extracted in each voxel of the planning target volume across all CBCT images from days 1, 2, 3, 7 and 14. For each patient and in each voxel, four regression models were fitted to voxel intensity; applying linear, Gaussian, quadratic and cubic methods. These models predicted the intensity value for each voxel in weeks 3 and 4, and the tumour volume found by thresholding. Each model was evaluated by computing the root mean square error in pixel value and structural similarity index metric (SSIM) for all patients. Finally, the sensitivity and specificity to predict a 30% change in volume were calculated for each model. Main results. The linear, Gaussian, quadratic and cubic models achieved a comparable similarity score, the average SSIM for all patients was 0.94, 0.94, 0.90, 0.83 in week 3, respectively. At week 3, a sensitivity of 84%, 53%, 90% and 88%, and specificity of 99%, 100%, 91% and 42% were observed for the linear, Gaussian, quadratic and cubic models respectively. Overall, the linear model performed best at predicting those patients that will benefit from RT adaptation. The linear model identified 21% and 23% of patients in our cohort with more than 30% tumour volume reduction to benefit from treatment adaptation in weeks 3 and 4 respectively. Significance. We have shown that it is feasible to predict the shape and volume of NSCLC tumours from routine CBCTs and effectively identify patients who will respond to treatment early. [ABSTRACT FROM AUTHOR]

This study aimed to investigate if a commercial, knowledge-based tool for radiotherapy planning could be used to estimate the amount of sparing in organs at risk (OARs) in the re-planning strategy for adaptive radiotherapy (ART). Eighty head and neck (HN) VMAT Pareto plans from our institute's database were used to train a knowledge-based planning (KBP) model. An evaluation set of another 20 HN patients was randomly selected. For each patient in the evaluation set, the planning computed tomography (CT) and 2 sets of on-board cone-beam CT, corresponding to the middle and second half of the radiotherapy treatment course, were extracted. The original plan was re-calculated on a daily deformed CT (delivered dose‐volume histogram (DVH)) and compared with the KBP DVH predictions and with the final KBP DVH after optimisation of the plan, which was performed on the same image sets. To evaluate the feasibility of this method, the range of KBP DVH uncertainties was compared with the gains obtained from re-planning. DVH differences and receiver operating characteristic (ROC) curve analysis were used for this purpose. On average, final KBP uncertainties were smaller than the gain in re-planning. Statistical tests confirmed significant differences between the two groups. ROC analysis showed KBP performance in terms of area under the curve values higher than 0.7, which confirmed a good accuracy in predicted values. Overall, for 48% of cases, KBP predicted a desirable outcome from re-planning, and the final dose confirmed an effective gain in 47% of cases. We have established a systematic workflow to identify effective OAR sparing in re-planning based on KBP predictions that can be implemented in an on-line, ART process. [ABSTRACT FROM AUTHOR]

Purpose. A radiotherapy system with a fixed treatment beam and a rotating patient positioning system could be smaller, more robust and more cost effective compared to conventional rotating gantry systems. However, patient rotation could cause anatomical deformation and compromise treatment delivery. In this work, we demonstrate an image-guided treatment workflow with a fixed beam prototype system that accounts for deformation during rotation to maintain dosimetric accuracy. Methods. The prototype system consists of an Elekta Synergy linac with the therapy beam orientated downward and a custom-built patient rotation system (PRS). A phantom that deforms with rotation was constructed and rotated within the PRS to quantify the performance of two image guidance techniques: motion compensated cone-beam CT (CBCT) for pre-treatment volumetric imaging and kilovoltage infraction monitoring (KIM) for real-time image guidance. The phantom was irradiated with a 3D conformal beam to evaluate the dosimetric accuracy of the workflow. Results. The motion compensated CBCT was used to verify pre-treatment position and the average calculated position was within −0.3 ± 1.1 mm of the phantom's ground truth position at 0°. KIM tracked the position of the target in real-time as the phantom was rotated and the average calculated position was within −0.2 ± 0.8 mm of the phantom's ground truth position. A 3D conformal treatment delivered on the prototype system with image guidance had a 3%/2 mm gamma pass rate of 96.3% compared to 98.6% delivered using a conventional rotating gantry linac. Conclusions. In this work, we have shown that image guidance can be used with fixed-beam treatment systems to measure and account for changes in target position in order to maintain dosimetric coverage during horizontal rotation. This treatment modality could provide a viable treatment option when there insufficient space for a conventional linear accelerator or where the cost is prohibitive. [ABSTRACT FROM AUTHOR]

Cone-beam computed tomography (CBCT) images acquired during radiotherapy may allow early response assessment. Previous studies have reported inconsistent findings on an association of CBCT-measured tumor volume changes with clinical outcomes. The purpose of this pilot study was twofold: (1) to characterize changes in CBCT-based radiomics features during treatment; and (2) to quantify the potential association of CBCT-based delta-radiomics features with overall survival in locally advanced lung cancer. We retrospectively identified 23 patients and calculated 658 radiomics features from each of 11 CBCT images per patient. Feature selection was performed based on repeatability, robustness against contouring uncertainties, and non-redundancy. We calculated the coefficient of determination (R2) for the relationship between the actual feature value at the end of treatment and predicted value based on linear models fitted using features between the first and kth fractions. We also quantified the predictive ability for survival with two methods by: (1) comparing delta-radiomics features (defined as the mean change between the first and kth fractions) between two groups of patients divided by a cutoff survival time of 18 months using the t-test or Wilcoxon rank-sum test; and (2) quantifying univariate discrimination of two groups divided by the median of delta-radiomics feature. All selected seven radiomics features during treatment (as early as the 10th fraction) were predictive of those at the end of treatment (R2 > 0.64). Three delta-radiomics features demonstrated significant differences (q < 0.05, as early as the 10th fraction) between the two groups of patients divided by the cutoff survival time. Two of those three features were also predictive of survival according to the log-rank statistics. We provided the first demonstration of a potential association of CBCT-based delta-radiomics features early during treatment with overall survival in locally advanced lung cancer. Our preliminary findings should be validated for a larger cohort of patients. [ABSTRACT FROM AUTHOR]

This study aims to utilize a deep convolutional neural network (DCNN) for synthesized CT image generation based on cone-beam CT (CBCT) and to apply the images to dose calculations for nasopharyngeal carcinoma (NPC). An encoder-decoder 2D U-Net neural network was produced. A total of 70 CBCT/CT paired images of NPC cancer patients were used for training (50), validation (10) and testing (10) datasets. The testing datasets were treated with the same prescription dose (70 Gy to PTVnx70, 68 Gy to PTVnd68, 62 Gy to the PTV62 and 54 Gy to the PTV54). The mean error (ME) and mean absolute error (MAE) for the true CT images were calculated for image quality evaluation of the synthesized CT. The dose-volume histogram (DVH) dose metric difference and 3D gamma pass rate for the true CT images were calculated for dose analysis, and the results were compared with those for the CBCT images (original CBCT images without any correction) and a patient-specific calibration (PSC) method. Compared with CBCT, the range of the MAE for synthesized CT images improved from (60, 120) to (6, 27) Hounsfield units (HU), and the ME improved from (−74, 51) to (−26, 4) HU. Compared with the true CT method, the average DVH dose metric differences for the CBCT, PSC and synthesized CT methods were 0.8%  ±  1.9%, 0.4%  ±  0.7% and 0.2%  ±  0.6%, respectively. The 1%/1 mm gamma pass rates within the body for the CBCT, PSC and synthesized CT methods were 90.8%  ±  6.2%, 94.1%  ±  4.4% and 95.5%  ±  1.6%, respectively, and the rates within the PTVnx70 were 80.3%  ±  16.6%, 87.9%  ±  19.7%, 98.6%  ±  2.9%, respectively. The DCNN model can generate high-quality synthesized CT images from CBCT images and be used for accurate dose calculations for NPC patients. This finding has great significance for the clinical application of adaptive radiotherapy for NPC. [ABSTRACT FROM AUTHOR]

Throughout the course of delivering a radiation therapy treatment, which may take several weeks, a patient’s anatomy may change drastically, and adaptive radiation therapy (ART) may be needed. Cone-beam computed tomography (CBCT), which is often available during the treatment process, can be used for both patient positioning and ART re-planning. However, due to the prominent amount of noise, artifacts, and inaccurate Hounsfield unit (HU) values, the dose calculation based on CBCT images could be inaccurate for treatment planning. One way to solve this problem is to convert CBCT images to more accurate synthesized CT (sCT) images. In this work, we have developed a cycle-consistent generative adversarial network framework (CycleGAN) to synthesize CT images from CBCT images. This model is capable of image-to-image translation using unpaired CT and CBCT images in an unsupervised learning setting. The sCT images generated from CBCT through this CycleGAN model are visually and quantitatively similar to real CT images with decreased mean absolute error (MAE) from 69.29 HU to 29.85 HU for head-and-neck (H&N) cancer patients. The dose distributions calculated on the sCT by CycleGAN demonstrated a higher accuracy than those on CBCT in a 3D gamma index analysis with increased gamma index pass rate from 86.92% to 96.26% under 1 mm/1% criteria, when using the deformed planning CT image (dpCT) as the reference. We also compared the CycleGAN model with other unsupervised learning methods, including deep convolutional generative adversarial networks (DCGAN) and progressive growing of GANs (PGGAN), and demonstrated that CycleGAN outperformed the other two models. A phantom study has been conducted to compare sCT with dpCT, and the increase of structural similarity index from 0.91 to 0.93 shows that CycleGAN performed better than DIR in terms of preserving anatomical accuracy. [ABSTRACT FROM AUTHOR]

Radiation-associated toxicities due to sophisticated ocular anatomy and shape variability of organs at risk (OARs) are major concerns during external beam radiation therapy (EBRT) of patients with intraocular cancer. A novel two-step image registration strategy between optical coherence tomography (OCT) and computed tomography (CT) images was proposed and validated to precisely localize the target in the EBRT of patients with intraocular cancer. Specifically, multiple features from OCT and CT images were extracted automatically, then spatial transformation based on thin-plate spline function was performed iteratively to achieve feature alignment between the CT and OCT images. Finally, an exclusive OR (XOR) algorithm was applied for precise 3D registration using a 3D-mesh model generated from OCT and CT volumes. The accuracy of the proposed novel registration strategy was validated and tested in a schematic-eye phantom with an artificially introduced tumor and in ten patients with confirmed primary and/or secondary intraocular cancer. There was an average registration error and computational time of 0.21  ±  0.05° and 259  ±  5 s, together with an average Dice similarity coefficient and Hausdorff distance of 88.4  ±  0.65 and 0.89  ±  0.09, respectively. The preliminary experimental results demonstrated that the proposed novel strategy to overcome current limitations on eye modeling and to localize precisely the tumor target during EBRT of intraocular cancer is promising. [ABSTRACT FROM AUTHOR]

Major advance was done in preclinical radiotherapy thanks to the development of image guided micro-irradiator. Nevertheless, some applications still can benefit of improvements, such as the irradiation of mobile tumors. This preclinical radiotherapy case presents increased difficulties compared to clinical practice because of the waveform of small animals breathing cycle, its frequency and amplitude. To answer this issue, we developed a specific beam shutter and implemented respiratory gating on the X-RAD 225Cx preclinical irradiator. In the first step of this study, the shutter was accurately characterized. Opening and closing speed of 1.28 and 0.33 mm ms−1 were respectively measured, and a transmission of 0.7% of the beam was measured with the shutter fully closed. Beam-on times were also determined for various gating parameters and highlighted a difference of 57 ms between the beam delivery duration and the gate width. This discrepancy was compensated during the respiratory monitoring adjustment. In a second step, a respiratory protocol was evaluated with two vertical beams of 2.5 and 5 mm diameters, for motion amplitudes ranging from 0.5 to 4 mm. This evaluation demonstrated the effectiveness of our set up to perform motion compensation for amplitude as small as 0.5 mm despite a dose gradient of 1.47 cGy mm−1 observed with the 5 mm irradiation field, due to the shutter opening and closing durations. We also investigated the efficiency of a scintillating fiber dosimeter, adapted to small beams and providing real-time dose rate measurements. This detector showed very good performances to detect motion in small irradiation fields and would be very suitable to monitor the number of delivered gates until the planned delivered dose is achieved. This study presented a new respiratory gating set up and showed that very efficient motion compensation could be achieved in small animal radiotherapy. [ABSTRACT FROM AUTHOR]

Radiotherapy dose calculations based on cone-beam CT (CBCT) images can be inaccurate due to unreliable Hounsfield units (HU) in the CBCT. Deformable image registration of planning CT images to CBCT, and direct correction of CBCT image values are two methods proposed to allow heterogeneity corrected dose calculations based on CBCT. In this paper we compare the accuracy and robustness of these two approaches. CBCT images for 44 patients were used including pelvis, lung and head & neck sites. CBCT HU were corrected using a ‘shading correction’ algorithm and via deformable registration of planning CT to CBCT using either Elastix or Niftyreg. Radiotherapy dose distributions were re-calculated with heterogeneity correction based on the corrected CBCT and several relevant dose metrics for target and OAR volumes were calculated. Accuracy of CBCT based dose metrics was determined using an ‘override ratio’ method where the ratio of the dose metric to that calculated on a bulk-density assigned version of the same image is assumed to be constant for each patient, allowing comparison to the patient’s planning CT as a gold standard. Similar performance is achieved by shading corrected CBCT and both deformable registration algorithms, with mean and standard deviation of dose metric error less than 1% for all sites studied. For lung images, use of deformed CT leads to slightly larger standard deviation of dose metric error than shading corrected CBCT with more dose metric errors greater than 2% observed (7% versus 1%). [ABSTRACT FROM AUTHOR]

A systematic method is presented for determining optimal parameters in variable-kernel deformable image registration of cone beam CT and CT images, in order to improve accuracy and convergence for potential use in online adaptive radiotherapy. Assessed conditions included the noise constant (symmetric force demons), the kernel reduction rate, the kernel reduction percentage, and the kernel adjustment criteria. Four such parameters were tested in conjunction with reductions of 5, 10, 15, 20, 30, and 40%. Noise constants ranged from 1.0 to 1.9 for pelvic images in ten prostate cancer patients. A total of 516 tests were performed and assessed using the structural similarity index. Registration accuracy was plotted as a function of iteration number and a least-squares regression line was calculated, which implied an average improvement of 0.0236% per iteration. This baseline was used to determine if a given set of parameters under- or over-performed. The most accurate parameters within this range were applied to contoured images. The mean Dice similarity coefficient was calculated for bladder, prostate, and rectum with mean values of 98.26%, 97.58%, and 96.73%, respectively; corresponding to improvements of 2.3%, 9.8%, and 1.2% over previously reported values for the same organ contours. This graphical approach to registration analysis could aid in determining optimal parameters for Demons-based algorithms. It also establishes expectation values for convergence rates and could serve as an indicator of non-physical warping, which often occurred in cases  >0.6% from the regression line. [ABSTRACT FROM AUTHOR]

Four-dimensional cone-beam computed tomography (4D-CBCT) of the free-breathing thorax is a valuable tool in image-guided radiation therapy of the thorax and the upper abdomen. It allows the determination of the position of a tumor throughout the breathing cycle, while only its mean position can be extracted from three-dimensional CBCT. The classical approaches are not fully satisfactory: respiration-correlated methods allow one to accurately locate high-contrast structures in any frame, but contain strong streak artifacts unless the acquisition is significantly slowed down. Motion-compensated methods can yield streak-free, but static, reconstructions. This work proposes a 4D-CBCT method that can be seen as a trade-off between respiration-correlated and motion-compensated reconstruction. It builds upon the existing reconstruction using spatial and temporal regularization (ROOSTER) and is called motion-aware ROOSTER (MA-ROOSTER). It performs temporal regularization along curved trajectories, following the motion estimated on a prior 4D CT scan. MA-ROOSTER does not involve motion-compensated forward and back projections: the input motion is used only during temporal regularization. MA-ROOSTER is compared to ROOSTER, motion-compensated Feldkamp–Davis–Kress (MC-FDK), and two respiration-correlated methods, on CBCT acquisitions of one physical phantom and two patients. It yields streak-free reconstructions, visually similar to MC-FDK, and robust information on tumor location throughout the breathing cycle. MA-ROOSTER also allows a variation of the lung tissue density during the breathing cycle, similar to that of planning CT, which is required for quantitative post-processing. [ABSTRACT FROM AUTHOR]

We present a study of extracting respiratory signals from cone beam computed tomography (CBCT) projections within the framework of the Amsterdam Shroud (AS) technique. Acquired prior to the radiotherapy treatment, CBCT projections were preprocessed for contrast enhancement by converting the original intensity images to attenuation images with which the AS image was created. An adaptive robust z-normalization filtering was applied to further augment the weak oscillating structures locally. From the enhanced AS image, the respiratory signal was extracted using a two-step optimization approach to effectively reveal the large-scale regularity of the breathing signals. CBCT projection images from five patients acquired with the Varian Onboard Imager on the Clinac iX System Linear Accelerator (Varian Medical Systems, Palo Alto, CA) were employed to assess the proposed technique. Stable breathing signals can be reliably extracted using the proposed algorithm. Reference waveforms obtained using an air bellows belt (Philips Medical Systems, Cleveland, OH) were exported and compared to those with the AS based signals. The average errors for the enrolled patients between the estimated breath per minute (bpm) and the reference waveform bpm can be as low as  −0.07 with the standard deviation 1.58. The new algorithm outperformed the original AS technique for all patients by 8.5% to 30%. The impact of gantry rotation on the breathing signal was assessed with data acquired with a Quasar phantom (Modus Medical Devices Inc., London, Canada) and found to be minimal on the signal frequency. The new technique developed in this work will provide a practical solution to rendering markerless breathing signal using the CBCT projections for thoracic and abdominal patients. [ABSTRACT FROM AUTHOR]

Optical cone-beam computed tomographic (CBCT) scanning of 3D radiochromic dosimeters may provide a practical method for 3D dose verification in radiation therapy. However, in cone-beam geometry stray light contaminates the projection images, degrading the accuracy of reconstructed linear attenuation coefficients. Stray light was measured using a beam pass aperture array (BPA) and structured illumination methods. The stray-to-primary ray ratio (SPR) along the central axis was found to be 0.24 for a 5% gelatin hydrogel, representative of radiochromic hydrogels. The scanner was modified by moving the spectral filter from the detector to the source, changing the light’s spatial fluence pattern and lowering the acceptance angle by extending distance between the source and object. These modifications reduced the SPR significantly from 0.24 to 0.06. The accuracy of the reconstructed linear attenuation coefficients for uniform carbon black liquids was compared to independent spectrometer measurements. Reducing the stray light increased the range of accurate transmission readings. In order to evaluate scanner performance for the more challenging application to small field dosimetry, a carbon black finger gel phantom was prepared. Reconstructions of the phantom from CBCT and fan-beam CT scans were compared. The modified source resulted in improved agreement. Subtraction of residual stray light, measured with BPA or structured illumination from each projection further improved agreement. Structured illumination was superior to BPA for measuring stray light for the smaller 1.2 and 0.5 cm diameter phantom fingers. At the costs of doubling the scanner size and tripling the number of scans, CBCT reconstructions of low-scattering hydrogel dosimeters agreed with those of fan-beam CT scans. [ABSTRACT FROM AUTHOR]

Characterization of anatomical change and other differences is important in sequential computed tomography (CT) imaging, where a high-fidelity patient-specific prior image is typically present, but is not used, in the reconstruction of subsequent anatomical states. Here, we introduce a penalized likelihood (PL) method called reconstruction of difference (RoD) to directly reconstruct a difference image volume using both the current projection data and the (unregistered) prior image integrated into the forward model for the measurement data. The algorithm utilizes an alternating minimization to find both the registration and reconstruction estimates. This formulation allows direct control over the image properties of the difference image, permitting regularization strategies that inhibit noise and structural differences due to inconsistencies between the prior image and the current data. Additionally, if the change is known to be local, RoD allows local acquisition and reconstruction, as opposed to traditional model-based approaches that require a full support field of view (or other modifications). We compared the performance of RoD to a standard PL algorithm, in simulation studies and using test-bench cone-beam CT data. The performances of local and global RoD approaches were similar, with local RoD providing a significant computational speedup. In comparison across a range of data with differing fidelity, the local RoD approach consistently showed lower error (with respect to a truth image) than PL in both noisy data and sparsely sampled projection scenarios. In a study of the prior image registration performance of RoD, a clinically reasonable capture ranges were demonstrated. Lastly, the registration algorithm had a broad capture range and the error for reconstruction of CT data was 35% and 20% less than filtered back-projection for RoD and PL, respectively. The RoD has potential for delivering high-quality difference images in a range of sequential clinical scenarios including image-guided surgeries and treatments where accurate and quantitative assessments of anatomical change is desired. [ABSTRACT FROM AUTHOR]

Kilo-voltage (KV) cone-beam computed tomography (CBCT) unit mounted onto a linear accelerator treatment system, often referred to as on-board imager (OBI), plays an increasingly important role in image-guided radiation therapy. While the FDK algorithm is currently used for reconstructing images from clinical OBI data, optimization-based reconstruction has also been investigated for OBI CBCT. An optimization-based reconstruction involves numerous parameters, which can significantly impact reconstruction properties (or utility). The success of an optimization-based reconstruction for a particular class of practical applications thus relies strongly on appropriate selection of parameter values. In the work, we focus on tailoring the constrained-TV-minimization-based reconstruction, an optimization-based reconstruction previously shown of some potential for CBCT imaging conditions of practical interest, to OBI imaging through appropriate selection of parameter values. In particular, for given real data of phantoms and patient collected with OBI CBCT, we first devise utility metrics specific to OBI-quality-assurance tasks and then apply them to guiding the selection of parameter values in constrained-TV-minimization-based reconstruction. The study results show that the reconstructions are with improvement, relative to clinical FDK reconstruction, in both visualization and quantitative assessments in terms of the devised utility metrics. [ABSTRACT FROM AUTHOR]

This study was designed to evaluate contouring variability of human-and deformable-generated contours on planning CT (PCT) and CBCT for ten patients with low-or intermediate-risk prostate cancer. For each patient in this study, five radiation oncologists contoured the prostate, bladder, and rectum, on one PCT dataset and five CBCT datasets. Consensus contours were generated using the STAPLE method in the CERR software package. Observer contours were compared to consensus contour, and contour metrics (Dice coefficient, Hausdorff distance, Contour Distance, Center-of-Mass [COM] Deviation) were calculated. In addition, the first day CBCT was registered to subsequent CBCT fractions (CBCTn: CBCT2–CBCT5) via B-spline Deformable Image Registration (DIR). Contours were transferred from CBCT1 to CBCTn via the deformation field, and contour metrics were calculated through comparison with consensus contours generated from human contour set. The average contour metrics for prostate contours on PCT and CBCT were as follows: Dice coefficient—0.892 (PCT), 0.872 (CBCT-Human), 0.824 (CBCT-Deformed); Hausdorff distance—4.75 mm (PCT), 5.22 mm (CBCT-Human), 5.94 mm (CBCT-Deformed); Contour Distance (overall contour)—1.41 mm (PCT), 1.66 mm (CBCT-Human), 2.30 mm (CBCT-Deformed); COM Deviation—2.01 mm (PCT), 2.78 mm (CBCT-Human), 3.45 mm (CBCT-Deformed). For human contours on PCT and CBCT, the difference in average Dice coefficient between PCT and CBCT (approx. 2%) and Hausdorff distance (approx. 0.5 mm) was small compared to the variation between observers for each patient (standard deviation in Dice coefficient of 5% and Hausdorff distance of 2.0 mm). However, additional contouring variation was found for the deformable-generated contours (approximately 5.0% decrease in Dice coefficient and 0.7 mm increase in Hausdorff distance relative to human-generated contours on CBCT). Though deformable contours provide a reasonable starting point for contouring on CBCT, we conclude that contours generated with B-Spline DIR require physician review and editing if they are to be used in the clinic. [ABSTRACT FROM AUTHOR]

The purpose of this work is to develop an automated landmark-guided deformable image registration (LDIR) algorithm between the planning CT and daily cone-beam CT (CBCT) with low image quality. This method uses an automated landmark generation algorithm in conjunction with a local small volume gradient matching search engine to map corresponding landmarks between the CBCT and the planning CT. The landmarks act as stabilizing control points in the following Demons deformable image registration. LDIR is implemented on graphics processing units (GPUs) for parallel computation to achieve ultra fast calculation. The accuracy of the LDIR algorithm has been evaluated on a synthetic case in the presence of different noise levels and data of six head and neck cancer patients. The results indicate that LDIR performed better than rigid registration, Demons, and intensity corrected Demons for all similarity metrics used. In conclusion, LDIR achieves high accuracy in the presence of multimodality intensity mismatch and CBCT noise contamination, while simultaneously preserving high computational efficiency. [ABSTRACT FROM AUTHOR]

Onboard cone-beam CT (CBCT) has been widely used in image guided radiation therapy. However, the longitudinal coverage is only 15.5 cm in the pelvis scan mode. As a result, a single CBCT scan cannot cover the planning target volume in the longitudinal direction for over 80% of the patients. The common approach is to use double- or multiple-circular scans and then combine multiple CBCT volumes after reconstruction. However it raises concerns regarding doubled imaging dose at the imaging beam junctions due to beam divergence. In this work, we present a new method, DSCS (Dual Scan with Complementary Shifts), to address the CBCT coverage problem using a pair of complementary circular scans. In DSCS, two circular scans were performed at 39.5 cm apart longitudinally. In the superior scan, the detector panel was offset by 16 cm to the left, 15 cm to the inferior. In the inferior scan, the detector panel was shifted 16 cm to the right and 15 cm to the superior. The effective imaging volume is 39.5 cm longitudinally with a 45 cm lateral field-of-view (FOV). Half beam blocks were used to confine the imaging radiation inside the volume of interest. A new image reconstruction algorithm was developed, based on the Feldkamp–Davis–Kress cone-beam CT reconstruction algorithm, to support the DSCS scanning geometry. Digital phantom simulations were performed to demonstrate the feasibility of DSCS. Physical phantom studies were performed using an anthropomorphic phantom on a commercial onboard CBCT system. With basic scattering corrections, the reconstruction results were acceptable. Other issues, including the discrepancy in couch vertical at different couch longitudinal positions, and the inaccuracy in couch table longitudinal movement, were manually corrected during the reconstruction process. In conclusion, the phantom studies showed that, using DSCS, a 39.5 cm longitudinal coverage with a 45 cm FOV was accomplished. The efficiency of imaging dose usage was near 100%. This proposed method could be potentially useful for image guidance and subsequent treatment plan adaptation. [ABSTRACT FROM AUTHOR]

Modern external beam radiotherapy requires detailed verification and quality assurance so that confidence can be placed on both the delivery of a single treatment fraction and on the consistency of delivery throughout the treatment course. To verify dose distributions, a comparison between prediction and measurement must be made. Comparisons between two dose distributions are commonly performed using a Gamma evaluation which is a calculation of two quantities on a pixel by pixel basis; the dose difference, and the distance to agreement. By providing acceptance criteria (e.g. 3%, 3 mm), the function will find the most appropriate match within its two degrees of freedom. For complex dynamic treatments such as IMRT or VMAT it is important to verify the dose delivery in a time dependent manner and so a gamma evaluation that includes a degree of freedom in the time domain via a third parameter, time to agreement, is presented here.A C++ (mex) based gamma function was created that could be run on either CPU and GPU computing platforms that would allow a degree of freedom in the time domain. Simple test cases were created in both 2D and 3D comprising of simple geometrical shapes with well-defined boundaries varying over time. Changes of varying magnitude in either space or time were introduced and repeated gamma analyses were performed varying the criteria. A clinical VMAT case was also included, artificial air bubbles of varying size were introduced to a patient geometry, along with shifts of varying magnitude in treatment time.For all test cases where errors in distance, dose or time were introduced, the time dependent gamma evaluation could accurately highlight the errors.The time dependent gamma function presented here allows time to be included as a degree of freedom in gamma evaluations. The function allows for 2D and 3D data sets which are varying over time to be compared using appropriate criteria without penalising minor offsets of subsequent radiation fields in time. [ABSTRACT FROM AUTHOR]

With the aim of maximally reducing imaging dose while meeting requirements for adaptive radiation therapy (ART), we propose in this paper a new cone beam CT (CBCT) acquisition and reconstruction method that delivers images with a low noise level inside a region of interest (ROI) and a relatively high noise level outside the ROI. The acquired projection images include two groups: densely sampled projections at a low exposure with a large field of view (FOV) and sparsely sampled projections at a high exposure with a small FOV corresponding to the ROI. A new algorithm combining the conventional filtered back-projection algorithm and the tight-frame iterative reconstruction algorithm is also designed to reconstruct the CBCT based on these projection data. We have validated our method on a simulated head-and-neck (HN) patient case, a semi-real experiment conducted on a HN cancer patient under a full-fan scan mode, as well as a Catphan phantom under a half-fan scan mode. Relative root-mean-square errors (RRMSEs) of less than 3% for the entire image and ~1% within the ROI compared to the ground truth have been observed. These numbers demonstrate the ability of our proposed method to reconstruct high-quality images inside the ROI. As for the part outside ROI, although the images are relatively noisy, it can still provide sufficient information for radiation dose calculations in ART. Dose distributions calculated on our CBCT image and on a standard CBCT image are in agreement, with a mean relative difference of 0.082% inside the ROI and 0.038% outside the ROI. Compared with the standard clinical CBCT scheme, an imaging dose reduction of approximately 3–6 times inside the ROI was achieved, as well as an 8 times outside the ROI. Regarding computational efficiency, it takes 1–3 min to reconstruct a CBCT image depending on the number of projections used. These results indicate that the proposed method has the potential for application in ART. [ABSTRACT FROM AUTHOR]