Recovers saturated spectral signals by modeling data as a low-rank bilinear system, flagging clipped channels as missing (NaNs), estimating scores and loadings with Orthogonalized Alternating Least Squares (O-ALS) at a known rank, and reconstructing the saturated regions from those estimated scores and loadings.
