High resolution imaging with focused kV x-rays for small animal radio-neuromodulation.

BACKGROUND: High precision radiotherapy with small irradiator size has potential in many treatment applications involving small shallow targets, with small animal radio-neuromodulation as an intriguing example. A focused kV technique based on novel usage of polycapillary x-ray lenses can focus x-ray beams to <0.2 mm in diameter, which is ideal for such uses. PURPOSE: Such an application also requires high resolution CT images for treatment planning and setup. In this work, we demonstrate the feasibility of using a virtual focal spot generated with an x-ray lens to perform high-resolution CBCT acquisition. METHOD: The experiment with x-ray lens was set up on an x-ray tabletop system to generate a virtual focal spot. The flood field images with and without the x-ray lens were first compared. A pinhole image was acquired for the virtual focal spot and compared with the one acquired with the conventional focal spot without the lens. The planar imaging resolution with and without the lens were evaluated using a line pair resolution phantom. The spatial resolution of the two settings were estimated by reconstructing a 0.15-mm wire phantom and comparing its full width half maximum (FWHM). A CBCT scan of a rodent head was also acquired to further demonstrate the improved resolution using the x-ray lens. RESULT: The proposed imaging setup with x-ray lens had a limited exposure area of 5 cm by 5 cm on the detector, which was suitable for guiding radio-neuromodulation to a small target in rodent brain. Compared to conventional imaging acquisition with a measured x-ray focal spot of 0.395 mm FWHM, the virtual focal spot size was measured at 0.175 mm. The reduction in focal spot size with lens leads to an almost doubled planar imaging resolution and a 26% enhancement in 3D spatial resolution. A realistic CBCT acquisition of a rodent head mimicked the imaging acquisition step for radio-neuromodulation and further showed the improved visualization for fine structures. CONCLUSION: This work demonstrated that the focused kV x-ray technique was capable of generating small focal spot size of <0.2 mm, which substantially improved x-ray imaging resolution for small animal imaging.

Single-Shot Quantitative X-ray Imaging Using a Primary Modulator and Dual-Layer Detector: Simulation and Phantom Studies.

Conventional x-ray imaging provides little quantitative information due to scatter, beam hardening, and overlaying tissues. A single-shot quantitative x-ray imaging (SSQI) method was previously developed to quantify material-specific densities in x-ray imaging by combining the use of a primary modulator (PM) and dual-layer (DL) detector. The feasibility of this concept was demonstrated with simulations using an iterative patch-based method. In this work, we propose a new algorithm pipeline for SSQI that enables accurate quantification and high computational efficiency. The DL images contain four measurements that are obtained behind the unattenuated and partially attenuated regions of the PM of each layer. Using the low-frequency property of scatter and a pre-calibrated material decomposition (MD), four unknowns (i.e., two scatter images and two material-specific images) are jointly recovered by directly solving four equations given by the four measurements. We tested this algorithm in simulations and further demonstrated its efficacy on chest phantom experiments. Through simulation, we show that the new method for MD is robust against scatter. Its performance improves with smaller PM pitch size and smaller focal spot blur. The RMSE in material-specific images compared to ground truth reduces by 52%-84% versus without scatter correction. For our experimental study, we successfully separated soft tissue and bone. The computational time for processing each view was ~8 s without optimization. The reported results further strengthen the potential of SSQI for widespread adoption, leading to quantitative imaging not only for x-ray imaging but also for real-time image guidance or cone-beam CT.

Characterization of x-ray focal spots using a rotating edge.

