The paper is devoted to an approach for image inpainting developed on the basis of neurogeometry of vision and sub-Riemannian geometry. Inpainting is realized by completing damaged isophotes (level lines of brightness) by optimal curves for the left-invariant sub-Riemannian problem on the group of rototranslations (motions) of a plane SE(2). The approach is considered as anthropomorphic inpainting since these curves satisfy the variational principle discovered by neurogeometry of vision. A parallel algorithm and software to restore monochrome binary or halftone images represented as series of isophotes were developed. The approach and the algorithm for computation of completing arcs are presented in detail.

In this paper, an extremal eigenvalue problem to the Sturm-Liouville equations with discontinuous coefficients and volume constraint is investigated. Liouville transformation is applied to change the problem into an equivalent minimization problem. Finite element method is proposed and the convergence for the finite element solution is established. A monotonic decreasing algorithm is presented to solve the extremal eigenvalue problem. A global convergence for the algorithm in the continuous case is proved. A few numerical results are given to depict the efficiency of the method.

Aiming at the isoparametric bilinear finite volume element scheme, we initially derive an asymptotic expansion and a high accuracy combination formula of the derivatives in the sense of pointwise by employing the energy-embedded method on uniform grids. Furthermore, we prove that the approximate derivatives are convergent of order two. Finally, numerical examples verify the theoretical results.

The z-transform is introduced to analyze a full discretization method for a partial integro-differential equation (PIDE) with a weakly singular kernel. In this method, spectral collocation is used for the spatial discretization, and, for the time stepping, the finite difference method combined with the convolution quadrature rule is considered. The global stability and convergence properties of complete discretization are derived and numerical experiments are reported.

This article describes the implementation of a simple wavelet-based optical-flow motion estimator dedicated to continuous motions such as fluid flows. The wavelet representation of the unknown velocity field is considered. This scale-space representation, associated to a simple gradient-based optimization algorithm, sets up a well-defined multiresolution framework for the optical flow estimation. Moreover, a very simple closure mechanism, approaching locally the solution by high-order polynomials is provided by truncating the wavelet basis at fine scales. Accuracy and efficiency of the proposed method is evaluated on image sequences of turbulent fluid flows.

This paper is concerned with obtaining an approximate solution and an approximate derivative of the solution for neutral Volterra integro-differential equation with a weakly singular kernel. The solution of this equation, even for analytic data, is not smooth on the entire interval of integration. The Jacobi collocation discretization is proposed for the given equation. A rigorous analysis of error bound is also provided which theoretically justifies that both the error of approximate solution and the error of approximate derivative of the solution decay exponentially in L∞ norm and weighted L2 norm. Numerical results are presented to demonstrate the effectiveness of the spectral method.

We consider the linear and non-linear enhancement of diffusion weighted magnetic resonance images (DW-MRI) to use contextual information in denoising and inferring fiber crossings. We describe the space of DW-MRI images in a moving frame of reference, attached to fiber fragments which allows for convection-diffusion along the fibers. Because of this approach, our method is naturally able to handle crossings in data. We will perform experiments showing the ability of the enhancement to infer information about crossing structures, even in diffusion tensor images (DTI) which are incapable of representing crossings themselves. We will present a novel non-linear enhancement technique which performs better than linear methods in areas around ventricles, thereby eliminating the need for additional preprocessing steps to segment out the ventricles. We pay special attention to the details of implementation of the various numeric schemes.

One of the classical optimization models for image segmentation is the well known Markov Random Fields (MRF) model. This model is a discrete optimization problem, which is shown here to formulate many continuous models used in image segmentation. In spite of the presence of MRF in the literature, the dominant perception has been that the model is not effective for image segmentation. We show here that the reason for the non-effectiveness is due to the lack of access to the optimal solution. Instead of solving optimally, heuristics have been engaged. Those heuristic methods cannot guarantee the quality of the solution nor the running time of the algorithm. Worse still, heuristics do not link directly the input functions and parameters to the output thus obscuring what would be ideal choices of parameters and functions which are to be selected by users in each particular application context.

