arr = np.array([1, 2, 3, 15, 16, 17, 17, 18, 19, 30, 41, 42, 43, 55, 56])
bin_edges = arr[np.where(np.diff(arr) > 10)[0]]
indices = np.digitize(arr, bin_edges, right=True)
print(indices)
Увеличьте элементы массива numpy на 1 по определенным индексам (для использования с группировкой таблицы astropy)