Efficient total variation algorithm for fetal MRI reconstruction

Fetal MRI reconstruction aims at finding a high-resolution image given a small set of low-resolution images. It is usually modeled as an inverse problem where the regularization term plays a central role in the reconstruction quality. Literature has considered several regularization terms s.a. Dirichlet/Laplacian energy [1], Total Variation (TV)based energies [2,3] and more recently non-local means [4]. Although TV energies are quite attractive because of their ability in edge preservation, standard explicit steepest gradient techniques have been applied to optimize fetal-based TV energies. The main contribution of this work lies in the introduction of a well-posed TV algorithm from the point of view of convex optimization. Specifically, our proposed TV optimization algorithm for fetal reconstruction is optimal w.r.t. the asymptotic and iterative convergence speeds O(1/n(2)) and O(1/root epsilon), while existing techniques are in O(1/n) and O(1/epsilon). We apply our algorithm to (1) clinical newborn data, considered as ground truth, and (2) clinical fetal acquisitions. Our algorithm compares favorably with the literature in terms of speed and accuracy.

Sex ratio bias, male aggression, and population collapse in lizards.

The adult sex ratio (ASR) is a key parameter of the demography of human and other animal populations, yet the causes of variation in ASR, how individuals respond to this variation, and how their response feeds back into population dynamics remain poorly understood. A prevalent hypothesis is that ASR is regulated by intrasexual competition, which would cause more mortality or emigration in the sex of increasing frequency. Our experimental manipulation of populations of the common lizard (Lacerta vivipara) shows the opposite effect. Male mortality and emigration are not higher under male-biased ASR. Rather, an excess of adult males begets aggression toward adult females, whose survival and fecundity drop, along with their emigration rate. The ensuing prediction that adult male skew should be amplified and total population size should decline is supported by long-term data. Numerical projections show that this amplifying effect causes a major risk of population extinction. In general, such an "evolutionary trap" toward extinction threatens populations in which there is a substantial mating cost for females, and environmental changes or management practices skew the ASR toward males.

Operational sex ratio, sexual conflict and the intensity of sexual selection.

Modern sexual selection theory indicates that reproductive costs rather than the operational sex ratio predict the intensity of sexual selection. We investigated sexual selection in the polygynandrous common lizard Lacerta vivipara. This species shows male aggression, causing high mating costs for females when adult sex ratios (ASR) are male-biased. We manipulated ASR in 12 experimental populations and quantified the intensity of sexual selection based on the relationship between reproductive success and body size. In sharp contrast to classical sexual selection theory predictions, positive directional sexual selection on male size was stronger and positive directional selection on female size weaker in female-biased populations than in male-biased populations. Thus, consistent with modern theory, directional sexual selection on male size was weaker in populations with higher female mating costs. This suggests that the costs of breeding, but not the operational sex ratio, correctly predicted the strength of sexual selection.

Question Waves: A multicast query routing algorithm for social search

This paper describes Question Waves, an algorithm that can be applied to social search protocols, such as Asknext or Sixearch. In this model, the queries are propagated through the social network, with faster propagation through more trustable acquaintances. Question Waves uses local information to make decisions and obtain an answer ranking. With Question Waves, the answers that arrive first are the most likely to be relevant, and we computed the correlation of answer relevance with the order of arrival to demonstrate this result. We obtained correlations equivalent to the heuristics that use global knowledge, such as profile similarity among users or the expertise value of an agent. Because Question Waves is compatible with the social search protocol Asknext, it is possible to stop a search when enough relevant answers have been found; additionally, stopping the search early only introduces a minimal risk of not obtaining the best possible answer. Furthermore, Question Waves does not require a re-ranking algorithm because the results arrive sorted

Female common lizards (Lacerta vivipara) do not adjust their sex-biased investment in relation to the adult sex ratio.

Sex allocation theory predicts that facultative maternal investment in the rare sex should be favoured by natural selection when breeders experience predictable variation in adult sex ratios (ASRs). We found significant spatial and predictable interannual changes in local ASRs within a natural population of the common lizard where the mean ASR is female-biased, thus validating the key assumptions of adaptive sex ratio models. We tested for facultative maternal investment in the rare sex during and after an experimental perturbation of the ASR by creating populations with female-biased or male-biased ASR. Mothers did not adjust their clutch sex ratio during or after the ASR perturbation, but produced sons with a higher body condition in male-biased populations. However, this differential sex allocation did not result in growth or survival differences in offspring. Our results thus contradict the predictions of adaptive models and challenge the idea that facultative investment in the rare sex might be a mechanism regulating the population sex ratio.

Dismissal of the illusion of uncertainty in the assessment of a likelihood ratio

