сложение матрицы с пропуском значений NaN

Question

Я хотел бы выполнить поэлементное сложение матриц, пропуская НЭН ценности.MATLAB и Octave имеют нансум, но он выполняет сложение по столбцам внутри матрицы.

Позволять:

a = NaN * zeros(3)
b = ones(3)

Я хочу:

c = nan+(a, b)

и

c = b

и:

d = nan+(a,a)

и

d = a

Answer 1

Вы все равно можете использовать nansum, если объедините свои массивы n-d по n+1-му измерению.

Для 2D

% commands de-nested for readability. You can do this with a single line, of course
tmp = cat(3,a,b);
c = nansum(tmp,3);

Общий случай

function out = nansumByElement(A,B)
%NANSUMBYELEMENT performs an element-wise nansum on the n-D arrays A and B
% A and B have to have the same size

% test input
if nargin < 2 || isempty(A) || isempty(B) || ndims(A)~=ndims(B) || ~all(size(A)==size(B))
error('please pass two non-empty arrays of the same size to nansumByElement')
end

% calculate output

nd = ndims(A); % get number of dimensions
% catenate and sum along n+1st dimension
out = nansum(cat(nd+1,A,B),nd+1);

Answer 2

a_fixed = a;
a_fixed(isnan(a)) = 0;
b_fixed = b;
b_fixed(isnan(b)) = 0;
c = a_fixed.+b_fixed;