information geometry on the curved q exponential family with application to survival data analysis abstract information geometry and statistical manifold have attracted wide attention in the past few decades the amari chentsov structure plays a central role both in optimization and statistical inference because of the invariance under sufficient statistics this paper discusses the amari chentsov structure on statistical manifolds induced by the curved q exponential family with type i censored data the reliability function and cumulant generating function of the family are obtained where the cumulant generating function contains a random variable depends on the experimental design after the construction of the statistical manifold the geometric quantities i e fisher metric amari tensor and connection coefficients are investigated we find that the affine connection is not equal to 0 which is different from the existing results the relationship among the geometric quantities based on the q exponential family and curved q exponential family is discussed the relationship between the curvature on the manifold in the normal and tangent directions is obtained employing the maximum likelihood and bayesian methods the parameter vector and its dual are estimated the main results are illustrated with the two parameter q exponential distribution and a real data analysis based on a terminal human cancer data set

unsupervised variational image segmentation classification using a weibull observation model studies have shown that the weibull distribution can model accurately a wide variety of images its parameters index a family of distributions which includes the exponential and approximations of the gaussian and the raleigh models widely used in image segmentation this study investigates the weibull distribution in unsupervised image segmentation and classification by a variational method the data term of the segmentation functional measures the conformity of the image intensity in each region to a weibull distribution whose parameters are determined jointly with the segmentation minimization of the functional is implemented by active curves via level sets and consists of iterations of two consecutive steps curve evolution via euler lagrange descent equations and evaluation of the weibull distribution parameters experiments with synthetic and real images are described which verify the validity of method and its implementation

active contours without edges we propose a new model for active contours to detect objects in a given image based on techniques of curve evolution mumford shah 1989 functional for segmentation and level sets our model can detect objects whose boundaries are not necessarily defined by the gradient we minimize an energy which can be seen as a particular case of the minimal partition problem in the level set formulation the problem becomes a mean curvature flow like evolving the active contour which will stop on the desired boundary however the stopping term does not depend on the gradient of the image as in the classical active contour models but is instead related to a particular segmentation of the image we give a numerical algorithm using finite differences finally we present various experimental results and in particular some examples for which the classical snakes methods based on the gradient are not applicable also the initial curve can be anywhere in the image and interior contours are automatically detected

sobolev active contours all previous geometric active contour models that have been formulated as gradient flows of various energies use the same l 2 type inner product to define the notion of gradient recent work has shown that this inner product induces a pathological riemannian metric on the space of smooth curves however there are also undesirable features associated with the gradient flows that this inner product induces in this paper we reformulate the generic geometric active contour model by redefining the notion of gradient in accordance with sobolev type inner products we call the resulting flows sobolev active contours sobolev metrics induce favorable regularity properties in their gradient flows in addition sobolev active contours favor global translations but are not restricted to such motions they are also less susceptible to certain types of local minima in contrast to traditional active contours these properties are particularly useful in tracking applications we demonstrate the general methodology by reformulating some standard edge based and region based active contour models as sobolev active contours and show the substantial improvements gained in segmentation

integrability and the painleve property for low dimensional systems we examine a system described by two first order nonlinear differential equations from the point of view of integrability the singularity analysis in the complex time plane is used to investigate the painleve property which according to the ablowitz ramani segur conjecture is a prerequisite for integrability for infinite dimensional systems we show that for such low dimensional systems the painleve analysis is still a most useful guide but integrable cases also exist which do not possess the painleve property

on the non integrability of a family of duffing van der pol oscillators we investigate the non integrability of a family of duffing van der pol oscillators x alpha x x2 1 x beta x3 gamma cos omega t by studying the analytic properties of the dynamics in complex time we find that the solutions of have no worse than algebraic singularities at t with only t t 1 2 terms present in their series expansions unlike for example the alpha 0 duffing case where typically log t t terms arise still when integrating around long enough contours a remarkably intricate pattern of square root singularities emerges on different sheets which appears to prevent solutions from ever returning to the original sheet such evidence of infinitely sheeted solutions termed the iss property has also been observed in a number of hamiltonian systems and is illustrated here on a simple example of a single first order differential equation we suggest that the iss property is a necessary condition for non integrability i e non existence of a complete set of analytic single valued constants of the motion which would permit the complete integration of a dynamical system in terms of quadratures

localizing region based active contours in this paper we propose a natural framework that allows any region based segmentation energy to be re formulated in a local way we consider local rather than global image statistics and evolve a contour based on local information localized contours are capable of segmenting objects with heterogeneous feature profiles that would be difficult to capture correctly using a standard global method the presented technique is versatile enough to be used with any global region based active contour energy and instill in it the benefits of localization we describe this framework and demonstrate the localization of three well known energies in order to illustrate how our framework can be applied to any energy we then compare each localized energy to its global counterpart to show the improvements that can be achieved next an in depth study of the behaviors of these energies in response to the degree of localization is given finally we show results on challenging images to illustrate the robust and accurate segmentations that are possible with this new class of active contour models