Purpose: The focal spot size and shape of an x-ray system are critical factors to the spatial resolution. Conventional approaches to characterizing the focal spot use specialized tools that usually require careful calibration. We propose an alternative to characterize the x-ray source's focal spot, simply using a rotating edge and flat-panel detector. Methods: An edge is moved to the beam axis, and an edge spread function (ESF) is obtained at a specific angle. Taking the derivative of the ESF provides the line spread function, which is the Radon transform of the focal spot in the direction parallel to the edge. By rotating the edge about the beam axis for 360 deg, we obtain a complete Radon transform, which is used for reconstructing the focal spot. We conducted a study on a clinical C-arm system with three focal spot sizes (0.3, 0.6, and 1.0 mm nominal size), then compared the focal spot imaged using the proposed method against the conventional pinhole approach. The full width at half maximum (FWHM) of the focal spots along the width and height of the focal spot were used for quantitative comparisons. Results: Using the pinhole method as ground truth, the proposed method accurately characterized the focal spot shapes and sizes. Quantitatively, the FWHM widths were 0.37, 0.65, and 1.14 mm for the pinhole method and 0.33, 0.60, and 1.15 mm for the proposed method for the 0.3, 0.6, and 1.0 mm nominal focal spots, respectively. Similar levels of agreement were found for the FWHM heights. Conclusions: The method uses a rotating edge to characterize the focal spot and could be automated in the future using a system's built-in collimator. The method could be included as part of quality assurance tests of image quality and tube health.

Single-pass metal artifact reduction using a dual-layer flat panel detector.

PURPOSE: Metal artifact remains a challenge in cone-beam CT images. Many image domain-based segmentation methods have been proposed for metal artifact reduction (MAR), which require two-pass reconstruction. Such methods first segment metal from a first-pass reconstruction and then forward-project the metal mask to identify them in projections. These methods work well in general but are limited when the metal is outside the scan field-of-view (FOV) or when the metal is moving during the scan. In the former, even reconstructing with a larger FOV does not guarantee a good estimate of metal location in the projections; and in the latter, the metal location in each projection is difficult to identify due to motion. Single-pass methods that detect metal in single-energy projections have also been developed, but often have imperfect metal detection that leads to residual artifacts. In this work, we develop a MAR method using a dual-layer (DL) flat panel detector, which improves performance for single-pass reconstruction. METHODS: In this work, we directly detect metal objects in projections using dual-energy (DE) imaging that generates material-specific images (e.g., soft tissue and bone), where the metal stands out in bone images when nonuniform soft tissue background is removed. Metal is detected via simple thresholding, and entropy filtration is further applied to remove false-positive detections. A DL detector provides DE images with superior temporal and spatial registration and was used to perform the task. Scatter correction was first performed on DE raw projections to improve the accuracy of material decomposition. One phantom mimicking a liver biopsy setup and a cadaver head were used to evaluate the metal reduction performance of the proposed method and compared with that of a standard two-pass reconstruction, a previously published sinogram-based method using a Markov random field (MRF) model, and a single-pass projection-domain method using single-energy imaging. The phantom has a liver steering setup placed in a hollow chest phantom, with embedded metal and a biopsy needle crossing the phantom boundary. The cadaver head has dental fillings and a metal tag attached to its surface. The identified metal regions in each projection were corrected by interpolation using surrounding pixels, and the images were reconstructed using filtered backprojection. RESULTS: Our current approach removes metal from the projections, which is robust to FOV truncation during imaging acquisition. In case of FOV truncation, the method outperformed the two-pass reconstruction method. The proposed method using DE renders better accuracy in metal segmentation than the MRF method and single-energy method, which were prone to false-positive errors that cause additional streaks. For the liver steering phantom, the average spatial nonuniformity was reduced from 0.127 in uncorrected images to 0.086 using a standard two-pass reconstruction and to 0.077 using the proposed method. For the cadaver head, the average standard deviation within selected soft tissue regions ( σ s ) was reduced from 209.1 HU in uncorrected images to 69.1 HU using a standard two-pass reconstruction and to 46.8 HU using our proposed method. The proposed method reduced the processing time by 31% as compared with the two-pass method. CONCLUSIONS: We proposed a MAR method that directly detects metal in the projection domain using DE imaging, which is robust to truncation and superior to that of single-energy imaging. The method requires only a single-pass reconstruction that substantially reduces processing time compared with the standard two-pass metal reduction method.

Deep learning-based reconstruction of interventional tools and devices from four X-ray projections for tomographic interventional guidance.

