Curve segmentation using directional information, relation to pattern detection

©2005 IEEE. Personal use of this material is permitted. However, permission to reprint/republish this material for advertising or promotional purposes or for creating new collective works for resale or distribution to servers or lists, or to reuse any copyrighted component of this work in other works must be obtained from the IEEE. This material is presented to ensure timely dissemination of scholarly and technical work. Copyright and all rights therein are retained by authors or by other copyright holders. All persons copying this information are expected to adhere to the terms and constraints invoked by each author's copyright. In most cases, these works may not be reposted without the explicit permission of the copyright holder.Presented at the 2005 International Conference on Image Processing (ICIP)September 11-14, 2005, Genova, Italy.DOI: 10.1109/ICIP.2005.1530175We propose an extension of the conformal (or geodesic) active contour framework in which the conformal factor depends not only on the position of the curve but also on the direction of its tangent. We describe several properties for variational curve segmentation schemes that justify the construction of optimal conformal factors (i.e., learning) in strong connection with pattern matching. The determination of optimal curves (i.e., segmentation) can be performed using either the calculus of variations or dynamic programming. The technique is illustrated on a road detection problem for different signal to noise ratios

Finsler Active Contours

©2008 IEEE. Personal use of this material is permitted. However, permission to reprint/republish this material for advertising or promotional purposes or for creating new collective works for resale or distribution to servers or lists, or to reuse any copyrighted component of this work in other works must be obtained from the IEEE. This material is presented to ensure timely dissemination of scholarly and technical work. Copyright and all rights therein are retained by authors or by other copyright holders. All persons copying this information are expected to adhere to the terms and constraints invoked by each author's copyright. In most cases, these works may not be reposted without the explicit permission of the copyright holder.DOI: 10.1109/TPAMI.2007.70713In this paper, we propose an image segmentation technique based on augmenting the conformal (or geodesic) active contour framework with directional information. In the isotropic case, the euclidean metric is locally multiplied by a scalar conformal factor based on image information such that the weighted length of curves lying on points of interest (typically edges) is small. The conformal factor that is chosen depends only upon position and is in this sense isotropic. Although directional information has been studied previously for other segmentation frameworks, here, we show that if one desires to add directionality in the conformal active contour framework, then one gets a well-defined minimization problem in the case that the factor defines a Finsler metric. Optimal curves may be obtained using the calculus of variations or dynamic programming-based schemes. Finally, we demonstrate the technique by extracting roads from aerial imagery, blood vessels from medical angiograms, and neural tracts from diffusion-weighted magnetic resonance imagery

On the filamentary environment of galaxies

The correlation between the large-scale distribution of galaxies and their spectroscopic properties at z=1.5 is investigated using the Horizon MareNostrum cosmological run. We have extracted a large sample of 10^5 galaxies from this large hydrodynamical simulation featuring standard galaxy formation physics. Spectral synthesis is applied to these single stellar populations to generate spectra and colours for all galaxies. We use the skeleton as a tracer of the cosmic web and study how our galaxy catalogue depends on the distance to the skeleton. We show that galaxies closer to the skeleton tend to be redder, but that the effect is mostly due to the proximity of large haloes at the nodes of the skeleton, rather than the filaments themselves. This effects translate into a bimodality in the colour distribution of our sample. The origin of this bimodality is investigated and seems to follow from the ram pressure stripping of satellite galaxies within the more massive clusters of the simulation. The virtual catalogues (spectroscopical properties of the MareNostrum galaxies at various redshifts) are available online at http://www.iap.fr/users/pichon/MareNostrum/cataloguesComment: 18 pages, 27 figures, accepted for publication in MNRA

Initial Conditions for Large Cosmological Simulations

This technical paper describes a software package that was designed to produce initial conditions for large cosmological simulations in the context of the Horizon collaboration. These tools generalize E. Bertschinger's Grafic1 software to distributed parallel architectures and offer a flexible alternative to the Grafic2 software for ``zoom'' initial conditions, at the price of large cumulated cpu and memory usage. The codes have been validated up to resolutions of 4096^3 and were used to generate the initial conditions of large hydrodynamical and dark matter simulations. They also provide means to generate constrained realisations for the purpose of generating initial conditions compatible with, e.g. the local group, or the SDSS catalog.Comment: 12 pages, 11 figures, submitted to ApJ

Non Gaussian extrema counts for CMB maps

In the context of the geometrical analysis of weakly non Gaussian CMB maps, the 2D differential extrema counts as functions of the excursion set threshold is derived from the full moments expansion of the joint probability distribution of an isotropic random field, its gradient and invariants of the Hessian. Analytic expressions for these counts are given to second order in the non Gaussian correction, while a Monte Carlo method to compute them to arbitrary order is presented. Matching count statistics to these estimators is illustrated on fiducial non-Gaussian "Planck" data.Comment: 4 pages, 1 figur

