Next Article in Journal
Similarity Analysis of Methods for Objective Determination of Weights in Multi-Criteria Decision Support Systems
Next Article in Special Issue
On Nonlinear Forced Impulsive Differential Equations under Canonical and Non-Canonical Conditions
Previous Article in Journal
Disconnection between Rat’s Left and Right Hemisphere Impairs Short-Term Memory but Not Long-Term Memory
Previous Article in Special Issue
Group Analysis of the Boundary Layer Equations in the Models of Polymer Solutions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

CT Image Reconstruction via Nonlocal Low-Rank Regularization and Data-Driven Tight Frame

1
School of Mathematics and Big Data, Dezhou University, Dezhou 253023, China
2
Data Recovery Key Laboratory of Sichuan Province, College of Mathematics and Information Science, Neijiang Normal University, Neijiang 641100, China
3
Financial Department of Dezhou University, Dezhou 253023, China
*
Author to whom correspondence should be addressed.
Symmetry 2021, 13(10), 1873; https://doi.org/10.3390/sym13101873
Submission received: 29 August 2021 / Revised: 26 September 2021 / Accepted: 28 September 2021 / Published: 4 October 2021

Abstract

:
X-ray computed tomography (CT) is widely used in medical applications, where many efforts have been made for decades to eliminate artifacts caused by incomplete projection. In this paper, we propose a new CT image reconstruction model based on nonlocal low-rank regularity and data-driven tight frame (NLR-DDTF). Unlike the Spatial-Radon domain data-driven tight frame regularization, the proposed NLR-DDTF model uses an asymmetric treatment for image reconstruction and Radon domain inpainting, which combines the nonlocal low-rank approximation method for spatial domain CT image reconstruction and data-driven tight frame-based regularization for Radon domain image inpainting. An alternative direction minimization algorithm is designed to solve the proposed model. Several numerical experiments and comparisons are provided to illustrate the superior performance of the NLR-DDTF method.

1. Introduction

