Assessing cardiac function from total-variation-regularized 4D C-arm CT in the presence of angular undersampling.

Time-resolved tomographic cardiac imaging using an angiographic C-arm device may support clinicians during minimally invasive therapy by enabling a thorough analysis of the heart function directly in the catheter laboratory. However, clinically feasible acquisition protocols entail a highly challenging reconstruction problem which suffers from sparse angular sampling of the trajectory. Compressed sensing theory promises that useful images can be recovered despite massive undersampling by means of sparsity-based regularization. For a multitude of reasons-most notably the desired reduction of scan time, dose and contrast agent required-it is of great interest to know just how little data is actually sufficient for a certain task. In this work, we apply a convex optimization approach based on primal-dual splitting to 4D cardiac C-arm computed tomography. We examine how the quality of spatially and temporally total-variation-regularized reconstruction degrades when using as few as [Formula: see text] projection views per heart phase. First, feasible regularization weights are determined in a numerical phantom study, demonstrating the individual benefits of both regularizers. Secondly, a task-based evaluation is performed in eight clinical patients. Semi-automatic segmentation-based volume measurements of the left ventricular blood pool performed on strongly undersampled images show a correlation of close to 99% with measurements obtained from less sparsely sampled data.

Image artefact propagation in motion estimation and reconstruction in interventional cardiac C-arm CT.

The acquisition of data for cardiac imaging using a C-arm computed tomography system requires several seconds and multiple heartbeats. Hence, incorporation of motion correction in the reconstruction step may improve the resulting image quality. Cardiac motion can be estimated by deformable three-dimensional (3D)/3D registration performed on initial 3D images of different heart phases. This motion information can be used for a motion-compensated reconstruction allowing the use of all acquired data for image reconstruction. However, the result of the registration procedure and hence the estimated deformations are influenced by the quality of the initial 3D images. In this paper, the sensitivity of the 3D/3D registration step to the image quality of the initial images is studied. Different reconstruction algorithms are evaluated for a recently proposed cardiac C-arm CT acquisition protocol. The initial 3D images are all based on retrospective electrocardiogram (ECG)-gated data. ECG-gating of data from a single C-arm rotation provides only a few projections per heart phase for image reconstruction. This view sparsity leads to prominent streak artefacts and a poor signal to noise ratio. Five different initial image reconstructions are evaluated: (1) cone beam filtered-backprojection (FDK), (2) cone beam filtered-backprojection and an additional bilateral filter (FFDK), (3) removal of the shadow of dense objects (catheter, pacing electrode, etc) before reconstruction with a cone beam filtered-backprojection (cathFDK), (4) removal of the shadow of dense objects before reconstruction with a cone beam filtered-backprojection and a bilateral filter (cathFFDK). The last method (5) is an iterative few-view reconstruction (FV), the prior image constrained compressed sensing combined with the improved total variation algorithm. All reconstructions are investigated with respect to the final motion-compensated reconstruction quality. The algorithms were tested on a mathematical phantom data set with and without a catheter and on two porcine models using qualitative and quantitative measures. The quantitative results of the phantom experiments show that if no dense object is present within the scan field of view, the quality of the FDK initial images is sufficient for motion estimation via 3D/3D registration. When a catheter or pacing electrode is present, the shadow of these objects needs to be removed before the initial image reconstruction. An additional bilateral filter shows no major improvement with respect to the final motion-compensated reconstruction quality. The results with respect to image quality of the cathFDK, cathFFDK and FV images are comparable. In conclusion, in terms of computational complexity, the algorithm of choice is the cathFDK algorithm.

Residual motion compensation in ECG-gated interventional cardiac vasculature reconstruction.

