Total Fractional-Order Variation-Based Constraint Image Deblurring Problem

.


Introduction
In image processing, image deblurring is an attractive topic due to its practical applications in robot vision [1], remote sensing [2], medical image processing [3], virtual reality [4], astronomical imaging [5], and many other fields.The mathematical relationship between the original u and blurry z images is follows: where ε denotes a noise function and K denotes the blurring operator: where k(x, y) = φ(x − y) is referred to as a translation-invariant kernel or a point spread function (PSF).Therefore, the task of recovering the u and K from z is called the deconvolution problem.If the blurring operator K is given, then the corresponding approach is referred to as non-blind deconvolution [6][7][8][9][10][11][12][13].However, when the blurring operator is unknown, the corresponding approach is referred to as blind deconvolution [14][15][16][17][18][19][20].In this paper, our primary focus is on non-blind deconvolution.K is the compact operator.However, the recovering of u from z poses challenges as it transforms the problem (1) into an unstable inverse problem [21][22][23].To address this issue, researchers have extensively explored the potential of energy minimization models to solve image deblurring problems, which have attracted significant attention over the past few decades.
where C represents a constrained set, R(u) is a regularization functional, and α > 0 is a smoothing parameter that determines the balance between the data fitting and smoothing terms, and * is a 2-D convolution operator.When applying these techniques to noisy and blurry photos, researchers must overcome two major challenges.The first challenge is dealing with non-linearity, while the second challenge involves resolving the massive matrix system involved.

Related Works
Non-blind deconvolution poses significant challenges as an ill-posed inverse problem.Numerous techniques for deconvolution have been developed to address these challenges by incorporating different image priors to regularize the solution.One such example is the utilization of the Tikhonov regularization model [21,22], where u Tik = Ω |u|dΩ.The Tikhonov model involves least-squares estimation, which often leads to excessively smoothed image reconstructions.As a result of its edge-preserving property, the total variation (TV) model [23][24][25] has become the most widely recognized non-linear energy minimization image deblurring model used for image deblurring.
where u TV,β = Ω | u| β dΩ and | u| β = u 2 x + u 2 y + β.Here, β > 0 is used to make functional u TV,β differentiable at zero.The TV model possesses numerous advantageous features; however, it does have one significant flaw.One notable drawback of the TV model is its tendency to transform smooth functions into piecewise constant functions, resulting in staircase effects in the restored images.The repaired photos appear blocky as a result.To mitigate the staircase effects in restored images, one solution is to employ total fractional-order variation (TFOV)-based models [26][27][28][29][30][31][32].These models have been proposed as an alternative approach to address the limitations of the TV model and reduce the undesirable staircase artifacts.
where α represents order of fractional derivative.u TFOV,β = Ω | ∇ α u | 2 +β 2 dΩ, and , where D α x and D α y are the fractional derivative operators along the x and y directions, respectively.The x-direction derivative is also denoted by D α a,x .Here, a and x are the lower and upper bounds of the integrals, and α represents the order of fractional derivative.In this notation, 0 < n − 1 < α < n, where n = [α] + 1 and [•] denotes the greatest integer function.Various definitions have been proposed to define a fractional-order derivative [33][34][35][36].The regularization models based on TFOV are known for their exceptional efficiency.These models preserve edges in the recovered images while simultaneously eliminating the undesirable staircase effect.
In recent years, the utilization of TFOV-based image processing methods has attracted increasing attention.Some basic conclusions have been drawn in the areas of image denoising, edge detection, and reconstruction [31,37,38].Mathieu et al. [37] proposed an edge detection method based on fractional differential, which effectively enhances image details such as edges and textures.Tian et al. [38] proposed a fractional-order adaptive regularization primal-dual algorithm for image denoising.Furthermore, Zhang et al. [31] proposed a TFOV model for image restoration, demonstrating its efficacy in suppressing the staircase effect.These studies have shown that, compared to first-order and second-order total variation methods, TFOV can more accurately and delicately represent image textures.More recently, Fairag et al. [39] and Guo et al. [30] have incorporated the TFOV model within the framework of image deblurring problems, further highlighting the applicability and potential of the TFOV-based approaches.

Scope of the Paper
Particularly in astronomical images, image deblurring frequently requires that the restored image has precisely non-negative intensities [40][41][42][43][44].However, it has been noted that solutions using current techniques may not always produce favorable outcomes.Images with many pixels having intensity values equal to or close to zero are known as images with negative intensities or black space.In this research, we provide a model for TFOV-based image deblurring that guarantees strictly positive outcomes for image intensities.The suggested model additionally restricts the image intensity values, maintaining them within a specified range.The removal of negative intensities or their confinement within the prescribed range also contributes to improving the quality of deblurred images.The main idea behind this paper is to convert the TFOV-based constrained image deblurring problem into an unconstrained one and then introduce the Lagrange multiplier.The optimization issues in computer vision and image processing have been successfully addressed by augmented Lagrangian methods [45][46][47].Augmented Lagrangian methods have demonstrated superior speed compared to other numerical techniques.It has been demonstrated that the original nontrivial minimization problem can be broken down into a number of straightforward and quick-to-solve subproblems using augmented Lagrangian methods.Some of them have closed forms of solutions, while others can be quickly solved using tools such as the fast Fourier transform (FFT).In our augmented Lagrangian method, the solution of one of the subproblems requires the conjugate gradient (CG) method.However, the CG method exhibits slow convergence due to an ill-conditioned matrix system.To overcome the slow convergence problem of CG, we introduce a new preconditioned matrix in this paper.
The main contributions of this paper are as follows: (i) it presents the one-sided and two-sided constraint methods for the TFOV-based constrained image deblurring problem; (ii) the proposed methods limit the upper boundary of the image intensity values, maintaining them within a specified range, while also guaranteeing strictly positive results; (iii) it presents a new circulant preconditioned matrix to improve the convergence of the CG method within the augmented Lagrangian method; and (iv) the proposed methods generate high-quality restored images compared to the most recent existing TFOV-based image deblurring methods.
This paper is divided into different sections.We discuss one-sided and two-sided constraint problems in Section 2. Section 3 presents Euler-Lagrange equations.Section 4 presents the cell discretization and the matrix-system of the model.The proposed preconditioned matrix is also presented in Section 4. The numerical application of our approaches is presented in Section 5. Section 6 contains the conclusions regarding the suggested methods and the Appendix A.

Constraint Image Deblurring Problem
In this section, we present a model for TFOV-based image deblurring that guarantees strictly positive image intensities as an outcome.The proposed model also imposes restrictions on the upper boundary of image intensity values, constraining them within a specified range.The main idea of this model is to convert the TFOV-based constrained image deblurring problem into an unconstrained problem and subsequently introduce a Lagrange multiplier.

One-Sided Constraint Problem
Consider the constrained image deblurring problem subject to : u ≥ 0 We convert the inequality constrained (8) into an equality constrained by introducing a function γ −u Now, Equations ( 7) and ( 8) become subject to : We have that u * is a local (global) minimum of Equations ( 7) and ( 8) if and only if (u * , γ * ), where γ * = √ u * is local (global) minimum of Equations ( 10) and (11).Now, consider the augmented Lagrangian functional for Equations (10) and (11) defined for positive penalty parameter c > 0 and multiplier λ by We are interested in minimizing augmented Lagrangian Equation (12) with respect to (u, γ) for different λ and c.Observe that the minimization of f c (u, γ, λ) with respect to γ can be found explicitly for each fixed u [48].The minimization of the problem above with respect to γ is equivalent to The above integrand is quadratic in w.The unconstrained (global) minimum at which the derivative is zero is w * .We have Therefore, the solution to Equation ( 13) is Substituting w * into the functional (12) gives Now, the method of multiplier [49,50] can be described as follows: For a given multiplier λ (k) and the penalty parameter c (k) , we minimize f c (k) (u, γ, λ (k) ) by obtaining u (k)  and γ (k) ; subsequently, we set In order to minimize the functional (12), we first pick a value for c and a function λ, then we compute w using Equation (15).Next, we compute u by minimizing the functional (12).This suggests following the one-sided constraint method (OSCM) Algorithm 1. Set: Set: Set: For m = 1, 2, ...

