Substract nd array slices

Question

Orangelynx 2017-11-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/366568-substract-nd-array-slices

评论： Orangelynx 2017-11-13

Hi,

I have a n-dimensional array which may contain singleton dimensions and I would like to substract the second page of dimension i from the first page of dimension i. My current solution looks something like this:

function A = substract_page( A, dim )
  orig_size = size(A);    % store the original shape
  A = shiftdim(A, dim);   % shift the dimension to be substracted to the end
  shifted_size = size(A); % aliter: shifted_size = circshift orig_size(dim); ?
  A = reshape(A, [], shifted_size(end));  % collapse first dimensions for easier substraction
  A = A(:, 2) - A(:, 1);
  shifted_size(end) = shifted_size(end) - 1;  % substraction leads to dimension loss
  A = reshape(A, shifted_size);               % restore shape
  lost_dimensions = length(orig_size) - ndims(A);
  A = shiftdim(A, -lost_dimensions);         % make sure that  singular dimensions  in the end don't get lost
  A = shiftdim(A, length(orig_size) - dim + lost_dimensions);    % restore order
end

But it seems way to complicated and I feel there is a much easier and more elegant solution.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Jan 2017-11-12

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/366568-substract-nd-array-slices#answer_290721

编辑：Jan 2017-11-12

在 MATLAB Online 中打开

Does your code perform this:

dim = 3;
A   = rand(2,3,4,5);
A   = A(:, :, 2, :) - A(:, :, 1, :)

? Then:

function A = substract_page(A, dim)
C1      = repmat({':'}, 1, ndims(A));
C2      = C1;
C1{dim} = 1;
C2{dim} = 2;
A       = A(C2{:}) - A(C1{:});
end

3 个评论
显示 1更早的评论隐藏 1更早的评论

Jan 2017-11-13

@Orangelynx: When the dimension to operate on has the length 2, diff is perfect. My code would be useful only, if the input has e.g. the length 4 as in my example, and only the difference between the first 2 slices should be calculated ignoring the rest.

Orangelynx 2017-11-13

good point, I guess I didn't specify that clearly in my original question.

请先登录，再进行评论。

Answer 2

Matt J 2017-11-12

3
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/366568-substract-nd-array-slices#answer_290702

在 MATLAB Online 中打开

Not sure, but maybe you are looking for

A=diff(A,1,dim);

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Orangelynx 2017-11-13

wow.. I knew about diff, but for some reason my head told me it's not what I was looking for, eventhough with the dim parameter, it's exactly what I wanted... thanks.

请先登录，再进行评论。