Shape from Contour

The problem of using image contours to infer the shapes and orientations of surfaces is treated as a problem of statistical estimation. The basis for solving this problem lies in an understanding of the geometry of contour formation, coupled with simple statistical models of the contour generating process. This approach is first applied to the special case of surfaces known to be planar. The distortion of contour shape imposed by projection is treated as a signal to be estimated, and variations of non-projective origin are treated as noise. The resulting method is then extended to the estimation of curved surfaces, and applied successfully to natural images. Next, the geometric treatment is further extended by relating countour curvature to surface curvature, using cast shadows as a model for contour generation. This geometric relation, combined with a statistical model, provides a measure of goodness-of-fit between a surface and an image contour. The goodness-of-fit measure is applied to the problem of establishing registration between an image and a surface model. Finally, the statistical estimation strategy is experimentally compared to human perception of orientation: human observers' judgements of tilt correspond closely to the estimates produced by the planar strategy.

Comparison of blinkers and explosive events: A case study

Bewsher, D., Innes, D.E., Parnell, C.E. and Brown, D.S., 2005, Comparison of blinkers and explosive events, Astronomy and Astrophysics, 432, 307. Sponsorship: PPARC

Quantum behavior in mesoscopic systems

In this thesis I present the work done during my PhD. The Thesis is divided into two parts; in the first one I present the study of mesoscopic quantum systems whereas in the second one I address the problem of the definition of Markov regime for quantum system dynamics. The first work presented is the study of vortex patterns in (quasi) two dimensional rotating Bose Einstein condensates (BECs). I consider the case of an anisotropy trapping potential and I shall show that the ground state of the system hosts vortex patterns that are unstable. In a second work I designed an experimental scheme to transfer entanglement from two entangled photons to two BECs. This work is meant to propose a feasible experimental set up to bring entanglement from microscopic to macroscopic systems for both the study of fundamental questions (quantum to classical transition) and technological applications. In the last work of the first part another experimental scheme is presented in order to detect coherences of a mechanical oscillator which is assumed to have been previously cooled down to the quantum regime. In this regime in fact the system can rapidly undergo decoherence so that new techniques have to be employed in order to detect and manipulate their states. In the scheme I propose a micro-mechanical oscillator is coupled to a BEC and the detection is performed by monitoring the BEC with a negligible back-action on the cantilever. In the second part of the thesis I give a definition of Markov regime for open quantum dynamics. The importance of such definition comes from both the mathematical description of the system dynamics and from the understanding of the role played by the environment in the evolution of an open system. In the Markov regime the mathematical description can be simplified and the role of the environment is a passive one.

Functional annotation signatures of disease susceptibility loci improve SNP association analysis.

BACKGROUND: Genetic association studies are conducted to discover genetic loci that contribute to an inherited trait, identify the variants behind these associations and ascertain their functional role in determining the phenotype. To date, functional annotations of the genetic variants have rarely played more than an indirect role in assessing evidence for association. Here, we demonstrate how these data can be systematically integrated into an association study's analysis plan. RESULTS: We developed a Bayesian statistical model for the prior probability of phenotype-genotype association that incorporates data from past association studies and publicly available functional annotation data regarding the susceptibility variants under study. The model takes the form of a binary regression of association status on a set of annotation variables whose coefficients were estimated through an analysis of associated SNPs in the GWAS Catalog (GC). The functional predictors examined included measures that have been demonstrated to correlate with the association status of SNPs in the GC and some whose utility in this regard is speculative: summaries of the UCSC Human Genome Browser ENCODE super-track data, dbSNP function class, sequence conservation summaries, proximity to genomic variants in the Database of Genomic Variants and known regulatory elements in the Open Regulatory Annotation database, PolyPhen-2 probabilities and RegulomeDB categories. Because we expected that only a fraction of the annotations would contribute to predicting association, we employed a penalized likelihood method to reduce the impact of non-informative predictors and evaluated the model's ability to predict GC SNPs not used to construct the model. We show that the functional data alone are predictive of a SNP's presence in the GC. Further, using data from a genome-wide study of ovarian cancer, we demonstrate that their use as prior data when testing for association is practical at the genome-wide scale and improves power to detect associations. CONCLUSIONS: We show how diverse functional annotations can be efficiently combined to create 'functional signatures' that predict the a priori odds of a variant's association to a trait and how these signatures can be integrated into a standard genome-wide-scale association analysis, resulting in improved power to detect truly associated variants.

