Simplifying for-find loop functions to speed up processing

2024 1 26
2 个回答
回答已采纳
更新时间：2024 1 26
6 次查看（30 天）

0 个投票

Hi all,
A pretty basic question, but I'm trying to find a more elegant solution to search through a 750,000x5 list to remove entries from a corresponding list when more than three value in any respective row are above a threshold of 10. The long way of doing so that I have is
% mD is a 750000 x 5 matrix containing distances to the nearest 5
% neighbors of the r,c,v point (from knnsearch) in each row
for i = size(r,1):-1:1
if size(find(mD(i,:)>10),2) > 3
r(i) = [];
c(i) = [];
v(i) = [];
end
end

4 个评论
显示 2更早的评论 隐藏 2更早的评论

Dyuman Joshi
Dyuman Joshi 2024-1-26
How does the condition checked depend on the for loop index?
If it does not, then you are effectively deleting the elements (with indices - size(r,1):-1:1)) or not, which can be done directly or not done at all.
Cameron
Cameron 2024-1-26
This code goes line by to check on a per-point basis which I'm trying to speed up, it goes in -1 steps to not run out of indexes since it deletes as it goes. It depends on the for loop because the corresponding values in the distance list and the r,c,v list are the same index. In pseudo-code terms, its:
Are there more than 3 values in mD row (i) above 10?
If yes, delete row(i) from r, c, v
Torsten
Torsten 2024-1-26
Are there more than 3 values in mD row (i) above 10?
But you don't refer to row i of mD in your loop - you refer to the complete matrix mD with your find-command.
Cameron
Cameron 2024-1-26
Oops, good catch, that's my mistake, that should be mD(i,:)

请先登录，再进行评论。

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

 采纳的回答

Dyuman Joshi
Dyuman Joshi 2024-1-26
编辑：Dyuman Joshi 2024-1-26

1 个投票

I assume r, c and v have the same number of rows -
%Check which rows from the given range in mD have more than 3 values greater than 10
idx = sum(mD(1:size(r,1),:)>10, 2)>3
%perform deletion
r(idx) = [];
c(idx) = [];
v(idx) = [];
If all the variables have the same number or rows then you can remove the indices used and just use mD.

更多回答（1 个）

Cameron
Cameron 2024-1-26
编辑：Cameron 2024-1-26

0 个投票

I've developed a better solution to the above, which is
[~,mD] = knnsearch(rcv,rcv,'K',5);
mD_thresh = find(mD>10);
[row,~] = ind2sub(size(mD),mD_thresh);
for i = 1:size(row,1)
rcv(i,:) = [];
end
This only takes a few seconds to run, but it sacrifices testing whether more than 3 values in the row are over the 10 threshold, and just deletes any row with a single value over the limit (preferable to a 200 hour run though). I also still feel there's a way to eliminate that for loop.

类别

在 帮助中心 和 File Exchange 中查找有关 Loops and Conditional Statements 的更多信息

产品

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by