3D cone-beam CT reconstruction for circular trajectories

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.

Density estimations and comparisons of a fragmented single fiber using X-ray computed tomography

A commercial X-ray computed tomography (CT) apparatus using a quasi-monochromatic beam was utilized for density estimations and comparisons of a fragmented single fiber. The validation of quasi-monochromaticity of the X-ray source was investigated by radiograph measurements. For the case of a transmittance higher than 50%, the contribution of Cu Kα characteristic X-rays was dominant. To realize a sufficient statistical quality, an attempt to increase the number of averaged voxels was demonstrated using the neighboring slices of the 3D-CT image. A minimum value of the coefficient of variation (CV) was achieved using multiple images rather than using a single image. The observed values of the inverse of the transmitted X-ray intensity (CT value) of the polymers showed a fairly good relationship with their density. An analytical curve derived from measurements of reference samples of known densities could provide the relative density of an unknown fragmented fiber down to the size of 30 μm in diameter and 35 μm in length. The CV of the estimated density was from 1.5 to 2%, which was estimated from the CV of CT values. Moreover, the correlation of CT values was improved with the linear absorption coefficient than the density. A better performance of discrimination of polymers including fibers might be realized with the difference of linear absorption coefficients for X-rays.

Artifacts in optical projection tomography due to refractive-index mismatch: model and correction

Optical projection tomography (OPT) is a powerful tool for three-dimensional (3D) imaging of mesoscopic samples. While it is able to achieve resolution of a few tens of microns over a sample volume of several cubic centimeters, the reconstructed images often suffer from artifacts caused by inaccurate calibration. In this work, we focus on the refractive-index mismatch between the sample and the surrounding medium. We derive a 3D cone-beam forward model of OPT that approximates the effect of refractive-index mismatch. We then implement a fast and efficient reconstruction method to correct for the induced seagull-shaped artifacts on experimental images of fluorescent beads.

Reciprocal-FDK reconstruction for x-ray diffraction computed tomography

Objective. X-ray diffraction (XRD) technology uses x-ray small-angle scattering signal for material analysis, which is highly sensitive to material inter-molecular structure. To meet the high spatial resolution requirement in applications such as medical imaging, XRD computed tomography (XRDCT) has been proposed to provide XRD intensity with improved spatial resolution from point-wise XRD scan. In XRDCT, 2D spatial tomography corresponds to a 3D reconstruction problem with the third dimension being the XRD spectrum dimension, i.e. the momentum transfer dimension. Current works in the field have studied reconstruction methods for either angular-dispersive XRDCT or energy-dispersive XRDCT for small samples. The approximations used are only suitable for regions near the XRDCT iso-center. A new XRDCT reconstruction method is needed for more general imaging applications. Approach. We derive a new FDK-type reconstruction method (Reciprocal-FDK) for XRDCT without limitation on object size. By introducing a set of reciprocal variables, the XRDCT model is transformed into a classical cone-parallel CT model, which is an extension of a circular-trajectory cone-beam CT model, after which the FDK method is applied for XRDCT reconstruction. Main results. Both analytical simulation and Monte Carlo simulation experiments are conducted to validate the XRDCT reconstruction method. The results show that when compared to existing analytical reconstruction methods, there are improvements in the proposed Reciprocal-FDK method with regard to relative structure reconstruction and XRD pattern peak reconstruction. Since cone-parallel CT does not satisfy the data completeness condition, cone-angle effect affects the reconstruction accuracy of XRDCT. The property of cone-angle effect in XRDCT is also analyzed with ablation studies. Significance. We propose a general analytical reconstruction method for XRDCT without constraint on object size. Reciprocal-FDK provides a complete derivation and theoretical support for XRDCT reconstruction by analogy to the well-studied cone-parallel CT model. In addition, the intrinsic problem with the XRDCT data model and the corresponding reconstruction error are discussed for the first time.

Non-destructive analysis of hollow-shaped single fibers using X-ray computed tomography

A commercial high-resolution X-ray computed tomography (CT) was utilized for non-destructive analysis of single fibers. The micro-CT apparatus was employed because it is applicable to both colored and colorless fibers. A sample preparation using adhesive sheets was demonstrated, and the method is similar to typical tape-lift sample collection method in crime cases. Different cross-sectional shapes of nylon and polyester single fibers were non-destructively distinguished, and the method is applicable to all types of fibers. Cross-sectional areas, aperture ratios, and volumes of individual fibers were directly and automatically measured using the open-source software. The observed parameters were within a coefficient of variation of 3%. In addition, a mass of a single fragment of a fiber can be estimated when the local density is given.

Feasibility Studies of X-Ray Computed Tomography for Forensic Examination of Single Fibers

Non-destructive observations of cross-sectional shapes of acrylic single fibers using X-ray computed tomography (CT) were performed. A commercial X-ray CT apparatus (Rigaku nano3DX) was employed because of its micrometer-scale spatial resolution and remarkable image contrast for soft matter. We assessed two types of sample support, a paraffin strip and a nylon string, for single fiber samples in terms of easy handling and sample recovery. Fixed individual single fibers were loaded into a narrow polyimide tube in both cases, and the tube ensured that the sample remained in the field of view during the CT measurements. In both cases, the cross-sectional shapes of individual single fibers could be distinguished, with a circular shape for one sample and a triangular shape for the other. However, the support using a nylon string was found to be more suitable for further analysis. The cross-sectional profile of the obtained tomographic image showed a clear difference between polyimide and nylon. The intensity ratio or the image contrast corresponded to that of the local densities. It was also found that the effect of the artifact appeared at around the boundaries of the objects, but the local density could still be utilized for examining individual single fibers.

Suppression of Cone-Beam Artefacts with Direct Iterative Reconstruction Computed Tomography Trajectories (DIRECTT)

The reconstruction of cone-beam computed tomography data using filtered back-projection algorithms unavoidably results in severe artefacts. We describe how the Direct Iterative Reconstruction of Computed Tomography Trajectories (DIRECTT) algorithm can be combined with a model of the artefacts for the reconstruction of such data. The implementation of DIRECTT results in reconstructed volumes of superior quality compared to the conventional algorithms.

Efficient high cone-angle artifact reduction in circular cone-beam CT using deep learning with geometry-aware dimension reduction