PURPOSE: Image guidance for minimally invasive interventions is usually performed by acquiring fluoroscopic images using a monoplanar or a biplanar C-arm system. However, the projective data provide only limited information about the spatial structure and position of interventional tools and devices such as stents, guide wires, or coils. In this work, we propose a deep learning-based pipeline for real-time tomographic (four-dimensional [4D]) interventional guidance at conventional dose levels. METHODS: Our pipeline is comprised of two steps. In the first one, interventional tools are extracted from four cone-beam CT projections using a deep convolutional neural network. These projections are then Feldkamp reconstructed and fed into a second network, which is trained to segment the interventional tools and devices in this highly undersampled reconstruction. Both networks are trained using simulated CT data and evaluated on both simulated data and C-arm cone-beam CT measurements of stents, coils, and guide wires. RESULTS: The pipeline is capable of reconstructing interventional tools from only four X-ray projections without the need for a patient prior. At an isotropic voxel size of 100 µ m , our methods achieve a precision/recall within a 100 µ m environment of the ground truth of 93%/98%, 90%/71%, and 93%/76% for guide wires, stents, and coils, respectively. CONCLUSIONS: A deep learning-based approach for 4D interventional guidance is able to overcome the drawbacks of today's interventional guidance by providing full spatiotemporal (4D) information about the interventional tools at dose levels comparable to conventional fluoroscopy.

An analysis of scatter characteristics in x-ray CT spectral correction.

X-ray scatter remains a major physics challenge in volumetric computed tomography (CT), whose physical and statistical behaviors have been commonly leveraged in order to eliminate its impact on CT image quality. In this work, we conduct an in-depth derivation of how the scatter distribution and scatter to primary ratio (SPR) will change during the spectral correction, leading to an interesting finding on the property of scatter. Such a characterization of scatter's behavior provides an analytic approach of compensating for the SPR as well as approximating the change of scatter distribution after spectral correction, even though both of them might be significantly distorted as the linearization mapping function in spectral correction could vary a lot from one detector pixel to another. We conduct an evaluation of SPR compensations (SPRCs) on a Catphan phantom and an anthropomorphic chest phantom to validate the characteristics of scatter. In addition, this scatter property is also directly adopted into CT imaging using a spectral modulator with flying focal spot technology (SMFFS) as an example to demonstrate its potential in practical applications. For cone-beam CT (CBCT) scans at both 80 and 120 kVp, CT images with accurate CT numbers can be achieved after spectral correction followed by the appropriate SPRC based on our presented scatter property. In the case of the SMFFS based CBCT scan of the Catphan phantom at 120 kVp, after a scatter correction using an analytic algorithm derived from the scatter property, CT image quality was significantly improved, with the averaged root mean square error reduced from 297.9 to 6.5 Hounsfield units.

Densely sampled spectral modulation for x-ray CT using a stationary modulator with flying focal spot: a conceptual and feasibility study of scatter and spectral correction.

PURPOSE: Modulation of the x-ray source in computed tomography (CT) by a designated filter to achieve a desired distribution of photon flux has been greatly advanced in recent years. In this work, we present a densely sampled spectral modulation (DSSM) as a promising low-cost solution to quantitative CT imaging in the presence of scatter. By leveraging a special stationary filter (namely a spectral modulator) and a flying focal spot, DSSM features a strong correlation in the scatter distributions across focal spot positions and sees no substantial projection sparsity or misalignment in data sampling, making it possible to simultaneously correct for scatter and spectral effects in a unified framework. METHODS: The concept of DSSM is first introduced, followed by an analysis of the design and benefits of using the stationary spectral modulator with a flying focal spot (SMFFS) that dramatically changes the data sampling and its associated data processing. With an assumption that the scatter distributions across focal spot positions have strong correlation, a scatter estimation and spectral correction algorithm from DSSM is then developed, where a dual-energy modulator along with two flying focal spot positions is of interest. Finally, a phantom study on a tabletop cone-beam CT system is conducted to understand the feasibility of DSSM by SMFFS, using a copper modulator and by moving the x-ray tube position in the X direction to mimic the flying focal spot. RESULTS: Based on our analytical analysis of the DSSM by SMFFS, the misalignment of low- and high-energy projection rays can be reduced by a factor of more than 10 when compared with a stationary modulator only. With respect to modulator design, metal materials such as copper, molybdenum, silver, and tin could be good candidates in terms of energy separation at a given attenuation of photon flux. Physical experiments using a Catphan phantom as well as an anthropomorphic chest phantom demonstrate the effectiveness of DSSM by SMFFS with much better CT number accuracy and less image artifacts. The root mean squared error was reduced from 297.9 to 6.5 Hounsfield units (HU) for the Catphan phantom and from 409.3 to 39.2 HU for the chest phantom. CONCLUSIONS: The concept of DSSM using a SMFFS is proposed. Phantom results on its scatter estimation and spectral correction performance validate our main ideas and key assumptions, demonstrating its potential and feasibility for quantitative CT imaging.