exploiting information geometry to improve the convergence of nonparametric active contours this paper presents a fast converging riemannian steepest descent method for nonparametric statistical active contour models with application to image segmentation unlike other fast algorithms the proposed method is general and can be applied to any statistical active contour model from the exponential family which comprises most of the models considered in the literature this is achieved by first identifying the intrinsic statistical manifold associated with this class of active contours and then constructing a steepest descent on that manifold a key contribution of this paper is to derive a general and tractable closed form analytic expression for the manifold s riemannian metric tensor which allows computing discrete gradient flows efficiently the proposed methodology is demonstrated empirically and compared with other state of the art approaches on several standard test images a phantom positron emission tomography scan and a b mode echography of in vivo human dermis

on the analytic structure of the driven pendulum the analytic structure of the solution of the driven pendulum is investigated through painleve analysis in the complex time plane the existence is pointed out of a two armed infinite sheeted riemann structure of the singularities after an exponential transformation

singularity clustering in the duffing oscillator an asymptotic and numerical study is made of the singularity structure in the complex t plane of the duffing oscillator the presence of logarithmic terms in the local psi series expansion of the form t4 ln t leads to a multisheeted singularity structure of great complexity this structure is built recursively from an elemental pattern which takes the form of four armed stars of singularities this construction is deduced analytically from the properties of the mapping z t4 ln t and is confirmed quite accurately numerically a systematic resummation of the psi series in terms of lame functions is developed this series exhibits the same analytic structure at all orders and provides a semi local analytical representation of the solution which is apparently valid even in the chaotic regime

the analytic structure of dynamical systems and self similar natural boundaries abstract in this paper we investigate the analytic complex time structure of the movable singularities for several dynamical systems in general it is found that there exists a direct connection between the occurencce of a certain type of multiple valuedness of the singularities and the existence of a class of remarkable self similar natural boundaries for these systems an asymptotic description of the distribution of singularities in the natural boundary is developed this provides a description of the fine scale structure of these natural boundaries that agrees closely with the numerical calculations

on the singularity analysis of intersecting separatrices in near integrable dynamical systems abstract it has been recently proved by s l ziglin that transversal intersections of separatrices invariant manifolds in near integrable hamiltonian systems of two degrees of freedom imply the existence of multi valued solutions with infinitely many riemann sheets ziglin s theorem is illustrated here on a simple example and then extended and applied to non hamiltonian analytic flows x f x eg x t with x x y and g x t g x t 2 which for e 0 possess the painleve property on the other hand the theoretically expected logarithmic singularities for e 0 are obtained explicitly in solutions near the intersecting separatrices thus we conjecture that dynamical systems with the painleve property can have no separatrix intersections and hence no strange attractors etc these singularities are then numerically located and found to form certain very interesting chimney patterns in the complex t plane on which they accumulate densely the upper part of the chimneys away from the re t axis is asymptotically quite insensitive to changes in parameter values or initial conditions the singularity pattern itself however becomes denser and each chimney is seen to gradually collapse towards the re t axis as the amplitude of the driving term increases and the motion becomes more chaotic

integrating the nonintegrable analytic structure of the lorenz system revisited abstract a study of the complex time analytic structure of the lorenz system in nonintegrable parameter regimes reveals the special sets of parameter values for which one time dependent integral of motion exists furthermore the analysis yields the exact form of the part of the integral with the highest homogeneous weight and a method to construct the rest of the integral recursive clustering of singularities in the chaotic regimes of the system is observed in computer studies and explained by a simple analytic argument the analytic techniques used in these studies a systematic resummation of a logarithmic psi series appears to be quite general and can provide explicit representations of a solution even in the chaotic regimes in the neighborhood of a given movable singularity furthermore we suggest that this technique provides a type of renormalization program to study a wide class of nonintegrable systems

image segmentation using active contours calculus of variations or shape gradients we consider the problem of segmenting an image through the minimization of an energy criterion involving region and boundary functionals we show that one can go from one class to the otherby solving poisson s orhelmholtz s equation with well chosen boundar y conditions using this equivalence we study the case of a large class of region functionals by standard methods of the calculus of variations and derive the corresponding euler lagrange equations we revisit this problem using the notion of a shape derivative and show that the same equations can be elegantly derived without going through the unnatural step of converting the region integrals into boundary integrals we also define a larger class of region functionals based on the estimation and comparison to a prototype of the probability density distribution of image features and show how the shape derivative tool allows us to easily compute the corresponding gateaux derivatives and euler lagrange equations finally we apply this new functional to the problem of regions segmentation in sequences of color images we briefly describe our numerical scheme and show some experimental results

estimating the granularity coefficient of a potts markov random field within a markov chain monte carlo algorithm this paper addresses the problem of estimating the potts parameter jointly with the unknown parameters of a bayesian model within a markov chain monte carlo mcmc algorithm standard mcmc methods cannot be applied to this problem because performing inference on requires computing the intractable normalizing constant of the potts model in the proposed mcmc method the estimation of is conducted using a likelihood free metropolis hastings algorithm experimental results obtained for synthetic data show that estimating jointly with the other unknown parameters leads to estimation results that are as good as those obtained with the actual value of on the other hand choosing an incorrect value of can degrade estimation performance significantly to illustrate the interest of this method the proposed algorithm is successfully applied to real bidimensional sar and tridimensional ultrasound images

