基于PICCS的非局部TGV能谱CT重建算法

能谱CT将宽谱划分为窄谱,导致通道内光子数目明显减少,加大了噪声影响,故从噪声投影中重建出高质量图像是能谱CT的一个研究热点.传统全变分(total variational,TV)容易造成重建图像中出现块状伪影等问题,总广义全变分(total generalized variation,TGV)算法可以逼近任意阶函数,再结合非局部均值算法的思想,同时考虑到不同能谱通道下重建图像的相关性,将高质量全能谱重建图像作为先验图像指导能谱CT重建,提出了基于先验图像约束压缩感知(prior image constrained compressed sensing,PICCS)的非局部TGV重建算法.实验结果表明,所提算法在抑制噪声的同时能够有效复原图像细节及边缘信息,且收敛速度快.     

A variational proximal alternating linearized minimization in a given metric for limited-angle CT image reconstruction

Due to the restriction of computed tomography (CT) scanning environment, the acquired projection data may be incomplete for exact CT reconstruction. Though some convex optimization methods, such as total variation minimization based method, can be used for incomplete data reconstruction, the edge of reconstruction image may be partly distorted for limited-angle CT reconstruction. To promote the quality of reconstruction image for limited-angle CT imaging, in this paper, a nonconvex and nonsmooth optimization model was investigated. To solve the model, a variational proximal alternating linearized minimization (VPALM) method based on proximal mapping in a given metric was proposed. The proposed method can avoid computing the inverse of a huge system matrix thus can be used to deal with the larger-scale inverse problems. What’s more, we show that each bounded sequence generated by VPALM globally converges to a critical point based on the Kurdyka–Lojasiewicz property. Real data experiments are used to demonstrate the viability and effectiveness of VPALM method, and the results show that the proposed method outperforms two classical CT reconstruction methods.     

Image restoration with a high-order total variation minimization method

In this paper, we propose a fast and efficient way to restore blurred and noisy images with a high-order total variation minimization technique. The proposed method is based on an alternating technique for image deblurring and denoising. It starts by finding an approximate image using a Tikhonov regularization method. This corresponds to a deblurring process with possible artifacts and noise remaining. In the denoising step, a high-order total variation algorithm is used to remove noise in the deblurred image. We see that the edges in the restored image can be preserved quite well and the staircase effect is reduced effectively in the proposed algorithm. We also discuss the convergence of the proposed regularization method. Some numerical results show that the proposed method gives restored images of higher quality than some existing total variation restoration methods such as the fast TV method and the modified TV method with the lagged diffusivity fixed-point iteration.     

Image decoding optimization based on compressive sensing

Transform-based image codec follows the basic principle: the reconstructed quality is decided by the quantization level. Compressive sensing (CS) breaks the limit and states that sparse signals can be perfectly recovered from incomplete or even corrupted information by solving convex optimization. Under the same acquisition of images, if images are represented sparsely enough, they can be reconstructed more accurately by CS recovery than inverse transform. So, in this paper, we utilize a modified TV operator to enhance image sparse representation and reconstruction accuracy, and we acquire image information from transform coefficients corrupted by quantization noise. We can reconstruct the images by CS recovery instead of inverse transform. A CS-based JPEG decoding scheme is obtained and experimental results demonstrate that the proposed methods significantly improve the PSNR and visual quality of reconstructed images compared with original JPEG decoder.     

Image decoding optimization based on compressive sensing

Transform-based image codec follows the basic principle: the reconstructed quality is decided by the quantization level. Compressive sensing (CS) breaks the limit and states that sparse signals can be perfectly recovered from incomplete or even corrupted information by solving convex optimization. Under the same acquisition of images, if images are represented sparsely enough, they can be reconstructed more accurately by CS recovery than inverse transform. So, in this paper, we utilize a modified TV operator to enhance image sparse representation and reconstruction accuracy, and we acquire image information from transform coefficients corrupted by quantization noise. We can reconstruct the images by CS recovery instead of inverse transform. A CS-based JPEG decoding scheme is obtained and experimental results demonstrate that the proposed methods significantly improve the PSNR and visual quality of reconstructed images compared with original JPEG decoder.     

Speckle noise removal in ultrasound images by first- and second-order total variation

