体积上的厄米对称

时间:2018-08-12 03:17:55

标签: matlab image-processing signal-processing fft ifft

我正在用实际值在3D图像上测试Hermitian对称性。但是,输出似乎受到一些伪影的影响。我的代码如下。如果有人可以帮助我找到问题,我将不胜感激。谢谢!

% A: Reference 3D image    
sizA = size(A);
Ak = fftshift(fftn(A)); %FFT of reference image
figure,imshow3Dfull(log(abs(Ak)+1),[]);

% Bk has only upper half of the reference image
Bk = zeros(sizA,'single'); 
Bk(1:round(sizA(1)/2),:,:)= Ak(1:round(sizA(1)/2),:,:);
figure,imshow3Dfull(log(abs(Bk)+1),[]);

% Applying Hermitian symmetry
Bk_C = Bk;
BkF = conj(flip(flip(flip(Bk,3),2),1));
Bk_C(1+round(sizA(1)/2):end,:,:) = BkF(1+round(sizA(1)/2):end,:,:);
figure,imshow3Dfull(log(abs(Bk_C)+1),[]);

B = real(ifftn(fftshift(Bk_C)));
figure,imshow3Dfull(B,[]);

0 个答案:

没有答案