We introduce a new algorithm for bidirectional
deconvolution. In our method, we
estimate the causal filters and anti-causal filters simultaneously
instead of alternately. We test three data examples (1D synthetic, 2D synthetic and 2D field data). The results show that the wavelet can
be compressed almost into a spike using our method. The two filters can be
estimated equally when we are dealing with a zero-phase wavelet. In addition, our method has a lower computational cost and faster convergence rate than the method discussed by
Zhang and Claerbout (2010).