1

次のコードをベクトル化しようとしています。

% code before 
% code before 
% a lot of code before we got to the current comment 
% 
% houghMatrix holds some values 
for i=1:n
    for j=1:m
        for k = 1:maximalRadius            
            % get the maximal threshold 
            if houghMatrix(i,j,k) > getMaximalThreshold(k)                           
                lhs = [j i k];

                % verify that the new circle is not listed 
                isCircleExist = verifyCircleExists(circles,lhs,circleCounter);

                % not listed - then we put it in the circles vector 
                if isCircleExist == 0
                    circles(circleCounter,:) = [j i k];                    
                    fprintf('Circle % d: % d, % d, % d \n', circleCounter, j, i, k);
                    circleCounter = circleCounter + 1;                    
                end
            end
        end
    end
end

tic tac を使用すると、以下の出力が得られました。

>> x = findCircles(ii);
Circle  1:  38,  38,  35 
Circle  2:  89,  51,  34 
Circle  3:  72,  66,  11 
Circle  4:  33,  75,  30 
Circle  5:  90,  81,  31 
Circle  6:  54,  96,  26 

    Elapsed time is 3.111176 seconds.
>> x = findCircles(ii);
Circle  1:  38,  38,  35 
Circle  2:  89,  51,  34 
Circle  3:  72,  66,  11 
Circle  4:  33,  75,  30 
Circle  5:  90,  81,  31 
Circle  6:  54,  96,  26 

    Elapsed time is 3.105642 seconds.
>> x = findCircles(ii);
Circle  1:  38,  38,  35 
Circle  2:  89,  51,  34 
Circle  3:  72,  66,  11 
Circle  4:  33,  75,  30 
Circle  5:  90,  81,  31 
Circle  6:  54,  96,  26 

    Elapsed time is 3.135818 seconds.

意味 - 平均 3.1 秒。

コードをベクトル化しようとしましたが、問題は、内部(3 番目)i,j,kの本体でインデックスを使用する必要があることです。forfor

コードをベクトル化する方法についての提案は大歓迎です

ありがとう

編集 ： 

% -- function [circleExists] = verifyCircleExists(circles,lhs,total) --
%
%
function [circleExists] = verifyCircleExists(circles,lhs,total)

    MINIMUM_ALLOWED_THRESHOLD = 2;

    circleExists = 0;
    for index = 1:total-1                
        rhs = circles(index,:);
        absExpr = abs(lhs - rhs);
        maxValue = max( absExpr );
        if  maxValue <= MINIMUM_ALLOWED_THRESHOLD + 1
            circleExists = 1;
            break
        end
    end

end
4

1 に答える 1

0

やりたいことは次のとおりです。有効な座標トリプレットごとに、近くのトリプレットがすでにあるかどうかを確認し、そうでない場合はリストに追加します。「連鎖」の可能性がない場合、つまり、可能な候補ボクセルの各クラスターが 1 つの中心しか収容できない場合、この操作は完全にベクトル化できます。この場合、次を使用するだけです。

%# create a vector of thresholds
maximalThreshold = getMaximalThreshold(1:maximalRadius);

%# make it 1-by-1-by-3
maximalThreshold = reshape(maximalThreshold,1,1,[]);

%# create a binary array the size of houghMatrix with 1's 
%# wherever we have a candidate circle center
validClusters = bsxfun(@gt, houghMatrix, maximalThreshold);

%# get the centroids of all valid clusters
stats = regionprops(validClusters,'Centroid');

%# collect centroids, round to get integer pixel values
circles = round(cat(1,stats.Centroid));

または、有効な円を選択するスキームに従いたい場合は、次のように ijk インデックスを取得できvalidClustersます。

[potentialCircles(:,1),potentialCircles(:,2), potentialCircles(:,3)]= ...  
       sub2ind(size(houghMatrix),find(validClusters));

nPotentialCircles = size(potentialCircles,1);


for iTest = 2:nPotentialCircles
    absDiff = abs(bsxfun(@minus,potentialCircles(1:iTest-1,:),potentialCircles(iTest,:)));

    if any(absDiff(:) <= MINIMUM_ALLOWED_THRESHOLD + 1)
       %# mask the potential circle
       potentialCircles(iTest,:) = NaN;
    end
end

circles = potentialCircles(isfinite(potentialCircles(:,1)),:);
于 2013-01-05T21:36:47.633 に答える