High cone-angle artifacts (HCAAs) appear frequently in circular cone-beam computed tomography (CBCT) images and can heavily affect diagnosis and treatment planning. To reduce HCAAs in CBCT scans, we propose a novel deep learning approach that reduces the three-dimensional (3D) nature of HCAAs to two-dimensional (2D) problems in an efficient way. Specifically, we exploit the relationship between HCAAs and the rotational scanning geometry by training a convolutional neural network (CNN) using image slices that were radially sampled from CBCT scans. We evaluated this novel approach using a dataset of input CBCT scans affected by HCAAs and high-quality artifact-free target CBCT scans. Two different CNN architectures were employed, namely U-Net and a mixed-scale dense CNN (MS-D Net). The artifact reduction performance of the proposed approach was compared to that of a Cartesian slice-based artifact reduction deep learning approach in which a CNN was trained to remove the HCAAs from Cartesian slices. In addition, all processed CBCT scans were segmented to investigate the impact of HCAAs reduction on the quality of CBCT image segmentation. We demonstrate that the proposed deep learning approach with geometry-aware dimension reduction greatly reduces HCAAs in CBCT scans and outperforms the Cartesian slice-based deep learning approach. Moreover, the proposed artifact reduction approach markedly improves the accuracy of the subsequent segmentation task compared to the Cartesian slice-based workflow.

Multi-pass approach to reduce cone-beam artifacts in a circular orbit cone-beam CT system

We propose a multi-pass approach to reduce cone-beam artifacts in a circular orbit cone-beam computed tomography (CT) system. Employing a large 2D detector array reduces the scan time but produces cone-beam artifacts in the Feldkamp, Davis, and Kress (FDK) reconstruction because of insufficient sampling for exact reconstruction. While the two-pass algorithm proposed by Hsieh is effective at reducing cone-beam artifacts, the correction performance is degraded when the bone density is moderate and the cone angle is large. In this work, we treated the cone-beam artifacts generated from bone and soft tissue as if they were from less dense bone objects and corrected them iteratively. The proposed method was validated using a numerical Defrise phantom, XCAT phantom data, and experimental data from a pediatric phantom followed by image quality assessment for FDK, the two-pass algorithm, the proposed method, and the total variation minimization-based iterative reconstruction (TV-IR). The results show that the proposed method was superior to the two-pass algorithm in cone-beam artifact reduction and effectively reduced the overcorrection by the two-pass algorithm near bone regions. It can also be observed that the proposed method produced better correction performance with fewer iterations than the TV-IR algorithm. A qualitative evaluation with mean-squared error, structural similarity, and structural dissimilarity demonstrated the effectiveness of the proposed method.

Multicolour imaging with spectral photon-counting CT: a phantom study

Background

To evaluate the feasibility of multicolour quantitative imaging with spectral photon-counting computed tomography (SPCCT) of different mixed contrast agents.

Methods

Phantoms containing eleven tubes with mixtures of varying proportions of two contrast agents (i.e. two selected from gadolinium, iodine or gold nanoparticles) were prepared so that the attenuation of each tube was about 280 HU. Scans were acquired at 120 kVp and 100 mAs using a five-bin preclinical SPCCT prototype, generating conventional, water, iodine, gadolinium and gold images. The correlation between prepared and measured concentrations was assessed using linear regression. The cross-contamination was measured for each material as the root mean square error (RMSE) of its concentration in the other material images, where no signal was expected. The contrast-to-noise ratio (CNR) relative to a phosphate buffered saline tube was calculated for each contrast agent.

Results

The solutions had similar attenuations (279 ± 10 HU, mean ± standard deviation) and could not be differentiated on conventional images. However, a distinction was observed in the material images within the same samples, and the measured and prepared concentrations were strongly correlated (R² ≥ 0.97, 0.81 ≤ slope ≤ 0.95, -0.68 ≤ offset ≤ 0.89 mg/mL). Cross-contamination in the iodine images for the mixture of gold and gadolinium contrast agents (RMSE = 0.34 mg/mL) was observed. CNR for 1 mg/mL of contrast agent was better for the mixture of iodine and gadolinium (CNR_iodine = 3.20, CNR_gadolinium = 2.80) than gold and gadolinium (CNR_gadolinium = 1.67, CNR_gold = 1.37).

Conclusions

SPCCT enables multicolour quantitative imaging. As a result, it should be possible to perform imaging of multiple uptake phases of a given tissue/organ within a single scan by injecting different contrast agents sequentially.

Dual scan CT image recovery from truncated projections

There are computerized tomography (CT) scanners available commercially for imaging small objects and they are often categorized as mini-CT X-ray machines. One major limitation of these machines is their inability to scan large objects with good image quality because of the truncation of projection data. An algorithm is proposed in this work which enables such machines to scan large objects while maintaining the quality of the recovered image.

The Effect of Surgical Technique and Spacer Texture on Bone Regeneration: A Caprine Study Using the Masquelet Technique

BACKGROUND

The Masquelet-induced-membrane technique is a commonly used method for treating segmental bone defects. However, there are no established clinical standards for management of the induced membrane before grafting.

QUESTIONS/PURPOSES

Two clinically based theories were tested in a chronic caprine tibial defect model: (1) a textured spacer that increases the induced-membrane surface area will increase bone regeneration; and (2) surgical scraping to remove a thin tissue layer of the inner induced-membrane surface will enhance bone formation.

METHODS

Thirty-two skeletally mature female goats were assigned to four groups: smooth spacer with or without membrane scraping and textured spacer with or without membrane scraping. During an initial surgical procedure (unilateral, left tibia), a defect was created excising bone (5 cm), periosteum (9 cm), and muscle (10 g). Segments initially were stabilized with an intramedullary rod and an antibiotic-impregnated polymethylmethacrylate spacer with a smooth or textured surface. Four weeks later, the spacer was removed and the induced-membrane was either scraped or left intact before bone grafting. Bone formation was assessed using micro-CT (total bone volume in 2.5-cm central defect region) as the primary outcome; radiographs and histologic analysis as secondary outcomes, with the reviewer blinded to the treatment groups of the samples being assessed 12 weeks after grafting. All statistical tests were performed using a linear mixed effects model approach.

RESULTS

Micro-CT analysis showed greater bone formation in defects with scraped induced membrane (mean, 3034.5 mm³; median, 1928.0 mm³; quartile [Q]1-Q3, 273.3-2921.1 mm³) compared with defects with intact induced membrane (mean, 1709.5 mm³; median, 473.8 mm³; Q1-Q3, 132.2-1272.3 mm³; p = 0.034). There was no difference in bone formation between textured spacers (mean, 2405.5 mm³; median, 772.7 mm³; Q1-Q3, 195.9-2743.8 mm³) and smooth spacers (mean, 2473.2 mm³; median, 1143.6 mm³; Q1-Q3, 230.2-451.1 mm³; p = 0.917).

