高斯模糊处理（图像复原）源代码

#define EXTRA_NAME "@gausssmooth."
#include "loadbmp.h"

#define In(x,y) lpInput[(x)+(y)*nWidth]
#define Out(x,y) lpOutput[(x)+(y)*nWidth]
#define Mediate(x,y) lpMediate[(x)+(y)*nWidth]
#define Point(x,y) lpPoints[(x)+(y)*nWidth]
#define sigma 1

void Gauss()
{
 int x,y,xx,xk,yk;
 BYTE *lpInput=new BYTE[nWidth*nHeight];
 BYTE *lpOutput=new BYTE[nWidth*nHeight];
 GetPoints(lpInput);
 int radius,window;
 double dev_inv=0.5/(sigma*sigma);
 radius = (int)ceil(3*sigma);
 window = radius*2+1;
 memcpy(lpOutput,lpInput,nWidth*nHeight);
 double* lpMediate=new double[nWidth*nHeight];

 double* mask=new double[window];
 double sum=0;
 double temp;
 double weight;
 for(x=0;x<radius;x++)
 {
  xx=(x-radius)*(x-radius);
  mask[x]=exp(-xx*dev_inv);
  mask[window-1-x]=mask[x];
  sum+=2*mask[x];
 }
 mask[radius]=1;
 sum+=1;
 for(x=0;x<window;x++)
 {
  mask[x]/=sum;
 }
 for(y=0;y<nHeight;y++)
 {
  for(x=radius;x<nWidth-radius;x++)
  {
   temp=0;
   for(xk=-radius;xk<radius+1;xk++)
    temp+=In(x+xk,y)*mask[xk+radius];
   Mediate(x,y)=temp;
  }
  for(x=0;x<radius;x++)
  {
   temp=0;
   weight=0;
   for(xk=-x;xk<radius+1;xk++)
   {
    temp+=In(x+xk,y)*mask[xk+radius];
    weight+=mask[xk+radius];
   }
   Mediate(x,y)=temp/weight;
  }
  for(x=nWidth-radius;x<nWidth;x++)
  {
   temp=0;
   weight=0;
   for(xk=-radius;xk<nWidth-x;xk++)
   {
    temp+=In(x+xk,y)*mask[xk+radius];
    weight+=mask[xk+radius];
   }
   Mediate(x,y)=temp/weight;
  }
 }
 double* Column = new double[nHeight];
 for(x=0;x<nWidth;x++)
 {
  for(y=radius;y<nHeight-radius;y++)
  {
   temp=0;
   for(yk=-radius;yk<radius+1;yk++)
    temp+=Mediate(x,y+yk)*mask[yk+radius];
   Column[y]=temp;
  }
  for(y=0;y<radius;y++)
  {
   temp=0;
   weight=0;
   for(yk=-y;yk<radius+1;yk++)
   {
    temp+=Mediate(x,y+yk)*mask[yk+radius];
    weight+=mask[yk+radius];
   }
   Column[y]=temp/weight;
  }
        for(y=nHeight-radius;y<nHeight;y++)
  {
   temp=0;
   weight=0;
   for(yk=-radius;yk<nHeight-y;yk++)
   {
    temp+=Mediate(x,y+yk)*mask[yk+radius];
    weight+=mask[yk+radius];
   }
   Column[y]=temp/weight;
  }

  for(y=0;y<nHeight;y++)
  {
   Out(x,y)=(int)Column[y];
  }
 }
 delete []Column;
 delete []mask;
 delete []lpMediate;
 PutPoints(lpOutput);
 delete lpInput;
 delete lpOutput;
}

void main(int argc, char *argv[])
{
 if(argc==2)
  FileName=argv[1];
 else
  return; 
 OpenFile();
 Gauss();
 SaveAs();
}