Medical imaging applications of X-ray computed tomography (CT) include cranial, chest, cardiac abdominal and pelvic imaging. The wide use of CT is due to the ability of viewing interior structures without destroying the surface of organs or subjects. However, the radiation caused by X-ray during imaging does harm to patients’ health. Thus, the approximately lossless image reconstruction based on low-dose X-ray, such as reducing the number of projections, is continuously considered by scientists and doctors. The basic model of the X-ray reconstruction problem can be represented as a linear inverse problem,
P u = f
where P R m × n is a projection operator representing the collection of discrete line integrations at different projection angles and along different beamlets. Numerically P is realized by Siddon’s algorithm [1,2].
When small number of projections are used, i.e., m < n , the matrix P has a rank deficiency. As a result, Equation (1) will have an infinite amount of solutions. There are many methods for deriving a useful result from them, such as the solutions by the filtered back-projection (FBP) algorithm [3] and algebraic reconstruction technique (ART) [4]; however, they usually contain artifacts and thus result in unreliable reconstructions.
In practice, the projection image f usually contains noise. Gaussian noise [1,5,6,7,8,9,10,11,12,13,14] and Poisson noise [15,16] both are considered to simulate a realistic environment. Thus, the basic model (1) is modified by
f = P u + ϵ .
In order to suppress noise and artifacts (local irregularities) during image reconstruction, it should usually contain tiny structural or high-frequency information to maintain fidelity and effectiveness. Various differential operator-based regularization methods, such as the total variation (TV)-based regularization approach (see [6,7,15,17] and references therein), low-rank models (see [8,18]) and wavelet frame-based methods (see [9,10,11,12]), have been proposed and received acceptable results for both, in vision and computational error. A standard form of the TV-based image reconstruction model can be written as
min u 1 2 P u f 2 2 + λ u 1 .
It is important to mention that the authors of [9] established a rigorous connection between a special model of wavelet frame-based approach, called the analysis-based approach, and variational models.
Two adjacent columns in the projected image f represent the information collected from adjacent projection angles. Therefore, f must be column-dependent locally. In order to exploit this kind of prior knowledge of f, the authors of [1] proposed an image inpainting method on the projection image f, following with another inverse problem for image restoration in the spatial domain. This wavelet frame-based regularity with Radon domain inpainting was shown to be useful for reconstructing high quality images from a very small number of projections,
min f , u 1 2 R Λ c ( P u f ) 2 2 + 1 2 R Λ ( P u ) f 0 2 2 + κ 2 R Λ f f 0 2 2 + λ 1 W 1 f 1 , 2 + λ 2 W 2 u 1 , 2 .
where R Λ C denotes the reconstruction on Ω \ Λ , and R Λ denotes the restriction on Λ . Here, f 0 is the projection image defined on the grid Λ of size N 0 × N p , where N 0 is the total number of detectors and N p is the number of angular projections. The unknown projection image f is defined on a grid Ω Λ of size N 0 × N ˜ p , with N ˜ p = 2 N p is selected in this paper. The first three terms of Equation (4) represent the data fidelity, making sure that f is consistent with f 0 on Λ and P f f . W 1 and W 2 are two different tight wavelet frame transforms. A fast algorithm was introduced to solve the model (4) based on the split Bregman algorithm [17,19,20] and the augmented Lagrangian method [21,22]. For more details about model and algorithm, we refer to [1].
Recently, it is known that in area of image restoration, data-driven tight frames or bi-frames generally outperform the regular pre-constructed wavelet frames [5,13,14]. The author of [5] proposed the SRD-DDTF model using data-driven tight frames as their sparsity priors for both u and f,
min f , u , v 1 , W 1 , v 2 , W 2 1 2 R Λ c ( P u f ) 2 2 + 1 2 R Λ ( P u ) f 0 2 2 + κ 2 R Λ f f 0 2 2 + λ 1 v 1 0 + μ 1 2 W 1 f v 1 2 2 + λ 2 v 2 0 + μ 2 2 W 2 u v 2 2 2 , s . t . W i T W i = I , i = 1 , 2 .
where tight frames W 1 and W 2 and corresponding coefficients v 1 and v 2 , are all treated as unknowns for enforcing the sparsity approximation to u and f. Here, 0 norm v i 0 stands for the number of non-zero elements in v i , i = 1 , 2 . For solving the additional W i and v i subproblems such as Equation (7), an alternative optimization algorithm was proposed in [13] based on Singular Value Decomposition (SVD) of a matrix associated with the image u (or f) after block rearrangement. As usual, the 0 regularization of v i can be solved by applying hard-thresholding to the singular values. A short review is stated in Section 2.1, for more details and convergence analysis, we refer to [5,13,14].
It is desired that there exist lots of similar structures both in natural and medical images. This property of nonlocal self-similarity was introduced by Dong et al. in [23,24] for exploiting structured sparsity in compressive sensing of both photographic and CT/MRI images. Random sampling and pseudo-radial sub-sampling are used for under sampling. However, they did not use any inpainting method in the transform domain for reconstruction. Xie et al. [16] used the nonlocal self-similar constraint in the position emission tomography (PET) image reconstruction with Total Variation(TV) method. The nonlocal low-rank approximation and TV model were also used for reducing noise in the denoising problem of low-dose CT image [25]. Xu et al. [26] proposed a model based on nonlocal low-rank and prior images under a given wavelet framework for reducing limited-angle artifacts in CT image restoration. There are also many other applications of nonlocal low-rank methods in the field of image denoising, one may see [27] for details.
In this paper, the nonlocal patch-based method will be used for restoration of CT image u. Roughly speaking, we shall use low-rank approximation for several image block patches after block matching. The rest of this paper is organized as follows. In Section 2, we will review image restoration methods based on data-driven tight frame and nonlocal low-rank regularization. In Section 3, we shall introduce our CT image restoration model combining data-driven tight frames for f and nonlocal low-rank regularization for u. In Section 4, several experiments will be provided for validating the merits of our proposed method, and the conclusion is discussed in Section 5.

2. Reviews and Preliminaries

2.1. Data-Driven Tight Frames

In [1,5], the sparse approximation methods were proposed for both the CT image u and projection image f, which are based on the pre-constructed wavelet frames or data-driven tight frames. When implemented, the data-driven tight frames are constructed based on the SVD of a matrix associated with the image u (or f) after block rearrangement.
It is well known that compared with a basis, using a frame results in a more robust representation of signals [28,29]. In discrete setting, the fast decomposition transform W and fast reconstruction transform W T , two finite matrix operators with entries from a filter { a j } j = 0 m associated to a translation-invariant tight wavelet transform, can be formed by convolution operators. Actually,
W T = [ S a 0 T ( · ) , S a 1 T ( · ) , , S a m T ( · ) ] T and W = [ S a 0 , S a 1 , , S a m ]
with convolution operator S a : 2 ( Z ) 2 ( Z ) ,
[ S a u ] ( n ) = [ a u ] ( n ) = k Z a ( n k ) u ( k ) for u 2 ( Z ) .
The tight frame property is guaranteed by W W T = I as usual, see [5]. The data-driven tight frame method takes the tight frame W as an unknown, which is also determined stepwisely by solving the following optimization problem:
min v , W λ v 0 + W f v 2 2 , W W T = I .
We also use a small-size W, thus, the phantom image f should be reshaped as the matrix F R ( N 1 N 2 ) × p . Here, each column is a vectorization of N 1 × N 2 patches extracted from f. The filters { a j } j = 0 m used for W are selected as the columns of the matrix D R ( N 1 N 2 ) × m . By denoting V R ( N 1 N 2 ) × p as the tight frame coefficients and taking m = N 1 N 2 , Equation (7) can be rewritten as
min V , D λ V 0 + D T F V 2 2 , D T D = I .
In the sequel, we follow the alternative optimization algorithm for solving Equation (8), see [5,13] for details,
D k + 1 = X Y T V k + 1 = T λ ( D k + 1 ) T F
with X and Y obtained by taking SVD of F ( V k ) T , and T λ being the hard-thresholding operator defined by
( T λ V ) [ i , j ] = 0 , if | x | < λ { 0 , V [ i , j ] } , if | x | = λ V [ i , j ] , otherwise .

