The feasibility of an inverse geometry CT system with stationary source arrays

Generalized-equiangular geometry CT: Concept and shift-invariant FBP algorithms

BACKGROUND

With advanced x-ray source and detector technologies being continuously developed, non-traditional CT geometries have been widely explored. Generalized-Equiangular Geometry CT (GEGCT) architecture, in which an x-ray source might be positioned radially far away from the focus of arced detector array that is equiangularly spaced, is of importance in many novel CT systems and designs.

PURPOSE

GEGCT, unfortunately, has no theoretically exact and shift-invariant analytical image reconstruction algorithm in general. In this study, to obtain fast and accurate reconstruction from GEGCT and to promote its system design and optimization, an in-depth investigation on a group of approximate Filtered Back-Projection (FBP) algorithms with a variety of weighting strategies has been conducted.

METHODS

The architecture of GEGCT is first presented and characterized by using a normalized-radial-offset distance (NROD). Next, shift-invariant weighted FBP-type algorithms are derived in a unified framework, with pre-filtering, filtering, and post-filtering weights, for both fixed and dynamic NROD configurations. Three viable weighting strategies are then presented including a classic one developed by Besson in the literature and two new ones generated from a curvature fitting and from an empirical formula, where all of the three weights can be expressed as certain functions of NROD. After that, an analysis of reconstruction accuracy is conducted with a wide range of NROD. Finally, the weighted FBP algorithm for GEGCT is extended to a three-dimensional form in the case of cone-beam scan with a cylindrical detector array.

RESULTS

Theoretical analysis and numerical study show that weights in the shift-invariant FBP algorithms can guarantee highly accurate reconstruction for GEGCT. A simulation of Shepp-Logan phantom and a GEGCT scan of lung mimicked by using a clinical lung CT dataset both demonstrate that FBP reconstructions with Besson and polynomial weights can achieve excellent image quality, with Peak Signal to Noise Ratio and Structural Similarity being at the same level as that from the standard equiangular fan-beam CT scan. Reconstruction of a cylinder object with multiple contrasts from simulated GEGCT scan with dynamic NROD is also highly consistent with fixed ones when using the Besson and polynomial weights, with root mean square error less than 7 hounsfield units, demonstrating the robustness and flexibility of the presented FBP algorithms. In terms of resolution, the direct FBP methods for GEGCT could achieve 1.35 lp/mm of spatial resolution at 10% modulation transfer functions point, higher than that of the rebinning method which can only reach 1.14 lp/mm. Moreover, 3D reconstructions of a disc phantom reveal that a greater value of NROD for GEGCT will bring less cone beam artifacts as expected.

CONCLUSIONS

We propose the concept of GEGCT and investigate the feasibility of using shift-invariant weighted FBP-type algorithms for reconstruction from GEGCT data without rebinning. A comprehensive analysis and phantom studies have been conducted to validate the effectiveness of proposed weighting strategies in a wide range of NROD for GEGCT with fixed and dynamic NROD.

Four-dimensional inverse-geometry computed tomography: a preliminary study

This study introduces and evaluates respiratory-correlated four-dimensional (4D) inverse geometry computed tomography (IGCT). The projection data of the IGCT were acquired in a single gantry rotation over 120 s. Three virtual phantoms-static Defrise, 4D Shepp-Logan, and 4D extended cardiac-torso (XCAT)-were used to obtain projection data for the IGCT and cone-beam computed tomography (CBCT). The projection acquisition parameters were determined to eliminate vacancies in the Radon space for an accurate rebinning process. Phase-based sorting was conducted within 10 phase bins, and the sorted projection data were binned into a cone beam geometry. Finally, Feldkamp-Davis-Kress reconstruction was conducted independently at each phase. The reconstructed images were compared using the structural similarity index measure (SSIM) and root mean square error (RMSE). The vertical profile of the Defrise reconstruction image was uniform, and the cone beam artefact was reduced in the IGCT image. Under an ideal projection acquisition condition, the mean coronal plane SSIMs of the Shepp-Logan and 4D XCAT phantoms were 0.899 and 0.706, respectively, which were higher than those of the CBCT (0.784 and 0.623, respectively). Similarly, the mean RMSEs of the coronal plane IGCT (0.036 and 0.158) exhibited an improvement over those of the CBCT (0.165 and 0.261, respectively). The mean standard deviations of the SSIM and RMSE were lower for IGCT than for CBCT. In particular, the SSIM and RMSE of the sagittal and coronal planes of the Shepp-Logan IGCT images were stable in all phase bins; however, those of the CBCT changed depending on the phase bins. Poor image quality was observed for IGCT under inappropriate conditions. This was caused by a vacancy in the Radon space, owing to an inappropriate scan setting. Overall, the proposed 4D IGCT exhibited better image quality than conventional CBCT.

