codebase / stimlib

function new_mask = rtsCropMask(mask, max_pixels)

%convert to logical
%mask = mask > 0;

if nnz(mask) < max_pixels
    new_mask = mask;
else
    j = 0;
    for i = 1:length(mask(:))
        if mask(i) > 0
            j = j + 1;
            if j == max_pixels
                last_pixel = i;
            end
        end
    end
    new_mask = zeros(size(mask, 1), size(mask, 2));
    new_mask(1:last_pixel) = mask(1:last_pixel);
end