Region-of-interest reconstruction of motion-contaminated data using a weighted backprojection filtration algorithm

Nonrigid registration-based coronary artery motion correction for cardiac computed tomography

PURPOSE

X-ray computed tomography angiography (CTA) is the modality of choice to noninvasively monitor and diagnose heart disease with coronary artery health and stenosis detection being of particular interest. Reliable, clinically relevant coronary artery imaging mandates high spatiotemporal resolution. However, advances in intrinsic scanner spatial resolution (CT scanners are available which combine nearly 900 detector columns with focal spot oversampling) can be tempered by motion blurring, particularly in patients with unstable heartbeats. As a result, recently numerous methods have been devised to improve coronary CTA imaging. Solutions involving hardware, multisector algorithms, or β-blockers are limited by cost, oversimplifying assumptions about cardiac motion, and populations showing contraindications to drugs, respectively. This work introduces an inexpensive algorithmic solution that retrospectively improves the temporal resolution of coronary CTA without significantly affecting spatial resolution.

METHODS

Given the goal of ruling out coronary stenosis, the method focuses on "deblurring" the coronary arteries. The approach makes no assumptions about cardiac motion, can be used on exams acquired at high heart rates (even over 75 beats/min), and draws on a fast and accurate three-dimensional (3D) nonrigid bidirectional labeled point matching approach to estimate the trajectories of the coronary arteries during image acquisition. Motion compensation is achieved by employing a 3D warping of a series of partial reconstructions based on the estimated motion fields. Each of these partial reconstructions is created from data acquired over a short time interval. For brevity, the algorithm "Subphasic Warp and Add" (SWA) reconstruction.

RESULTS

The performance of the new motion estimation-compensation approach was evaluated by a systematic observer study conducted using nine human cardiac CTA exams acquired over a range of average heart rates between 68 and 86 beats/min. Algorithm performance was based-lined against exams reconstructed using standard filtered-backprojection (FBP). The study was performed by three experienced reviewers using the American Heart Association's 15-segment model. All vessel segments were evaluated to quantify their viability to allow a clinical diagnosis before and after motion estimation-compensation using SWA. To the best of the authors' knowledge this is the first such observer study to show that an image processing-based software approach can improve the clinical diagnostic value of CTA for coronary artery evaluation.

CONCLUSIONS

Results from the observer study show that the SWA method described here can dramatically reduce coronary artery motion and preserve real pathology, without affecting spatial resolution. In particular, the method successfully mitigated motion artifacts in 75% of all initially nondiagnostic coronary artery segments, and in over 45% of the cases this improvement was enough to make a previously nondiagnostic vessel segment clinically diagnostic.

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.

Diffuse lung disease: CT of the chest with adaptive statistical iterative reconstruction technique

PURPOSE

To compare visualization of subtle normal and abnormal findings at computed tomography (CT) of the chest for diffuse lung disease with images reconstructed with filtered back projection and adaptive statistical iterative reconstruction (ASIR) techniques.

MATERIALS AND METHODS

In this HIPAA-compliant, institutional review board-approved study, 24 patients underwent 64-section multi-detector row CT of the chest for evaluation of diffuse lung disease. Scanning parameters included a pitch of 0.984:1 and 120 kVp in thin-section mode, with 2496 views per rotation compared with 984 views acquired for normal mode. The 0.625-mm-thick images were reconstructed with filtered back projection, ASIR, and ASIR high-definition (ASIR-HD) kernels. Two thoracic radiologists independently assessed the filtered back projection, ASIR, and ASIR-HD images for small anatomic details (interlobular septa, centrilobular region, and small bronchi and bronchioles), abnormal findings (reticulation, tiny nodules, altered attenuation, bronchiectasis), image quality (graded by using a six-point scale, where 1 = excellent image quality, and 5 = interpretation impossible), image noise, and artifacts. Data were tabulated for statistical testing.

RESULTS

For visualization of normal and pathologic structures, CT image series reconstructed with ASIR-HD were rated substantially better than those reconstructed with filtered back projection and ASIR (P < .001). ASIR-HD images were superior to filtered back projection images in 15 of 24 (62%) patients for visualization of normal structures and in 24 of 24 (100%) patients for pathologic findings. ASIR-HD was superior to ASIR in three of 24 (12%) images for normal anatomic findings and in seven of 24 (29%) images for pathologic evaluation. None of the images in the three groups were rated as unacceptable for noise (P < .001).

