Comparative Study of Weighted and Coupled Second- and Fourth-Order PDEs for Image Despeckling in Grayscale, Color, SAR, and Ultrasound

Abstract: Partial Differential Equation (PDE)-based approaches have gained significant attention in image despeckling due to their strong capability to preserve structural details while suppressing noise. However, conventional second-order PDE models tend to generate blocky artifacts, whereas higher-order models often introduce speckle patterns. To resolve it, this paper proposes and comparatively analyzes two advanced PDE-based frameworks designed for speckle noise suppression while preserving the fine edges. The first model introduces a novel weighted formulation that combines second and fourth-order PDEs through a weighting parameter. The second-order diffusion coefficient employs grayscale and gradient-based indicators, while the fourth-order term is guided solely by a Laplacian-based indicator. The second model constructs a coupled PDE framework, where independent fourth and second-order components are explicitly solved in an iterative manner. In this coupled structure, each diffusion coefficient is defined separately to enhance adaptability in varying image regions. Both models are implemented using the explicit finite difference method. The proposed techniques are extensively evaluated on a variety of datasets, including standard grayscale, color, Synthetic Aperture Radar (SAR), and ultrasound images. Comparative experiments with the existing Telegraph Diffusion Model (TDM) and Fourth-Order Telegraph Diffusion Model (TDFM) demonstrate the superiority of the proposed approaches in reducing speckle noise while effectively preserving fine image structures and edges. Quantitative evaluations using PSNR, SSIM and Speckle Index metrics confirm that the proposed models produce higher image quality and enhanced visual perception. Overall, the presented PDE-based formulations provide a reliable and efficient framework for image despeckling in both natural and medical imaging.

• The paper introduces two novel PDE-based frameworks that balance noise suppression and structure preservation using weighted and coupled second-/fourth-order approaches.
• It employs explicit finite-difference schemes and adaptive parameter tuning to optimize PSNR, MSSIM, and speckle index across varied imaging modalities.
• Results show that Model 1 excels at preserving fine image structures while Model 2 achieves superior intensity recovery, supporting robust real-time despeckling.

Comparative Analysis of Weighted and Coupled Second- and Fourth-Order PDEs for Image Despeckling

Introduction and Problem Context

Speckle noise is a fundamental challenge in imaging systems employing coherent modalities, including synthetic aperture radar (SAR) and medical ultrasound. Characterized by its multiplicative nature, speckle substantially degrades visual interpretability and quantitative analyses. Traditional spatial domain techniques, transform-based approaches, and recent advances in deep learning have seen broad application; however, each suffers from trade-offs relating to structure preservation, artifact introduction, computational burden, or robustness to unseen noise statistics.

Partial Differential Equation (PDE)-based methods, especially those leveraging anisotropic and adaptive diffusion, have demonstrated compelling structural preservation properties. Yet, second-order PDEs (notably Perona-Malik anisotropic diffusion) are prone to staircase artifacts, while fourth-order models can induce unnatural smoothing or residual speckle structures, especially on strong edges or in highly textured regions.

This paper presents a comparative study of two advanced PDE-based frameworks: a weighted combination of second- and fourth-order PDEs (Model 1) and a coupled iterative approach solving second- and fourth-order terms independently (Model 2). Both aim to optimize the trade-off between effective despeckling and fidelity to fine image structures across grayscale, color, SAR, and ultrasound datasets.

PDE-Based Despeckling Model Formulations

The proposed denoising architectures are rooted in the hyperbolic telegraph PDE structure, augmenting classical parabolic diffusion with second-order temporal derivatives and customizable damping, thus enabling a more flexible balance of spatial and temporal restoration dynamics.

Model 1 introduces a weighted sum of the second- and fourth-order diffusion operators. The principal innovation is the weighting parameter $A$ ($0 \leq A \leq 1$), which directly governs the relative impact of edge-aware (second-order) and texture-preserving (fourth-order) smoothing:

\begin{equation} \frac{\partial^2 I}{\partial t^2} + \gamma \frac{\partial I}{\partial t} = A \nabla \cdot [C_1(I_\xi) C_3(|\nabla I_\xi|) \nabla I]

• (1-A)\Delta[C_2(|\Delta I_\xi|)\Delta I]
• \lambda \left(\frac{I-f}{I}\right)^2 \end{equation}

Here, $C_1$, $C_2$, and $C_3$ are locally adaptive diffusion coefficients constructed from smoothed intensity and gradient/Laplacian statistics.

Model 2 instead deploys a coupled system, iteratively applying second- and fourth-order PDEs with independently optimized coefficient functions for each diffusion order. This design offers greater adaptability to heterogeneous image regions by decoupling the regularization pathways.

Both models leverage explicit finite-difference discretizations in space and time, offering efficient and transparent numerical implementation. The color extension is achieved by applying the framework channel-wise on RGB decompositions and aggregating results post-denoising.

Parameter Selection and Optimization