Characterization and potential applications of a dual-layer flat-panel detector.

PURPOSE: Dual-energy (DE) x-ray imaging has many clinical applications in radiography, fluoroscopy, and CT. This work characterizes a prototype dual-layer (DL) flat-panel detector (FPD) and investigates its DE imaging capabilities for applications in two-dimensional (2D) radiography/fluoroscopy and quantitative three-dimensional (3D) cone-beam CT. Unlike other DE methods like kV switching, a DL FPD obtains DE images from a single exposure, making it robust against patient and system motion. METHODS: The DL FPD consists of a top layer with a 200 µm-thick CsI scintillator coupled to an amorphous silicon (aSi) FPD of 150 µm pixel size and a bottom layer with a 550 µm thick CsI scintillator coupled to an identical aSi FPD. The two layers are separated by a 1-mm Cu filter to increase spectral separation. Images (43 × 43 cm2 active area) can be readout in 2 × 2 binning mode (300 µm pixels) at up to 15 frames per second. Detector performance was first characterized by measuring the MTF, NPS, and DQE for the top and bottom layers. For 2D applications, a qualitative study was conducted using an anthropomorphic thorax phantom containing a porcine heart with barium-filled coronary arteries (similar to iodine). Additionally, fluoroscopic lung tumor tracking was investigated by superimposing a moving tumor phantom on the thorax phantom. Tracking accuracies of single-energy (SE) and DE fluoroscopy were compared against the ground truth motion of the tumor. For 3D quantitative imaging, a phantom containing water, iodine, and calcium inserts was used to evaluate overall DE material decomposition capabilities. Virtual monoenergetic (VM) images ranging from 40 to 100 keV were generated, and the optimal VM image energy which achieved the highest image uniformity and maximum contrast-to-noise ratio (CNR) was determined. RESULTS: The spatial resolution of the top layer was substantially higher than that of the bottom layer (top layer 50% MTF = 2.2 mm-1 , bottom layer = 1.2 mm-1 ). A substantial increase in NNPS and reduction in DQE were observed for the bottom layer mainly due to photon loss within the top layer and Cu filter. For 2D radiographic and fluoroscopic applications, the DL FPD was capable of generating high-quality material-specific images separating soft tissue from bone and barium. For lung tumor tracking, DE fluoroscopy yielded more accurate results than SE fluoroscopy, with an average reduction in the root mean square error (RMSE) of over 10×. For the DE-CBCT studies, accurate basis material decompositions were obtained. The estimated material densities were 294.68  ±  17.41 and 92.14  ±  15.61 mg/ml for the 300 and 100 mg/ml calcium inserts, respectively, and 8.93  ±  1.45, 4.72  ±  1.44, and 2.11  ±  1.32 mg/ml for the 10, 5, and 2 mg/ml iodine inserts, respectively, with an average error of less than 5%. The optimal VM image energy was found to be 60 keV. CONCLUSIONS: We characterized a prototype DL FPD and demonstrated its ability to perform accurate single-exposure DE radiography/fluoroscopy and DE-CBCT. The merits of the DL detector approach include superior spatial and temporal registration between its constituent images, and less complicated acquisition sequences.

Comparative Study of Dual Energy Cone-Beam CT using a Dual-Layer Detector and kVp Switching for Material Decomposition.

