Вы можете использовать numpy.r_
для объединения диапазонов со скалярами.Единственное осложнение заключается в том, что вам нужно использовать pd.DataFrame.iloc
вместо этого, но это можно сделать с помощью df.columns.get_loc
.
Вот демонстрационная версия:
import pandas as pd
import numpy as np
df = pd.DataFrame(columns=['column'+str(i) for i in range(1, 82)])
colidx = df.columns.get_loc
res = df.iloc[:, np.r_[colidx('column1'):colidx('column5'), colidx('column80')]]
print(res.columns)
Index(['column1', 'column2', 'column3', 'column4', 'column80'], dtype='object')