Three-dimensional reconstruction of cardiac vasculature from angiographic C-arm CT (rotational angiography) data is a major challenge. Motion artefacts corrupt image quality, reducing usability for diagnosis and guidance. Many state-of-the-art approaches depend on retrospective ECG-gating of projection data for image reconstruction. A trade-off has to be made regarding the size of the ECG-gating window. A large temporal window is desirable to avoid undersampling. However, residual motion will occur in a large window, causing motion artefacts. We present an algorithm to correct for residual motion. Our approach is based on a deformable 2D-2D registration between the forward projection of an initial, ECG-gated reconstruction, and the original projection data. The approach is fully automatic and does not require any complex segmentation of vasculature, or landmarks. The estimated motion is compensated for during the backprojection step of a subsequent reconstruction. We evaluated the method using the publicly available CAVAREV platform and on six human clinical datasets. We found a better visibility of structure, reduced motion artefacts, and increased sharpness of the vessels in the compensated reconstructions compared to the initial reconstructions. At the time of writing, our algorithm outperforms the leading result of the CAVAREV ranking list. For the clinical datasets, we found an average reduction of motion artefacts by 13 ± 6%. Vessel sharpness was improved by 25 ± 12% on average.

Interventional 4D motion estimation and reconstruction of cardiac vasculature without motion periodicity assumption.

Anatomical and functional information of cardiac vasculature is a key component in the field of interventional cardiology. With the technology of C-arm CT it is possible to reconstruct static intraprocedural 3D images from angiographic projection data. Current approaches attempt to add the temporal dimension (4D). In the assumption of periodic heart motion, ECG-gating techniques can be used. However, arrhythmic heart signals and slight breathing motion are degrading image quality frequently. To overcome those problems, we present a reconstruction method based on a 4D time-continuous B-spline motion field. The temporal component of the motion field is parameterized by the acquisition time and does not assume a periodic heart motion. The analytic dynamic FDK-reconstruction formula is used directly for the motion estimation and image reconstruction. In a physical phantom experiment two vessels of size 3.1mm and 2.3mm were reconstructed using the proposed method and an algorithm with periodicity assumption. For a periodic motion both methods obtained an error of 0.1mm. For a non-periodic motion the proposed method was superior, obtaining an error of 0.3mm/0.2mm in comparison to 1.2mm/1.0mm for the algorithm with periodicity assumption. For a clinical test case of a left coronary artery it could be further shown that the method is capable to produce diameter measurements with an absolute error of 0.1mm compared to state-of-the-art measurement tools from orthogonal coronary angiography. Further, it is shown for three different clinical cases (left/right coronary artery, coronary sinus) that the proposed method is able to handle a large variability of vascular structures and motion patterns. The complete algorithm is hardware-accelerated using the GPU requiring a computation time of less than 3min for typical clinical scenarios.

New X-ray imaging modalities and their integration with intravascular imaging and interventions.

During recent years various techniques emerged providing more detailed images and insights in the cardiovascular system. C-Arm computed tomography is currently introduced in cardiac imaging offering the potential of three dimensional imaging of the coronary arteries, the cardiac chambers, venous system and a variety of anatomic anomalies inside the interventional environment. Furthermore it might enable perfusion imaging during percutaneous coronary intervention (PCI). Intravascular ultrasound (IVUS) and optical coherence tomography (OCT) are meanwhile established tools for detailed assessment of the coronary arteries. Their use might further increase with automated tissue characterization, three dimensional reconstruction, integration in angiography systems, and new emerging techniques. Parameters of fluid tissue interactions are important factors in the pathogenesis of atherosclerosis. These parameters can be calculated using computational fluid dynamics based on three dimensional models of the coronary vessels which can be derived from various sources including multislice computed tomography (MSCT), C-Arm CT or 3D reconstructed IVUS or OCT. Their use in the clinical setting has yet to be determined especially with regard to their ability in increasing treatment efficiency and clinical outcome.

[Application of C-arm computed tomography in cardiology]. / Kardiale Anwendung der C-Arm-Computertomographie.

C-arm computed tomography is currently being introduced into cardiac imaging and offers the potential for three-dimensional imaging of the cardiac anatomy within the interventional environment. This detailed view is necessary to support complex interventional strategies, such as transcutaneous valve replacement, interventional therapy of atrial fibrillation, implantation of biventricular pacemakers and assessment of myocardial perfusion. Currently, the major limitation of this technology is its insufficient temporal resolution which limits the visualization of fast moving parts of the heart.

