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

avatar
作者
筋斗云
阅读量:2

图像增强算法主要目的是改善图像的视觉效果或为后续处理(如图像分析、识别等)提供更好的图像质量。以下是一些常见的图像增强算法及其原理:

  1. 直方图均衡化:通过拉伸图像的灰度直方图,使得直方图更加均匀分布,从而增加图像的对比度。

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

  3. 幂律(伽马)变换:通过调整伽马值,可以对图像的暗部和亮部进行不同程度的增强。

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

  5. 双边滤波:结合空间邻近度和像素相似度,对图像进行平滑处理同时保持边缘清晰。

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

  7. 自适应直方图均衡化:根据图像的局部区域特征进行直方图均衡化,以避免全局直方图均衡化可能导致的过度增强问题。

下面是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)进行图像超分辨率。

示例:使用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'); 

 

广告一刻

为您即时展示最新活动产品广告消息,让您随时掌握产品活动新动态!