주파수 도메인에서 컨볼 루션을 수행 할 클래스에 대해 C++ 프로그램을 작성 중이며 최종 결과에 오류가있는 것으로 나타났습니다. 그래서 나는 MATLAB에서 그것을 시도하고 똑같은 결과를 얻었다. 예 :주파수 영역의 이미지 컨볼 루션
http://engronline.ee.memphis.edu/eece7214/images/Downlodable.htm
에서 카메라맨을 사용하여 나는
a = imread('cameraman.pgm');
h = ones(25,25)/25/25;
a(512,512) = 0;
h(512,512) = 0;
c = ifft2(fft2(a).*fft2(h))/256;
c = c(1:256, 1:256);
c = real(c);
imwrite(c,'test2.png')
주요 질문과 무관하게이 MATLAB 코드에서 256으로 나누어야하는 이유를 알고 싶습니다. 내 C++ 코드에서는 결과의 크기를 조정할 필요가 없었으며이 MATLAB 코드와 동일한 대답을 얻었습니다.
편집 : 1 차원 벡터 (conv 및 ifft (fft * fft) 수행)를 사용하여 조금 놀았으며 '오류'는 출력에서 '전체'컨볼 루션을 맨 위에 표시합니다 '같은'회선 대신에 왼쪽 구석. 그러나 그 경우에도 결정론적인 방식으로 코드를 작성하는 방법을 잘 모르겠습니다. "전체"의 왼쪽 상단 256x256 부분 대신 '같은 부분'만 가져 오기 "
편집 : 더 많은 인터넷 검색 http://jeremy.fix.free.fr/IMG/pdf/fftconvolution.pdf을 통해 가능한 해결책을 찾았습니다. 그것은 이전에 보지 못했던 많은 수학 기호를 가지고 있습니다. 그러나 여러분이 nxn과 mxm을 convolving하는 경우 m : (m + n-1)을 추출하여 '같은'컨볼 루션을 얻을 수 있습니다. fft 근사. 나는 여전히 내가하는 것보다 전문가 인 사람으로부터 소식을 듣고 싶습니다. 따라서이 업데이트를 기반으로 의견을 말하지 않기로 선택하지 마십시오!