void FT::LowpassFilter(int ** InputImage, int ** OutputImage, double ** FreqReal, double ** FreqImag, int h, int w)
{
	double filterPartA,filter;
	double n=5,d=20;
	for (int r = 0; r < h; r++)
	{
		for (int c = 0; c < w; c++)
		{
			filterPartA = pow(pow(r - h/2, 2) + pow(c - w/2, 2), 0.5);
			filter = 1 / (1 + pow((filterPartA / d), 2 * n));
			FreqReal[r][c] = FreqReal[r][c]*filter;
			FreqImag[r][c] = FreqImag[r][c]*filter;
			OutputImage[r][c] = sqrt(pow(FreqReal[r][c], (double) 2.0) + pow(FreqImag[r][c], (double) 2.0))*h;
		}
	}
}