Fission dynamics for capture reactions in 58,64ni + 208pb systems: New results in terms of thermal energy and neutron multiplicity correlated distributions

The neutron multidetector DéMoN has been used to investigate the symmetric splitting dynamics in the reactions 58.64Ni + 208Pb with excitation energies ranging from 65 to 186 MeV for the composite system. An analysis based on the new backtracing technique has been applied on the neutron data to determine the two-dimensional correlations between the parent composite system initial thermal energy (EthCN) and the total neutron multiplicity (νtot), and between pre- and post-scission neutron multiplicities (νpre and νpost, respectively). The νpre distribution shape indicates the possible coexistence of fast-fission and fusion-fission for the system 58Ni + 208Pb (Ebeam = 8.86 A MeV). The analysis of the neutron multiplicities in the framework of the combined dynamical statistical model (CDSM) gives a reduced friction coefficient β = 23 ± 2512 × 1021 s-1, above the one-body dissipation limit. The corresponding fission time is τf = 40 ± 4620 × 10-21 s. © 1999 Elsevier Science B.V. All rights reserved.

Potential Climate Change Effects on the Habitat of Antarctic Krill in the Weddell Quadrant of the Southern Ocean

Antarctic krill is a cold water species, an increasingly important fishery resource and a major prey item for many fish, birds and mammals in the Southern Ocean. The fishery and the summer foraging sites of many of these predators are concentrated between 0 degrees and 90 degrees W. Parts of this quadrant have experienced recent localised sea surface warming of up to 0.2 degrees C per decade, and projections suggest that further widespread warming of 0.27 degrees to 1.08 degrees C will occur by the late 21st century. We assessed the potential influence of this projected warming on Antarctic krill habitat with a statistical model that links growth to temperature and chlorophyll concentration. The results divide the quadrant into two zones: a band around the Antarctic Circumpolar Current in which habitat quality is particularly vulnerable to warming, and a southern area which is relatively insensitive. Our analysis suggests that the direct effects of warming could reduce the area of growth habitat by up to 20%. The reduction in growth habitat within the range of predators, such as Antarctic fur seals, that forage from breeding sites on South Georgia could be up to 55%, and the habitat's ability to support Antarctic krill biomass production within this range could be reduced by up to 68%. Sensitivity analysis suggests that the effects of a 50% change in summer chlorophyll concentration could be more significant than the direct effects of warming. A reduction in primary production could lead to further habitat degradation but, even if chlorophyll increased by 50%, projected warming would still cause some degradation of the habitat accessible to predators. While there is considerable uncertainty in these projections, they suggest that future climate change could have a significant negative effect on Antarctic krill growth habitat and, consequently, on Southern Ocean biodiversity and ecosystem services.

Risk adjustment for hospital use using social security data: cross sectional small area analysis

Objectives: To identify demographic and socioeconomic determinants of need for acute hospital treatment at small area level. To establish whether there is a relation between poverty and use of inpatient services. To devise a risk adjustment formula for distributing public funds for hospital services using, as far as possible, variables that can be updated between censuses. Design: Cross sectional analysis. Spatial interactive modelling was used to quantify the proximity of the population to health service facilities. Two stage weighted least squares regression was used to model use against supply of hospital and community services and a wide range of potential needs drivers including health, socioeconomic census variables, uptake of income support and family credit, and religious denomination. Setting: Northern Ireland. Main outcome measure: Intensity of use of inpatient services. Results: After endogeneity of supply and use was taken into account, a statistical model was produced that predicted use based on five variables: income support, family credit, elderly people living alone, all ages standardised mortality ratio, and low birth weight. The main effect of the formula produced is to move resources from urban to rural areas. Conclusions: This work has produced a population risk adjustment formula for acute hospital treatment in which four of the five variables can be updated annually rather than relying on census derived data. Inclusion of the social security data makes a substantial difference to the model and to the results produced by the formula.

