Eigenvectors and eigenvalues of discrete Laplacians are often used for manifold learning and nonlinear dimensionality reduction. Graph Laplacian is one widely used discrete laplacian on point cloud. It was previously proved by Belkin and Niyogi [4] that the eigenvectors and eigenvalues of the graph Laplacian converge to the eigenfunctions and eigenvalues of the Laplace-Beltrami operator of the manifold in the limit of innitely many data points sampled independently from the uniform distribution over the manifold. Recently, we introduced Point Integral method (PIM) [15] to solve elliptic equations and corresponding eigenvalue problem on point clouds. In this paper, we prove that the eigenvectors and eigenvalues obtained by PIM converge in the limit of innitely many random samples. Moreover, one estimate of the rate of the convergence is given.

Image segmentation is an important task in the domain of computer vision and medical imaging. In natural and medical images, intensity inhomogeneity, i.e. the varying image intensity, occurs often and it poses considerable challenges for image segmentation. In this paper, we propose an efficient variational method for segmenting images with intensity inhomogeneity. The method is inspired by previous works on two-stage segmentation and variational Retinex. Our method consists of two stages. In the first stage, we decouple the image into reflection and illumination parts by solving a convex energy minimization model with either total variation or tight-frame regularisation. In the second stage, we segment the original image by thresholding on the reflection part, and the inhomogeneous intensity is estimated by the smoothly varying illumination part. We adopt a primal dual algorithm to solve the convex model in the first stage, and the convergence is guaranteed. Numerical experiments clearly show that our method is robust and efficient to segment both natural and medical images.

We use a positive parameter to develop a dierentiable perturbed reconstruction model to solve the L1–L2 magnetic resonance image (MRI) reconstruction problem. We use Bregman iterative formulation to solve the dierentiable perturbed L1–L2 model, and lagged diusivity xed-point iteration to solve the minimization problem in the Bregman iteration. Two Fourier transforms and an inverse Fourier transform are used to accelerate L1–L2 MRI reconstruction. Real MR images are used to test the method in numerical experiments. The results demonstrate that the proposed method is very ecient for L1–L2MRI reconstruction.

No abstract uploaded!

In this paper we study the stochastic Galerkin approximation for the linear transport equation with random inputs and diffusive scal- ing. We first establish uniform (in the Knudsen number) stability results in the random space for the transport equation with uncertain scattering coefficients, and then prove the uniform spectral conver- gence (and consequently the sharp stochastic Asymptotic-Preserving property) of the stochastic Galerkin method. A micro-macro decom- position based fully discrete scheme is adopted for the problem and proved to have a uniform stability. Numerical experiments are conducted to demonstrate the stability and asymptotic properties of the method.