Filtering point spread function in backprojection cone-beam CT and its applications in long object imaging.

In backprojection cone-beam CT the cone-beam projection images are first filtered, then 3D backprojected into the object space. In this paper the point spread function (PSF) for the filtering operation is studied. For the cases where the normalization matrix is a constant, i.e. all integration planes intersect the scan path the same number of times, the derivation of the PSF is extended to the general case of limited angular range for the Radon line integrals. It is found that the 2D component of the PSF can be reduced to the form of space-variant 1D Hilbert transforms. The application of the PSF to a number of aspects in long object imaging will be discussed.

Exact radon rebinning algorithm for the long object problem in helical cone-beam CT.

This paper addresses the long object problem in helical cone-beam computed tomography. We present the PHI-method, a new algorithm for the exact reconstruction of a region-of-interest (ROI) of a long object from axially truncated data extending only slightly beyond the ROI. The PHI-method is an extension of the Radon-method, published by Kudo, Noo, and Defrise in issue 43 of journal Physics in Medicine and Biology. The key novelty of the PHI-method is the introduction of a virtual object fpsi(x) for each value of the azimuthal angle psi in the image space, with each virtual object having the property of being equal to the true object f(x) in some ROI omegam. We show that, for each psi, one can calculate exact Radon data corresponding to the two-dimensional (2-D) parallel-beam projection of fpsi(x) onto the meridian plane of angle psi. Given an angular range of length pi of such parallel-beam projections, the ROI omegam can be exactly reconstructed because f(x) is identical to fpsi(x) in Omegam. Simulation results are given for both the Radon-method and the PHI-method indicating that 1) for the case of short objects, the Radon- and PHI-methods produce comparable image quality, 2) for the case of long objects, the PHI-method delivers the same image quality as in the short object case, while the Radon-method fails, and 3) the image quality produced by the PHI-method is similar for a large range of pitch values.

Exact (spiral + circles) scan region-of-interest cone beam reconstruction via backprojection.

We present a (spiral + circles) scan cone beam reconstruction algorithm in which image reconstruction proceeds via backprojection in the object space. In principle, the algorithm can reconstruct sectional region-of-interest (ROI) in a long object. The approach is a generalization of the cone beam backprojection technique developed by Kudo and Saito in two aspects: the resource-demanding normalization step in the Kudo and Saito's algorithm is eliminated through the technique of data combination that we published earlier, and the elimination of the restriction that the detector be big enough to capture the entire cone beam projection of the ROI. Restricting the projection data to the appropriate angular range required by data combination can be accomplished by a masking process. Because of the simplification resulting from the elimination of the normalization step, the most time-consuming operations of the algorithm can be approximated by the efficient step of line-by-line ramp filtering the cone beam image in the direction of the scan path, plus a correction image. The correction image, which can be computed exactly, is needed because data combination is not properly matched at the mask boundary when ramp filtering is involved. Empirical two-dimensional (2-D) point spread function (PSF) is developed to improve matching with the correction image which is computed with finite samplings. The use of transition region to further improve matching is introduced. The results of testing the algorithm on simulated phantoms are presented.

Do cyanobacteria contain "mammalian-type" cytochrome oxidase?

The cytochrome content of membranes isolated from seven species of cyanobacteria was investigated in terms of conventional difference spectra, carbon monoxide difference spectra, photoaction spectra and photodissociation spectra, and by extraction of acid-labile heme followed by spectral identification. In addition, the effect of various inhibitors and activators on the oxidation of horse heart cytochrome c by the membrane was studied. Both the spectral features and the properties of the cytochrome oxidase reaction catalysed by the membranes suggested the presence of a terminal oxidase strikingly similar to mitochondrial ferrocytochrome c:oxygen oxidoreductase (EC. 1.9.3.1).