Skip to content
Snippets Groups Projects
Commit 858de8d6 authored by Robin Stampa's avatar Robin Stampa
Browse files

deleted unneeded code

parent d5e4a723
No related branches found
No related tags found
No related merge requests found
...@@ -70,7 +70,7 @@ int main() ...@@ -70,7 +70,7 @@ int main()
int phase_ERR = 0; int phase_ERR = 0;
long phase_width; long phase_width;
long phase_height; long phase_height;
unsigned char* phase_data = readpng("Images/input.png", &phase_ERR, &phase_width, &phase_height); unsigned char* phase_data = readpng("Images/input_phase.png", &phase_ERR, &phase_width, &phase_height);
if (!phase_data) if (!phase_data)
{ {
bAmplitudeWithPhase = false; bAmplitudeWithPhase = false;
...@@ -82,7 +82,6 @@ int main() ...@@ -82,7 +82,6 @@ int main()
bAmplitudeWithPhase = false; bAmplitudeWithPhase = false;
} }
fftw_complex* signal_s = new fftw_complex[image_length]; fftw_complex* signal_s = new fftw_complex[image_length];
fftw_complex* signal_t = new fftw_complex[image_length]; fftw_complex* signal_t = new fftw_complex[image_length];
fftw_plan forward_plan = fftw_plan_dft_2d(image_height, image_width, signal_s, signal_t, FFTW_FORWARD, FFTW_MEASURE); fftw_plan forward_plan = fftw_plan_dft_2d(image_height, image_width, signal_s, signal_t, FFTW_FORWARD, FFTW_MEASURE);
...@@ -102,23 +101,11 @@ int main() ...@@ -102,23 +101,11 @@ int main()
} }
fftw_execute(forward_plan); fftw_execute(forward_plan);
//find maximum
double max = 0;
for (long i = 0; i < image_length; i++)
{
if (max < own_abs_2(signal_t[i]) / image_length && i != 0)
max = own_abs_2(signal_t[i]) / image_length;
}
unsigned char* result_data = new unsigned char[image_length]; unsigned char* result_data = new unsigned char[image_length];
for (long i = 0; i < image_length; i++) for (long i = 0; i < image_length; i++)
{ {
//result_data[i] = own_abs_2(signal_t[image_switch_index(i, image_width, image_height)])/image_length*255/max;
//result_data[i] = d2uc_converter(own_abs_2(signal_t[image_switch_index(i, image_width, image_height)]) / image_length);
result_data[i] = own_abs_2(signal_t[image_switch_index(i, image_width, image_height)]) / image_length; result_data[i] = own_abs_2(signal_t[image_switch_index(i, image_width, image_height)]) / image_length;
if (result_data[i] > 0 && false)
printf("%d at %d, %d\n", result_data[i], i%image_width, i / image_width);
} }
writepng("Images/result.png", result_data, image_width, image_height); writepng("Images/result.png", result_data, image_width, image_height);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment