does this function really need to be part of dipy? This doesn't seem general enough to be used across many users and it seems like this function could be replaced by a fairly start-forward, albeit long, one liner, `np.all(((cfa >= threshold[0::2]) & (cfa < threshold[1::2]) & ROI[..., None]), -1)`.