3D Image Reconstruction From Truncated Helical Cone Beam Projection Data - A Linear Prediction Approach

With the introduction of 2D flat-panel X-ray detectors, 3D image reconstruction using helical cone-beam tomography is fast replacing the conventional 2D reconstruction techniques. In 3D image reconstruction, the source orbit or scanning geometry should satisfy the data sufficiency or completeness condition for exact reconstruction. The helical scan geometry satisfies this condition and hence can give exact reconstruction. The theoretically exact helical cone-beam reconstruction algorithm proposed by Katsevich is a breakthrough and has attracted interest in the 3D reconstruction using helical cone-beam Computed Tomography.In many practical situations, the available projection data is incomplete. One such case is where the detector plane does not completely cover the full extent of the object being imaged in lateral direction resulting in truncated projections. This result in artifacts that mask small features near to the periphery of the ROI when reconstructed using the convolution back projection (CBP) method assuming that the projection data is complete. A number of techniques exist which deal with completion of missing data followed by the CBP reconstruction. In 2D, linear prediction (LP)extrapolation has been shown to be efficient for data completion, involving minimal assumptions on the nature of the data, producing smooth extensions of the missing projection data.In this paper, we propose to extend the LP approach for extrapolating helical cone beam truncated data. The projection on the multi row flat panel detectors has missing columns towards either ends in the lateral direction in truncated data situation. The available data from each detector row is modeled using a linear predictor. The available data is extrapolated and this completed projection data is backprojected using the Katsevich algorithm. Simulation results show the efficacy of the proposed method.

Reconstruction from incomplete data in cone-beam tomography

A method for reconstruction of an object f(x) x=(x,y,z) from a limited set of cone-beam projection data has been developed. This method uses a modified form of convolution back-projection and projection onto convex sets (POCS) for handling the limited (or incomplete) data problem. In cone-beam tomography, one needs to have a complete geometry to completely reconstruct the original three-dimensional object. While complete geometries do exist, they are of little use in practical implementations. The most common trajectory used in practical scanners is circular, which is incomplete. It is, however, possible to recover some of the information of the original signal f(x) based on a priori knowledge of the nature of f(x). If this knowledge can be posed in a convex set framework, then POCS can be utilized. In this report, we utilize this a priori knowledge as convex set constraints to reconstruct f(x) using POCS. While we demonstrate the effectiveness of our algorithm for circular trajectories, it is essentially geometry independent and will be useful in any limited-view cone-beam reconstruction.

Image reconstruction with laterally truncated projections in helical cone-beam CT: Linear prediction based projection completion techniques

Lateral or transaxial truncation of cone-beam data can occur either due to the field of view limitation of the scanning apparatus or iregion-of-interest tomography. In this paper, we Suggest two new methods to handle lateral truncation in helical scan CT. It is seen that reconstruction with laterally truncated projection data, assuming it to be complete, gives severe artifacts which even penetrates into the field of view. A row-by-row data completion approach using linear prediction is introduced for helical scan truncated data. An extension of this technique known as windowed linear prediction approach is introduced. Efficacy of the two techniques are shown using simulation with standard phantoms. A quantitative image quality measure of the resulting reconstructed images are used to evaluate the performance of the proposed methods against an extension of a standard existing technique.

Estimation of missing data using windowed linear prediction in laterally truncated projection & in cone beam CT

In this paper, we address the reconstruction problem from laterally truncated helical cone-beam projections. The reconstruction problem from lateral truncation, though similar to that of interior radon problem, is slightly different from it as well as the local (lambda) tomography and pseudo-local tomography in the sense that we aim to reconstruct the entire object being scanned from a region-of-interest (ROI) scan data. The method proposed in this paper is a projection data completion approach followed by the use of any standard accurate FBP type reconstruction algorithm. In particular, we explore a windowed linear prediction (WLP) approach for data completion and compare the quality of reconstruction with the linear prediction (LP) technique proposed earlier.

Correlation of Calcification on Excised Aortic Valves by Micro-Computed Tomography with Severity of Aortic Stenosis