Stationary computed tomography with source and detector in linear symmetric geometry: Direct filtered backprojection reconstruction

PURPOSE

Inverse-geometry computed tomography (IGCT) could have great potential in medical applications and security inspections, and has been actively investigated in recent years. In this work, we explore a special architecture of IGCT in a stationary configuration: symmetric-geometry computed tomography (SGCT), where the x-ray source and detector are linearly distributed in a symmetric design. A direct filtered backprojection (FBP)-type algorithm is developed to analytically reconstruct images from the SGCT projections.

METHODS

In our proposed SGCT system, a big number of x-ray source points equally distributed along a straight-line trajectory will sequentially fire in an ultra-fast manner in one side, and an equispaced detector whose total length is comparable to that of the source will continuously collect data in the opposite side, as the object to be scanned moves into the imaging plane. We firstly present the overall design of SGCT. An FBP-type reconstruction algorithm is then derived for this unique imaging configuration. With finite length of x-ray source and detector arrays, projection data from one segment of SGCT scan are insufficient for an exact reconstruction. As a result, in practical applications, dual-SGCT scan whose detector segments are placed perpendicular to each other, is of particular interest and is proposed. Two segments of SGCT together can make sure that the passing rays cover at least 180 degrees for each and every point if carefully designed. In general, however, there exists a data redundancy problem for a dual-SGCT. So a weighting strategy is developed to maximize the use of projection data collected while avoid image artifacts. In addition, we further extend the fan-beam SGCT to cone beam and obtain a Feldkamp-Davis-Kress (FDK)-type reconstruction algorithm. Finally, we conduct a set of experimental studies both in simulation and on a prototype SGCT system and validate our proposed methods.

RESULTS

A simulation study using the Shepp-Logan head phantom confirms that CT images can be exactly reconstructed from dual-SGCT scan and that our proposed weighting strategy is able to handle the data redundancy properly. Compared with the rebinning-to-parallel-beam method using the forward projection of an abdominal CT dataset, our proposed method is seen to be less sensitive to data truncation. Our algorithm can achieve 10.64 lp/cm of spatial resolution at 50% modulation transfer functions point, higher than that of the rebinning method which can only reach at 9.42 lp/cm even with extremely fine interpolation. Real experiments of a cylindrical object on a prototype SGCT further prove the effectiveness and practicability of the direct FBP method proposed, with similar level of noise performance to rebinning algorithm.

CONCLUSIONS

A new concept of SGCT with linearly distributed source and detector is investigated in this work, in which spinning of sources and detectors is no longer needed during data acquisition, simplifying its system design, development, and manufacturing. A direct FBP-type algorithm is developed for analytical reconstruction from SGCT projection data. Numerical and real experiments validate our method and show that exact CT image can be reconstructed from dual-SGCT scan, where data redundancy problem can be solved by our proposed weighting function.

Reduced-Dose Imageless Needle and Patient Tracking in Interventional CT Procedures

This paper describes a new method for imageless needle and patient tracking in interventional CT procedures based on fractional CT scanning. Our method accurately locates a needle with a spherical marker attached to it at a known distance from the tip with respect to the patient in the CT scanner coordinate frame with online sparse scan sampling and without reconstructing the CT image. The key principle of our method is to detect the needle and attached spherical marker in projection (sinogram) space based on the strongly attenuated X-ray signal due to the metallic composition of the needle and the needle's thin cylindrical geometry, and based on the marker's spherical geometry. A transformation from projection space to physical space uniquely determines the location and orientation of the needle and the needle tip position. Our method works directly in projection space and simultaneously performs patient registration and needle localization for every fractional CT scanning acquisition using the same sparse set of views. We performed registration and needle tip localization in five abdomen phantom scans using a rigid needle, and obtained a voxel-size tip localization error. Our experimental results indicate a voxel-sized deviation of the localization from a comparable method in 3-D image space, with the benefit of allowing X-ray dose reduction via fractional scanning at each localization. This benefit enables more frequent tip localizations during needle insertion for a similar total dose, or a reduced total dose for the same frequency of tip localization.

