图像的高斯平滑C++实现

#include<iostream>#include<highgui\highgui.hpp>#include<core/core.hpp>  using namespace cv;  using namespace std;  void GetGuassFilter(float* pfGaussFilter, const int iFilterHeight, const int iFilterWidth, const float fSigma){    cout<<"pfGaussFilter"<<endl;    for(int i=0;i<iFilterHeight;i++)    {        for(int j=0;j<iFilterWidth;j++)        {            pfGaussFilter[i*iFilterHeight+j]=1.0/exp(((i-iFilterHeight/2)^2+(j-iFilterWidth)^2));            //pfGaussFilter[i*iFilterWidth+j]=1.0/25.0;            //cout<<pfGaussFilter[i*iFilterWidth+j]<<" ";        }        cout<<endl;    }    float sum=0;    for(int j=0;j<iFilterWidth*iFilterHeight;j++)    {        sum+=pfGaussFilter[j];    }    for(int j=0;j<iFilterHeight*iFilterWidth;j++)    {        pfGaussFilter[j]=pfGaussFilter[j]/sum;    }}void GaussFilt2D(uchar* pSrc, uchar* pDst, const int iHeight, const int iWidth,                 const int iFilterHeight, const int iFilterWidth, const float fSigma){    float* pfGaussFilter=new float[iFilterHeight*iFilterWidth];    GetGuassFilter(pfGaussFilter,iFilterHeight,iFilterWidth,fSigma);    float sum=0;    int index=0;    float current_src=0;    for(int i=-iFilterHeight;i<iHeight-1;i++)        for(int j=-iFilterWidth;j<iWidth-1;j++)        {            sum=0;            index=0;            for(int m=i-iFilterHeight/2;m<i+iFilterHeight/2+1;m++)            {                for(int n=j-iFilterWidth/2;n<j+iFilterWidth/2+1;n++)                {                    if(m<0||n<0||m>iHeight-1||n>iWidth-1)                        continue;                    current_src=(float)pSrc[m*iWidth+n];                    sum+=current_src*pfGaussFilter[index++];                }            }            if(i+iFilterHeight/2>0&&i+iFilterHeight/2<iHeight-1&&j+iFilterWidth/2>0&&j+iFilterWidth/2<iWidth-1)                pDst[(i+iFilterHeight/2)*iWidth+j+iFilterWidth/2]=(uchar)sum;            //pDst[i*iWidth+j]=(uchar)sum;            if(pDst[i*iWidth+j]<0)                pDst[i*iWidth+j]=0;            if(pDst[i*iWidth+j]>255)                pDst[i*iWidth+j]=255;        }}void GaussFilt2DFull(uchar* pSrc, uchar* pDst, const int iHeight, const int iWidth,                 const int iFilterHeight, const int iFilterWidth, const float fSigma){    float* pfGaussFilter=new float[iFilterHeight*iFilterWidth];    GetGuassFilter(pfGaussFilter,iFilterHeight,iFilterWidth,fSigma);    float sum=0;    int index=0;    float current_src=0;    for(int i=0;i<iHeight-1;i++)        for(int j=0;j<iWidth-1;j++)        {            sum=0;            index=0;            for(int m=i-iFilterHeight/2;m<i+iFilterHeight/2;m++)                for(int n=j-iFilterWidth/2;n<j+iFilterWidth/2;n++)                {                    if(m<0||n<0||m>iHeight-1||n>iWidth-1)                        continue;                    current_src=(float)pSrc[m*iWidth+n];                    sum+=current_src*pfGaussFilter[index++];                }            //if(i+iFilterHeight/2>0&&i+iFilterHeight/2<iHeight-1&&j+iFilterWidth/2>0&&j+iFilterWidth/2<iWidth-1)            pDst[i*iWidth+j]=(uchar)sum;            //pDst[i*iWidth+j]=(uchar)sum;            if(pDst[i*iWidth+j]<0)                pDst[i*iWidth+j]=0;            if(pDst[i*iWidth+j]>255)                pDst[i*iWidth+j]=255;        }}int main()  {  //读入图片，注意图片路径      //Mat image=imread("D:\\Picture\\lena.jpg");      //Mat image=imread("E:\\Projects\\MyOpenCV\\MyOpenCV\\1.png");     Mat image=imread("C:\\迅雷下载\\图像处理\\Projects\\MyOpenCV\\MyOpenCV\\DIP3ECH06\\Fig0638(a)(lenna_RGB).tif",0);     Mat image_Dst=imread("C:\\迅雷下载\\图像处理\\Projects\\MyOpenCV\\MyOpenCV\\DIP3ECH06\\Fig0648(b)(lenna-noise-G-gauss-mean0-var800).tif",0);     //Mat image=imread("E:\\Projects\\MyOpenCV\\MyOpenCV\\DIP3ECH06\\Fig0635(middle_row_left_chalk ).tif",0);    //Mat image_Dst=imread("E:\\Projects\\MyOpenCV\\MyOpenCV\\DIP3ECH06\\Fig0635(middle_row_left_chalk ).tif",0);    const int iHeight=image.rows;    const int iWidth=image.cols;    uchar* pSrc=image.data;//new uchar[iHeight*iWidth];    uchar* pDst=image_Dst.data;//new uchar[iHeight*iWidth];    //GaussFilt2D(pSrc,pDst,iHeight,iWidth,5,5,1);    GaussFilt2DFull(pSrc,pDst,iHeight,iWidth,25,25,1);    namedWindow("Origin",1);    imshow("Origin",image);    //创建一个名字为“Lena”的图像显示窗口，（不提前声明也可以）      namedWindow("Lena",1);      //显示图像      imshow("Lena",image_Dst);      //等待按键      waitKey();      return 0;  }