【正文】 ILE *fpr,*fpg,*fpb。 if((fpg=fopen(,wb))==NULL) exit(0)。 fseek(fpb,0,SEEK_SET)。width,4,1,fpg)。height,4,1,fpb)。height,4,1,fpr)。 fclose(fpr)。 if((fp=fopen(,rb))==NULL) exit(0)。 fclose(fp)。 for(h=0。i++) { q_count=0。kq_width。sg_width*g_height)//是否越介 data_buf[g_count]=((float)graph_data[s])*quad_data[q_count]+data_buf[g_count]。ig_height*g_width。ig_height*g_width。} //Sobel or Priwitt operator functionvoid Sobel (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。 for(h=0。i++) { q_count=0。k3。sg_width*g_height)//是否越介 x_buf[g_count]=((float)graph_data[s])*quadx[q_count]+x_buf[g_count]。 for(h=0。i++) { q_count=0。k3。sg_width*g_height)//是否越介 y_buf[g_count]=((float)graph_data[s])*quady[q_count]+y_buf[g_count]。ig_height*g_width。i++) { if(x_buf[i]maxdata) maxdata=x_buf[i]。i++) { x_buf[i]=(x_buf[i]mindata)*250/(maxdatamindata)。} //Robertsvoid 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。 for(h=0。i++) { q_count=0。k2。sg_width*g_height)//是否越介 x_buf[g_count]=((float)graph_data[s])*quadx[q_count]+x_buf[g_count]。 for(h=0。i++) { q_count=0。k2。sg_width*g_height)//是否越介 y_buf[g_count]=((float)graph_data[s])*quady[q_count]+y_buf[g_count]。ig_height*g_width。i++) { if(x_buf[i]maxdata) maxdata=x_buf[i]。i++) { x_buf[i]=(x_buf[i]mindata)*250/(maxdatamindata)。} //canny算子void Canny (unsigned char *gray_buf,unsigned long width,unsigned long height) { unsigned long i=0。iwidth*height。 for(i=0。 get_graph(gray_buf,width,height)。hheight。 maxdata=0。k3。swidth*height)//是否越介 if (gray_buf[s]=maxdata) maxdata=gray_buf[s]。 } }}void get_graph(unsigned char * gray_buf,unsigned long width,unsigned long height){ unsigned long i=0,j=0,x=0。i++) { for(j=0。 x++。 unsigned long gray_histogram[256]。i++)//初始化 { gray_histogram[i]=0。i++) //計(jì)算直方圖 gray_histogram[gray_buf[i]]++。i++) gray_p[i]=(float)gray_histogram[i]/((float)width*(float)height)+gray_p[i1]。 gray_p[0]=0。i++) { gray_b=(float)(gray_buf[i])。 gray_buf[i]=(unsigned char)gray_b。 int c=0,z=0,b=0。i++)//初始化 gray_histogram[i]=0。jwidth。i256。 r=r*j。x++) { if(gray_histogram[x]j) b++。x) { if(gray_histogram[x]j) z++。 } c=255zb。iheight。 x_buf=r*(x_bufb)+。 } } }29。 gray_buf[k]=(unsigned char)x_buf。jwidth。 k=0。 } if(z==b) if(z==256) { z=0。 } for (x=255。 for (x=0。 } r=。 k++。iheight。 //直方圖緩存 for (i=0。 long x=0。 gray_b=c*gray_b。 for (i=0。i256。 for (i=1。 } for (i=0。 for (i=0。 float gray_b=0。j++) { if(gray_buf[x]=155) gray_buf[x]= for(i=0。 } } if(maxdatagray_buf[g_count]) gray_buf[g_count]=0。 if(s=0amp。j3。iwidth。 unsigned char maxdata=0,mindata=255。i++) gray_buf[i]=255gray_buf[i]。 Sobel (quad_data6,quad_data7,gray_buf,width,height)。 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, 1 ,1 }, quad_data11[]=//Roberts Y { 1 ,1, 1 , 1 }。 } free(x_buf)。 } for(i=0。 for(i=0。 } } g_count++。 if(s=0amp。j2。h++) { for(i=0。 } } g_count++。 if(s=0amp。j2。h++) { for(i=0。 if((x_buf=calloc(g_width*g_height+4,4))==NULL) exit(0)。