Multisource inverse-geometry CT. Part I. System concept and development

PURPOSE

This paper presents an overview of multisource inverse-geometry computed tomography (IGCT) as well as the development of a gantry-based research prototype system. The development of the distributed x-ray source is covered in a companion paper [V. B. Neculaes et al., "Multisource inverse-geometry CT. Part II. X-ray source design and prototype," Med. Phys. 43, 4617-4627 (2016)]. While progress updates of this development have been presented at conferences and in journal papers, this paper is the first comprehensive overview of the multisource inverse-geometry CT concept and prototype. The authors also provide a review of all previous IGCT related publications.

METHODS

The authors designed and implemented a gantry-based 32-source IGCT scanner with 22 cm field-of-view, 16 cm z-coverage, 1 s rotation time, 1.09 × 1.024 mm detector cell size, as low as 0.4 × 0.8 mm focal spot size and 80-140 kVp x-ray source voltage. The system is built using commercially available CT components and a custom made distributed x-ray source. The authors developed dedicated controls, calibrations, and reconstruction algorithms and evaluated the system performance using phantoms and small animals.

RESULTS

The authors performed IGCT system experiments and demonstrated tube current up to 125 mA with up to 32 focal spots. The authors measured a spatial resolution of 13 lp/cm at 5% cutoff. The scatter-to-primary ratio is estimated 62% for a 32 cm water phantom at 140 kVp. The authors scanned several phantoms and small animals. The initial images have relatively high noise due to the low x-ray flux levels but minimal artifacts.

CONCLUSIONS

IGCT has unique benefits in terms of dose-efficiency and cone-beam artifacts, but comes with challenges in terms of scattered radiation and x-ray flux limits. To the authors' knowledge, their prototype is the first gantry-based IGCT scanner. The authors summarized the design and implementation of the scanner and the authors presented results with phantoms and small animals.

Multi-Mounted X-Ray Computed Tomography

Most existing X-ray computed tomography (CT) techniques work in single-mounted mode and need to scan the inspected objects one by one. It is time-consuming and not acceptable for the inspection in a large scale. In this paper, we report a multi-mounted CT method and its first engineering implementation. It consists of a multi-mounted scanning geometry and the corresponding algebraic iterative reconstruction algorithm. This approach permits the CT rotation scanning of multiple objects simultaneously without the increase of penetration thickness and the signal crosstalk. Compared with the conventional single-mounted methods, it has the potential to improve the imaging efficiency and suppress the artifacts from the beam hardening and the scatter. This work comprises a numerical study of the method and its experimental verification using a dataset measured with a developed multi-mounted X-ray CT prototype system. We believe that this technique is of particular interest for pushing the engineering applications of X-ray CT.

A novel concept for CT with fixed anodes (FACT): Medical imaging based on the feasibility of thermal load capacity

Focussing primarily on thermal load capacity, we describe the performance of a novel fixed anode CT (FACT) compared with a 100 kW reference CT. Being a fixed system, FACT has no focal spot blurring of the X-ray source during projection. Monte Carlo and finite element methods were used to determine the fluence proportional to thermal capacity. Studies of repeated short-time exposures showed that FACT could operate in pulsed mode for an unlimited period. A virtual model for FACT was constructed to analyse various temporal sequences for the X-ray source ring, representing a circular array of 1160 fixed anodes in the gantry. Assuming similar detector properties at a very small integration time, image quality was investigated using an image reconstruction library. Our model showed that approximately 60 gantry rounds per second, i.e. 60 sequential targetings of the 1160 anodes per second, were required to achieve a performance level equivalent to that of the reference CT (relative performance, RP = 1) at equivalent image quality. The optimal projection duration in each direction was about 10 μs. With a beam pause of 1 μs between projections, 78.4 gantry rounds per second with consecutive source activity were thermally possible at a given thermal focal spot. The settings allowed for a 1.3-fold (RP = 1.3) shorter scan time than conventional CT while maintaining radiation exposure and image quality. Based on the high number of rounds, FACT supports a high image frame rate at low doses, which would be beneficial in a wide range of diagnostic and technical applications.

A dynamic attenuator improves spectral imaging with energy-discriminating, photon counting detectors

