PubMed Health⌕ Search

Biomedical subjects

Seung Wook Lee

Publications and source records attributed to Seung Wook Lee.

3 recordsLinked to original sources

Grangeat-type helical half-scan computerized tomography algorithm for reconstruction of a short object.

Currently, cone-beam computerized tomography (CT) and micro-CT scanners are under rapid development for major biomedical applications. Half-scan cone-beam image reconstruction algorithms assume only part of a scanning turn, and are advantageous in terms of temporal resolution and image artifacts. While the existing half-scan cone-beam algorithms are in the Feldkamp framework, we have published a half-scan algorithm in the Grangeat framework for a circular trajectory [Med. Phys. 30, 689-700 (2003)]. In this paper, we extend our previous work to a helical case without data truncation. We modify the Grangeat's formula for utilization and estimation of Radon data. Specifically, we categorize each characteristic point in the Radon space into singly, doubly, triply sampled, and shadow regions, respectively. A smooth weighting strategy is designed to compensate for data redundancy and inconsistency. In the helical half-scan case, the concepts of projected trajectories and transition points on meridian planes are introduced to guide the design of weighting functions. Then, the shadow region is recovered via linear interpolation after smooth weighting. The Shepp-Logan phantom is used to verify the correctness of the formulation, and demonstrate the merits of the Grangeat-type half-scan algorithm. Our Grangeat-type helical half-scan algorithm is not only valuable for quantitative and/or dynamic biomedical applications of CT and micro-CT, but also serves as an intermediate step towards solving the long object problem.

Algorithms↗

A Grangeat-type half-scan algorithm for cone-beam CT.

Modern CT and micro-CT scanners are rapidly moving from fan-beam toward cone-beam geometry. Half-scan CT algorithms are advantageous in terms of temporal resolution, and widely used in fan-beam and cone-beam geometry. While existing half-scan algorithms for cone-beam CT are in the Feldkamp framework, in this paper we compensate missing data explicitly in the Grangeat framework, and formulate a half-scan algorithm in the circular scanning case. The half-scan spans 180 degrees plus two cone angles that guarantee sufficient data for reconstruction of the midplane defined by the source trajectory. The smooth half-scan weighting functions are designed for the suppression of data inconsistency. Numerical simulation results are reported for verification of our formulas and programs. This Grangeat-type half-scan algorithm produces excellent image quality, without off-mid-plane artifacts associated with Feldkamp-type half-scan algorithms. The Grangeat-type half-scan algorithm seems promising for quantitative and dynamic biomedical applications of CT and micro-CT.

Algorithms↗

Artifacts associated with implementation of the Grangeat formula.

To compensate for image artifacts introduced in approximate cone-beam reconstruction, exact cone-beam reconstruction algorithms are being developed for medical x-ray CT. Although the exact cone-beam approach is theoretically error-free, it is subject to image artifacts due to the discrete nature of numerical implementation. We report a study on image artifacts associated with the Grangeat algorithm as applied to a circular scanning locus. Three types of artifacts are found, which are thorn, wrinkle, and V-shaped artifacts. The thorn pattern is created by inappropriate extrapolation into the shadow zone in the radon domain. If the shadow zone is filled in with continuous data, the thorn artifacts along the boundary of the shadow zone can be removed. The wrinkle appearance arises if interpolated first derivatives of the radon data are not smooth between adjacent detector planes. In particular, the nearest-neighbor interpolation method should not be used. If the number of projections is not small, the bilinear interpolation method is effective to suppress the wrinkle artifacts. The V-shaped artifacts on the meridian plane come from the line integrations through the transition zones where derivative data change abruptly. Two remedies are to increase the sampling rate and suppress data noise.

Algorithms↗