Advertisement

【图像增强】matlab实现几种典型的图像增强算法

阅读量:

图像增强算法旨在优化图像的视觉表现并为其后续处理(如图像分析和识别)提供更高质量的支持。以下是一些常见的图像增强算法及其工作原理:

直方图均衡化 :通过调整图像的灰度直方图使其趋于均匀分布,并显著提升了图像的空间对比度。

对数变换 :对数变换可以增强低亮度区域,使图像的整体对比度得到提升。

Gamma变换:通过设定不同的Gamma值来提升图像对比度的方法。

锐化滤波器 :使用高通滤波器(如拉普拉斯滤波器)增强图像的边缘和细节。

双边滤波:基于空间邻域关系和像素值相似程度,在图像去噪的同时有效保护了边缘细节。

高频增强 :通过提取图像的高频分量并加以增强,可以提升图像的细节。

自适应直方图均衡化(Adaptive Histogram Equalization, AHE):该方法基于图像的不同区域特性实施直方图均衡化操作,旨在有效规避全局直方图均衡化可能引入的过度增强现象。

下面是MATLAB中实现图像直方图均衡化的一个例子:

复制代码
 % 读取图像

    
 I = imread('cameraman.tif');
    
  
    
 % 将图像转换为双精度浮点型并进行直方图均衡化
    
 I_eq = histeq(double(I));
    
  
    
 % 显示原始图像和增强后的图像
    
 subplot(1, 2, 1), imshow(I), title('Original Image');
    
 subplot(1, 2, 2), imshow(I_eq), title('Enhanced Image');

以下是MATLAB中实现双边滤波的一个例子:

复制代码
 % 读取图像

    
 I = imread('cameraman.tif');
    
  
    
 % 定义双边滤波的参数
    
 sigma_s = 50; % 空间域标准差
    
 sigma_r = 10; % 范围域标准差
    
  
    
 % 进行双边滤波
    
 I_bilateral = imfilter(double(I), fspecial('motion', 9), 'replicate');
    
 I_bilateral = imfilter(I_bilateral, fspecial('gaussian', [2*sigma_r+1, 2*sigma_r+1], sigma_r), 'same', 'replicate');
    
  
    
 % 显示原始图像和双边滤波后的图像
    
 subplot(1, 2, 1), imshow(I), title('Original Image');
    
 subplot(1, 2, 2), imshow(I_bilateral), title('Bilateral Filtered Image');

对数变换增强

复制代码
 % 读取图像

    
 I = imread('cameraman.tif');
    
  
    
 % 将图像转换为双精度浮点型
    
 I_double = double(I);
    
  
    
 % 对数变换
    
 I_log = log(1 + I_double);
    
  
    
 % 显示原始图像和增强后的图像
    
 subplot(1, 2, 1), imshow(I), title('Original Image');
    
 subplot(1, 2, 2), imshow(I_log), title('Log Transformation');

幂律(伽马)变换增强

复制代码
 % 读取图像

    
 I = imread('cameraman.tif');
    
  
    
 % 定义伽马值
    
 gamma = 0.4;
    
  
    
 % 幂律变换
    
 I_gamma = impow(double(I), gamma);
    
  
    
 % 显示原始图像和增强后的图像
    
 subplot(1, 2, 1), imshow(I), title('Original Image');
    
 subplot(1, 2, 2), imshow(I_gamma), title(['Gamma Correction (γ=' num2str(gamma) ')']);

拉普拉斯锐化滤波器

复制代码
 % 读取图像

    
 I = imread('cameraman.tif');
    
  
    
 % 创建拉普拉斯滤波器核
    
 laplacianFilter = fspecial('laplacian');
    
  
    
 % 应用拉普拉斯滤波器
    
 I_sharpened = imfilter(double(I), laplacianFilter, 'same');
    
  
    
 % 显示原始图像和增强后的图像
    
 subplot(1, 2, 1), imshow(I), title('Original Image');
    
 subplot(1, 2, 2), imshow(I_sharpened), title('Laplacian Sharpening');

