Неточность в подсчете fftw - PullRequest
       17

Неточность в подсчете fftw

1 голос
/ 05 октября 2011

Я вычисляю fft изображения, пиксели которого в два раза больше, с библиотекой FFTW, для этой цели я использую библиотеку fftw (http://www.fftw.org/) библиотека в Visual Studio 2008 и использую ITK (http://www.itk.org/).

* 1004)* Когда я вычисляю fft по двойному указателю с помощью метода 'fftw_plan_dft_c2r_2d', он выполняет вычисления, но вводит неточность, это составляет около + - 0,5 ...

Я действительно не знаю, что не так, может бытьрасстояние между пикселями? Тип используемых данных? У кого-то была эта проблема или что-то подобное раньше?

Большое спасибо в продвинутом.

Антонио

КОД:

Код здесь:

// FFT CALCULATION
    // Inizialization of the neccesary elemnets for calculating the FFT. 
    fftw_plan p1;  //variable for storing the FFT 
    int N_fft= ancho*alto;  //number of points of the image
    fftw_complex *F1 =(fftw_complex*) fftw_malloc(sizeof(fftw_complex)*alto*((ancho/2)+1));  // result pointer
    p1 = fftw_plan_dft_r2c_2d(alto,ancho, f1, F1, FFTW_ESTIMATE);  // FFT planning
    fftw_execute(p1);  // FFT calculation
    fftw_destroy_plan(p1); // plan is destroyed 

Я хочу сделать БПФ для двойного указателя с именем f1.

...