Decomposition-Guided Curvelet Thresholding for Sharp-to-Soft CT Kernel Conversion
Authors: Mahmoud Nasr, Jan K. Argasinski, Krzysztof Brzostowski, Adam Piorkowski
Organizations: Department of Biocybernetics and Biomedical Engineering, AGH University of Krakow, 30-059 Krakow, Poland. · Sano Centre for Computational Medicine, Czarnowiejska 36/C5, Kraków, 30-054, Poland. · Faculty of Physics, Astronomy and Applied Computer Science, Jagiellonian University, Łojasiewicza 11, Kraków, 30-348, Poland. · Faculty of Information and Communication Technology, Wroclaw University of Science and Technology, Wyb. Wyspiańskiego 27, 50-370 Wroclaw, Poland.
Image denoising is a crucial task in image processing, focused on improving image quality by minimizing noise while maintaining essential structural elements. This study presents a hybrid denoising framework that combines several decomposition techniques, including empirical mode decomposition (EMD), variational mode decomposition (VMD), multichannel EMD (MEMD), and bidimensional EMD (BEMD), with curvelet transform thresholding. Each decomposition mode undergoes processing through both soft and hard thresholding, and the denoised modes are combined to rebuild the final image. Comprehensive evaluations of standard CT image datasets reconstructed with various kernels (B50, B46, B41, B36) reveal substantial enhancements in denoising efficacy. VMD consistently achieves the highest peak signal-to-noise ratio (PSNR) and structural similarity index (SSIM), signifying exceptional noise reduction and feature preservation. The study analyses the trade-offs between soft and hard thresholding: soft thresholding maintains intricate visual details, whilst harsh thresholding provides enhanced noise reduction. The suggested method surpasses traditional techniques in both reference and non-reference quality criteria, indicating its potential for broader application in medical imaging and future incorporation with adaptive thresholding algorithms.
Figures & tables
Figure 1: This figure shows the structure of the proposed method based on the 3 steps for image denoising and enhancing image quality.
Figure 2: Samples of the denoised images of an input reconstructed image B50 (sharp Kernel) by applying the soft threshold.
Metric
Curvelet
EMD + Curvelet
VMD + Curvelet
MEMD + Curvelet
BEMD + Curvelet
SSIM
0.931
0.903
0.949
0.954
0.949
PSNR
32.344
32.548
33.894
33.094
32.816
SNR
19.079
19.283
20.629
19.829
19.551
NIQE Score
4.656
9.543
5.396
5.040
5.056
BRISQUE Score
40.575
42.032
44.369
42.191
42.331
PIQE Score
65.059
66.000
61.523
75.042
70.478
Table 1: Performance Metrics for Soft Thresholding (B50)
Metric
Curvelet
EMD + Curvelet
VMD + Curvelet
MEMD + Curvelet
BEMD + Curvelet
SSIM
0.949
0.915
0.954
0.957
0.953
PSNR
37.856
36.516
38.475
38.187
37.796
SNR
24.523
23.183
25.143
24.854
24.463
NIQE Score
4.686
9.888
5.061
5.294
5.340
BRISQUE Score
42.536
41.963
40.965
43.431
43.758
PIQE Score
71.384
71.945
73.732
76.598
71.873
Table 2: Performance Metrics for Soft Thresholding (B46)
Metric
Curvelet
EMD + Curvelet
VMD + Curvelet
MEMD + Curvelet
BEMD + Curvelet
SSIM
0.949
0.913
0.957
0.958
0.953
PSNR
38.515
37.009
39.528
39.358
38.803
SNR
26.551
25.044
27.564
27.394
26.839
NIQE Score
4.937
11.547
5.510
5.515
5.537
BRISQUE Score
43.902
42.462
43.255
43.603
43.923
PIQE Score
78.073
75.470
77.409
80.670
78.187
Table 3: Performance Metrics for Soft Thresholding (B41)
Metric
Curvelet
EMD + Curvelet
VMD + Curvelet
MEMD + Curvelet
BEMD + Curvelet
SSIM
0.964
0.915
0.957
0.959
0.958
PSNR
43.112
38.321
41.343
41.700
41.550
SNR
30.687
25.896
28.917
29.274
29.124
NIQE Score
5.018
11.728
5.171
5.322
5.226
BRISQUE Score
42.755
42.630
44.310
44.577
44.366
PIQE Score
78.889
75.834
80.792
83.809
82.571
Table 4: Performance Metrics for Soft Thresholding (B36)
Figure 3: Samples of the denoised images of an input reconstructed image B50 (sharp Kernel) by applying the hard threshold.
Metric
Curvelet
EMD + Curvelet
VMD + Curvelet
MEMD + Curvelet
BEMD + Curvelet
SSIM
0.945
0.900
0.940
0.952
0.946
PSNR
35.905
35.650
37.429
35.909
35.731
SNR
22.640
22.385
24.164
22.643
22.466
NIQE Score
4.486
7.464
5.509
5.117
5.057
BRISQUE Score
42.684
36.714
39.820
40.952
44.160
PIQE Score
46.854
26.838
36.184
48.065
49.961
Table 5: Performance Metrics for Hard Thresholding (B50)
Metric
Curvelet
EMD + Curvelet
VMD + Curvelet
MEMD + Curvelet
BEMD + Curvelet
SSIM
0.971
0.941
0.972
0.972
0.970
PSNR
42.329
40.082
42.587
41.754
41.914
SNR
28.997
26.750
29.254
28.421
28.581
NIQE Score
4.708
7.997
4.785
5.184
5.056
BRISQUE Score
40.812
38.801
45.586
39.776
40.756
PIQE Score
59.926
38.451
52.565
65.817
62.978
Table 6: Performance Metrics for Hard Thresholding (B46)
Metric
Curvelet
EMD + Curvelet
VMD + Curvelet
MEMD + Curvelet
BEMD + Curvelet
SSIM
0.975
0.949
0.978
0.977
0.977
PSNR
42.930
40.783
43.812
43.178
43.184
SNR
30.966
28.819
31.848
31.214
31.220
NIQE Score
5.114
8.948
5.128
5.582
5.387
BRISQUE Score
45.910
41.408
47.861
43.284
44.929
PIQE Score
66.931
46.319
53.507
75.715
69.959
Table 7: Performance Metrics for Hard Thresholding (B41)
Metric
Curvelet
EMD + Curvelet
VMD + Curvelet
MEMD + Curvelet
BEMD + Curvelet
SSIM
0.979
0.947
0.974
0.974
0.975
PSNR
48.418
42.252
45.942
46.410
46.836
SNR
35.992
29.827
33.516
33.984
34.410
NIQE Score
5.308
9.469
5.020
5.348
5.172
BRISQUE Score
40.779
40.778
40.380
42.512
41.557
PIQE Score
67.228
51.531
69.299
77.144
76.257
Table 8: Performance Metrics for Hard Thresholding (B36)
The quality of computed tomography (CT) images is significantly affected by the selection of reconstruction kernels: sharp kernels improve spatial resolution but increase noise, whereas soft kernels diminish noise at the expense of edge clarity. This study presents an innovative enhancement framework utilising Bidimensional Empirical Mode Decomposition in conjunction with Quaternion Bilateral Filtering (BEMD--QBF) to convert sharp-kernel CT images into representations resembling soft-kernels, while maintaining critical anatomical structures. The technique disaggregates each image into intrinsic mode functions via BEMD and analyzes them inside a cohesive quaternion framework to attain efficient noise reduction and structural integrity. The proposed methodology is evaluated using several reconstruction kernels (B50, B46, B41, B36, B35, B31) and compared with recognised filtering strategies, including Non-Local Means, Anisotropic Diffusion, Bilateral Filtering, and Quaternion Bilateral Filtering. Quantitative evaluations of the Structural Similarity Index (SSIM) and Peak Signal-to-Noise Ratio (PSNR) indicate that BEMD-QBF consistently attains superior structural fidelity and competitive noise reduction across all evaluated kernels. The results underscore the efficacy of the proposed strategy as a viable approach to enhancing post-reconstruction CT images, yielding superior image quality without requiring access to raw projection data.
Mahmoud Nasr, Jan K. Argasinski, Krzysztof Brzostowski +1
Sano Centre for Computational Medicine, Czarnowiejska 36/C5, Kraków, 30-054, Poland · Department of Biocybernetics and Biomedical Engineering, AGH University of Krakow, 30-059 Krakow, Poland · Faculty of Physics, Astronomy and Applied Computer Science, Jagiellonian University, Łojasiewicza 11, Kraków, 30-348, Poland +1
Computed tomography (CT) plays a crucial role in medical diagnosis, but minimizing radiation exposure while maintaining image quality remains a critical challenge. Low-dose CT (LDCT) protocols reduce radiation risks but inevitably suffer from severe noise and artifacts that compromise diagnostic accuracy. While existing deep learning methods have achieved promising results, there remains a continuous quest for generative paradigms that intrinsically capture global-to-local structural dependencies to better preserve fine anatomical details. To this end, we propose DeVAR, a novel generative framework that applies visual autoregressive modeling (VAR) to LDCT denoising for the first time. Conditioned on global context provided by LDCT prefix tokens, DeVAR progressively generates discrete token maps of the target normal-dose CT (NDCT) via next-scale prediction. Because quantization inherently discards high-frequency information, we introduce a residual refiner to capture subtle anatomical structures beyond the capacity of a discrete codebook. Finally, empowered by a dual-representation hybrid training strategy, our hybrid NDCT decoder seamlessly integrates continuous and discrete latents to reconstruct high-fidelity, detail-preserved images. Extensive experiments on two public datasets demonstrate that DeVAR consistently achieves superior qualitative and quantitative performance compared to state-of-the-art LDCT denoising methods.
Xizhuo Zhang, Yannian Gu, Zhongzhen Huang +2
Shanghai Jiao Tong University, China · SenseTime Research, China · Shanghai Innovation Institute, China
The total variation (TV) method is an image denoising technique that aims to reduce noise by minimizing the total variation of the image, which measures the variation in pixel intensities. The TV method has been widely applied in image processing and computer vision for its ability to preserve edges and enhance image quality. In this paper, we propose a Mixed-norm TV (MixTV) model for image denoising and the associated numerical algorithm to carry out the procedure, which is particularly effective in removing several types of noise and their combinations. Our MixTV admits a unique solution and the associated numerical algorithm guarantees convergence. Numerical experiments are demonstrated to show improved effectiveness and denoising quality compared to other TV models. Such encouraging results further enhance the utility of the TV method in image processing. Our project page is available at https://jing-en-huang.github.io/MixTV.
Jing-En Huang, Jia-Wei Liao, Ku-Te Lin +2
National Yang Ming Chiao Tung University · National Taiwan University · National Central University +1