CONCLUSION

ASIR-HD reconstruction results in superior visualization of subtle and tiny anatomic structures and lesions in diffuse lung disease compared with ASIR and filtered back projection reconstructions.

Region-of-interest image reconstruction in circular cone-beam microCT

Cone-beam microcomputed tomography (microCT) is one of the most popular choices for small animal imaging which is becoming an important tool for studying animal models with transplanted diseases. Region-of-interest (ROI) imaging techniques in CT, which can reconstruct an ROI image from the projection data set of the ROI, can be used not only for reducing imaging-radiation exposure to the subject and scatters to the detector but also for potentially increasing spatial resolution of the reconstructed images. Increasing spatial resolution in microCT images can facilitate improved accuracy in many assessment tasks. A method proposed previously for increasing CT image spatial resolution entails the exploitation of the geometric magnification in cone-beam CT. Due to finite detector size, however, this method can lead to data truncation for a large geometric magnification. The Feldkamp-Davis-Kress (FDK) algorithm yields images with artifacts when truncated data are used, whereas the recently developed backprojection filtration (BPF) algorithm is capable of reconstructing ROI images without truncation artifacts from truncated cone-beam data. We apply the BPF algorithm to reconstructing ROI images from truncated data of three different objects acquired by our circular cone-beam microCT system. Reconstructed images by use of the FDK and BPF algorithms from both truncated and nontruncated cone-beam data are compared. The results of the experimental studies demonstrate that, from certain truncated data, the BPF algorithm can reconstruct ROI images with quality comparable to that reconstructed from nontruncated data. In contrast, the FDK algorithm yields ROI images with truncation artifacts. Therefore, an implication of the studies is that, when truncated data are acquired with a configuration of a large geometric magnification, the BPF algorithm can be used for effective enhancement of the spatial resolution of a ROI image.

Temporally targeted imaging method applied to ECG-gated computed tomography: preliminary phantom and in vivo experience

RATIONALE AND OBJECTIVES

Existing cardiac imaging methods do not allow for improved temporal resolution when considering a targeted region of interest (ROI). The imaging method presented here enables improved temporal resolution for ROI imaging (namely, a reconstruction volume smaller than the complete field of view). Clinically, temporally targeted reconstruction would not change the primary means of reconstructing and evaluating images, but rather would enable the adjunct technique of ROI imaging, with improved temporal resolution compared with standard reconstruction ( approximately 20% smaller temporal scan window). In gated cardiac computed tomography (CT) scans improved temporal resolution directly translates into a reduction in motion artifacts for rapidly moving objects such as the coronary arteries.

MATERIALS AND METHODS

Retrospectively electrocardiogram gated coronary angiography data from a 64-slice CT system were used. A motion phantom simulating the motion profile of a coronary artery was constructed and scanned. Additionally, an in vivo study was performed using a porcine model. Comparisons between the new reconstruction technique and the standard reconstruction are given for an ROI centered on the right coronary artery, and a pulmonary ROI.

RESULTS

In both a well-controlled motion model and a porcine model results show a decrease in motion induced artifacts including motion blur and streak artifacts from contrast enhanced vessels within the targeted ROIs, as assessed through both qualitative and quantitative observations.

CONCLUSION

Temporally targeted reconstruction techniques demonstrate the potential to reduce motion artifacts in coronary CT. Further study is warranted to demonstrate the conditions under which this technique will offer direct clinical utility. Improvement in temporal resolution for gated cardiac scans has implications for improving: contrast enhanced CT angiography, calcium scoring, and assessment of the pulmonary anatomy.

Approximate and exact cone-beam reconstruction with standard and non-standard spiral scanning

The long object problem is practically important and theoretically challenging. To solve the long object problem, spiral cone-beam CT was first proposed in 1991, and has been extensively studied since then. As a main feature of the next generation medical CT, spiral cone-beam CT has been greatly improved over the past several years, especially in terms of exact image reconstruction methods. Now, it is well established that volumetric images can be exactly and efficiently reconstructed from longitudinally truncated data collected along a rather general scanning trajectory. Here we present an overview of some key results in this area.