CONCLUSIONS

Scraping the induced-membrane surface to remove the innermost layer of the induced-membrane increased bone regeneration. A textured spacer that increased the induced-membrane surface area had no effect on bone regeneration.

CLINICAL RELEVANCE

Scraping the induced membrane during the second stage of the Masquelet technique may be a rapid and simple means of improving healing of segmental bone defects, which needs to be confirmed clinically.

Light field reconstruction from projection modeling of focal stack

This paper aims to reconstruct the object-side light field from the focal stack focusing on different imaging planes. In the forward problem, the focal stack was modeled as projections of the light field. Based on projection modeling, both the filtered-backprojection(FBP) method and the Landweber iterative scheme of solving the inverse problem regarding light field reconstruction from focal stack were derived by applying the methods of computerized tomography(CT). The experimental results show that the high-precision light field can be reconstructed via FBP and Simultaneous-Algebraic-Reconstruction-Technique(SART) algorithm, and depth and surface of the scene can be reconstructed via the reconstructed light field.

Study of Image Qualities From 6D Robot-Based CBCT Imaging System of Small Animal Irradiator

Purpose:

To assess the quality of cone beam computed tomography images obtained by a robotic arm-based and image-guided small animal conformal radiation therapy device.

Method and Materials:

The small animal conformal radiation therapy device is equipped with a 40 to 225 kV X-ray tube mounted on a custom made gantry, a 1024 × 1024 pixels flat panel detector (200 μm resolution), a programmable 6 degrees of freedom robot for cone beam computed tomography imaging and conformal delivery of radiation doses. A series of 2-dimensional radiographic projection images were recorded in cone beam mode by placing and rotating microcomputed tomography phantoms on the “palm’ of the robotic arm. Reconstructed images were studied for image quality (spatial resolution, image uniformity, computed tomography number linearity, voxel noise, and artifacts).

Results:

Geometric accuracy was measured to be 2% corresponding to 0.7 mm accuracy on a Shelley microcomputed tomography QA phantom. Qualitative resolution of reconstructed axial computed tomography slices using the resolution coils was within 200 μm. Quantitative spatial resolution was found to be 3.16 lp/mm. Uniformity of the system was measured within 34 Hounsfield unit on a QRM microcomputed tomography water phantom. Computed tomography numbers measured using the linearity plate were linear with material density (R² > 0.995). Cone beam computed tomography images of the QRM multidisk phantom had minimal artifacts.

Conclusion:

Results showed that the small animal conformal radiation therapy device is capable of producing high-quality cone beam computed tomography images for precise and conformal small animal dose delivery. With its high-caliber imaging capabilities, the small animal conformal radiation therapy device is a powerful tool for small animal research.

[Evaluation for Basic Image Qualities Dependence on the Position in XYZ Directions and Acquisition Parameters of the Cone Beam CT for Angiography System with Flat Panel Detector]

The purpose of this study was to investigate the effect of the position in XYZ directions and acquisition parameters on the basic image qualities of for cone beam computed tomography (CBCT) in an angiography system with flat panel detector. The resolution property (modulation transfer function: MTF) and the noise property (Wiener spectrum: WS) of CBCT images in X-Y plane were measured with different acquisition parameters (scan matrix number and projection number) and the effect of the position in XYZ directions. The MTFs with 1024×1024 matrix were higher than those of 512×512 matrix and decreased in the peripheral areas due to the reduction of projection number. The highest and the lowest MTFs were measured at the X-ray tube side and on the detector side of the position in X-Y plane, respectively. The WS-doubled projection number showed about 50% lesser noise level. There were differences in the Wiener spectra (WS) at the position in XYZ directions. We conclude that the resolution and the noise property of CBCT image in X-Y plane showed dependences on the position in XYZ directions and acquisition parameters of the CBCT.

A temporal interpolation approach for dynamic reconstruction in perfusion CT

This article presents a dynamic CT reconstruction algorithm for objects with time dependent attenuation coefficient. Projection data acquired over several rotations are interpreted as samples of a continuous signal. Based on this idea, a temporal interpolation approach is proposed which provides the maximum temporal resolution for a given rotational speed of the CT scanner. Interpolation is performed using polynomial splines. The algorithm can be adapted to slow signals, reducing the amount of data acquired and the computational cost. A theoretical analysis of the approximations made by the algorithm is provided. In simulation studies, the temporal interpolation approach is compared with three other dynamic reconstruction algorithms based on linear regression, linear interpolation, and generalized Parker weighting. The presented algorithm exhibits the highest temporal resolution for a given sampling interval. Hence, our approach needs less input data to achieve a certain quality in the reconstruction than the other algorithms discussed or, equivalently, less x-ray exposure and computational complexity. The proposed algorithm additionally allows the possibility of using slow rotating scanners for perfusion imaging purposes.

Interesting detector shapes for third generation CT scanners

PURPOSE

To discuss options in designing detector shapes in third generation CT and to quantify potential cost savings for compact third generation CT systems, and to extend the work from two-dimensional fan-beam CT to three-dimensional cone-beam CT for circular, sequential, and spiral scan trajectories.

METHODS

Third generation CT scanners typically comprise detectors which are flat or whose shape is the segment of a cylinder or a sphere that is focused onto the focal spot of the x-ray source. There appear to be two design criteria that favor this choice of detector shape. One is the possibility of performing fan-beam and cone-beam filtered backprojection in the native geometry (without rebinning) and the other criterion could be to enable the early use of focused antiscatter grids. It is less known, however, that other detector shapes may also have these properties. While these designs have been evaluated for 2D CT from a theoretical standpoint more than one decade ago the authors revisit and generalize these considerations, extend them to 3D circular, sequential, and spiral cone-beam CT and propose an optimal design in terms of detector costs while keeping image quality constant. Their considerations and conclusions are based on considering the sampling density of the x-rays, including the effects of finite focal spot and finite detector element size. Proposing image reconstruction algorithms or numerically evaluating the results by reconstructing simulated projection data is not within the scope of this work.

RESULTS

