Created
November 5, 2017 23:53
-
-
Save alimanfoo/c5977e87111abe8127453b21204c1065 to your computer and use it in GitHub Desktop.
Find runs of consecutive items in a numpy array.
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
import numpy as np | |
def find_runs(x): | |
"""Find runs of consecutive items in an array.""" | |
# ensure array | |
x = np.asanyarray(x) | |
if x.ndim != 1: | |
raise ValueError('only 1D array supported') | |
n = x.shape[0] | |
# handle empty array | |
if n == 0: | |
return np.array([]), np.array([]), np.array([]) | |
else: | |
# find run starts | |
loc_run_start = np.empty(n, dtype=bool) | |
loc_run_start[0] = True | |
np.not_equal(x[:-1], x[1:], out=loc_run_start[1:]) | |
run_starts = np.nonzero(loc_run_start)[0] | |
# find run values | |
run_values = x[loc_run_start] | |
# find run lengths | |
run_lengths = np.diff(np.append(run_starts, n)) | |
return run_values, run_starts, run_lengths |
Glad it's useful!
Btw, if anyone needs an artihmetic RLE library, check out https://github.com/pyranges/pyrle :)
Btw, if anyone needs an artihmetic RLE library, check out https://github.com/pyranges/pyrle :)
👍
Thank you :)
Thanks! This is exactly what I needed!
Thanks! Nice work, just what I want!
膜拜大佬
Very good! Thanks.
Thanks finally a solution that works as intended
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
This is great! Just what I needed, thank you. Included in my research code with attribution.