The weighting parameter $A$ in Model 1 is empirically tuned ($A=0.7$ yields optimal results across datasets) to provide maximal noise suppression consistent with minimal blurring (see Figure 1).

Figure 1: Effect of varying the weighting parameter $A$ on image restoration results.

A systematic grid search is carried out for all core parameters, including $\gamma$ (damping), $\lambda$ (fidelity), $0 \leq A \leq 1$0, and $0 \leq A \leq 1$1, with best-fit settings specific to modality and noise severity. This ensures calibration against the distinctive statistical profiles typical of natural, SAR, and ultrasound images.

Quantitative and Qualitative Evaluation

Performance is evaluated via PSNR (reconstruction fidelity), MSSIM (perceptual structure), and speckle index (S.I.; relevant to SAR/ultrasound), all computed on both synthetic and real-world datasets.

Representative input images are shown in Figures 1 and 2.

Figure 2: Images: (a) Parrots, (b) Texture.

Figure 3: Color images: (a) caps, (b) baboon.

Comprehensive restoration results across a wide range of noise intensities (parameterized by the number of looks $0 \leq A \leq 1$2) are detailed in Figures 4–7 for the typical grayscale and color datasets. These illustrate the progressive performance gains from TDM and TDFM (prior state-of-the-art telegraph-based PDEs) to the proposed Model 1 and Model 2 at each noise level.

Figure 4: Look=1.

Figure 5: Look=1.

Figure 6: Look=1.

Figure 7: Look=1.

The SAR and ultrasound case studies further demonstrate the relevance in practical, reference-free settings (see Figure 8).

Figure 8: Comparison of SAR Image Restoration. Each row: (1) Noisy Image, (2) Model 1 Restored, (3) Model 2 Restored.

The metrics are synthesized in Figures 10 and 11, showing PSNR and MSSIM trends for all methods across increasing noise, delineating the regimes in which each model excels.

Figure 9: Comparison of PSNR and MSSIM for Proposed Model 1, 2, and state-of-art TDM and TDFM for the gray parrots image.

Figure 10: Comparison of PSNR and MSSIM for Proposed Model 1, 2, and state-of-art TDM and TDFM for the color caps image.

Model 1 consistently exhibits the highest MSSIM in grayscale (parrots/textures), especially at low and medium noise, reflecting optimal preservation of fine image structures, while Model 2 dominates in PSNR, corresponding to best squared-error intensity recovery. In the context of complex, richly textured images (baboon), the superiority of Model 1 is magnified across both metrics and all noise intensities.

Notably, S.I. reduction for SAR and ultrasound underscores robust speckle suppression for both new schemes, with Model 1 achieving the minimal residual index.

These findings are complemented by ROI visualizations (see figure below) contrasting Model 1/2 with TDFM and confirming improved edge and detail preservation by both new formulations.

Figure 11: TDFM.

Implications, Theoretical and Practical

These results have far-reaching implications for both foundational image restoration theory and real-world imaging pipelines. The explicit coupling of second and fourth-order regularization in both weighted and decoupled forms provides a tunable pathway to optimize the bias-variance frontier for different application domains and image characteristics. The models’ reliance on interpretable and physically motivated regularization contrasts with black-box deep learning, offering a transparent approach favored in critical contexts (medical, remote sensing) where model trust is essential.

Practically, the efficient explicit scheme and modest parameterization facilitate easy adaptation for real-time processing or automated quality-driven parameter selection. The color extension is straightforward, indicating applicability for natural imaging tasks. Furthermore, the SAR and ultrasound experiments position these PDE designs as highly competitive alternatives amid growing reliance on computational imaging in operational settings.

Theoretically, the formulation enriches the class of PDE-based denoising approaches, motivating future analysis of well-posedness, convergence, and optimality of such hybrid-order schemes and their relation to variational models.

Future Directions

A key limitation is the reliance on first-order explicit time integration, which may constrain stability and time-accuracy; higher-order or adaptive time-stepping could further reduce artifacts and improve convergence. Exploring spatially and temporally adaptive mixing of second/fourth-order terms and integrating learned priors could yield even higher adaptability. Deeper theoretical investigation into the optimal mixing parameterization—particularly in the regime of strongly non-stationary speckle—remains open. On the application side, the coupling framework may naturally extend to multimodal image fusion, 3D volumetric imaging, or unsupervised learning-based parameter optimization.

Conclusion

This comparative study demonstrates that weighted and coupled combinations of second- and fourth-order PDEs markedly advance the trade-off between speckle noise suppression and intricate structure retention across varied modalities. Model 1, with its weighted formulation, and Model 2, with its explicit coupling strategy, both consistently surpass recent state-of-the-art PDE models on PSNR, MSSIM, and S.I. Model 1 is especially effective for perceptual and structural fidelity. The results substantiate the value of advanced PDE blending strategies for general-purpose, interpretable, and high-fidelity image restoration in challenging noise environments, positioning these approaches as robust baselines for future investigations into both theory and deployment of despeckling algorithms.