目录
13.1 数字图像概述
13.2 图像的基本空域处理
Matlab代码示例：均值滤波与锐化滤波
13.3 图像的频域处理
Matlab代码示例：傅里叶变换与频域滤波
13.4 图像的边缘检测
Matlab代码示例：Sobel与Canny边缘检测
13.5 图像的形态学处理
Matlab代码示例：形态学处理
习题 13
总结
数字图像处理是通过计算机对数字图像进行分析和处理的一门技术，广泛应用于医学图像、工业检测、遥感、计算机视觉等领域。通过图像处理技术，可以对图像进行增强、压缩、分割、特征提取等操作，以便于人类或机器进一步理解和利用。本章将介绍数字图像处理的基本概念、常见的图像处理技术，以及在Matlab中的实现。
13.1 数字图像概述
数字图像是由离散的像素组成的，每个像素通常表示为一个数值，代表其灰度级或颜色信息。数字图像处理的核心在于对这些像素值进行各种数学运算，从而达到增强图像质量、提取有用信息等目的。
-
灰度图像：灰度图像是只有一个通道的图像，每个像素的值表示亮度，通常在0到255之间。
-
彩色图像：彩色图像通常由三个通道组成（如RGB），每个通道表示不同的颜色分量。
-
分辨率：图像的分辨率表示图像的清晰度，通常以像素数来衡量。
13.2 图像的基本空域处理
空域处理是直接对图像像素进行操作的方法，通常用于图像增强，包括平滑、锐化等操作。
-
平滑滤波：平滑滤波用于降低图像中的噪声，使得图像更加柔和。常见的方法有均值滤波和高斯滤波。
-
锐化滤波：锐化滤波用于突出图像中的边缘信息，使得图像更加清晰。常见的方法有拉普拉斯滤波和Sobel滤波。
Matlab代码示例：均值滤波与锐化滤波
% 读取图像
img = imread('example.jpg');
% 转换为灰度图像
grayImg = rgb2gray(img);
% 进行均值滤波
h = fspecial('average', [5 5]);
smoothedImg = imfilter(grayImg, h);
% 进行锐化处理
sharpenedImg = imsharpen(grayImg);
% 显示结果
figure;
subplot(1, 3, 1), imshow(grayImg), title('原始灰度图像');
subplot(1, 3, 2), imshow(smoothedImg), title('均值滤波后的图像');
subplot(1, 3, 3), imshow(sharpenedImg), title('锐化后的图像');在上述代码中，我们读取了一幅图像并将其转换为灰度图像，然后使用fspecial和imfilter函数对图像进行均值滤波，同时使用imsharpen函数对图像进行锐化处理。
13.3 图像的频域处理
频域处理是将图像从空域变换到频域，通过傅里叶变换来分析图像的频率特性。频域处理常用于图像的增强和滤波。
-
傅里叶变换：傅里叶变换将图像从空域转换为频域，得到图像的频谱信息。
-
频域滤波：在频域中对图像进行滤波，可以实现低通滤波（去除高频噪声）或高通滤波（增强边缘）。
Matlab代码示例：傅里叶变换与频域滤波
% 读取图像并转换为灰度图像
grayImg = rgb2gray(imread('example.jpg'));
% 进行傅里叶变换
F = fft2(grayImg);
Fshift = fftshift(F);
% 计算频谱
g = log(1 + abs(Fshift));
% 显示频谱图
figure;
imshow(mat2gray(g)), title('频谱图');
% 设计一个理想低通滤波器
[M, N] = size(grayImg);
D0 = 50;
[U, V] = meshgrid(1:N, 1:M);
D = sqrt((U - N/2).^2 + (V - M/2).^2);
H = double(D <= D0);
% 进行滤波
G = Fshift .* H;
Gshift = ifftshift(G);
filteredImg = real(ifft2(Gshift));
% 显示滤波结果
figure;
imshow(mat2gray(filteredImg)), title('低通滤波后的图像');在该代码中，我们使用fft2函数对图像进行了傅里叶变换，并设计了一个理想低通滤波器对图像进行频域滤波，最后使用ifft2函数将图像转换回空域。
13.4 图像的边缘检测
边缘检测是数字图像处理中非常重要的操作，用于提取图像中显著的边缘信息，通常用于对象识别和图像分割。
-
Sobel算子：通过计算图像在水平方向和垂直方向上的梯度来检测边缘。
-
Canny算子：一种多级边缘检测算法，通过计算图像梯度、应用非极大值抑制和双阈值处理，得到更精确的边缘。
Matlab代码示例：Sobel与Canny边缘检测
% 读取灰度图像
grayImg = rgb2gray(imread('example.jpg'));
% Sobel边缘检测
sobelEdges = edge(grayImg, 'sobel');
% Canny边缘检测
cannyEdges = edge(grayImg, 'canny');
% 显示结果
figure;
subplot(1, 2, 1), imshow(sobelEdges), title('Sobel边缘检测');
subplot(1, 2, 2), imshow(cannyEdges), title('Canny边缘检测');该代码展示了如何使用edge函数进行Sobel和Canny边缘检测，提取图像中的边缘信息。
13.5 图像的形态学处理
形态学处理是一种基于形状的图像处理方法，主要用于二值图像的分析，如图像中的物体分割、边缘提取等。常见的形态学操作包括腐蚀、膨胀、开运算和闭运算。
-
膨胀与腐蚀：膨胀使图像中的物体变大，填补空洞；腐蚀则使物体变小，去除噪声。
-
开运算与闭运算：开运算用于去除小的物体，保留主要结构；闭运算用于填补小的孔洞。
Matlab代码示例：形态学处理
% 读取二值图像
bwImg = imread('binary_example.png');
% 进行膨胀操作
se = strel('disk', 5);
dilatedImg = imdilate(bwImg, se);
% 进行腐蚀操作
erodedImg = imerode(bwImg, se);
% 显示结果
figure;
subplot(1, 3, 1), imshow(bwImg), title('原始二值图像');
subplot(1, 3, 2), imshow(dilatedImg), title('膨胀后的图像');
subplot(1, 3, 3), imshow(erodedImg), title('腐蚀后的图像');在该代码中，我们对二值图像进行了膨胀和腐蚀处理，展示了形态学操作在图像分析中的应用。
习题 13
在第十三章结束后，提供了一些相关的习题，帮助读者深入理解数字图像处理的基本方法。习题13包括：
图像平滑与锐化：对一幅图像进行平滑和锐化处理，观察处理效果的变化。
频域滤波：使用不同的滤波器对一幅图像进行频域滤波，比较不同滤波器的效果。
边缘检测：对不同类型的图像使用Sobel和Canny算子进行边缘检测，比较其边缘提取效果。
形态学操作：对一幅二值图像进行开运算和闭运算，分析其对图像中噪声和孔洞的影响。
通过这些习题，读者可以进一步掌握数字图像处理的各种技术，以及如何利用Matlab工具进行图像的处理与分析。
总结
第十三章介绍了数字图像处理的基本概念和常见方法，包括空域处理、频域处理、边缘检测和形态学操作等。数字图像处理在许多领域有着广泛的应用，通过对图像进行处理，可以提取有用的信息、提高图像质量，帮助我们更好地理解和利用图像数据。通过本章的学习，读者可以掌握常见的图像处理技术，并利用Matlab对图像进行各种操作和分析。
