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