图像二值化.docx

实验五 图像二值化实验目的：掌握在计算机上进行灰度图象二值化处理，包括直方图法、微分直方图方法等；实验要求：1编写图像二值化的C语言程序以及相应的显示程序。（直方图法和微分直方图法）2对指定图象进行直方图法二值化处理(T=整幅图像灰度值的平均值)，将原始图象及二值化后的图象都显示于屏幕上，查看图像分割的效果。3对指定图象进行直方图法二值化处理(T取大于或小于整幅图像灰度值的平均值的某一个值)，将原始图象及二值化后的图象都显示于屏幕上，查看图像分割的效果，与2的结果进行比较。 4、对指定图象进行微分直方图法二值化处理，将原始图象及二值化后的图象都显示于屏幕上，查看图像分割的效果，并与2的结果进行比较。实验内容：（1）直方图法二值化处理(T=图像A2灰度值的平均值)核心C语言程序：enchence() unsigned int x,y,a,b,c,d; long data1256,data=0; for(c=0;c256;c+) data1c=0; for(b=0;b64;b+) for(a=0;a64;a+) data1*(buf+b*64+a)+; for(d=0;d256;d+) data=data+data1d*d; data=data/4096; /*以图灰度值和的平均值作为阈值进行图象二值化*/ for(b=0;b64;b+) for(a=0;a64;a+) if(*(buf+b*64+a)data) *(buf+b*64+a)=125; else *(buf+b*64+a)=255; for(y=0;y64;y+) for(x=0;x64;x+) putpixel2(x*2+10,y*2+150,*(buf+(long)y*64+x); putpixel2(x*2+11,y*2+150,*(buf+(long)y*64+x); putpixel2(x*2+10,y*2+151,*(buf+(long)y*64+x); putpixel2(x*2+11,y*2+151,*(buf+(long)y*64+x); 运行结果： （2）直方图法二值化处理(T取小于整幅图像灰度值的平均值的某个值)核心程序：enchence() unsigned int x,y,a,b,c,d; long data1256,data=0; for(c=0;c256;c+) data1c=0; for(b=0;b64;b+) for(a=0;a64;a+) data1*(buf+b*64+a)+; for(d=0;d256;d+) data=data+data1d*d; data=data/4096-10; /*以一个小于图灰度值和的平均值作为阈值进行图象二值化*/ for(b=0;b64;b+) for(a=0;a64;a+) if(*(buf+b*64+a)data) *(buf+b*64+a)=0; else *(buf+b*64+a)=255; for(y=0;y64;y+) for(x=0;x64;x+) putpixel2(x*2+10,y*2+150,*(buf+(long)y*64+x); putpixel2(x*2+11,y*2+150,*(buf+(long)y*64+x); putpixel2(x*2+10,y*2+151,*(buf+(long)y*64+x); putpixel2(x*2+11,y*2+151,*(buf+(long)y*64+x); 运行结果：（3）直方图法二值化处理(T取大于整幅图像灰度值的平均值的某个值)核心程序：enchence() unsigned int x,y,a,b,c,d; long data1256,data=0; for(c=0;c256;c+) data1c=0; for(b=0;b64;b+) for(a=0;a64;a+) data1*(buf+b*64+a)+; for(d=0;d256;d+) data=data+data1d*d; data=data/4096+20; /*以一个大于图灰度值和的平均值作为阈值进行图象二值化*/ for(b=0;b64;b+) for(a=0;a64;a+) if(*(buf+b*64+a)data) *(buf+b*64+a)=0; else *(buf+b*64+a)=255; for(y=0;y64;y+) for(x=0;x64;x+) putpixel2(x*2+10,y*2+150,*(buf+(long)y*64+x); putpixel2(x*2+11,y*2+150,*(buf+(long)y*64+x); putpixel2(x*2+10,y*2+151,*(buf+(long)y*64+x); putpixel2(x*2+11,y*2+151,*(buf+(long)y*64+x);运行结果：（4） 微分直方图法二值化处理核心程序：enchence() unsigned int x,y,a,b,c,d,m,k; long data3256,data2=0,data44096; for(c=0;c256;c+) data3c=0 ; /*微分累加和存放数组初始化*/ c=0;for(b=0;b64;b+) for(a=0;a64;a+) data4b*64+a=*(buf+b*64+a); /*开辟一个数组存放各个灰度级数*/for(b=0;b64;b+)for(a=0;a64;a+) data3*(buf+b*64+a)=data3*(buf+b*64+a)+(abs(*(buf+b*64+a)-*(buf+(b-1)*64+a-1)+abs(*(buf+b*64+a)-*(buf+(b-1)*64+a)+abs(*(buf+b*64+a)-*(buf+b*64+a-1)+abs(*(buf+b*64+a)-*(buf+b*64+a+1)+abs(*(buf+b*64+a)-*(buf+(b-1)*64+a+1)+abs(*(buf+b*64+a)-*(buf+(b+1)*64+a-1)+abs(*(buf+b*64+a)-*(buf+(b+1)*64+a)+abs(*(buf+b*64+a)-*(buf+(b+1)*64+a+1)/8; /*微分累加和*/ d=0; data2=data30; for(k=1;k256;k+) if(data2data3k) data2=data3k; d=k; /*取最大微分和的灰度值为梯度值*/ for(b=0;b64;b+) for(a=0;a64;a+) if(*(buf+b*64+a)d) *(buf+b*64+a)=0; else *(buf+b*64+a)=255; /*进行二值化分割*/ for(y=0;y64;y+) for(x=0;x64;x+) putpixel2(x*2+108,y*2+12,*(buf+(long)y*64+x); putpixel2(x*2+109,y*2+12,*(buf+(long)y*64+x); putpixel2(x*2+108,y*2+13,*(buf+(long)y*64+x); putpixel2(x*2+109,y*2+13,*(buf+(long)y*64+x); 运行结果：实验小结：对比实验要求2、3的结果可以看出，对于