Using the Split Bregman Algorithm to Solve the Self-repelling Snakes Model

Preserving contour topology during image segmentation is useful in many practical scenarios. By keeping the contours isomorphic, it is possible to prevent over-segmentation and under-segmentation, as well as to adhere to given topologies. The Self-repelling Snakes model (SR) is a variational model that preserves contour topology by combining a non-local repulsion term with the geodesic active contour model. The SR is traditionally solved using the additive operator splitting (AOS) scheme. In our paper, we propose an alternative solution to the SR using the Split Bregman method. Our algorithm breaks the problem down into simpler sub-problems to use lower-order evolution equations and a simple projection scheme rather than re-initialization. The sub-problems can be solved via fast Fourier transform or an approximate soft thresholding formula which maintains stability, shortening the convergence time, and reduces the memory requirement. The Split Bregman and AOS algorithms are compared theoretically and experimentally.

Efficient Iterative Regularization Method for Total Variation-Based Image Restoration

Total variation (TV) regularization has received much attention in image restoration applications because of its advantages in denoising and preserving details. A common approach to address TV-based image restoration is to design a specific algorithm for solving typical cost function, which consists of conventional ℓ2 fidelity term and TV regularization. In this work, a novel objective function and an efficient algorithm are proposed. Firstly, a pseudoinverse transform-based fidelity term is imposed on TV regularization, and a closely-related optimization problem is established. Then, the split Bregman framework is used to decouple the complex inverse problem into subproblems to reduce computational complexity. Finally, numerical experiments show that the proposed method can obtain satisfactory restoration results with fewer iterations. Combined with the restoration effect and efficiency, this method is superior to the competitive algorithm. Significantly, the proposed method has the advantage of a simple solving structure, which can be easily extended to other image processing applications.

A Two-Staged Feature Extraction Method Based on Total Variation for Hyperspectral Images

Effective feature extraction (FE) has always been the focus of hyperspectral images (HSIs). For aerial remote-sensing HSIs processing and its land cover classification, in this article, an efficient two-staged hyperspectral FE method based on total variation (TV) is proposed. In the first stage, the average fusion method was used to reduce the spectral dimension. Then, the anisotropic TV model with different regularization parameters was utilized to obtain featured blocks of different smoothness, each containing multi-scale structure information, and we stacked them as the next stage’s input. In the second stage, equipped with singular value transformation to reduce the dimension again, we followed an isotropic TV model based on split Bregman algorithm for further detail smoothing. Finally, the feature-extracted block was fed to the support vector machine for classification experiments. The results, with three hyperspectral datasets, demonstrate that our proposed method can competitively outperform state-of-the-art methods in terms of its classification accuracy and computing time. Also, our proposed method delivers robustness and stability by comprehensive parameter analysis.

Tensor Dictionary Learning with an Enhanced Sparsity Constraint for Sparse-View Spectral CT Reconstruction

Spectral computed tomography (CT) can divide collected photons into multi-energy channels and gain multi-channel projections synchronously by using photon-counting detectors. However, reconstructed images usually contain severe noise due to the limited number of photons in the corresponding energy channel. Tensor dictionary learning (TDL)-based methods have achieved better performance, but usually lose image edge information and details, especially from an under-sampling dataset. To address this problem, this paper proposes a method termed TDL with an enhanced sparsity constraint for spectral CT reconstruction. The proposed algorithm inherits the superiority of TDL by exploring the correlation of spectral CT images. Moreover, the method designs a regularization using the L0-norm of the image gradient to constrain images and the difference between images and a prior image in each energy channel simultaneously, further improving the ability to preserve edge information and subtle image details. The split-Bregman algorithm has been applied to address the proposed objective minimization model. Several numerical simulations and realistic preclinical mice are studied to assess the effectiveness of the proposed algorithm. The results demonstrate that the proposed method improves the quality of spectral CT images in terms of noise elimination, edge preservation, and image detail recovery compared to the several existing better methods.

Fast Terahertz Imaging Model Based on Group Sparsity and Nonlocal Self-Similarity

In order to solve the problems of long-term image acquisition time and massive data processing in a terahertz time domain spectroscopy imaging system, a novel fast terahertz imaging model, combined with group sparsity and nonlocal self-similarity (GSNS), is proposed in this paper. In GSNS, the structure similarity and sparsity of image patches in both two-dimensional and three-dimensional space are utilized to obtain high-quality terahertz images. It has the advantages of detail clarity and edge preservation. Furthermore, to overcome the high computational costs of matrix inversion in traditional split Bregman iteration, an acceleration scheme based on conjugate gradient method is proposed to solve the terahertz imaging model more efficiently. Experiments results demonstrate that the proposed approach can lead to better terahertz image reconstruction performance at low sampling rates.

Strip-Map SAR Image Formulation Based on the Modified Alternating Split Bregman Method

Conventional compressive sensing (CS)-based imaging methods allow images to be reconstructed from a small amount of data, while they suffer from high computational burden even for a moderate scene. To address this problem, this paper presents a novel two-dimensional (2D) CS imaging algorithm for strip-map synthetic aperture radars (SARs) with zero squint angle. By introducing a 2D separable formulation to model the physical procedure of the SAR imaging, we separate the large measurement matrix into two small ones, and then the induced algorithm can deal with 2D signal directly instead of converting it into 1D vector. As a result, the computational load can be reduced significantly. Furthermore, thanks to its superior performance in maintaining contour information, the gradient space of the SAR image is exploited and the total variation (TV) constraint is incorporated to improve resolution performance. Due to the non-differentiable property of the TV regularizer, it is difficult to directly solve the induced TV regularization problem. To overcome this problem, an improved split Bregman method is presented by formulating the TV minimization problem into a sequence of unconstrained optimization problem and Bregman updates. It yields an accurate and simple solution. Finally, the synthesis and real experiment results demonstrate that the proposed algorithm remains competitive in terms of high resolution and high computational efficiency.

Hybrid Sparsity Model for Fast Terahertz Imaging

In order to shorten the long-term image acquisition time of the terahertz time domain spectroscopy imaging system while ensuring the imaging quality, a hybrid sparsity model (HSM) is proposed for fast terahertz imaging in this paper, which incorporates both intrinsic sparsity prior and nonlocal self-similarity constraints in a unified statistical model. In HSM, a weighted exponentiation shift-invariant wavelet transform is introduced to enhance the sparsity of the terahertz image. Simultaneously, the nonlocal self-similarity by means of the three-dimensional sparsity in the transform domain is exploited to ensure high-quality terahertz image reconstruction. Finally, a new split Bregman-based iteration algorithm is developed to solve the terahertz imaging model more efficiently. Experiments are presented to verify the effectiveness of the proposed approach.