numpy.ma.take

ma.take(a, indices, axis=None, out=None, mode='raise')