Scattered light noise affects the sensitivity of gravitational waves detectors. The characterization of such noise is needed to mitigate it. The time-varying filter empirical mode decomposition algorithm is suitable for identifying signals with time-dependent frequency such as scattered light noise (or scattering). We present a fully automated pipeline based on the pytvfemd library, a python implementation of the tvf-EMD algorithm, to identify objects inducing scattering in the gravitational-wave channel with their motion. The pipeline application to LIGO Livingston O3 data shows that most scattering noise is due to the penultimate mass at the end of the X-arm of the detector (EXPUM) and with a motion in the micro-seismic frequency range.