Four-dimensional Cone Beam CT Reconstruction and Enhancement using a Temporal Non-Local Means Method

Four-dimensional Cone Beam Computed Tomography (4D-CBCT) has been developed to provide respiratory phase resolved volumetric imaging in image guided radiation therapy (IGRT). Inadequate number of projections in each phase bin results in low quality 4D-CBCT images with obvious streaking artifacts. In this work, we propose two novel 4D-CBCT algorithms: an iterative reconstruction algorithm and an enhancement algorithm, utilizing a temporal nonlocal means (TNLM) method. We define a TNLM energy term for a given set of 4D-CBCT images. Minimization of this term favors those 4D-CBCT images such that any anatomical features at one spatial point at one phase can be found in a nearby spatial point at neighboring phases. 4D-CBCT reconstruction is achieved by minimizing a total energy containing a data fidelity term and the TNLM energy term. As for the image enhancement, 4D-CBCT images generated by the FDK algorithm are enhanced by minimizing the TNLM function while keeping the enhanced images close to the FDK results. A forward-backward splitting algorithm and a Gauss-Jacobi iteration method are employed to solve the problems. The algorithms are implemented on GPU to achieve a high computational efficiency. The reconstruction algorithm and the enhancement algorithm generate visually similar 4D-CBCT images, both better than the FDK results. Quantitative evaluations indicate that, compared with the FDK results, our reconstruction method improves contrast-to-noise-ratio (CNR) by a factor of 2.56~3.13 and our enhancement method increases the CNR by 2.75~3.33 times. The enhancement method also removes over 80% of the streak artifacts from the FDK results. The total computation time is ~460 sec for the reconstruction algorithm and ~610 sec for the enhancement algorithm on an NVIDIA Tesla C1060 GPU card.Comment: 20 pages, 3 figures, 2 table

Improvements in Cardiac Spect/CT for the Purpose of Tracking Transplanted Cells

Regenerative therapy via stem cell transplantation has received increased attention to help treat the myocardial injury associated with heart disease. Currently, the hybridisation of SPECT with X-ray CT is expanding the utility of SPECT. This thesis compared two SPECT/CT systems for attenuation correction using slow or fast-CT attenuation maps (mu-maps). We then developed a method to localize transplanted cells in relation to compromised blood flow in the myocardium following a myocardial infarction using SPECT/CT. Finally, a method to correct for image truncation was studied for a new SPECT/CT design that incorporated small field-of-view (FOV) detectors. Computer simulations compared gated-SPECT reconstructions using slow-CT and fast-CT mu-maps with gated-CT mu-maps. Using fast-CT mu-maps improved the Root Mean Squared (RMS) error from 4.2% to 4.0%. Three canine experiments were performed comparing SPECT/CT reconstruction using the Infinia/Hawkeye-4 (slow-CT) and Symbia T6 (fast-CT). Canines were euthanized prior to imaging, and then ventilated. The results showed improvements in both RMS errors and correlation coefficients for all canines. A first-pass contrast CT imaging technique can identify regions of myocardial infarction and can be fused with SPECT. Ten canines underwent surgical ligation of the left-anterior-descending artery. Cells were labeled with 111In-tropolone and transplanted into the myocardium. SPECT/CT was performed on day of transplantation, 4, and 10 days post-transplantation. For each imaging session first-pass perfusion CT was performed and successfully delineated the infarct zone. Delayed-enhanced MRI was performed and correlated well with first-pass CT. Contrast-to-noise ratios were calculated for 111In-SPECT and suggested that cells can be followed for 11 effective half-lives. We evaluated a modified SPECT/CT acquisition and reconstruction method for truncated SPECT. Cardiac SPECT/CT scans were acquired in 14 patients. The original projections were truncated to simulate a small FOV acquisition. Data was reconstructed in three ways: non-truncated and standard reconstruction (NTOSEM), which was our gold-standard; truncated and standard reconstruction (TOSEM); and truncated and a modified reconstruction (TMOSEM). Compared with NTOSEM, small FOV imaging incurred an average cardiac count ratio error greater than 100% using TOSEM and 8.9% using TMOSEM. When we plotted NTOSEM against TOSEM and TMOSEM the correlation coefficient was 0.734 and 0.996 respectively

Improvements In computed tomography perfusion output using complex singular value decomposition and the maximum slope algorithm