If the detector arc is curved to be nearly concentric with the circle describing the edge of the field of measurement significantly less detector area and detector pixels are required compared to today's third generation CT systems where the detector arc is centered about the focal spot. Combined with a detector that just covers the spiral Tam window cost savings of 60% or more are possible in compact CT systems. In terms of practicability the new designs appear to be nearly as easy to realize as today's third generation systems.

CONCLUSIONS

Compact CT systems, which require the focal spot to be mounted close to the edge of the field of measurement, may significantly benefit from using detector shapes other than the typical equiangular detector that is focused onto the focal spot.

Generation of voxelized breast phantoms from surgical mastectomy specimens

PURPOSE

In the research and development of dedicated tomographic breast imaging systems, digital breast object models, also known as digital phantoms, are useful tools. While various digital breast phantoms do exist, the purpose of this study was to develop a realistic high-resolution model suitable for simulating three-dimensional (3D) breast imaging modalities. The primary goal was to design a model capable of producing simulations with realistic breast tissue structure.

METHODS

The methodology for generating an ensemble of digital breast phantoms was based on imaging surgical mastectomy specimens using a benchtop, cone-beam computed tomography system. This approach allowed low-noise, high-resolution projection views of the mastectomy specimens at each angular position. Reconstructions of these projection sets were processed using correction techniques and diffusion filtering prior to segmentation into breast tissue types in order to generate phantoms.

RESULTS

Eight compressed digital phantoms and 20 uncompressed phantoms from which an additional 96 pseudocompressed digital phantoms with voxel dimensions of 0.2 mm(3) were generated. Two distinct tissue classification models were used in forming breast phantoms. The binary model classified each tissue voxel as either adipose or fibroglandular. A multivalue scaled model classified each tissue voxel as percentage of adipose tissue (range 1%-99%). Power spectral analysis was performed to compare simulated reconstructions using the breast phantoms to the original breast specimen reconstruction, and fits were observed to be similar.

CONCLUSIONS

The digital breast phantoms developed herein provide a high-resolution anthropomorphic model of the 3D uncompressed and compressed breast that are suitable for use in evaluating and optimizing tomographic breast imaging modalities. The authors believe that other research groups might find the phantoms useful, and therefore they offer to make them available for wider use.

FDK-Type Algorithms with No Backprojection Weight for Circular and Helical Scan CT

We develop two Feldkamp-type reconstruction algorithms with no backprojection weight for circular and helical trajectory with planar detector geometry. Advances in solid-state electronic detector technologies lend importance to CT systems with the equispaced linear array, the planar (flat panel) detectors, and the corresponding algorithms. We derive two exact Hilbert filtered backprojection (FBP) reconstruction algorithms with no backprojection weight for 2D fan-beam equispace linear array detector geometry (complement of the equi-angular curved array detector). Based on these algorithms, the Feldkamp-type algorithms with no backprojection weight for 3D reconstruction are developed using the standard heuristic extension of the divergent beam FBP algorithm. The simulation results show that the axial intensity drop in the reconstructed image using the FDK algorithms with no backprojection weight with circular trajectory is similar to that obtained by using Hu's and T-FDK, algorithms. Further, we present efficient algorithms to reduce the axial intensity drop encountered in the standard FDK reconstructions in circular cone-beam CT. The proposed algorithms consist of mainly two steps: reconstruction of the object using FDK algorithm with no backprojection weight and estimation of the missing term. The efficient algorithms are compared with the FDK algorithm, Hu's algorithm, T-FDK, and Zhu et al.'s algorithm in terms of axial intensity drop and noise. Simulation shows that the efficient algorithms give similar performance in axial intensity drop as that of Zhu et al.'s algorithm while one of the efficient algorithms outperforms Zhu et al.'s algorithm in terms of computational complexity.

Improved compressed sensing-based cone-beam CT reconstruction using adaptive prior image constraints

Volumetric cone-beam CT (CBCT) images are acquired repeatedly during a course of radiation therapy and a natural question to ask is whether CBCT images obtained earlier in the process can be utilized as prior knowledge to reduce patient imaging dose in subsequent scans. The purpose of this work is to develop an adaptive prior image constrained compressed sensing (APICCS) method to solve this problem. Reconstructed images using full projections are taken on the first day of radiation therapy treatment and are used as prior images. The subsequent scans are acquired using a protocol of sparse projections. In the proposed APICCS algorithm, the prior images are utilized as an initial guess and are incorporated into the objective function in the compressed sensing (CS)-based iterative reconstruction process. Furthermore, the prior information is employed to detect any possible mismatched regions between the prior and current images for improved reconstruction. For this purpose, the prior images and the reconstructed images are classified into three anatomical regions: air, soft tissue and bone. Mismatched regions are identified by local differences of the corresponding groups in the two classified sets of images. A distance transformation is then introduced to convert the information into an adaptive voxel-dependent relaxation map. In constructing the relaxation map, the matched regions (unchanged anatomy) between the prior and current images are assigned with smaller weight values, which are translated into less influence on the CS iterative reconstruction process. On the other hand, the mismatched regions (changed anatomy) are associated with larger values and the regions are updated more by the new projection data, thus avoiding any possible adverse effects of prior images. The APICCS approach was systematically assessed by using patient data acquired under standard and low-dose protocols for qualitative and quantitative comparisons. The APICCS method provides an effective way for us to enhance the image quality at the matched regions between the prior and current images compared to the existing PICCS algorithm. Compared to the current CBCT imaging protocols, the APICCS algorithm allows an imaging dose reduction of 10-40 times due to the greatly reduced number of projections and lower x-ray tube current level coming from the low-dose protocol.

Virtual temporal bone dissection system: OSU virtual temporal bone system: development and testing

OBJECTIVES/HYPOTHESIS

The objective of this project was to develop a virtual temporal bone dissection system that would provide an enhanced educational experience for the training of otologic surgeons.

STUDY DESIGN

A randomized, controlled, multi-institutional, single-blinded validation study.

METHODS

The project encompassed four areas of emphasis: structural data acquisition, integration of the system, dissemination of the system, and validation.

RESULTS

Structural acquisition was performed on multiple imaging platforms. Integration achieved a cost-effective system. Dissemination was achieved on different levels including casual interest, downloading of software, and full involvement in development and validation studies. A validation study was performed at eight different training institutions across the country using a two-arm randomized trial where study subjects were randomized to a 2-week practice session using either the virtual temporal bone or standard cadaveric temporal bones. Eighty subjects were enrolled and randomized to one of the two treatment arms; 65 completed the study. There was no difference between the two groups using a blinded rating tool to assess performance after training.

