Reconstruction of coronary arteries from X-ray angiography: A review.

Despite continuous progress in X-ray angiography systems, X-ray coronary angiography is fundamentally limited by its 2D representation of moving coronary arterial trees, which can negatively impact assessment of coronary artery disease and guidance of percutaneous coronary intervention. To provide clinicians with 3D/3D+time information of coronary arteries, methods computing reconstructions of coronary arteries from X-ray angiography are required. Because of several aspects (e.g. cardiac and respiratory motion, type of X-ray system), reconstruction from X-ray coronary angiography has led to vast amount of research and it still remains as a challenging and dynamic research area. In this paper, we review the state-of-the-art approaches on reconstruction of high-contrast coronary arteries from X-ray angiography. We mainly focus on the theoretical features in model-based (modelling) and tomographic reconstruction of coronary arteries, and discuss the evaluation strategies. We also discuss the potential role of reconstructions in clinical decision making and interventional guidance, and highlight areas for future research

3D reconstruction of coronary artery using Feldkamp-Davis-Kress algorithm

An important cause of death in industrialized countries is coronary heart diseases. To treat those pathologies, a percutaneous intervention that consists in inserting a catheter in the femoral artery is performed. The instrument is directed to the affected arteries, and coronary angiography is used to lead the surgeon in an interventional context. However, 2D angiography which is frequently used during an intervention, does not consider depth, resulting in high doses of contrast agent and an extended exposure to X-ray. To mitigate the impact of these problems, medical imaging techniques such as 3D coronary artery imaging are used to assist surgeons during the intervention. Many imaging modalities are used to acquire the sequences, but the rotational angiography is favored due to its lower contrast agent use and its ease of use in an interventional context. This imaging technique allows the surgeon to guide the catheter in 3D in a clear manner, and limit the use of X-rays and contrast agent by reducing the duration of the intervention. In this thesis, we present a flexible algorithm, Feldkamp-Davis-Kress (FDK), to reconstruct 3D model of coronary artery in multiple angle views. The dual-axis rotational coronary artery angiography is proposed to use along with this algorithm. The cameras parameters are first calibrated by a nonlinear optimization where the reprojection error is minimized. Then the optimal working view is calculated to avoiding the vessel overlap and foreshortening effects. To reduce the cardiac motion effect, ECG-gated is applied into the reconstruction algorithm. The proposed method can be used in the framework to improve 3D navigation guidance in surgery. It could be a good tool for clinicians in coronary artery disease

Motion compensated iterative reconstruction for cardiac X-ray tomography

Within this Ph.D. project, three-dimensional reconstruction methods for moving objects (with a focus on the human heart) from cone-beam X-ray projections using iterative reconstruction algorithms were developed and evaluated. This project was carried in collaboration with the Digital Imaging Group of Philips Research Europe – Hamburg. In cardiac cone-beam computed tomography (CT) a large effort is continuously dedicated to increase scanning speed in order to minimize patient or organ motion during acquisition. In particular, motion causes severe artifacts such as blurring and streaks in tomographic images. While for a large class of applications the current scanning speed is sufficient, in cardiac CT image reconstruction improvements are still required. Whereas it is currently feasible to achieve stable image quality in the resting phases of the cardiac cycle, in the phase of fast motion data acquisition is too slow. A variety of algorithms to reduce or compensate for motion artifacts have been proposed in literature. Most of the correction methods address the calculation of consistent projection data belonging to the same motion state (gated CT reconstruction). Even if gated CT leads to better results, not only with respect to the processing time but also regarding the image quality, it is also limited in its temporal and spatial resolution due to the mechanical movement of the gantry. This can lead to motion blurring, especially in the phases of fast cardiac motion during the RR interval. A motion-compensated reconstruction method for CT can be used to improve the resolution of the reconstructed image and to suppress motion blurring. Iterative techniques are a promising approach to solve this problem, since no direct inversion methods are known for arbitrarily moving objects. In this work, we therefore introduced motion compensation into image reconstruction. In order to determine the unknown cardiac motion, 3 different cardiac-motion estimation methodologies were implemented. Visual and quantitative assessment of the method in a number of applications, including: phantoms; cardiac CT reconstructions; Region of Interest (ROI) CT reconstructions of left and right coronaries of several clinical patients, confirmed its potential

The Estimation and Correction of Rigid Motion in Helical Computed Tomography