The use of the Bayes factor (BF) or likelihood ratio as a metric to assess the probative value of forensic traces is largely supported by operational standards and recommendations in different forensic disciplines. However, the progress towards more widespread consensus about foundational principles is still fragile as it raises new problems about which views differ. It is not uncommon e.g. to encounter scientists who feel the need to compute the probability distribution of a given expression of evidential value (i.e. a BF), or to place intervals or significance probabilities on such a quantity. The article here presents arguments to show that such views involve a misconception of principles and abuse of language. The conclusion of the discussion is that, in a given case at hand, forensic scientists ought to offer to a court of justice a given single value for the BF, rather than an expression based on a distribution over a range of values.

Unbreakable PT symmetry of solitons supported by inhomogeneous defocusing nonlinearity

Peer Reviewed

Parental body condition does not correlate with offspring sex-ratio in Cory's shearwater

We analyzed offspring sex ratio variation in Mediterranean Cory's Shearwater (Calonectris d. diomedea) during two consecutive breeding seasons in two colonies. We test for differential breeding conditions between years and colonies looking at several breeding parameters and parental condition. We then explored the relationship between offspring sex ratio and parental condition and breeding parameters. This species is sexually dimorphic with males larger and heavier than females; consequently we expected differential parental cost in rearing sexes, or a greater sensitivity of male chicks to adverse conditions, which may lead to biased sex ratios. Chicks were sexed molecularly by the amplification of the CHD genes. Offspring sex ratio did not differ from parity, either at hatching or fledging, regardless of the colony or year. However, parental body condition and breeding parameters such as egg size and breeding success were different between years and colonies. Nevertheless, neither nestling mortality nor body condition at fledging varied between years or colonies, suggesting that male and female chicks were probably not differentially affected by variability in breeding conditions.

Pose-based SLAM with probabilistic scan matching algorithm using a mechanical scanned imaging sonar

This paper proposes a pose-based algorithm to solve the full SLAM problem for an autonomous underwater vehicle (AUV), navigating in an unknown and possibly unstructured environment. The technique incorporate probabilistic scan matching with range scans gathered from a mechanical scanning imaging sonar (MSIS) and the robot dead-reckoning displacements estimated from a Doppler velocity log (DVL) and a motion reference unit (MRU). The proposed method utilizes two extended Kalman filters (EKF). The first, estimates the local path travelled by the robot while grabbing the scan as well as its uncertainty and provides position estimates for correcting the distortions that the vehicle motion produces in the acoustic images. The second is an augment state EKF that estimates and keeps the registered scans poses. The raw data from the sensors are processed and fused in-line. No priory structural information or initial pose are considered. The algorithm has been tested on an AUV guided along a 600 m path within a marina environment, showing the viability of the proposed approach

A concurrent region growing algorithm guided by circumscribed contours

Image segmentation of natural scenes constitutes a major problem in machine vision. This paper presents a new proposal for the image segmentation problem which has been based on the integration of edge and region information. This approach begins by detecting the main contours of the scene which are later used to guide a concurrent set of growing processes. A previous analysis of the seed pixels permits adjustment of the homogeneity criterion to the region's characteristics during the growing process. Since the high variability of regions representing outdoor scenes makes the classical homogeneity criteria useless, a new homogeneity criterion based on clustering analysis and convex hull construction is proposed. Experimental results have proven the reliability of the proposed approach

Design and evaluation of learning algorithms for dynamic resource management in virtual networks

Network virtualisation is considerably gaining attentionas a solution to ossification of the Internet. However, thesuccess of network virtualisation will depend in part on how efficientlythe virtual networks utilise substrate network resources.In this paper, we propose a machine learning-based approachto virtual network resource management. We propose to modelthe substrate network as a decentralised system and introducea learning algorithm in each substrate node and substrate link,providing self-organization capabilities. We propose a multiagentlearning algorithm that carries out the substrate network resourcemanagement in a coordinated and decentralised way. The taskof these agents is to use evaluative feedback to learn an optimalpolicy so as to dynamically allocate network resources to virtualnodes and links. The agents ensure that while the virtual networkshave the resources they need at any given time, only the requiredresources are reserved for this purpose. Simulations show thatour dynamic approach significantly improves the virtual networkacceptance ratio and the maximum number of accepted virtualnetwork requests at any time while ensuring that virtual networkquality of service requirements such as packet drop rate andvirtual link delay are not affected.

Studies of Eu- and Yb-induced Reconstructions on the Si(100) Surface

