沿3维对3D数组排序

Question

我有一个144x91x92矩阵存储在o3{1} 。 我想浏览第三个维度（代表一个季节中的天数）并对其进行排序。 然后，我想提取值的前10个百分位数和后10个百分位数，包括这些值的索引。 这将找到每个网格单元的前10天和最后10％天（每个网格单元会有所不同）。

我正在尝试以下方法：

[Y,I] = sort(o3{1},3); % sort o3 along 3rd dimension
o3_sorted = o3{1}(I);
ind_top10 = o3_sorted(90);
ind_bot10 = o3_sorted(10);    

但是我知道我并没有提出正确的前10个百分点和前10个百分点。 另外，这种方式不会告诉我前10天和下10天百分位的索引（每个144x91网格单元的索引不同）。 我希望最终得到144x91x10矩阵，以表示前10％的天数，后10％的天数以及每个索引。

Answer 1

像这样尝试：

[~, I] = sort(o3{1},3); %// sort o3 along 3rd dimension
ind_top10 = I(:,:,end-8:end);
ind_bot10 = I(:,:,1:9);

I3 = cat(3, ind_top10, ind_bot10); %// you might want to skip this but and just work with the top and bottom separately from here on

[I1, I2, ~] = ndgrid(1:size(o3{1},1), 1:size(o3{1},2), 1:size(I3,3));
ind = sub2ind(size(o3{1}),I1,I2,I3)

现在

o3{2}(ind)
o3{3}(ind)
%// etc...

Answer 2

这与Dan提出的想法略有不同：

%% // Init
clear variables; clc;

%% // Generate some data:
o3 = cell(3,1);
for indO = 1:3
  o3{indO} = randi(intmax('uint16'),144,91,92,'uint16');
end

%% // Find 10 & 90 percentiles:
percentiles = cat(3,prctile(o3{1},10,3),prctile(o3{1},90,3));

%% // Find indices of relevant values
select_bot10 = bsxfun(@ge,o3{1},percentiles(:,:,1)); %// replace @ge with @gt if needed
select_top10 = bsxfun(@le,o3{1},percentiles(:,:,2)); %// replace @le with @lt if needed
%// Another optional way to index the values if required:
[rb,cb,vb] = ind2sub(size(select_bot10),find(select_bot10));
[rt,ct,vt] = ind2sub(size(select_top10),find(select_top10));

%% // Get values from o3{1..3} etc.
bot10 = o3{1}(select_bot10);
top10 = o3{1}(select_top10);
%// etc.

此解决方案可能不适合您的特定需求，但改编应该很简单。 还要注意，由于采用了精确的百分位数，因此bot10和top10之间的元素数量可能会有所不同。

3D查找的功劳归功于 Kenneth Eaton / gnogno 。