codebase / stim-matlab

  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