영상 파일 퓨리에 변환 하는 것 좀 도와주세요.
보미나
질문 제목 : 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;
}
번호 | 제 목 | 글쓴이 | 날짜 |
---|---|---|---|
2655328 | 함수에서 포인터 (3) | 모아 | 2024-05-18 |
2655299 | 아이콘바꾸는법..? | 소 | 2024-05-18 |
2655271 | scanf 입력받으면 쭉 실행되는.. (5) | ComeOn | 2024-05-18 |
2655184 | 이 세개 프로그램좀 해결해주실.... 고수님들???? (3) | 말달리자 | 2024-05-17 |
2655129 | 다차원 배열 원소들 주소계산하는거 질문이요~ | 터전 | 2024-05-17 |
2655103 | AVRstudio MoveAxis. MoveLine 함수질문 | 앨버트 | 2024-05-16 |
2655072 | 콘솔창의 위치와 크기를 정하고 싶습니다. | 딥레드 | 2024-05-16 |
2655044 | 자료구조랑 online judge site에 대해서 질문있습니다... (1) | 뽀송뽀송 | 2024-05-16 |
2655016 | 별찍기 문제 중 이해안되는 부분이 있어요... | 도전중 | 2024-05-16 |
2654986 | 입력받은 성적을 배열에 저장하기 (3) | 안녕나야 | 2024-05-15 |
2654959 | ===========-이런거 미리 지정하는법 (1) | 김자영 | 2024-05-15 |
2654930 | 도와주십시요 ㅜㅜ 자료구조 ㅜㅜ | 도널드 | 2024-05-15 |
2654881 | PID 모듈 검색 부분 질문합니다 . | 이름없음 | 2024-05-14 |
2654850 | 조건에 만족하는 다음 수를 찾고 싶습니다. | 초엘 | 2024-05-14 |
2654790 | 두파일간 byte단위 txt파일 복사 中 | 슬찬 | 2024-05-14 |
2654768 | 사전만들기 입니다. | 여우By | 2024-05-13 |
2654716 | c언어 고수님들 답좀적어주세요 급합니다 ㅠㅠ | 한말글 | 2024-05-13 |
2654688 | c언어 ㅠㅠ 너 무 어려 워요 ㅠ 도와주세요!! | 별빛 | 2024-05-13 |
2654659 | 구구단 소스 궁금한 점 도움 부탁드립니다! (2) | 미련곰팅이 | 2024-05-12 |
2654633 | 오류가 무었인가요..? (2) | 단순랩퍼 | 2024-05-12 |