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.

Fast unsupervised extraction of endmembers spectra from hyperspectral data

Linear unmixing decomposes an hyperspectral image into a collection of re ectance spectra, called endmember signatures, and a set corresponding abundance fractions from the respective spatial coverage. This paper introduces vertex component analysis, an unsupervised algorithm to unmix linear mixtures of hyperpsectral data. VCA exploits the fact that endmembers occupy vertices of a simplex, and assumes the presence of pure pixels in data. VCA performance is illustrated using simulated and real data. VCA competes with state-of-the-art methods with much lower computational complexity.

Hyperspectral unmixing algorithm via dependent component analysis

This paper introduces a new method to blindly unmix hyperspectral data, termed dependent component analysis (DECA). This method decomposes a hyperspectral images into a collection of reflectance (or radiance) spectra of the materials present in the scene (endmember signatures) and the corresponding abundance fractions at each pixel. DECA assumes that each pixel is a linear mixture of the endmembers signatures weighted by the correspondent abundance fractions. These abudances are modeled as mixtures of Dirichlet densities, thus enforcing the constraints on abundance fractions imposed by the acquisition process, namely non-negativity and constant sum. The mixing matrix is inferred by a generalized expectation-maximization (GEM) type algorithm. This method overcomes the limitations of unmixing methods based on Independent Component Analysis (ICA) and on geometrical based approaches. The effectiveness of the proposed method is illustrated using simulated data based on U.S.G.S. laboratory spectra and real hyperspectral data collected by the AVIRIS sensor over Cuprite, Nevada.

On the analysis of data emerging in non-linear and complex systemscomparison of two X-Ray prony spectra

New arguments proving that successive (repeated) measurements have a memory and actually remember each other are presented. The recognition of this peculiarity can change essentially the existing paradigm associated with conventional observation in behavior of different complex systems and lead towards the application of an intermediate model (IM). This IM can provide a very accurate fit of the measured data in terms of the Prony's decomposition. This decomposition, in turn, contains a small set of the fitting parameters relatively to the number of initial data points and allows comparing the measured data in cases where the “best fit” model based on some specific physical principles is absent. As an example, we consider two X-ray diffractometers (defined in paper as A- (“cheap”) and B- (“expensive”) that are used after their proper calibration for the measuring of the same substance (corundum a-Al2O3). The amplitude-frequency response (AFR) obtained in the frame of the Prony's decomposition can be used for comparison of the spectra recorded from (A) and (B) - X-ray diffractometers (XRDs) for calibration and other practical purposes. We prove also that the Fourier decomposition can be adapted to “ideal” experiment without memory while the Prony's decomposition corresponds to real measurement and can be fitted in the frame of the IM in this case. New statistical parameters describing the properties of experimental equipment (irrespective to their internal “filling”) are found. The suggested approach is rather general and can be used for calibration and comparison of different complex dynamical systems in practical purposes.

Application of Infrared Thermal Imaging in a Violinist with Temporomandibular Disorder

Temporomandibular disorders (TMD) consist of a group of pathologies that affect the masticatory muscles, temporomandibular joints (TMJ), and/or related structures. String instrumentalists, like many orchestra musicians, can spend hours with head postures that may influence the biomechanical behavior of the TMJ and the muscles of the craniocervicomandibular complex (CCMC). The adoption of abnormal postures acquired during performance by musicians can lead to muscular hyperactivity of the head and cervical muscles, with the possible appearance of TMD. Medical infrared thermography is a non-invasive procedure that can monitor the changes in the superficial tissue related to blood circulation and may serve as a complement to the clinical examination. The objective of this study was to use infrared thermography to evaluate, in one subject, the cutaneous thermal changes adjacent to the CCMC that occur before, during, and after playing a string instrument.

Formation of a stable cyano-bridged dinuclear iron cluster following oxidation of the superoxide reductases from Treponema pallidum and Desulfovibrio vulgaris with K(3)Fe(CN)(6)

Inorg. Chem., 2003, 42 (4), pp 938–940 DOI: 10.1021/ic0262886

X-ray crystal structure and EPR spectra of "arsenite-inhibited" Desulfovibriogigas aldehyde dehydrogenase: a member of the xanthine oxidase family

J. Am. Chem. Soc., 2004, 126 (28), pp 8614–8615 DOI: 10.1021/ja0490222

Electron impact ionization cross sections and analysis of X-ray spectra

A thesis submitted for the degree of Ph. D. in Physics

Thermal performance of cool facades: evaluation by Infrared Thermography

