Results 1  10
of
136
Topologically Adaptable Snakes
 Medical Image Analysis
, 1995
"... This paper presents a topologically adaptable snakes model for image segmentation and object representation. The model is embedded in the framework of domain subdivision using simplicial decomposition. This framework extends the geometric and topological adaptability of snakes while retaining all of ..."
Abstract

Cited by 195 (5 self)
 Add to MetaCart
This paper presents a topologically adaptable snakes model for image segmentation and object representation. The model is embedded in the framework of domain subdivision using simplicial decomposition. This framework extends the geometric and topological adaptability of snakes while retaining all of the features of traditionalsnakes, such as user interaction, and overcoming many of the limitations of traditionalsnakes. By superposing a simplicial grid over the image domain and using this grid to iteratively reparameterize the deforming snakes model, the model is able to flow into complex shapes, even shapes with significant protrusions or branches, and to dynamically change topology as necessitated by the data. Snakes can be created and can split into multiple parts or seamlessly merge into other snakes. The model can also be easily converted to and from the traditional parametric snakes model representation. We apply a 2D model to various synthetic and real images in order to segment ...
SOLVING SYSTEMS OF POLYNOMIAL EQUATIONS
, 2002
"... These are the lecture notes for ten lectures to be given at the CBMS ..."
Abstract

Cited by 150 (12 self)
 Add to MetaCart
These are the lecture notes for ten lectures to be given at the CBMS
Numerical methods for image registration
, 2004
"... In this paper we introduce a new framework for image registration. Our formulation is based on consistent discretization of the optimization problem coupled with a multigrid solution of the linear system which evolve in a GaussNewton iteration. We show that our discretization is helliptic independ ..."
Abstract

Cited by 148 (25 self)
 Add to MetaCart
In this paper we introduce a new framework for image registration. Our formulation is based on consistent discretization of the optimization problem coupled with a multigrid solution of the linear system which evolve in a GaussNewton iteration. We show that our discretization is helliptic independent of parameter choice and therefore a simple multigrid implementation can be used. To overcome potential large nonlinearities and to further speed up computation, we use a multilevel continuation technique. We demonstrate the efficiency of our method on a realistic highly nonlinear registration problem. 1 Introduction and problem setup Image registration is one of today’s challenging image processing problems. Given a socalled reference R and a socalled template image T, the basic idea is to find a “reasonable ” transformation such that a transformed version of the template image becomes “similar ” to the reference image. Image registration
Parallel LagrangeNewtonKrylovSchur methods for PDEconstrained optimization. Part I: The KrylovSchur solver
 SIAM J. Sci. Comput
, 2000
"... Abstract. Large scale optimization of systems governed by partial differential equations (PDEs) is a frontier problem in scientific computation. The stateoftheart for such problems is reduced quasiNewton sequential quadratic programming (SQP) methods. These methods take full advantage of existin ..."
Abstract

Cited by 72 (11 self)
 Add to MetaCart
Abstract. Large scale optimization of systems governed by partial differential equations (PDEs) is a frontier problem in scientific computation. The stateoftheart for such problems is reduced quasiNewton sequential quadratic programming (SQP) methods. These methods take full advantage of existing PDE solver technology and parallelize well. However, their algorithmic scalability is questionable; for certain problem classes they can be very slow to converge. In this twopart article we propose a new method for steadystate PDEconstrained optimization, based on the idea of full space SQP with reduced space quasiNewton SQP preconditioning. The basic components of the method are: Newton solution of the firstorder optimality conditions that characterize stationarity of the Lagrangian function; Krylov solution of the KarushKuhnTucker (KKT) linear systems arising at each Newton iteration using a symmetric quasiminimum residual method; preconditioning of the KKT system using an approximate state/decision variable decomposition that replaces the forward PDE Jacobians by their own preconditioners, and the decision space Schur complement (the reduced Hessian) by a BFGS approximation or by a twostep stationary method. Accordingly, we term the new method LagrangeNewtonKrylov Schur (LNKS). It is fully parallelizable, exploits the structure of available parallel algorithms for the PDE forward problem, and is locally quadratically convergent. In the first part of the paper we investigate the effectiveness of the KKT linear system solver. We test the method on two optimal control problems in which the flow is described by the steadystate Stokes equations. The
Numerical Decomposition of the Solution Sets of Polynomial Systems into Irreducible Components
, 2001
"... In engineering and applied mathematics, polynomial systems arise whose solution sets contain components of different dimensions and multiplicities. In this article we present algorithms, based on homotopy continuation, that compute much of the geometric information contained in the primary decomposi ..."
Abstract

Cited by 56 (26 self)
 Add to MetaCart
In engineering and applied mathematics, polynomial systems arise whose solution sets contain components of different dimensions and multiplicities. In this article we present algorithms, based on homotopy continuation, that compute much of the geometric information contained in the primary decomposition of the solution set. In particular, ignoring multiplicities, our algorithms lay out the decomposition of the set of solutions into irreducible components, by finding, at each dimension, generic points on each component. As byproducts, the computation also determines the degree of each component and an upper bound on itsmultiplicity. The bound issharp (i.e., equal to one) for reduced components. The algorithms make essential use of generic projection and interpolation, and can, if desired, describe each irreducible component precisely as the common zeroesof a finite number of polynomials.
Sparse Elimination and Applications in Kinematics
, 1994
"... This thesis proposes efficient algorithmic solutions to problems in computational algebra and computational algebraic geometry. Moreover, it considers their application to different areas where algebraic systems describe kinematic and geometric constraints. Given an arbitrary system of nonlinear mul ..."
Abstract

Cited by 49 (11 self)
 Add to MetaCart
This thesis proposes efficient algorithmic solutions to problems in computational algebra and computational algebraic geometry. Moreover, it considers their application to different areas where algebraic systems describe kinematic and geometric constraints. Given an arbitrary system of nonlinear multivariate polynomial equations, its resultant serves in eliminating variables and reduces root finding to a linear eigenproblem. Our contribution is to describe the first efficient and general algorithms for computing the sparse resultant. The sparse resultant generalizes the classical homogeneous resultant and exploits the structure of the given polynomials. Its size depends only on the geometry of the input Newton polytopes. The first algorithm uses a subdivision of the Minkowski sum and produces matrix...
PHoM  a Polyhedral Homotopy Continuation Method for Polynomial Systems
 Computing
, 2003
"... PHoM is a software package in C++ for finding all isolated solutions of polynomial systems using a polyhedral homotopy continuation method. Among three modules constituting the package, the first module StartSystem constructs a family of polyhedrallinear homotopy functions, based on the polyhedral ..."
Abstract

Cited by 31 (13 self)
 Add to MetaCart
PHoM is a software package in C++ for finding all isolated solutions of polynomial systems using a polyhedral homotopy continuation method. Among three modules constituting the package, the first module StartSystem constructs a family of polyhedrallinear homotopy functions, based on the polyhedral homotopy theory, from input data for a given system of polynomial equations f (x) = 0. The second module CMPSc traces the solution curves of the homotopy equations to compute all isolated solutions of f (x) = 0. The third module Verify checks whether all isolated solutions of f (x) = 0 have been approximated correctly. We describe numerical methods used in each module and the usage of the package. Numerical results to demonstrate the performance of PHoM include some large polynomial systems that have not been solved previously.