当
卷积的计算步骤:(动态演示)
-
对h(n)绕纵轴折叠,得h(-n);
-
对h(-m)移位得h(n-m);
-
将x(m)和h(n-m)所有对应项相乘之后相加得离散卷积结果y(n)。
说明:
因此,x(m)与h(n-m)的位置可对调(即输入为x(n)、单位脉冲响应为h(n)的线性时不变系统与输入为h(n)、单位脉冲响应为x(n)的线性时不变系统具有同样的输出)。
系统的稳定性与因果性
时,该线性时不变系统是稳定的。
这种稳定因果系统既是可实现的又是稳定工作的,这种系统是最主要的系统。
http://zlgc.seu.edu.cn/jpkc2/ipkc/signal/new/course/one/1_3_2.htm
这里讨论利用输入图像中像素的小邻域来产生输出图像的方法,在信号处理中这种方法称为滤波(filtering)。其中,最常用的是线性滤波:输出像素是输入邻域像素的加权和。
1.相关算子(Correlation Operator)
步骤:
例:
A = [17
24 1 8 15 h = [8 1 6
23 5 7 14 16 3 5 7
4 6 13 20 22 4 9 2]
10 12 19 21 3
11 18 25 2 9] Matlab 函数:imfilter(A,h)
2.卷积算子(Convolution)
步骤:
Matlab 函数:Matlab 函数:imfilter(A,h,'conv')% imfilter默认是相关算子,因此当进行卷积计算时需要传入参数'conv'
3.边缘效应
当对图像边缘的进行滤波时,核的一部分会位于图像边缘外面。
常用的策略包括:
1)使用常数填充:imfilter默认用0填充,这会造成处理后的图像边缘是黑色的。
2)复制边缘像素:I3 = imfilter(I,h,'replicate');
4.常用滤波
fspecial函数可以生成几种定义好的滤波器的相关算子的核。
例:unsharp masking 滤波
12345I = imread(
'moon.tif'
);
h = fspecial(
'unsharp'
);
I2 = imfilter(I,h);
imshow(I), title(
'Original Image'
)
figure, imshow(I2), title(
'Filtered Image'
)