CONCLUSIONS

A virtual temporal bone dissection system has been developed and compared to cadaveric temporal bones for practice using a multicenter trial. There was no statistical difference between practice on the current simulator compared to practice on human cadaveric temporal bones. Further refinements in structural acquisition and interface design have been identified, which can be implemented prior to full incorporation into training programs and used for objective skills assessment.

Inter-plane artifact suppression in tomosynthesis using 3D CT image data

BACKGROUND

Despite its superb lateral resolution, flat-panel-detector (FPD) based tomosynthesis suffers from low contrast and inter-plane artifacts caused by incomplete cancellation of the projection components stemming from outside the focal plane. The incomplete cancellation of the projection components, mostly due to the limited scan angle in the conventional tomosynthesis scan geometry, often makes the image contrast too low to differentiate the malignant tissues from the background tissues with confidence.

METHODS

In this paper, we propose a new method to suppress the inter-plane artifacts in FPD-based tomosynthesis. If 3D whole volume CT images are available before the tomosynthesis scan, the CT image data can be incorporated into the tomosynthesis image reconstruction to suppress the inter-plane artifacts, hence, improving the image contrast. In the proposed technique, the projection components stemming from outside the region-of-interest (ROI) are subtracted from the measured tomosynthesis projection data to suppress the inter-plane artifacts. The projection components stemming from outside the ROI are calculated from the 3D whole volume CT images which usually have lower lateral resolution than the tomosynthesis images. The tomosynthesis images are reconstructed from the subtracted projection data which account for the x-ray attenuation through the ROI. After verifying the proposed method by simulation, we have performed both CT scan and tomosynthesis scan on a phantom and a sacrificed rat using a FPD-based micro-CT.

RESULTS

We have measured contrast-to-noise ratio (CNR) from the tomosynthesis images which is an indicator of the residual inter-plane artifacts on the focal-plane image. In both cases of the simulation and experimental imaging studies of the contrast evaluating phantom, CNRs have been significantly improved by the proposed method. In the rat imaging also, we have observed better visual contrast from the tomosynthesis images reconstructed by the proposed method.

CONCLUSIONS

The proposed tomosynthesis technique can improve image contrast with aids of 3D whole volume CT images. Even though local tomosynthesis needs extra 3D CT scanning, it may find clinical applications in special situations in which extra 3D CT scan is already available or allowed.

Cone-beam artifact evaluation of the factorization method

PURPOSE

The authors investigate the CB artifact behavior of the factorization approach recently suggested for image reconstruction in circular cone-beam computed tomography. This investigation is carried out in a typical C-arm geometry and involves simulated data and for the first time also phantom and clinical CB data acquired with a commercially available angiographic system.

METHODS

The CB artifact level is first measured using quantitative figures-of-merit that are computed from the reconstructions of the mathematical FORBILD head phantom and of a modified disk phantom. The authors then show reconstructions from a physical thorax phantom and clinical head data sets for a visual assessment of image quality. The performance of the factorization method is primarily compared to that of short-scan FDK, but the authors also show the results obtained with the full-scan FDK and the virtual PI-line BPF method for the simulation studies, as a benchmark.

RESULTS

Quantitatively, the FORBILD head phantom reconstructions of both FDK methods show a spatially averaged bias of up to 1.2% in the axial slices about 9 cm away from the plane of the scan, which is placed 4 cm below the central slice through the phantom. The artifact level for the short-scan FDK method and the virtual PI-line BPF method noticeably depends on the scan orientation. The factorization approach can significantly reduce both, this dependency as well as the reconstruction bias. It also shows visually an improved quality of the clinical images compared to short-scan FDK, particularly close to the spine and in the subcranial regions of the clinical data sets.

CONCLUSIONS

The factorization approach comes with noticeably lower reconstruction bias than the FDK methods and is least sensitive to the scan orientation among all considered short-scan methods. The data inconsistencies contained in the real data sets, such as scatter, beam hardening, or data truncation, show only little impact on the factorization results. Hence, in both, reconstructions from real and simulated data, the factorization method yields better image quality than short-scan FDK, albeit at the cost of some slight, directed high-frequency artifacts that are mostly visible in axial slices.

Dimensional stability in composite cone beam computed tomography

An automated increase in the field of view (FOV) for multipurpose cone beam CT (CBCT) by "stitching" (joining) up to three component volumes to yield a larger composite volume must still ensure dimensional stability, especially if the image is to form the basis for a surgical splint. Dimensional stability, image discrepancies and the influence of movement artefacts between exposures were evaluated. The first consumer installation of the Kodak 9000 three-dimensional (3D) extraoral imaging system with stitching software was used for the evaluation of a human mandible with three endodontic instruments as markers. The distances between several reproducible points were measured directly and the results compared with the values measured on screen. Displacements of the mandible along all axes between exposures as well as angular displacements were conducted to test the capability of the system. The standard deviations (SD) of the results for the vertical distances varied between 0.212 mm and 0.409 mm (approximately 1-2 voxels; range, 0.6-1.3 mm) and may be considered the systematic error. The SD of the results for the horizontal and diagonal distances varied between 0.195 mm and 0.571 mm (approximately 1-3 voxels; range, 0.6-1.7 mm) if the group with overall horizontal angulations of 10° and a central rotation of 20° was omitted. In conclusion, the evaluated stitching software is a useful tool to expand the options of combined CBCT with an initial small FOV by allowing a merger of up to three component volumes to yield a larger FOV of about 80 × 80 × 37 mm. The dimensional stability was acceptable when seen in relation to the induced disturbance. Further evaluation of this composite CBCT/digital imaging and communications in medicine system for subsequent splint fabrication may yield promising results.

Application of polychromatic µCT for mineral density determination

Accurate assessment of mineral density (MD) provides information critical to the understanding of mineralization processes of calcified tissues, including bones and teeth. High-resolution three-dimensional assessment of the MD of teeth has been demonstrated by relatively inaccessible synchrotron radiation microcomputed tomography (SRµCT). While conventional desktop µCT (CµCT) technology is widely available, polychromatic source and cone-shaped beam geometry confound MD assessment. Recently, considerable attention has been given to optimizing quantitative data from CµCT systems with polychromatic x-ray sources. In this review, we focus on the approaches that minimize inaccuracies arising from beam hardening, in particular, beam filtration during the scan, beam-hardening correction during reconstruction, and mineral density calibration. Filtration along with lowest possible source voltage results in a narrow and near-single-peak spectrum, favoring high contrast and minimal beam-hardening artifacts. More effective beam monochromatization approaches are described. We also examine the significance of beam-hardening correction in determining the accuracy of mineral density estimation. In addition, standards for the calibration of reconstructed grey-scale attenuation values against MD, including K(2)PHO(4) liquid phantom, and polymer-hydroxyapatite (HA) and solid hydroxyapatite (HA) phantoms, are discussed.