Development of sampling methods for Raman analysis of solid dosage forms of therapeutic and illicit drugs

The results of a study aimed at determining the most important experimental parameters for automated, quantitative analysis of solid dosage form pharmaceuticals (seized and model 'ecstasy' tablets) are reported. Data obtained with a macro-Raman spectrometer were complemented by micro-Raman measurements, which gave information on particle size and provided excellent data for developing statistical models of the sampling errors associated with collecting data as a series of grid points on the tablets' surface. Spectra recorded at single points on the surface of seized MDMA-caffeine-lactose tablets with a Raman microscope (lambda(ex) = 785 nm, 3 mum diameter spot) were typically dominated by one or other of the three components, consistent with Raman mapping data which showed the drug and caffeine microcrystals were ca 40 mum in diameter. Spectra collected with a microscope from eight points on a 200 mum grid were combined and in the resultant spectra the average value of the Raman band intensity ratio used to quantify the MDMA: caffeine ratio, mu(r), was 1.19 with an unacceptably high standard deviation, sigma(r), of 1.20. In contrast, with a conventional macro-Raman system (150 mum spot diameter), combined eight grid point data gave mu(r) = 1.47 with sigma(r) = 0.16. A simple statistical model which could be used to predict sigma(r) under the various conditions used was developed. The model showed that the decrease in sigma(r) on moving to a 150 mum spot was too large to be due entirely to the increased spot diameter but was consistent with the increased sampling volume that arose from a combination of the larger spot size and depth of focus in the macroscopic system. With the macro-Raman system, combining 64 grid points (0.5 mm spacing and 1-2 s accumulation per point) to give a single averaged spectrum for a tablet was found to be a practical balance between minimizing sampling errors and keeping overhead times at an acceptable level. The effectiveness of this sampling strategy was also tested by quantitative analysis of a set of model ecstasy tablets prepared from MDEA-sorbitol (0-30% by mass MDEA). A simple univariate calibration model of averaged 64 point data had R-2 = 0.998 and an r.m.s. standard error of prediction of 1.1% whereas data obtained by sampling just four points on the same tablet showed deviations from the calibration of up to 5%.

Influence of mix proprtions on rheolgy of cement grouts containing limestone powder

In this paper the parameters of cement grout affecting rheological behaviour and compressive strength are investigated. Factorial experimental design was adopted in this investigation to assess the combined effects of the following factors on fluidity, rheological properties, induced bleeding and compressive strength: water/binder ratio (W/B), dosage of superplasticiser (SP), dosage of viscosity agent (VA), and proportion of limestone powder as replacement of cement (LSP). Mini-slump test, Marsh cone, Lombardi plate cohesion meter, induced bleeding test, coaxial rotating cylinder viscometer were used to evaluate the rheology of the cement grout and the compressive strengths at 7 and 28 days were measured. A two-level fractional factorial statistical model was used to model the influence of key parameters on properties affecting the fluidity, the rheology and compressive strength. The models are valid for mixes with 0.35-0.42 W/B, 0.3-1.2% SP, 0.02-0.7% VA (percentage of binder) and 12-45% LSP as replacement of cement. The influences of W/B, SP, VA and LSP were characterised and analysed using polynomial regression which can identify the primary factors and their interactions on the measured properties. Mathematical polynomials were developed for mini-slump, plate cohesion meter, inducing bleeding, yield value, plastic viscosity and compressive strength as function of W/B, SP, VA and proportion of LSP. The statistical approach used highlighted the limestone powder effect and the dosage of SP and VA on the various rheological characteristics of cement grout

Uranium halide complexes in ionic liquids

