Unmixing hyperspectral data: independent and dependent component analysis

The development of high spatial resolution airborne and spaceborne sensors has improved the capability of ground-based data collection in the fields of agriculture, geography, geology, mineral identification, detection [2, 3], and classification [4–8]. The signal read by the sensor from a given spatial element of resolution and at a given spectral band is a mixing of components originated by the constituent substances, termed endmembers, located at that element of resolution. This chapter addresses hyperspectral unmixing, which is the decomposition of the pixel spectra into a collection of constituent spectra, or spectral signatures, and their corresponding fractional abundances indicating the proportion of each endmember present in the pixel [9, 10]. Depending on the mixing scales at each pixel, the observed mixture is either linear or nonlinear [11, 12]. The linear mixing model holds when the mixing scale is macroscopic [13]. The nonlinear model holds when the mixing scale is microscopic (i.e., intimate mixtures) [14, 15]. The linear model assumes negligible interaction among distinct endmembers [16, 17]. The nonlinear model assumes that incident solar radiation is scattered by the scene through multiple bounces involving several endmembers [18]. Under the linear mixing model and assuming that the number of endmembers and their spectral signatures are known, hyperspectral unmixing is a linear problem, which can be addressed, for example, under the maximum likelihood setup [19], the constrained least-squares approach [20], the spectral signature matching [21], the spectral angle mapper [22], and the subspace projection methods [20, 23, 24]. Orthogonal subspace projection [23] reduces the data dimensionality, suppresses undesired spectral signatures, and detects the presence of a spectral signature of interest. The basic concept is to project each pixel onto a subspace that is orthogonal to the undesired signatures. As shown in Settle [19], the orthogonal subspace projection technique is equivalent to the maximum likelihood estimator. This projection technique was extended by three unconstrained least-squares approaches [24] (signature space orthogonal projection, oblique subspace projection, target signature space orthogonal projection). Other works using maximum a posteriori probability (MAP) framework [25] and projection pursuit [26, 27] have also been applied to hyperspectral data. In most cases the number of endmembers and their signatures are not known. Independent component analysis (ICA) is an unsupervised source separation process that has been applied with success to blind source separation, to feature extraction, and to unsupervised recognition [28, 29]. ICA consists in finding a linear decomposition of observed data yielding statistically independent components. Given that hyperspectral data are, in given circumstances, linear mixtures, ICA comes to mind as a possible tool to unmix this class of data. In fact, the application of ICA to hyperspectral data has been proposed in reference 30, where endmember signatures are treated as sources and the mixing matrix is composed by the abundance fractions, and in references 9, 25, and 31–38, where sources are the abundance fractions of each endmember. In the first approach, we face two problems: (1) The number of samples are limited to the number of channels and (2) the process of pixel selection, playing the role of mixed sources, is not straightforward. In the second approach, ICA is based on the assumption of mutually independent sources, which is not the case of hyperspectral data, since the sum of the abundance fractions is constant, implying dependence among abundances. This dependence compromises ICA applicability to hyperspectral images. In addition, hyperspectral data are immersed in noise, which degrades the ICA performance. IFA [39] was introduced as a method for recovering independent hidden sources from their observed noisy mixtures. IFA implements two steps. First, source densities and noise covariance are estimated from the observed data by maximum likelihood. Second, sources are reconstructed by an optimal nonlinear estimator. Although IFA is a well-suited technique to unmix independent sources under noisy observations, the dependence among abundance fractions in hyperspectral imagery compromises, as in the ICA case, the IFA performance. Considering the linear mixing model, hyperspectral observations are in a simplex whose vertices correspond to the endmembers. Several approaches [40–43] have exploited this geometric feature of hyperspectral mixtures [42]. Minimum volume transform (MVT) algorithm [43] determines the simplex of minimum volume containing the data. The MVT-type approaches are complex from the computational point of view. Usually, these algorithms first find the convex hull defined by the observed data and then fit a minimum volume simplex to it. Aiming at a lower computational complexity, some algorithms such as the vertex component analysis (VCA) [44], the pixel purity index (PPI) [42], and the N-FINDR [45] still find the minimum volume simplex containing the data cloud, but they assume the presence in the data of at least one pure pixel of each endmember. This is a strong requisite that may not hold in some data sets. In any case, these algorithms find the set of most pure pixels in the data. Hyperspectral sensors collects spatial images over many narrow contiguous bands, yielding large amounts of data. For this reason, very often, the processing of hyperspectral data, included unmixing, is preceded by a dimensionality reduction step to reduce computational complexity and to improve the signal-to-noise ratio (SNR). Principal component analysis (PCA) [46], maximum noise fraction (MNF) [47], and singular value decomposition (SVD) [48] are three well-known projection techniques widely used in remote sensing in general and in unmixing in particular. The newly introduced method [49] exploits the structure of hyperspectral mixtures, namely the fact that spectral vectors are nonnegative. The computational complexity associated with these techniques is an obstacle to real-time implementations. To overcome this problem, band selection [50] and non-statistical [51] algorithms have been introduced. This chapter addresses hyperspectral data source dependence and its impact on ICA and IFA performances. The study consider simulated and real data and is based on mutual information minimization. Hyperspectral observations are described by a generative model. This model takes into account the degradation mechanisms normally found in hyperspectral applications—namely, signature variability [52–54], abundance constraints, topography modulation, and system noise. The computation of mutual information is based on fitting mixtures of Gaussians (MOG) to data. The MOG parameters (number of components, means, covariances, and weights) are inferred using the minimum description length (MDL) based algorithm [55]. We study the behavior of the mutual information as a function of the unmixing matrix. The conclusion is that the unmixing matrix minimizing the mutual information might be very far from the true one. Nevertheless, some abundance fractions might be well separated, mainly in the presence of strong signature variability, a large number of endmembers, and high SNR. We end this chapter by sketching a new methodology to blindly unmix hyperspectral data, where abundance fractions are modeled as a mixture of Dirichlet sources. This model enforces positivity and constant sum sources (full additivity) constraints. The mixing matrix is inferred by an expectation-maximization (EM)-type algorithm. This approach is in the vein of references 39 and 56, replacing independent sources represented by MOG with mixture of Dirichlet sources. Compared with the geometric-based approaches, the advantage of this model is that there is no need to have pure pixels in the observations. The chapter is organized as follows. Section 6.2 presents a spectral radiance model and formulates the spectral unmixing as a linear problem accounting for abundance constraints, signature variability, topography modulation, and system noise. Section 6.3 presents a brief resume of ICA and IFA algorithms. Section 6.4 illustrates the performance of IFA and of some well-known ICA algorithms with experimental data. Section 6.5 studies the ICA and IFA limitations in unmixing hyperspectral data. Section 6.6 presents results of ICA based on real data. Section 6.7 describes the new blind unmixing scheme and some illustrative examples. Section 6.8 concludes with some remarks.