Region of Interest Imaging for a General Trajectory with the Rebinned BPF Algorithm

The back-projection-filtration (BPF) algorithm has been applied to image reconstruction for cone-beam configurations with general source trajectories. The BPF algorithm can reconstruct 3-D region-of-interest (ROI) images from data containing truncations. However, like many other existing algorithms for cone-beam configurations, the BPF algorithm involves a back-projection with a spatially varying weighting factor, which can result in the non-uniform noise levels in reconstructed images and increased computation time. In this work, we propose a BPF algorithm to eliminate the spatially varying weighting factor by using a rebinned geometry for a general scanning trajectory. This proposed BPF algorithm has an improved noise property, while retaining the advantages of the original BPF algorithm such as minimum data requirement.

Low-noise dynamic reconstruction for x-ray tomographic perfusion studies using low sampling rates

Functional imaging based on tomographic X-ray imaging relies on the reconstruction of a temporal sequence of images which accurately reproduces the time attenuation curves of the tissue. The main constraints of these techniques are temporal resolution and dose. Using current techniques the data acquisition has to be performed fast so that the dynamic attenuation values can be regarded as static during the scan. Due to the relatively high number of repeated scans the dose per single scan has to be low yielding a poor signal-to-noise ratio (SNR) in the reconstructed images. In a previous publication a temporal interpolation scheme in the projection data space was relaxing the temporal resolution constraint. The aim of this contribution is the improvement of the SNR. A temporal smoothing term is introduced in the temporal interpolation scheme such that only the physiologic relevant bandwidth is considered. A significant increase of the SNR is achieved. The obtained level of noise only depends on the total dose applied and is independent of the number of scans and the SNR of a single reconstructed image. The approach might be the first step towards using slowly rotating CT systems for perfusion imaging like C-arm or small animal CT scanners.

Evaluation of the effects of sagging shifts on isocenter accuracy and image quality of cone-beam CT from kV on-board imagers