generalized gradients priors on minimization flows this paper tackles an important aspect of the variational problem underlying active contours optimization by gradient flows classically the definition of a gradient depends directly on the choice of an inner product structure this consideration is largely absent from the active contours literature most authors explicitely or implicitely assume that the space of admissible deformations is ruled by the canonical l 2 inner product the classical gradient flows reported in the literature are relative to this particular choice here we investigate the relevance of using i other inner products yielding other gradient descents and ii other minimizing flows not deriving from any inner product in particular we show how to induce different degrees of spatial consistency into the minimizing flow in order to decrease the probability of getting trapped into irrelevant local minima we report numerical experiments indicating that the sensitivity of the active contours method to initial conditions which seriously limits its applicability and efficiency is alleviated by our application specific spatially coherent minimizing flows we show that the choice of the inner product can be seen as a prior on the deformation fields and we present an extension of the definition of the gradient toward more general priors

ergodicity of billiard flows and quadratic differentials on considere un systeme billard de 2 objets de masses m 1 et m 2 on montre que pour un ensemble dense de paires m 1 m 2 ce systeme est ergodique

image segmentation with a sobolev gradient method abstract the most effective methods for finding object boundaries in a digital image involve minimizing a functional over a set of curves or surfaces where the functional includes internal energy terms for regularization and external energy terms that align the curves or surfaces with object boundaries current practice is to seek critical points of the energy functional by what amounts to a steepest descent iteration with the discretized l 2 gradient since the functional involves derivatives a descent method with a discretized sobolev gradient is likely to be much more efficient we demonstrate this with test results for an implementation of a variational level set method for edge based segmentation with active contours in two dimensions

closed trajectories for quadratic differentials with an application to billiards on considere les orbites periodiques pour le systeme dynamique d une table polygonale avec des angles multiples rationnels de le billard rationnel

on the complete and partial integrability of non hamiltonian systems abstract the methods of singularity analysis are applied to several third order non hamiltonian systems of physical significance including the lotka volterra equations the three wave interaction and the rikitake dynamo model complete integrability is defined and new completely integrable systems are discovered by means of the painleve property in all these cases we obtain integrals which reduce the equations either to a final quadrature or to an irreducible second order ordinary differential equation ode solved by painleve transcendents relaxing the painleve property we find many partially integrable cases whose movable singularities are poles at leading order with in t t 0 terms entering at higher orders in an n th order generalized rossler model a precise relation is established between the partial fulfillment of the painleve conditions and the existence of n 2 integrals of the motion

new possibilities with sobolev active contours recently the sobolev metric was introduced to define gradient flows of various geometric active contour energies it was shown that the sobolev metric outperforms the traditional metric for the same energy in many cases such as for tracking where the coarse scale changes of the contour are important some interesting properties of sobolev gradient flows include that they stabilize certain unstable traditional flows and the order of the evolution pdes are reduced when compared with traditional gradient flows of the same energies in this paper we explore new possibilities for active contours made possible by sobolev metrics the sobolev method allows one to implement new energy based active contour models that were not otherwise considered because the traditional minimizing method render them ill posed or numerically infeasible in particular we exploit the stabilizing and the order reducing properties of sobolev gradients to implement the gradient descent of these new energies we give examples of this class of energies which include some simple geometric priors and new edge based energies we also show that these energies can be quite useful for segmentation and tracking we also show that the gradient flows using the traditional metric are either ill posed or numerically difficult to implement and then show that the flows can be implemented in a stable and numerically feasible manner using the sobolev gradient

asymptotic analysis of time singularities for a class of time dependent hamiltonians we analyse the local and global structure of time singularities for a class of quasi integrable hamiltonian systems in the arnold liouville sense we show that there is good agreement between the numerically observed local behaviour of the solutions and the perturbative scheme we produce using asymptotic approximations of the solution around the singularities we also prove the convergence of the psi series associated to the movable singularities of the systems considered we also propose a simple model in order to analyse the global structure of the singularities in the directions of exponential growth of the potential in time

a second order shape optimization approach for image segmentation the problem of segmentation of a given image using the active contour technique is considered an abstract calculus to find appropriate speed functions for active contour models in image segmentation or related problems based on variational principles is presented the speed method from shape sensitivity analysis is used to derive speed functions which correspond to gradient or newton type directions for the underlying optimization problem the newton type speed function is found by solving an elliptic problem on the current active contour in every time step numerical experiments comparing the classical gradient method with newton s method are presented

teichmuller curves in moduli space eisenstein series and an application to triangular billiards there exists a teichmuller disc n containing the riemann surface ofy2 x n 1 in the genus n 1 2 teichmuller space such that the stabilizer of n in the mapping class group has a fundamental domain of finite poincare volume in n application is given to an asymptotic formula for the length spectrum of the billiard in isosceles triangles with angles n n n 2 n and to the uniform distribution of infinite billiard trajectories in the same triangles