第七章图像的逻辑运算.doc

第10章 图像的合成本章要点： 图像的差影法 图像的逻辑运算10.1图像的差影法本小节介绍： 图像的差影法基本理论 图像加运算 图像减运算10.1.1图像的差影法基本理论差影法实际是对图像进行代数运算的一种不同的叫法，代数运算是指对两幅输入图像进行点对点的加、减、乘、除计算而得到输出图像的运算。而对于相加和相乘的情况则可能不止有两幅图像参加运算。下面我们来看一看这四种代数运算的数学表达式：加运算 C（x，y）=A（x，y）B（x，y） （101）减运算 C（x，y）=A（x，y）B（x，y） （102）乘运算 C（x，y）=A（x，y）B（x，y） （103）除运算 C（x，y）=A（x，y）B（x，y） （104）其中A（x，y）和B（x，y）为输入图像，而C（x，y）为输出图像。以上四种运算在实际中有很重要的作用。图像相加的一个重要应用就是对同一场景的多幅图像求平均值，它可以有效地降低随机噪声的影响，这是因为对于一幅有噪声的图像S（x，y），可以看成是由原始无噪声的图像F（x，y）和噪声G（x，y）叠加而成的，即 S（x，y）=G（x，y）+F（x，y） （105）如果叠加在图像上的噪声G（x，y）是非相关、具有零均值的随机噪声时，那么，把针对同一目标物在相同条件下，做M次重复摄取的图像相加，取平均值作为输出图像，即： （106）这样输出的平均图像（x ,y）中随机噪声的含量，比单幅图像s（x，y）的噪声含量大大削弱了，这样，便达到了使图像变平滑的处理目的。图像相减可用于去除一幅图像中所不需要的图案，也可以用于检测同一场景的两幅图之间的变化，例如检测物体运动。对于同一场景的两幅图S(x,y)、F(x,y)，因为是同一场景，所以目标图像的背景是大致相同的，当两幅图进行相减时，相同位置的背景点因为灰度值相同，结果变成了黑点，而目标图像因为位置不同，所以相减后值不为0，这样一来，处理后的图像就只在两个目标图像的位置有像素点，背景变为全黑，此时，我们计算两个目标之间的距离就非常简单了，这也就达到了检测物体运动的目的。乘法和除法在数字图像处理中一般应用得不多，但它们的用途也很重要。数字化对一幅图像各点的敏感程度可能有变化，乘和除运算有可能纠正这种影响。除运算可产生对颜色和多光谱图像分析十分重要的比率图像，而用一幅图像乘某一图像可以遮住该图像中的某些部分，仅留下感兴趣的物体。10.1.2图像加运算1. 理论基础现在来具体分析一下加运算为什么能消减噪声。假定同一目标物在相同条件下，做M次重复摄取的图像S（x，y），其中含有的噪声G（x，y）是非相关、具有零均值的随机噪声，那么图像S（x，y）可以看成是由原始无噪声的图像F（x，y）和噪声G（x，y）叠加而成的，即 S（x，y）=G（x，y）+F（x，y）对于图像中任意点，我们定义其功率信噪比为P（x，y）= ，表示期望算符。因为我们假定每幅图的噪声都是互不相干的、噪声均值等于0的随机噪声，这就意味着全部满足： （107） （ij） （108）和 （ij） （109）如果对M幅图像作加法后平均化，可得： （1010）则功率信噪比为： （1011） 根据公式108、109可以将公式1011变形为： （1012）现在我们注意到式（1012）的分母的第二项等于0（因为107）而且，因为M 个噪声样本来自相同的样本集，分母的第一项中求和项都相同。所以有： （1013）因此，对M幅图像进行平均，使图像中每一点的功率信噪比提高了M倍。幅度信噪比是功率信噪比的开平方，即： （1014）随着图像的数目增多，其消除噪声的效果也就越来越好。下面介绍一个简单的加运算的思路，这里只计算两幅图的加法。首先将两幅图的像素读到内存缓冲区中，将两幅图的对应位置相同的像素点的灰度值相加，如果结果大于255，则置为255，否则将结果存在内存中，当所有像素点都处理完后，将内存中的结果复制回原图像。 2实现步骤(1)获取第一幅图的首地址，及图像的高和宽。(2)开辟一块内存缓冲区，并初始化为255。(3)选取进行加运算的第二幅图并获取第二幅图的首地址，及图像的高和宽。(4)将两幅图的对应像素点进行加运算，结果大于255，则置为255，不大于将结果保存下来。(5)将内存中的结果复制回第一幅图的数据区。3. 编程代码/* *函数名称：Add(LPBYTE p_data, LPBYTE p_dataBK,int wide,int height) *函数类型：void *参数说明：p_data 原图像首地址* p_dataBK 背景图像首地址* wide，height 原图像的高和宽 *功能：对图像进行加运算。 */ void TuXiangHeChengDib:Add(LPBYTE p_data, LPBYTE p_dataBK,int wide,int height) int i,j; /循环变量LPBYTE lpSrc; /原图像像素指针 LPBYTE lpDst; /目标像素指针 LPBYTE lpSrcBK; /背景像素指针 LPBYTE temp=new BYTEwide*height; /新图像缓冲区/初始化新分配的内存为255memset(temp,255,wide*height);/进行每个像素的加运算 for(j=1;jheight-1;j+)for(i=1;i255) *lpDst=255; else *lpDst=(*lpSrc+*lpSrcBK); memcpy(p_data, temp,wide*height); /复制处理后的图像delete temp; 4. 处理效果图 （a）原图 （b）背景 （c）相加效果图图10-1相加效果图由以上三幅图可以看出，加运算将两幅不同的图重叠起来，并使原来图中较暗的部分变亮了。10.1.3图像减运算1. 理论基础减运算的目的是为了从图像中去除不需要的图像，只保留所关心的图像，比如医生想观察药物在人身体中的流动情况，但是由于脊椎等的影响，图像很模糊，为了获得清晰的图像，我们可以先获得一张病人未吃药的图像，用带有药物的图像减去没有药物的图像，因为是一个人的影像，所以脊椎等干扰图像是相同的，减处理后，就只剩下药物的影像了。利用这种原理，也可以用来识别物体在特定区域的运动情况。本节介绍的是减运算的基本原理，运算的两幅图要求完全对准，这在实际中是很难达到的。在这里我们就不讨论了。 2. 实现步骤(1)获取需要提取图像的数据首地址，及图像的高和宽。(2)开辟一块内存缓冲区，初始化为255。(3)选择背景图像，并获取其首地址，及图像的高和宽。(4)用第一幅图像的每一个像素点的灰度值减去背景图像的相应像素的灰度值，结果小于0，则该点置为0，不小于0，则保存在缓冲区中。(5)将内存中的结果复制回第一幅图的数据区。3. 编程代码 /* *函数名称：Sub(LPBYTE p_data, LPBYTE p_dataBK,int wide,int height) *函数类型：void *参数说明：p_data 原图像首地址* p_dataBK 背景图像首地址* wide，height 原图像的高和宽 *功能：对图像进行减运算。 */ void TuXiangHeChengDib:Sub(LPBYTE p_data, LPBYTE p_dataBK, int wide, int height) int i,j; /循环变量 LPBYTE lpSrc; /原图像像素指针LPBYTE lpDst; /目标像素指针LPBYTE lpSrcBK; /背景像素指针LPBYTE temp=new BYTEwide*height; /新图像缓冲区/初始化新分配的内存为255memset(temp,255,wide*height);/进行每个像素的减运算for(j=1;jheight-1;j+)for(i=1;iwide-1;i+)/获得原像素指针 lpSrc=p_data+wide*j+i; /获得目标像素指针 lpDst=temp+wide*j+i; /获得背景像素指针 lpSrcBK=p_dataBK+wide*j+i; /进行减运算 if(*lpSrc-*lpSrcBK)0) *lpDst=0; else *lpDst=(*lpSrc-*lpSrcBK);memcpy(p_data, temp,wide*height); /复制处理后的图像delete temp; 4. 处理效果图 (a)原图 (b)背景 (c)相减效果图图10-2 相减效果图很明显减运算将图像进行了分离。10.2图像的逻辑运算逻辑运算是数学中关于集合的一类运算，其基本运算是交、并、补三种。与运算的法则是当两个值都是1时两个数与的结果等于1，其它全为0。或运算的法则是当两个值都是0时两个数或的结果等于0，其它全为1。非运算的法则是当数值为0结果为1，当数值为1结果为0。与或运算是将两个数先进行与运算，再进行或运算，与非方法类似。异或运算的法则是当运算的两个数取值不同，结果为1，相同则为0。本节介绍六种运算是用来处理两幅二值图: 与运算 或运算 非运算 或非运算 与非运算 异或运算10.2.1与运算1.实现步骤(1)把二值图像转化为0、1值；(2)循环两幅图像的像素点；(3)像素点相与；(4)把0、1值转化为二值图像数据，为0的像素置0，为1的像素置255。2. 编程代码/* *函数名称：Yuyunsuan(LPBYTE p_data, LPBYTE p_dataBK,int wide,int height) *函数类型：void *参数说明：p_data 原图像首地址* p_dataBK 背景图像首地址* wide，height 原图像的高和宽 *功能：对图像进行与运算。 */void TuXiangHeChengDib:Yuyunsuan(LPBYTE p_data, LPBYTE p_dataBK, int wide, int height)int i,j; /循环变量LPBYTE temp1=new BYTEwide*height; /新图像缓冲区 LPBYTE temp2=new BYTEwide*height; memset( temp1,0,wide*height); memset( temp2,0,wide*height);/对原图二值化for(j=0;jheight;j+)for(i=0;i128) temp1wide*j+i=1; /对背景二值化for(j=0;jheight;j+)for(i=0;i128) temp2wide*j+i=1; /进行每个像素的与运算 for(j=1;jheight-1;j+)for(i=1;iwide-1;i+) if( temp2wide*j+i=1& temp1wide*j+i=1) temp1wide*j+i=255;else temp1wide*j+i=0; memcpy(p_data, temp1,wide*height); /复制处理后的图像delete temp1;delete temp2;4.处理效果图 (a)原图 (b)背景 (c) 与运算效果图图10-3与运算效果图可以看出与运算在对图像进行二值化的同时保留了原图和背景中相同的部分，去除了不同的部分。10.2.2或运算1.实现步骤(1)把二值图像转化为0、1值；(2)循环两幅图像的像素点；(3)两像素点同时为0置0,否则置1；(4)把0、1值转化为二值图像数据，为0的像素置0，为1的像素置255。2.编程代码/* *函数名称：Huoyunsuan(LPBYTE p_data, LPBYTE p_dataBK,int wide,int height) *函数类型：void *参数说明：p_data 原图像首地址* p_dataBK 背景图像首地址* wide，height 原图像的高和宽 *功能：对图像进行或运算。 */ void TuXiangHeChengDib:Huoyunsuan(LPBYTE p_data, LPBYTE p_dataBK, int wide, int height) int i,j; /循环变量LPBYTE temp1=new BYTEwide*height; /新图像缓冲区 LPBYTE temp2=new BYTEwide*height; memset( temp1,0,wide*height); memset( temp2,0,wide*height);/对原图二值化for(j=0;jheight;j+)for(i=0;i128) temp1wide*j+i=1; /对背景二值化for(j=0;jheight;j+)for(i=0;i128) temp2wide*j+i=1; /进行每个像素的或运算 for(j=1;jheight-1;j+) for(i=1;iwide-1;i+) if( temp2wide*j+i=0& temp1wide*j+i=0) temp1wide*j+i=0;else temp1wide*j+i=255; memcpy(p_data, temp1,wide*height); /复制处理后的图像delete temp1;delete temp2;4处理效果图 (a)原图 (b)背景 (c) 或运算效果图图10-4 或运算效果图通过观察可以看出或运算的效果和加运算有些近似，只是或运算对图像进行了二值化。10.2.3非运算1.实现步骤(1)把二值图像转化为0、1值；(2)循环两幅图像的像素点；(3)像素点值为1置0,为0置1；(4)把0、1值转化为二值图像数据，为0的像素置0，为1的像素置255。2.编程代码/* *函数名称：Feiyunsuan() *函数类型：void *功能：对图像进行非运算。 */ void TuXiangHeChengDib:Feiyunsuan()int i,j; /循环变量BYTE* p_data;/图像指针int wide,height; /原图的高和宽p_data=GetData();wide=GetWidth();height=GetHeight();LPBYTE p_temp1=new BYTEwide*height; /新图像缓冲区 memset(p_temp1,0,wide*height); /初始化为0/对原图二值化for(j=0;jheight;j+)for(i=0;i128)p_temp1wide*j+i=1; /进行每个像素的非运算 for(j=1;jheight-1;j+) for(i=1;iwide-1;i+) if(p_temp1wide*j+i=0 ) p_temp1wide*j+i=255;else p_temp1wide*j+i=0; memcpy(p_data,p_temp1,wide*height); /复制处理后的图像delete p_temp1;4处理效果图 (a)原图 (b) 非运算效果图 图10-5非运算效果图非运算实际上是对图像二值化后进行反色变换的结果。10.2.4或非运算1.实现步骤(1)把二值图像转化为0、1值；(2)循环两幅图像的像素点；(3)两像素点同时为0置1，否则置0；(4)把0、1值转化为二值图像数据，为0的像素置0，为1的像素置255。2.编程代码/* *函数名称：Huofei(LPBYTE p_data, LPBYTE p_dataBK,int wide,int height) *函数类型：void *参数说明：p_data 原图像首地址* p_dataBK 背景图像首地址* wide，height 原图像的高和宽 *功能：对图像进行或非运算。 */ void TuXiangHeChengDib:Huofei(LPBYTE p_data, LPBYTE p_dataBK, int wide, int height)int i,j; /循环变量 LPBYTE temp1=new BYTEwide*height; /新图像缓冲区 LPBYTE temp2=new BYTEwide*height; memset( temp1,0,wide*height); memset( temp2,0,wide*height);/对原图二值化for(j=0;jheight;j+)for(i=0;i128) temp1wide*j+i=1; /对背景二值化for(j=0;jheight;j+)for(i=0;i128) temp2wide*j+i=1; /进行每个像素的与运算 for(j=1;jheight-1;j+)for(i=1;iwide-1;i+) if( temp2wide*j+i=0& temp1wide*j+i=0) temp1wide*j+i=255;else temp1wide*j+i=0; memcpy(p_data, temp1,wide*height); /复制处理后的图像delete temp1;delete temp2;4.处理效果图 (a)原图 (b)背景 (c) 或非运算效果图图10-6或非运算效果图10.2.5与非运算1.实现步骤(1)把二值图像转化为0、1值；(2)循环两幅图像的像素点；(3)两像素点同时为1置0,否则置1；(4)把0、1值转化为二值图像数据，为0的像素置0，为1的像素置255。2.编程代码/* *函数名称：Yufei(LPBYTE p_data, LPBYTE p_dataBK,int wide,int height) *函数类型：void *参数说明：p_data 原图像首地址* p_dataBK 背景图像首地址* wide，height 原图像的高和宽 *功能：对图像进行与非运算。 */void TuXiangHeChengDib:Yufei(LPBYTE p_data, LPBYTE p_dataBK, int wide, int height)int i,j; /循环变量LPBYTE temp1=new BYTEwide*height; /新图像缓冲区 LPBYTE temp2=new BYTEwide*height; memset( temp1,0,wide*height); memset( temp2,0,wide*height);/对原图二值化for(j=0;jheight;j+)for(i=0;i128) temp1wide*j+i=1; /对背景二值化for(j=0;jheight;j+)for(i=0;i128) temp2wide*j+i=1; /进行每个像素的与非运算 for(j=1;jheight-1;j+) for(i=1;iwide-1;i+) if( temp2wide*j+i=1& temp1wide*j+i=1) temp1wide*j+i=0;else temp1wide*j+i=255; memcpy(p_data, temp1,wide*height); /复制处理后的图像delete temp1;delete temp2;3效果图 (a)原图 (b)背景 (c) 与非运算效果图图10-7与非运算效果图10.2.6异或运算1.实现步骤(1)把二值图像转化为0、1值；(2)循环两幅图像的像素点；(3)两像素点相同置0,不同置1；(4)把0、1值转化为二值图像数据，为0的像素置0，为1的像素置255。2.编程代码/* *函数名称：Yihuo(LPBYTE p_data, LPBYTE p_dataBK,int wide,int height) *函数类型：v