void bilateralFilter(unsigned char* pSrc, unsigned char* pDest, int width, int height, int radius)

{

	float delta = 0.1f;

	float eDelta = 1.0f / (2 * delta * delta);

	int colorDistTable[256 * 256] = { 0 };

	for (int x = 0; x < 256; x++)

	{

		int  * colorDistTablePtr = colorDistTable + (x * 256);

		for (int y = 0; y < 256; y++)

		{

			int  mod = ((x - y) * (x - y))*(1.0f / 256.0f);

			colorDistTablePtr[y] = 256 * exp(-mod * eDelta);

		}

	}

	for (int Y = 0; Y < height; Y++)

	{

		int Py = Y * width;

		unsigned char* LinePD = pDest + Py;

		unsigned char* LinePS = pSrc + Py;

		for (int X = 0; X < width; X++)

		{

			int sumPix = 0;

			int sum = 0;

			int factor = 0;

			for (int i = -radius; i <= radius; i++)

			{

				unsigned char* pLine = pSrc + ((Y + i + height) % height)* width;

				int cPix = 0;

				int  * colorDistPtr = colorDistTable + LinePS[X] * 256;

				for (int j = -radius; j <= radius; j++)

				{

					cPix = pLine[ (X + j+width)%width];

					factor = colorDistPtr[cPix];

					sum += factor;

					sumPix += (factor *cPix);

				}

			}

			LinePD[X] = (sumPix / sum);

		}

	}

}