Background and aim of the study: The quantification of incidentally found aortic valve calcification on computed tomography (CT) is not performed routinely, as data relating to the accuracy of aortic valve calcium for estimating the severity of aortic stenosis (AS) is neither consistent nor validated. As aortic valve calcium quantification by CT is confounded by wall and coronary ostial calcification, as well as motion artifact, the ex-vivo micro-computed tomography (micro-CT) of stenotic aortic valves allows a precise measurement of the amounts of calcium present. The study aim, using excised aortic valves from patients with confirmed AS, was to determine if the amount of calcium on micro-CT correlated with the severity of AS. Methods: Each of 35 aortic valves that had been excised from patients during surgical valve replacement were examined using micro-CT imaging. The amount of calcium present was determined by absolute and proportional values of calcium volume in the specimen. Subsequently, the correlation between calcium volume and preoperative mean aortic valve gradient (MAVG), peak transaortic velocity (V-max), and aortic valve area (AVA) on echocardiography, was evaluated. Results: The mean calcium volume across all valves was 603.2 +/- 398.5 mm(3), and the mean ratio of calcium volume to total valve volume was 0.36 +/- 0.16. The mean aortic valve gradient correlated positively with both calcium volume and ratio (r = 0.72, p <0.001). V-max also correlated positively with the calcium volume and ratio (r = 0.69 and 0.76 respectively; p <0.001). A logarithmic curvilinear model proved to be the best fit to the correlation. A calcium volume of 480 mm(3) showed sensitivity and specificity of 0.76 and 0.83, respectively, for a diagnosis of severe AS, while a calcium ratio of 0.37 yielded sensitivity and specificity of 0.82 and 0.94, respectively. Conclusion: A radiological estimation of calcium amount by volume, and its proportion to the total valve volume, were shown to serve as good predictive parameters for severe AS. An estimation of the calcium volume may serve as a complementary measure for determining the severity of AS when aortic valve calcification is identified on CT imaging. The Journal of Heart Valve Disease 2012;21:320-327

Biocompatibility property of 100% strontium-substituted SiO2-Al2O3-P2O5-CaO-CaF2 glass ceramics over 26 weeks implantation in rabbit model: Histology and micro-Computed Tomography analysis

One of the desired properties for any new biomaterial composition is its long-term stability in a suitable animal model and such property cannot be appropriately assessed by performing short-term implantation studies. While hydroxyapatite (HA) or bioglass coated metallic biomaterials are being investigated for in vivo biocompatibility properties, such study is not extensively being pursued for bulk glass ceramics. In view of their inherent brittle nature, the implant stability as well as impact of long-term release of metallic ions on bone regeneration have been a major concern. In this perspective, the present article reports the results of the in vivo implantation experiments carried out using 100% strontium (Sr)-substituted glass ceramics with the nominal composition of 4.5 SiO2-3Al(2)O(3)-1.5P(2)O(5)-3SrO-2SrF(2) for 26 weeks in cylindrical bone defects in rabbit model. The combination of histological and micro-computed tomography analysis provided a qualitative and quantitative understanding of the bone regeneration around the glass ceramic implants in comparison to the highly bioactive HA bioglass implants (control). The sequential polychrome labeling of bone during in vivo osseointegration using three fluorochromes followed by fluorescence microscopy observation confirmed homogeneous bone formation around the test implants. The results of the present study unequivocally confirm the long-term implant stability as well as osteoconductive property of 100% Sr-substituted glass ceramics, which is comparable to that of a known bioactive implant, that is, HA-based bioglass. (c) 2014 Wiley Periodicals, Inc. J Biomed Mater Res Part B: Appl Biomater, 103B: 1168-1179, 2015.

A Method for Delineation of Bone Surfaces in Photoacoustic Computed Tomography of the Finger

Photoacoustic (PA) imaging of interphalangeal peripheral joints is of interest in the context of using the synovial membrane as a surrogate marker of rheumatoid arthritis. Previous work has shown that ultrasound (US) produced by absorption of light at the epidermis reflects on the bone surfaces within the finger. When the reflected signals are backprojected in the region of interest, artifacts are produced, confounding interpretation of the images. In this work, we present an approach where the PA signals known to originate from the epidermis are treated as virtual US transmitters, and a separate reconstruction is performed as in US reflection imaging. This allows us to identify the bone surfaces. Furthermore, the identification of the joint space is important as this provides a landmark to localize a region-of-interest in seeking the inflamed synovial membrane. The ability to delineate bone surfaces allows us to identify not only the artifacts but also the interphalangeal joint space without recourse to new US hardware or a new measurement. We test the approach on phantoms and on a healthy human finger.

