Recycling Krylov subspaces for efficient large-scale electrical impedance tomography

Electrical impedance tomography (EIT) captures images of internal features of a body. Electrodes are attached to the boundary of the body, low intensity alternating currents are applied, and the resulting electric potentials are measured. Then, based on the measurements, an estimation algorithm obtains the three-dimensional internal admittivity distribution that corresponds to the image. One of the main goals of medical EIT is to achieve high resolution and an accurate result at low computational cost. However, when the finite element method (FEM) is employed and the corresponding mesh is refined to increase resolution and accuracy, the computational cost increases substantially, especially in the estimation of absolute admittivity distributions. Therefore, we consider in this work a fast iterative solver for the forward problem, which was previously reported in the context of structural optimization. We propose several improvements to this solver to increase its performance in the EIT context. The solver is based on the recycling of approximate invariant subspaces, and it is applied to reduce the EIT computation time for a constant and high resolution finite element mesh. In addition, we consider a powerful preconditioner and provide a detailed pseudocode for the improved iterative solver. The numerical results show the effectiveness of our approach: the proposed algorithm is faster than the preconditioned conjugate gradient (CG) algorithm. The results also show that even on a standard PC without parallelization, a high mesh resolution (more than 150,000 degrees of freedom) can be used for image estimation at a relatively low computational cost. (C) 2010 Elsevier B.V. All rights reserved.

Methodological analysis of gamma tomography system for large random packed columns

Gamma ray tomography experiments have been carried out to detect spatial patterns in the porosity in a 0.27 m diameter column packed with steel Rashig rings of different sizes: 12.6, 37.9, and 76 mm. using a first generation CT system (Chen et al., 1998). A fast Fourier transform tomographic reconstruction algorithm has been used to calculate the spatial variation over the column cross section. Cross-sectional gas porosity and solid holdup distribution were determinate. The values of cross-sectional average gas porosity were epsilon=0.849, 0.938 and 0.966 for the 12.6, 37.9, and 76 mm rings, respectively. Radial holdup variation within the packed bed has been determined. The variation of the circumferentially averaged gas holdup in the radial direction indicates that the porosity in the column wall region is a somewhat higher than that in the bulk region, due to the effect of the column wall. (C) 2009 Elsevier Ltd. All rights reserved.

Three Dimensional Finite Element Analyses of Oral Structures by Computerized Tomography

Our aim was to document the benefits of three dimensional finite element model generations from computed tomography data as well as the realistic creation of all oral structures in a patient. The stresses resulting from the applied load in our study did not exceed the structure limitations, suggesting a clinically acceptable physiological condition.

Prospective ultramorphological characterization of human hair by optical coherence tomography

Background/purpose The continuous advancement in cosmetic science has led to an increasing demand for the development of non-invasive, reliable scientific techniques directed toward claim substantiation, which is of utmost relevance, to obtain data regarding the efficacy and safety of cosmetic products. Methods In this work, we used the optical coherence tomography (OCT) technique to produce in vitro transversal section-images of human hair. We also compared the OCT signal before and after chemical treatment with an 18% w/w ammonium thioglycolate solution. Results The mean diameter of the medulla was 29 +/- 7 mu m and the hair diameter was 122 +/- 16 mu m in our samples of standard Afro-ethnic hair. A three-dimensional (3D) image was constructed starting from 601 cross-sectional images (slices). Each slice was taken in steps of 6.0 mu m at eight frames per second, and the entire 3D image was constructed in 60 s. Conclusion It was possible to identify, using the A-scan protocol, the principal structures: the cuticle, cortex and medulla. After chemical treatment, it was not possible to identify the main structures of hair fiber due to index matching promoted by deleterious action of the chemical agent.

Refractive index tomography of turbid media by bifocal optical coherence refractometry

We demonstrate tomographic imaging of the refractive index of turbid media using bifocal optical coherence refractometry (BOCR). The technique, which is a variant of optical coherence tomography, is based on the measurement of the optical pathlength difference between two foci simultaneously present in a medium of interest. We describe a new method to axially shift the bifocal optical pathlength that avoids the need to physically relocate the objective lens or the sample during an axial scan, and present an experimental realization based on an adaptive liquid-crystal lens. We present experimental results, including video clips, which demonstrate refractive index tomography of a range of turbid liquid phantoms, as well as of human skin in vivo.

X-ray computed tomography to quantify tree rooting spatial distributions

