Vectorize a for loop

4 ビュー (過去 30 日間)
Ramon Ticzon
Ramon Ticzon 2019 年 3 月 19 日
コメント済み: Ramon Ticzon 2019 年 3 月 19 日
Hi. How can I speed up the process of my for loop code.
function [missingCase] = CaseChecking(bimage)
missingCase = '';
for img = 1:66
mapList = xlsread('map.xlsx',img);
counter=0;
for i = 1:200
if (~bimage(mapList(i,2),mapList(i,1)))
counter=counter+1;
end
end
if counter<50
missingCase = img;
break
end
end

回答 (2 件)

Andrei Bobrov
Andrei Bobrov 2019 年 3 月 19 日
Please attach your data (bimage,map.xlsx - small example) as mat-file and xlsx-file.
ml = zeros(200,2,66);
for ii = 1:size(ml,3)
ml(:,:,ii) = xlsread('map.xlsx',ii);
end
mlw = permute(ml(:,end:-1:1,:),[1,3,2]);
out = find(~bimage(sub2ind(size(bimage),mlw(:,:,1),mlw(:,:,2))) < 50,'first',1);
  1 件のコメント
Ramon Ticzon
Ramon Ticzon 2019 年 3 月 19 日
These are the sample files. Thank you

サインインしてコメントする。


Andrei Bobrov
Andrei Bobrov 2019 年 3 月 19 日
Solution for attached data:
[~,sheets] = xlsfinfo('MapList.xlsx');
n = numel(sheets);
ml = cell(n,1);
for ii = 1:n
ml{ii} = xlsread('MapList.xlsx',ii);
end
jj = repelem((1:n)',cellfun('size',ml,1));
A = cat(1,ml{:});
k = accumarray(jj,(1:size(A,1))',[],...
@(x) sum(~bimage(sub2ind(size(bimage),A(x,2),A(x,1)))));
out = find(k < 50,'first',1);
  2 件のコメント
Ramon Ticzon
Ramon Ticzon 2019 年 3 月 19 日
where is the variable 'missingCase'? also the 'img' variable.
Ramon Ticzon
Ramon Ticzon 2019 年 3 月 19 日
Hi. Im having an error.
Error in CaseChecking (line 14)
@(x) sum(~bimage(sub2ind(size(bimage),A(x,2),A(x,1)))));

サインインしてコメントする。

カテゴリ

Help Center および File ExchangeLoops and Conditional Statements についてさらに検索

Community Treasure Hunt

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

Start Hunting!

Translated by