Exploring the knowledge contained in neuroimages: Statistical discriminant analysis and automatic segmentation of the most significant changes

Objective: The aim of this article is to propose an integrated framework for extracting and describing patterns of disorders from medical images using a combination of linear discriminant analysis and active contour models. Methods: A multivariate statistical methodology was first used to identify the most discriminating hyperplane separating two groups of images (from healthy controls and patients with schizophrenia) contained in the input data. After this, the present work makes explicit the differences found by the multivariate statistical method by subtracting the discriminant models of controls and patients, weighted by the pooled variance between the two groups. A variational level-set technique was used to segment clusters of these differences. We obtain a label of each anatomical change using the Talairach atlas. Results: In this work all the data was analysed simultaneously rather than assuming a priori regions of interest. As a consequence of this, by using active contour models, we were able to obtain regions of interest that were emergent from the data. The results were evaluated using, as gold standard, well-known facts about the neuroanatomical changes related to schizophrenia. Most of the items in the gold standard was covered in our result set. Conclusions: We argue that such investigation provides a suitable framework for characterising the high complexity of magnetic resonance images in schizophrenia as the results obtained indicate a high sensitivity rate with respect to the gold standard. (C) 2010 Elsevier B.V. All rights reserved.

Pulmonary function testing in quadriplegic subjects

Study design: Cross-sectional study. Objective: Pulmonary functional capacity in 23 Brazilian quadriplegic subjects (ASIA A), aged 30 (9.5) years, weight 66 (10.75) kg, height 176 (7) cm, was investigated at 42 ( 64) months postinjury. Setting: University Hospital-UNICAMP, Campinas, Brazil. Method: Subjects performed forced vital capacity ( FVC) and maximal voluntary ventilation (MVV) tests while seated in their standard wheelchairs. Forced Expired Volume after 1 s (FEV1) and FVC/FEV1 ratio were calculated from these tests. Values obtained were compared to three prediction equations from the literature that are used specifically for spinal cord subjects and include different variables in their formulae, such as age, gender, height, postinjury time and injury level. Data are expressed as median (interquartile interval). Differences between values were demonstrated by median confidence interval with significance level set at a 0.05. Results: Obtained data were statistically different from prediction equation results, with FVC 3.11 ( 0.81), 4.46 (0.28), 4.16 (0.33), 4.26 (0.42); FEV1 2.77 (1.03), 3.67 (0.21), 3.66 (0.30), 3.45 (0.39) and MVV 92 (27), 154.2 (11.9), 156.6 (14),157.3 (16.8), where the first value is obtained experimentally and the second, third and fourth values correspond to predicted values. The results obtained from spirometry test in this study differed significantly from the results obtained when prediction equations were used. Conclusion: The use of prediction equations developed to estimate pulmonary function in wheelchair users significantly overestimates pulmonary function of quadriplegic individuals with complete lesions (ASIA group A), in comparison to measured values.

Human Fatalities from Cyanobacteria: Chemical and Biological Evidence for Cyanotoxins

An outbreak of acute liver failure occurred at a dialysis center in Caruaru, Brazil (8 degrees 17 'S, 35 degrees 58 'W), 134 km from Recife, the state capital of Pernambuco. At the clinic, 116 (89%) of 131 patients experienced visual disturbances, nausea, and vomiting after routine hemodialysis treatment on 13-20 February 1996. Subsequently, 100 patients developed acute liver failure, and of these 76 died. As of December 1996, 52 of the deaths could be attributed to a common syndrome now called Caruaru syndrome. Examination of phytoplankton from the dialysis clinic's water source, analyses of the clinic's water treatment system, plus serum and liver tissue of clinic patients led to the identification of two groups of cyanobacterial toxins, the hepatotoxic cyclic peptide microcystins and the hepatotoxic alkaloid cylindrospermopsin. Comparison of victims' symptoms and pathology using animal studies of these two cyanotoxins leads us to conclude that the major contributing factor to death of the dialyses patients was intravenous exposure to microcystins, specifically microcystin-YR, -LR, and -AR. From liver concentrations and exposure volumes, it was estimated that 19.5 mug/L microcystin was in the water used for dialysis treatments. This is 19.5 times the level set as a guideline for safe drinking water supplies by the World. Health Organization.