2.2. Nonlocal Low-Rank Regularization

The nonlocal low-rank regularization is used for compressed sensing [24], PET image reconstruction [16] and CT image restoration from limited-angle projection data [26]. All these problems are handled by the alternative direction multiplier method technique, which is converted into several sub-problems and solved iteratively and alternately. When implemented, the rank-minimization sub-problem is also solved by the singular value thresholding (SVT) method, and the image reconstruction sub-problem is solved by different regularization methods.
In this article, the nonlocal low-rank regularization is used for CT image restoration, for the purpose of exploiting the nonlocal property in the image. After extracting n 1 × n 1 patches at position i in CT image u, denoted by x i R n 1 , m 1 -nearest patches based on Euclidean distance are used for obtaining the matrix X i = [ x i , 0 , x i , 1 , , x i , m 1 1 ] . Here, each x i , r is its reordered position index, satisfying S i = { r | x i x i , r T i } and T i is the pre-defined threshold.
Since in the natural image there exist lots of nonlocal self-similar structures, the formed data X i can be split into two parts, the low-rank matrix L i and the Gaussian noise matrix W i , i.e., X i = L i + W i . Then, L i can be recovered by the low-rank regularization model,
L i = argmin L i Rank ( L i ) , s . t . X i L i 2 2 σ W 2 ,
where σ W 2 denotes the variance of Gaussian noise. In order to solve Equation (11), as in [24], we consider an approximate version,
L i k + 1 = argmin L i X i L i 2 2 + λ 2 ϕ ( L i , σ ϵ k ) ,
where ϕ ( L i , σ ϵ k ) = r = 0 n 0 σ ϵ , r k σ r denotes the weighted nuclear norm with σ ϵ , r k = 1 σ r k + ϵ and n 0 = min { n 1 , m 1 } . The singular values σ r and σ r k corresponding to L i and L i k , respectively, are ordered in descending order. We adopt the weighted singular value thresholding method to solve Equation (12) efficiently and effectively,
L i k + 1 = S σ ϵ k ( X i ) = U Σ λ 2 diag ( σ ϵ k ) + V T ,
where U Σ V T is the SVD of X i , ( x ) + = max { x , 0 } and σ ϵ k = 1 σ 1 k + ϵ , 1 σ 2 k + ϵ , , 1 σ n 0 k + ϵ .

3. Models and Algorithm

3.1. CT Image Reconstruction Model