Energy-discriminating, photon counting (EDPC) detectors have high potential in spectral imaging applications but exhibit degraded performance when the incident count rate approaches or exceeds the characteristic count rate of the detector. In order to reduce the requirements on the detector, we explore the strategy of modulating the X-ray flux field using a recently proposed dynamic, piecewise-linear attenuator. A previous paper studied this modulation for photon counting detectors but did not explore the impact on spectral applications. In this work, we modeled detection with a bipolar triangular pulse shape (Taguchi et al., 2011) and estimated the Cramer-Rao lower bound (CRLB) of the variance of material selective and equivalent monoenergetic images, assuming deterministic errors at high flux could be corrected. We compared different materials for the dynamic attenuator and found that rare earth elements, such as erbium, outperformed previously proposed materials such as iron in spectral imaging. The redistribution of flux reduces the variance or dose, consistent with previous studies on benefits with conventional detectors. Numerical simulations based on DICOM datasets were used to assess the impact of the dynamic attenuator for detectors with several different characteristic count rates. The dynamic attenuator reduced the peak incident count rate by a factor of 4 in the thorax and 44 in the pelvis, and a 10 Mcps/mm (2) EDPC detector with dynamic attenuator provided generally superior image quality to a 100 Mcps/mm (2) detector with reference bowtie filter for the same dose. The improvement is more pronounced in the material images.

Studies of a prototype linear stationary x-ray source for tomosynthesis imaging

A prototype linear x-ray source to implement stationary source-stationary detector tomosynthesis (TS) imaging has been studied. Potential applications include human breast and small animal imaging. The source is comprised of ten x-ray source elements each consisting of a field emission cathode, electrostatic lens, and target. The electrostatic lens and target are common to all elements. The source elements form x-ray focal spots with minimum diameters of 0.3-0.4 mm at electron beam currents of up to 40 mA with a beam voltage of 40 kV. The x-ray flux versus time was quantified from each source. X-ray bremsstrahlung spectra from tungsten targets were produced using electron beam energies from 35 to 50 keV. The half-value layer was measured to be 0.8, 0.9, and 1.0 mm, respectively, for the 35, 40, and 45 kV tube potentials using the tungsten target. The suppression of voltage breakdown events, particularly during source operation, and the use of a modified form of the standard cold-cathode geometry, enhanced source reliability. The prototype linear source was used to collect tomographic data sets of a mouse phantom using digital TS reconstruction methods and demonstrated a slice-sensitivity profile with a full-width-half-maximum of 1.3 mm. Lastly, preliminary studies of tomographic imaging of flow through the mouse phantom were performed.

The piecewise-linear dynamic attenuator reduces the impact of count rate loss with photon-counting detectors

Photon counting x-ray detectors (PCXDs) offer several advantages compared to standard energy-integrating x-ray detectors, but also face significant challenges. One key challenge is the high count rates required in CT. At high count rates, PCXDs exhibit count rate loss and show reduced detective quantum efficiency in signal-rich (or high flux) measurements. In order to reduce count rate requirements, a dynamic beam-shaping filter can be used to redistribute flux incident on the patient. We study the piecewise-linear attenuator in conjunction with PCXDs without energy discrimination capabilities. We examined three detector models: the classic nonparalyzable and paralyzable detector models, and a 'hybrid' detector model which is a weighted average of the two which approximates an existing, real detector (Taguchi et al 2011 Med. Phys. 38 1089-102). We derive analytic expressions for the variance of the CT measurements for these detectors. These expressions are used with raw data estimated from DICOM image files of an abdomen and a thorax to estimate variance in reconstructed images for both the dynamic attenuator and a static beam-shaping ('bowtie') filter. By redistributing flux, the dynamic attenuator reduces dose by 40% without increasing peak variance for the ideal detector. For non-ideal PCXDs, the impact of count rate loss is also reduced. The nonparalyzable detector shows little impact from count rate loss, but with the paralyzable model, count rate loss leads to noise streaks that can be controlled with the dynamic attenuator. With the hybrid model, the characteristic count rates required before noise streaks dominate the reconstruction are reduced by a factor of 2 to 3. We conclude that the piecewise-linear attenuator can reduce the count rate requirements of the PCXD in addition to improving dose efficiency. The magnitude of this reduction depends on the detector, with paralyzable detectors showing much greater benefit than nonparalyzable detectors.