我有一个名为grnPixels
的大小为(1 x 40)
的单元格数组,其中每个单独的单元格都有一个M x 1
向量数组,其中M
是可变的。我还有一个名为redCentroid
的矢量数组N x 1
。
我想检查redCentroid
中的值是否与grnPixels
中的任何值相对应。我已经制作了一个代码但是在这个Matlab代码中非常慢。我怎样才能改善这个?
nRedCells = length(propsRed);
nGrnCells = length(propsGrn);
grnPixels = cell(1,nGrnCells);
redCentroid = zeros(nRedCells,1);
matchMemory = zeros(nRedCells,1);
for j = 1:nRedCells
for i = 1:nGrnCells
for n = 1:length(grnPixels{i})
matchment = ismember(redCentroid(j),grnPixels{i}(n));
if matchment == 1
matchMemory(j,1:2) = [j i];
end
continue
end
end
end
样本数据
redCentroid
51756
65031
100996
118055
122055
169853
197175
233860
244415
253822
grnPixels{1}
142
143
100996
167
168
grnPixels{2}
537
538
539
540
541
542
233860
244415
545
546
547
548
答案 0 :(得分:1)
ismember
可以接受第一个或第二个输入的矩阵,因此不需要外部循环或最内部循环。
matchMemory = zeros(numel(redCentroid), 2);
for k = 1:numel(grnPixels)
% Check which of the centroids are in grnpixels
isPresent = ismember(redCentroid, grnPixels{k});
% If they were present fill in the first column with the index in red
matchMemory(isPresent, 1) = find(isPresent);
% Fill in the second column with the index in green
matchMemory(isPresent, 2) = k;
end
答案 1 :(得分:1)
如果您想查找任何匹配而不管订单
这应该在2 * M * log(M)+ 2 * M时间内运行。
如果要查找原始矩阵中与匹配项对应的索引,只需将收集器数组中的每个元素与两个矩阵的元素进行比较,每当找到匹配项时记录索引,并继续直到达到端。
如果元素按特定顺序(如坐标),只需将第一组中的元素1与第二组中的元素1进行比较。