Speckle noise contamination is a common issue in ultrasound imaging system. Due to the edge-preserving feature, total variation (TV) regularization-based techniques have been extensively utilized for speckle noise removal. However, TV regularization sometimes causes staircase artifacts as it favors solutions that are piecewise constant. In this paper, we propose a new model to overcome this deficiency. In this model, the regularization term is represented by a combination of total variation and high-order total variation, while the data fidelity term is depicted by a generalized Kullback-Leibler divergence. The proposed model can be efficiently solved by alternating direction method with multipliers (ADMM). Compared with some state-of-the-art methods, the proposed method achieves higher quality in terms of the peak signal to noise ratio (PSNR) and the structural similarity index (SSIM). Numerical experiments demonstrate that our method can remove speckle noise efficiently while suppress staircase effects on both synthetic images and real ultrasound images.     

Image reconstruction model for the exterior problem of computed tomography based on weighted directional total variation

The exterior problem of computed tomography (CT) is a special imaging modality and it is an important research issue in Non-Destructive Testing (NDT). In the exterior problem, the missing projection data vary with the position of the pixels. Furthermore, amounts of theoretical and experimental results showed that the edges that are tangent to the radial direction are much harder to be reconstructed than edges that are tangent to the angular direction. In this paper, a weighted directional total variation (WDTV) based regularization model was proposed to better deal with the exterior problem. By introducing the WDTV regularization term, our model can reconstruct high quality images. First, WDTV of an image also describes the sparsity of image gradient magnitude which can enforce the reconstructed result to be a nearly flat image. Second, our model can preserve edges and reduce artifacts near edges that are tangent to the radial direction. Third, a numerical implementation called SART+WDTV algorithm was developed to solve our model. Simulated and real experiments demonstrated that our model was more capable of suppressing artifacts and preserving edges.     

High-order total variation-based Poissonian image deconvolution with spatially adapted regularization parameter

Images captured by image acquisition systems using photon-counting devices such as astronomical imaging, positron emission tomography and confocal microscopy imaging, are often contaminated by Poisson noise. Total variation (TV) regularization, which is a classic regularization technique in image restoration, is well-known for recovering sharp edges of an image. Since the regularization parameter is important for a good recovery, Chen and Cheng (2012) proposed an effective TV-based Poissonian image deblurring model with a spatially adapted regularization parameter. However, it has drawbacks since the TV regularization produces staircase artifacts. In this paper, in order to remedy the shortcoming of TV of their model, we introduce an extra high-order total variation (HTV) regularization term. Furthermore, to balance the trade-off between edges and the smooth regions in the images, we also incorporate a weighting parameter to discriminate the TV and the HTV penalty. The proposed model is solved by an iterative algorithm under the framework of the well-known alternating direction method of multipliers. Our numerical results demonstrate the effectiveness and efficiency of the proposed method, in terms of signal-to-noise ratio (SNR) and relative error (RelRrr).     

基于分数阶微积分正则化的图像处理

全变分正则化方法已被广泛地应用于图像处理,利用此方法可以较好地去除噪声,并保持图像的边缘特征,但得到的优化解会产生"阶梯"效应.为了克服这一缺点,本文通过分数阶微积分正则化方法,建立了一个新的图像处理模型.为了克服此模型中非光滑项对求解带来的困难,本文研究了基于不动点方程的迫近梯度算法.最后,本文利用提出的模型与算法进行了图像去噪、图像去模糊与图像超分辨率实验,实验结果表明分数阶微积分正则化方法能较好的保留图像纹理等细节信息.     

Total image constrained diffusion tensor for spectral computed tomography reconstruction

Photon counting detector (PCD)-based spectral computed tomography (CT) is a promising imaging technique that enables high energy resolution imaging with narrow energy bins. However, the image quality is degraded because the number of photons in each energy bin is less than the number of photons in the full spectrum. To reconstruct high quality spectral CT images with narrow energy bins, we developed a total image constrained diffusion tensor (TICDT) for statistical iterative reconstruction (SIR) based on a penalized weighted least-squares (PWLS) principle, which is called “PWLS-TICDT.” Specifically, TICDT uses supplementary information from a high-quality total image as a structural prior for SIR, so that the narrow energy bin image can be enhanced, while some primary features are preserved. We also developed an alternating minimization algorithm to solve the associated objective function. We conducted qualitative and quantitative studies to validate and evaluate the PWLS-TICDT method using digital phantoms and preclinical data. Results from both numerical simulation and real PCD data studies show that the proposed PWLS-TICDT method achieves noticeable gains over competing methods in terms of suppressing noise, detecting low contrast objects, and preserving resolution. More importantly, the multi-energy images reconstructed by PWLS-TICDT method can generate more accurate basis material decomposition results than the other methods.     