高频增强

复制代码
 % 读取图像

    
 I = imread('cameraman.tif');
    
  
    
 % 将图像转换为双精度浮点型
    
 I_double = double(I);
    
  
    
 % 高频增强
    
 I_highfreq = imhighfreq(I_double);
    
  
    
 % 显示原始图像和增强后的图像
    
 subplot(1, 2, 1), imshow(I), title('Original Image');
    
 subplot(1, 2, 2), imshow(I_highfreq), title('High Frequency Enhancement');

深度学习增强算法

深度学习技术得到了广泛应用,并且卷积神经网络(CNN)在该领域取得了显著成果。基于生成对抗网络的技术已成为当前最引人注目的模型之一,在图像处理方面展现出独特优势。

示例:使用GAN进行图像超分辨率

复制代码
 % 加载预训练的SRGAN模型

    
 load('trainedSRGANModel.mat');
    
  
    
 % 读取一个低分辨率图像
    
 I_lr = imread('lena_50x50.png');
    
  
    
 % 将图像转换为单通道
    
 I_lr_singleChannel = rgb2gray(I_lr);
    
  
    
 % 进行超分辨率重建
    
 I_sr = srgan(double(I_lr_singleChannel));
    
  
    
 % 显示原始低分辨率图像和超分辨率重建图像
    
 figure;
    
 subplot(1, 2, 1), imshow(I_lr_singleChannel), title('Low Resolution Image');
    
 subplot(1, 2, 2), imshow(mat2gray(I_sr)), title('Super Resolution Image');

多尺度变换增强算法

在图像处理中,多尺度变换方法包括小波变换等技术,并进而对不同频带中的信号特征进行增强处理。

示例:使用小波变换进行图像增强

复制代码
 % 读取图像

    
 I = imread('cameraman.tif');
    
  
    
 % 进行二维小波分解
    
 [C, S] = wavedec2(double(I), 3, 'haar');
    
  
    
 % 对高频子带进行增强
    
 C = wfilters(C, 'enhance');
    
  
    
 % 进行二维小波重构
    
 I_enhanced = waverec2(C, S, 'haar');
    
  
    
 % 显示原始图像和增强后的图像
    
 subplot(1, 2, 1), imshow(I), title('Original Image');
    
 subplot(1, 2, 2), imshow(I_enhanced), title('Wavelet Enhanced Image');

图像金字塔增强算法

金字塔状的图像通过将原始图象分解为多个不同分辨率层次来实现信息表征,在具体应用中不仅能够放大并增强图象中的细节信息,在处理局部特征时也展现出良好的适应性

示例:使用图像金字塔进行细节增强

复制代码
 % 读取图像

    
 I = imread('peppers.png');
    
  
    
 % 创建高斯图像金字塔
    
 [LP, HP] = impyramid(I, 'gaussian', 4);
    
  
    
 % 对高分辨率金字塔层进行增强
    
 I_hp = imresize(HP{4}, 2);
    
  
    
 % 重建图像金字塔
    
 I_enhanced = LP{3} + imresize(I_hp, 2);
    
  
    
 % 显示原始图像和增强后的图像
    
 subplot(1, 2, 1), imshow(I), title('Original Image');
    
 subplot(1, 2, 2), imshow(I_enhanced), title('Laplacian Pyramid Enhanced Image');

自适应直方图均衡化

复制代码
 % 读取图像

    
 I = imread('cameraman.tif');
    
  
    
 % 创建自适应直方图均衡化对象
    
 claheObj = vision.ContrastAdjuster('Method', 'adhisteq');
    
  
    
 % 应用自适应直方图均衡化
    
 I_adhisteq = step(claheObj, I);
    
  
    
 % 显示原始图像和增强后的图像
    
 subplot(1, 2, 1), imshow(I), title('Original Image');
    
 subplot(1, 2, 2), imshow(I_adhisteq), title('Adaptive Histogram Equalization');

全部评论 (0)

还没有任何评论哟~