View difference between Paste ID: WkS5ShMf and 5KzhvYbU
SHOW: | | - or go back to the newest paste.
1
    [rows, cols, scls] = size(scalespace);
2
    scalespace_feats = cell(scls,1);
3
    [XX, YY] = meshgrid(1:cols, 1:rows);
4
    for rr=1:rows
5
        for cc=1:cols
6
            distance = sqrt( (YY-rr).^2 + (XX-cc).^2  );
7
8
            for ss=1:scls
9
                % scalespace_radii contaings the radius associated to a given scale, i.e.: radii = scale * sqrt(2)
10
                mask_feat_radii = (distance <= scalespace_radii(ss));
11
                currscale       = scalespace(:,:,ss);
12
                responses       = currscale(mask_feat_radii);
13
14
                scalespace_feats{ss}(rr,cc,:) = [mean(responses), std(responses), median(responses)];                
15
            end
16
            fprintf('\n');
17
        end
18
    end