Parallel method for sparse semisupervised hyperspectral unmixing

Parallel hyperspectral unmixing problem is considered in this paper. A semisupervised approach is developed under the linear mixture model, where the abundance's physical constraints are taken into account. The proposed approach relies on the increasing availability of spectral libraries of materials measured on the ground instead of resorting to endmember extraction methods. Since Libraries are potentially very large and hyperspectral datasets are of high dimensionality a parallel implementation in a pixel-by-pixel fashion is derived to properly exploits the graphics processing units (GPU) architecture at low level, thus taking full advantage of the computational power of GPUs. Experimental results obtained for real hyperspectral datasets reveal significant speedup factors, up to 164 times, with regards to optimized serial implementation.

Parallel hyperspectral unmixing method on GPU

Many Hyperspectral imagery applications require a response in real time or near-real time. To meet this requirement this paper proposes a parallel unmixing method developed for graphics processing units (GPU). This method is based on the vertex component analysis (VCA), which is a geometrical based method highly parallelizable. VCA is a very fast and accurate method that extracts endmember signatures from large hyperspectral datasets without the use of any a priori knowledge about the constituent spectra. Experimental results obtained for simulated and real hyperspectral datasets reveal considerable acceleration factors, up to 24 times.

Parallel sparse unmixing of hyperspectral data

In this paper, a new parallel method for sparse spectral unmixing of remotely sensed hyperspectral data on commodity graphics processing units (GPUs) is presented. A semi-supervised approach is adopted, which relies on the increasing availability of spectral libraries of materials measured on the ground instead of resorting to endmember extraction methods. This method is based on the spectral unmixing by splitting and augmented Lagrangian (SUNSAL) that estimates the material's abundance fractions. The parallel method is performed in a pixel-by-pixel fashion and its implementation properly exploits the GPU architecture at low level, thus taking full advantage of the computational power of GPUs. Experimental results obtained for simulated and real hyperspectral datasets reveal significant speedup factors, up to 1 64 times, with regards to optimized serial implementation.

Algebraic structure for interaction on mixed models

Binary operations on commutative Jordan algebras, CJA, can be used to study interactions between sets of factors belonging to a pair of models in which one nests the other. It should be noted that from two CJA we can, through these binary operations, build CJA. So when we nest the treatments from one model in each treatment of another model, we can study the interactions between sets of factors of the first and the second models.

From weak Allee effect to no Allee effect in Richards’ growth models

Population dynamics have been attracting interest since many years. Among the considered models, the Richards’ equations remain one of the most popular to describe biological growth processes. On the other hand, Allee effect is currently a major focus of ecological research, which occurs when positive density dependence dominates at low densities. In this chapter, we propose the dynamical study of classes of functions based on Richards’ models describing the existence or not of Allee effect. We investigate bifurcation structures in generalized Richards’ functions and we look for the conditions in the (β, r) parameter plane for the existence of a weak Allee effect region. We show that the existence of this region is related with the existence of a dovetail structure. When the Allee limit varies, the weak Allee effect region disappears when the dovetail structure also disappears. Consequently, we deduce the transition from the weak Allee effect to no Allee effect to this family of functions. To support our analysis, we present fold and flip bifurcation curves and numerical simulations of several bifurcation diagrams.

Symbolic dynamics and big bang bifurcation in Weibull-Gompertz-Fréchet's growth models

In this paper, motivated by the interest and relevance of the study of tumor growth models, a central point of our investigation is the study of the chaotic dynamics and the bifurcation structure of Weibull-Gompertz-Fréchet's functions: a class of continuousdefined one-dimensional maps. Using symbolic dynamics techniques and iteration theory, we established that depending on the properties of this class of functions in a neighborhood of a bifurcation point PBB, in a two-dimensional parameter space, there exists an order regarding how the infinite number of periodic orbits are born: the Sharkovsky ordering. Consequently, the corresponding symbolic sequences follow the usual unimodal kneading sequences in the topological ordered tree. We verified that under some sufficient conditions, Weibull-Gompertz-Fréchet's functions have a particular bifurcation structure: a big bang bifurcation point PBB. This fractal bifurcations structure is of the so-called "box-within-a-box" type, associated to a boxe ω1, where an infinite number of bifurcation curves issues from. This analysis is done making use of fold and flip bifurcation curves and symbolic dynamics techniques. The present paper is an original contribution in the framework of the big bang bifurcation analysis for continuous maps.

Explosion birth and extinction: Double big bang bifurcations and Allee effect in Tsoularis-Wallace's growth models

This work concerns dynamics and bifurcations properties of a new class of continuous-defined one-dimensional maps: Tsoularis-Wallace's functions. This family of functions naturally incorporates a major focus of ecological research: the Allee effect. We provide a necessary condition for the occurrence of this phenomenon of extinction. To establish this result we introduce the notions of Allee's functions, Allee's effect region and Allee's bifurcation curve. Another central point of our investigation is the study of bifurcation structures for this class of functions, in a three-dimensional parameter space. We verified that under some sufficient conditions, Tsoularis-Wallace's functions have particular bifurcation structures: the big bang and the double big bang bifurcations of the so-called "box-within-a-box" type. The double big bang bifurcations are related to the existence of flip codimension-2 points. Moreover, it is verified that these bifurcation cascades converge to different big bang bifurcation curves, where for the corresponding parameter values are associated distinct kinds of boxes. This work contributes to clarify the big bang bifurcation analysis for continuous maps and understand their relationship with explosion birth and extinction phenomena.

