1. A method for noise reduction of seismic data, the method comprising: obtaining an input seismic data set acquired using one or more seismic sources and receivers; conducting a first inversion on a first frequency filtered input seismic data set across all slowness values in the plurality of slowness values using an inversion matrix with dipole modulation comprising a plurality of Tau-p operators and a plurality of dipole terms to transform the first frequency filtered input seismic data set from a time domain to a Tau-p domain comprising a plurality of first inversion Tau-p coefficients, each Tau-p coefficient associated with one of the plurality of slowness values; using the first inversion Tau-p coefficients to identify a subset of slowness values from the plurality of slowness values; performing iterative sparse Tau-p inversion using the identified subset of slowness values; and modifying the input seismic data set to generate a noise attenuated seismic data set based on the iterative sparse Tau-p inversion.