We report here the syntheses, characterisation and electrochemistry of some 1-ethyl-3-methylimidazolium, [emim], uranium halide salts. The electrochemistry of the uranium halide salts were investigated in both basic and acidic haloaluminate ionic liquids (ILs). The solid state structures of the uranium chloride salts have previously been reported, but have now been re-evaluted using a new statistical model to determine the presence or absence of weak hydrogen bonding interactions in the crystalline state.

Cold-Atom-Induced Control of an Optomechanical Device

We consider a cavity with a vibrating end mirror and coupled to a Bose-Einstein condensate. The cavity field mediates the interplay between mirror and collective oscillations of the atomic density. We study the implications of this dynamics and the possibility of an indirect diagnostic. Our predictions can be observed in a realistic setup that is central to the current quest for mesoscopic quantumness.

The top-down mechanism for body-mass-abundance scaling

Scaling relationships between mean body masses and abundances of species in multitrophic communities continue to be a subject of intense research and debate. The top-down mechanism explored in this paper explains the frequently observed inverse linear relationship between body mass and abundance (i.e., constant biomass) in terms of a balancing of resource biomasses by behaviorally and evolutionarily adapting foragers, and the evolutionary response of resources to this foraging pressure. The mechanism is tested using an allometric, multitrophic community model with a complex food web structure. It is a statistical model describing the evolutionary and population dynamics of tens to hundreds of species in a uniform way. Particularities of the model are the detailed representation of the evolution and interaction of trophic traits to reproduce topological food web patterns, prey switching behavior modeled after experimental observations, and the evolutionary adaptation of attack rates. Model structure and design are discussed. For model states comparable to natural communities, we find that (1) the body-mass-abundance scaling does not depend on the allometric scaling exponent of physiological rates in the form expected from the energetic equivalence rule or other bottom-up theories; (2) the scaling exponent of abundance as a function of body mass is approximately -1, independent of the allometric exponent for physiological rates assumed; (3) removal of top-down control destroys this pattern, and energetic equivalence is recovered. We conclude that the top-down mechanism is active in the model, and that it is a viable alternative to bottom-up mechanisms for controlling body-mass-abundance relations in natural communities.

Entanglement detection in hybrid optomechanical systems

We study a device formed by a Bose-Einstein condensate (BEC) coupled to the field of a cavity with a moving end mirror and find a working point such that the mirror-light entanglement is reproduced by the BEC-light quantum correlations. This provides an experimentally viable tool for inferring mirror-light entanglement with only a limited set of assumptions. We prove the existence of tripartite entanglement in the hybrid device, persisting up to temperatures of a few milli-Kelvin, and discuss a scheme to detect it.

Probing mechanical quantum coherence with an ultracold-atom meter

We propose a scheme to probe quantum coherence in the state of a nanocantilever based on its magnetic coupling (mediated by a magnetic tip) with a spinor Bose Einstein condensate (BEC). By mapping the BEC into a rotor, its coupling with the cantilever results in a gyroscopic motion whose properties depend on the state of the cantilever: the dynamics of one of the components of the rotor angular momentum turns out to be strictly related to the presence of quantum coherence in the state of the cantilever. We also suggest a detection scheme relying on Faraday rotation, which produces only a very small back-action on the BEC and is thus suitable for a continuous detection of the cantilever's dynamics.

2D silhouette and 3D skeletal models for human detection and tracking

In this paper we propose a statistical model for detection and tracking of human silhouette and the corresponding 3D skeletal structure in gait sequences. We follow a point distribution model (PDM) approach using a Principal Component Analysis (PCA). The problem of non-lineal PCA is partially resolved by applying a different PDM depending of pose estimation; frontal, lateral and diagonal, estimated by Fisher's linear discriminant. Additionally, the fitting is carried out by selecting the closest allowable shape from the training set by means of a nearest neighbor classifier. To improve the performance of the model we develop a human gait analysis to take into account temporal dynamic to track the human body. The incorporation of temporal constraints on the model increase reliability and robustness.