Two-Sided Constraint Problem
Next, we take the case where pixel values of digital images must lie in a specific interval [a 1 , a 2 ].For instance, for 8-bit images, the interval is [a 1 , a 2 ] = [0, 255].We consider solving the constrained model: subject to : First, we convert the inequality (18) into two equalities Then, Equations ( 17) and ( 18) become subject to : Let us consider the augmented Lagrangian functional for Equations ( 20)-( 22) defined for a positive penalty parameter c > 0 and multipliers λ 1 , λ 2 by Now, we want to minimize augmented Lagrangian (23) with respect to (u, and c.
Similar to the one-side constraint case, minimization of g c with respect to γ 1 and γ 2 can be explicitly found for each fixed u.Minimization with respect to γ 1 and γ 2 are equivalent to min min Hence, the solutions of Equations ( 24) and ( 25) are The method of the multiplier can be described as follows: Given multipliers λ 2 , and a penalty parameter c (k) , we minimize g c (k) by obtaining u (k) , γ (k) In order to minimize the functional (23), we first select a value for c and a function λ 1 .Then, we compute w * 1 using Equation (26).After that we choose a function λ 2 and compute w * 2 using Equation (27).Next, we compute u by minimizing the functional ( 23).This two-sided constraint method (TSCM) is explained in Algorithm 2.

Algorithm 2 TWO-SIDED CONSTRAINT METHOD function:
Set: Set: Set: Test: Stopping criteria 10. Set: 11. Set: In both algorithms (OSCM and TSCM), we have to compute u by minimizing the functionals ( 12) and (23), respectively.In both cases, we require Euler-Lagrange equations to compute u.The Euler-Lagrange equations for u are the same for both cases.

Euler-Lagrange Equations
This section includes the presentation of the Euler-Lagrange equations connected with TFOV for the image deblurring problem.Theorem 1.For the functional given in Equation ( 12) and α ∈ (1, 2), the Euler-Lagrange equations are where η is the unit outward normal vector.The K * is an adjoint operator of the integral operator K and non-linear differential operator L α (u) is defined as follows: Proof.The proof is given in Appendix A.
Note that Equation ( 30) can be written as follows: with the dual, or flux, variable We apply the Galerkin method to Equations ( 32) and ( 33) together with midpoint quadrature for the integral term and the cell-centered finite difference method for the derivative part.

Numerical Implementation
First, we will present the discretization of our proposed model.The computational domain Ω = (0, 1) × (0, 1) is divided into N 2 equal squares (cells), where N represents the number of equispaced partitions in the x or y direction.We proceed the same discretization approach in [31,51].Next, let (x k , y l ), k, l = 0, 1, ..., N + 1 be discrete points for the image domain Ω.We assume that u satisfies homogenous Dirichlet boundary condition.To discretize the fractional derivative of order α at the inner point (x k , y l ) (for k, l = 0, 1, ..., N) in the x-direction, we employ the shifted Gr ünwald approximation approach [52].
which is applicable to both Riemann-Liouville and Caputo derivatives [53,54].Here, f l s = f s,l and ω α j = (−1) j α j j = 0, 1, ..., N and ω α 0 = 1, ω α j = (1 − 1+α j )ω α j−1 for j > 0. From Equation (35), one can observe that the first-order approximation of with fixed y l .After using the homogenous boundary condition in the matrix estimation of the fractional-order derivative, all N equations of the fractional derivatives in the x-direction in Equation ( 35) can be expressed as follows: .
By the definition of fractional derivative (35), for any 1 < α < 2, the coefficients ω k α possess the properties given below: (1) By applying the Gershgorin circle theorem, it can be concluded that the matrix B α N is a symmetric and negative-definite Toeplitz matrix (i.e., −B α N is the positive definite Toeplitz matrix).Let U ∈ R N×N represent the solution matrix at all nodes (kh x ; lh y ), k, l = 1, ..., N corresponding to the spatial discretization nodes in the x and y directions.The ordered solution vector of U is denoted by u ∈ R N 2 ×1 .The discrete and direct analogue to differentiation for an arbitrary order α derivative is In the same way, all values in the y-direction having order α derivative of u(x; y) for these nodes are estimated using where u = u 11 , u 12 , ..., u NN and ⊗ represents the Kronecker product.The αth-order derivative of u x α of u(x; y) along all x-direction nodes in Ω can be represented by the matrix B α N U.For further details on the discretization, we recommend reading References [35,54].The fractional discretization mentioned above utilizes cell-centered finite difference method (CCFDM) and takes advantage of the fact that [(−1) n ∇ α •] is the adjoint operator of the operator ∇ α .Consequently, Equations ( 32) and (33) yield the following system where N F is the number of the Fixed Point Iterations (FPI) used to linearize the non-linear term in the square root in (34).The matrix K h is obtained by using the midpoint quadrature for the integral operator as follows: with entries [K h U] ij,lm = h 2 k(x i − x j , y l − y m ), where the lexico-graphical order is used, K h is a block Toeplitz with Toeplitz block (BTTB) matrix.The discrete scheme of the matrix L α h U is given by (41)   where • is the point-wise multiplication, m is the m-th fixed-point iteration, and U is an N × N-sized reshaped matrix of the vector u.D 1 (U m ) and D 2 (U m ) are the diagonals of the Hadamard inverses of B x α (U m ) and B y α (U m ), respectively.I N is the identity matrix.Now, if we eliminate V from the system (39), then we have the following primal system of the TFOV-based image deblurring model: If we use a simple total variation (TV) regularization functional, then we have the following similar primal form: where The L TV h (U m ) is derived from the discretization of total variational functional.The details can be seen in Reference [22].The matrix B h has the following structure, where both G 1 and G 2 are of size n x (n x − 1) × n 2 x , and H h is a diagonal matrix obtained by the discretization of the expression | u| 2 + β 2 , which has the following structure: where H x is a size of (n x − 1) × n x , and H y is a size of n x × (n x − 1).
To obtain the value of our primal variable u, one needs to solve the matrix system (42).As for the remaining variables and Lagrange multipliers, we solve them directly after discretizing them at the grid points (x i , y j ).
As mentioned earlier, to compute the value of u, we need to solve the matrix system (42), which is a nonlinear system.The Hessian matrix Λ = K * h K h + αL α (U m ) of the system (42) is extremely large and tends to be ill-conditioned when α is small.This is primarily due to the clustering of eigenvalues of K h around zero [23], while K * h K h is a full matrix, the Fast Fourier transformation (FFT) can be used to evaluate K * h K h U in O(n x logn x ) operations [23] because the blurring kernel exhibits translation-invariant behavior.The advantageous aspect is that the Hessian matrix is symmetric positive definite (SPD).
Therefore, the conjugate gradient (CG) method is suitable for solving the system (42).However, the CG method can have a slow convergence rate for large and ill-conditioned systems.In order to achieve faster convergence, we use the preconditioned conjugate gradient (PCG) method [55][56][57].Here, we introduce our SPD circulant preconditioned matrix P of Strang-type [58].
where Kh is a circulant approximation of matrix K h .The diag(L TV h (U m )) is a diagonal structure of L TV h (U m ).This is summarized in Algorithm 3.

Algorithm 3 THE PCG METHOD function:
[ U ] = PCG(α, U, K, Z) 1. Set: Set: Set: b m = K * h Z, 6. Set: Test: Stopping criteria 8. end While applying PCG to Equation ( 42), we need to take the inverse of the preconditioned matrix P. For the inversion of the first term K * h Kh , we require O(n x logn x ) floating-point operations using FFTs [23].As for the second term in P, which is a diagonal matrix, inversion is not problematic.Now, let the eigenvalues of K * h K h , L α h and L TV h be λ K i , λ α i and λ TV i , respectively, such that λ K i → 0, and λ α i → ∞.Consequently, the eigenvalues of P −1 Λ are given by where λ dTV i are the eigenvalues of diag(L TV h (U m )).It is evident that Therefore, the spectrum of P −1 Λ is more favorable than that of the Hessian matrix Λ.The flowchart of our proposed method is illustrated in Figure 1.

Numerical Experiments
In this section, we use our algorithms to solve the unconstrained TFOV problem (5) and compare the results with the one-sided constrained problems ( 7) and ( 8) and twosided constrained problems (17) and (18).We conduct several sets of experiments using different digital images.The algorithm code is written in MATLAB and all computational experiments are performed on an Intel (R) Core (TM) i7-4510U CPU @2.60 GHz.To evaluate the quality of the restored image, we use the peak signal-to-noise ratio (PSNR) in decibels (dB) [59] and the structural similarity index measure (SSIM) [60].A higher PSNR and SSIM value indicates a high quality of the restored image.The degree to which the algorithms satisfy the constraints is measured by counting the number of pixels with negative values for the one-sided constraint method and the number of pixels outside the range [0, 255] for the two-sided constraint method.To observe the optimum values of our initial parameters c 0 and d, we performed computations on the Barbara image.We observed numerically (see Figure 2) that the optimum ranges for the initial parameters c (0) and d of the augmented Lagrangian method are c (0) ∈ [1 × 10 −5 , 1 × 10 −4 ] and d ∈ [4,10], respectively.For our experiments, we choose c (0) = 0.000085, d = 5.The λ (0) = λ (0) The values of α = 1.8, α = 1 × 10 −8 and β = 0.1 are chosen according to [39].In all experiments, the stopping criterion for the numerical iterations is defined as b − Ax k < tol b , where x k = (v k , u k ) is the solution vector in the k-th iteration.The results are presented in figures and tables.Example 1.Here, we have compared Algorithm 1 (OSCM) with an unconstrained TFOV-based deblurring technique.We used Barbara's image, which presents a challenge due to the combination of a large-scale cartoon element (the face) with a small-scale texture (the shirt).The restored images are shown in Figure 3, with each subfigure having a size of 256 × 256.Table 1 lists the PSNR, SSIM, and the number of negative pixels in the experiment.In this example, the test images are blurred using the PSF given in Figure 4, which is a circular Gaussian kernel.The criterion to stop the computational algorithm is based on a tolerance value of tol = 1 × 10 −7 .

1.
From Figure 3, it can be observed that the deblurred images produced by the OSCM exhibit significantly better quality compared to the unconstrained method.

2.
In Table 1, one can observe that the PSNR and SSIM values of the OSCM are considerably higher than the PSNR and SSIM values of the unconstrained method.The OSCM identifies negative pixels and reduces them as the iterations progress.Finally, it removes them in just 12 iterations, resulting a clear, blur-free image.Example 2. Here, we have compared Algorithm 2 (TSCM) with the unconstrained TFOV-based deblurring method.We used the Moon image, which is a non-texture image.The restored images are shown in Figure 5, each of size 256 × 256.Table 2 lists the PSNR, SSIM, and the number of pixels outside the interval [0, 255].In this example, we have used the PSF given in Figure 4.The criterion to stop the computational algorithm is based on a tolerance of tol = 1 × 10 −7 .

1.
From Figure 5, it can be seen that the deblurred images created by the TSCM are much better than those of the unconstrained method.This can also be verified with the data given in Table 2.

2.
In Table 2, one can observe that the PSNR and SSIM values of TSCM are considerably higher than the PSNR and SSIM values of the unconstrained method.The TSCM identifies pixels that are outside the given range of [0, 255] and modifies them as the number of iterations increases.

Example 3.
Here, we have compared our algorithms (OSCM and TSCM) with the TFOV-based methods by Fairag et al. [39], in which they created the one-level method (OLM) and the two-level method (TLM) for the TFOV-based image deblurring problem.Two different kinds of images were used in this investigation: Goldhills (actual) and Cameraman (complicated).Figures 6 and 7 display the various facets of each picture.Each subfigure has a size of 512 × 512.In this example, test images are blurred by a Gaussian kernel (PSF).Additionally, Gaussian noise with mean µ = 0.01 and σ 2 = 0.5 is added to the images.For the TFOV-based algorithms (OLM and TLM), we used α = 1.8, λ α L = 1 × 10 −16 , and β = 0.1.The Level-II parameter λ J in TLM is calculated according to the formula given in Reference [39].OLM and TLM also use the CG method for numerical solution.The stopping criterion of the computational algorithm is based on a tolerance tol = 1 × 10 −7 .Table 3 contains all the information related to this experiment.Remark 3.

1.
One can see from Figures 6 and 7 that our algorithms (OSCM and TSCM) produce results of slightly higher quality compared to other methods.2.
From Table 3, it can be observed that our algorithms (OSCM and TSCM) consistently achieve higher PSNR and SSIM values compared to other methods (OLM, and TLM) for all photos.
Although the TLM generates higher PSNR and SSIM values more quickly, the quality of the PSNR and SSIM is inferior to that of OSCM and TSCM.Despite taking less time, our algorithms produce significantly better quality compared to other methods.For example, for the Goldhills image, OLM and TLM require 1005.2589 and 526.5476 s, respectively, to achieve PSNR/SSIM values of 33.1589/0.7704and 33.1458/0.7690,respectively.However, OSCM and TSCM take 896.4058 and 909.5469 s, respectively, to achieve higher PSNR/SSIM values of 34.8945/0.7788and 34.8965/0.7759,respectively.Similar trends can be observed in the Cameraman's image.Therefore, our algorithms (OSCM and TSCM) produce high-quality deblurred images compared to other methods.Example 4. Here, we have also compared our algorithms (OSCM and TSCM) with the TFOV-based methods (OLM and TLM) proposed by Fairag et al. [39].We used a Pepper image (a non-texture image) for this comparison.Figure 8 displays various facets of the given image, with each subfigure having a size of 512 × 512.In this example, the test images are blurred by a motion kernel (PSF) with a motion length of l = 256 and an angle of motion θ = 15 0 .Salt and pepper noise with a small density of 0.01 is added to the blurry image.For the TFOV-based algorithms (OLM and TLM), we used α = 1.8, λ α L = 1 × 10 −12 , and β = 0.1.The Level-II parameter λ J in TLM is calculated according to the formula given in [39].The stopping criterion of the computational algorithm is based on a tolerance of tol = 1 × 10 −6 .Table 4 contains all the information related to this experiment.Remark 4.

1.
One can see from Figure 8 that our algorithms (OSCM and TSCM) produce results of slightly higher quality compared to other methods.2.
From Table 4, it can be observed that for all photos, our algorithms (OSCM and TSCM) exhibit higher PSNR values compared to other methods (OLM and TLM).Although the TLM generates faster PSNR and SSIM computation, its quality is inferior to that of OSCM and TSCM.Therefore, our algorithms (OSCM and TSCM) produce superior-quality deblurred images when compared to other methods.Example 5.In this example, we used two satellite images from Reference [61].The blurred images were corrupted by Poisson noise and blurring artifacts.The blurring process was conducted using the f special( gaussian , 9, sqrt(3)) kernel.Dealing with the addition of Poisson in the images proves to be particularly challenging for most deblurring methods.Imaging modalities like this often exhibit the presence of Poisson noise, primarily arising from photon counting.Simultaneously, blurring is an inevitable outcome resulting from the physical mechanism of an imaging system, accurately represented as the convolution of the image with a point spread function.For comparison, we used the non-blind fractional order total variation-based method (NFOV) [61] and Richardson-Lucy algorithm with total variation regularization (RLTV) [62].The restored images of Galaxy are shown in Figure 9, each with a size of 256 × 256.The restored images of Satel are shown in Figure 10, each with a size of 128 × 128.For the NFOV and RLTV methods, the parameters used are according to [61].Table 5 lists the information about this experiment.6, it is evident that our methods, OSCM and TSCM, consistently produce better values for PSNR and SSIM when compared to all other methods.These results demonstrate the strong performance of the OSCM and TSCM, which consistently generate high-quality images.A comparison of PSNR and SSIM values for different methods using Levin's dataset is depicted in Figure

Conclusions
For the image deblurring problem, we presented OSCM and TSCM using a TFOV regularization functional.In addition to guaranteeing strictly positive outcomes, both OSCM and TSCM impose upper limitations on image intensity levels, maintaining them within a predetermined range.We applied our proposed approaches to conduct numerical tests on various image types, including synthetic, real, complex, satellite and non-texture images.To evaluate our algorithms, we also used images from Levin's dataset [63].We compared the OSCM and TSCM with the most recent TFOV-based approaches, mainly TV, OLM, TLM, RLTV, and NFOV.The numerical tests demonstrated the efficiency of our proposed techniques.In the future, we plan to develop the OSCM and TSCM for the other computationally expensive regularization functionals, such as mean curvature functional will be pursued in the future.Additionally, we aim to design a constrained model within a similar framework for the blind image deblurring problem.It is worth noting that under specific conditions, the proposed techniques can be implemented in other image processing models.Then, it suffices to take ν ∈ C 1 0 (Ω, R) (the space of first-order differentiable functions vanishes at the boundary), which implies Hence, Equation (A2) reduces to Equation (30).Case-II: If ν ∈ W α 1 (Ω), then Therefore, the boundary terms in Equation (A3) can only become zero if This concludes the proof.

Figure 2 .
Figure 2. PSNR for Barbara image against different values of parameters.

Figure 3 .
Figure 3. (a) is an exact image.(b) is a blurry image.(c) is a deblurred image by the unconstrained method.The images from (d-o) are deblurred by the OSCM according to iterations from k = 1 to k = 12, respectively.

Figure 5 .
Figure 5. (a) is an exact image.(b) is a blurry image.(c) is a deblurred image by the unconstrained method.The images from (d-s) are deblurred by TSCM according to iterations k = 1 to k = 16, respectively.

Figure 6 .
Figure 6.(a) is an exact image, (b) is a blurry image, (c) is an image deblurred by OLM, (d) is an image deblurred by TLM (e), deblurred image by OSCM and (f) is a deblurred image by TSCM.

Figure 7 .
Figure 7. (a) is an exact image, (b) is a blurry image, (c) is an image deblurred by OLM, (d) is a an imaged deblurred by TLM (e) deblurred image by OSCM and (f) is a deblurred image by TSCM.

Figure 8 .
Figure 8.(a) is an exact image, (b) is a blurry image, (c) is a deblurred image by OLM, (d) is an image deblurred by image TLM, (e) is an image deblurred image by OSCM, and (f) is an image deblurred by TSCM.

Figure 9 .
Figure 9. Galaxy image: (a) The blurry image, (b) deblurred image by the RLTV method, (c) deblurred image by the NFOV method, (d) deblurred image by the OSCM and (e) deblurred image by the TSCM.

Figure 11 .
Figure 11.Image 1: (a) Exact image, (b) blurry image, (c) image deblurred by the TV method, (d) image deblurred by the OLM, (e) image deblurred by the TLM, (f) image deblurred by the RLTV method, (g) image deblurred by the NFOV method, (h) image deblurred by the OSCM and (i) deblurred image by the TSCM.

Figure 12 .
Figure 12.Image 2: (a) Exact image (b) blurry image, (c) deblurred image by the TV method, (d) deblurred image by the OLM, (e) deblurred image by the TLM, (f) deblurred image by the RLTV method, (g) deblurred image by the NFOV method, (h) deblurred image by OSCM and (i) deblurred image by the TSCM.

Figure 13 .
Figure 13.Image 3: (a) Exact image, (b) blurry image, (c) image deblurred by the TV method, (d) image deblurred by the OLM method, (e) image deblurred by the TLM, (f) image deblurred by the RLTV method, (g) image deblurred by the NFOV method, (h) image deblurred by the OSCM method, and (i) deblurred image by the TSCM.

Figure 14 .
Figure 14.Image 4: (a) Exact image, (b) blurry image, (c) image deblurred by the TV method, (d) image deblurred by the OLM method, (e) image deblurred by the TLM, (f) image deblurred by the RLTV method, (g) image deblurred by the NFOV method, (h) image deblurred by the OSCM method, and (i) deblurred image by the TSCM.

Figure 15 .
Figure 15.Comparison of PSNR and SSIM values for different methods using Levin's dataset.

Table 1 .
The PSNR, SSIM, and number of negative pixels for Example 1.

Table 3 .
PSNR, SSIM, and CPU-Time comparison of different methods for Example 3.

Table 4 .
PSNR, SSIM and CPU-Time Comparison of different methods for Example 4.

Table 5 .
PSNR and SSIM comparison of different methods for Example 5.

Table 5 ,
one can observe that all of the methods generated nearly identical results.However, OSCM and TSCM better PSNR and SSIM values compared to all other methods.This demonstrates the effectiveness of the OSCM and TSCM in generating high-quality images.
[63]ple 6.In this example, we used four different images from the dataset of Levin et al.[63].The ke g en(N, 100, 5) kernel was used for blurring.For comparison, we used TV, OLM, TLM, RLTV, NFOV, OSCM, and TSCM.Restored images are shown in Figures11-14.Each one is of the size 255 × 255.Table6lists the information of this experiment.

Table 6 .
PSNR and SSIM comparison of different methods for Example 6.