X-ray CT is a tomographic imaging tool used in medicine and industry. Although technological developments have significantly improved the performance of CT systems, the accuracy of images produced by state-of-the-art scanners is still often limited by artefacts due to object motion. To tackle this problem, a number of motion estimation and compensation methods have been proposed. However, no methods with the demonstrated ability to correct for rigid motion in helical CT scans appear to exist. The primary aims of this thesis were to develop and evaluate effective methods for the estimation and correction of arbitrary six degree-of-freedom rigid motion in helical CT. As a first step, a method was developed to accurately estimate object motion during CT scanning with an optical tracking system, which provided sub-millimetre positional accuracy. Subsequently a motion correction method, which is analogous to a method previously developed for SPECT, was adapted to CT. The principle is to restore projection consistency by modifying the source-detector orbit in response to the measured object motion and reconstruct from the modified orbit with an iterative reconstruction algorithm. The feasibility of this method was demonstrated with a rapidly moving brain phantom, and the efficacy of correcting for a range of human head motions acquired from healthy volunteers was evaluated in simulations. The methods developed were found to provide accurate and artefact-free motion corrected images with most types of head motion likely to be encountered in clinical CT imaging, provided that the motion was accurately known. The method was also applied to CT data acquired on a hybrid PET/CT scanner demonstrating its versatility. Its clinical value may be significant by reducing the need for repeat scans (and repeat radiation doses), anesthesia and sedation in patient groups prone to motion, including young children

Improving Statistical Image Reconstruction for Cardiac X-ray Computed Tomography.

Technological advances in CT imaging pose new challenges such as increased X-ray radiation dose and complexity of image reconstruction. Statistical image reconstruction methods use realistic models that incorporate the physics of the measurements and the statistical properties of the measurement noise, and they have potential to provide better image quality and dose reduction compared to the conventional filtered back-projection (FBP) method. However, statistical methods face several challenges that should be addressed before they can replace the FBP method universally. In this thesis, we develop various methods to overcome these challenges of statistical image reconstruction methods. Rigorous regularization design methods in Fourier domain were proposed to achieve more isotropic and uniform spatial resolution or noise properties. The design framework is general so that users can control the spatial resolution and the noise characteristics of the estimator. In addition, a regularization design method based on the hypothetical geometry concept was introduced to improve resolution or noise uniformity. Proposed designs using the new concept effectively improved the spatial resolution or noise uniformity in the reconstructed image. The hypothetical geometry idea is general enough to be applied to other scan geometries. Statistical weighting modification, based on how much each detector element affects insufficiently sampled region, was proposed to reduce the artifacts without degrading the temporal resolution within the region-of-interest (ROI). Another approach using an additional regularization term, that exploits information from the prior image, was investigated. Both methods effectively removed short-scan artifacts in the reconstructed image. We accelerated the family of ordered-subsets algorithms by introducing a double surrogate so that faster convergence speed can be achieved. Furthermore, we present a variable splitting based algorithm for motion-compensated image reconstruction (MCIR) problem that provides faster convergence compared to the conjugate gradient (CG) method. A sinogram-based motion estimation method that does not require any additional measurements other than the short-scan amount of data was introduced to provide decent initial estimates for the joint estimation. Proposed methods were evaluated using simulation and real patient data, and showed promising results for solving each challenge. Some of these methods can be combined to generate more complete solutions for CT imaging.PhDElectrical Engineering: SystemsUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/110319/1/janghcho_1.pd

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

Surrogate driven respiratory motion model derived from CBCT projection data

Cone Beam Computed Tomography (CBCT) is the most common imaging method for Image Guided Radiation Therapy (IGRT). However due to the slow rotating gantry, the image quality of CBCT can be adversely affected by respiratory motion, as it blurs the tumour and nearby organs at risk (OARs), which makes visualization of organ boundaries difficult, in particular for organs in the thoracic region. Currently one approach to tackle the problem of respiratory motion is the use of respiratory motion model to compensate for the motion during CBCT image reconstruction. The overall goal of this work is to estimate the 3D motion, including the breath-to-breath variability, on the day of treatment directly from the CBCT projection data, without requiring any external devices. The work presented here consist of two main parts: firstly, we introduce a novel data driven method based on Principal Component Analysis PCA, with the goal to extract a surrogate signal related to the internal anatomy from the CBCT projections. Secondly, using the extracted signals, we use surrogate-driven respiratory motion models to estimate the patient’s 3D respiratory motion. We utilized a recently developed generalized framework that unifies image registration and correspondence model fitting into a single optimization. This enables the model to be fitted directly to unsorted/unreconstructed data (CBCT projection data), thereby allowing an estimate of the patient’s respiratory motion on the day of treatment. To evaluate our methods, we have used an anthropomorphic software phantom combined with CBCT projection simulations. We have also tested the proposed method on clinical data with promising results obtained

Doctor of Philosophy

dissertationEach year in the United States, a quarter million cases of stroke are caused directly by atherosclerotic disease of the cervical carotid artery. This represents a significant portion of health care costs that could be avoided if high-risk carotid artery lesions could be detected early on in disease progression. There is mounting evidence that Magnetic Resonance Imaging of the carotid artery can better classify subjects who would benefit from interventions. Turbo Spin Echo sequences are a class of Magnetic Resonance Imaging sequences that provide a variety of tissue contrasts. While high resolution Turbo Spin Echo images have demonstrated important details of carotid artery morphology, it is evident that pulsatile blood and wall motion related to the cardiac cycle are still significant sources of image degradation. In addition, patient motion artifacts due to the relatively long scan times of Turbo Spin Echo sequences result in an unacceptable fraction of noninterpretable studies. This dissertation presents work done to detect and correct for types of voluntary and physiological patient motion