【正文】 width*g_height)//是否越介 y_buf[g_count]=((float)graph_data[s])*quady[q_count]+y_buf[g_count]。 q_count++。 } } g_count++。 } } for(i=0。ig_height*g_width。i++) x_buf[i]=pow(pow(x_buf[i],2)+pow(y_buf[i],2),)。 for(i=0。ig_height*g_width。i++) { if(x_buf[i]maxdata) maxdata=x_buf[i]。 if(x_buf[i]mindata) mindata=x_buf[i]。 } for(i=0。ig_height*g_width。i++) { x_buf[i]=(x_buf[i]mindata)*250/(maxdatamindata)。 graph_data[i]=(unsigned char)x_buf[i]。 } free(x_buf)。 free(y_buf)。 } //Roberts void Robert(float * quadx,float * quady,unsigned char *graph_data,unsigned long g_width,unsigned long g_height) { unsigned long h=0,i=0,j=0,k=0,q_count=0,g_count=0,s=0。 float *x_buf,*y_buf,maxdata=0,mindata=255。 if((x_buf=calloc(g_width*g_height+4,4))==NULL) exit(0)。 if((y_buf=calloc(g_width*g_height+4,4))==NULL) exit(0)。 for(h=0。hg_height。h++) 畢業(yè)設(shè)計 25 { for(i=0。ig_width。i++) { q_count=0。 for(j=0。j2。j++) { for(k=0。k2。k++) { s=g_width*(h1+j)+(i1+k)。 if(s=0amp。amp。sg_width*g_height)//是否越介 x_buf[g_count]=((float)graph_data[s])*quadx[q_count]+x_buf[g_count]。 q_count++。 } } g_count++。 } } g_count=0。 for(h=0。hg_height。h++) { for(i=0。ig_width。i++) { q_count=0。 for(j=0。j2。j++) { for(k=0。k2。k++) { s=g_width*(h1+j)+(i1+k)。 if(s=0amp。amp。sg_width*g_height)//是否越介 y_buf[g_count]=((float)graph_data[s])*quady[q_count]+y_buf[g_count]。 q_count++。 } } g_count++。 } } for(i=0。ig_height*g_width。i++) x_buf[i]=pow(pow(x_buf[i],2)+pow(y_buf[i],2),)。 for(i=0。ig_height*g_width。i++) { if(x_buf[i]maxdata) maxdata=x_buf[i]。 if(x_buf[i]mindata) mindata=x_buf[i]。 畢業(yè)設(shè)計 26 } for(i=0。ig_height*g_width。i++) { x_buf[i]=(x_buf[i]mindata)*250/(maxdatamindata)。 graph_data[i]=(unsigned char)x_buf[i]。 } free(x_buf)。 free(y_buf)。 } //canny 算子 void Canny (unsigned char *gray_buf,unsigned long width,unsigned long height) { unsigned long i=0。 float vaw=0。 float quad_data5[]= { , , 1, , ,1, 1 , 1, , 1, 1, ,1, , , 1, , }, quad_data6[]=//Sobel算子 X { 1, 2, 1, 0, 0, 0, 1, 2, 1 }, quad_data7[]=//Sobel算子 Y { 1, 0, 1, 2, 0, 2, 1, 0, 1 }, quad_data[]=//laplacian 算子 { 0, 1, 0, 1, 4, 1, 0, 1, 0 }, quad_data10[]=//Roberts X { 1 ,1, 畢業(yè)設(shè)計 27 1 ,1 }, quad_data11[]=//Roberts Y { 1 ,1, 1 , 1 }。 for(i=0。iwidth*height。i++) gray_buf[i]=255gray_buf[i]。 Sobel (quad_data6,quad_data7,gray_buf,width,height)。contrast (gray_buf,width,height)。 for(i=0。iwidth*height。i++) gray_buf[i]=255gray_buf[i]。 contrast (gray_buf,width,height)。 get_graph(gray_buf,width,height)。 } void nonmax_disable(unsigned char * gray_buf,unsigned long width,unsigned long height) { unsigned long h=0,i=0,j=0,k=0,q_count=0,g_count=0,s=0。 unsigned char maxdata=0,mindata=255。 for(h=0。hheight。h++) { for(i=0。iwidth。i++) { q_count=0。 maxdata=0。 for(j=0。j3。j++) { for(k=0。k3。k++) { s=width*(h1+j)+(i1+k)。 if(s=0amp。amp。swidth*height)//是否越介 if (gray_buf[s]=maxdata) maxdata=gray_buf[s]。 q_count++。 } } if(maxdatagray_buf[g_count]) gray_buf[g_count]=0。 g_count++。 } } } void get_graph(unsigned char * gray_buf,unsigned long width,unsigned long height) { 畢業(yè)設(shè)計 28 unsigned long i=0,j=0,x=0。 x=0。 for(i=0。iheight。i++) { for(j=0。jwidth。j++) { if(gray_buf[x]=155) gray_buf[x]= else gray_buf[x]=0。 x++。 } } } void equal (unsigned char *gray_buf,unsigned long width,unsigned long height) { unsigned long i=0,j=0,k=0,data_size=0。 float gray_b=0。 float c=0。 unsigned long gray_histogram[256]。 //直方圖緩存 float gray_p[256]。 for (i=0。i256。i++)//初始化 { gray_histogram[i]=0。 gray_p[i]=0。 } for (i=0。iheight*width。i++) //計算直方圖 gray_histogram[gray_buf[i]]++。 gray_p[0]=(float)gray_histogram[0]/(float)(width*height)。 for (i=1。i256。i++) gray_p[i]=(float)gray_histogram[i]/((float)width*(float)height)+gray_p[i1]。 for (i=1。i256。i++) gray_p[i]=gray_p[i]*255。 gray_p[0]=0。 k=0。 for (i=0。iheight*width。i++) { gray_b=(float)(gray_buf[i])。 c=gray_p[gray_buf[i]]/(float)(gray_buf[i])。 gray_b=c*gray_b。 if(gray_b255)//