Modern techniques for constraint solving the CASPER experience

Dissertação apresentada para obtenção do Grau de Doutor em Engenharia Informática, pela Universidade Nova de Lisboa, Faculdade de Ciências e Tecnologia

Tectónica da região de Sintra

Photo-interpretation of aerial stereopairs of the Sintra region on the approx. 1/32 000 scale together with field work allowed the production of the present Tectonic Map of the Sintra region. It is now possible to separate structures which resulted from two different tectonic events: one, corresponding to the intrusion of the Late Cretaceous Sintra igneous diapir, and the other the Miocene compressive event, the most important tectonic inversion phase of the Lusitanian Basin. The former are present to the south, southeast and east of the intrusion and within the intrusion itself, affecting the peripheral granites and their contacts with the gabbro-syenite core. These structures comprehend: i) faults and conical fractures striking parallel to the massif boundary, which were intruded by dykes, ii) vertical faults and fractures of two conjugate sets, dextral NNW-SSE and sinistral NNE-SSW. These faults are certainly associated with the E-W striking massif's northwards directed thrust and indicate a N-S oriented horizontal maximum compressive stress. The Miocene compressive event reactivated most of the inherited structures as follows. The NNWSSE faults located on the Sintra southern platform were reactivated as dextral strike slip faults and the E-W thrust along the northern boundary of the massif was also reactivated. This thrust propagated to the east. It also enhanced the asymmetry of the rim-syncline, uplifted the massif and reactivated the NNE-SSW faults as sinistral lateral ramps, which also accommodated vertical throw. The present Tectonic Map of Sintra together with the available geophysical data (MOREIRA, 1984, KULLBERG et al., 1991, SILVA & MIRANDA, 1994) allowed reassessment of the models proposed for the emplacement of the Sintra, Sines and Monchique igneous massifs, which intruded during Late Cretaceous times along the deep dextral NNW-SSE oriented strike slip fault (RIBEIRO et al., 1979; TERRINHA, 1998; TERRINHA & KULLBERG, 1998).

An integer programming framework for sequencing cutting patterns based on interval graph completion

We derived a framework in integer programming, based on the properties of a linear ordering of the vertices in interval graphs, that acts as an edge completion model for obtaining interval graphs. This model can be applied to problems of sequencing cutting patterns, namely the minimization of open stacks problem (MOSP). By making small modifications in the objective function and using only some of the inequalities, the MOSP model is applied to another pattern sequencing problem that aims to minimize, not only the number of stacks, but also the order spread (the minimization of the stack occupation problem), and the model is tested.

An integer programming model for the minimum interval graph completion problem

The minimum interval graph completion problem consists of, given a graph G = ( V, E ), finding a supergraph H = ( V, E ∪ F ) that is an interval graph, while adding the least number of edges |F| . We present an integer programming formulation for solving the minimum interval graph completion problem recurring to a characteri- zation of interval graphs that produces a linear ordering of the maximal cliques of the solution graph.

Multi-objective Parallel Particle Swarm Optimization for Day-ahead Vehicle-To-Grid Scheduling

This paper presents a methodology for multi-objective day-ahead energy resource scheduling for smart grids considering intensive use of distributed generation and Vehicle- To-Grid (V2G). The main focus is the application of weighted Pareto to a multi-objective parallel particle swarm approach aiming to solve the dual-objective V2G scheduling: minimizing total operation costs and maximizing V2G income. A realistic mathematical formulation, considering the network constraints and V2G charging and discharging efficiencies is presented and parallel computing is applied to the Pareto weights. AC power flow calculation is included in the metaheuristics approach to allow taking into account the network constraints. A case study with a 33-bus distribution network and 1800 V2G resources is used to illustrate the performance of the proposed method.

Anthropocentric-based robotic and autonomous systems: assessment for new organisational options

Text based on the paper presented at the Conference "Autonomous systems: inter-relations of technical and societal issues" held at Monte de Caparica (Portugal), Universidade Nova de Lisboa, November, 5th and 6th 2009 and organized by IET-Research Centre on Enterprise and Work Innovation

Increase of the Delivered Power Probability in Distribution Networks using Pareto DC Programming

A methodology to increase the probability of delivering power to any load point through the identification of new investments in distribution network components is proposed in this paper. The method minimizes the investment cost as well as the cost of energy not supplied in the network. A DC optimization model based on mixed integer non-linear programming is developed considering the Pareto front technique in order to identify the adequate investments in distribution networks components which allow increasing the probability of delivering power for any customer in the distribution system at the minimum possible cost for the system operator, while minimizing the energy not supplied cost. Thus, a multi-objective problem is formulated. To illustrate the application of the proposed methodology, the paper includes a case study which considers a 180 bus distribution network

Polyglutamine diseases (PolyQ): construction and characterization of yeast models

Dissertação apresentada para a obtenção do Grau de Mestre em Genética Molecular e Biomedicina, pela Universidade Nova de Lisboa, Faculdade de Ciências e Tecnologia