High reflective paints (cool paints) are used on flat roofs to reduce heat gains from the incidence of solar radiation and thus improve the thermal comfort and energy efficiency of buildings, especially in summer periods. Given the application potential of these paints on vertical surfaces, a research study has been developed to evaluate the thermal performance of reflective paints on walls under real exposure conditions. Accordingly, different reflective paints have been applied as the final coating of an ETICS type solution, on the facades of a full scale experimental cell built at LNEC campus. For being applied in an ETICS system a paint has to fulfill several requirements, whether aesthetic or functional (such as the adhesion between the coating layers or the durability of the insulation), essential for its efficient performance. Since this construction coating system is subject to a prolonged sun exposure, various problems may arise, such as paint degradation or deterioration of the thermal insulation properties, particularly when dark colors are applied. To evaluate the thermal performance of the chosen paints, the method of non-destructive analysis by Infrared Thermography was used. Thermography allows knowing the temperature distribution of facades by measuring the radiation emitted by their surfaces. To complement the thermographic diagnosis, thermocouples were placed between the insulation and the paint system of the experimental cell. Additional laboratory tests allowed the characterization of the optical properties (reflectance and emittance) of the different reflective paints used in this study. The comparative analysis of the thermal performance of reflective and conventional paints revealed that the reflective paint allows a reduction of the facade surface temperature, reducing the risk of loss of insulating properties of the ETICS system and thus ensuring its longevity and functionality. The color of the paint used affects, naturally, the reflective ability of the surface and may have an important role in energy balance of the building. This paper also showed the potential of infrared thermography in the evaluation of the thermal performance of reflective paints.

Automated quantitative analysis of MCC-IMS spectra

Ion Mobility Spectrometry coupled with Multi Capillary Columns (MCC -IMS) is a fast analytical technique working at atmospheric pressure with high sensitivity and selectivity making it suitable for the analysis of complex biological matrices. MCC-IMS analysis generates its information through a 3D spectrum with peaks, corresponding to each of the substances detected, providing quantitative and qualitative information. Sometimes peaks of different substances overlap, making the quantification of substances present in the biological matrices a difficult process. In the present work we use peaks of isoprene and acetone as a model for this problem. These two volatile organic compounds (VOCs) that when detected by MCC-IMS produce two overlapping peaks. In this work it’s proposed an algorithm to identify and quantify these two peaks. This algorithm uses image processing techniques to treat the spectra and to detect the position of the peaks, and then fits the data to a custom model in order to separate the peaks. Once the peaks are separated it calculates the contribution of each peak to the data.

Production of polyhydroxyalkanoates from oil-containing substrates

Different oil-containing substrates, namely, used cooking oil (UCO), fatty acids-byproduct from biodiesel production (FAB) and olive oil deodorizer distillate (OODD) were tested as inexpensive carbon sources for the production of polyhydroxyalkanoates (PHA) using twelve bacterial strains, in batch experiments. The OODD and FAB were exploited for the first time as alternative substrates for PHA production. Among the tested bacterial strains, Cupriavidus necator and Pseudomonas resinovorans exhibited the most promising results, producing poly-3-hydroxybutyrate, P(3HB), form UCO and OODD and mcl-PHA mainly composed of 3-hydroxyoctanoate (3HO) and 3-hydroxydecanoate (3HD) monomers from OODD, respectively. Afterwards, these bacterial strains were cultivated in bioreactor. C. necator were cultivated in bioreactor using UCO as carbon source. Different feeding strategies were tested for the bioreactor cultivation of C. necator, namely, batch, exponential feeding and DO-stat mode. The highest overall PHA productivity (12.6±0.78 g L-1 day-1) was obtained using DO-stat mode. Apparently, the different feeding regimes had no impact on polymer thermal properties. However, differences in polymer‟s molecular mass distribution were observed. C. necator was also tested in batch and fed-batch modes using a different type of oil-containing substrate, extracted from spent coffee grounds (SCG) by super critical carbon dioxide (sc-CO2). Under fed-batch mode (DO-stat), the overall PHA productivity were 4.7 g L-1 day-1 with a storage yield of 0.77 g g-1. Results showed that SCG can be a bioresource for production of PHA with interesting properties. Furthermore, P. resinovorans was cultivated using OODD as substrate in bioreactor under fed-batch mode (pulse feeding regime). The polymer was highly amorphous, as shown by its low crystallinity of 6±0.2%, with low melting and glass transition temperatures of 36±1.2 and -16±0.8 ºC, respectively. Due to its sticky behavior at room temperature, adhesiveness and mechanical properties were also studied. Its shear bond strength for wood (67±9.4 kPa) and glass (65±7.3 kPa) suggests it may be used for the development of biobased glues. Bioreactor operation and monitoring with oil-containing substrates is very challenging, since this substrate is water immiscible. Thus, near-infrared spectroscopy (NIR) was implemented for online monitoring of the C. necator cultivation with UCO, using a transflectance probe. Partial least squares (PLS) regression was applied to relate NIR spectra with biomass, UCO and PHA concentrations in the broth. The NIR predictions were compared with values obtained by offline reference methods. Prediction errors to these parameters were 1.18 g L-1, 2.37 g L-1 and 1.58 g L-1 for biomass, UCO and PHA, respectively, which indicates the suitability of the NIR spectroscopy method for online monitoring and as a method to assist bioreactor control. UCO and OODD are low cost substrates with potential to be used in PHA batch and fed-batch production. The use of NIR in this bioprocess also opened an opportunity for optimization and control of PHA production process.

