1. A method for deblending seismic data associated with a subsurface of the earth, the method comprising: receiving an input dataset generated by first and second sources S 1 and S 2 that are operating as simultaneous sources; arranging the input dataset based on the firing times of source S 1 ; applying with a computing system an annihilation filter to the arranged input dataset to estimate cross-talk noise; convolving the cross-talk noise estimate with an operator to form a signal estimate using the firing times of S 1 and S 2 ; and generating an image of the subsurface based on the signal estimate, wherein the annihilation filter removes coherent energy while a coherency filter removes incoherent energy.