数字图像处理及MATLAB实现3

数字图像处理武汉理工大学信息学院(BasicOperationinDigitalImageProcessing)3.1图像基本运算的概述(Introduction)3.2点运算(PointOperation)3.3代数与逻辑运算(AlgebraandLogicalOperation)3.4几何运算(GeometricOperation)第3章图像基本运算3.1图像基本运算的概述(Introduction)图像基本运算点运算(PointOperation)代数运算(AlgebraOperation)逻辑运算(LogicalOperation)几何运算(GeometricOperation)按图像处理运算的数学特征,图像基本运算可分为：图像基本运算的分类3.1图像基本运算的概述(Introduction)点运算点运算是指对一幅图像中每个像素点的灰度值进行计算的方法。代数运算、逻辑运算代数运算或逻辑运算是指将两幅或多幅图像通过对应像素之间的加、减、乘、除运算或逻辑与、或、非运算得到输出图像的方法。3.1图像基本运算的概述(Introduction)几何运算几何运算就是改变图像中物体对象（像素）之间的空间关系。从变换性质来分，几何变换可以分为图像的位置变换（平移、镜像、旋转）、形状变换（放大、缩小）以及图像的复合变换等。3.2点运算(PointOperation)1.点运算的定义设输入图像的灰度为f(x,y)，输出图像的灰度为g(x,y)，则点运算可以表示为：(,)[(,)]gxyTfxy其中T[]是对f在（x，y）点值的一种数学运算，即点运算是一种像素的逐点运算，是灰度到灰度的映射过程，故称T[]为灰度变换函数。灰度变换函数3.2点运算(PointOperation)若令f(x,y)和g(x,y)在任意点（x，y）的灰度级分别为r和s，则灰度变换函数可简化表示为：[]sTr3.1对比度增大255178482550sr3.2加亮、减暗图像2552550sr128218非线性灰度变换点运算可以改变图像数据所占据的灰度值范围，从而改善图像显示效果。3.2点运算(PointOperation)点运算又称为“对比度增强”、“对比度拉伸”、“灰度变换”等，按灰度变换函数T[]的性质，可将点运算分为：点运算灰度变换增强直方图增强(5.2.2基于直方图处理的图像增强)线性灰度变换（线性点运算）非线性灰度变换（非线性点运算）分段线性灰度变换（分段线性点运算）2.点运算的分类3.2.1线性点运算（LinearPointOperation）线性点运算的灰度变换函数形式可以采用线性方程描述，即sarbsarbr图3.3线性点运算1、线性点运算黑线：0,10ba红线：蓝线：0,1ba输出灰度压缩输出灰度扩展整体变亮0,1ba输出灰度不变绿线：0,10ba输出灰度压缩，整体变暗45º线性点运算的应用1）如果a1，输出图像的对比度增大（灰度扩展）3.4对比度增大变换前变换后sarb3.2.1线性点运算（LinearPointOperation）255178482550sr2）如果0a1,输出图像的对比度减小（灰度压缩）25525514203.5降低对比度变换前变换后3.2.1线性点运算（LinearPointOperation）3）如果a为负值，暗区域将变亮，亮区域将变暗0255255变换前变换后3.2.1线性点运算（LinearPointOperation）2、分段线性点运算将感兴趣的灰度范围线性扩展，相对抑制不感兴趣的灰度区域。设f(x,y)灰度范围为[0,Mf]，g(x,y)灰度范围为[0,Mg],ayxfyxfacbyxfacayxfabcdMyxfbdbyxfbMdMyxgffg),(0),(),(]),([),(]),([),(0f(x,y)g(x,y)abcdMfMg3.2.1线性点运算（LinearPointOperation）0f(x,y)g(x,y)abcdMfMg变换前变换后分段线性点运算的应用3.2.1线性点运算（LinearPointOperation）对数变换的一般表达式为：s=clog(1+r)其中C是一个常数。图3.9对数曲线图1、非线性点运算非线性点运算的输出灰度级与输入灰度级呈非线性关系，常见的非线性灰度变换为对数变换和幂次变换。1）、对数变换rss=log(1+r)低灰度区扩展，高灰度区压缩。图像加亮、减暗。非线性拉伸不是对图像的整个灰度范围进行扩展，而是有选择地对某一灰度值范围进行扩展，其他范围的灰度值则有可能被压缩。3.2.2非线性点运算(Non-LinearPointOperation)非线性点运算应用实例1对比度拉伸效果：图像加亮、减暗3.2.2非线性点运算(Non-LinearPointOperation)非线性点运算应用实例2：傅里叶频谱的显示s=clog(1+r)此时，C=1经对数灰度变换后的频谱图原始图像原始图像的傅里叶谱3.2.2非线性点运算(Non-LinearPointOperation)幂次变换的一般形式为：scr2）、幂次变换rs0.040.10.20.40.6711.52.551025其中C和为正常数。加亮、减暗图像加暗、减亮图像1013.2.2非线性点运算(Non-LinearPointOperation)=25.0=10.0=2.5=1=0.4=0.1=0.04L-10L/2L/2L-1输入灰度级r不同的s=cr曲线及图像变换结果输出灰度级s=1.5=0.66原始图像非线性点运算应用实例3加亮、减暗图像加暗、减亮图像3.2.2非线性点运算(Non-LinearPointOperation)加暗、减亮图像3.2.2非线性点运算(Non-LinearPointOperation)思考问题:1、点运算是否会改变图像内像素点之间的空间位置关系？2、对图像灰度的拉伸，非线性拉伸与分段线性拉伸的区别？非线性拉伸不是通过在不同灰度值区间选择不同的线性方程来实现对不同灰度值区间的扩展与压缩，而是在整个灰度值范围内采用统一的非线性变换函数，利用函数的数学性质实现对不同灰度值区间的扩展与压缩。点运算是一种像素的逐点运算，它与相邻的像素之间没有运算关系，点运算不会改变图像内像素点之间的空间位置关系。3.2.2非线性点运算(Non-LinearPointOperation)3.3代数运算与逻辑运算(AlgebraandLogicalOperation)1.概念代数运算是指两幅或多幅输入图像之间进行点对点的加、减、乘、除运算得到输出图像的过程。如果记输入图像为A(x,y)和B(x,y)，输出图像为C(x,y)，则有如下四种形式：),(),(),(yxByxAyxC),(),(),(yxByxAyxC),(),(),(yxByxAyxC),(),(),(yxByxAyxC代数运算的四种基本形式逻辑运算在进行图像理解与分析领域比较有用。运用这种方法可以为图像提供模板，与其他运算方法结合起来可以获得某种特殊的效果。3.3代数运算与逻辑运算(AlgebraandLogicalOperation)逻辑运算是指将两幅或多幅图像通过对应像素之间的与、或、非逻辑运算得到输出图像的方法。1、加法运算3.3.1加法运算(Addition)),(),(),(yxByxAyxC主要应用举例：去除“叠加性”随机噪音生成图像叠加效果去除“叠加性”噪音对于原图象f(x,y),有一个噪音图像集{gi(x,y)}i=1,2,...M其中：gi(x,y)=f(x,y)+ei(x,y)当：噪音ei(x,y)为互不相关，且均值为0时，上述图象均值将降低噪音的影响。M个图像的均值为：3.3.1加法运算(Addition)随机噪声原始图像混入噪声的图像yxeyxfyxg,,,MiiMiiiyxeMyxfyxeyxfMyxg11),(1),(),(),(1),(3.3.1加法运算(Addition)MiiyxgMEyxgE1,1,MiiyxgEM1,1MiiiyxeEyxfEM1,,1MiiyxfyxfM1,,1则是的无偏估计y,xgyxf,MiiMiiiyxeMyxfyxeyxfMyxg11),(1),(),(),(1),(利用同一景物的多幅图像取平均、消除噪声。取M个图像相加求平均得到1幅新图像，一般选8幅取平均。相加3.3.1加法运算(Addition)Addition:–averagingfornoisereductionM=1M=2M=4M=16生成图象叠加效果：可以得到各种图像合成的效果，也可以用于两张图片的衔接。3.3.1加法运算(Addition)3.3.2减法运算(Subtraction)减法运算),(),(),(yxByxAyxC主要应用举例：差影法(检测同一场景两幅图像之间的变化)混合图像的分离将同一景物在不同时间拍摄的图像或同一景物在不同波段的图像相减，这就是图像的减法运算。实际中常称为差影法。差值图像提供了图像间的差值信息，能用于指导动态监测、运动目标的检测和跟踪、图像背景的消除及目标识别等。检测同一场景两幅图像之间的变化设：时刻1的图像为T1(x,y)，时刻2的图像为T2(x,y)g(x,y)=T2(x,y)-T1(x,y)=-T1(x,y)T2(x,y)g(x,y)3.3.2减法运算(Subtraction)3.3.2减法运算(Subtraction)差影法在自动现场监测中的应用1、在银行金库内，摄像头每隔一固定时间拍摄一幅图像，并与上一幅图像做差影，如果图像差别超过了预先设置的阈值，则表明可能有异常情况发生，应自动或以某种方式报警；4、利用差值图像还能鉴别出耕地及不同的作物覆盖情况。2、用于遥感图像的动态监测，差值图像可以发现森林火灾、洪水泛滥，监测灾情变化等；3、也可用于监测河口、海岸的泥沙淤积及监视江河、湖泊、海岸等的污染；3.3.2减法运算(Subtraction)图3.6差影法进行混合图像的分离（a）混合图像（b）被减图像（c）差影图像混合图像的分离3.3.2减法运算(Subtraction)消除背景影响设：背景图像b(x,y)，前景背景混合图像f(x,y)g(x,y)=f(x,y)–b(x,y)g(x,y)为去除了背景图像即去除不需要的叠加性图案背景图像差值图像（a）从病人头顶向下拍摄的X光照片（b）碘元素注入后拍摄的X光照片与背景图像的差值3.3.3乘法运算(Multiplication)乘法运算),(),(),(yxByxAyxC主要应用举例：图像的局部显示改变图像的灰度级图像的局部显示3.3.3乘法运算(Multiplication)(a)原图(b)乘以1.2(c)乘以2图3.8乘法运算结果改变图像的灰度级3.3.4除法运算(Division)除法运算简单的除法运算可用于改变图像的灰度级，常用于遥感图像处理中。),(),(),(yxByxAyxC在四种算术运算中，减法与加法在图像增强处理中最为有用。3.3.5逻辑运算(LogicalOperation)图3.7图像的逻辑运算（a）A图（b）B图(c)A、B相与结果图(d)A、B相或结果图(e)A取反结果图“与”、“或”，“非”逻辑运算逻辑运算主要以像素对像素为基础在两幅或多幅图像间进行。3.3.5逻辑运算(LogicalOperation)“与”、“或”逻辑运算可以从一幅图像中提取子图像3.4几何运算(GeometricOperation)几何运算几何运算就是改变图像中物体对象（像素）之间的空间关系。从变换性质来分，几何变换可以分为图像的位置变换（平移、镜像、旋转）、