In this paper, we also focus on the column number of the unknown projection image satisfying N ˜ p = 2 N p , with N p being the column number of f 0 . f is considered to be restored from f 0 with doubled angular sampling resolution. Unlike the symmetric model (5), we use nonlocal low-rank regularization instead of the data-driven tight frame method when dealing with CT images. Our asymmetric CT image reconstruction model based on nonlocal low-rank regularity and data-driven tight frames reads as follows,
min u , f , W , v , L i 1 2 R Λ C ( P u f ) 2 2 + 1 2 R Λ P u f 0 2 2 + κ 2 R Λ f f 0 2 2 + λ 1 v 0 + μ 2 W f v 2 2 + η j 1 2 R ˜ i u L i 2 2 + λ 2 ϕ ( L i , σ ϵ k )
where R ˜ i u = R i , 0 u , R i , 1 u , , R i , m 1 1 u denotes the first m 1 closest patches under cosine similarity for every exemplar patch x i . In this section, we demonstrate that the proposed objective functional can be efficiently solved by the method of alternative minimization. After k-th iteration, the variables of step k + 1 are determined by the following sub-problems.
(1)
f sub-problem
If all the parameters but f are determined by k-th iteration, Equation (14) can be practically rewritten as
f k + 1 = argmin f 1 2 R Λ C ( P u f ) 2 2 + κ 2 R Λ f f 0 2 2 + μ 2 W k f V k 2 2 + a 2 f f k 2 2 .
Here, we add a special term f f k 2 2 in order to ensure that the new f k + 1 is not far from f k , which can theoretically justify the convergence of this algorithm. By differentiating the right side of Equation (15) with respect to f, we obtain the following closed-form solution
f k + 1 = R Λ C + κ R Λ + ( μ + a ) 1 ( R Λ C P u k + κ R Λ f 0 + μ W k T V k + a f k ) ,
where R Λ C + κ R Λ + ( μ + a ) is a diagonal matrix and, hence, Equation (16) can be computed easily.
(2)
W and v sub-problem
In this paper, we use the data-driven tight frame method explained in [5] or Section 2, for W and v sub-problems,
W k + 1 = argmin W T W = I μ 2 W f k + 1 v k 2 2 + c 2 W W k 2 2 v k + 1 = argmin v μ 2 W k + 1 f k + 1 v 2 2 + λ 1 v 0 + d 2 v v k 2 2
As discussed in the previous section, the variables f, v and W are reformulated in the practice calculation. To be more specific, we need to solve the following reformulated sub-problems:
D k + 1 = argmin D T D = I μ 2 D T F k + 1 V k 2 2 + c 2 D D k 2 2
V k + 1 = argmin V μ 2 ( D k + 1 ) T F k + 1 V 2 2 + λ 1 V 0 + d 2 V V k 2 2
Thus, to solve sub-problems (17) and (18), we can simply compute
D k + 1 = X Y T where X Σ Y T = F k + 1 ( V k + 1 ) T + c μ D k
and
V k + 1 = T λ 1 μ + d μ ( D k + 1 ) T F k + 1 + d V k μ + d ,
where T a ( · ) is the hard thresholding operator given in Equation (10).
(3)
L i and u sub-problems
Based on the augmented Lagrangian method, we formulate the L i sub-problem as,
L i k + 1 = argmin L i 1 2 R ˜ i u k L i 2 2 + λ 2 ϕ ( L i , σ ϵ k ) + e 2 L i L i k 2 2
with ϕ ( L i , σ ϵ k ) = r = 1 n 0 σ r σ r k + ϵ . Here, σ r and σ r k are the r-th ordered singular value with respect to L i and L i k . Applying the weighted singular value thresholding operator to Equation (21) yields the algorithm
L i k + 1 = S σ ϵ k R ˜ i u k + e L i k = U Σ ˜ λ 2 1 + e diag ( σ ϵ k ) + V T
where U Σ ˜ V T is the SVD of R ˜ i u k + e L i k .
After L i being updated, the current estimation of CT image u can be solved by our proposed model with other variables fixed,
u k + 1 = argmin u 1 2 R Λ C P u f k + 1 2 2 + 1 2 R Λ P u f 0 2 2 + η 2 k R ˜ i u L i k + 1 2 2 + b 2 u u k 2 2 ,
where the term u u k 2 2 stands for the similarity in two adjacent steps. It is a quadratic optimization problem admitting a closed-form solution,
u k + 1 = P T P + η i R ˜ i T R ˜ i + b I 1 P T R Λ C f k + 1 + P T R Λ C f 0 + η i R ˜ i T L i k + 1 + b u k ,
where R ˜ i T L i k + 1 = r = 0 m 1 R i , r T L i , r k + 1 and R ˜ i T R ˜ i = r = 0 m 1 R i , r T R i , r is a diagonal matrix. In order to derive a much faster algorithm for solving Equation (23), we also use the alternative direction multiplier method (ADMM) [30]. Thus, we obtain
u k + 1 = argmin u 1 2 R Λ C P u f k + 1 2 2 + 1 2 R Λ P u f 0 2 2 + b 2 u u k 2 2 + β 2 u z + μ β 2 2 + η 2 i R ˜ i z L i k + 1 2 2 ,
where z R n × n is an auxiliary variable, u R n × n is the Lagrangian multiplier, and β is a positive scalar parameter. The split iteration version of the optimization of Equation (24) can be stated as follows:
z k + 1 = argmin z β k 2 u k z + μ k β k 2 2 + η 2 i R ˜ i z L i k + 1 2 2 u k + 1 = argmin u 1 2 R Λ C P u f k 2 2 + 1 2 R Λ P u f 0 2 2 + b 2 u u k 2 2 + β k 2 u z k + 1 + μ k β k 2 2 μ k + 1 = μ k + β k ( μ k z k + 1 ) β k + 1 = 1.2 β k .
Here, z k + 1 admits a closed-form solution,
z k + 1 = η i R ˜ i T R ˜ i + β k I 1 β k u k + μ k + η i R ˜ i T L i k + 1 ,
i R ˜ i T R ˜ i is a diagonal matrix with each diagonal element being the number of overlapping patches that cover the current pixel location. i R ˜ i T L i k + 1 denotes the weighted low-rank approximation, i.e., gathering useful information from the collected similar patches.
For u k + 1 , we also use the conjugate gradient method to solve a linear equation,
P T P + ( b + β k ) I u = P T R Λ C f k + 1 + μ k + η i R ˜ i T L i k + 1
Collecting all together, we derive the following Algorithm 1.
Algorithm 1: CT image reconstruction via nonlocal low-rank regularity and data-driven tight frame (NLR-DDTF).
1Input: compute u 0 from (5) and set f 0 = P u 0 . Compute v 0 , w 0 from (4).
2Repeat:
3  (1): update f k + 1 by optimizing f using (16)
4  (2): update D k + 1 and V k + 1 by (19) and (20)
5  (3): update L i k + 1 ,for all i using (22)
6  (4): update u k + 1 by solving (24)
7Until: Relative error u k + 1 u k / u k + 1 < ϵ
8Output:  u k

