28
28
local_sum_A2 = dic_local_sum(double_sub_base.*double_sub_base,n,n);
29
29
denom_A = sqrt( max(( local_sum_A2 - (local_sum_A.^2)/mn ) / (mn-1), 0) );
30
30
i_nonzero = find(denom_A~=0);
31
normxcorr_hw(hwid, 2, icp, sub_base, single(local_sum_A), single(denom_A), uint16(i_nonzero));
32
[lsum, denom] = normxcorr_hw(hwid, 2, icp, base);