Driversb.rar_C++与图片相减_eye_平滑

#include <alloc.h> #include <stdio.h> #include <graphics.h> #include <math.h> int maxopenleft=30; int eyeclosenum=0; void save(unsigned char image_in[100][100]) {FILE *fp; char p; int m,n; if((fp=fopen("c:\\save.txt","w"))==NULL) {printf("can not open file\n"); return; } p='\n'; for(m=0;m<100;m++) for(n=0;n<100;n++) { if(fprintf(fp,"%d,",image_in[m][n])==NULL) printf("file write error\n"); if(n==99) fprintf(fp,"%c",p); } } void setbmp(FILE *fp,unsigned char image_in[100][100]) /*设置图像数据*/ { int a,m,n,num; unsigned char image[30000],k1,k2,k3,b; for(m=0;m<30000;m++) { fseek(fp,54+m,0); fread(&image[m],1,1,fp); } for(m=0;m<100;m++) for(n=0;n<100;n++) { a=m*300+n*3; k1=image[a]; a++; k2=image[a]; a++; k3=image[a]; image_in[m][n]=(unsigned char)(k3*0.299+k2*0.587 +k1*0.114); } for(m=0;m<50;m++) for(n=0;n<100;n++) { b=image_in[99-m][n]; image_in[99-m][n]=image_in[m][n]; image_in[m][n]=b; } } void sub(unsigned char image1_in[100][100],unsigned char image2_in[100][100],unsigned char image_in[100][100]) { int i,j; int a; for(i=0;i<100;i++) for(j=0;j<100;j++) {a=(int)image2_in[i][j]-(int)image1_in[i][j]; if(a<0) image_in[i][j]=0; else image_in[i][j]=(unsigned char)a; } } void smooth(unsigned char image_in[100][100]) /*高斯平滑*/ { int a,b,buf,c; unsigned char image_out[100][100]; for(a=1;a<99;a++) for(b=1;b<99;b++) { buf=(int)image_in[a-1][b-1] +2*(int)image_in[a-1][b] +(int)image_in[a-1][b+1] +2*(int)image_in[a][b-1] +4*(int)image_in[a][b] +2*(int)image_in[a][b+1] +(int)image_in[a+1][b-1] +2*(int)image_in[a+1][b] +(int)image_in[a+1][b+1]; image_out[a][b]=(unsigned char)(buf/16); } for(a=1;a<99;a++) for(b=1;b<99;b++) image_in[a][b]=image_out[a][b]; } void highlow(unsigned char image_in[100][100]) /*二值化*/ { int i,j; for(i=0;i<100;i++) for(j=0;j<100;j++) if(image_in[i][j]<=50) image_in[i][j]=0; else image_in[i][j]=1; } void arrayline(unsigned char image_in[100][100])/*行列算法*/ { int i,j,k=0,k1,k2; int p=0,p1,p2; /*对列进行处理�是阈值，超过4的像素归为0*/ for(j=0;j<100;j++) { for(i=0;i<100;i++) { if(image_in[i][j]==1) k=k+1; else if(image_in[i][j]==0) { if(k<=4) k=0; else if(k>4) { k2=i-1; for(k1=0;k1<k;k1++) { image_in[k2][j]=0; k2=k2-1; } k=0; } } } } /*对行进行处理，6是阈值，超过6的像素归为0*/ for(i=0;i<100;i++) { for(j=0;j<100;j++) { if(image_in[i][j]==1) p=p+1; else if(image_in[i][j]==0) { if(p<=6) p=0; else if(p>6) { p2=j-1; for(p1=0;p1<p;p1++) { image_in[i][p2]=0; p2=p2-1; } p=0; } } } } } void regionfind(unsigned char image_in[100][100],int *p1) { int r[100][100]; int i,j,m,n; int counter[100]; int maxlabel=0; int tempvalue=0; int maxtop=0; int maxbottom=0; int prepoint=0; int latterpoint=0; int maxopensize=0; int k=1; for(i=0;i<100;i++) /*初始化r数组*/ for(j=0;j<100;j++) r[i][j]=0; for(i=1;i<100;i++) for(j=1;j<100;j++) if(image_in[i][j]==1) { if(image_in[i-1][j]==0&&image_in[i][j-1]==0) { r[i][j]=k; k++; } else if(image_in[i-1][j]==1&&image_in[i][j-1]==0) r[i][j]=r[i-1][j]; else if(image_in[i-1][j]==0&&image_in[i][j-1]==1) r[i][j]=r[i][j-1]; else if(image_in[i-1][j]==1&&image_in[i][j-1]==1) { r[i][j]=r[i-1][j]; if(r[i-1][j]!=r[i][j-1]) {prepoint=r[i-1][j]; latterpoint=r[i][j-1]; for(m=1;m<=i;m++) {for(n=1;n<100;n++) {if(m==i&&n==j) break; else if(r[m][n]==latterpoint) r[m][n]=prepoint; } if(m==i&&n==j) break; } k--; } } } else r[i][j]=0; maxlabel=1; for(i=1;i<100;i++) counter[i]=0; for(i=0;i<100;i++) for(j=0;j<100;j++) counter[r[i][j]]++; tempvalue=counter[1]; for(i=2;i<=k;i++) if(counter[i]>tempvalue) { maxlabel=i; tempvalue=counter[i]; } for(i=1;i<100;i++) {for(j=1;j<100;j++) if(r[i][j]==maxlabel) {maxtop=i; break; } if(r[i][j]==maxlabel) break; } for(i=99;i>=1;i--) { for(j=1;j<100;j++) if(r[i][j]==maxlabel) { maxbottom=i; break; } if(r[i][j]==maxlabel) break; } maxopensize=maxbottom-maxtop; printf("%d,%d\n",maxtop,maxbottom); *p1=maxopensize; } void geteyeclosenum(int x) { float i; i=(float)(x)/(float)(maxopenleft); if(i<0.1) eyeclosenum++; } void getmaxopen(int openleft) { unsigned char image3[100][100],image4[100][100],imagesub2[100][100]; unsigned char image5[100][100],image6[100][100],imagesub3[100][100]; int x2=0,x3=0; int *p3,*p4; FILE *fp3,*fp4,*fp5,*fp6; fp3=fopen("c:\\3.bmp","rb"); fp4=fopen("c:\\4.bmp","rb"); fp5=fopen("c:\\5.bmp","rb"); fp6=fopen("c:\\6.bmp","rb"); if(fp3&&fp4&&fp5&&fp6) { setbmp(fp3,image3); setbmp(fp4,image4); sub(image3,image4,imagesub2); smooth(imagesub2); highlow(imagesub2); p3=&x2; regionfind(imagesub2,p3); x2=*p3; setbmp(fp5,image5); setbmp(fp6,image6); sub(image5,image6,imagesub3); smooth(imagesub3); highlow(imagesub3); p4=&x3; regionfind(imagesub3,p4); x3=*p4; if(openleft>maxopenleft) maxopenleft=openleft; if(x2>maxopenleft) maxopenleft=x2; if(x3>maxopenleft) maxopenleft=x3; geteyeclosenum(x2); geteyeclosenum(x3); } fclose(fp3); fclose(fp4); fclose(fp5); fclose(fp6); } void recognizeimage(int x) { int j; float i; j=x; geteyeclosenum(j); i=(float)eyeclosenum/(float)3; if(i<0.075) printf("status is 1"); else printf("status is 0"); } void main() { unsigned char image1[100][100],image2[100][100],imagesub[100][100]; /*存放图像数据*/ int openleft,m,n,i; int *p1,*p2; FILE *fp1,*fp2; if((fp1=fopen("c:\\3.bmp","rb"))==NULL) /*打开位图1*/ { printf("cannot open 1.bmp\n"); exit(1); } setbmp(fp1,image1); if((fp2=fopen("c:\\4.bmp","rb"))==NULL) /*打开位图2*/ { printf("cannot open 2.bmp\n"); exit(1); } setbmp(fp2,image2); sub(image1,image2,imagesub); save(imagesub); smooth(imagesub); highlow(imagesub);save(imagesub); arrayline(imagesub); save(imagesub); p1=&openleft; regionfind(imagesub,p1);/*区域划分，获得开度*/ openleft=*p1; i=openle