# 4. Median filter (removes any remaining salt-and-pepper noise) denoised = cv2.medianBlur((denoised * 255).astype(np.uint8), 3).astype(np.float32) / 255.0
# 3. Strong Bilateral filter (smooths while keeping edges) denoised = denoise_bilateral(denoised, sigma_color=0.3, sigma_spatial=5, multichannel=(image.ndim==3))
# Load a test image original = img_as_float(data.camera()) max denoise
import numpy as np import cv2 import pywt from skimage.restoration import denoise_nl_means, denoise_bilateral from skimage.util import random_noise def max_denoise(image, sigma=0.1, h=1.15, wavelet='db8'): """ Apply maximum-strength denoising using a cascade of methods.
# Display fig, axes = plt.subplots(1, 3, figsize=(12, 4)) axes[0].imshow(original, cmap='gray') axes[0].set_title('Original') axes[1].imshow(noisy, cmap='gray') axes[1].set_title('Noisy') axes[2].imshow(denoised, cmap='gray') axes[2].set_title('Max Denoised') for ax in axes: ax.axis('off') plt.tight_layout() plt.show() # Display fig, axes = plt
# Apply hard thresholding to detail coefficients def threshold_coeffs(coeff_list, thr): return [pywt.threshold(c, thr, mode='hard') for c in coeff_list]
# Apply maximal denoising denoised = max_denoise(noisy, sigma=0.2, h=1.5) # Display fig
return np.clip(denoised, 0, 1) if name == " main ": import matplotlib.pyplot as plt from skimage import data, img_as_float