Methods are presented for calibrating or characterizing a printing system with respect to at least one color, in which a toner reproduction curve (TRC) is measured and curve-fitted according to a least squares solution using a set of spline basis functions having high spline density in regions of high TRC curvature and/or high measurement noise, with the weights of the spline functions being restricted to positive values to maintain monotonicity of the TRC.