Brief note on the variation of constants formula for fuzzy differential equations

This note gives a theory of state transition matrices for linear systems of fuzzy differential equations. This is used to give a fuzzy version of the classical variation of constants formula. A simple example of a time-independent control system is used to illustrate the methods. While similar problems to the crisp case arise for time-dependent systems, in time-independent cases the calculations are elementary solutions of eigenvalue-eigenvector problems. In particular, for nonnegative or nonpositive matrices, the problems at each level set, can easily be solved in MATLAB to give the level sets of the fuzzy solution. (C) 2002 Elsevier Science B.V. All rights reserved.

PAH air pollution at a Portuguese urban area: carcinogenic risks and sources identification

This study aimed to characterize air pollution and the associated carcinogenic risks of polycyclic aromatic hydrocarbon (PAHs) at an urban site, to identify possible emission sources of PAHs using several statistical methodologies, and to analyze the influence of other air pollutants and meteorological variables on PAH concentrations.The air quality and meteorological data were collected in Oporto, the second largest city of Portugal. Eighteen PAHs (the 16 PAHs considered by United States Environment Protection Agency (USEPA) as priority pollutants, dibenzo[a,l]pyrene, and benzo[j]fluoranthene) were collected daily for 24 h in air (gas phase and in particles) during 40 consecutive days in November and December 2008 by constant low-flow samplers and using polytetrafluoroethylene (PTFE) membrane filters for particulate (PM10 and PM2.5 bound) PAHs and pre-cleaned polyurethane foam plugs for gaseous compounds. The other monitored air pollutants were SO2, PM10, NO2, CO, and O3; the meteorological variables were temperature, relative humidity, wind speed, total precipitation, and solar radiation. Benzo[a]pyrene reached a mean concentration of 2.02 ngm−3, surpassing the EU annual limit value. The target carcinogenic risks were equal than the health-based guideline level set by USEPA (10−6) at the studied site, with the cancer risks of eight PAHs reaching senior levels of 9.98×10−7 in PM10 and 1.06×10−6 in air. The applied statistical methods, correlation matrix, cluster analysis, and principal component analysis, were in agreement in the grouping of the PAHs. The groups were formed according to their chemical structure (number of rings), phase distribution, and emission sources. PAH diagnostic ratios were also calculated to evaluate the main emission sources. Diesel vehicular emissions were the major source of PAHs at the studied site. Besides that source, emissions from residential heating and oil refinery were identified to contribute to PAH levels at the respective area. Additionally, principal component regression indicated that SO2, NO2, PM10, CO, and solar radiation had positive correlation with PAHs concentrations, while O3, temperature, relative humidity, and wind speed were negatively correlated.

Numerische Simulation dendritischen Kristallwachstums

Dendritic Growth, Stefan-Problem, Finite-Element-Method, Level-Set-Method

Blood Pressure and Hemodynamic Adaptations after a Training Program in Young Individuals with Down Syndrome

Background:Cardiovascular diseases affect people worldwide. Individuals with Down Syndrome (DS) have an up to sixteen-time greater risk of mortality from cardiovascular diseases.Objective:To evaluate the effects of aerobic and resistance exercises on blood pressure and hemodynamic variables of young individuals with DS.Methods:A total of 29 young individuals with DS participated in the study. They were divided into two groups: aerobic training (AT) (n = 14), and resistance training (TR) (n = 15). Their mean age was 15.7 ± 2.82 years. The training program lasted 12 weeks, and had a frequency of three times a week for AT and twice a week for RT. AT was performed in treadmill/ bicycle ergometer, at an intensity between 50%-70% of the HR reserve. RT comprised nine exercises with three sets of 12 repetition-maximum. Systolic blood pressure (SBP), diastolic blood pressure (DBP), mean blood pressure (MBP) and hemodynamic variables were assessed beat-to-beat using the Finometer device before/after the training program. Descriptive analysis, the Shapiro-Wilk test to check the normality of data, and the two-way ANOVA for repeated measures were used to compare pre- and post-training variables. The Pearson’s correlation coefficient was calculated to correlate hemodynamic variables. The SPSS version 18.0 was used with the significance level set at p < 0.05.Results:After twelve weeks of aerobic and/or resistance training, significant reductions in variables SBP, DBP and MBP were observed.Conclusion:This study suggests a chronic hypotensive effect of moderate aerobic and resistance exercises on young individuals with DS.

