高斯滤波函数(二维)my.docVIP

void Snake_Fun::Gauss_Image(WORD*m_Image_Data,int m_Image_Colunms,int m_Image_Rows,int m_Image_Layers,int kerner ,double sigma) { int colunm=m_Image_Colunms; //输入图像数据 int row=m_Image_Rows; int layer=m_Image_Layers; int plane_rc_data=row*colunm; int plane_tol_data=layer*plane_rc_data; int half_kerner=(kerner+1)/2-1; //模板步长； WORD * G_Image=new WORD[plane_tol_data]; //高斯滤波后的图像 for (int i=0;iplane_tol_data;i++) { G_Image[i]=0; } double **Gauss_temp=new double*[kerner]; //高斯矩阵 for (int i=0;ikerner;i++) { Gauss_temp[i]=new double[kerner]; } double *Gauss_One=new double[kerner*kerner]; //一维高斯数组 for (int i=0;ikerner*kerner;i++) { Gauss_One[i]=0; } double ** gau_x=new double *[kerner]; double ** gau_y=new double *[kerner]; double** gau_xx=new double*[kerner]; double **gau_yy=new double *[kerner]; for (int i=0;ikerner;i++) { gau_x[i]=new double [kerner]; gau_y[i]=new double [kerner]; gau_xx[i]=new double[kerner]; gau_yy[i]=new double [kerner]; } for (int i=0;ikerner;i++) { int row_offset=kerner-i-1; for (int j=0;jkerner;j++) { int colunm_offset=kerner-j-1; gau_x[row_offset][colunm_offset]=i-(kerner-1)/2; } } for (int i=0;ikerner;i++) { for (int j=0;jkerner;j++) { gau_y[j][i]=gau_x[i][j]; } } for (int i=0;ikerner;i++) { for (int j=0;jkerner;j++) { gau_xx[i][j]=gau_x[i][j]*gau_x[i][j]; gau_yy[i][j]=gau_y[i][j]*gau_y[i][j]; } } for (int i=0;ikerner;i++) { for (int j=0;jkerner;j++) { Gauss_temp[i][j]=exp(-(gau_xx[i][j]+gau_yy[i][j])/(2*sigma)); } } double max_ker=0; for (int i=0;ikerner;i++) //二维模板转化为一维数组 { for (int j=0;jkerner;j++) { Gauss_One[i*kerner+j]=Gauss_temp[i][j]; max_ker+=Gauss_One[i*kerner+j]; } } for (int i=0;ikerner*kerner;i++) { Gauss_One[i]=Gauss_One[i]/max_ker; } Conv2D(m_Image_Data,G_Image,Gauss_One,m_Image_Colunms,m_Image_Row