In this section, numerical examples will be provided to demonstrate the effectiveness of the proposed model and algorithm. All experiments were conducted on a Lenovo laptop equipped with an Intel 2.5 GHz CPU and 4 GB of memory.
6.2. Experimental Results
In
Figure 2,
Figure 3,
Figure 4 and
Figure 5, we compared the proposed method with several classic competing methods, FastTVMM [
2] and WFBM [
1].
Figure 2 and
Figure 3 display the recovery of the Lena image, previously impaired by the implementation of Scenario I and superimposed with Gaussian noise characterized by zero mean and a standard deviation of 10
−3. In
Figure 4 and
Figure 5, our restoration work encompassed the Chimpanzee image, degraded as a result of Scenario II and the imposition of Gaussian noise characterized by zero mean and a standard deviation of 10
−3. For example, the eyes and eyelashes of Lena in
Figure 2 and
Figure 3, or the hair on the face of the chimpanzee in
Figure 4 and
Figure 5. Its recovery results often exhibit an oversmooth phenomenon, and readers can observe Lena’s eye corner area in
Figure 2 and
Figure 3 (the face area of the circle in
Figure 4 and
Figure 5).
FastTVMM [
2] is widely recognized as a traditional and effective method for addressing image blurring issues. It performs admirably in the restoration of smooth regions, achieving high precision, but its performance in preserving image details is not satisfactory when restoring images. WFBM [
1] is a new image deblurring model based on wavelet frames, which explicitly treats images as piecewise smoothing functions. Contrasting with existing methods, our novel approach boasts a notable advantage in precisely restoring subtle edge details, such as the complex texture of feathers on hats and the minute strands of eyelashes, as evidenced in
Figure 2 and
Figure 3. To put it differently, our methodology avoids the common pitfall of producing artifacts in areas containing edges, setting it apart from competing approaches.
For Lena images (see
Figure 2 and
Figure 3), the proposed method is significantly superior to the other two methods. Contrasting with existing methods, our novel approach boasts a notable advantage in precisely restoring subtle edge details, such as the complex texture of feathers on hats and the minute strands of eyelashes, as evidenced in
Figure 2 and
Figure 3. For Chimpanzee images (see
Figure 4 and
Figure 5), the subjective quality of the deblurring results obtained by the proposed method is significantly better than other competing schemes. The fundamental reason for this phenomenon is that the integer derivative operator is a local operator that only considers the relationship between pixels in the image and their neighboring pixels while ignoring the correlation of global information in the image. The texture features of images often exhibit obvious non-locality and self-similarity. Therefore, as a non-local extension of integer order calculus, the proposed fractional order model has a good texture preservation function. Precisely, the non-local feature of the fractional-order derivative in the model is instrumental in safeguarding texture integrity.
Table 2 provides the parallel comparison results of the corresponding PSNR and SSIM values in
Figure 2,
Figure 3,
Figure 4 and
Figure 5. As evidenced by the data in
Table 2, the proposed methodology shows strong capabilities, particularly in the aspects measured by PSNR and SSIM. From the last row of
Table 2, it can be seen that the proposed method has the highest average values in both PSNR and SSIM metrics. Specifically, the proposed method improved PSNR values by 0.9733 (i.e., 0.9733 is obtained by subtracting 32.3494 from 31.3761) and SSIM values by 0.0111 (i.e., 0.0111 is obtained by subtracting 0.8752 from 0.8641) compared to the FastTVMM method.
The effectiveness of the proposed model in addressing noise and blur image restoration is evaluated and compared in
Figure 6 and
Figure 7.
Figure 6a is damaged by Scenario III and zero-mean Gaussian noise with a standard deviation of 2.
Figure 7a is damaged by Scenario IV and zero-mean Gaussian noise with a standard deviation of 2. As shown in
Figure 6b and
Figure 7b, despite the impressive denoising prowess displayed by the NNTGV [
3] technique, a substantial amount of detail is sacrificed, prominently seen in the absence of intricate features, including the house located farthest to the right in the visual frame, and the restoration results show an oversmooth phenomenon (such as the lawn below the image being almost smoothed out). Compared with the NNTGV [
3] method, although the QDSV [
4] method has made significant improvements in detail preservation, like the NNTGV [
3] method, the QDSV [
4] method still lacks satisfactory deblurring results. Readers wishing for a more detailed analysis are invited to closely inspect the representation of the grassland as shown in
Figure 6c and
Figure 7c. In contrast, although the proposed method generates some residual noise in the results, it has clearer visual results. For more details, readers can observe the grassland in
Figure 6d and
Figure 7d. The reason for the presence of noise in
Figure 6d and
Figure 7d is that the proposed model requires both denoising and deblurring, and the blurring operator is ill-conditioned. In general, when dealing with simultaneous deblurring and denoising tasks, the images requiring restoration are commonly impaired by the presence of Gaussian noise with a mean value of zero and a standard deviation set at 0.001. Therefore, the Gaussian noise interfering with
Figure 6a and
Figure 7a is very large. It is for this reason that the results illustrated in
Figure 6d and
Figure 7d exhibit a notable presence of noise.
Overall, the NNTGV [
3] method has impressive denoising performance, but suffers from severe loss of details. Contrasted with the NNTGV [
3] technique, the QDSV [
4] approach has indeed made strides in preserving intricate details; however, much like NNTGV [
3], QDSV [
4] remains deficient in achieving commendable deblurring effects. In comparison, our method can achieve a compromise between denoising and preserving details. That is, in the presence of very small noise, images that can or have rich details.
Figure 8 shows the convergence behavior of the iterative support shrinking algorithm. Specifically,
Figure 8a shows the history of the energy function value
changing with the number of outer iteration number;
Figure 8b shows the history of the
value changing with the number of outer iteration number. We can clearly observe that as the outer iteration number
k increases,
decreases and converges, which validates the theoretical result in Lemma 4. Meanwhile,
converging to 0 corroborates the theoretical finding in Lemma 5.
To demonstrate that the proposed model has good restoration results for images contaminated by different deblurring scenes, we present the restoration results of the proposed method for images contaminated by Scenario V and Scenario VI in
Figure 9. Based on the results of restoring images contaminated by different scenes, we found that the proposed method applies to a wide range of images and blur kernels. In fact, the proposed method has extensive utility across various fields in the real world. Specifically, the proposed method not only performs well in natural image (see
Figure 2,
Figure 3,
Figure 4,
Figure 5,
Figure 6 and
Figure 7) restoration, but also performs well in restoring piecewise constant images (see
Figure 9). In
Figure 10,
Figure 11,
Figure 12,
Figure 13,
Figure 14,
Figure 15,
Figure 16 and
Figure 17, we show the restoration results of our method on two test images and six blurring kernels. From
Figure 10,
Figure 11,
Figure 12,
Figure 13,
Figure 14,
Figure 15,
Figure 16 and
Figure 17, it further verified that our method is suitable for a wide range of images and blurred kernels. This is because the variable exponential regularization term in the model allows our scheme to be completely independent of any assumptions about the spatial distribution of the kernel. However, the proposed method still faces some challenges, namely, its effectiveness is not ideal for situations with high noise levels. The reason lies in the dual requirement of the proposed model to address both noise and blur, making the blurring operator inherently difficult to manage. Moreover, we also provide the values of SNR and RMSE at the annotations in
Figure 10,
Figure 11,
Figure 12,
Figure 13,
Figure 14,
Figure 15,
Figure 16 and
Figure 17. We found that most of our results have good results in terms of SNR and RMSE values. This further validates the advantages of the proposed method.
In addition, we also compared our methods with those of the past two years. The experimental results are displayed in
Figure 18 and
Figure 19. From
Figure 18c and
Figure 19c and the corresponding enlarged images, it can be clearly seen that the TV
p,q(x) method produces results with significant loss of texture details. The visual quality of the deblurring results obtained using D-TGV methods is too smooth and accompanied by obvious artifacts. The specific results can be seen in
Figure 18d and
Figure 19d. In contrast, the proposed method achieves higher visual quality, including preserving details, preventing oversmoothing, and avoiding artifacts.
To better demonstrate the effectiveness of the proposed method, we provide a computational complexity metrics. The specific values of each metrics are recorded in
Table 3. From
Table 3, we can clearly see that the single iteration time of the three images with a size of
(i.e.,
Figure 10a–
Figure 12a) exceeds 10 s, the single iteration time of the three images with a size of
(i.e.,
Figure 13a and
Figure 14a and
Figure 17a) exceeds 1 s, and the single iteration time of the two smaller images (i.e.,
Figure 15a and
Figure 16a) is less than 1 s. Similarly, the total iteration time and total number of iterations follow the same pattern. In addition, we have provided the time required to obtain restoration results using the proposed method and other methods. The data we tested is the time required for the restoration of
Figure 10a–
Figure 17a. The specific numerical results are recorded in
Table 4 below. From
Table 4, we clearly find that, except for the FastTVMM method, the proposed method requires the least amount of time to restore blurry images. The FastTVVM method requires less time because it uses an alternating iteration algorithm. However, the restoration results of the FastTVMM method exhibit oversmoothness, resulting in the loss of texture details. This is not the result we want to obtain.
Finally, we also present the deblurring results obtained using the proposed method with different fractional-orders in
Figure 20. From
Figure 20, we can clearly see that the proposed method maintains different levels of detail in deblurring images using different fractional-orders. The reason for this phenomenon is determined by the internal structure of the image itself.
6.3. Discussion
In this paper, we propose a novel fractional-order non-convex TVα,p model in image deblurring. The fractional-order derivatives we introduce have non-locality and self-similarity, which can effectively simulate the texture information of images. Therefore, it can produce results that preserve the texture of the image. On the contrary, integer-order derivative operators are local operators that only consider the relationship between pixels in the image and their neighboring pixels, while ignoring the correlation of global information in the image. Therefore, the method of modeling based on integer-order gradient generally produces image processing results with oversmoothing, resulting in the inability to restore texture and other details in the image. Compared with several classic methods and several latest methods, the experimental results have verified the effectiveness of the proposed method.
However, there are also some drawbacks. Firstly, the proposed model should not handle blurry images with too much noise. If the blurred image contains too much noise, the proposed model’s deblurring results will also have significant noise. This is because deblurring is essentially an inverse problem that involves solving a ill-conditioned equation. This inverse problem is numerically unstable, meaning that even small errors can be significantly amplified. In addition, in the presence of a large amount of noise, deblurring operations may amplify these noises, resulting in noticeable noise traces in the restored image. In order to more intuitively demonstrate the impact of different noise levels on the deblurring effect of the proposed method, we conducted a set of experiments in
Figure 21. From
Figure 21, we can clearly see that when the noise level
, the proposed method produces unsatisfactory results in deblurring. Secondly, the proposed model is experimented on specific types of datasets (i.e., natural scenery), which may limit its direct application in other fields (e.g., medical imaging, satellite imagery). Moreover, for images under extreme conditions (such as extreme low light, high dynamic range scenes), the recovery results generated by the proposed model are not ideal. In order to more intuitively demonstrate the impact of extreme conditions (taking darkness levels as an example) on the deblurring effect of the proposed method, we conducted a set of experiments in
Figure 22. From
Figure 22, we can clearly see that when the darkness level exceeds 3 (including when the darkness level is equal to 3), the proposed method produces unsatisfactory results in deblurring.
It is worth mentioning that the proposed model also has a wide range of applications, including virtual reality and augmented reality (VR/AR), nature conservation and documentation, artistic and commercial uses, and wildlife photography. High-quality, deblurred nature images can be used to create more immersive VR/AR experiences, such as virtual tours of natural landscapes, which can be beneficial for education and tourism. By providing clearer and more detailed images, our model can support the documentation of biodiversity, assist in species identification, and help in creating high-resolution databases for ecological research. Deblurred nature images can also find applications in the creation of high-quality prints, digital art, and commercial photography, enhancing the visual appeal and marketability of nature-related media. The proposed method can improve the clarity of images taken in natural settings, where motion blur is common due to the unpredictable movements of animals. This can aid in wildlife monitoring, research, and conservation efforts.
In our future work, we plan to investigate adaptive fractional-order models. Additionally, we will also study the characteristics of images under extreme conditions (such as extreme low light, high dynamic range scenes) and design reasonable methods to handle different blurs and noise. Moreover, it is essential to develop more efficient algorithms to enhance the computational performance of the model.