Cone-beam CT (CBCT) is widely used in diagnostic imaging and image-guided procedures, leading to an increasing need for advanced CBCT techniques, such as dual energy (DE) imaging. Previous studies have shown that DE-CBCT can perform quantitative material decomposition, including quantification of contrast agents, electron density, and virtual monoenergetic images. Currently, most CBCT systems perform DE imaging using a kVp switching technique. However, the disadvantages of this method are spatial and temporal misregistration as well as total scan time increase, leading to errors in the material decomposition. DE-CBCT with a dual layer flat panel detector potentially overcomes these limitations by acquiring the dual energy images simultaneously. In this work, we investigate the DE imaging performance of a prototype dual layer detector by evaluating its material decomposition capability and comparing its performance to that of the kVp switching method. Two sets of x-ray spectra were used for kVp switching: 80/120 kVp and 80/120 kVp + 1 mm Cu filtration. Our results show the dual layer detector outperforms kVp switching at 80/120 kVp with matched dose. The performance of kVp switching was better by adding 1 mm copper filtration to the high energy images (80/120 kVp + 1 mm Cu), though the dual layer detector still provided comparable performance for material decomposition tasks. Overall, both the dual layer detector and kVp switching methods provided quantitative material decomposition images in DE-CBCT, with the dual layer detector having additional potential advantages.

Reconstruction of x-ray focal spot distribution using a rotating edge.

The size and shape of an x-ray source's focal spot is a critical factor in the imaging system's overall spatial resolution. The conventional approach to imaging the focal spot uses a pinhole camera, but this requires careful, manual measurements. Instead, we propose a novel alternative, simply using the collimator available on many x-ray systems. After placing the edge of a collimator blade in the center of the beam, we can obtain an image of its edge spread function (ESF). Each ESF provides information about the focal spot distribution - specifically, the parallel projection of the focal spot in the direction parallel to the edge. If the edge is then rotated about the beam axis, each image provides a different parallel projection of the focal spot until a complete Radon transform of the focal spot distribution is obtained. The focal spot can then be reconstructed by the inverse Radon transform, or parallel-beam filtered backprojection. We conducted a study on a clinical C-arm system with 3 focal spot sizes (0.3, 0.6, 1.0 mm nominal size), comparing the focal spot obtained using the rotating edge method against the conventional pinhole approach. Our results demonstrate accurate characterization of the size and shape of the focal spot.

Novel motor design for rotating anode x-ray tubes operating in the fringe field of a magnetic resonance imaging system.

PURPOSE: Using hybrid x-ray∕MR (XMR) systems for image guidance during interventional procedures could enhance the diagnosis and treatment of neurologic, oncologic, cardiovascular, and other disorders. The authors propose a close proximity hybrid system design in which a C-arm fluoroscopy unit is placed immediately adjacent to the solenoid magnet of a MR system with a minimum distance of 1.2 m between the x-ray and MR imaging fields of view. Existing rotating anode x-ray tube designs fail within MR fringe field environments because the magnetic fields alter the electron trajectories in the x-ray tube and act as a brake on the induction motor, reducing the rotation speed of the anode. In this study the authors propose a novel motor design that avoids the anode rotation speed reduction. METHODS: The proposed design replaces the permanent magnet stator found in brushed dc motors with the radial component of the MR fringe field. The x-ray tube is oriented such that the radial component of the MR fringe field is orthogonal to the cathode-anode axis. Using a feedback position sensor and the support bearings as electrical slip rings, the authors use electrical commutation to eliminate the need for mechanical brushes and commutators. A vacuum compatible prototype of the proposed motor design was assembled, and its performance was evaluated at various operating conditions. The prototype consisted of a 3.1 in. diameter anode rated at 300 kHU with a ceramic rotor that was 5.6 in. in length and had a 2.9 in. diameter. The material chosen for all ceramic components was MACOR, a machineable glass ceramic developed by Corning Inc. The approximate weight of the entire assembly was 1750 g. The maximum rotation speed, angular acceleration, and acceleration time of the motor design were investigated, as well as the dependence of these parameters on rotor angular offset, magnetic field strength, and field orientation. The resonance properties of the authors' assembly were also evaluated to determine its stability during acceleration, and a pulse width modulation algorithm was implemented to control the rotation speed of the motor. RESULTS: At a magnetic flux density of 41 mT orthogonal to the axis of rotation (on the lower end of the expected flux density in the MR suite) the maximum speed of the motor was found to be 5150 revolutions per minute (rpm). The acceleration time necessary to reach 3000 rpm was found to be approximately 10 s at 59 mT. The resonance frequency of the assembly with the anode attached was 1310 rpm (21.8 Hz) which is far below the desired operating speeds. Pulse width modulation provides an effective method to control the speed of the motor with a resolution of 100 rpm. CONCLUSIONS: The proposed design can serve as a direct replacement to the conventional induction motor used in rotating anode x-ray tubes. It does not suffer from a reduced rotation speed when operating in a MR environment. The presence of chromic steel bearings in the prototype prevented testing at the higher field strengths, and future iterations of the design could eliminate this shortcoming. The prototype assembly demonstrates proof of concept of the authors' design and overcomes one of the major obstacles for a MR compatible rotating anode x-ray tube.