We describe here how MRF can model and solve efficiently several known continuous models for image segmentation and describe briefly a very efficient polynomial time algorithm, which is provably fastest possible, to solve optimally the MRF problem. The MRF algorithm is enhanced here compared to the algorithm in Hochbaum (2001) by allowing the set of assigned labels to be any discrete set. Other enhancements include dynamic features that permit adjustments to the input parameters and solves optimally for these changes with minimal computation time. Several new theoretical results on the properties of the algorithm are proved here and are demonstrated for images in the context of medical and biological imaging. An interactive implementation tool for MRF is described, and its performance and flexibility in practice are demonstrated via computational experiments.

We conclude that many continuous models common in image segmentation have discrete analogs to various special cases of MRF and as such are solved optimally and efficiently, rather than with the use of continuous techniques, such as PDE methods, that restrict the type of functions used and furthermore, can only guarantee convergence to a local minimum.

In this paper, we explore the use of the diffusion geometry framework for the fusion of geometric and photometric information in local and global shape descriptors. Our construction is based on the definition of a diffusion process on the shape manifold embedded into a high-dimensional space where the embedding coordinates represent the photometric information. Experimental results show that such data fusion is useful in coping with different challenges of shape analysis where pure geometric and pure photometric methods fail.

Partial similarity of shapes is a challenging problem arising in many important applications in computer vision, shape analysis, and graphics, e.g. When one has to deal with partial information and acquisition artifacts. The problem is especially hard when the underlying shapes are non-rigid and are given up to a deformation. Partial matching is usually approached by computing local descriptors on a pair of shapes and then establishing a point-wise non-bijective correspondence between the two, taking into account possibly different parts. In this paper, we introduce an alternative correspondence-less approach to matching fragments to an entire shape undergoing a non-rigid deformation. We use region-wise local descriptors and optimize over the integration domains on which the integral descriptors of the two parts match. The problem is regularized using the Mumford-Shah functional. We show an efficient discretization based on the Ambrosio-Tortorelli approximation generalized to triangular point clouds and meshes, and present experiments demonstrating the success of the proposed method.

Detecting similarity between non-rigid shapes is one of the fundamental problems in computer vision. In order to measure the similarity the shapes must first be aligned. As opposite to rigid alignment that can be parameterized using a small number of unknowns representing rotations, reflections and translations, non-rigid alignment is not easily parameterized. Majority of the methods addressing this problem boil down to a minimization of a certain distortion measure. The complexity of a matching process is exponential by nature, but it can be heuristically reduced to a quadratic or even linear for shapes which are smooth two-manifolds. Here we model the shapes using both local and global structures, employ these to construct a quadratic dissimilarity measure, and provide a hierarchical framework for minimizing it to obtain sparse set of corresponding points. These correspondences may serve as an initialization for dense linear correspondence search.

Image deconvolution problems with a symmetric point-spread function arise in many areas of science and engineering. These problems often are solved by the Richardson-Lucy method, a nonlinear iterative method. We first show a convergence result for the Richardson-Lucy method. The proof sheds light on why the method may converge slowly. Subsequently, we describe an iterative active set method that imposes the same constraints on the computed solution as the Richardson-Lucy method. Computed examples show the latter method to yield better restorations than the Richardson-Lucy method and typically require less computational effort.

Total variation regularization has good performance in noise removal and edge preservation but lacks in texture restoration. Here we present a texture-preserving strategy to restore images contaminated by blur and noise. According to a texture detection strategy, we apply spatially adaptive fractional order diffusion. A fast algorithm based on the half-quadratic technique is used to minimize the resulting objective function. Numerical results show the effectiveness of our strategy.

In this paper, we present a surface reconstruction via 2D strokes and a vector field on the strokes based on a two-step method. In the first step, from sparse strokes drawn by artists and a given vector field on the strokes, we propose a nonlinear vector interpolation combining total variation (TV) and H1 regularization with a curl-free constraint for obtaining a dense vector field. In the second step, a height map is obtained by integrating the dense vector field in the first step. Jump discontinuities in surface and discontinuities of surface gradients can be well reconstructed without any surface distortion. We also provide a fast and efficient algorithm for solving the proposed functionals. Since vectors on the strokes are interpreted as a projection of surface gradients onto the plane, different types of strokes are easily devised to generate geometrically crucial structures such as ridge, valley, jump, bump, and dip on the surface. The stroke types help users to create a surface which they intuitively imagine from 2D strokes. We compare our results with conventional methods via many examples.