Upsampled image FFT contains NaNs.

4 次查看（过去 30 天）
raheem mian
raheem mian 2020-1-27
评论： 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
Matt J 2020-1-27
But there are no FFT operations in the code that you have posted...??
raheem mian
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.

请先登录，再进行评论。

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

采纳的回答

Matt J
Matt J 2020-1-27
编辑：Matt J 2020-1-27
You have a division by zero in this line,
res = y./scale;
That leads to NaNs in cp_vol and hence also in shift_vol.
  2 个评论
raheem mian
raheem mian 2020-1-27
Yeah I understand that but I need to do that calculation for phase correlation. I was wondering if my upsampling looks correct.
Matt J
Matt J 2020-1-27
I don't see anything in your upsample() function that would generate NaNs, though that would be easy to check.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Correlation and Convolution 的更多信息

标签

产品


版本

R2019b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by