This thesis presents experimental studies of rare earth (RE) metal induced structures on Si(100) surfaces. Two divalent RE metal adsorbates, Eu and Yb, are investigated on nominally flat Si(100) and on vicinal, stepped Si(100) substrates. Several experimental methods have been applied, including scanning tunneling microscopy/spectroscopy (STM/STS), low energy electron diffraction (LEED), synchrotron radiation photoelectron spectroscopy (SR-PES), Auger electron spectroscopy (AES), thermal desorption spectroscopy (TDS), and work function change measurements (Δφ). Two stages can be distinguished in the initial growth of the RE/Si interface: the formation of a two-dimensional (2D) adsorbed layer at submonolayer coverage and the growth of a three-dimensional (3D) silicide phase at higher coverage. The 2D phase is studied for both adsorbates in order to discover whether they produce common reconstructions or reconstructions common to the other RE metals. For studies of the 3D phase Yb is chosen due to its ability to crystallize in a hexagonal AlB2 type lattice, which is the structure of RE silicide nanowires, therefore allowing for the possibility of the growth of one-dimensional (1D) wires. It is found that despite their similar electronic configuration, Eu and Yb do not form similar 2D reconstructions on Si(100). Instead, a wealth of 2D structures is observed and atomic models are proposed for the 2×3-type reconstructions. In addition, adsorbate induced modifications on surface morphology and orientational symmetry are observed. The formation of the Yb silicide phase follows the Stranski-Krastanov growth mode. Nanowires with the hexagonal lattice are observed on the flat Si(100) substrate, and moreover, an unexpectedly large variety of growth directions are revealed. On the vicinal substrate the growth of the silicide phase as 3D islands and wires depends drastically on the growth conditions. The conditions under which wires with high aspect ratio and single orientation parallel to the step edges can be formed are demonstrated.

Fragmentation of massive dense cores down to <~ 1000 AU: relation between fragmentation and density structure

In order to shed light on the main physical processes controlling fragmentation of massive dense cores, we present a uniform study of the density structure of 19 massive dense cores, selected to be at similar evolutionary stages, for which their relative fragmentation level was assessed in a previous work. We inferred the density structure of the 19 cores through a simultaneous fit of the radial intensity profiles at 450 and 850 μm (or 1.2 mm in two cases) and the spectral energy distribution, assuming spherical symmetry and that the density and temperature of the cores decrease with radius following power-laws. Even though the estimated fragmentation level is strictly speaking a lower limit, its relative value is significant and several trends could be explored with our data. We find a weak (inverse) trend of fragmentation level and density power-law index, with steeper density profiles tending to show lower fragmentation, and vice versa. In addition, we find a trend of fragmentation increasing with density within a given radius, which arises from a combination of flat density profile and high central density and is consistent with Jeans fragmentation. We considered the effects of rotational-to-gravitational energy ratio, non-thermal velocity dispersion, and turbulence mode on the density structure of the cores, and found that compressive turbulence seems to yield higher central densities. Finally, a possible explanation for the origin of cores with concentrated density profiles, which are the cores showing no fragmentation, could be related with a strong magnetic field, consistent with the outcome of radiation magnetohydrodynamic simulations.

New Distance Transforms for Gray Level Image Compression

This thesis deals with distance transforms which are a fundamental issue in image processing and computer vision. In this thesis, two new distance transforms for gray level images are presented. As a new application for distance transforms, they are applied to gray level image compression. The new distance transforms are both new extensions of the well known distance transform algorithm developed by Rosenfeld, Pfaltz and Lay. With some modification their algorithm which calculates a distance transform on binary images with a chosen kernel has been made to calculate a chessboard like distance transform with integer numbers (DTOCS) and a real value distance transform (EDTOCS) on gray level images. Both distance transforms, the DTOCS and EDTOCS, require only two passes over the graylevel image and are extremely simple to implement. Only two image buffers are needed: The original gray level image and the binary image which defines the region(s) of calculation. No other image buffers are needed even if more than one iteration round is performed. For large neighborhoods and complicated images the two pass distance algorithm has to be applied to the image more than once, typically 3 10 times. Different types of kernels can be adopted. It is important to notice that no other existing transform calculates the same kind of distance map as the DTOCS. All the other gray weighted distance function, GRAYMAT etc. algorithms find the minimum path joining two points by the smallest sum of gray levels or weighting the distance values directly by the gray levels in some manner. The DTOCS does not weight them that way. The DTOCS gives a weighted version of the chessboard distance map. The weights are not constant, but gray value differences of the original image. The difference between the DTOCS map and other distance transforms for gray level images is shown. The difference between the DTOCS and EDTOCS is that the EDTOCS calculates these gray level differences in a different way. It propagates local Euclidean distances inside a kernel. Analytical derivations of some results concerning the DTOCS and the EDTOCS are presented. Commonly distance transforms are used for feature extraction in pattern recognition and learning. Their use in image compression is very rare. This thesis introduces a new application area for distance transforms. Three new image compression algorithms based on the DTOCS and one based on the EDTOCS are presented. Control points, i.e. points that are considered fundamental for the reconstruction of the image, are selected from the gray level image using the DTOCS and the EDTOCS. The first group of methods select the maximas of the distance image to new control points and the second group of methods compare the DTOCS distance to binary image chessboard distance. The effect of applying threshold masks of different sizes along the threshold boundaries is studied. The time complexity of the compression algorithms is analyzed both analytically and experimentally. It is shown that the time complexity of the algorithms is independent of the number of control points, i.e. the compression ratio. Also a new morphological image decompression scheme is presented, the 8 kernels' method. Several decompressed images are presented. The best results are obtained using the Delaunay triangulation. The obtained image quality equals that of the DCT images with a 4 x 4