Central and Cortical Gray Mater Segmentation of Magnetic Resonance Images of the Fetal Brain

Motivation. The study of human brain development in itsearly stage is today possible thanks to in vivo fetalmagnetic resonance imaging (MRI) techniques. Aquantitative analysis of fetal cortical surfacerepresents a new approach which can be used as a markerof the cerebral maturation (as gyration) and also forstudying central nervous system pathologies [1]. However,this quantitative approach is a major challenge forseveral reasons. First, movement of the fetus inside theamniotic cavity requires very fast MRI sequences tominimize motion artifacts, resulting in a poor spatialresolution and/or lower SNR. Second, due to the ongoingmyelination and cortical maturation, the appearance ofthe developing brain differs very much from thehomogenous tissue types found in adults. Third, due tolow resolution, fetal MR images considerably suffer ofpartial volume (PV) effect, sometimes in large areas.Today extensive efforts are made to deal with thereconstruction of high resolution 3D fetal volumes[2,3,4] to cope with intra-volume motion and low SNR.However, few studies exist related to the automatedsegmentation of MR fetal imaging. [5] and [6] work on thesegmentation of specific areas of the fetal brain such asposterior fossa, brainstem or germinal matrix. Firstattempt for automated brain tissue segmentation has beenpresented in [7] and in our previous work [8]. Bothmethods apply the Expectation-Maximization Markov RandomField (EM-MRF) framework but contrary to [7] we do notneed from any anatomical atlas prior. Data set &Methods. Prenatal MR imaging was performed with a 1-Tsystem (GE Medical Systems, Milwaukee) using single shotfast spin echo (ssFSE) sequences (TR 7000 ms, TE 180 ms,FOV 40 x 40 cm, slice thickness 5.4mm, in plane spatialresolution 1.09mm). Each fetus has 6 axial volumes(around 15 slices per volume), each of them acquired inabout 1 min. Each volume is shifted by 1 mm with respectto the previous one. Gestational age (GA) ranges from 29to 32 weeks. Mother is under sedation. Each volume ismanually segmented to extract fetal brain fromsurrounding maternal tissues. Then, in-homogeneityintensity correction is performed using [9] and linearintensity normalization is performed to have intensityvalues that range from 0 to 255. Note that due tointra-tissue variability of developing brain someintensity variability still remains. For each fetus, ahigh spatial resolution image of isotropic voxel size of1.09 mm is created applying [2] and using B-splines forthe scattered data interpolation [10] (see Fig. 1). Then,basal ganglia (BS) segmentation is performed on thissuper reconstructed volume. Active contour framework witha Level Set (LS) implementation is used. Our LS follows aslightly different formulation from well-known Chan-Vese[11] formulation. In our case, the LS evolves forcing themean of the inside of the curve to be the mean intensityof basal ganglia. Moreover, we add local spatial priorthrough a probabilistic map created by fitting anellipsoid onto the basal ganglia region. Some userinteraction is needed to set the mean intensity of BG(green dots in Fig. 2) and the initial fitting points forthe probabilistic prior map (blue points in Fig. 2). Oncebasal ganglia are removed from the image, brain tissuesegmentation is performed as described in [8]. Results.The case study presented here has 29 weeks of GA. Thehigh resolution reconstructed volume is presented in Fig.1. The steps of BG segmentation are shown in Fig. 2.Overlap in comparison with manual segmentation isquantified by the Dice similarity index (DSI) equal to0.829 (values above 0.7 are considered a very goodagreement). Such BG segmentation has been applied on 3other subjects ranging for 29 to 32 GA and the DSI hasbeen of 0.856, 0.794 and 0.785. Our segmentation of theinner (red and blue contours) and outer cortical surface(green contour) is presented in Fig. 3. Finally, torefine the results we include our WM segmentation in theFreesurfer software [12] and some manual corrections toobtain Fig.4. Discussion. Precise cortical surfaceextraction of fetal brain is needed for quantitativestudies of early human brain development. Our workcombines the well known statistical classificationframework with the active contour segmentation forcentral gray mater extraction. A main advantage of thepresented procedure for fetal brain surface extraction isthat we do not include any spatial prior coming fromanatomical atlases. The results presented here arepreliminary but promising. Our efforts are now in testingsuch approach on a wider range of gestational ages thatwe will include in the final version of this work andstudying as well its generalization to different scannersand different type of MRI sequences. References. [1]Guibaud, Prenatal Diagnosis 29(4) (2009). [2] Rousseau,Acad. Rad. 13(9), 2006, [3] Jiang, IEEE TMI 2007. [4]Warfield IADB, MICCAI 2009. [5] Claude, IEEE Trans. Bio.Eng. 51(4) (2004). [6] Habas, MICCAI (Pt. 1) 2008. [7]Bertelsen, ISMRM 2009 [8] Bach Cuadra, IADB, MICCAI 2009.[9] Styner, IEEE TMI 19(39 (2000). [10] Lee, IEEE Trans.Visual. And Comp. Graph. 3(3), 1997, [11] Chan, IEEETrans. Img. Proc, 10(2), 2001 [12] Freesurfer,http://surfer.nmr.mgh.harvard.edu.

Neurally adjusted ventilatory assist (NAVA) improves patient-ventilator interaction during non-invasive ventilation delivered by face mask.

PURPOSE: To determine if, compared to pressure support (PS), neurally adjusted ventilatory assist (NAVA) reduces patient-ventilator asynchrony in intensive care patients undergoing noninvasive ventilation with an oronasal face mask. METHODS: In this prospective interventional study we compared patient-ventilator synchrony between PS (with ventilator settings determined by the clinician) and NAVA (with the level set so as to obtain the same maximal airway pressure as in PS). Two 20-min recordings of airway pressure, flow and electrical activity of the diaphragm during PS and NAVA were acquired in a randomized order. Trigger delay (T(d)), the patient's neural inspiratory time (T(in)), ventilator pressurization duration (T(iv)), inspiratory time in excess (T(iex)), number of asynchrony events per minute and asynchrony index (AI) were determined. RESULTS: The study included 13 patients, six with COPD, and two with mixed pulmonary disease. T(d) was reduced with NAVA: median 35 ms (IQR 31-53 ms) versus 181 ms (122-208 ms); p = 0.0002. NAVA reduced both premature and delayed cyclings in the majority of patients, but not the median T(iex) value. The total number of asynchrony events tended to be reduced with NAVA: 1.0 events/min (0.5-3.1 events/min) versus 4.4 events/min (0.9-12.1 events/min); p = 0.08. AI was lower with NAVA: 4.9 % (2.5-10.5 %) versus 15.8 % (5.5-49.6 %); p = 0.03. During NAVA, there were no ineffective efforts, or late or premature cyclings. PaO(2) and PaCO(2) were not different between ventilatory modes. CONCLUSION: Compared to PS, NAVA improved patient ventilator synchrony during noninvasive ventilation by reducing T(d) and AI. Moreover, with NAVA, ineffective efforts, and late and premature cyclings were absent.

Harmonic active contours.

We propose a segmentation method based on the geometric representation of images as 2-D manifolds embedded in a higher dimensional space. The segmentation is formulated as a minimization problem, where the contours are described by a level set function and the objective functional corresponds to the surface of the image manifold. In this geometric framework, both data-fidelity and regularity terms of the segmentation are represented by a single functional that intrinsically aligns the gradients of the level set function with the gradients of the image and results in a segmentation criterion that exploits the directional information of image gradients to overcome image inhomogeneities and fragmented contours. The proposed formulation combines this robust alignment of gradients with attractive properties of previous methods developed in the same geometric framework: 1) the natural coupling of image channels proposed for anisotropic diffusion and 2) the ability of subjective surfaces to detect weak edges and close fragmented boundaries. The potential of such a geometric approach lies in the general definition of Riemannian manifolds, which naturally generalizes existing segmentation methods (the geodesic active contours, the active contours without edges, and the robust edge integrator) to higher dimensional spaces, non-flat images, and feature spaces. Our experiments show that the proposed technique improves the segmentation of multi-channel images, images subject to inhomogeneities, and images characterized by geometric structures like ridges or valleys.

