2013-08-21 4 views
1

2 개의 실제 2D 배열을 FFTW (r2c 및 c2r 계획 사용)와 결합하는 간단한 프로그램 (ANSI C)을 작성했습니다. 모든 것이 (규모에 대한 나의 이해를 제외하고는) 잘 작동하고 결과물은 괜찮아 보입니다. (. : 이산 푸리에 변환 I :. 역 DFT) 최종 결과의 규모가 올바른지 fftw와의 회선 규모

문제

내가 처리 후 IDFT[DFT(A)DFT(B)] 있도록 두 배열의 크기를 조절하는 방법에 약간의 혼란 스러워요 것입니다. 이미지 A는 너비가 a 픽셀이고 컨볼 루션 커널이 b 픽셀 인 사각형이라고 가정 해보십시오. 그런 다음 폭이 p=a+b 인 두 개의 채워진 정사각형 배열을 만들고 내 AB 이미지를 배치하고 컨볼 루션을 수행합니다.

최종 출력의 모든 픽셀에 1/p^2을 곱하고 있습니까? 미리 감사드립니다.

http://fftw.org/doc/What-FFTW-Really-Computes.html (웹 adress은 충분히 명시 적, 내가 생각) : 당신이 대답으로 그것을 원하기 때문에

+1

이 페이지가 도움이 될 수 있습니다. http://www.fftw.org/doc/What-FFTW-Really-Computes.html –

+0

감사합니다. 이해합니다. 내 추측이 맞았습니다. 나는 DFT의 절반 만 사용하는 r2c와 c2r 루틴에 똑같이 적용되는지 확실하지 않았습니다. 그러나 당신이 링크 한 기사는 정확한 답을줍니다. 동의 할 수 있도록 답변으로 추가해주세요. 다시 감사합니다. – makhlaghi

답변

1

.

일반적으로 FFTW의 문서는 매우 유용하지만 사용자가 알아야 할 세부 정보가 충분하지 않아 어수선하게 정리되어 있습니다. 예를 들어 프로그램을 닫기 전에 일부 메모리 누수를 방지하려면 fftw_cleanup()을 호출해야합니다.