We present techniques developed for numerical modeling of wave propagation, and source-signature removal in seismic imaging, based on a class of linear operators known as Gabor multipliers. These operators are localized Fourier multipliers, whose actions is selectively localized by an element of a partition of unity. We discuss boundedness and stability properties for these operators, approximations to PDEs and pseudodifferential operators, and an approximate functional calculus.