Geometry-based demosaicking

Demosaicking is a particular case of interpolation problems where, from a scalar image in which each pixel has either the red, the green or the blue component, we want to interpolate the full-color image. State-of-the-art demosaicking algorithms perform interpolation along edges, but these edges are estimated locally. We propose a level-set-based geometric method to estimate image edges, inspired by the image in-painting literature. This method has a time complexity of O(S) , where S is the number of pixels in the image, and compares favorably with the state-of-the-art algorithms both visually and in most relevant image quality measures.

A post-processing pipeline to reconstruct the developing fetal brain using low-resolution MRI

Introduction. Development of the fetal brain surfacewith concomitant gyrification is one of the majormaturational processes of the human brain. Firstdelineated by postmortem studies or by ultrasound, MRIhas recently become a powerful tool for studying in vivothe structural correlates of brain maturation. However,the quantitative measurement of fetal brain developmentis a major challenge because of the movement of the fetusinside the amniotic cavity, the poor spatial resolution,the partial volume effect and the changing appearance ofthe developing brain. Today extensive efforts are made todeal with the âeurooepost-acquisitionâeuro reconstruction ofhigh-resolution 3D fetal volumes based on severalacquisitions with lower resolution (Rousseau, F., 2006;Jiang, S., 2007). We here propose a framework devoted tothe segmentation of the basal ganglia, the gray-whitetissue segmentation, and in turn the 3D corticalreconstruction of the fetal brain. Method. Prenatal MRimaging was performed with a 1-T system (GE MedicalSystems, Milwaukee) using single shot fast spin echo(ssFSE) sequences in fetuses aged from 29 to 32gestational weeks (slice thickness 5.4mm, in planespatial resolution 1.09mm). For each fetus, 6 axialvolumes shifted by 1 mm were acquired (about 1 min pervolume). First, each volume is manually segmented toextract fetal brain from surrounding fetal and maternaltissues. Inhomogeneity intensity correction and linearintensity normalization are then performed. A highspatial resolution image of isotropic voxel size of 1.09mm is created for each fetus as previously published byothers (Rousseau, F., 2006). B-splines are used for thescattered data interpolation (Lee, 1997). Then, basalganglia segmentation is performed on this superreconstructed volume using active contour framework witha Level Set implementation (Bach Cuadra, M., 2010). Oncebasal ganglia are removed from the image, brain tissuesegmentation is performed (Bach Cuadra, M., 2009). Theresulting white matter image is then binarized andfurther given as an input in the Freesurfer software(http://surfer.nmr.mgh.harvard.edu/) to provide accuratethree-dimensional reconstructions of the fetal brain.Results. High-resolution images of the cerebral fetalbrain, as obtained from the low-resolution acquired MRI,are presented for 4 subjects of age ranging from 29 to 32GA. An example is depicted in Figure 1. Accuracy in theautomated basal ganglia segmentation is compared withmanual segmentation using measurement of Dice similarity(DSI), with values above 0.7 considering to be a verygood agreement. In our sample we observed DSI valuesbetween 0.785 and 0.856. We further show the results ofgray-white matter segmentation overlaid on thehigh-resolution gray-scale images. The results arevisually checked for accuracy using the same principlesas commonly accepted in adult neuroimaging. Preliminary3D cortical reconstructions of the fetal brain are shownin Figure 2. Conclusion. We hereby present a completepipeline for the automated extraction of accuratethree-dimensional cortical surface of the fetal brain.These results are preliminary but promising, with theultimate goal to provide âeurooemovieâeuro of the normal gyraldevelopment. In turn, a precise knowledge of the normalfetal brain development will allow the quantification ofsubtle and early but clinically relevant deviations.Moreover, a precise understanding of the gyraldevelopment process may help to build hypotheses tounderstand the pathogenesis of several neurodevelopmentalconditions in which gyrification have been shown to bealtered (e.g. schizophrenia, autismâeuro¦). References.Rousseau, F. (2006), 'Registration-Based Approach forReconstruction of High-Resolution In Utero Fetal MR Brainimages', IEEE Transactions on Medical Imaging, vol. 13,no. 9, pp. 1072-1081. Jiang, S. (2007), 'MRI of MovingSubjects Using Multislice Snapshot Images With VolumeReconstruction (SVR): Application to Fetal, Neonatal, andAdult Brain Studies', IEEE Transactions on MedicalImaging, vol. 26, no. 7, pp. 967-980. Lee, S. (1997),'Scattered data interpolation with multilevel B-splines',IEEE Transactions on Visualization and Computer Graphics,vol. 3, no. 3, pp. 228-244. Bach Cuadra, M. (2010),'Central and Cortical Gray Mater Segmentation of MagneticResonance Images of the Fetal Brain', ISMRM Conference.Bach Cuadra, M. (2009), 'Brain tissue segmentation offetal MR images', MICCAI.

Vessel centerline tracking and boundary sgementation in coronary MRA with minimal manual interaction.

Magnetic resonance angiography (MRA) provides a noninvasive means to detect the presence, location and severity of atherosclerosis throughout the vascular system. In such studies, and especially those in the coronary arteries, the vessel luminal area is typically measured at multiple cross-sectional locations along the course of the artery. The advent of fast volumetric imaging techniques covering proximal to mid segments of coronary arteries necessitates automatic analysis tools requiring minimal manual interactions to robustly measure cross-sectional area along the three-dimensional track of the arteries in under-sampled and non-isotropic datasets. In this work, we present a modular approach based on level set methods to track the vessel centerline, segment the vessel boundaries, and measure transversal area using two user-selected endpoints in each coronary of interest. Arterial area and vessel length are measured using our method and compared to the standard Soap-Bubble reformatting and analysis tool in in-vivo non-contrast enhanced coronary MRA images.

Morphoquantitative description of bovine digital cushion

Abstract The digital cushion is characterized as a modified subcutaneous tissue that absorbs the shock during gait, assists venous return of the hoof and supports a considerable part of body weight. Digital cushions have particular importance in the pathogenesis of the hoof, since they need to properly work in order to prevent compression and traumas in soft tissues. This study aimed to measure and determine how is the arrangement of these structures, and for this it was established the proportions of connective, adipose, vascular tissues and collagen fibers and collagen types found in palmar and plantar digital cushion of bovine using fore and hindlimbs of twelve adult zebu cattle of both sexes, 11 male and one female, with 269kg average carcass weight and without limb disorders. Fragments of cushions were subjected to conventional histology, cut to a thickness of 4µm and stained with Red Picrosirius. With digital optical microscope, the quantification of the connective tissue and differentiation of types of collagen used the Image Pro Plus® software, and of adipose and vascular tissue, the test point system. The mean and standard error were estimated with the GraphPad Prism 5.0 software, and then data were subjected to Kolmogorov-Smirnov normality test and Student's t-test with significance level set at 5% for determining the amount of different tissues between fore and hindlimbs of studied animals. In forelimbs the mean and standard error of the connective tissue proportion was 50.10%+1.54, of the adipose tissue was 21.34%+1.44, and of vascular tissue was 3.43%+0.28. Hindlimbs presented a proportion of connective tissue of 61.61%+1.47, 20.66%+1.53 of adipose tissue, and 3.06%+0.20 of vascular tissue. A significant difference (p<0.001) was detected in the connective tissue proportion between fore and hindlimbs. Types I and II collagen fibers have presented, respectively, a proportion of 31.89% and 3.9% in forelimbs and 34.05% and 1.78% in hindlimbs. According to the used methodology, digital cushions had a clear differentiation relative to adipose tissue between fore and hindlimbs.

Simulation numérique de feux de forêt avec réinitialisation et contournement d’obstacles

Ce travail présente une technique de simulation de feux de forêt qui utilise la méthode Level-Set. On utilise une équation aux dérivées partielles pour déformer une surface sur laquelle est imbriqué notre front de flamme. Les bases mathématiques de la méthode Level-set sont présentées. On explique ensuite une méthode de réinitialisation permettant de traiter de manière robuste des données réelles et de diminuer le temps de calcul. On étudie ensuite l’effet de la présence d’obstacles dans le domaine de propagation du feu. Finalement, la question de la recherche du point d’ignition d’un incendie est abordée.