Massive spheroids can form in single minor mergers

Accepted for publication in MNRAS, 12 pages, 6 figuresUnderstanding how rotationally supported discs transform into dispersion-dominated spheroids is central to our comprehension of galaxy evolution. Morphological transformation is largely merger-driven. While major mergers can efficiently create spheroids, recent work has highlighted the significant role of other processes, like minor mergers, in driving morphological change. Given their rich merger histories, spheroids typically exhibit large fractions of ‘ex situ’ stellar mass, i.e. mass that is accreted, via mergers, from external objects. This is particularly true for the most massive galaxies, whose stellar masses typically cannot be attained without a large number of mergers. Here, we explore an unusual population of extremely massive (M ∗ > 10 11M) spheroids, in the Horizon-AGN simulation, which exhibit anomalously low ex situ mass fractions, indicating that they form without recourse to significant merging. These systems form in a single minor-merger event (with typical merger mass ratios of 0.11–0.33), with a specific orbital configuration, where the satellite orbit is virtually co-planar with the disc of the massive galaxy. The merger triggers a catastrophic change in morphology, over only a few hundred Myr, coupled with strong in situ star formation. While this channel produces a minority (∼5 per cent) of such galaxies, our study demonstrates that the formation of at least some of the most massive spheroids need not involve major mergers – or any significant merging at all – contrary to what is classically believed.Peer reviewedFinal Accepted Versio

Generation of Vorticity and Velocity Dispersion by Orbit Crossing

We study the generation of vorticity and velocity dispersion by orbit crossing using cosmological numerical simulations, and calculate the backreaction of these effects on the evolution of large-scale density and velocity divergence power spectra. We use Delaunay tessellations to define the velocity field, showing that the power spectra of velocity divergence and vorticity measured in this way are unbiased and have better noise properties than for standard interpolation methods that deal with mass weighted velocities. We show that high resolution simulations are required to recover the correct large-scale vorticity power spectrum, while poor resolution can spuriously amplify its amplitude by more than one order of magnitude. We measure the scalar and vector modes of the stress tensor induced by orbit crossing using an adaptive technique, showing that its vector modes lead, when input into the vorticity evolution equation, to the same vorticity power spectrum obtained from the Delaunay method. We incorporate orbit crossing corrections to the evolution of large scale density and velocity fields in perturbation theory by using the measured stress tensor modes. We find that at large scales (k~0.1 h/Mpc) vector modes have very little effect in the density power spectrum, while scalar modes (velocity dispersion) can induce percent level corrections at z=0, particularly in the velocity divergence power spectrum. In addition, we show that the velocity power spectrum is smaller than predicted by linear theory until well into the nonlinear regime, with little contribution from virial velocities.Comment: 27 pages, 14 figures. v2: reorganization of the material, new appendix. Accepted by PR

Crustal structure and the mechanism of extension in the northern Red Sea: Constraints from gravity anomalies

Gravity data from detailed marine surveys were used to investigate the density structure of the crust and upper mantle beneath the northern Red Sea, a continental rift considered by many to be in the process of establishing a mid-ocean ridge type spreading center. A gravity and bathymetry profile across the sea near 27°N was analyzed in detail. Model constraints include single-channel and multichannel seismic reflection lines coincident with the gravity profile, which provide detailed information on sedimentary layers, a series of ESP seismic experiments centered nearly on the gravity line which constrain the sediment and crustal thickness, and a series of 58 closely spaced heat flow measurements nearly coincident with the gravity line which closely constrain the mantle thermal structure and thus density variations within the mantle. The calculated gravity anomalies based on the known crustal and mantle structure fail to match the observed data. Specifically, the model anomalies become steadily more negative from the margins toward the center of the sea with the result that the calculated gravity is approximately 25 mGal more negative near the axis relative to the marginal areas than are the observed anomalies. Two possible solutions were investigated. One is to introduce a lower crustal layer resulting from crustal underplating within what was originally considered to be the upper mantle beneath the marginal areas. The other is to introduce an area of slightly denser crust resulting from intrusion of new mafic crustal material within the axial depression. Both solutions can reproduce the observed gravity anomalies within the Red Sea. However, the complete set of geophysical data, and in particular the heat flow in the Red Sea and gravity and seismic refraction data on land, suggests that the most satisfactory model is that the crust of the axial depression is denser than that of the marginal areas. This implies that crustal extension in the northern Red Sea is presently occurring in large part by intrusion of new mafic material