4. Experiments

In this section, we report the experimental results of the proposed asymmetric low-rank approximation and data-driven tight frame-based asymmetric CT image recovery method. The proposed joint regulation method is actually based on the SRD-DDTF model [5]. The major difference is the optimization method for updating u k . Thus, we select the same warm-started initial value u 0 as in [5], given by analysis of the model in Section 3.1. In the following, we consider three phantoms, head image with smaller details from [5], head image and brain image with much more details from [24]. All experiments are conducted on a 3.20 GHZ Intel(R) Core(TM) i7-8700cpu with 16 GB memory. The PSNR, relative error and correlation results of the reconstructed images are included in Table 1, additionally, the visual effects are reported in Figure 1, Figure 2 and Figure 3. Both of them reveal that the nonlocal low-rank method does preserve more information from the underground truth other than noise, which contains less structural information.
The degrade projection image f 0 is constructed by the every-other-angular extraction method from synthesized data, based on the Radon transform of real data u and the Monte Carlo simulation, f 0 = R Λ ( P u + ϵ ) , where ϵ is some Gaussian white noise and the projection number N p is considered to be 15 , 30 , 45 , 60 in all examples. We choose the standard deviation of noise as max ( | P u | ) / 300 in our experiments. The projection image f is calculated based on the photon numbers the detectors received. Gaussian white noise is also generated when the corresponding optical equipment counts photons, which is related to the uncertainty of photons. For an initial value of f, we directly use P u 0 as in [5].
The patch size for u is 6 × 6 and the patch size for f is 8 × 2 to adapt to the different shape of u and f. For optimization of the W and v subproblems, we use μ = 8400 λ 1 and λ 1 is selected the same as in the SRD-DDTF model (5). When splitting the image u into overlapped small patches, we extract a paragon patch in every N s pixels along both horizontal and vertical directions. Here, we use N s = 5 and N s = 3 for realizing our proposed NLR-DDTF model. For each patch of u, m 1 = 41 similar patches are selected when the low-rank approximation is executed to derive L i . In our application, we choose η = 0.8 for smaller details in the image “head1” and η = 1.6 for the many more details in images “head2” and “brain”. The regularization parameter λ 2 is selected to be 1.
For quantitative accuracy evaluation, the relative error and correlation for the reconstructed u corresponding to the ground truth image u t are defined as:
e r r ( u t , u ) = u u t 2 u t 2 and c o r r ( u t , u ) = ( u u ¯ ) ( u t u ¯ t ) u u ¯ 2 u t u ¯ t 2
with u ¯ and u ¯ t being the mean of u and u t . The maximum iteration number of all our experiments is set to 1000 with the following stopping criterion: whenever the relative change between two consecutive iterations is below 10 3 (or 10 4 for “head1” with higher N p than 15).

5. Conclusions

In this work, we introduced a novel spatial-Radon domain CT image reconstruction model, which combines the spatial domain inverse problem based on the nonlocal low-rank approximation method and the Radon domain inpainting model using a data-driven tight frame-based regularization. The proposed model reconstructs high quality and high resolution CT images even when the images contain more details than the observed low-resolution projection images. We used the split Bregman algorithm for numerical simulations of the proposed model.
The success of the regularization method in CT image restoration, no matter TV-based or wavelet frame-based, lies in the property of local smoothness or structure sparsity of the ground truth image. It is shown in [5] that using data-driven tight frame as sparsity priors for Radon domain projection image preformed better than using pre-determined wavelet frame systems. The nonlocal low-rank prior can reflect the group sparsity characteristic of similar patches in CT images well. Our proposed CT image restoration method combines the advantages of the above two methods. The PSNR and subjective quality comparison shown in Table 1 clearly reveal that the proposed NLR-DDTF method performs better than just the tight frame in several cases. Compared with the SRD-DDTF method, our proposed NLR-DDFT model’s shortcoming is the slightly higher time consumption. When N s = 5 , the average time required for all examples is 42 % longer, and the average time consumption is 120 % times higher if N s = 3 .

Author Contributions

