C-arm CT imaging with the extended line-ellipse-line trajectory: first implementation on a state-of-the-art robotic angiography system.

Three-dimensional cone-beam imaging has become valuable in interventional radiology. Currently, this tool, referred to as C-arm CT, employs a circular short-scan for data acquisition, which limits the axial volume coverage and yields unavoidable cone-beam artifacts. To improve flexibility in axial coverage and image quality, there is a critical need for novel data acquisition geometries and related image reconstruction algorithms. For this purpose, we previously introduced the extended line-ellipse-line trajectory, which allows complete scanning of arbitrary volume lengths in the axial direction together with adjustable axial beam collimation, from narrow to wide depending on the targeted application. A first implementation of this trajectory on a state-of-the-art robotic angiography system is reported here. More specifically, an assessment of the quality of this first implementation is presented. The assessment is in terms of geometric fidelity and repeatability, complemented with a first visual inspection of how well the implementation enables imaging an anthropomorphic head phantom. The geometric fidelity analysis shows that the ideal trajectory is closely emulated, with only minor deviations that have no impact on data completeness and clinical practicality. Also, mean backprojection errors over short-term repetitions are shown to be below the detector pixel size at field-of-view center for most views, which indicates repeatability is satisfactory for clinical utilization. These repeatability observations are further supported by values of the Structural Similarity Index Metric above 94% for reconstructions of the FORBILD head phantom from computer-simulated data based on repeated data acquisition geometries. Last, the real data experiment with the anthropomorphic head phantom shows that the high contrast features of the phantom are well reconstructed without distortions as well as without breaks or other disturbing transition zones, which was not obvious given the complexity of the data acquisition geometry and the major variations in axial coverage that occur over the scan.

Coping with real world data: Artifact reduction and denoising for motion-compensated cardiac C-arm CT.

PURPOSE: Detailed analysis of cardiac motion would be helpful for supporting clinical workflow in the interventional suite. With an angiographic C-arm system, multiple heart phases can be reconstructed using electrocardiogram gating. However, the resulting angular undersampling is highly detrimental to the quality of the reconstructed images, especially in nonideal intraprocedural imaging conditions. Motion-compensated reconstruction has previously been shown to alleviate this problem, but it heavily relies on a preliminary reconstruction suitable for motion estimation. In this work, the authors propose a processing pipeline tailored to augment these initial images for the purpose of motion estimation and assess how it affects the final images after motion compensation. METHODS: The following combination of simple, direct methods inspired by the core ideas of existing approaches proved beneficial: (a) Streak reduction by masking high-intensity components in projection domain after filtering. (b) Streak reduction by subtraction of estimated artifact volumes in reconstruction domain. (c) Denoising in spatial domain using a joint bilateral filter guided by an uncompensated reconstruction. (d) Denoising in temporal domain using an adaptive Gaussian smoothing based on a novel motion detection scheme. RESULTS: Experiments on a numerical heart phantom yield a reduction of the relative root-mean-square error from 89.9% to 3.6% and an increase of correlation with the ground truth from 95.763% to 99.995% for the motion-compensated reconstruction when the authors' processing is applied to the initial images. In three clinical patient data sets, the signal-to-noise ratio measured in an ideally homogeneous region is increased by 37.7% on average. Overall visual appearance is improved notably and some anatomical features are more readily discernible. CONCLUSIONS: The authors' findings suggest that the proposed sequence of steps provides a clear advantage over an arbitrary sequence of individual image enhancement methods and is fit to overcome the issue of lacking image quality in motion-compensated C-arm imaging of the heart. As for future work, the obtained results pave the way for investigating how accurately cardiac functional motion parameters can be determined with this modality.

Extended ellipse-line-ellipse trajectory for long-object cone-beam imaging with a mounted C-arm system.

Recent reports show that three-dimensional cone-beam (CB) imaging with a floor-mounted (or ceiling-mounted) C-arm system has become a valuable tool in interventional radiology. Currently, a circular short scan is used for data acquisition, which inevitably yields CB artifacts and a short coverage in the direction of the patient table. To overcome these two limitations, a more sophisticated data acquisition geometry is needed. This geometry should be complete in terms of Tuy's condition and should allow continuous scanning, while being compatible with the mechanical constraints of mounted C-arm systems. Additionally, the geometry should allow accurate image reconstruction from truncated data. One way to ensure such a feature is to adopt a trajectory that provides full R-line coverage within the field-of-view (FOV). An R-line is any segment of line that connects two points on a source trajectory, and the R-line coverage is the set of points that belong to an R-line. In this work, we propose a novel geometry called the extended ellipse-line-ellipse (ELE) for long-object imaging with a mounted C-arm system. This trajectory is built from modules consisting of two elliptical arcs connected by a line. We demonstrate that the extended ELE can be configured in many ways so that full R-line coverage is guaranteed. Both tight and relaxed parametric settings are presented. All results are supported by extensive mathematical proofs provided in appendices. Our findings make the extended ELE trajectory attractive for axially-extended FOV imaging in interventional radiology.

Contrast-Enhanced C-arm Computed Tomography Imaging of Myocardial Infarction in the Interventional Suite.

OBJECTIVES: Cardiac C-arm computed tomography (CT) uses a standard C-arm fluoroscopy system rotating around the patient to provide CT-like images during interventional procedures without moving the patient to a conventional CT scanner. We hypothesized that C-arm CT can be used to visualize and quantify the size of perfusion defects and late enhancement resulting from a myocardial infarction (MI) using contrast-enhanced techniques similar to previous CT and magnetic resonance imaging studies. MATERIALS AND METHODS: A balloon occlusion followed by reperfusion in a coronary artery was used to study acute and subacute MI in 12 swine. Electrocardiographically gated C-arm CT images were acquired the day of infarct creation (n = 6) or 4 weeks after infarct creation (n = 6). The images were acquired immediately after contrast injection, then at 1 minute, and every 5 minutes up to 30 minutes with no additional contrast. The volume of the infarct as measured on C-arm CT was compared against pathology. RESULTS: The volume of acute MI, visualized as a combined region of hyperenhancement with a hypoenhanced core, correlated well with pathologic staining (concordance correlation, 0.89; P < 0.0001; mean [SD] difference, 0.67 [2.98]cm3). The volume of subacute MI, visualized as a region of hyperenhancement, correlated well with pathologic staining at imaging times 5 to 15 minutes after contrast injection (concordance correlation, 0.82; P < 0.001; mean difference, -0.64 [1.94]cm3). CONCLUSIONS: C-arm CT visualization of acute and subacute MI is possible in a porcine model, but improvement in the imaging technique is important before clinical use. Visualization of MI in the catheterization laboratory may be possible and could provide 3-dimensional images for guidance during interventional procedures.

Axially extended-volume C-arm CT using a reverse helical trajectory in the interventional room.

C-arm computed tomography (CT) is an innovative technique that enables a C-arm system to generate 3-D images from a set of 2-D X-ray projections. This technique can reduce treatment-related complications and may improve interventional efficacy and safety. However, state-of-the-art C-arm systems rely on a circular short scan for data acquisition, which limits coverage in the axial direction. This limitation was reported as a problem in hepatic vascular interventions. To solve this problem, as well as to further extend the value of C-arm CT, axially extended-volume C-arm CT is needed. For example, such an extension would enable imaging the full aorta, the peripheral arteries or the spine in the interventional room, which is currently not feasible. In this paper, we demonstrate that performing long object imaging using a reverse helix is feasible in the interventional room. This demonstration involved developing a novel calibration method, assessing geometric repeatability, implementing a reconstruction method that applies to real reverse helical data, and quantitatively evaluating image quality. Our results show that: 1) the reverse helical trajectory can be implemented and reliably repeated on a multiaxis C-arm system; and 2) a long volume can be reconstructed with satisfactory image quality using reverse helical data.

Interventional heart wall motion analysis with cardiac C-arm CT systems.

Today, quantitative analysis of three-dimensional (3D) dynamics of the left ventricle (LV) cannot be performed directly in the catheter lab using a current angiographic C-arm system, which is the workhorse imaging modality for cardiac interventions. Therefore, myocardial wall analysis is completely based on the 2D angiographic images or pre-interventional 3D/4D imaging. In this paper, we present a complete framework to study the ventricular wall motion in 4D (3D+t) directly in the catheter lab. From the acquired 2D projection images, a dynamic 3D surface model of the LV is generated, which is then used to detect ventricular dyssynchrony. Different quantitative features to evaluate LV dynamics known from other modalities (ultrasound, magnetic resonance imaging) are transferred to the C-arm CT data. We use the ejection fraction, the systolic dyssynchrony index a 3D fractional shortening and the phase to maximal contraction (ϕi, max) to determine an indicator of LV dyssynchrony and to discriminate regionally pathological from normal myocardium. The proposed analysis tool was evaluated on simulated phantom LV data with and without pathological wall dysfunctions. The LV data used is publicly available online at https://conrad.stanford.edu/data/heart. In addition, the presented framework was tested on eight clinical patient data sets. The first clinical results demonstrate promising performance of the proposed analysis tool and encourage the application of the presented framework to a larger study in clinical practice.

Evaluation of interpolation methods for surface-based motion compensated tomographic reconstruction for cardiac angiographic C-arm data.

PURPOSE: For interventional cardiac procedures, anatomical and functional information about the cardiac chambers is of major interest. With the technology of angiographic C-arm systems it is possible to reconstruct intraprocedural three-dimensional (3D) images from 2D rotational angiographic projection data (C-arm CT). However, 3D reconstruction of a dynamic object is a fundamental problem in C-arm CT reconstruction. The 2D projections are acquired over a scan time of several seconds, thus the projection data show different states of the heart. A standard FDK reconstruction algorithm would use all acquired data for a filtered backprojection and result in a motion-blurred image. In this approach, a motion compensated reconstruction algorithm requiring knowledge of the 3D heart motion is used. The motion is estimated from a previously presented 3D dynamic surface model. This dynamic surface model results in a sparse motion vector field (MVF) defined at control points. In order to perform a motion compensated reconstruction, a dense motion vector field is required. The dense MVF is generated by interpolation of the sparse MVF. Therefore, the influence of different motion interpolation methods on the reconstructed image quality is evaluated. METHODS: Four different interpolation methods, thin-plate splines (TPS), Shepard's method, a smoothed weighting function, and a simple averaging, were evaluated. The reconstruction quality was measured on phantom data, a porcine model as well as on in vivo clinical data sets. As a quality index, the 2D overlap of the forward projected motion compensated reconstructed ventricle and the segmented 2D ventricle blood pool was quantitatively measured with the Dice similarity coefficient and the mean deviation between extracted ventricle contours. For the phantom data set, the normalized root mean square error (nRMSE) and the universal quality index (UQI) were also evaluated in 3D image space. RESULTS: The quantitative evaluation of all experiments showed that TPS interpolation provided the best results. The quantitative results in the phantom experiments showed comparable nRMSE of ≈0.047 ± 0.004 for the TPS and Shepard's method. Only slightly inferior results for the smoothed weighting function and the linear approach were achieved. The UQI resulted in a value of ≈ 99% for all four interpolation methods. On clinical human data sets, the best results were clearly obtained with the TPS interpolation. The mean contour deviation between the TPS reconstruction and the standard FDK reconstruction improved in the three human cases by 1.52, 1.34, and 1.55 mm. The Dice coefficient showed less sensitivity with respect to variations in the ventricle boundary. CONCLUSIONS: In this work, the influence of different motion interpolation methods on left ventricle motion compensated tomographic reconstructions was investigated. The best quantitative reconstruction results of a phantom, a porcine, and human clinical data sets were achieved with the TPS approach. In general, the framework of motion estimation using a surface model and motion interpolation to a dense MVF provides the ability for tomographic reconstruction using a motion compensation technique.

Simulation tools for two-dimensional experiments in x-ray computed tomography using the FORBILD head phantom.

Mathematical phantoms are essential for the development and early stage evaluation of image reconstruction algorithms in x-ray computed tomography (CT). This note offers tools for computer simulations using a two-dimensional (2D) phantom that models the central axial slice through the FORBILD head phantom. Introduced in 1999, in response to a need for a more robust test, the FORBILD head phantom is now seen by many as the gold standard. However, the simple Shepp-Logan phantom is still heavily used by researchers working on 2D image reconstruction. Universal acceptance of the FORBILD head phantom may have been prevented by its significantly higher complexity: software that allows computer simulations with the Shepp-Logan phantom is not readily applicable to the FORBILD head phantom. The tools offered here address this problem. They are designed for use with Matlab®, as well as open-source variants, such as FreeMat and Octave, which are all widely used in both academia and industry. To get started, the interested user can simply copy and paste the codes from this PDF document into Matlab® M-files.

Accurate image reconstruction using real C-arm data from a Circle-plus-arc trajectory.

OBJECTIVE: Developing an efficient tool for accurate three-dimensional imaging from projections measured with C-arm systems. MATERIAL AND METHODS: A circle-plus-arc trajectory, which is complete and thus amenable to accurate reconstruction, is used. This trajectory is particularly attractive as its implementation does not require moving the patient. For reconstruction, we use the "M-line method", which allows processing the data in the efficient filtered backprojection mode. This method also offers the advantage of not requiring an ideal data acquisition geometry, i.e., the M-line algorithm can account for known deviations in the scanning geometry, which is important given that sizeable deviations are generally encountered in C-arm imaging. RESULTS: A robust implementation scheme of the "M-line method" that applies straightforwardly to real C-arm data is presented. In particular, a numerically stable technique to compute the view-dependent derivative with respect to the source trajectory parameter is applied, and an efficient way to compute the π-line backprojection intervals via a polygonal weighting mask is presented. Projection data of an anthropomorphic thorax phantom were acquired on a medical C-arm scanner and used to demonstrate the benefit of using a complete data acquisition geometry with an accurate reconstruction algorithm versus using a state-of-the-art implementation of the conventional Feldkamp algorithm with a circular short scan of cone-beam data. A significant image quality improvement based on visual assessment is shown in terms of cone-beam artifacts.

Line plus arc source trajectories and their R-line coverage for long-object cone-beam imaging with a C-arm system.

Cone-beam imaging with C-arm systems has become a valuable tool in interventional radiology. Currently, a simple circular trajectory is used, but future applications should use more sophisticated source trajectories, not only to avoid cone-beam artifacts but also to allow extended volume imaging. One attractive strategy to achieve these two goals is to use a source trajectory that consists of two parallel circular arcs connected by a line segment, possibly with repetition. In this work, we address the question of R-line coverage for such a trajectory. More specifically, we examine to what extent R-lines for such a trajectory cover a central cylindrical region of interest (ROI). An R-line is a line segment connecting any two points on the source trajectory. Knowledge of R-line coverage is crucial because a general theory for theoretically exact and stable image reconstruction from axially truncated data is only known for the points in the scanned object that lie on R-lines. Our analysis starts by examining the R-line coverage for the elemental trajectories consisting of (i) two parallel circular arcs and (ii) a circular arc connected orthogonally to a line segment. Next, we utilize our understanding of the R-lines for the aforementioned elemental trajectories to determine the R-line coverage for the trajectory consisting of two parallel circular arcs connected by a tightly fit line segment. For this trajectory, we find that the R-line coverage is insufficient to completely cover any central ROI. Because extension of the line segment beyond the circular arcs helps to increase the R-line coverage, we subsequently propose a trajectory composed of two parallel circular arcs connected by an extended line. We show that the R-lines for this trajectory can fully cover a central ROI if the line extension is long enough. Our presentation includes a formula for the minimum line extension needed to achieve full R-line coverage of an ROI with a specified size, and also includes a preliminary study on the required detector size, showing that the R-lines added by the line extension are not constraining.