Low-dose spectral CT reconstruction using image gradient ℓ0–norm and tensor dictionary

Spectral computed tomography (CT) has a great superiority in lesion detection, tissue characterization and material decomposition. To further extend its potential clinical applications, in this work, we propose an improved tensor dictionary learning method for low-dose spectral CT reconstruction with a constraint of image gradient ℓ₀-norm, which is named as ℓ₀TDL. The ℓ₀TDL method inherits the advantages of tensor dictionary learning (TDL) by employing the similarity of spectral CT images. On the other hand, by introducing the ℓ₀-norm constraint in gradient image domain, the proposed method emphasizes the spatial sparsity to overcome the weakness of TDL on preserving edge information. The split-bregman method is employed to solve the proposed method. Both numerical simulations and real mouse studies are perform to evaluate the proposed method. The results show that the proposed ℓ₀TDL method outperforms other competing methods, such as total variation (TV) minimization, TV with low rank (TV+LR), and TDL methods.     

Anisotropic Total Variation Filtering

Total variation regularization and anisotropic filtering have been established as standard methods for image denoising because of their ability to detect and keep prominent edges in the data. Both methods, however, introduce artifacts: In the case of anisotropic filtering, the preservation of edges comes at the cost of the creation of additional structures out of noise; total variation regularization, on the other hand, suffers from the stair-casing effect, which leads to gradual contrast changes in homogeneous objects, especially near curved edges and corners. In order to circumvent these drawbacks, we propose to combine the two regularization techniques. To that end we replace the isotropic TV semi-norm by an anisotropic term that mirrors the directional structure of either the noisy original data or the smoothed image. We provide a detailed existence theory for our regularization method by using the concept of relaxation. The numerical examples concluding the paper show that the proposed introduction of an anisotropy to TV regularization indeed leads to improved denoising: the stair-casing effect is reduced while at the same time the creation of artifacts is suppressed.     

An ADMM algorithm for second-order TV-based MR image reconstruction

In this paper, we propose a new model for MR image reconstruction based on second order total variation ( \(\text {TV}^{2}\) ) regularization and wavelet, which can be considered as requiring the image to be sparse in both the spatial finite differences and wavelet transforms. Furthermore, by applying the variable splitting technique twice, augmented Lagrangian method and the Barzilai-Borwein step size selection scheme, an ADMM algorithm is designed to solve the proposed model. It reduces the reconstruction problem to several unconstrained minimization subproblems, which can be solved by shrinking operators and alternating minimization algorithms. The proposed algorithm needs not to solve a fourth-order PDE but to solve several second-order PDEs so as to improve calculation efficiency. Numerical results demonstrate the effectiveness of the presented algorithm and illustrate that the proposed model outperforms some reconstruction models in the quality of reconstructed images.     

基于字典学习方法的CT不完全投影图像重建算法

对于不完全投影角度的重建研究是CT图像重建中一个重要的问题.将压缩感知中字典学习的方法与CT重建算法ART迭代算法相结合.字典学习方法中字典更新采用K-SVD(K-奇异值分解)算法,稀疏编码采用OMP(正交匹配追踪)算法.最后通过对标准Head头部模型进行仿真实验,验证了字典学习方法在CT图像重建中对于提高图像的重建质量和提高信噪比的可行性与有效性.另外还研究了字典学习中图像块大小和滑动距离对重建图像的影响     

A multi-scale image reconstruction algorithm for electrical capacitance tomography