Models, methodology and formal analysis, Y.S. and S.S.; conceptualization and data curation, F.X.; convergence and supervision, Y.L.; programming and experiments, Y.S. and F.X.; writing—original draft preparation, Y.S. and S.S.; writing—review and editing, Y.L., X.Y. and X.Z.; project administration, F.X., Y.L. and X.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Science and Technology Research Program for Universities in Shandong Province, China (Grant No. J16LI58 and J18KB099), Natural Science Foundation of Shandong Province (Grant No. ZR2018LA003), Open Research Fund Program of Data Recovery Key Laboratory of Sichuan Province (Grant No. DRN19020), National Natural Science Foundation of China (No. 11801060, 11501082) and Research and Development Project of Dezhou City in China.

Acknowledgments

The authors would like to thank B. Dong and W. Dong for providing the test images and source code used in the simulation experiments.

Conflicts of Interest

The authors declare no conflict of interest. The funders S.Y., L.Y. and Y.X. had a role in the design of the study, in the analyses and interpretation of data, in the writing of the manuscript, and in the decision to publish the results.

References

  1. Dong, B.; Li, J.; Shen, Z. X-ray ct image reconstruction via wavelet frame based regularization and radon domain inpainting. J. Sci. Comput. 2013, 54, 333–349. [Google Scholar] [CrossRef]
  2. Siddon, R. Fast calculation of the exact radiological path for a three-dimensional CT array. Med. Phys. 1985, 12, 252. [Google Scholar] [CrossRef] [PubMed]
  3. Feldkamp, L.; Davis, L.; Kress, J. Practical cone-beam algorithm. J. Opt. Soc. Am. A 1984, 1, 612–619. [Google Scholar] [CrossRef] [Green Version]
  4. Gordon, R.; Bender, R.; Herman, G. Algebraic reconstruction techniques (ART) for three-dimensional electron microscopy and X-ray photography. J. Theor. Biol. 1970, 29, 471–481. [Google Scholar] [CrossRef]
  5. Zhan, R.; Dong, B. CT image reconstruction by spatial-Radon domain data-driven tight frame regulatization. SIAM J. Imaging Sci. 2016, 9, 1063–1083. [Google Scholar] [CrossRef] [Green Version]
  6. Rudin, L.; Osher, S.; Fatemi, E. Nonlinear total variation based noise removal algorithms. Phys. D Nonlinear Phenom. 1992, 60, 259–268. [Google Scholar] [CrossRef]
  7. Sidky, E.; Pan, X. Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization. Phys. Med. Biol. 2008, 53, 4777. [Google Scholar] [CrossRef] [Green Version]
  8. Cai, J.; Candes, E.; Shen, Z. A singular value thresholding algorithm for matrix completion. SIAM J. Optim. 2010, 20, 1956–1982. [Google Scholar] [CrossRef]
  9. Cai, J.; Dong, B.; Osher, S.; Shen, Z. Image restorations: Total variation, wavelet frames and beyond. J. Am. Math. Soc. 2012, 25, 1033–1089. [Google Scholar] [CrossRef] [Green Version]
  10. Cai, J.; Dong, B.; Shen, Z. Image restorations: A wavelet frame based model for piecewise smooth functions and beyond. Appl. Comput. Harmon. Anal. 2016, 41, 94–138. [Google Scholar] [CrossRef] [Green Version]
  11. Dong, B.; Jiang, Q.; Shen, Z. Image Restoration: Wavelet Frame Shrinkage, Nonlinear Evolution PDEs, and Beyond. Multiscale Model. Simul. 2017, 15, 606–660. [Google Scholar] [CrossRef] [Green Version]
  12. Dong, B.; Shen, Z.; Xie, P. Image restoration: A general wavelet frame based model and its asymptotic analysis. SIAM J. Math. Anal. 2017, 49, 421–445. [Google Scholar] [CrossRef] [Green Version]
  13. Cai, J.; Ji, H.; Shen, Z.; Ye, G. Data-driven tight frame construction and image denoising. Appl. Comput. Harmon. Anal. 2014, 37, 89–105. [Google Scholar] [CrossRef]
  14. Bao, C.; Ji, H.; Shen, Z. Convergence analysis for iterative data-driven tight frame construction scheme. Appl. Comput. Harmon. Anal. 2015, 38, 510–523. [Google Scholar] [CrossRef]
  15. Chen, Z.; Jin, X.; Li, L.; Wang, G. A limited-angle ct reconstruction method based on anisotropic tv minimization. Phys. Med. Biol. 2013, 58, 2019. [Google Scholar] [CrossRef]
  16. Xie, N.; Chen, Y.; Liu, H. Nonlocal low-rank and total variation constrained PET image reconstruction. In Proceedings of the 2018 24th International Conference on Pattern Recognition (ICPR), Beijing, China, 20–24 August 2018; pp. 3874–3879. [Google Scholar] [CrossRef]
  17. Zhang, H.; Wang, L.; Yan, B.; Li, L.; Xi, X.; Lu, L. Image reconstruction based on total-variation minimization and alternating direction method in linear scan computed tomography. Chin. Phys. B 2013, 22, 078701. [Google Scholar] [CrossRef]
  18. Cai, J.; Jia, X.; Jiang, H.G.S.; Shen, Z.; Zhao, H. Cine cone beam ct reconstruction using low-rank matrix factorization: Algorithm and a proof-of-principle study. IEEE Trans. Med. Imaging 2014, 33, 1581–1591. [Google Scholar] [CrossRef] [PubMed]
  19. Goldstein, T.; Osher, S. The split Bregman algorithm for L1 regularized problems. SIAM J. Imaging Sci. 2009, 2, 323–343. [Google Scholar] [CrossRef]
  20. Cai, J.; Osher, S.; Shen, Z. Split Bregman methods and frame based image restoration. Multiscale Model. Simulation A SIAM Interdiscip. J. 2009, 8, 337–369. [Google Scholar] [CrossRef]
  21. Esser, E. Applications of Lagrangian-based alternating direction methods and connections to split Bregman. CAM Rep. 2009, 9, 31. [Google Scholar]
  22. Tai, X.; Wu, C. Augmented Lagrangian method, dual methods and split Bregman iteration for ROF model. Scale Space Var. Methods Comput. Vis. 2009, 5567, 502–513. [Google Scholar]
  23. Dong, W.; Li, X.; Zhang, L.; Shi, G. Sparsity-based image denoising via dictionary learning and structural clustering. In Proceedings of the 2011 IEEE Conference on Computer Vision and Pattern Recognition(CVPR), Colorado Springs, CO, USA, 20–25 June 2011; pp. 457–464. [Google Scholar] [CrossRef]
  24. Dong, W.; Shi, G.; Li, X.; Ma, Y.; Huang, F. Compressive sensing bia nonlocal low-rank regularizaiton. IEEE Trans. Image Process. 2014, 23, 3618–3632. [Google Scholar] [CrossRef] [PubMed]
  25. Sagheer, S.; George, S. Denoising of low-dose CT images via low-rank tensor modelling and total variation regularization. Artif. Intell. Med. 2019, 94, 1–17. [Google Scholar] [CrossRef] [PubMed]
  26. Xu, N.; Hou, H.; Cheng, Z.; Wang, M.; Wang, Y.; Wang, G. Nonlocal low-rank and prior image-based reconstruction in a wavelet tight frame using limited-angle projection data. IEEE Access 2021, 9, 24616–24628. [Google Scholar] [CrossRef]
  27. Hu, Z.; Nie, F.; Wang, R.; Li, X. Low rank regularization: A review. Neural Netw. 2021, 136, 218–232. [Google Scholar] [CrossRef] [PubMed]
  28. Li, Y.; Han, D. Framelet perturbation and application to nouniform sampling approximation for Sobolev space. arXiv 2017, arXiv:1707.01325. [Google Scholar]
  29. Shen, Y.; Yang, S.; Yuan, D. d-Refinable (dual) pseudo-splines and their regularities. Int. J. Wavelets Multiresol. Inf. Process. 2017, 15, 1750002. [Google Scholar] [CrossRef]
  30. Lin, Z.; Chen, M.; Wu, L.; Ma, Y. The Augmented Lagrange Multiplier Method for Exact Recovery of Corrupted Low-Rank Matrices; Tech. Rep. UILU-ENG-09-2215; Department of Electrical and Computer Engineering, University of Illinois Urbana-Champaign: Urbana, IL, USA, 2009. [Google Scholar]