Poor root development due to constraining soil conditions could be an important factor influencing health of urban trees. Therefore, there is a need for efficient techniques to analyze the spatial distribution of tree roots. An analytical procedure for describing tree rooting patterns from X-ray computed tomography (CT) data is described and illustrated. Large irregularly shaped specimens of undisturbed sandy soil were sampled from Various positions around the base of trees using field impregnation with epoxy resin, to stabilize the cohesionless soil. Cores approximately 200 mm in diameter by 500 mm in height were extracted from these specimens. These large core samples were scanned with a medical X-ray CT device, and contiguous images of soil slices (2 mm thick) were thus produced. X-ray CT images are regarded as regularly-spaced sections through the soil although they are not actual 2D sections but matrices of voxels similar to 0.5 mm x 0.5 mm x 2 mm. The images were used to generate the equivalent of horizontal root contact maps from which three-dimensional objects, assumed to be roots, were reconstructed. The resulting connected objects were used to derive indices of the spatial organization of roots, namely: root length distribution, root length density, root growth angle distribution, root spatial distribution, and branching intensity. The successive steps of the method, from sampling to generation of indices of tree root organization, are illustrated through a case study examining rooting patterns of valuable urban trees. (C) 1999 Elsevier Science B.V. All rights reserved.

3D reconstruction and quantification of macropores using X-ray computed tomography and image analysis

Axial X-ray Computed tomography (CT) scanning provides a convenient means of recording the three-dimensional form of soil structure. The technique has been used for nearly two decades, but initial development has concentrated on qualitative description of images. More recently, increasing effort has been put into quantifying the geometry and topology of macropores likely to contribute to preferential now in soils. Here we describe a novel technique for tracing connected macropores in the CT scans. After object extraction, three-dimensional mathematical morphological filters are applied to quantify the reconstructed structure. These filters consist of sequences of so-called erosions and/or dilations of a 32-face structuring element to describe object distances and volumes of influence. The tracing and quantification methodologies were tested on a set of undisturbed soil cores collected in a Swiss pre-alpine meadow, where a new earthworm species (Aporrectodea nocturna) was accidentally introduced. Given the reduced number of samples analysed in this study, the results presented only illustrate the potential of the method to reconstruct and quantify macropores. Our results suggest that the introduction of the new species induced very limited chance to the soil structured for example, no difference in total macropore length or mean diameter was observed. However. in the zone colonised by, the new species. individual macropores tended to have a longer average length. be more vertical and be further apart at some depth. Overall, the approach proved well suited to the analysis of the three-dimensional architecture of macropores. It provides a framework for the analysis of complex structures, which are less satisfactorily observed and described using 2D imaging. (C) 2002 Elsevier Science B.V. All rights reserved.

Quantum process tomography of a controlled-NOT gate

We demonstrate complete characterization of a two-qubit entangling process-a linear optics controlled-NOT gate operating with coincident detection-by quantum process tomography. We use a maximum-likelihood estimation to convert the experimental data into a physical process matrix. The process matrix allows an accurate prediction of the operation of the gate for arbitrary input states and a calculation of gate performance measures such as the average gate fidelity, average purity, and entangling capability of our gate, which are 0.90, 0.83, and 0.73, respectively.

Thoracic aorta calcification detected by electron beam tomography predicts all-cause mortality

Background: The presence of coronary artery calcium (CAC) is an independent marker of increased risk of cardiovascular disease (CVD) events and mortality. However, the predictive value of thoracic aorta calcification (TAC), which can be additionally identified without further scanning during assessment of CAC, is unknown. Methods: We followed a cohort of 8401 asymptomatic individuals (mean age: 53 +/- 10 years, 69% men) undergoing cardiac risk factor evaluation and TAC and CAC testing with electron beam computed tomography. Multivariable Cox proportional hazards models were developed to predict all-cause mortality based on the presence of TAC. Results: During a median follow-up period of 5 years, 124 (1.5%) deaths were observed. Overall survival was 96.9% and 98.9% for those with and without detectable TAC, respectively (p < 0.0001). Compared to those with no TAC, the hazard ratio for mortality in the presence of TAC was 3.25 (95% CI: 2.28-4.65, p < 0.0001) in unadjusted analysis. After adjusting for age, gender, hypertension, dyslipidemia, diabetes mellitus, smoking and family history of premature coronary artery disease, and presence of CAC the relationship remained robust (HR 1.61, 95% CI: 1.10-2.27, p = 0.015). Likelihood ratio chi(2) statistics demonstrated that the addition of TAC contributed significantly in predicting mortality to traditional risk factors alone (chi(2) = 13.62, p = 0.002) as well as risk factors + CAC (chi(2) = 5.84, p = 0.02) models. Conclusion: In conclusion, the presence of TAC was associated with all-cause mortality in our study; this relationship was independent of conventional CVD risk factors as well as the presence of CAC. (C) 2009 Elsevier Ireland Ltd. All rights reserved.

Bedside estimation of recruitable alveolar collapse and hyperdistension by electrical impedance tomography

To present a novel algorithm for estimating recruitable alveolar collapse and hyperdistension based on electrical impedance tomography (EIT) during a decremental positive end-expiratory pressure (PEEP) titration. Technical note with illustrative case reports. Respiratory intensive care unit. Patients with acute respiratory distress syndrome. Lung recruitment and PEEP titration maneuver. Simultaneous acquisition of EIT and X-ray computerized tomography (CT) data. We found good agreement (in terms of amount and spatial location) between the collapse estimated by EIT and CT for all levels of PEEP. The optimal PEEP values detected by EIT for patients 1 and 2 (keeping lung collapse < 10%) were 19 and 17 cmH(2)O, respectively. Although pointing to the same non-dependent lung regions, EIT estimates of hyperdistension represent the functional deterioration of lung units, instead of their anatomical changes, and could not be compared directly with static CT estimates for hyperinflation. We described an EIT-based method for estimating recruitable alveolar collapse at the bedside, pointing out its regional distribution. Additionally, we proposed a measure of lung hyperdistension based on regional lung mechanics.

Characterization of blood-filled cyst by contrast echocardiography and computed tomography

Contrast echocardiography has been demonstrated useful for left ventricular opacification and improvement of endocardial border delineation. Another important clinical application of this technique refers to the better characterization of cardiac tumors and masses. We here described an asymptomatic patient with cystic mass attached to submitral valve apparatus in which contrast echocardiography was performed after intravenous injection of lipid-encapsulated microbubbles. It resulted in enhancement of the cystic borders and allowed for better definition of its diagnosis. Multislice computed tomography confirmed the echocardiographic findings. This case illustrates the potential of contrast echocardiography to improve the anatomic evaluation of cardiac masses.

Three-dimensional electrical impedance, tomography: A topology optimization approach

Electrical impedance tomography is a technique to estimate the impedance distribution within a domain, based on measurements on its boundary. In other words, given the mathematical model of the domain, its geometry and boundary conditions, a nonlinear inverse problem of estimating the electric impedance distribution can be solved. Several impedance estimation algorithms have been proposed to solve this problem. In this paper, we present a three-dimensional algorithm, based on the topology optimization method, as an alternative. A sequence of linear programming problems, allowing for constraints, is solved utilizing this method. In each iteration, the finite element method provides the electric potential field within the model of the domain. An electrode model is also proposed (thus, increasing the accuracy of the finite element results). The algorithm is tested using numerically simulated data and also experimental data, and absolute resistivity values are obtained. These results, corresponding to phantoms with two different conductive materials, exhibit relatively well-defined boundaries between them, and show that this is a practical and potentially useful technique to be applied to monitor lung aeration, including the possibility of imaging a pneumothorax.

Relation of Aortic Valve Calcium Detected by Cardiac Computed Tomography to All-Cause Mortality

Aortic valve calcium (AVC) can be quantified on the same computed tomographic scan as coronary artery calcium (CAC). Although CAC is an established predictor of cardiovascular events, limited evidence is available for an independent predictive value for AVC. We studied a cohort of 8,401 asymptomatic subjects (mean age 53 10 years, 69% men), who were free of known coronary heart disease and were undergoing electron beam computed tomography for assessment of subclinical atherosclerosis. The patients were followed for a median of 5 years (range 1 to 7) for the occurrence of mortality from any cause. Multivariate Cox regression models were developed to predict all-cause mortality according to the presence of AVC. A total of 517 patients (6%) had AVC on electron beam computed tomography. During follow-up, 124 patients died (1.5%), for an overall survival rate of 96.1% and 98.7% for those with and without AVC, respectively (hazard ratio 3.39, 95% confidence interval 2.09 to 5.49). After adjustment for age, gender, hypertension, dyslipidemia, diabetes mellitus, smoking, and a family history of premature coronary heart disease, AVC remained a significant predictor of mortality (hazard ratio 1.82, 95% confidence interval 1.11 to 2.98). Likelihood ratio chi-square statistics demonstrated that the addition of AVC contributed significantly to the prediction of mortality in a model adjusted for traditional risk factors (chi-square = 5.03, p = 0.03) as well as traditional risk factors plus the presence of CAC (chi-square = 3.58, p = 0.05). In conclusion, AVC was associated with increased all-cause mortality, independent of the traditional risk factors and the presence of CAC. (C) 2010 Published by Elsevier Inc. (Am J Cardiol 2010;106:1787-1791)

Dipyridamole Stress and Rest Myocardial Perfusion by 64-Detector Row Computed Tomography in Patients With Suspected Coronary Artery Disease

Recently, stress myocardial computed tomographic perfusion (CTP) was shown to detect myocardial ischemia. Our main objective was to evaluate the feasibility of dipyridamole stress CTP and compare it to single-photon emission computed tomography (SPECT) to detect significant coronary stenosis using invasive conventional coronary angiography (CCA; stenosis >70%) as the reference method. Thirty-six patients (62 +/- 8 years old, 20 men) with previous positive results with SPECT (<2 months) as the primary inclusion criterion and suspected coronary artery disease underwent a customized multidetector-row CT protocol with myocardial perfusion evaluation at rest and during stress and coronary CT angiography (CTA). Multidetector-row computed tomography was performed in a 64-slice scanner with dipyridamole stress perfusion acquisition before a second perfusion/CT angiographic acquisition at rest. Independent blinded observers performed analysis of images from CTP, CTA, and CCA. All 36 patients completed the CT protocol with no adverse events (mean radiation dose 14.7 +/- 3.0 mSv) and with interpretable scans. CTP results were positive in 27 of 36 patients (75%). From the 9 (25%) disagreements, 6 patients had normal coronary arteries and 2 had no significant stenosis (8 false-positive results with SPECT, 22%). The remaining patient had an occluded artery with collateral flow confirmed by conventional coronary angiogram. Good agreement was demonstrated between CTP and SPECT on a per-patient analysis (kappa 0.53). In 26 patients using CCA as reference, sensitivity, specificity, and positive and negative predictive values were 88.0%, 79.3%, 66.7%, and 93.3% for CTP and 68.8, 76.1%, 66.7%, and 77.8%, for SPECT, respectively (p = NS). In conclusion, dipyridamole CT myocardial perfusion at rest and during stress is feasible and results are similar to single-photon emission CT scintigraphy. The anatomical-perfusion information provided by this combined CT protocol may allow identification of false-positive results by SPECT. (C) 2010 Elsevier Inc. All rights reserved. (Am J Cardiol 2010;106:310-315)

A computerized tomography scan method for calculating the hernia sac and abdominal cavity volume in complex large incisional hernia with loss of domain

Preoperative progressive pneumoperitoneum (PPP) is a safe and effective procedure in the treatment of large incisional hernia (size > 10 cm in width or length) with loss of domain (LIHLD). There is no consensus in the literature on the amount of gas that must be insufflated in a PPP program or even how long it should be maintained. We describe a technique for calculating the hernia sac volume (HSV) and abdominal cavity volume (ACV) based on abdominal computerized tomography (ACT) scanning that eliminates the need for subjective criteria for inclusion in a PPP program and shows the amount of gas that must be insufflated into the abdominal cavity in the PPP program. Our technique is indicated for all patients with large or recurrent incisional hernias evaluated by a senior surgeon with suspected LIHLD. We reviewed our experience from 2001 to 2008 of 23 consecutive hernia surgical procedures of LIHLD undergoing preoperative evaluation with CT scanning and PPP. An ACT was required in all patients with suspected LIHLD in order to determine HSV and ACV. The PPP was performed only if the volume ratio HSV/ACV (VR = HSV/ACV) was a parts per thousand yen25% (VR a parts per thousand yen 25%). We have performed this procedure on 23 patients, with a mean age of 55.6 years (range 31-83). There were 16 women and 7 men with an average age of 55.6 years (range 31-83), and a mean BMI of 38.5 kg/m(2) (range 23-55.2). Almost all patients (21 of 23 patients-91.30%) were overweight; 43.5% (10 patients) were severely obese (obese class III). The mean calculated volumes for ACV and HSV were 9,410 ml (range 6,060-19,230 ml) and 4,500 ml (range 1,850-6,600 ml), respectively. The PPP is performed by permanent catheter placed in a minor surgical procedure. The total amount of CO(2) insufflated ranged from 2,000 to 7,000 ml (mean 4,000 ml). Patients required a mean of 10 PPP sessions (range 4-18) to achieve the desired volume of gas (that is the same volume that was calculated for the hernia sac). Since PPP sessions were performed once a day, 4-18 days were needed for preoperative preparation with PPP. The mean VR was 36% (ranged from 26 to 73%). We conclude that ACT provides objective data for volume calculation of both hernia sac and abdominal cavity and also for estimation of the volume of gas that should be insufflated into the abdominal cavity in PPP.