Upsampled image FFT contains NaNs.

Question

raheem mian 2020-1-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/502212-upsampled-image-fft-contains-nans

评论： Matt J 2020-1-27

I upsampled images and when I take the FFT I get NaN values. How can I prevent NaN values from showing when I take the fourier transform ?

Heres my code to upsample. Upsampling frames in a volume:

ps. ref is 3D

vol = upsample(ref);
pad = GetPad(fftr);
fftr = fftn(vol, pad);
fftt = fftn(tar, pad);
% calculate the power spectrum to remove the magnitude
cp_vol = CrossPowerSpectrum(fftr, fftt);
% go back to the time domain and find the peak
shift_vol = ifftn(cp_vol);
function [res] = CrossPowerSpectrum(A, B)
y = A.*conj(B);
scale = abs(y);
res = y./scale;
end
function [new_dim] = GetPad(vol)
dim = nextpow2(size(ref));
new_dim = 2.^dim;
end
function new_vol = upsample(vol)
m = 2;
[x,y,z] = size(vol);
new_vol = zeros(2*x,2*y,z);
for i = 1 : size(vol,3)
    new_vol(1:m:end,1:m:end,i) = vol(:,:,i);
end
sz = m*m;
H = fspecial('average',[sz sz]);
for i = 1 : size(vol,3)
    new_vol(:,:,i) = conv2(new_vol(:,:,i), H, 'same');
end
end

2 个评论
显示无隐藏无

Matt J 2020-1-27

But there are no FFT operations in the code that you have posted...??

raheem mian 2020-1-27

编辑：raheem mian 2020-1-27

My mistake. Added more code. I see Nan values in shift_vol. I see one NaN in cp_vol. I see a 0 in fftt but no Nans. I think the problem is the 0 in fftt. But I dont have this problem when I dont upsample the image. So I thought maybe I am upsampling incorrectly.

请先登录，再进行评论。

请先登录，再回答此问题。