Alterações da substância branca cerebral relacionadas com o envelhecimento

INTRODUÇÃO: Estudos prévios, com técnicas de imagem, documentam de forma consistente a existência de alterações da substância branca cerebral relacionadas com o envelhecimento (ASBRE). Tais alterações poderão ter um papel importante no declínio funcional do idoso, reflectindo‐se sobretudo no desempenho motor e cognitivo, com repercussão evidente na prática clínica. Apesar disso, a caracterização em definitivo dos fenótipos clínicos e da evolução das ASBRE continua por esclarecer, possivelmente pelas dificuldades metodológicas de que se reveste o seu estudo, incluindo: a adequação das baterias neuropsicológicas, a utilização de amostras de doentes com diferentes graus de severidade e de envolvimento regional, as limitações das diferentes escalas e a sensibilidade dos diferentes métodos de imagem. A Ressonância Magnética (RM) de difusão tem revelado grande sensibilidade para as alterações isquémicas, admitindo‐se que poderá permitir uma melhor caracterização das ASBRE e deste modo possibilitar uma correlação mais precisa com as variáveis cognitivas e motoras, permitindo avaliar ainda a substância branca aparentemente normal (SBAN). OBJECTIVOS: Descrever a evolução imagiológica das ASBRE no intervalo de um ano e analisar a sua expressão clínica e impacto funcional; identificar factores preditivos de progressão das ASBRE e de declínio funcional associado. Descrever a expressão clínica e perfil evolutivo dos doentes com ASBRE com envolvimento preferencial da região parieto‐occipital; comparar este grupo de doentes com os doentes com ASBRE, sem envolvimento preferencial desta região. Medir os coeficientes de difusão aparente (CDA), utilizando regiões de interesse (RDI), em diferentes localizações da substância branca, incluindo substância branca lesada e SBAN, descrever sua evolução temporal no intervalo de um ano e determinar suas correlações clínicas e imagiológicas. MÉTODOS: Utilizando uma amostra de conveniência, foram estudados 30 doentes, com mais de 65 anos, sem incapacidade funcional ou com incapacidade mínima, avaliada pela escala de actividades instrumentais da vida diária (IADL), apresentando ASBRE em TC. Foi utilizado um protocolo exaustivo de avaliação clínica (com particular destaque para as funções motoras e cognitivas) e imagiológica, em dois momentos de avaliação separados por um ano de intervalo (t0 e t1). As ASBRE foram avaliadas com escalas visuais, escala ARWMC e escala de Fazekas, e os doentes foram estudados em função do grau de severidade (ligeiro versus moderado a grave na escala de Fazekas) e de um envolvimento preferencial posterior (definido como 2 ou mais pontos na escala ARWMC na região parieto‐occipital por comparação com a região frontal). Os CDA foram avaliados mediante estudo de RDI, na substância branca frontal lesada (SBFL) e SBAN frontal, parieto‐occipital e dos pedúnculos cerebelosos. Para verificar diferenças na ordem de distribuição das variáveis foi usado o teste de Mann‐Whitney e para comparação de proporções, o teste exacto de Fisher. Na comparação entre a avaliação em t0 e t1 foi usado o teste Wilcoxon Signed Ranks na comparação da distribuição da ordem das variáveis e o teste McNemar na análise de frequências. Na análise correlacional foram utilizados os testes de T para variáveis emparelhadas e as correlações entre estas foram efectuadas com o coeficiente de correlação de Spearman ou de Pearson. O trabalho foi aprovado pela Comissão de Ética do hospital onde foi realizado e todos os doentes incluídos assinaram um consentimento informado. RESULTADOS: A idade média da população estudada foi 72,5 anos (17 doentes eram do sexo masculino). No final de um ano, 1 doente tinha falecido e 3 doentes não completaram a avaliação imagiológica. Registou‐se uma progressão significativa das ASBRE segundo a escala ARWMC (t0: 8,37 / t1: 9,65 ; p<0,001). Na análise funcional, motora e cognitiva, não houve um agravamento significativo. Avaliando os doentes em t0 e t1 segundo o grau de severidade das ASBRE, o grupo com atingimento moderado a grave (ASBRE2) comparado com o grupo com atingimento ligeiro (ASBRE1) apresentava: maior extensão de lesão da substância branca (ARWMC t0: 11,9 / 4,8 ; p<0.001 ; t1: 14,0 / 5,9 ; p<0,001); tendência a pior desempenho funcional (IADL t0: 90,7 / 99,2 ; p=0,023; t1: 86,4 / 96,7 ; p=n.s.) e motor (SPPB t0: 9,8 / 10,3 ; p=n.s. ; t1: 9,5 / 10,5 ; p=0,058); tendência a maior compromisso do humor (Escala Cornell t0: 6,7 / 3,5 ; p=0,037; t1: 6,2 / 4,5 ; p=n.s.). Analisando a evolução, de t0 para t1, de cada um dos grupos (ASBRE2 e ASBRE1) registou‐se: aumento da extensão da lesão da substância branca em ambos (ASBRE2: 12,0 / 14,0;z=‐2,687 ; p=0,007; ASBR1: 4,8 / 5,9 ; z=‐2,724 ; p=0,006); variação não significativa funcional e motora; tendência ao agravamento em ambos na prova de Cancelamento de dígitos (ASBRE2: 17,5 / 17,4 ; p=n.s. ; ASBRE1: 19,9 / 16,9 ; z=‐2,096 ; p=0,036);tendência à melhoria em ambos no MMS (ASBRE2: 25,7 / 27,5 ; z=‐2,155 ; p=0,031; ASBRE1: 27,5 / 28,2 ; p=n.s). Avaliando os doentes em t0 e t1 em função do padrão de distribuição das ASBRE, os doentes com um envolvimento preferencial posterior (ASBREP) comparados com os restantes (ASBREnP), apresentavam: maior extensão da lesão (ARWMC t0: 10,8 / 6,9 ; p=0,025; t1: 12,9 / 7,6 ; p=0,011); diferenças não significativas no desempenho motor; tendência a melhor desempenho na prova dos Labirintos (t0: 8,1 / 11,8 ; p=0,06; t1: 8,7 / 9,5 ; p=n.s.) e Cancelamento de dígitos (t0: 20,9 / 17,4 ; p=0,045; t1: 18,5 / 16,3 ; p=n.s.); tendência a maior compromisso depressivo na GDS (t0: 5,0 / 3,68 ; p=n.s. ; t1: 5,7 / 3,3 p=0,033). Analisando o perfil evolutivo de t0 para t1, registou‐se: aumento da extensão da lesão nos dois grupos (ASBREP: 10,8 / 12,9 ; z=‐2,555 ; P=0,011; ASBREnP: 6,4 / 7,6 ; z=‐2,877 ; p=0,04); variação em sentidos diferentes com melhoria funcional no grupo ASBREP (91,0 / 95,5 ; z=‐0,926 ; p=0,036) e agravamento no grupo ASBREnP (96,7 / 89,8 ; z=‐2,032 ; p=0,042); variação sem sentidos diferentes, com agravamento significativo no grupo ASBREnP no item estação de pé do SPPB (ASBREP 3,8/3,9 p=n.s.; ASBREnP 3,9/3,6; z=‐2,236 ; p=0,025); tendência à melhoria nos dois grupos no MMS (ASBREP: 27,2 / 28,2 ; p=n.s.; ASBREnP: 26,3 / 27,7 ; z=‐2,413 ; p=0,016) e tendência em sentidos diferentes no Trail Making, com eventual melhoria no grupo ASBREP (113,9 / 91,6 ; p=n.s.) e agravamento no grupo ASBREnP (113,7 / 152,0 ; z=‐2,155 ; p=0,031). Na análise da imagem, utilizando a escala ARWMC e o estudo dos CDA, na avaliação transversal na inclusão, a comparação entre as pontuações médias da escala ARWML nas diferentes regiões mostrava diferenças significativas (F=39,54 , p<0,0001). A análise comparativa post‐hoc de Bonferroni mostrou valores significativamente mais altos para as regiões frontais e parieto‐occipitais (p<0,0001). Os valores médios dos CDA eram significativamente diferentes entre regiões (F=44,56; p<0,0001), sendo mais altos na SBFL (p<0,0001). Não existia diferença significativa entre os valores registados na SBAN nas regiões frontais e parieto‐occipitais. As pontuações regionais da escala ARWMC e os valores médios dos CDA correlacionavam‐se todos de forma positiva. A pontuação da escala ARWMC na região frontal correlacionava‐se significativamente com os valores do CDA da SBFL (r=0,467 ; p=0,012). Existia tendência para uma correlação positiva entre as pontuações da escala ARWMC na região frontal e os valores médios dos CDA na SBAN frontal (r=0,276 ; p=0,155). As pontuações da escala ARWMC e os CDA correlacionavam‐se de forma positiva com a idade e com a tensão arterial (TA). Foram encontradas correlações significativas entre: idade e SBAN frontal (r=0,440 ; p=0,019); TA diastólica e SBFL (r=0,386 ; p=0,034); TA sistólica e SBAN Parieto‐occipital (r=0,407 ; P=0,032). Na avaliação motora e cognitiva, dado elevado número de variáveis, foi efectuada uma análise de factor principal. Registou‐se uma tendência global negativa na correlação entre as pontuações da escala visual na região frontal, os valores dos CDA, e o desempenho motor e cognitivo. Na análise evolutiva, (n=19), registou‐se variação significativa dos CDA, com aumento na SBFL (Direita: z=‐2,875 ; p=0,004 ; Esquerda: z=‐2,113 ; p=0,035) e diminuição na SBAN dos pedúnculos cerebelosos (Direita: z=‐2,094 ; p=0,036 ; Esquerda: z=‐1,989 ; p=0,047). Foi observada uma correlação negativa entre a variação do CDA na SBAN dos pedúnculos cerebelosos e na SBFL contralateral (SBAN pedúnculo cerebeloso Esquerdo / SBFL Direita: r=‐0,133 ; p=n.s.; SBAN pedúnculo cerebeloso Direito / SBFL Esquerda: r=‐0,561 ; p=0,012). Os valores dos CDA à direita correlacionavam‐se de forma positiva com a velocidade da marcha (r=0,562 ; p=0,012). CONCLUSÕES: A progressão das ASBRE pode ser observada com uma escala visual detalhada no intervalo de um ano. Contudo, o eventual agravamento da incapacidade funcional, motora e cognitiva, não parece ser apreciável em igual intervalo de tempo. A maior severidade das ASBRE associa‐se a uma tendência para um maior compromisso funcional, motor e possivelmente do humor. A questão da progressão em escalas simplificadas, de um estádio ligeiro para um estádio moderado a grave, não é elucidada pelos resultados do presente trabalho. Os doentes com um envolvimento preferencial da região parieto‐occipital poderão constituir um subgrupo distinto que, apesar de ter maior extensão de lesão, parece ter um melhor desempenho motor e cognitivo. O perfil evolutivo destes doentes parece igualmente ser distinto, não se observando a tendência ao agravamento funcional, motor e cognitivo (sobretudo em provas de função executiva) que se encontra nos restantes doentes. A análise transversal na inclusão, utilizando uma escala visual e o estudo dos CDA, sugere que a severidade das ASBRE se correlaciona com o compromisso motor e cognitivo, bem como com a idade e com a TA. Uma maior vulnerabilidade da substância branca frontal à lesão vascular parece ter um papel importante no compromisso motor e na disfunção executiva, (essencialmente à custa do compromisso da atenção), possivelmente associada à desconexão dos circuitos fronto‐subcorticais. A análise dos CDA sugere que isso é válido igualmente para a SBAN e sublinha que, as imagens de RM convencional poderão não traduzir a verdadeira extensão da lesão e consequentemente do compromisso motor e cognitivo. A relação entre a progressão da doença vascular em lesões frontais constituídas e a redução do CDA no pedúnculo cerebeloso contralateral poderá estar associada a um pior desempenho motor. A disrupção dos circuitos fronto‐cerebelosos, determinando hipometabolismo e diminuição da perfusão no cerebelo, poderá ser responsável pela diminuição do CDA no cerebelo. ABSTRACT INTRODUCTION: Previous studies, with new imaging techniques, have consistently documented the presence of age‐related white matter lesions (ARWML), emphasizing their role in agerelated functional decline, mainly related to motor and cognitive impairment, and inherent consequences in clinical practice. However clinical significance of ARWML remains to be elucidated, probably on account of methodological difficulties such as: specific neuropsychological batteries, utilization of samples with different degrees of severity and regional involvement, utilization of different imaging scales and different sensitivity of imaging techniques. Recently, Diffusion Weighted Magnetic Ressonance imaging (DWI) has shown a higher sensitivity to ischemic lesions, suggesting it might be superior for characterization of ARWML, allowing more precise correlation with motor and cognitive variables, and evaluating also normal appearing white matter (NAWM). OBJECTIVES: To describe imagiologic evolution of ARWML within one year interval and to analyse its clinical and functional significance. To identify predictors of ARWML progression and associated functional impairment. To describe clinical characteristics and evolution profile of patients with predominantly posterior lesions; to compare this group of patients with patients without predominantly posterior lesions. To study average Apparent Diffusion Coeficcients (ADC) in different white matter regions using regions of interest (ROI); to analyse their evolution profile and to determine their clinical and imagiologic correlations. METHODS: A sample of 30 patients older than 65 years, without functional impairment or with minimal impairment, according to the Instrumental Activities of Daily Lliving scale, with ARWML on CT scan, were studied in a cross‐sectional design. An extensive clinical(with detailed motor and cognitive evaluation) and imagiologic protocol was applied in two one‐year interval separate moments (t0 and t1). ARWML were studied using visual scales, ARWMC and Fazekas’s scale, and patients were studied according to degree of severity (Fazekas scale mild versus moderate / severe) and preferential involvement of the posterior region (defined as 2 or more points in the ARWMC scale in the parietooccipital region compared with frontal region). Evaluation of ADC was performed using ROI in frontal lesioned white matter (FLWM) and NAWM (frontal, parieto‐occipital and cerebellar regions). To study differences in the distribution of variables the Mann‐Whitney test was used and to compare proportions the exact Fisher Test was used. To compare temporal evolution profile between t0 and t1, the Wilcoxon Signed ranks Test was used to analyse the distribution of variables and the Mc Nemar Test to analyse frequencies. Correlation analysis was performed using Spearman or Pearson tests. The study was approved by the local Ethics Committee and all patients signed an informed consent. RESULTS: Mean age was 72.5 years (17 patients were male). By the end of the study, one patient was dead and 3 patients did not undergo brain imaging. There was a higher extent of ARWML evaluated with the ARWMC scale (t0: 8.37 / t1: 9.65 ; p<0.001). Functional, motor and cognitive performance did not progress significantly. Evaluating patients in t0 and t1 according to the degree of severity (Fazekas scale), the moderate / severe group of patients (WML2), compared with the mild group (WML1), showed: higher extent of lesion (ARWMC scale t0: 11.9 / 4.8 ; p<0.001 ; t1: 14.0 / 5.9 ; p<0.001); tendency to worse functional (IADL t0: 90.7 / 99.2 ; p=0.023; t1: 86.4 / 96.7 ; p=n.s.) and motor (SPPB t0: 9.8 / 10.3 ; p=n.s. ; t1: 9.5 / 10.5 ; p=0.058) performance; tendency to higher depressive scores (Cornell Scale t0: 6.7 / 3.5 ; p=0.037; t1: 6.2 / 4.5; p=n.s.). Analysing the evolution profile from t0 to t1 of each group (WML2 and WML1), there was a higher extent of lesion (ARWMC scale) in both (WML2: 12.0 / 14.0; z=‐2.687 ; p=0.007; WML1: 4.8 / 5.9 ; z=‐2.724 ; p=0.006); non‐significant variation in functional and motor performances; tendency to worse performance on the Digit Cancelling (WML2: 17.5 / 17.4 ; p=n.s. ; WML1: 19.9 / 16.9 ; z=‐2.096 ; p=0,036) and to better performance on the MMS (WML2: 25.7 / 27.5 ; z=‐2.155 ; p=0.031; WML1: 27.5/ 28.2 ; p=n.s). Evaluating patients in t0 and t1 according to the regional distribution of ARWML, patients with predominantly posterior lesions (WMLP) compared with the rest of the group (WMLnP), showed: higher extent of lesion (ARWMC scale t0: 10.8 / 6.9 ; p=0.025; t1:12.9 / 7.6 ; p=0.011); non significant differences on motor evaluation; tendency to a better performance on Maze (t0: 8.1 / 11.8 ; p=0.06; t1: 8.7 / 9.5 ; p=n.s.) and Digit cancelling (t0: 20.9 / 17.4 ; p=0.045; t1: 18.5 / 16.3 ; p=n.s.) tests;tendency to higher scores on GDS (t0: 5.0 / 3.68 ; p=n.s. ; t1: 5.7 / 3.3 p=0.033). Analysing the evolution profile from t0 to t1 of each group (WMLP and WMLnP), there was: higher extent of lesion (ARWMC scale) in both groups (WMLP: 10.8 / 12.9 ;z=‐2,555 ; P=0,011; WMLnP: 6.4 / 7.6 ; z=‐2.877; p=0.04); variation in different directions with better functional performance in the group WMLP (91.0 / 95.5 ;z=‐0.926 ; p=0.036) and worse in WMLnP (96.7 / 89.8 ; z=‐2.032 ; p=0.042); variation in different directions with worse motor performance in one SPPB item (total stands) in the group WMLnP (WMLP 3.8/3.9 p=n.s.; ASBREnP 3.9/3.6; z=‐2.236 ; p=0.025);tendency to improvement in both groups in MMS (WMLP: 27.2 / 28.2 ; p=n.s.; WMLnP:26.3 / 27.7 ; z=‐2.413 ; p=0.016); tendency to a variation in different directions in the Trail Making Test, with possible improvement in the group WMLP (113.9 / 91.6 ;p=n.s.) and worsening in the group WMLnP (113.7 / 152.0 ; z=‐2.155 ; p=0.031). Imaging analysis in the inclusion, using the ARWMC scale and ADC evaluation, showed significant differences in different regions (F=39.54, p<0.0001). Comparative post‐hoc Bonferroni analysis showed significantly higher scores in the frontal and parieto‐occipital regions (p<0.0001. ADC values were significantly different between regions (F=44.56; p<0.0001), being higher in FLWM (p<0‐0001). There was no significant difference between ADC in NAWM in frontal and parieto‐occipital regions. ARWMC scores and ADC values correlated positively. Significant correlations were found between frontal ARWMC score and FLWM ADC values (r=0.467 ; p=0.012). ARWMC scores and ADC values correlated positively with age and blood pressure. Significant correlations were: age and frontal NAWM (r=0.440 ; p=0.019); Diastolic blood pressure and FLWM (r=0.386 ; p=0.034); sistolic blood pressure and parietooccipital NAWM (r=0.407 ; P=0.032). Due to the higher number of motor and cognitive variables a preliminary study was done, using principal component analysis. A global tendency to a negative correlation was found between ARWMC scores, ADC values and motor and cognitive performances. Evolutive analysis of ADC (n=19), showed a significant variation, with higher values in t1 in FLWM (Right: z=‐2.875 ; p=0.004 ; Left: z=‐2.113 ; p=0.035) and lower values in t1 in cerebellar NAWM (Right: z=‐2.094 ; p=0.036 ; Left: z=‐1.989 ; p=0.047). A negative correlation was found between ADC variation in cerebellar NAWM and contralateral FLWM (Left cerebellar NAWM / Right FLWM: r=‐0.133 ; p=n.s.; Right cerebellar NAWM/ Left FLWM: r=‐0.561 ; p=0.012). ADC values on the right correlated positively with walking speed (r=0,562 ; p=0,012). CONCLUSIONS: Progression of ARWML can be documented with a detailed visual scale in a one year interval. However, functional, motor and cognitive impairment, do not seem to progress significantly within the same period. A higher severity of ARWML is associated with a tendency to a worse functional and motor performance (and possibly to higher scores in depression scales). The issue of progression in a simplified visual scale from a mild to a moderate / severe degree of ARWML is not further elucidated. Patients with predominantly posterior lesions may be a subset of ARWML patients, with a different profile, that despite higher extent of lesion, seem to fair better than the rest of the group, namely with better performance on motor and cognitive tests. Evolution profile of this subset of patients also seems to be different, without a clearcut tendency to worsening functional, motor and cognitive (particularly for executive function tests) performance that is observed in the rest of the group. Imaging analysis, with a visual scale and ADC evaluation, suggests that severity of ARWML correlates negatively with cognitive and motor performance and positively with age and blood pressure. A higher vulnerability of frontal white matter to vascular disease seems to play an important role in motor and cognitive dysfunction, mainly determined by impairment of attention skills associated with frontal‐subcortical disconnection. DWI results, suggest that this may also be true for NAWM, underlining that conventional MR images may not represent the true extent of cognitive decline. The relation between vascular disease progression inside frontal lesions and ADC reduction in contralateral cerebellar peduncles, may be associated with a worse motor performance. Disruption of fronto‐cerebellar cicuits, with associated regional hypometabolism, may be responsible for the reduction of cerebellar ADC.

Vertex component analysis: a geometric-based approach to unmix hyperspectral data

Hyperspectral remote sensing exploits the electromagnetic scattering patterns of the different materials at specific wavelengths [2, 3]. Hyperspectral sensors have been developed to sample the scattered portion of the electromagnetic spectrum extending from the visible region through the near-infrared and mid-infrared, in hundreds of narrow contiguous bands [4, 5]. The number and variety of potential civilian and military applications of hyperspectral remote sensing is enormous [6, 7]. Very often, the resolution cell corresponding to a single pixel in an image contains several substances (endmembers) [4]. In this situation, the scattered energy is a mixing of the endmember spectra. A challenging task underlying many hyperspectral imagery applications is then decomposing a mixed pixel into a collection of reflectance spectra, called endmember signatures, and the corresponding abundance fractions [8–10]. Depending on the mixing scales at each pixel, the observed mixture is either linear or nonlinear [11, 12]. Linear mixing model holds approximately when the mixing scale is macroscopic [13] and there is negligible interaction among distinct endmembers [3, 14]. If, however, the mixing scale is microscopic (or intimate mixtures) [15, 16] and the incident solar radiation is scattered by the scene through multiple bounces involving several endmembers [17], the linear model is no longer accurate. Linear spectral unmixing has been intensively researched in the last years [9, 10, 12, 18–21]. It considers that a mixed pixel is a linear combination of endmember signatures weighted by the correspondent abundance fractions. Under this model, and assuming that the number of substances and their reflectance spectra are known, hyperspectral unmixing is a linear problem for which many solutions have been proposed (e.g., maximum likelihood estimation [8], spectral signature matching [22], spectral angle mapper [23], subspace projection methods [24,25], and constrained least squares [26]). In most cases, the number of substances and their reflectances are not known and, then, hyperspectral unmixing falls into the class of blind source separation problems [27]. Independent component analysis (ICA) has recently been proposed as a tool to blindly unmix hyperspectral data [28–31]. ICA is based on the assumption of mutually independent sources (abundance fractions), which is not the case of hyperspectral data, since the sum of abundance fractions is constant, implying statistical dependence among them. This dependence compromises ICA applicability to hyperspectral images as shown in Refs. [21, 32]. In fact, ICA finds the endmember signatures by multiplying the spectral vectors with an unmixing matrix, which minimizes the mutual information among sources. If sources are independent, ICA provides the correct unmixing, since the minimum of the mutual information is obtained only when sources are independent. This is no longer true for dependent abundance fractions. Nevertheless, some endmembers may be approximately unmixed. These aspects are addressed in Ref. [33]. Under the linear mixing model, the observations from a scene are in a simplex whose vertices correspond to the endmembers. Several approaches [34–36] have exploited this geometric feature of hyperspectral mixtures [35]. Minimum volume transform (MVT) algorithm [36] determines the simplex of minimum volume containing the data. The method presented in Ref. [37] is also of MVT type but, by introducing the notion of bundles, it takes into account the endmember variability usually present in hyperspectral mixtures. The MVT type approaches are complex from the computational point of view. Usually, these algorithms find in the first place the convex hull defined by the observed data and then fit a minimum volume simplex to it. For example, the gift wrapping algorithm [38] computes the convex hull of n data points in a d-dimensional space with a computational complexity of O(nbd=2cþ1), where bxc is the highest integer lower or equal than x and n is the number of samples. The complexity of the method presented in Ref. [37] is even higher, since the temperature of the simulated annealing algorithm used shall follow a log( ) law [39] to assure convergence (in probability) to the desired solution. Aiming at a lower computational complexity, some algorithms such as the pixel purity index (PPI) [35] and the N-FINDR [40] still find the minimum volume simplex containing the data cloud, but they assume the presence of at least one pure pixel of each endmember in the data. This is a strong requisite that may not hold in some data sets. In any case, these algorithms find the set of most pure pixels in the data. PPI algorithm uses the minimum noise fraction (MNF) [41] as a preprocessing step to reduce dimensionality and to improve the signal-to-noise ratio (SNR). The algorithm then projects every spectral vector onto skewers (large number of random vectors) [35, 42,43]. The points corresponding to extremes, for each skewer direction, are stored. A cumulative account records the number of times each pixel (i.e., a given spectral vector) is found to be an extreme. The pixels with the highest scores are the purest ones. N-FINDR algorithm [40] is based on the fact that in p spectral dimensions, the p-volume defined by a simplex formed by the purest pixels is larger than any other volume defined by any other combination of pixels. This algorithm finds the set of pixels defining the largest volume by inflating a simplex inside the data. ORA SIS [44, 45] is a hyperspectral framework developed by the U.S. Naval Research Laboratory consisting of several algorithms organized in six modules: exemplar selector, adaptative learner, demixer, knowledge base or spectral library, and spatial postrocessor. The first step consists in flat-fielding the spectra. Next, the exemplar selection module is used to select spectral vectors that best represent the smaller convex cone containing the data. The other pixels are rejected when the spectral angle distance (SAD) is less than a given thresh old. The procedure finds the basis for a subspace of a lower dimension using a modified Gram–Schmidt orthogonalizati on. The selected vectors are then projected onto this subspace and a simplex is found by an MV T pro cess. ORA SIS is oriented to real-time target detection from uncrewed air vehicles using hyperspectral data [46]. In this chapter we develop a new algorithm to unmix linear mixtures of endmember spectra. First, the algorithm determines the number of endmembers and the signal subspace using a newly developed concept [47, 48]. Second, the algorithm extracts the most pure pixels present in the data. Unlike other methods, this algorithm is completely automatic and unsupervised. To estimate the number of endmembers and the signal subspace in hyperspectral linear mixtures, the proposed scheme begins by estimating sign al and noise correlation matrices. The latter is based on multiple regression theory. The signal subspace is then identified by selectin g the set of signal eigenvalue s that best represents the data, in the least-square sense [48,49 ], we note, however, that VCA works with projected and with unprojected data. The extraction of the end members exploits two facts: (1) the endmembers are the vertices of a simplex and (2) the affine transformation of a simplex is also a simplex. As PPI and N-FIND R algorithms, VCA also assumes the presence of pure pixels in the data. The algorithm iteratively projects data on to a direction orthogonal to the subspace spanned by the endmembers already determined. The new end member signature corresponds to the extreme of the projection. The algorithm iterates until all end members are exhausted. VCA performs much better than PPI and better than or comparable to N-FI NDR; yet it has a computational complexity between on e and two orders of magnitude lower than N-FINDR. The chapter is structure d as follows. Section 19.2 describes the fundamentals of the proposed method. Section 19.3 and Section 19.4 evaluate the proposed algorithm using simulated and real data, respectively. Section 19.5 presents some concluding remarks.

Development of Secondary Resistance to Fluconazole in Cryptococcus neoformans Isolated from a Patient with AIDS

Cryptococcus neoformans is the fifth most common opportunistic agent of infection in patients with AIDS in the USA, exceeded only by Candida species, Pneumocystis carinii, cytomegalovirus and Mycobacterium avium1, 2, 6, 10, 11. In Brazil is the sixth, exceeded by Candida species, P. carinii, Mycobacterium species, Toxoplasma gondii, and herpes simplex virus (AIDS, Boletim Epidemiológico, set/nov 96, Ministério da Saúde, Brasil). During 30 years, the treatment of C. neoformans meningitis was based on the use of amphotericin B with or without flucytosine13. Nowadays, with the immunodepression caused by human immunodeficiency virus (HIV) infection and the availability of new antifungal drugs as the triazoles, the concept related to cure and relapses of cryptococcosis has been altered7, 20. Patients are treated with amphotericin B with or without flucytosine as initial therapy, but maintenance therapy is always necessary in AIDS patients with C. neoformans infections

Robustness of corroded reinforced concrete structures. A structural performance approach

Structure and Infrastructure Engineering, 1-17

Infective endocarditis (IE) first diagnosed at autopsy: analysis of 31 cases in Ribeirão Preto, Brazil

Thirty one infective endocarditis (IE) fatal cases whose diagnosis was first obtained at autopsy were studied. The clinical data of these patients (Group 1) showed significant differences compared to other 141 IE cases (Group 2). The average age of 53 years in Group 1 patients was 18 years higher than that of Group 2. The Group 1 patients had a low frequency of IE predisposing heart disease. Both patient groups presented fever (about 87%), but a significant low frequency of cardiac murmur (25.8%) was observed in Group 1 patients and echocardiography tests were performed in only 16.1%, suggesting that IE diagnosis was not suspected. Likewise, although most Group 1 patients appeared with severe acute illness, they did not present the classic IE clinical presentation. Blood cultures were performed in only 64.5% of the Group 1 patients. However, bacteria were isolated in 70% of these blood cultures and Staphylococcus aureus was isolated in 71.4%. The bacteria attacked mitral and aortic valves. Complications such as embolizations and cardiac failure occurred in almost half of the cases and they also presented with infections of the lungs, urinary tract, and central nervous system. Medical procedures were performed in practically all fatal cases whose diagnosis was first obtained at autopsy. Sepsis occurred in about half of the patients and it was followed by shock in more than 25%. This form of IE must be suspected in mature and in old febrile hospitalized patients having infection predisposing diseases, embolization, and suffering medical procedures.

Analysis of proteins from membrane and soluble fractions of Giardia duodenalis trophozoites of two Brazilian axenic strains

In the present study, we have analyzed by sodium docecyl sulphate - polyacrilamide gel electrophoresis (SDS-PAGE), immunoblotting and Concanavalin A blotting (Con A blotting) proteins of membrane fractions and soluble fractions obtained from Giardia duodenalis trophozoites of two axenic strains isolated in Brazil from a symptomatic (BTU-11) and an asymptomatic patient (BTU-10), as compared to the reference strain Portland 1. Both Brazilian strains showed a complex and homogeneous electrophoretic pattern of proteins, but some differences could be observed. Several glycoproteins were detected, particularly the proteins of 81, 72, 59 kDa and the protein of 62 kDa in the membrane proteins and cytosol, respectively. Many antigenic components were revealed by anti-Giardia rabbit IgG antibodies in the immunoblotting analysis. Among these components, the membrane protein of 32 kDa and the cytosol protein of 30 kDa could be related to giardin, as previously demonstrated.

Analysis of the clonal relationship among clinical isolates of Salmonella enterica serovar Infantis by different typing methods

Salmonella Infantis has been the second most common serovar in Argentina in the last two years, being isolated mostly from paediatric hospitalised patients. In order to determine the clonal relationship among Salmonella Infantis strains, we examined 15 isolates from paediatric patient faeces in Argentina (12 geographically related and 3 geographically non-related) by using antimicrobial susceptibility, plasmid profiling, repetitive extragenic palindromic (REP) PCR, enterobacterial repetitive intergenic consensus (ERIC) PCR, and low-frequency restriction analysis of chromosomal DNA by pulsed field gel electrophoresis (PFGE). Four Spanish strains were included as controls of clonal diversity in molecular techniques. Antibiotype and plasmid profile was not useful as epidemiological tools. PFGE and REP-PCR were able to discriminate between Argentinean and Spanish isolates of Salmonella Infantis allowing to detect genetically related strains in three different cities. This finding indicates that a possible spread of a clone of this serovar in the North-eastern Region of Argentina has taken place in 1998.

Novel LTCC-potentiometric microfluidic device for biparametric analysis of organic compounds carrying plastic antibodies as ionophores: Application to sulfamethoxazole and trimethoprim

Monitoring organic environmental contaminants is of crucial importance to ensure public health. This requires simple, portable and robust devices to carry out on-site analysis. For this purpose, a low-temperature co-fired ceramics (LTCC) microfluidic potentiometric device (LTCC/μPOT) was developed for the first time for an organic compound: sulfamethoxazole (SMX). Sensory materials relied on newly designed plastic antibodies. Sol–gel, self-assembling monolayer and molecular-imprinting techniques were merged for this purpose. Silica beads were amine-modified and linked to SMX via glutaraldehyde modification. Condensation polymerization was conducted around SMX to fill the vacant spaces. SMX was removed after, leaving behind imprinted sites of complementary shape. The obtained particles were used as ionophores in plasticized PVC membranes. The most suitable membrane composition was selected in steady-state assays. Its suitability to flow analysis was verified in flow-injection studies with regular tubular electrodes. The LTCC/μPOT device integrated a bidimensional mixer, an embedded reference electrode based on Ag/AgCl and an Ag-based contact screen-printed under a micromachined cavity of 600 μm depth. The sensing membranes were deposited over this contact and acted as indicating electrodes. Under optimum conditions, the SMX sensor displayed slopes of about −58.7 mV/decade in a range from 12.7 to 250 μg/mL, providing a detection limit of 3.85 μg/mL and a sampling throughput of 36 samples/h with a reagent consumption of 3.3 mL per sample. The system was adjusted later to multiple analyte detection by including a second potentiometric cell on the LTCC/μPOT device. No additional reference electrode was required. This concept was applied to Trimethoprim (TMP), always administered concomitantly with sulphonamide drugs, and tested in fish-farming waters. The biparametric microanalyzer displayed Nernstian behaviour, with average slopes −54.7 (SMX) and +57.8 (TMP) mV/decade. To demonstrate the microanalyzer capabilities for real applications, it was successfully applied to single and simultaneous determination of SMX and TMP in aquaculture waters.

Dynamical Analysis and Visualization of Tornadoes Time Series

In this paper we analyze the behavior of tornado time-series in the U.S. from the perspective of dynamical systems. A tornado is a violently rotating column of air extending from a cumulonimbus cloud down to the ground. Such phenomena reveal features that are well described by power law functions and unveil characteristics found in systems with long range memory effects. Tornado time series are viewed as the output of a complex system and are interpreted as a manifestation of its dynamics. Tornadoes are modeled as sequences of Dirac impulses with amplitude proportional to the events size. First, a collection of time series involving 64 years is analyzed in the frequency domain by means of the Fourier transform. The amplitude spectra are approximated by power law functions and their parameters are read as an underlying signature of the system dynamics. Second, it is adopted the concept of circular time and the collective behavior of tornadoes analyzed. Clustering techniques are then adopted to identify and visualize the emerging patterns.

Cutaneous anthrax in Lima, Peru: retrospective analysis of 71 cases, including four with a meningoencephalic complication

Anthrax is a zoonosis produced by Bacillus anthracis, and as an human infection is endemic in several areas in the world, including Peru. More than 95% of the reported naturally acquired infections are cutaneous, and approximately 5% of them can progress to meningoencephalitis. In this study we review the clinical and epidemiological characteristics of the patients with diagnosis of cutaneous anthrax evaluated between 1969 and 2002 at the Hospital Nacional Cayetano Heredia (HNCH) and the Instituto de Medicina Tropical Alexander von Humboldt in Lima, Peru. Seventy one patients were included [49/71 (69%) of them men], with a mean age of 37 years. The diagnoses were classified as definitive (44%) or probable (56%). The most common occupation of the patients was agriculture (39%). The source of infection was found in 63 (88.7%) patients. All the patients had ulcerative lesions, with a central necrosis. Most of the patients (65%) had several lesions, mainly located in the upper limbs (80%). Four patients (5.6%) developed meningoencephalitis, and three of them eventually died. In conclusion, considering its clinical and epidemiological characteristics, cutaneous anthrax must be included in the differential diagnosis of skin ulcers. A patient with clinical suspicion of the disease should receive effective treatment soon, in order to avoid neurological complications which carry a high fatality rate.

Clinical and epidemiological analysis of patients with HIV/AIDS admitted to a reference hospital in the northeast region of Brazil

The AIDS epidemic has become a worldwide phenomenon of enormous magnitude and extension, deeply transforming medical practices and public health initiatives. This retrospective survey aimed to analyze clinical and epidemiological characteristics of patients with HIV/AIDS admitted to the Institute of Tropical Diseases Natan Portella, Teresina, Piauí, Brazil, from January, 2001 through December, 2004. Of the 828 patients, 43% were from other states and 71.3% were men. Average patient age was 35.4 ± 11.5 years-old and 85.5% were illiterate or had primary education. The main form of exposure to HIV was heterosexual behavior (54.1%), while injectable drug use was confirmed by only 2.7% of registered cases. The most frequent infectious complications were candidiasis (42.4%) and pneumocystosis (22.2%). Sixty-eight cases (8.2%) of visceral leishmaniasis were registered. Using multivariate analysis, individuals aged over 40 years-old, patients with active tuberculosis, Pneumocystis carinii pneumonia and central nervous system cryptococcosis showed increased risk of death. In this study, young male adults with low educational levels predominated and the most frequent opportunistic infections were candidiasis and pneumocystosis.

Schistosomiasis haematobia: histopathological course determined by cystoscopy in a patient in whom praziquantel treatment failed

Schistosomiasis haematobia or urinary schistosomiasis is one of the main public health problems in Africa and the Middle East. A single dose of 40 mg praziquantel per kg body weight continues to be the treatment of choice for this infection. The aims of this follow-up were to study the post-treatment course of a patient infected with S. haematobium and not submitted to re-exposure, and to identify complications of the disease and/or therapeutic failure after praziquantel treatment by histopathological analysis. Treatments were repeated under medical supervision to ensure the correct use of the drug. In view of the suspicion of lesions in cystoscopy, the patient was submitted to bladder biopsy. The histopathological characteristics observed in biopsies obtained, after each treatment, indicated viability of parasite eggs and activity of granulomas.

Municipalities of higher vulnerability to Sylvatic Yellow Fever occurrence in the São Paulo State, Brazil

Until 1999 the endemic cases of Sylvatic Yellow Fever were located in the states of northern, midwestern and pre-Amazon regions. Since then, the disease progressively expanded its territory of occurrence, cases being registered beyond the traditional boundaries of endemism. The São Paulo State is considered to be part of this context, since after decades without registration of autochthonous cases of the disease, it reported, in 2000 and 2008-2009, epizootic occurrence in non-human primates and 30 cases in humans. Facts like these, added to the increase in incidences of serious adverse effects resulting from the Yellow Fever vaccination, have highlighted the importance of defining priority municipalities for vaccination against the disease in the state. Two groups of municipalities, some affected and some non-affected by YF, were compared for environmental variables related to the eco-epidemiology of the disease according to literature. The Multiple Correspondence Analysis (MCA) was used to pinpoint the factor able to differentiate the two groups of municipalities and define the levels of risk. The southeast region of the São Paulo State was considered to be the area with a higher number of municipalities classified as high risk and should be considered a priority for the application of prevention measures against Yellow Fever.

Infection and immune-mediated meningococcal-associated arthritis: combination features in the same patient

We present a case of a 16-year-old male patient with sudden-onset, rash, arthritis and meningitis by Neisseria meningitidis one week after an acute upper respiratory infection. On the 10th day of treatment followed by neurological and arthritis clinical improvement, he presented once again a tender and swollen left knee with a moderate effusion, and active and passive range of motion was severely limited secondary to pain, and when he was submitted to surgical drainage and synovial fluid analysis he showed inflammatory characteristics. A non-steroidal anti-inflammatory drug was taken for five days with complete improvement of symptoms. The case is notable for its combination of features of septic and immune-mediated arthritis, which has rarely been reported in the same patient.