To investigate the effects of sagging shifts of three on-board kV imaging systems (OBI) on the isocenter positioning accuracy and image quality of cone-beam CT (CBCT). A cubical phantom having a metal marker in the center that can be aligned with the radiation isocenter was used to measure sagging shifts and their variation with gantry angle on three Varian linacs with kV on-board imaging systems. A marker-tracking algorithm was applied to detect the shadow of the metal marker and localize its center in the two-dimensional cone-beam radiographic projections. This tracking algorithm is based on finding the position of maximum cross-correlation between a region-of-interest from a template image (including the metal marker) and the projections containing the shadow of the metal marker. Sagging shifts were corrected by mapping the center of the metal marker to a reference position for all projections acquired over a full gantry rotation (0-360 degrees). The sag-corrected radiographic projections were then used to reconstruct CBCT using Feldkamp back-projection. A standard quality assurance phantom was used to evaluate the image quality of CBCT before and after sagging correction. Sagging affects both the positioning accuracy of the OBI isocenter and the CBCT image quality. For example, on one linac, the position of the marker on the cone-beam radiographic projections depends on the angular view and has maximal shifts of about 2 mm along the imager x-direction (patient's cross-plane). Sagging produces systematic shifts of the OBI isocenter as large as 1 mm posterior and 1 mm left in patient coordinates relative to the radiation isocenter. Further, it causes spatial distortion and blurring in CBCT image reconstructed from radiographic projections that are not corrected for OBI sagging. CBCT numbers vary by about 1% in full-fan scans and up to 3.5% in half-fan scans because of sagging. In order to achieve better localization accuracy in image-guided radiation therapy, sagging shifts of the kV OBI need to be corrected. In addition, correction of sagging improves image and provides better visualization of internal structures. Frequent quality assurance is required to monitor and maintain standards of variations in the mechanical accuracy of isocenter and image quality of CBCT because of sagging shifts.

[Combinations of scan parameters and image quality at C-arm CT for abdominal imaging]

We measured the contrast-to-noise ratio (CNR) and evaluated low-contrast images and streak artifacts to optimize abdominal C-arm CT imaging, and we investigated the view number, acquisition matrix, and pixel depth. To measure CNR, we filled 0.125-1.0-inch cavities in an American Association of Physicists in Medicine (AAPM) CT performance phantom with a sodium chloride solution. Five radiological technologists visually evaluated the noise, signal conspicuity at low-and high-signal density, and the overall image quality using paired comparisons based on Thurstone's law. In a given acquisition matrix, the total view number had the greatest effect on the image noise, artifacts, and signal detectability on C-arm CT images. For a given incident dose per view on the flat-panel detector (FPD), fewer images with noise and streak artifacts resulted when a larger view number was selected.

Cone-beam technique for 64-MDCT of lung: image quality comparison with stepwise (step-and-shoot) technique

OBJECTIVE

The purpose of this study was to use phantom and patient data acquired with 64-MDCT to compare the image quality and characteristics of helical high-resolution CT images obtained with cone-beam reconstruction with those of stepwise high-resolution CT images obtained with fan-beam reconstruction.

SUBJECTS AND METHODS

We reconstructed helical high-resolution CT images with cone-beam technique and stepwise high-resolution CT images with fan-beam technique. In the phantom study, we measured high-contrast spatial resolution and image noise using a phantom. Streak artifact was evaluated by five radiologists using the phantom. In the clinical phase of the study, two radiologists independently evaluated high-resolution helical and stepwise CT images of the lung fields of 30 patients with diffuse lung disease. Using a 3-point ordinal scale, the radiologists assessed the sharpness of peripheral vessels and interlobular fissures, artifacts, and graininess in the lung fields; overall image quality; and the sharpness of the contour of the left ventricle.

RESULTS

In high-contrast spatial resolution, the contrast curves in each spatial frequency were similar on the helical and stepwise images. In the clinical study, there was no statistically significant difference between helical and stepwise images with respect to sharpness of the contour of the left ventricle, peripheral vessels, or interlobular fissures (p>0.05). With respect to streak artifacts and graininess in the lung fields, helical images received a significantly higher quality grade than did stepwise images (p<0.05).

CONCLUSION

Our phantom and clinical evaluation showed that the quality of high-resolution CT images of the lung obtained with helical scanning was comparable with the quality of stepwise scans.

An efficient estimation method for reducing the axial intensity drop in circular cone-beam CT

Reconstruction algorithms for circular cone-beam (CB) scans have been extensively studied in the literature. Since insufficient data are measured, an exact reconstruction is impossible for such a geometry. If the reconstruction algorithm assumes zeros for the missing data, such as the standard FDK algorithm, a major type of resulting CB artifacts is the intensity drop along the axial direction. Many algorithms have been proposed to improve image quality when faced with this problem of data missing; however, development of an effective and computationally efficient algorithm remains a major challenge. In this work, we propose a novel method for estimating the unmeasured data and reducing the intensity drop artifacts. Each CB projection is analyzed in the Radon space via Grangeat's first derivative. Assuming the CB projection is taken from a parallel beam geometry, we extract those data that reside in the unmeasured region of the Radon space. These data are then used as in a parallel beam geometry to calculate a correction term, which is added together with Hu's correction term to the FDK result to form a final reconstruction. More approximations are then made on the calculation of the additional term, and the final formula is implemented very efficiently. The algorithm performance is evaluated using computer simulations on analytical phantoms. The reconstruction comparison with results using other existing algorithms shows that the proposed algorithm achieves a superior performance on the reduction of axial intensity drop artifacts with a high computation efficiency.

Experimental study of intracranial hematoma detection with flat panel detector C-arm CT

BACKGROUND AND PURPOSE

Intracranial hemorrhage is a commonly acknowledged complication of interventional neuroradiology procedures, and the ability to image hemorrhage at the time of the procedure would be very beneficial. A new C-arm system with 3D functionality extends the capability of C-arm imaging to include soft-tissue applications by facilitating the detection of low-contrast objects. We evaluated its ability to detect small intracranial hematomas in a swine model.

MATERIALS AND METHODS

Intracranial hematomas were created in 7 swine by autologous blood injection of various hematocrits (19%-37%) and volumes (1.5-5 mL). Four animals received intravascular contrast before obtaining autologous blood (group 1), and 3 did not (group 2). We scanned each animal by using the C-arm CT system, acquiring more than 500 images during a 20-second rotation through more than 200 degrees . Multiplanar reformatted images with isotropic resolution were reconstructed on the workstation by using product truncation, scatter, beam-hardening, and ring-artifact correction algorithms. The brains were harvested and sliced for hematoma measurement and compared with imaging findings.

RESULTS

Five intracranial hematomas were created in group 1 animals, and all were visualized. Six were created in group 2, and 3 were visualized. One nonvisualized hematoma was not confirmed at necropsy. All the others in both groups were confirmed. In group 1 (with contrast), small hematomas were detectable even when the hematocrit was 19%-20%. In group 2 (without contrast) C-arm CT was able to detect small hematomas (<1.0 cm(2)) created with hematocrits of 29%-37%. The area of hematoma measured from the C-arm CT data was, on average, within 15% of the area measured from harvested brain.

CONCLUSIONS

The image quality obtained with this implementation of C-arm CT was sufficient to detect experimentally created small intracranial hematomas. This capability should provide earlier detection of hemorrhagic complications that may occur during neurointerventional procedures.

A short-scan reconstruction for cone-beam CT using shift-invariant FBP and equal weighting

A 3D reconstruction formula has been derived for a circular cone-beam (CB) short scan using ID shift-invariant filtering, CB backprojection, and equal weighting. By first converting the divergent projections to parallel projections, we analyze the circular CB data using the classic central slice theorem. The sampling density in Fourier space is investigated and 1D shift-invariant filtering before backprojection can be used to compensate for the nonuniformity. The final formula consists of a conventional FDK reconstruction and a correction term using differential backprojection and the 1D Hilbert transform in the image domain. On a full scan, the approach reduces to the FDK algorithm, while for a short scan, the CB artifacts are suppressed by the second term. This algorithm outperforms the modified FDK algorithm with Parker's weighting, as illustrated by computer simulations and experimental results. Due to its shift-invariant filtered-backprojection structure, the proposed algorithm is implemented efficiently, and requires a simple adaptation of the FDK algorithm.

In Vivo Assessment of Bone Graft/Endplate Contact Pressure in a Caprine Interbody Pseudarthrosis Model: A Preliminary Biomechanical Characterization of the Fusion Process for the Development of a Microelectromechanical Systems (MEMS) Biosensor

Background

In this preliminary study we used a goat model to quantify pressure at an interbody bone graft interface. Although the study was designed to assess fusion status, the concept behind the technology could lead to early detection of implant failure and potential hazardous complications related to motion-preservation devices. The purpose of this study was to investigate the feasibility of in vivo pressure monitoring as a strategy to determine fusion status.

Methods

Telemetric pressure transducers were implanted, and pressure at the bone graft interfaces of cervical interbody fusion autografts placed into living goats (Groups A and B) was evaluated. Group A constituted the 4-month survival group and Group B the 6-month survival group. One goat served as the study control (Group C) and was not implanted with a pressure transducer. An additional six cadaveric goat cervical spines (Group D) were obtained from a local slaughterhouse and implanted with bone grafts and ventral plates and used for in vitro biomechanical comparison to the specimens from Groups A and B.

Results

All goats demonstrated an increase in interface pressure within the first 10 days postoperatively, with the largest relative change in pressure occurring between the sixth and ninth days. The goats from Groups A and B had a 200% to 400% increase in relative pressure.

Conclusions

Although this was a pilot study to assess pressure as an indicator for a fusion or pseudarthrosis, the preliminary data suggest that early bone healing is detectable by an increase in pressure. Thus, pressure may serve as an indicator of fusion status by detecting altered biomechanical parameters.

Iliac bars lever reduction and fixation system used in the treatment of spondylolisthesis

Background

The purpose of the current study was to use the Iliac Bars Lever Reduction and Fixation System (IBLRFS) for Grades 1 and 2 spondylolytic spondylolisthesis, evaluate its stability and reductive efficacy, and examine the complications.

Methods

Between April 2005 and August 2006, 44 patients with Grades 1 and 2 spondylolytic spondylolisthesis were treated surgically: 21 patients underwent posterior Iliac Bars Lever Reduction and Fixation (IBRLFS), 23 patients were treated with traditional stabilization and reduction systems (SRS). The follow-up periods ranged from 1 to 2 years (mean, 1 year and 2 months). The clinical outcome, fusion rate, average percentile degree of displacement, displacement angle, sacral inclination, ratio of intervertebral height, and complications were evaluated. Operating time, blood loss, and duration of hospital stay were compared.

Results

There were no statistically significant differences between the 2 groups in blood loss, recovery rate, and radiographic results. However, there were statistically significant differences in operating time (P < .05), duration of hospital stay (P < .05). There were no cases of nonunion in the two groups. In the IBLRFS group, preoperatively, the average percentile degree of displacement, displacement angle, sacral inclination, and ratio of intervertebral height were 23.48% ± 5.36%, 2.2° ± 1.1°, 29.4° ± 6.5°, and 0.68 ± 0.21, respectively. Postoperatively, the respective measurements were 6.47% ± 1.49%, 10.3° ± 3.3°, 42.6° ± 8.1°, and 0.85 ± 0.12. No patients experienced major complications.

In the SRS group, preoperatively, the average percentile degree of displacement, displacement angle, sacral inclination, and ratio of intervertebral height were 21.78% ± 5.16%, 2.3° ± 1.0°, 26.4° ± 8.5°, and 0.62 ± 0.25, respectively. Postoperatively, the respective measurements were 6.34% ± 2.01%, 9.8° ± 2.1°, 44.1° ± 7.6°, and 0.79 ± 0.23. One patient experienced a badly placed screw in the right pedicle of lumbar 4.

Conclusions

This kind of new fixation system (IBLRFS) was shown to be useful in the treatment of spondylolisthesis, and its use was associated with minimal complications after 14 months of mean follow-up.

Level of Evidence

Therapeutic, case studies (level 4).

A factorization approach for cone-beam reconstruction on a circular short-scan

In this paper, we introduce a new algorithm for 3-D image reconstruction from cone-beam (CB) projections acquired along a partial circular scan. Our algorithm is based on a novel, exact factorization of the initial 3-D reconstruction problem into a set of independent 2-D inversion problems, each of which corresponds to finding the object density on one, single plane. Any such 2-D inversion problem is solved numerically using a projected steepest descent iteration scheme. We present a numerical evaluation of our factorization algorithm using computer-simulated CB data, without and with noise, of the FORBILD head phantom and of a disk phantom. First, we study quantitatively the impact of the reconstruction parameters on the algorithm performance. Next, we present reconstruction results for visual assessment of the achievable image quality and provide, for comparison, results obtained with two other state-of-the-art reconstruction algorithms for the circular short-scan.

Resolution and noise trade-off analysis for volumetric CT

Until recently, most studies addressing the trade-off between spatial resolution and quantum noise were performed in the context of single-slice CT. In this study, we extend the theoretical framework of previous works to volumetric CT and further extend it by taking into account the actual shapes of the preferred reconstruction kernels. In the experimental study, we also attempt to explore a three-dimensional approach for spatial resolution measurement, as opposed to the conventional two-dimensional approaches that were widely adopted in previously published studies. By scanning a finite-sized sphere phantom, the MTF was measured from the edge profile along the spherical surface. Cases of different resolutions (and noise levels) were generated by adjusting the reconstruction kernel. To reduce bias, the total photon fluxes were matched: 120 kVp, 200 mA, and 1 s per gantry rotation. All data sets were reconstructed using a modified FDK algorithm under the same condition: Scan field-of-view (SFOV) = 10 cm, and slice thickness = 0.625 mm. The theoretical analysis indicated that the variance of noise is proportional to > 4th power of the spatial resolution. Our experimental results supported this conclusion by showing the relationship is 4.6th (helical) or 5th (axial) power.

Bone formation following OP-1 implantation is improved by addition of autogenous bone marrow cells in a canine femur defect model

Osteogenic Protein-1 (OP-1, BMP-7) acts locally on connective tissue progenitors (CTPs) to induce bone formation. The response to OP-1 and similar agents is potentially limited by the number of local CTPs. This study tested the hypothesis that supplementing local CTPs using autogenous bone marrow will enhance bone formation at an OP-1 implant. Four 1.0-cm diameter unicortical cylindrical defects in the left proximal femur were grafted in each of seven dogs. Radial ingrowth of new bone formation was assessed at 4 weeks using micro CT. The OP-1 (3.5 mg rhOP-1 in 1 g bovine collagen I matrix) was implanted in each site combined with either clotted blood or aspirated bone marrow (BM). Bone formation was increased in the group augmented with transplanted marrow. These data suggest that increasing the local population of cells and CTPs using aspirated bone marrow can enhance the performance of OP-1, but may not eliminate the effects of site variation on the response to OP-1 and similar agents. The canine multiple femoral defect model defined in this study is well suited to quantitatively evaluate strategies for augmenting bone repair using local cell targeting and cell transplantation strategies.

GPU-based streaming architectures for fast cone-beam CT image reconstruction and demons deformable registration

This paper shows how to significantly accelerate cone-beam CT reconstruction and 3D deformable image registration using the stream-processing model. We describe data-parallel designs for the Feldkamp, Davis and Kress (FDK) reconstruction algorithm, and the demons deformable registration algorithm, suitable for use on a commodity graphics processing unit. The streaming versions of these algorithms are implemented using the Brook programming environment and executed on an NVidia 8800 GPU. Performance results using CT data of a preserved swine lung indicate that the GPU-based implementations of the FDK and demons algorithms achieve a substantial speedup--up to 80 times for FDK and 70 times for demons when compared to an optimized reference implementation on a 2.8 GHz Intel processor. In addition, the accuracy of the GPU-based implementations was found to be excellent. Compared with CPU-based implementations, the RMS differences were less than 0.1 Hounsfield unit for reconstruction and less than 0.1 mm for deformable registration.

Development of a non-linear weighted hybrid cone-beam CT reconstruction for circular trajectories

We investigated an image reconstruction algorithm to reduce cone-beam artifacts in cone-beam CT. Our new algorithm to reduce such artifacts features: (1) a change in weighting with respect to projection data obtained at different projection angles; (2) distribution of correction coefficients so that they are larger near the center of the detector, while taking individual channel data for the detector into account, and smaller near the edges; (3) three-dimensional back-projection of corrected projection data. These findings confirmed that this algorithm reduces cone-beam artifacts and generates high-quality reconstruction images.