Electro-Optical/Infrared sensor turret integration on an aircraft - structural im-pact on LOCKHEED MARTIN C-130H and design methodology

The growing need to patrol and survey large maritime and terrestrial areas increased the need to integrate external sensors on aircraft in order to accomplish those patrols at increasingly higher altitudes, longer range and not depending upon vehicle type. The main focus of this work is to elaborate a practical, simple, effective and efficient methodology for the aircraft modification procedure resulting from the integration of an Elec-tro-Optical/Infra-Red (EO/IR) turret through a support structure. The importance of the devel-opment of a good methodology relies on the correct management of project variables as time, available resources and project complexity. The key is to deliver a proper tool for a project de-sign team that will be used to create a solution that fulfils all technical, non-technical and certi-fication requirements present in this field of transportation. The created methodology is inde-pendent of two main inputs: sensor model and aircraft model definition, and therefore it is in-tended to deliver the results for different projects besides the one that was presented in this work as a case study. This particular case study presents the development of a structure support for FLIR STAR SAPHIRE III turret integration on the front lower fuselage bulkhead (radome) of the LOCKHEED MARTIN C-130 H. Development of the case study focuses on the study of local structural analysis through the use of Finite Element Method (FEM). Development of this Dissertation resulted in a cooperation between Faculty of Science and Technology - Universidade Nova de Lisboa and the company OGMA - Indústria Aeronáutica de Portugal

Bioactive dihydroxyfuranonaphthoquinones from the bark of Tabebuia incana A.H. Gentry (Bignoniaceae) and HPLC analysis of commercial pau d'arco and certified T.incana bark infusions

Tabebuia incana A.H. Gentry (Bignoniaceae) is a tree from the Brazilian Amazon having medicinal uses and is one several Tabebuia spp. known as pau d'arco or palo de arco in this region. Fractionation of the bark ethanolic extract afforded a mixture of 5 and 8-hydroxy-2-(1-hydroxyethyl)naphtho[2,3-b]furan-4,9-diones (1 and 2, respectively) identified on the basis of nuclear magnetic resonance (NMR), infrared (IR) and mass (MS) spectra, whose in vitro antimalarial and antitumor activity have been shown previously. This is the first study on T. incana bark, and 2 are described in this species for the first time. Also, high performance liquid chromatography (HPLC) analysis of T. incana bark tea revealed the presence of the 1 + 2 mixture peak corresponding to a concentration in the range 10-6-10-5 M. The chromatograms of teas prepared from commercial pau d' arco and T. incana bark were also studied and the presence of the 1 + 2 peak has potential for quality control of commercial plant materials.

Surface plasmon resonance-coupled photoluminescence and resistive switching behavior of pulsed laser-deposited Ag:SiC nanocermet thin films

In this study, Ag:SiC nanocermets were prepared via rapid thermal annealing (RTA) of pulsed laser-deposited SiC/Ag/SiC trilayers grown on Si substrate. Atomic force microscope images show that silver nanoparticles (Ag NPs) are formed after RTA, and the size of NPs increases with increasing Ag deposition time (t Ag). Sharp dip observed in the reflectance spectra confirmed the existence of Ag surface plasmons (SPs). The infrared transmission spectra showed an intense and broad absorption band around 780–800 cm−1 that can be assigned to Si-C stretching vibration mode. Influence of t Ag on the spectral characteristics of SP-enhanced photoluminescence (PL) and electrical properties of silicon carbide (SiC) films has been investigated. The maximum PL enhancement by 5.5 times for Ag:SiC nanocermets is achieved when t Ag ≈ 50 s. This enhancement is due to the strong resonant coupling between SiC and the SP oscillations of the Ag NPs. Presence of Ag NPs in SiC also induces a forming-free resistive switching with switching ratio of 2 × 10−2. The analysis of I–V curves demonstrates that the trap-controlled space-charge-limited conduction with filamentary model is the governing mechanism for the resistive switching in nanocerment thin films.

Measurement of charged-particle spectra in Pb+Pb collisions at √sNN = 2.76TeV with the ATLAS detector at the LHC

Charged-particle spectra obtained in 0.15nb−1 of Pb+Pb interactions at sNN−−−√=2.76TeV and 4.2pb−1 of pp interactions at s√=2.76TeV with the ATLAS detector at the LHC are presented in a wide transverse momentum (0.5spectra are presented as a function of collision centrality, which is determined by the response of the forward calorimeter located on both sides of the interaction point. The nuclear modification factors RAA and RCP are presented in detail as function of centrality, pT and η. They show a distinct pT-dependence with a pronounced minimum at about 7GeV. Above 60GeV, RAA is consistent with a plateau at a centrality-dependent value, within the uncertainties. The value is 0.55±0.01(stat.)±0.04(syst.) in the most central collisions. The RAA distribution is consistent with flat |η| dependence over the whole transverse momentum range in all centrality classes.