diff --git a/jarvis/core/utils.py b/jarvis/core/utils.py index 69e0298b..97d65e85 100644 --- a/jarvis/core/utils.py +++ b/jarvis/core/utils.py @@ -362,7 +362,18 @@ def baseline_als(y, lam, p, niter=10): z = spsolve(Z, w*y) w = p * (y > z) + (1-p) * (y < z) return z - + +def recast_array(x_original = [],y_original = [],x_new = np.arange(0, 90, 1),tol=0.1): + x_original=np.array(x_original) + # Initialize the new y array with NaNs or a default value + y_new = np.full_like(x_new, 0, dtype=np.float64) + + # Fill the corresponding bins + for x_val, y_val in zip(x_original, y_original): + closest_index = np.abs(x_new - x_val).argmin() # Find the closest x_new index + y_new[closest_index] = y_val + #y_new[y_new