Electrical capacitance tomography (ECT) is considered as a promising process tomography (PT) technology, and its successful applications depend mainly on the precision and speed of the image reconstruction algorithms. In this paper, based on the wavelet multi-scale analysis method, an efficient image reconstruction algorithm is presented. The original inverse problem is decomposed into a sequence of inverse problems, which are solved successively from the largest scale to the smallest scale. At different scales, the inverse problem is solved by a generalized regularized total least squares (TLS) method, which is developed using a combinational minimax estimation method and an extended stabilizing functional, until the solution of the original inverse problem is found. The homotopy algorithm is employed to solve the objective functional. The proposed algorithm is tested by the noise-free capacitance data and the noise-contaminated capacitance data, and excellent numerical performances and satisfactory results are observed. In the cases considered in this paper, the reconstruction results show remarkable improvement in the accuracy. The spatial resolution of the reconstructed images by the proposed algorithm is enhanced and the artifacts in the reconstructed images can be eliminated effectively. As a result, a promising algorithm is introduced for ECT image reconstruction.     

Efficient Threshold Selection for Multivariate Total Variation Denoising

Total variation (TV) denoising is still attracting attention with theoretical and computational motivations, for its conceptual simplicity of solving a lasso-like convex problem and its good properties for preserving sharp edges and contours in objects with spatial structures like natural images, although more modern and recent techniques specifically tailored to image processing have been developed. TV induces variation-sparsity in the sense that the reconstruction is piecewise constant with a small number of jumps. A threshold parameter λ controls the number of jumps and the quality of the estimation. Since calculation of the TV estimate in high dimension is computationally intensive for a given λ, we propose to calculate the TV estimate for only two sequential λ’s. Our adaptive procedure is based on large deviation of stochastic processes and extreme value theory. We also show that TV can perform exact segmentation in dimension one, under an alternating sign condition for some prescribed threshold. We apply our procedure to denoise a collection of 1D and 2D test signals verifying empirically the effectiveness of our approach. Codes are given to reproduce our results in a provided PURL.     

A Subpixel Image Restoration Algorithm

Abstract

In statistical image reconstruction, data are often recorded on a regular grid of squares, known as pixels, and the reconstructed image is defined on the same pixel grid. Thus, the reconstruction of a continuous planar image is piecewise constant on pixels, and boundaries in the image consist of horizontal and vertical edges lying between pixels. This approximation to the true boundary can result in a loss of information that may be quite noticeable for small objects, only a few pixels in size. Increasing the resolution of the sensor may not be a practical alternative. If some prior assumptions are made about the true image, however, reconstruction to a greater accuracy than that of the recording sensor's pixel grid is possible. We adopt a Bayesian approach, incorporating prior information about the true image in a stochastic model that attaches higher probability to images with shorter total edge length. In reconstructions, pixels may be of a single color or split between two colors. The model is illustrated using both real and simulated data.     

Continuous-time image reconstruction for binary tomography

Binary tomography is the process of reconstructing a binary image from a finite number of projections. We present a novel method for solving binary tomographic inverse problems using a continuous-time image reconstruction (CIR) system described by nonlinear differential equations based on the minimization of a double Kullback–Leibler divergence. We prove theoretically that the divergence measure monotonically decreases in time. Moreover, we demonstrate numerically that the quality of the reconstructed images of the nonlinear CIR system is better than those from an iterative reconstruction method.     

三维加权PFDK重建算法研究

基于滤波反投影的Feldkamp-Davis-Kress(FDK)算法,具有数学形式简单、容易实现和计算速度快等优点,在医疗和工业等领域得到了广泛的应用.平行重排(PF DK)算法是FDK算法的一种推广,针对PFDK算法重建出的图像受锥角的影响加大的问题,给出一种三维加权PFDK图像重建算法,并研究了重排过程中径向插值间隔对重建图像质量的影响,分别采用三种不同插值总数(插值间隔分别是1单位,0.5单位,0.25单位)重排数据.实验结果表明给出的三维加权PFDK算法可有效减少锥角的影响,且当采用2倍插值总数时重建结果较好.     

Algorithms and software for total variation image reconstruction via first-order methods

This paper describes new algorithms and related software for total variation (TV) image reconstruction, more specifically: denoising, inpainting, and deblurring. The algorithms are based on one of Nesterov’s first-order methods, tailored to the image processing applications in such a way that, except for the mandatory regularization parameter, the user needs not specify any parameters in the algorithms. The software is written in C with interface to Matlab (version 7.5 or later), and we demonstrate its performance and use with examples.