DSA image registration using non-uniform MRF model and pivotal control points

In order to reduce the motion artifacts in DSA, non-rigid image registration is commonly used before subtracting the mask from the contrast image. Since DSA registration requires a set of spatially non-uniform control points, a conventional MRF model is not very efficient. In this paper, we introduce the concept of pivotal and non-pivotal control points to address this, and propose a non-uniform MRF for DSA registration. We use quad-trees in a novel way to generate the non-uniform grid of control points. Our MRF formulation produces a smooth displacement field and therefore results in better artifact reduction than that of registering the control points independently. We achieve improved computational performance using pivotal control points without compromising on the artifact reduction. We have tested our approach using several clinical data sets, and have presented the results of quantitative analysis, clinical assessment and performance improvement on a GPU. (C) 2013 Elsevier Ltd. All rights reserved.

Formation flying of small satellites using suboptimal MPSP guidance

Using the recently developed model predictive static programming (MPSP), a suboptimal guidance logic is presented in this paper for formation flying of small satellites. Due to the inherent nature of the problem formulation, MPSP does not require the system dynamics to be linearized. The proposed guidance scheme is valid both for high eccentricity chief satellite orbits as well as large separation distance between chief and deputy satellites. Moreover, since MPSP poses the desired conditions as a set of `hard constraints', the final accuracy level achieved is very high. The proposed guidance scheme has been tested successfully for a variety of initial conditions and for a variety of formation commands as well. Comparison with standard Linear Quadratic Regulator (LQR) solution (which serves as a guess solution for MPSP) and another nonlinear controller, State Dependent Riccati Equation (SDRE) reveals that MPSP guidance achieves the objective with higher accuracy and with lesser amount of control usage as well.

Effects of elastase and collagenase on the nonlinearity and anisotropy of porcine aorta

We use enzymatic manipulation methods to investigate the individual and combined roles of elastin and collagen on arterial mechanics. Porcine aortic tissues were treated for differing amounts of time using enzymes elastase and collagenase to cause degradation in substrate proteins elastin and collagen and obtain variable tissue architecture. We use equibiaxial mechanical tests to quantify the material properties of control and enzyme treated tissues and histological methods to visualize the underlying tissue microstructure in arterial tissues. Our results show that collagenase treated tissues were more compliant in the longitudinal direction as compared to control tissues. Collagenase treatment also caused a decrease in the tissue nonlinearity as compared to the control samples in the study. A one hour collagenase treatment was sufficient to cause fragmentation and degradation of the adventitial collagen. In contrast, elastase treatment leads to significantly stiffer tissue response associated with fragmented and incomplete elastin networks in the tissue. Thus, elastin in arterial walls distributes tensile stresses whereas collagen serves to reinforce the vessel wall in the circumferential direction and also contributes to tissue anisotropy. A microstructurally motivated strain energy function based on circumferentially oriented medial fibers and helically oriented collagen fibers in the adventitia is useful in describing these experimental results.

Optical properties of an inhomogeneously broadened Lambda V-system with multiple excited states

We present a theoretical model using a density matrix approach to show the influence of multiple excited states on the optical properties of an inhomogeneously broadened Lambda V-system of the Rb-87 D2 line. These closely spaced multiple excited states cause asymmetry in absorption and dispersion profiles. We observe the reduced absorption profiles, due to dressed state interactions of the applied electromagnetic fields, which results the Mollow sideband-like transparency windows. In a room temperature vapor, we obtain a narrow enhanced absorption and steep positive dispersion at the line center when the strengths of control and pump fields are equal. Here, we show how the probe transmittance varies when it passes through the atomic medium. We also discuss the transient behavior of our system which agrees well with the corresponding absorption and dispersion profiles. This study has potential applications in controllability of group velocity, and for optical and quantum information processing.

Experimental assessment of concrete damage due to exposure to high temperature and efficacy of the repair system

The present study experimentally evaluates the performance of control (standard cylinder specimen), damaged (mechanical loading after thermal exposure) and repaired / retrofitted normal plain concrete cylinders using different repair schemes such as on use of FRP wraps, Geo-polymers, etc., to restore the capacity of damaged structural concrete elements. The control-companion specimen in the series provides the reference frame against which both, specimen damage levels were quantified and the benefits of a specimen repaired subsequent to damage were assessed.

Wall-mode instability in plane shear flow of viscoelastic fluid over a deformable solid

The linear stability analysis of a plane Couette flow of an Oldroyd-B viscoelastic fluid past a flexible solid medium is carried out to investigate the role of polymer addition in the stability behavior. The system consists of a viscoelastic fluid layer of thickness R, density rho, viscosity eta, relaxation time lambda, and retardation time beta lambda flowing past a linear elastic solid medium of thickness HR, density rho, and shear modulus G. The emphasis is on the high-Reynolds-number wall-mode instability, which has recently been shown in experiments to destabilize the laminar flow of Newtonian fluids in soft-walled tubes and channels at a significantly lower Reynolds number than that for flows in rigid conduits. For Newtonian fluids, the linear stability studies have shown that the wall modes become unstable when flow Reynolds number exceeds a certain critical value Re c which scales as Sigma(3/4), where Reynolds number Re = rho VR/eta, V is the top-plate velocity, and dimensionless parameter Sigma = rho GR(2)/eta(2) characterizes the fluid-solid system. For high-Reynolds-number flow, the addition of polymer tends to decrease the critical Reynolds number in comparison to that for the Newtonian fluid, indicating a destabilizing role for fluid viscoelasticity. Numerical calculations show that the critical Reynolds number could be decreased by up to a factor of 10 by the addition of small amount of polymer. The critical Reynolds number follows the same scaling Re-c similar to Sigma(3/4) as the wall modes for a Newtonian fluid for very high Reynolds number. However, for moderate Reynolds number, there exists a narrow region in beta-H parametric space, corresponding to very dilute polymer solution (0.9 less than or similar to beta < 1) and thin solids (H less than or similar to 1.1), in which the addition of polymer tends to increase the critical Reynolds number in comparison to the Newtonian fluid. Thus, Reynolds number and polymer properties can be tailored to either increase or decrease the critical Reynolds number for unstable modes, thus providing an additional degree of control over the laminar-turbulent transition.

Experimental studies on the influence of intermediate principal stress and inclination on the mechanical behaviour of angular sands

A comprehensive experimental study has been made on angular sand to investigate various aspects of mechanical behavior. A hollow cylinder torsion testing apparatus is used in this program to apply a range of stress conditions on this angular quartzitic fine sand under monotonic drained shear. The effect of the magnitude and inclination of the principal stresses on an element of sand is studied through these experiments. This magnitude and inclination of the principal stresses are presented as an ``ensemble measure of fabric in sands''. This ensemble measure of fabric in the sands evolves through the shearing process, and reaches the final state, which indeed has a unique fabric. The sand shows significant variation in strength with changing inclination of the principal stresses. The locus of the final stress state in principal stress space is also mapped from these series of experiments. Additional aspects of non-coaxiality, a benchmarking exercise with a few constitutive models is presented here. This experimental approach albeit indirect shows that a unique state which is dependent on the fabric, density and confining stress exists. This suite of experiments provides a well-controlled data set for a clear understanding on the mechanical behavior of sands.

The Rise and Fall of Water Net (Hydrodictyon reticulatum) in New Zealand

During the late 1980s to early 1990s a range of aquatic habitats in the central North Island of New Zealand were invaded by the filamentous green alga, water net Hydrodictyon reticulatum (Linn. Lagerheim). The alga caused significant economic and recreational impacts at major sites of infestation, but it was also associated with enhanced invertebrate numbers and was the likely cause of an improvement in the trout fishery. The causes of prolific growth of water net and the range of control options pursued are reviewed. The possible causes of its sudden decline in 1995 are considered, including physical factors, increase in grazer pressure, disease, and loss of genetic vigour.

Open Personalization: Involving Third Parties in Improving the User Experience of Websites

Traditional software development captures the user needs during the requirement analysis. The Web makes this endeavour even harder due to the difficulty to determine who these users are. In an attempt to tackle the heterogeneity of the user base, Web Personalization techniques are proposed to guide the users’ experience. In addition, Open Innovation allows organisations to look beyond their internal resources to develop new products or improve existing processes. This thesis sits in between by introducing Open Personalization as a means to incorporate actors other than webmasters in the personalization of web applications. The aim is to provide the technological basis that builds up a trusty environment for webmasters and companion actors to collaborate, i.e. "an architecture of participation". Such architecture very much depends on these actors’ profile. This work tackles three profiles (i.e. software partners, hobby programmers and end users), and proposes three "architectures of participation" tuned for each profile. Each architecture rests on different technologies: a .NET annotation library based on Inversion of Control for software partners, a Modding Interface in JavaScript for hobby programmers, and finally, a domain specific language for end-users. Proof-of-concept implementations are available for the three cases while a quantitative evaluation is conducted for the domain specific language.

Validation of daily growth increments and estimation of growth rates of larval and early juvenile black skipjack, Euthynnus lineatus, using otoliths

ENGLISH: Increments in otoliths (sagittae) were examined, using light and scanning electron microscopy, to determine ages and estimate growth rates of larval and early-juvenile black skipjack, Euthynnus lineatus. Larvae and juveniles were collected between 1987 and 1989 from coastal waters of Panama in the eastern Pacific Ocean. Results from a laboratory experiment indicated that immersion for 6 and 12 hours in a 200 mg/L solution of tetracycline hydrochloride adequately marks otoliths and that increments are formed daily in the sagittae of postflexion larvae and early juveniles. Further, survival rates of tetracycline-treated fish were not significantly different from those of control fish. Growth rates were derived from length-age relationships of 218 field-collected specimens ranging in size from 5.7 to 20.3 mm SL. A growth rate of 0.70 mm/d was estimated from the weighted regression of standard length on age for all specimens. This rate lies within the range reported for larvae and early juveniles of other species of subtropical and tropical scombrids. Growth rates of postflexion larvae and early juveniles were not significantly different between the rainy season in July-August 1988 and the dry, upwelling season in January-February 1989. Growth was, however, significantly more variable for older individuals in July-August than in January-February, and may correspond, in part, to seasonal patchiness of prey. The growth rates of the otoliths relative to fish length were also not significantly different between seasons; however, the otoliths were larger relative to the lengths of fish collected in the rainy season, which may reflect slower growth during earlier larval stages. SPANISH: Se examinaron incrementos en otolitos (ságitas), usando microscopia de luz y de barrido electrónico, a fin de determinar la edad y estimar las tasas de crecimiento de barriletes negros, Euthynnus lineatus, larvales y juveniles tempranos. Entre 1987 y 1989 se capturaron larvas y juveniles en las aguas costeras de Panamá en el Océano Pacífico oriental. Los resultados de un experimento de laboratorio indicaron que una inmersión de 6 a 12 horas de duración en una solución de 200 mg/L de hidrocloro de tetraciclina marca los otolitos adecuadamente y que los incrementos se forman a diario en las ságitas de larvas en postflexión y juveniles tempranos. Además, las tasas de supervivencia de los peces tratados con tetraciclina no fueron significativamente diferentes a aquellas de los peces de control. Se calcularon las tasas de crecimiento a partir de las relaciones de talla-edad de 218 especímenes de TE entre 5.7 y 20.3 mm capturados en el mar. Se estimó.una tasa de crecimiento de 0.70 mm/día a partir de la regresión ponderada de talla estándar sobre edad para todos los especímenes. Esta tasa cae dentro del rango reportado para larvas y juveniles tempranos de otras especies de escómbridos subtropicales y tropicales. Las tasas de crecimiento de larvas en postflexión y juveniles tempranos no fueron significativamente diferentes entre la temporada de lluvias en julio-agosto de 1988 y la temporada de sequía y afloramiento en enero-febrero de 1989. Sin emoargo, el crecimiento fue significativamente más variable para los individuos de mayor edad en julio-agosto que en enero-febrero, y quizás corresponda parcialmente a la irregularidad temporal de la abundancia de presas. Las tasas de crecimiento de los otolitos en relación a la talla de los peces tampoco fueron significativamente diferentes entre temporadas; sin embargo, los otolitos eran más grandes en relación a la talla en peces capturados en la temporada de lluvias, lo cual podría reflejar crecimiento más lento durante las etapas larvales más tempranas. (PDF contains 42 pages.)

Lake Kainji water hyacinth infestation and the control strategies

The paper documents the aquatic vascular plants situation of Lake Kainji, particularly with the invasion of water hyacinth (Eichhormia crassipes) hitherto alien to the lake system. The frequency of occurrence of some aquatic plants on Lake Kainji have increased whilst a few decreased over the period between 1984 and 1995. More than 53.6% of the lake surface has now been covered by aquatic vascular plants of which water hyacinth alone covered 30.6% (16.4% of the system). Apart from the characteristic rapid multiplication and biomass build-up of the plant, more than 9.1 clumps with a diameter range between 0.8-4cm enter the lake per minute. Thus more than 16.4ha per day or annually 6000ha of water hyacinth are added as drifts. Because of the peculiar hydrological cycle of the lake and the plant aggressiveness and invasiveness, NIFFR (with the support of the GTZ) had developed an approach of control which combines manual, biological, natural, boom construction and chemical measures. These are discussed

Neural routing circuits for forming invariant representations of visual objects

This thesis presents a biologically plausible model of an attentional mechanism for forming position- and scale-invariant representations of objects in the visual world. The model relies on a set of control neurons to dynamically modify the synaptic strengths of intra-cortical connections so that information from a windowed region of primary visual cortex (Vl) is selectively routed to higher cortical areas. Local spatial relationships (i.e., topography) within the attentional window are preserved as information is routed through the cortex, thus enabling attended objects to be represented in higher cortical areas within an object-centered reference frame that is position and scale invariant. The representation in V1 is modeled as a multiscale stack of sample nodes with progressively lower resolution at higher eccentricities. Large changes in the size of the attentional window are accomplished by switching between different levels of the multiscale stack, while positional shifts and small changes in scale are accomplished by translating and rescaling the window within a single level of the stack. The control signals for setting the position and size of the attentional window are hypothesized to originate from neurons in the pulvinar and in the deep layers of visual cortex. The dynamics of these control neurons are governed by simple differential equations that can be realized by neurobiologically plausible circuits. In pre-attentive mode, the control neurons receive their input from a low-level "saliency map" representing potentially interesting regions of a scene. During the pattern recognition phase, control neurons are driven by the interaction between top-down (memory) and bottom-up (retinal input) sources. The model respects key neurophysiological, neuroanatomical, and psychophysical data relating to attention, and it makes a variety of experimentally testable predictions.

Earthquakes of the Nepal Himalaya : towards a physical model of the seismic cycle

Home to hundreds of millions of souls and land of excessiveness, the Himalaya is also the locus of a unique seismicity whose scope and peculiarities still remain to this day somewhat mysterious. Having claimed the lives of kings, or turned ancient timeworn cities into heaps of rubbles and ruins, earthquakes eerily inhabit Nepalese folk tales with the fatalistic message that nothing lasts forever. From a scientific point of view as much as from a human perspective, solving the mysteries of Himalayan seismicity thus represents a challenge of prime importance. Documenting geodetic strain across the Nepal Himalaya with various GPS and leveling data, we show that unlike other subduction zones that exhibit a heterogeneous and patchy coupling pattern along strike, the last hundred kilometers of the Main Himalayan Thrust fault, or MHT, appear to be uniformly locked, devoid of any of the “creeping barriers” that traditionally ward off the propagation of large events. The approximately 20 mm/yr of reckoned convergence across the Himalaya matching previously established estimates of the secular deformation at the front of the arc, the slip accumulated at depth has to somehow elastically propagate all the way to the surface at some point. And yet, neither large events from the past nor currently recorded microseismicity nearly compensate for the massive moment deficit that quietly builds up under the giant mountains. Along with this large unbalanced moment deficit, the uncommonly homogeneous coupling pattern on the MHT raises the question of whether or not the locked portion of the MHT can rupture all at once in a giant earthquake. Univocally answering this question appears contingent on the still elusive estimate of the magnitude of the largest possible earthquake in the Himalaya, and requires tight constraints on local fault properties. What makes the Himalaya enigmatic also makes it the potential source of an incredible wealth of information, and we exploit some of the oddities of Himalayan seismicity in an effort to improve the understanding of earthquake physics and cipher out the properties of the MHT. Thanks to the Himalaya, the Indo-Gangetic plain is deluged each year under a tremendous amount of water during the annual summer monsoon that collects and bears down on the Indian plate enough to pull it away from the Eurasian plate slightly, temporarily relieving a small portion of the stress mounting on the MHT. As the rainwater evaporates in the dry winter season, the plate rebounds and tension is increased back on the fault. Interestingly, the mild waggle of stress induced by the monsoon rains is about the same size as that from solid-Earth tides which gently tug at the planets solid layers, but whereas changes in earthquake frequency correspond with the annually occurring monsoon, there is no such correlation with Earth tides, which oscillate back-and-forth twice a day. We therefore investigate the general response of the creeping and seismogenic parts of MHT to periodic stresses in order to link these observations to physical parameters. First, the response of the creeping part of the MHT is analyzed with a simple spring-and-slider system bearing rate-strengthening rheology, and we show that at the transition with the locked zone, where the friction becomes near velocity neutral, the response of the slip rate may be amplified at some periods, which values are analytically related to the physical parameters of the problem. Such predictions therefore hold the potential of constraining fault properties on the MHT, but still await observational counterparts to be applied, as nothing indicates that the variations of seismicity rate on the locked part of the MHT are the direct expressions of variations of the slip rate on its creeping part, and no variations of the slip rate have been singled out from the GPS measurements to this day. When shifting to the locked seismogenic part of the MHT, spring-and-slider models with rate-weakening rheology are insufficient to explain the contrasted responses of the seismicity to the periodic loads that tides and monsoon both place on the MHT. Instead, we resort to numerical simulations using the Boundary Integral CYCLes of Earthquakes algorithm and examine the response of a 2D finite fault embedded with a rate-weakening patch to harmonic stress perturbations of various periods. We show that such simulations are able to reproduce results consistent with a gradual amplification of sensitivity as the perturbing period get larger, up to a critical period corresponding to the characteristic time of evolution of the seismicity in response to a step-like perturbation of stress. This increase of sensitivity was not reproduced by simple 1D-spring-slider systems, probably because of the complexity of the nucleation process, reproduced only by 2D-fault models. When the nucleation zone is close to its critical unstable size, its growth becomes highly sensitive to any external perturbations and the timings of produced events may therefore find themselves highly affected. A fully analytical framework has yet to be developed and further work is needed to fully describe the behavior of the fault in terms of physical parameters, which will likely provide the keys to deduce constitutive properties of the MHT from seismological observations.

Binding site size limitations of imidazole-pyrrole polyamides for recognition in the minor groove of DNA

The discovery that the three ring polyamide Im-Py-Py-Dp containing imidazole (Im) and pyrrole (Py) carboxamides binds the DNA sequence 5'-(A,T)G(A,T)C(A,T)-3' as an antiparallel dimer offers a new model for the design of ligands for specific recognition of sequences in the minor groove containing both G,C and A,T base pairs. In Chapter 2, experiments are described in which the sequential addition of five N- methylpyrrolecarboxamides to the imidazole-pyrrole polyamide Im-Py-Py-Dp affords a series of six homologous polyamides, Im-(Py)_2-7-Dp, that differ in the size of their binding site, apparent first order binding affinity, and sequence specificity. These results demonstrate that DNA sequences up to nine base pairs in length can be specifically recognized by imidazole-pyrrole polyamides containing three to seven rings by 2:1 polyamide-DNA complex formation in the minor groove. Recognition of a nine base pair site defines the new lower limit of the binding site size that can be recognized by polyamides containing exclusively imidazole and pyrrolecarboxamides. The results of this study should provide useful guidelines for the design of new polyamides that bind longer DNA sites with enhanced affinity and specificity.

In Chapter 3 the design and synthesis of the hairpin polyamide Im-Py-Im-Py-γ-Im- Py-Im-Py-Dp is described. Quantitative DNase I footprint titration experiments reveal that Im-Py-Im-Py-γ-Im-Py-Im-Py-Dp binds six base pair 5'-(A,T)GCGC(A,T)-3' sequences with 30-fold higher affinity than the unlinked polyamide Im-Py-Im-Py-Dp. The hairpin polyamide does not discriminate between A•T and T•A at the first and sixth positions of the binding site as three sites 5'-TGCGCT-3', 5'-TGCGCA-3', and 5 'AGCGCT- 3' are bound with similar affinity. However, Im-Py-Im-Py-γ-Im-Py-Im-PyDp is specific for and discriminates between G•C and C•G base pairs in the 5'-GCGC-3' core as evidenced by lower affinities for the mismatched sites 5'-AACGCA-3', 5'- TGCGTT-3', 5'-TGCGGT-3', and 5'-ACCGCT-3'.

In Chapter 4, experiments are described in which a kinetically stable hexa-aza Schiff base La³⁺ complex is covalently attached to a Tat(49-72) peptide which has been shown to bind the HIV-1 TAR RNA sequence. Although these metallo-peptides cleave TAR site-specifically in the hexanucleotide loop to afford products consistent with hydrolysis, a series of control experiments suggests that the observed cleavage is not caused by a sequence-specifically bound Tat(49-72)-La(L)³⁺ peptide.

Efficient methods for stochastic optimal control

The Hamilton Jacobi Bellman (HJB) equation is central to stochastic optimal control (SOC) theory, yielding the optimal solution to general problems specified by known dynamics and a specified cost functional. Given the assumption of quadratic cost on the control input, it is well known that the HJB reduces to a particular partial differential equation (PDE). While powerful, this reduction is not commonly used as the PDE is of second order, is nonlinear, and examples exist where the problem may not have a solution in a classical sense. Furthermore, each state of the system appears as another dimension of the PDE, giving rise to the curse of dimensionality. Since the number of degrees of freedom required to solve the optimal control problem grows exponentially with dimension, the problem becomes intractable for systems with all but modest dimension.

In the last decade researchers have found that under certain, fairly non-restrictive structural assumptions, the HJB may be transformed into a linear PDE, with an interesting analogue in the discretized domain of Markov Decision Processes (MDP). The work presented in this thesis uses the linearity of this particular form of the HJB PDE to push the computational boundaries of stochastic optimal control.

This is done by crafting together previously disjoint lines of research in computation. The first of these is the use of Sum of Squares (SOS) techniques for synthesis of control policies. A candidate polynomial with variable coefficients is proposed as the solution to the stochastic optimal control problem. An SOS relaxation is then taken to the partial differential constraints, leading to a hierarchy of semidefinite relaxations with improving sub-optimality gap. The resulting approximate solutions are shown to be guaranteed over- and under-approximations for the optimal value function. It is shown that these results extend to arbitrary parabolic and elliptic PDEs, yielding a novel method for Uncertainty Quantification (UQ) of systems governed by partial differential constraints. Domain decomposition techniques are also made available, allowing for such problems to be solved via parallelization and low-order polynomials.

The optimization-based SOS technique is then contrasted with the Separated Representation (SR) approach from the applied mathematics community. The technique allows for systems of equations to be solved through a low-rank decomposition that results in algorithms that scale linearly with dimensionality. Its application in stochastic optimal control allows for previously uncomputable problems to be solved quickly, scaling to such complex systems as the Quadcopter and VTOL aircraft. This technique may be combined with the SOS approach, yielding not only a numerical technique, but also an analytical one that allows for entirely new classes of systems to be studied and for stability properties to be guaranteed.

The analysis of the linear HJB is completed by the study of its implications in application. It is shown that the HJB and a popular technique in robotics, the use of navigation functions, sit on opposite ends of a spectrum of optimization problems, upon which tradeoffs may be made in problem complexity. Analytical solutions to the HJB in these settings are available in simplified domains, yielding guidance towards optimality for approximation schemes. Finally, the use of HJB equations in temporal multi-task planning problems is investigated. It is demonstrated that such problems are reducible to a sequence of SOC problems linked via boundary conditions. The linearity of the PDE allows us to pre-compute control policy primitives and then compose them, at essentially zero cost, to satisfy a complex temporal logic specification.