Macropores generated during shrinkage in two paddy soils using X-ray micro-computed tomography

Soil shrinkage curve represents a decrease of total porosity or an increase of bulk density with water loss. However, our knowledge of the dynamics of pores and their geometry during soil shrinkage is scarce, partially due to lack of reliable methods for determining soil pores in relation to change in soil water. This study aimed to investigate the dynamics of macropores (>30 mu m) of paddy soils during shrinkage. Two, paddy soils, which were sampled from one paddy field cultivated for 20 years (YPF) and the other one for over 100 years (OPF), represented difference in crack geometry in the field. Macropore parameters (volume, connectivity, and orientation of pores) and soil shrinkage parameters were determined on the same undisturbed soil cores by X-ray microtomography and shrinkage curve, respectively. Macroporosity was on average four times larger in the YPF than in the OPF whereas the shrinkage capacity was lower in the YPF as compared to the OPF (0.09 vs. 0.15 COLE). Soil shrinkage increased the volume of pores by 3.7% in the YPF and by 1.6% in the OPF as well as their connectivity. The formation of macropores occurred mostly in the proportional shrinkage phase. As a result, the slope of the proportional shrinkage phase was smaller in the YPF (0.65) than in the OPF (0.89). New macropores were cracks and extended pre-existing pores in the range of 225-1215 pm size without any preferential orientation. This work provides image evidences that in paddy soils with high shrinkage capacity more macropores are generated in the soil presenting a smaller proportional shrinkage slope. (C) 2015 Elsevier B.V. All rights reserved.

Helical FDK Algorithms With No Backprojection Weight

In this paper, we present two new filtered backprojection (FBP) type algorithms for cylindrical detector helical cone-beam geometry with no position dependent backprojection weight. The algorithms are extension of the recent exact Hilbert filtering based 2D divergent beam reconstruction with no backprojection weight to the FDK type algorithm for reconstruction in 3D helical trajectory cone-beam tomography. The two algorithms named HFDK-W1 and HFDK-W2 result in better image quality, noise uniformity, lower noise and reduced cone-beam artifacts.

A MatLAB Based Virtual Phantom for 2D Electrical Impedance Tomography (MatVP2DEIT): Studying the Medical Electrical Impedance Tomography Reconstruction in Computer

