수다닷컴

질문 제목 : 512*512 짜리 gray 그림 파일을 이산퓨리에 변환하는 과제인데 도저히 하기가 힘드네요.

현재 퓨리에 변환을 하는 것까지되는 것 같습니다. 그런데 진폭을 구해서 image 파일로 어떻게 저장하는 지를 모르겠습니다. 그리고 512*512영상이라서 시간이 꽤 걸려서 첫번째 for문을 구간을 정해서 했는데 전혀 이상하게 나오는것 같습니다.

질문 내용 : 각각 주파수 좌표에 맞게 영상의 값을 복소 지수함수로 곱한 값을 영상의 모든점에 맞게 더하면 주파수 영역을 구할수 있다고 생각해서 코드를 짰습니다.
각각 더한 값을 영상크기로 나눈뒤 진폭을 sqrt(real제곱 + image제곱) 으로 구해서
영상으로 출력할수 있게 unsigned char 배열에 넣었습니다.
출력영상은 사각형의 꼭지점부분이 하얀색을 띄고 가운데는 검정색으로 나올거라고 생각했는데
전혀 그렇게 안나오고 뒤죽 박죽이네요
시간이 오래걸려서 첫번째 for문을 100~110 이런식으로 해도 전혀 형태가 안나옵니다.
조언좀 부탁드립니다.

// dft
#include stdio.h
#include math.h
#include stdlib.h
#define hor_size 512
#define ver_size 512
#define image_size 512*512
#define pi 3.14159265
int main(void)
{
int i, j, x, y;
double magnitude = 0.0, magnitude1 = 0.0;;
unsigned char inputimagebuf[image_size] = {0,};
unsigned char outputimagebuf[image_size] = {0,};
double real = 0.0;
double image = 0.0;

file* in = fopen(lena.raw,rb);
file* out = fopen(lena_output.raw,wb);

if(in==null)
{
printf(파일오픈에러!\n);
exit(0);
}

fread(inputimagebuf, sizeof(char), image_size,in); //파일을 불러오는 부분입니다.
fclose(in);

// dft 연산을 위한 부분입니다.
for(i=0; iver_seze; i++)
{
for(j=0; jhor_size; j++)
{
real = 0.0;//복소지수함수를 실수와 허수값으로 나누어서 계산하는데 초기화를 해주는 부분입니다.
image = 0.0;

for(x=0; xver_size; x++)
{
for(y=0; yhor_size; y++)
{
//각 좌표의 값과 복소지수함수를 곱해서 더하는 부분입니다.
real += inputimagebuf[x*hor_size+y]*cos(2*pi*(i*x+j*y)/hor_size);
image += -1*inputimagebuf[x*hor_size+y]*sin(2*pi*(i*x+j*y)/hor_size);
}
}
//이미지의 크기로 나누어주는 부분입니다.
real = real /image_size;
image = image/image_size;

magnitude = hypot(real, image); //진폭을 구하는 부분입니다.
outputimagebuf[i*hor_size+j] = magnitude; //진폭을 출력을 할 배열에 저장하는 부분입니다.
}
printf(%d / %d \n, i+1, ver_size);
}

//저장하는 부분입니다.
fwrite(outputimagebuf, sizeof(char), image_size, out);
return 0;
}

신청하기