Scatter correction method for x-ray CT using primary modulation: phantom studies.

PURPOSE: Scatter correction is a major challenge in x-ray imaging using large area detectors. Recently, the authors proposed a promising scatter correction method for x-ray computed tomography (CT) using primary modulation. Proof of concept was previously illustrated by Monte Carlo simulations and physical experiments on a small phantom with a simple geometry. In this work, the authors provide a quantitative evaluation of the primary modulation technique and demonstrate its performance in applications where scatter correction is more challenging. METHODS: The authors first analyze the potential errors of the estimated scatter in the primary modulation method. On two tabletop CT systems, the method is investigated using three phantoms: A Catphan 600 phantom, an anthropomorphic chest phantom, and the Catphan 600 phantom with two annuli. Two different primary modulators are also designed to show the impact of the modulator parameters on the scatter correction efficiency. The first is an aluminum modulator with a weak modulation and a low modulation frequency, and the second is a copper modulator with a strong modulation and a high modulation frequency. RESULTS: On the Catphan 600 phantom in the first study, the method reduces the error of the CT number in the selected regions of interest (ROIs) from 371.4 to 21.9 Hounsfield units (HU); the contrast to noise ratio also increases from 10.9 to 19.2. On the anthropomorphic chest phantom in the second study, which represents a more difficult case due to the high scatter signals and object heterogeneity, the method reduces the error of the CT number from 327 to 19 HU in the selected ROIs and from 31.4% to 5.7% on the overall average. The third study is to investigate the impact of object size on the efficiency of our method. The scatter-to-primary ratio estimation error on the Catphan 600 phantom without any annulus (20 cm in diameter) is at the level of 0.04, it rises to 0.07 and 0.1 on the phantom with an elliptical annulus (30 cm in the minor axis and 38 cm in the major axis) and with a circular annulus (38 cm in diameter). CONCLUSIONS: On the three phantom studies, good scatter correction performance of the proposed method has been demonstrated using both image comparisons and quantitative analysis. The theory and experiments demonstrate that a strong primary modulation that possesses a low transmission factor and a high modulation frequency is preferred for high scatter correction accuracy.

Characterization of a novel anthropomorphic plastinated lung phantom.

Phantoms are widely used during the development of new imaging systems and algorithms. For development and optimization of new imaging systems such as tomosynthesis, where conventional image quality metrics may not be applicable, a realistic phantom that can be used across imaging systems is desirable. A novel anthropomorphic lung phantom was developed by plastination of an actual pig lung. The plastinated phantom is characterized and compared with reference to in vivo images of the same tissue prior to plastination using high resolution 3D CT. The phantom is stable over time and preserves the anatomical features and relative locations of the in vivo sample. The volumes for different tissue types in the phantom are comparable to the in vivo counterparts, and CT numbers for different tissue types fall within a clinically useful range. Based on the measured CT numbers, the phantom cardiac tissue experienced a 92% decrease in bulk density and the phantom pulmonary tissue experienced a 78% decrease in bulk density compared to their in vivo counterparts. By-products in the phantom from the room temperature vulcanizing silicone and plastination process are also identified. A second generation phantom, which eliminates most of the by-products, is presented. Such anthropomorphic phantoms can be used to evaluate a wide range of novel imaging systems.

Inverse-geometry volumetric CT system with multiple detector arrays for wide field-of-view imaging.