Electrical Impedance Tomography (EIT) is a computerized medical imaging technique which reconstructs the electrical impedance images of a domain under test from the boundary voltage-current data measured by an EIT electronic instrumentation using an image reconstruction algorithm. Being a computed tomography technique, EIT injects a constant current to the patient's body through the surface electrodes surrounding the domain to be imaged (Omega) and tries to calculate the spatial distribution of electrical conductivity or resistivity of the closed conducting domain using the potentials developed at the domain boundary (partial derivative Omega). Practical phantoms are essentially required to study, test and calibrate a medical EIT system for certifying the system before applying it on patients for diagnostic imaging. Therefore, the EIT phantoms are essentially required to generate boundary data for studying and assessing the instrumentation and inverse solvers a in EIT. For proper assessment of an inverse solver of a 2D EIT system, a perfect 2D practical phantom is required. As the practical phantoms are the assemblies of the objects with 3D geometries, the developing of a practical 2D-phantom is a great challenge and therefore, the boundary data generated from the practical phantoms with 3D geometry are found inappropriate for assessing a 2D inverse solver. Furthermore, the boundary data errors contributed by the instrumentation are also difficult to separate from the errors developed by the 3D phantoms. Hence, the errorless boundary data are found essential to assess the inverse solver in 2D EIT. In this direction, a MatLAB-based Virtual Phantom for 2D EIT (MatVP2DEIT) is developed to generate accurate boundary data for assessing the 2D-EIT inverse solvers and the image reconstruction accuracy. MatVP2DEIT is a MatLAB-based computer program which simulates a phantom in computer and generates the boundary potential data as the outputs by using the combinations of different phantom parameters as the inputs to the program. Phantom diameter, inhomogeneity geometry (shape, size and position), number of inhomogeneities, applied current magnitude, background resistivity, inhomogeneity resistivity all are set as the phantom variables which are provided as the input parameters to the MatVP2DEIT for simulating different phantom configurations. A constant current injection is simulated at the phantom boundary with different current injection protocols and boundary potential data are calculated. Boundary data sets are generated with different phantom configurations obtained with the different combinations of the phantom variables and the resistivity images are reconstructed using EIDORS. Boundary data of the virtual phantoms, containing inhomogeneities with complex geometries, are also generated for different current injection patterns using MatVP2DEIT and the resistivity imaging is studied. The effect of regularization method on the image reconstruction is also studied with the data generated by MatVP2DEIT. Resistivity images are evaluated by studying the resistivity parameters and contrast parameters estimated from the elemental resistivity profiles of the reconstructed phantom domain. Results show that the MatVP2DEIT generates accurate boundary data for different types of single or multiple objects which are efficient and accurate enough to reconstruct the resistivity images in EIDORS. The spatial resolution studies show that, the resistivity imaging conducted with the boundary data generated by MatVP2DEIT with 2048 elements, can reconstruct two circular inhomogeneities placed with a minimum distance (boundary to boundary) of 2 mm. It is also observed that, in MatVP2DEIT with 2048 elements, the boundary data generated for a phantom with a circular inhomogeneity of a diameter less than 7% of that of the phantom domain can produce resistivity images in EIDORS with a 1968 element mesh. Results also show that the MatVP2DEIT accurately generates the boundary data for neighbouring, opposite reference and trigonometric current patterns which are very suitable for resistivity reconstruction studies. MatVP2DEIT generated data are also found suitable for studying the effect of the different regularization methods on reconstruction process. Comparing the reconstructed image with an original geometry made in MatVP2DEIT, it would be easier to study the resistivity imaging procedures as well as the inverse solver performance. Using the proposed MatVP2DEIT software with modified domains, the cross sectional anatomy of a number of body parts can be simulated in PC and the impedance image reconstruction of human anatomy can be studied.

Fully 3D refraction correction dosimetry system

The irradiation of selective regions in a polymer gel dosimeter results in an increase in optical density and refractive index (RI) at those regions. An optical tomography-based dosimeter depends on rayline path through the dosimeter to estimate and reconstruct the dose distribution. The refraction of light passing through a dose region results in artefacts in the reconstructed images. These refraction errors are dependant on the scanning geometry and collection optics. We developed a fully 3D image reconstruction algorithm, algebraic reconstruction technique-refraction correction (ART-rc) that corrects for the refractive index mismatches present in a gel dosimeter scanner not only at the boundary, but also for any rayline refraction due to multiple dose regions inside the dosimeter. In this study, simulation and experimental studies have been carried out to reconstruct a 3D dose volume using 2D CCD measurements taken for various views. The study also focuses on the effectiveness of using different refractive-index matching media surrounding the gel dosimeter. Since the optical density is assumed to be low for a dosimeter, the filtered backprojection is routinely used for reconstruction. We carry out the reconstructions using conventional algebraic reconstruction (ART) and refractive index corrected ART (ART-rc) algorithms. The reconstructions based on FDK algorithm for cone-beam tomography has also been carried out for comparison. Line scanners and point detectors, are used to obtain reconstructions plane by plane. The rays passing through dose region with a RI mismatch does not reach the detector in the same plane depending on the angle of incidence and RI. In the fully 3D scanning setup using 2D array detectors, light rays that undergo refraction are still collected and hence can still be accounted for in the reconstruction algorithm. It is found that, for the central region of the dosimeter, the usable radius using ART-rc algorithm with water as RI matched medium is 71.8%, an increase of 6.4% compared to that achieved using conventional ART algorithm. Smaller diameter dosimeters are scanned with dry air scanning by using a wide-angle lens that collects refracted light. The images reconstructed using cone beam geometry is seen to deteriorate in some planes as those regions are not scanned. Refraction correction is important and needs to be taken in to consideration to achieve quantitatively accurate dose reconstructions. Refraction modeling is crucial in array based scanners as it is not possible to identify refracted rays in the sinogram space.

