/************************************************************************/
/* double sigma = 1.0;
nR = 3
*/
/************************************************************************/
void Gauss_Smooth(uchar *pSrc, int nWidth, int nHeight,uchar* pDst)
{
// 循环控制变量
int x,y;
// 一维高斯数据滤波器
int pK[2] = {70,116};
int *pTemp = (int *)malloc(nWidth*nHeight*sizeof(int));
int Sum = 0;
int nW1 = nWidth-1;
// x方向进行滤波
for(y=0; y<nHeight; y++)
{
// 第一行
pTemp[y*nWidth] = ((int)pSrc[y*nWidth]) * pK[1]+((int)pSrc[y*nWidth+1] + 1 ) * pK[0];
for(x=1; x<nWidth-1; y++)
{
pTemp[y*nWidth + x] = (((int)pSrc[y*nWidth + x-1] + (int)pSrc[y*nWidth + x+1]) * pK[0]
+ (int)pSrc[y*nWidth + x]* pK[1])>>16;
}
// 最后一行
pTemp[y*nWidth+nW1] = ((int)pSrc[y*nWidth+nW1]) * pK[1] + ((int)pSrc[y*nWidth+nW1-1] + 1) * pK[0];
}
// y方向进行滤波
int nhW = (nHeight-1)*nWidth;
for(x=0; x<nWidth; x++)
{
// 第一列
pDst[x] = ( pTemp[x]* (pK[1]+pK[0]) +pTemp[x+nWidth]* pK[0])>>16;
for(y=1; y<nHeight-1; y++)
{
pDst[y*nWidth + x] = ((pTemp[(y-1)*nWidth + x] + pTemp[(y+1)*nWidth + x]) * pK[0]
+pTemp[y*nWidth + x]* pK[1])>>16;
}
// 最后一列
pDst[nhW+x] = (pTemp[nhW+x] * (pK[0]+pK[1])+pTemp[nhW+x-nWidth] * pK[1])>>16;
}
free(pTemp);
}