OBJECTIVE: Determine if complex singular value decomposition (cSVD) used as preprocessing in the maximum slope algorithm reduces image noise of resultant physiologic parametric images. Noise will be decreased in the parametric maps of cerebral blood flow (CBF), cerebral blood volume (CBV) as compared to the same algorithm and data set with no cSVD applied. MATERIALS AND METHODS: A set of 10 patients (n=15) underwent a total combined 15 CT perfusion studies upon presenting with stroke symptoms. It was determined these patients suffered from occlusions resulting in a prolonged arrival time of blood to the brain. DICOM data files of these patients scans were selected based on this increased arrival delay. We compared the output of estimation calculations for cerebral blood flow (CBF), and cerebral blood volume (CBV), using preprocessing cSVD against the same scan data with no preprocessing cSVD. Image noise was assessed through the calculation of the standard deviation within specific regions of interest copied to specific areas of grey and white matter as well as CSF space. A decrease in the standard deviation values will indicate improvement in the noise level of the resultant images.. Results for the mean value within the regions of interest are expected to be similar between the groups calculated using cSVD and those calculated under the standard method. This will indicate the presence of minimal bias. RESULTS: Between groups of the standard processing method and the cSVD method standard deviation (SD) reductions were seen in both CBF and CBV values across all three ROIs. In grey matter measures of CBV, SD was reduced an average of 0.0034 mL/100g while measures of CBF saw SD reduced by an average of 0.073 mL/100g/min. In samples of white matter, standard deviations of CBV values were reduced on average by 0.0041mL/100g while CBF SD's were reduced by 0.073 mL/100g/min. CSF ROIs in CBV calculations saw SD reductions averaging 0.0047 mL/100g and reductions of 0.074 mL/100g/min in measures of CBF. Bias within CBV calculations was at most minimal as determined by no significant changes in mean calculated values. Calculations of CBF saw large downward bias in the mean values. CONCLUSIONS: The application of the cSVD method to preprocessing of CT perfusion imaging studies produces an effective method of noise reduction. In calculations of CBV, cSVD noise reduction results in overall improvement. In calculations of CBF, cSVD, while effective in noise reduction, caused mean values to be statistically lower than the standard method. It should be noted that there is currently no evaluation of which values can be considered more accurate physiologically. Simulations of the effect of noise on CBF showed a positive correlation suggesting that the CBF algorithm itself is sensitive to the level of noise

"Dose of the day" based on cone beam computed tomography and deformable image registration for lung cancer radiotherapy.

PURPOSE:Adaptive radiotherapy (ART) has potential to reduce toxicity and facilitate safe dose escalation. Dose calculations with the planning CT deformed to cone beam CT (CBCT) have shown promise for estimating the "dose of the day". The purpose of this study is to investigate the "dose of the day" calculation accuracy based on CBCT and deformable image registration (DIR) for lung cancer radiotherapy. METHODS:A total of 12 lung cancer patients were identified, for which daily CBCT imaging was performed for treatment positioning. A re-planning CT (rCT) was acquired after 20 Gy for all patients. A virtual CT (vCT) was created by deforming initial planning CT (pCT) to the simulated CBCT that was generated from deforming CBCT to rCT acquired on the same day. Treatment beams from the initial plan were copied to the vCT and rCT for dose calculation. Dosimetric agreement between vCT-based and rCT-based accumulated doses was evaluated using the Bland-Altman analysis. RESULTS:Mean differences in dose-volume metrics between vCT and rCT were smaller than 1.5%, and most discrepancies fell within the range of ± 5% for the target volume, lung, esophagus, and heart. For spinal cord Dmax , a large mean difference of -5.55% was observed, which was largely attributed to very limited CBCT image quality (e.g., truncation artifacts). CONCLUSION:This study demonstrated a reasonable agreement in dose-volume metrics between dose accumulation based on vCT and rCT, with the exception for cases with poor CBCT image quality. These findings suggest potential utility of vCT for providing a reasonable estimate of the "dose of the day", and thus facilitating the process of ART for lung cancer

Recent advances in x-ray cone-beam computed laminography

X-ray computed tomography is a well established volume imaging technique used routinely in medical diagnosis, industrial non-destructive testing, and a wide range of scientific fields. Traditionally, computed tomography uses scanning geometries with a single axis of rotation together with reconstruction algorithms specifically designed for this setup. Recently there has however been increasing interest in more complex scanning geometries. These include so called X-ray computed laminography systems capable of imaging specimens with large lateral dimensions, or large aspect ratios, neither of which are well suited to conventional CT scanning procedures. Developments throughout this field have thus been rapid, including the introduction of novel system trajectories, the application and refinement of various reconstruction methods, and the use of recently developed computational hardware and software techniques to accelerate reconstruction times. Here we examine the advances made in the last several years and consider their impact on the state of the art