Structure and deformation correlation of closed-cell aluminium foam subject to uniaxial compression

We report the results of an experimental and numerical study conducted on a closed-cell aluminium foam that was subjected to uniaxial compression with lateral constraint. X-ray computed tomography was utilized to gain access into the three-dimensional (3-D) structure of the foam and some aspects of the deformation mechanisms. A series of advanced 3-D image analyses are conducted on the 3-D images aimed at characterizing the strain localization regions. We identify the morphological/geometrical features that are responsible for the collapse of the cells and the strain localization. A novel mathematical approach based on a Minkowski tensor analysis along with the mean intercept length technique were utilized to search for signatures of anisotropy across the foam sample and its evolution as a function of loading. Our results show that regions with higher degrees of anisotropy in the undeformed foam have a tendency to initiate the onset of cell collapse. Furthermore, we show that strain hardening occurs predominantly in regions with large cells and high anisotropy. We combine the finite element method with the tomographic images to simulate the mechanical response of the foam. We predict further deformation in regions where the foam is already deformed. Crown Copyright (C) 2012 Published by Elsevier Ltd. on behalf of Acta Materialia Inc. All rights reserved.

Cryogenically cured hydroxyapatite-gelatin nanobiocomposite for bovine serum albumin protein adsorption and release

This research paper presents the first results on the protein adsorption and release kinetics and in vitro biodegradability of cryogenically cured hydroxyapatite-gelatin based micro/macroporous scaffolds (CHAMPS). While the adsorption and release of bovine serum albumin (BSA) protein exhibits steady state behavior over an incubation period of up to 10 days, Fourier transform infrared (FT-IR) analysis importantly confirms the absence of any change in the secondary structure of BSA proteins due to interaction with the CHAMPS scaffold. The compression properties of the CHAMPS scaffold with interconnected porosity (pore size similar to 50-200 mm) is characterized by a non-linear stress-strain response with a strength close to 5 MPa and a maximum strain of up to 24%. The slow but systematic increase in weight loss over a period of 7 days as well as apatite layer formation indicates its good bioactivity. The extensive micro-computed tomography (micro-CT) analysis establishes cancellous bone-like highly interconnected and complex porous architecture of the CHAMPS scaffold. Importantly, the excellent adsorption (up to 50%) and release (up to 60% of adsorbed protein) of BSA has been uniquely attributed to the inherent porous microstructure of the CHAMPS scaffold. Overall, the present study provides an assessment of the interaction of protein with the gelatin-hydroxyapatite macroporous scaffold in vitro, as well as reporting for the first time the efficacy of such scaffolds to release 60% of BSA loaded onto the scaffold in vitro, which is significantly higher than earlier literature reports.

Early osseointegration of a strontium containing glass ceramic in a rabbit model

The most important property of a bone cement or a bone substitute in load bearing orthopaedic implants is good integration with host bone with reduced bone resorption and increased bone regeneration at the implant interface. Long term implantation of metal-based joint replacements often results in corrosion and particle release, initiating chronic inflammation leading onto osteoporosis of host bone. An alternative solution is the coating of metal implants with hydroxyapatite (HA) or bioglass or the use of bulk bioglass or HA-based composites. In the above perspective, the present study reports the in vivo biocompatibility and bone healing of the strontium (Sr)-stabilized bulk glass ceramics with the nominal composition of 4.5SiO(2)-3Al(2)O(3)-1.5P(2)O(5)-3SrO-2SrF(2) during short term implantation of up to 12 weeks in rabbit animal model. The progression of healing and bone regeneration was qualitatively and quantitatively assessed using fluorescence microscopy, histological analysis and micro-computed tomography. The overall assessment of the present study establishes that the investigated glass ceramic is biocompatible in vivo with regards to local effects after short term implantation in rabbit animal model. Excellent healing was observed, which is comparable to that seen in response to a commercially available implant of HA-based bioglass alone. (C) 2013 Elsevier Ltd. All rights reserved.