Current volumetric computed tomography (CT) methods require seconds to acquire a thick volume (>8 cm) with high resolution. Inverse-geometry CT (IGCT) is a new system geometry under investigation that is anticipated to be able to image a thick volume in a single gantry rotation with isotropic resolution and no cone-beam artifacts. IGCT employs a large array of source spots opposite a smaller detector array. The in-plane field of view (FOV) is primarily determined by the size of the source array, in much the same way that the FOV is determined by the size of the detector array in a conventional CT system. Thus, the size of the source array can be a limitation on the achievable FOV. We propose adding additional detector arrays, spaced apart laterally, to increase the in-plane FOV while still using a modestly sized source array. We determine optimal detector placement to maximize the FOV while obtaining relatively uniform sampling. We also demonstrate low wasted radiation of the proposed system through design and simulation of a pre-patient collimator. Reconstructions from simulated projection data show no artifacts when combining the data from the detector arrays. Finally, to demonstrate feasibility of the concept, an anthropomorphic thorax phantom containing a porcine heart was scanned on a prototype table-top system. The reconstructed axial images demonstrate a 45 cm in-plane FOV using a 23 cm source array.

Scatter correction method for X-ray CT using primary modulation: theory and preliminary results.

An X-ray system with a large area detector has high scatter-to-primary ratios (SPRs), which result in severe artifacts in reconstructed computed tomography (CT) images. A scatter correction algorithm is introduced that provides effective scatter correction but does not require additional patient exposure. The key hypothesis of the algorithm is that the high-frequency components of the X-ray spatial distribution do not result in strong high-frequency signals in the scatter. A calibration sheet with a checkerboard pattern of semitransparent blockers (a "primary modulator") is inserted between the X-ray source and the object. The primary distribution is partially modulated by a high-frequency function, while the scatter distribution still has dominant low-frequency components, based on the hypothesis. Filtering and demodulation techniques suffice to extract the low-frequency components of the primary and hence obtain the scatter estimation. The hypothesis was validated using Monte Carlo (MC) simulation, and the algorithm was evaluated by both MC simulations and physical experiments. Reconstructions of a software humanoid phantom suggested system parameters in the physical implementation and showed that the proposed method reduced the relative mean square error of the reconstructed image in the central region of interest from 74.2% to below 1%. In preliminary physical experiments on the standard evaluation phantom, this error was reduced from 31.8% to 2.3%, and it was also demonstrated that the algorithm has no noticeable impact on the resolution of the reconstructed image in spite of the filter-based approach. Although the proposed scatter correction technique was implemented for X-ray CT, it can also be used in other X-ray imaging applications, as long as a primary modulator can be inserted between the X-ray source and the imaged object.

A prototype table-top inverse-geometry volumetric CT system.

A table-top volumetric CT system has been implemented that is able to image a 5-cm-thick volume in one circular scan with no cone-beam artifacts. The prototype inverse-geometry CT (IGCT) scanner consists of a large-area, scanned x-ray source and a detector array that is smaller in the transverse direction. The IGCT geometry provides sufficient volumetric sampling because the source and detector have the same axial, or slice direction, extent. This paper describes the implementation of the table-top IGCT scanner, which is based on the NexRay Scanning-Beam Digital X-ray system (NexRay, Inc., Los Gatos, CA) and an investigation of the system performance. The alignment and flat-field calibration procedures are described, along with a summary of the reconstruction algorithm. The resolution and noise performance of the prototype IGCT system are studied through experiments and further supported by analytical predictions and simulations. To study the presence of cone-beam artifacts, a "Defrise" phantom was scanned on both the prototype IGCT scanner and a micro CT system with a +/-5 cone angle for a 4.5-cm volume thickness. Images of inner ear specimens are presented and compared to those from clinical CT systems. Results showed that the prototype IGCT system has a 0.25-mm isotropic resolution and that noise comparable to that from a clinical scanner with equivalent spatial resolution is achievable. The measured MTF and noise values agreed reasonably well with theoretical predictions and computer simulations. The IGCT system was able to faithfully reconstruct the laminated pattern of the Defrise phantom while the micro CT system suffered severe cone-beam artifacts for the same object. The inner ear acquisition verified that the IGCT system can image a complex anatomical object, and the resulting images exhibited more high-resolution details than the clinical CT acquisition. Overall, the successful implementation of the prototype system supports the IGCT concept for single-rotation volumetric scanning free from cone-beam artifacts.