《数字图像处理基础》由会员分享，可在线阅读，更多相关《数字图像处理基础（72页珍藏版）》请在金锄头文库上搜索。
1、数字图象处理北京大学计算机研究所陈晓鸥2021/6/161第二章数字图像处理基础2.1 图像运算算术运算逻辑运算2.2 空域变换几何变换非几何变换2.3 频域变换傅立叶变换导言傅立叶变换的特性快速傅立叶变换第第二二章章 数数字字图图象象处处理理基基础础2021/6/162第一节图象运算2.1.1 算术运算加法、减法乘法、除法2.1.2 逻辑运算求反异或、或与第第二二章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算2021/6/1632.1.1 图象运算：算术运算加法运算的定义C(x,y) = A(x,y) + B(x,y)主要应用举例去除“叠加性”噪音生成图象叠加效果第第二二
2、章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算2021/6/1642.1.1 图象运算：算术运算去除“叠加性”噪音对于原图象f(x,y),有一个噪音图象集 gi(x,y) i =1,2,.M其中：gi(x,y) = f(x,y) + h(x,y)iM个图象的均值定义为：g(x,y) = 1/M (g0(x,y)+g1(x,y)+ gM(x,y)当：噪音h(x,y)i为互不相关，且均值为0时，上述图象均值将降低噪音的影响。第第二二章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算2021/6/1652.1.1 图象运算：算术运算生成图象叠加效果对于两个图象f(
3、x,y)和h(x,y)的均值有：g(x,y) = 1/2f(x,y) + 1/2h(x,y)会得到二次暴光的效果。推广这个公式为：g(x,y) = f(x,y) + h(x,y)其中+= 1 我们可以得到各种图象合成的效果，也可以用于两张图片的衔接第第二二章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算2021/6/1662.1.1 图象运算：算术运算第第二二章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算2021/6/1672.1.1 图象运算：算术运算减法的定义C(x,y) = A(x,y) - B(x,y)主要应用举例去除不需要的叠加性图案检测同一场景
4、两幅图象之间的变化计算物体边界的梯度第第二二章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算2021/6/1682.1.1 图象运算：算术运算去除不需要的叠加性图案设：背景图象b(x,y)，前景背景混合图象f(x,y)g(x,y) = f(x,y) b(x,y)g(x,y) 为去除了背景的图象。电视制作的蓝屏技术就基于此第第二二章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算2021/6/1692.1.1 图象运算：算术运算检测同一场景两幅图象之间的变化设：时间1的图象为T1(x,y)，时间2的图象为T2(x,y) g(x,y) = T2 (x,y) - T
5、1(x,y)= =- -第第二二章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算2021/6/16102.1.1 图象运算：算术运算计算物体边界的梯度在一个图象内，寻找边缘时，梯度幅度（描绘变化陡峭程度的量）的近似计算 |Vf(x,y)| = max(f(x,y)f(x+1,y) ,f(x,y)f(x,y+1) 以后还会讲到第第二二章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算2021/6/16112.1.1 图象运算：算术运算第第二二章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算2021/6/16122.1.1 图象运算：算术运算乘法
6、的定义C(x,y) = A(x,y) * B(x,y) 主要应用举例图象的局部显示用二值蒙板图象与原图象做乘法第第二二章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算2021/6/16132.1.1 图象运算：算术运算第第二二章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算2021/6/16142.1.2 图象运算：逻辑运算2.1.2 逻辑运算求反异或、或与第第二二章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算2021/6/16152.1.1 图象运算：逻辑运算求反的定义g(x,y) = 255 - f(x,y)主要应用举例获得一个阴图象
7、获得一个子图像的补图像绘制区别于背景的、可恢复的图形第第二二章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算2021/6/16162.1.1 图象运算：逻辑运算获得一个阴图象第第二二章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算2021/6/16172.1.1 图象运算：逻辑运算获得一个子图像的补图像255-第第二二章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算= =2021/6/16182.1.1 图象运算：逻辑运算绘制区别于背景的、可恢复的图形第第二二章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算原图原图求反画图
8、求反画图求反恢复原图求反恢复原图2021/6/16192.1.1 图象运算：逻辑运算异或运算的定义g(x,y) = f(x,y) h(x,y)主要应用举例获得相交子图象绘制区别于背景的、可恢复的图形第第二二章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算2021/6/16202.1.1 图象运算：逻辑运算获得相交子图象第第二二章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算 = =2021/6/16212.1.1 图象运算：逻辑运算绘制区别于背景的、可恢复的图形第第二二章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算原图原图异或画图异或画图
9、异或恢复原图异或恢复原图2021/6/16222.1.1 图象运算：逻辑运算或运算的定义 g(x,y) = f(x,y) v h(x,y)主要应用举例合并子图像第第二二章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算 = =2021/6/16232.1.1 图象运算：逻辑运算与运算的定义 g(x,y) = f(x,y) h(x,y)主要应用举例求两个子图像的相交子图第第二二章章数数字字图图象象处处理理基基础础 第第一一节节图图象象运运算算 = =2021/6/1624第二节空域变换2.2.1 几何变换基本变换灰度级插值2.2.2 非几何变换模板运算灰度级变换直方图第第二二章章
10、数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16252.2.1 几何变换: 基本变换基本几何变换的定义常用的基本几何变换平移变换旋转变换镜像变换：水平镜像、垂直镜像放缩变换拉伸变换离散几何变换的计算第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16262.2.1 几何变换: 基本变换基本几何变换的定义对于原图象f(x,y)，坐标变换函数x = a(x,y); y = b(x,y)唯一确定了几何变换：g(x,y) = f(a(x,y), b(x,y);（1）g(x,y)是目标图象。表面看没有值的改变。第第二二章章数数字
11、字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16272.2.1 几何变换: 基本变换平移变换设: a(x,y) = x + x0; b(x,y) = y + y0; 用齐次矩阵表示：a(x,y)1 0 x0 xb(x,y)=0 1 y0 y 1 0 0 1 1第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16282.2.1 几何变换: 基本变换旋转变换:绕原点旋转 度设: a(x,y) = x * cos( ) - y * sin( ); b(x,y) = x * sin( ) + y * cos( ); 用齐次矩阵表示：
12、a(x,y)cos( ) 0x0 xb(x,y)=sin( ) 1 y0 y 1 0 0 1 1第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16292.2.1 几何变换: 基本变换0,0xy第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16302.2.1 几何变换: 基本变换水平镜像设: a(x,y) = -x; b(x,y) = y; 用齐次矩阵表示：a(x,y)-1 0 0 xb(x,y)=0 1 0 y 1 0 0 1 1第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2
13、021/6/16312.2.1 几何变换: 基本变换0,0xy第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16322.2.1 几何变换: 基本变换垂直镜像设: a(x,y) = x; b(x,y) = -y; 用齐次矩阵表示：a(x,y)1 0 0 xb(x,y)=0 -1 0 y 1 0 0 1 1第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16332.2.1 几何变换: 基本变换0,0xy第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16342.2.1 几
14、何变换: 基本变换放缩变换:x方向放缩c倍，y方向放缩d倍设: a(x,y) = x*c; b(x,y) = y*d; 用齐次矩阵表示：a(x,y)c 0 0 xb(x,y)=0 d 0 y 1 0 0 1 1第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16352.2.1 几何变换: 基本变换拉伸变换FDCBAFDCAB第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16362.2.1 几何变换:基本变换离散几何变换的计算问题向前映射法向后映射法第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空
15、空域域变变换换2021/6/16372.2.1 几何变换:基本变换离散几何变换的计算问题旋转、放缩、拉伸变换的漏点问题第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16382.2.1 几何变换:基本变换向前映射计算法g(x,y) = f(a(x,y), b(x,y);从原图象坐标计算出目标图象坐标镜像、平移变换使用这种计算方法第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16392.2.1 几何变换:基本变换向后映射计算法g(a(x,y), b(x,y) = f(x,y);从结果图象的坐标计算原图象的坐标旋
16、转、拉伸、放缩可以使用解决了漏点的问题，出现了马赛克第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16402.2.1 几何变换:灰度级插值灰度级插值最邻近插值法双线性插值（一阶插值）高阶插值第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16412.2.1 几何变换:灰度级插值最邻近插值法就是最临近点重复第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16422.2.1 几何变换:灰度级插值双线性插值（一阶插值）已知正方形的4个顶点，求正方形内部的点，有双线性方程：
17、f(x,y) = ax + by + cxy + d设4个顶点的坐标为：(x0,y0), (x1,y0), (x0,y1), (x1,y1)f(x, y0) = f(x0,y0)+xf(x1,y0)f(x0,y0) / (x1x0)f(x, y1) = f(x0,y1)+xf(x1,y1)f(x0,y1) / (x1x0).f(x, y) = f(x, y0) + yf(x, y1) f(x, y0) / (y1 y0)第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16432.2.1 几何变换:灰度级插值双线性插值（一阶插值）第第二二章章数数字字图图
18、象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16442.2.1 几何变换:灰度级插值高阶插值双线性插值的缺陷平滑作用使图象细节退化，尤其在放大时不连续性会产生不希望的结果高阶插值的实现用三次样条插值常用卷积来实现将大大增加计算量第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16452.2.2 非几何变换2.2.2 非几何变换非几何变换的定义模板运算灰度级变换直方图第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16462.2.2 非几何变换：非几何变换的定义非几何变换的定义对于原图象f(
19、x,y)，灰度值变换函数T(f(x,y)唯一确定了非几何变换：g(x,y) = T(f(x,y) （2）g(x,y)是目标图象没有几何位置的改变。第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16472.2.2 非几何变换：非几何变换的定义非几何变换的定义对于彩色原图象f(x,y)，颜色值变换函数Tr(f(x,y); Tg(f(x,y); Tb(f(x,y);唯一确定了非几何变换： gr(x,y) = Tr(f(x,y) gg(x,y) = Tg(f(x,y) gb(x,y) = Tb(f(x,y) 第第二二章章数数字字图图象象处处理理基基础础 第第
20、二二节节 空空域域变变换换2021/6/16482.2.2 非几何变换:模板运算模板运算模板的定义模板大小模板系数模板运算的定义模板运算举例：均值变换第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16492.2.2 非几何变换:模板运算模板的定义所谓模板就是一个系数矩阵模板大小：经常是奇数，如：3x3 5x5 7x7模板系数: 矩阵的元素w1 w2 w3w4 w5 w6w7 w8 w9第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16502.2.2 非几何变换:模板运算模板运算的定义对于某图象的子图像：z1
21、 z2 z3 z4 z5 z6 z7 z8 z9z5的模板运算公式为：R = w1z1 + w2z2 + . + w9z9第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16512.2.2 非几何变换:模板运算模板运算举例：均值变换模板系数：wi = 1/9计算公式： R = 1/9(w1z1 + w2z2 + . + w9z9)第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16522.2.2 非几何变换:灰度级变换灰度级变换灰度级变换的定义灰度级变换的实现灰度级变换举例图象求反对比度拉伸动态范围压缩灰度级切
22、片第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16532.2.2 非几何变换:灰度级变换灰度级变换（点运算）的定义(1)对于输入图象f(x,y)，灰度级变换T将产生一个输出图像g(x,y)，且g(x,y)的每一个像素值，都是由f(x,y)的对应输入像素点的值决定的。g(x,y) = T(f(x,y)第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16542.2.2 非几何变换：非几何变换的定义灰度级变换（点运算）的定义(2) 对于原图象f(x,y)，灰度值变换函数T(f(x,y) 由于灰度值总是有限个如：0
23、-255 非几何变换可定义为:R = T(r) 其中R,r在0-255之间取值第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16552.2.2 非几何变换：非几何变换的定义灰度级变换（点运算）的实现R = T(r) 定义了输入像素值与输出像素之间的映射关系，通常通过查表来实现。因此灰度级变换也被称为LUT（Look Up Table）变换。 0 1 2 3 4 5 6 7 8 9 250 251 252 253 254 255 0 3 5 7 9 11 13 15 17 19 252 253 254 254 254 255第第二二章章数数字字图图象象
24、处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16562.2.2 非几何变换:灰度级变换灰度级变换举例图象求反0255255第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16572.2.2 非几何变换:灰度级变换灰度级变换举例对比度拉伸第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16582.2.2 非几何变换:灰度级变换灰度级变换举例动态范围压缩0255255第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16592.2.2 非几何变换:灰度级变换
25、灰度级变换举例灰度级切片0255255第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16602.2.2 非几何变换:直方图直方图图象直方图的定义直方图应用举例直方图均衡化直方图匹配第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16612.2.2 非几何变换:直方图图象直方图的定义（1） 一个灰度级别在范围0，L-1的数字图象的直方图是一个离散函数 p(rk)= nk/n n 是图象的像素总数 nk是图象中第k个灰度级的像素总数 rk 是第k个灰度级，k = 0,1,2,L-1第第二二章章数数字字图图象象处处
26、理理基基础础 第第二二节节 空空域域变变换换2021/6/16622.2.2 非几何变换:直方图图象直方图的定义举例 p(rk) Nk第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换0.1 0.2 0.3 0.4 31 15 7 23 2021/6/16632.2.2 非几何变换:直方图图象直方图的定义（2） 一个灰度级别在范围0，L-1的数字图象的直方图是一个离散函数 p(rk)= nk k = 0,1,2,L-1 由于rk的增量是1，直方图可表示为：p(k)= nk即，图象中不同灰度级像素出现的次数第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空
27、域域变变换换2021/6/16642.2.2 非几何变换:直方图两种图象直方图定义的比较 p(rk)= nkp(rk)= nk/n使函数值正则化到0,1区间，成为实数函数函数值的范围与象素的总数无关给出灰度级rk在图象中出现的概率密度统计第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16652.2.2 非几何变换:直方图较暗图象的直方图 p(rk) nk第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16662.2.2 非几何变换:直方图较亮图象的直方图 p(rk) nk第第二二章章数数字字图图象象处处理理基
28、基础础 第第二二节节 空空域域变变换换2021/6/16672.2.2 非几何变换:直方图对比度较低图象的直方图 p(rk) nk第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16682.2.2 非几何变换:直方图对比度较高图象的直方图 p(rk) nk第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/16692.2.2 非几何变换:直方图直方图应用举例直方图均衡化一种自动调节图象对比度质量的算法使用的方法是灰度级变换：s = T(r) 基本思想是通过灰度级r的概率密度函数p(rk )，求出灰度级变换T(r) ，建立等值像素出现的次数与结果图象像素值之间的关系。第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/1670请提问第第二二章章数数字字图图象象处处理理基基础础 第第二二节节 空空域域变变换换2021/6/1671 结束语结束语若有不当之处，请指正，谢谢！若有不当之处，请指正，谢谢！
