Canny算子C++【DOC精选】.doc

Canny算子C++语言实现： 包含三个文件：EdgeDetect.cpp；bmp_io.cpp；bmp_io.h； 1. EdgeDetect.cpp #include stdio.h #include stdlib.h #include math.h #include bmp_io.h #define WINDOW 11 #define ST -5.0 #define VAR 1.0 #define H_THRESHOLD 100.0 #define L_THRESHOLD 15.0 void translate(FILE *f, double **A, BMP N, unsigned char **Q){ int i,j; double max=A[0][0],min=A[0][0]; writehead(f,N.Header); writepalette(f,N.Palette,N.Header.palette_size); for(i=(WINDOW-1)/2;iN.Header.height-(WINDOW-1)/2;i++) for(j=(WINDOW-1)/2;jN.Header.width-(WINDOW-1)/2;j++){ if(A[i][j]min) min=A[i][j]; if(A[i][j]max) max=A[i][j]; } for(i=0;iN.Header.height;i++) for(j=0;jN.Header.width;j++){ if(A[i][j]max) Q[i][j]=255; else if(A[i][j]min) Q[i][j]=0; else Q[i][j]=(unsigned char)(((A[i][j]-min)*254.0)/(max-min)+1.000000001); fprintf(f,%c,Q[i][j]); } return; } int main(void){ FILE *fp = fopen(Gray_ChienI.bmp,rb); FILE *fp2 = fopen(Chien_Edge_NT.bmp,wb); BMP pic; double **I; double **Ix; double **Iy; double **J; unsigned char **Q; int *Queue; double Gx[WINDOW][WINDOW],Gy[WINDOW][WINDOW]; int i,j,k,l,W,H,cnt,d[10][2],QPtr=0,DD[8][2]; double x=ST,y=ST; for(i=0;iWINDOW;i++){ y=ST; for(j=0;jWINDOW;j++){ Gx[i][j]=0.0; Gy[i][j]=0.0; Gx[i][j]=-x*exp(-(x*x+y*y)/(2*VAR)); Gy[i][j]=-y*exp(-(x*x+y*y)/(2*VAR)); y+=1.0; } x+=1.0; } readhead(fp,pic.Header); W=pic.Header.width; H=pic.Header.height; pic.p=new unsigned char [pic.Header.cpr_size]; Queue = new int [pic.Header.cpr_size]; I= new double *[H]; Ix= new double *[H]; Iy= new double *[H]; J= new double *[H]; Q = new unsigned char *[H]; for(i=0;iH;i++){ I[i] = new double [W]; Ix[i] = new double [W]; Iy[i] = new double [W]; J[i] = new double [W]; Q[i] = new unsigned char [W]; } readpalette(fp,pic.Palette,pic.Header.palette_size); for(i=0;ipic.Header.cpr_size;i++){ fscanf(fp,%c,pic.p[i]); I[i/W][i%W]=(double)pic.p[i]; } for(i=0;iH;i++) for(j=0;jW;j++) Ix[i]