Преобразование Фурье на C#

[MypkAnn]

Есть входные данные - либо синус, косинус, либо картинки.
Нужно сделать прямое и обратное преобразования Фурье (БПФ или ДПФ) и что бы после этих процедур было возможно вывести на экран результат ( тот же синус или картинку).
Вот прямое преобразование Фурье:

        const int N = 50;
        double[] InnerArray = new double[N];
         DSP.PF.Complexs[] Furie = new DSP.PF.Complexs[N];
         DSP.PF.Complexs[] OuterArray = new DSP.PF.Complexs[N];
        
        public class Complexs //Класс комплексных чисел
        {
            public double Re;//Реальная часть
            public double Im;//Мнимая
         }

            Complexs[] c=new Complexs[N];
            Complexs[] c0=new Complexs[N / 2];
            Complexs[] c2 = new Complexs[N / 2];
            Complexs[] Res = new Complexs[N];
            Complexs c1 = new Complexs();

            int n = (int)c.Length;
            if (n == 1) return;
            for (int i = 0, j = 0; i < n; i += 2, ++j)
            {
                if (j<n)
                c0[j] = c[i];
                c2[j] = c[i + 1];
            }

            double Arg = 2 * Math.PI / n * (invert ? -1 : 1);
            for (int i = 0; i < n / 2; ++i)
            {
                Inner[i] = (Math.Cos(2 * Math.PI * i / N) / 2) * 400 + 230;

                c[i].Re += Inner[i] * Math.Cos(Arg);
                c[i].Im -= Inner[i] * Math.Sin(Arg);
          
                Res[i].Re = c0[i].Re * c2[i].Re - c0[i].Im * c2[i].Im;
                Res[i].Im = c0[i].Re * c2[i].Im + c0[i].Im * c2[i].Re;

Как реализовать обратное преобразование Фурье что бы можно было вывести результаты, т.е. как избавиться от комплексных чисел?

[ds.Dante]

Я бы задал такой вопрос математикам, а не программистам. Нормальный программист скорее заюзает готовый код, а лучше - библиотеку.

[MypkAnn]

так в том то и дело, что кода нет и нужно его написать...

[Альтернатива]

public void FourierTransformnoise()//обратное дискретное преобразование Фурье
        {
            Complex ImOne = new Complex(0, 1);
            Complex complex;
            Complex sum;
            Complex sk;
            for (int n = 0; n < N; n++)
            {
                sum = new Complex(0, 0);
                for (int k = 0; k < K; k++)
                {
                    sk = new Complex(listNoise[k].Y, 0);
                    complex = new Complex(-2 * Math.PI * k * n / K, 0);
                    complex = Complex.Multiply(ImOne, complex);
                    complex = Complex.Exp(complex);
                    complex = Complex.Multiply(sk, complex);
                    sum = Complex.Add(sum, complex);
                }
                listSSHF.Add(n, Complex.Abs(sum));
            }
        }

P.S. если не актуально,то может пригодится кому)

[WorldMaster]

Есть чудесная библа AlgLib... если нужен именно метод то проще ее взять и не парить моцк.