Figure 1. The tomographic results for the image “head1”. (top left) True data and local detail area for comparison; (top right) the restored image by SRD-DDTF method [5]; (bottom left) the restored image by NLR-DDTF method with N s = 5 ; (bottom right) the restored image by NLR-DDTF method with N s = 3 .
Figure 1. The tomographic results for the image “head1”. (top left) True data and local detail area for comparison; (top right) the restored image by SRD-DDTF method [5]; (bottom left) the restored image by NLR-DDTF method with N s = 5 ; (bottom right) the restored image by NLR-DDTF method with N s = 3 .
Symmetry 13 01873 g001
Figure 2. The tomographic results for the image “head2”. (top left) True data and local detail area for comparison; (top right) the restored image by SRD-DDTF method [5]; (bottom left) the restored image by NLR-DDTF method with N s = 5 ; (bottom right) the restored image by NLR-DDTF method with N s = 3 .
Figure 2. The tomographic results for the image “head2”. (top left) True data and local detail area for comparison; (top right) the restored image by SRD-DDTF method [5]; (bottom left) the restored image by NLR-DDTF method with N s = 5 ; (bottom right) the restored image by NLR-DDTF method with N s = 3 .
Symmetry 13 01873 g002
Figure 3. The tomographic results for the image “brain”. (top left) True data and local detail area for comparison; (top right) the restored image by SRD-DDTF method [5]; (top right) the restored image by NLR-DDTF method with N s = 5 ; (bottom right) the restored image by NLR-DDTF method with N s = 3 .
Figure 3. The tomographic results for the image “brain”. (top left) True data and local detail area for comparison; (top right) the restored image by SRD-DDTF method [5]; (top right) the restored image by NLR-DDTF method with N s = 5 ; (bottom right) the restored image by NLR-DDTF method with N s = 3 .
Symmetry 13 01873 g003
Table 1. Comparison of relative errors (in percentage), correlations (in percentage), PSNR (in db) and running time (in seconds). For image “head1”, NLR-DDTF methods can improve PSNR 3.8–27.2% with consuming 15.7–240.2% more time; For image “head2”, NLR-DDTF methods can improve PSNR 0.2–24.6% with consuming 29.1–233.7% more time; for image “brain”, NLR-DDTF methods can improve PSNR 1.3–11.0% with consuming 18.9–225.7% more time.
Table 1. Comparison of relative errors (in percentage), correlations (in percentage), PSNR (in db) and running time (in seconds). For image “head1”, NLR-DDTF methods can improve PSNR 3.8–27.2% with consuming 15.7–240.2% more time; For image “head2”, NLR-DDTF methods can improve PSNR 0.2–24.6% with consuming 29.1–233.7% more time; for image “brain”, NLR-DDTF methods can improve PSNR 1.3–11.0% with consuming 18.9–225.7% more time.
“head1”
N p SRD-DDTFNLR-DDTF with N s = 5 NLR-DDTF with N s = 3
errcorrpsnrtimeerrcorrpsnrtimeerrcorrpsnrtime
1510.5099.0635.81470.608.9799.1337.16860.718.1099.4438.031600.87
305.4099.7541.57923.972.0399.9650.051251.051.9399.9750.501937.15
454.2899.8443.581367.351.2899.9654.051672.661.2299.98754.482368.40
603.8099.8844.611809.700.9899.9956.382095.440.94199.99256.742782.54
“head2”
N p SRD-DDTFNLR-DDTF with N s = 5 NLR-DDTF with N s = 3
errcorrpsnrtimeerrcorrpsnrtimeerrcorrpsnrtime
1529.3892.3422.47495.4129.2192.3922.52888.5129.2292.3922.521653.06
3020.7196.2825.51914.7615.6697.8727.931284.9515.4497.9328.051953.55
4516.9497.5427.251378.7010.0499.1331.801718.879.9099.1531.912443.62
6015.3797.9828.091773.867.0199.5834.912290.076.93999.5835.002861.41
“brain”
N p SRD-DDTFNLR-DDTF with N s = 5 NLR-DDTF with N s = 3
errcorrpsnrtimeerrcorrpsnrtimeerrcorrpsnrtime
1523.8795.5128.73476.9522.8795.8829.10909.4422.5096.0229.241553.29
3015.5798.1232.44944.7912.3498.8234.461284.8112.0998.8734.632059.53
4512.9798.7034.031390.109.3199.3336.911728.389.2299.3436.992453.39
6011.6098.9635.001777.777.5199.5638.772112.937.4599.5738.842838.36
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shen, Y.; Sun, S.; Xu, F.; Liu, Y.; Yin, X.; Zhou, X. CT Image Reconstruction via Nonlocal Low-Rank Regularization and Data-Driven Tight Frame. Symmetry 2021, 13, 1873. https://doi.org/10.3390/sym13101873

AMA Style

Shen Y, Sun S, Xu F, Liu Y, Yin X, Zhou X. CT Image Reconstruction via Nonlocal Low-Rank Regularization and Data-Driven Tight Frame. Symmetry. 2021; 13(10):1873. https://doi.org/10.3390/sym13101873

Chicago/Turabian Style

Shen, Yanfeng, Shuli Sun, Fengsheng Xu, Yanqin Liu, Xiuling Yin, and Xiaoshuang Zhou. 2021. "CT Image Reconstruction via Nonlocal Low-Rank Regularization and Data-Driven Tight Frame" Symmetry 13, no. 10: 1873. https://doi.org/10.3390/sym13101873

APA Style

Shen, Y., Sun, S., Xu, F., Liu, Y., Yin, X., & Zhou, X. (2021). CT Image Reconstruction via Nonlocal Low-Rank Regularization and Data-Driven Tight Frame. Symmetry, 13(10), 1873. https://doi.org/10.3390/sym13101873

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop