function T = quantile_landmarks(I, no_quants, t_range=[1, i_end]) Arguments: I is integrate(X) no_quants: number of quantiles/landmarks, always including 0% and 100%. Number of intervals will be one less than no_quants. t_range: for the final conversion of T's values to the proper scale. If not informed, it will be [1, size(I, 2)]