fft — C ++ Image 2D быстрое преобразование Фурье

Я должен реализовать 2D FFT-преобразование на изображении (я не могу использовать библиотеку, чтобы сделать это для меня — часть курса). я использую CImg загружать и сохранять изображения. Я сделал следующий код:

CImg<Complex> FastFourier(CImg<unsigned char> &originalImage)
{
//check size in the main.cpp
CImg<Complex> resultantImage = TransformToComplex(originalImage);
vector< vector< vector< Complex > > > vectorImage = imageToVector(resultantImage);
//cout << "Transform to complex" << endl;
int size = originalImage.width();

for(int i = 0; i < size; i++)
FastFourier1D(vectorImage[i], false);

vectorImage = rotateVector(vectorImage);

for(int i = 0; i < size; i++)
FastFourier1D(vectorImage[i], false);

vectorImage = rotateVector(vectorImage);

resultantImage = vectorToImage(vectorImage);

return resultantImage;
}

А также:

void FastFourier1D(vector< vector< Complex > > &input, bool inverse)
{
int size = input.size();
double angle;

if(size <= 1)
return;

int channels = input[0].size();
vector< vector< Complex > > even;
vector< vector< Complex > > odd;

for(int i = 0; i < size; i+=2)
{
vector< Complex > tempEven;
vector< Complex > tempOdd;
for(int channelIterator = 0; channelIterator < channels; channelIterator++)
{
tempEven.push_back(input[i][channelIterator]);
tempOdd.push_back(input[i + 1][channelIterator]);
}

even.push_back(tempEven);
odd.push_back(tempOdd);
}

FastFourier1D(even, inverse);
FastFourier1D(odd, inverse);

for(int channelIterator = 0; channelIterator < channels; channelIterator++)
{
for(int i = 0; i < size / 2; i++)
{
if(inverse == false)
angle = -2.0 * (double)PI * (double)i / (double)size;
else
angle = 2.0 * (double)PI * (double)i / (double)size;

double real = cos(angle);
double imaginary = sin(angle);

Complex W;
W.setRP(real);
W.setIP(imaginary);

W = W * odd[i][channelIterator];

input[i][channelIterator] = even[i][channelIterator] + W;
input[(size / 2) + i][channelIterator] = even[i][channelIterator] - W;
}
}
}

Однако результаты не являются хорошими. Входное изображение:

БПФ (без какого-либо преобразования):

Обратное БПФ:

Как вы видите, он имеет цвета лены, но не похож на лену. Не могли бы вы помочь мне? Есть ли какая-то ошибка?